Simple cation-pi interaction between a phenyl ring and a protonated amine stabilizes an alpha-helix in water

Cation-pi interactions have been proposed to be important contributors to protein structure and function. In particular, these interactions have been suggested to provide significant stability at the solvent-exposed surface of a protein. We have investigated the magnitude of cation-pi interactions between phenylalanine (Phe) and lysine (Lys), ornithine (Orn), and diaminobutanoic acid (Dab) in the context of an alpha-helix and have found that only the Phe...Orn interaction provides significant stability to the helix, stabilizing it by -0.4 kcal/mol. This interaction energy is in the same range as a salt bridge in an alpha-helix, and equivalent to the recently reported Trp...Arg interaction in an alpha-helix, despite the fact that Trp...guanidinium interactions have been proposed to be stronger than Phe...ammonium interactions. These results indicate that even the simplest cation-pi interaction can provide significant stability to protein structure and demonstrate the subtle factors that can influence the observed interaction energies in designed systems.     

Retention behavior of peptides in hydrophilic-interaction chromatography

Selected hydrophilic interaction chromatography (HILIC) columns packed with bare silica, bridge-ethyl hybrid silica, or an amide sorbent chemistry were utilized for an investigation of chromatographic behavior and separation selectivity of tryptic peptides. Retention model was proposed allowing for retention prediction of peptides with correlation coefficient R(2)~0.92-0.97 for various columns. The values of optimized amino acid retention coefficients were compared to those obtained for reversed-phase liquid chromatography (Gilar et al., Anal. Chem. 2010, 82, 265-275) and used to elucidate the impact of different amino acid on peptide HILIC retention. In contrast to reversed-phase chromatography, where presence of Phe, Trp, Ile, and Leu amino acid residues in sequence strongly promoted, and presence of hydrophilic His, Lys and Arg residues strongly reduced peptide retention, the effects of these amino acid residues in HILIC were opposite (His, Lys and Arg promote, Phe, Trp, Ile and Leu demote peptide retention in HILIC). Retention coefficient optimized for pH experiments illustrated the impact of silanols on HILIC retention.     

Putative One‐Pot Prebiotic Polypeptides with Ribonucleolytic Activity

KIA7, a peptide with a highly restricted set of amino acids (Lys, Ile, Ala, Gly and Tyr), adopts a specifically folded structure. Some amino acids, including Lys, Ile, Ala, Gly and His, form under the same putative prebiotic conditions, whereas different conditions are needed for producing Tyr, Phe and Trp. Herein, we report the 3D structure and conformational stability of the peptide KIA7H, which is composed of only Lys, Ile, Ala, Gly and His. When the imidazole group is neutral, this 20‐mer peptide adopts a four‐helix bundle with a specifically packed hydrophobic core. Therefore, one‐pot prebiotic proteins with well‐defined structures might have arisen early in chemical evolution. The Trp variant, KIA7W, was also studied. It adopts a 3D structure similar to that of KIA7H and its previously studied Tyr and Phe variants, but is remarkably more stable. When tested for ribonucleolytic activity, KIA7H, KIA7W and even short, unstructured peptides rich in His and Lys, in combination with Mg⁺⁺, Mn⁺⁺ or Ni⁺⁺ (but not Cu⁺⁺, Zn⁺⁺ or EDTA) specifically cleave the single‐stranded region in an RNA stem–loop. This suggests that prebiotic peptide–divalent cation complexes with ribonucleolytic activity might have co‐inhabited the RNA world.     

A two model receptor system of the alpha1D adrenergic receptor to describe interactions with epinephrine and BMY7378

In this study, we have developed a two receptor model system to describe the R and R states of G-protein coupled receptors, specifically the alpha(1D) adrenergic receptor. The two models interact with agonist (epinephrine) and antagonist (BMY7378) differently. The active model has increased interactions with epinephrine. The inactive model has increased interactions with BMY7378. We also explored the protonation state of the ligands. When the most basic amine was protonated, we found increased hydrogen bonding and increased aromatic interactions. Protonated epinephrine hydrogen bonds with Asp176 and has aromatic residues Trp172, Trp235, Trp361, and Phe388 within 3 Angstroms. Protonated BMY7378 hydrogen bonds with Trp172 and Lys236 and has aromatic residues Trp172, Trp254, Phe364, Phe384, and Phe388 within 3 Angstroms. We conclude that the two model system is required to represent the two states of the receptor and that protonation of the ligand is also critical.     

Patterns of cation binding to the aromatic amino acid R groups in Trp,Tyr, and Phe

Previous joint experimental and theoretical work demonstrates that typically soluble peptides will be rendered insoluble in the presence of saturated sodium ions in aqueous solution due to disruption of cation-π interactions between Trp and Lys. The present work utilizes quantum chemical methods including density functional theory, symmetry-adapted perturbation theory, and even coupled cluster theory to determine the strengths of cation-π interactions for the aromatic R groups of Trp, Tyr, and Phe (approximated as skatole, methyl phenol, and toluene) with both alkali and alkaline-Earth atomic cations and electron-accepting R groups from Lys, Arg, and His approximated as methyl ammonium, guanidinium, and imidazolium cations. This work shows that sodium ion is still the most likely disrupter of peptide folding built upon cation-π interactions, since Trp, Tyr, and Phe all bind more strongly to sodium ion than to any of the polyatomic cations. Additionally, the atomic cation complex binding energies decrease with an increase in partial charge on the atomic cation in the complex. However, as the average partial charge increases in the interacting hydrogen atoms in the polyatomic cations, the binding energy increases. The disruption of such peptide–peptide cation-π interactions is certainly relevant for peptide design in β-sheets or β-hairpin structures, but it could also have implications for astrobiology.     

Design of a bacterial host for site-specific incorporation of p-bromophenylalanine into recombinant proteins

Introduction of a yeast suppressor tRNA (ytRNA(Phe)(CUA)) and a mutant yeast phenylalanyl-tRNA synthetase (yPheRS (T415G)) into an Escherichia coli expression host allows in vivo incorporation of phenylalanine analogues into recombinant proteins in response to amber stop codons. However, high-fidelity incorporation of non-natural amino acids is precluded in this system by mischarging of ytRNA(Phe)(CUA) with tryptophan (Trp) and lysine (Lys). Here we show that ytRNA(Phe)(CUA) and yPheRS can be redesigned to achieve high-fidelity amber codon suppression through delivery of p-bromophenylalanine (pBrF). Two strategies were used to reduce misincorporation of Trp and Lys. First, Lys misincorporation was eliminated by disruption of a Watson-Crick base pair between nucleotides 30 and 40 in ytRNA(Phe)(CUA). Loss of this base pair reduces mischarging by the E. coli lysyl-tRNA synthetase. Second, the binding site of yPheRS was redesigned to enhance specificity for pBrF. Specifically, we used the T415A variant, which exhibits 5-fold higher activity toward pBrF as compared to Trp in ATP-PP(i) exchange assays. Combining mutant ytRNA(Phe)(CUA) and yPheRS (T415A) allowed incorporation of pBrF into murine dihydrofolate reductase in response to an amber codon with at least 98% fidelity.     

Differentiation of the local structure around tryptophan 51 and 64 in recombinant human erythropoietin by tryptophan phosphorescence

Abstract Recombinant human erythropoietin is a 4-helix bundle, glycosylated cytokine containing three tryptophan residues at positions 51, 64 and 88 whose phosphorescence emission may represent a sensitive probe of the structure at multiple sites near or at the protein surface. This report characterizes the phosphorescence properties (spectral energy, thermal spectral relaxation and phosphorescence lifetime), from low temperature glasses to ambient temperature, of the native protein plus that of three single point mutation analogs where each Trp was replaced by Phe. The structural information inferred from the phosphorescence parameters was essentially in good agreement with the structure of the Escherichia coli-produced nonglycosylated protein determined by nuclear magnetic resonance (Cheetham et al., Nat Struct Biol [1998] 5:861). The results showed that the fluorescence and phosphorescence spectra of the native protein were entirely due to independent contributions of Trp51 and Trp64 and that Trp88 was quenched under all conditions. The phosphorescence emissions of Trp51 and Trp64 were differentiated by their unique spectra at 77 K with Trp64 exhibiting an unusually blueshifted spectrum likely due to the attractive interaction of Arg110 and Lys116 with the ground state dipole of Trp64. In the native protein the room temperature phosphorescence lifetime of Trp64 was relatively short with a time of 1.62 ms whereas the lifetime of Trp51 was five-fold longer. Characterization of the single point mutation analogs showed that each lifetime was composed of multiple components revealing the presence of multiple stable conformations of the protein at these surface sites.     

The role of cation-pi interactions in biomolecular association. Design of peptides favoring interactions between cationic and aromatic amino acid side chains.

Cation-pi interactions between amino acid side chains are increasingly being recognized as important structural and functional features of proteins and other biomolecules. Although these interactions have been found in static protein structures, they have not yet been detected in dynamic biomolecular systems. We determined, by (1)H NMR spectroscopic titrations, the energies of cation-pi interactions of the amino acid derivative AcLysOMe (1) with AcPheOEt (2) and with AcTyrOEt (3) in aqueous and three organic solvents. The interaction energy is substantial; it ranges from -2.1 to -3.4 kcal/mol and depends only slightly on the dielectric constant of the solvent. To assess the effects of auxiliary interactions and structural preorganization on formation of cation-pi interactions, we studied these interactions in the association of pentapeptides. Upon binding of the positively-charged peptide AcLysLysLysLysLysNH(2) (5) to the negatively-charged partner AcAspAspXAspAspNH(2) (6), in which X is Leu (6a), Tyr (6b), and Phe (6c), multiple interactions occur. Association of the two pentapeptides is dynamic. Free peptides and their complex are in fast exchange on the NMR time-scale, and 2D (1)H ROESY spectra of the complex of the two pentapeptides do not show intermolecular ROESY peaks. Perturbations of the chemical shifts indicated that the aromatic groups in peptides 6b and 6c were affected by the association with 5. The association constants K(A) for 5 with 6a and with 6b are nearly equal, (4.0 +/- 0.7) x 10(3) and (5.0 +/- 1.0) x 10(3) M(-)(1), respectively, while K(A) for 5 with 6c is larger, (8.3 +/- 1.3) x 10(3) M(-)(1). Molecular-dynamics (MD) simulations of the pentapeptide pairs confirmed that their association is dynamic and showed that cation-pi contacts between the two peptides are stereochemically possible. A transient complex between 5 and 6 with a prominent cation-pi interaction, obtained from MD simulations, was used as a template to design cyclic peptides C(X) featuring persistent cation-pi interactions. The cyclic peptide C(X) had a sequence in which X is Tyr, Phe, and Leu. The first two peptides do, but the third does not, contain the aromatic residue capable of interacting with a cationic Lys residue. This covalent construct offered conformational stability over the noncovalent complexes and allowed thorough studies by 2D NMR spectroscopy. Multiple conformations of the cyclic peptides C(Tyr) and C(Phe) are in slow exchange on the NMR time-scale. In one of these conformations, cation-pi interaction between Lys3 and Tyr9/Phe9 is clearly evident. Multiple NOEs between the side chains of residues 3 and 9 are observed; chemical-shift changes are consistent with the placement of the side chain of Lys3 over the aromatic ring. In contrast, the cyclic peptide C(Leu) showed no evidence for close approach of the side chains of Lys3 and Leu9. The cation-pi interaction persists in both DMSO and aqueous solvents. When the disulfide bond in the cyclic peptide C(Phe) was removed, the cation-pi interaction in the acyclic peptide AC(Phe) remained. To test the reliability of the pK(a) criterion for the existence of cation-pi interactions, we determined residue-specific pK(a) values of all four Lys side chains in all three cyclic peptides C(X). While NOE cross-peaks and perturbations of the chemical shifts clearly show the existence of the cation-pi interaction, pK(a) values of Lys3 in C(Tyr) and in C(Phe) differ only marginally from those values of other lysines in these dynamic peptides. Our experimental results with dynamic peptide systems highlight the role of cation-pi interactions in both intermolecular recognition at the protein-protein interface and intramolecular processes such as protein folding.     

The recognition of nucleotides with model beta-hairpin receptors: investigation of critical contacts and nucleotide selectivity

[reaction: see text] We have investigated the factors that contribute to binding of ATP by a designed 12-residue beta-hairpin peptide, WKWK, and have determined its selectivity for binding to the naturally occurring nucleotide triphosphates. We have previously shown that WKWK creates an ATP binding pocket on one face of the beta-hairpin consisting of two Trp and two Lys residues. Mutation of the two Lys residues on the binding face of the beta-hairpin resulted in a lower affinity, indicating that each is involved in ATP binding and that each residue contributes approximately -1.5 kcal/mol to the energy of complexation. Replacement of either Trp residue of the ATP binding pocket with Phe or Leu destabilizes the complex formed with ATP by approximately 1 kcal/mol, indicating that both Trp residues participate in interactions with ATP. For binding to the nucleotide triphosphates, the order of binding affinity was shown to follow dTTP > GTP > ATP > CTP, with differences in binding energies spanning as much as 1.6 kcal/mol. NMR analysis demonstrates that both aromatic interactions with the Trp side chains and CH-pi interactions between the ribose protons and the Trp residues may contribute significantly to binding. The results from our model system provide useful thermodynamic information regarding protein-nucleic acid interactions that occur at the surface of a beta-sheet.     

Minimalist protein design: a beta-hairpin peptide that binds ssDNA

A 28-residue beta-hairpin dimer (WKWK)2 with two Trp and two Lys residues on one face of each beta-sheet was shown to form a complex with single-stranded oligonucleotides at low micromolar concentrations. Each beta-hairpin of the dimer contains a cross-strand Trp-Trp pair in a diagonal orientation which has previously been shown to create a cleft for the intercalation of aromatic guests such as adenine (J. Am. Chem. Soc. 2003, 125, 9580). The beta-hairpin dimer binds 5-residue ssDNA sequences 5'-AAAAA-3', 5'-TTTTT-3', and 5'-CCCCC-3' in water with dissociation constants in the range of 12-30 muM. A weak energetic preference for binding to sequence 5'-AAAAA-3' was observed, which is believed to result from stronger stacking interactions between Trp and the adenine base. The interaction of 5'-AAAAA-3' with the Lys and Trp residues of the peptide was evident by NMR, and a 1:1 association was demonstrated. The recognition of an 11-residue ssDNA sequence occurred with a dissociation constant of 3 muM under near-physiological ionic strength and pH, demonstrating that the beta-hairpin dimer binds ssDNA as strongly as many naturally occurring proteins. The salt dependence of the interaction of the 11-residue oligonucleotide with the peptide dimer indicates that Trp-nucleobase stacking interactions contribute about -4 kcal/mol to recognition, which is much greater than the contribution of nonionic interactions in unstructured peptides containing Trp. Moreover, recognition of the ssDNA demonstrated reduced salt dependence relative to the corresponding duplex, resulting in selectivity for ssDNA under high salt conditions. Peptide (WKWK)2 is a relevant mimic of OB-fold (oligonucleotide/oligosaccharide-binding) proteins which bind ssDNA on the surface of a beta-sheet.     

Interaction between PEO-PPO-PEO copolymers and a hexapeptide in aqueous solutions

Interaction between PEO-PPO-PEO copolymers and a hexapeptide, growth hormone releasing peptide-6 (GHRP-6), was investigated by NMR to study the potential use of the copolymers in peptide drug delivery. (1)H NMR and nuclear Overhauser effect spectroscopy (NOESY) measurements determined that PO methyl protons interacted with methyl protons of the Ala moiety, aromatic protons of the Trp moiety, and some of the Phe aromatic protons. The Lys moiety and part of the Phe moiety entered the hydrophilic EO environment via hydrogen bonding. PEO-PPO-PEO copolymers and the peptide formed a complex in 1:1 stoichiometry. Binding constants between copolymers and GHRP-6 were determined, and it was indicated that the copolymers containing more EO and PO units will lead to greater affinity with the peptide. Isothermal titration calorimetry (ITC) measurements confirmed the results of NMR experiments. This study indicates that PEO-PPO-PEO copolymers have great potential in delivering peptide drugs.     

Use of a retroinverso p53 peptide as an inhibitor of MDM2

An N-terminal helical region of the tumor suppressor p53 binds in a hydrophobic cleft of the oncoprotein MDM2. A retroinverso isomer of the natural N-terminal helical peptide was found to interact with MDM2 using the same hydrophobic residues, Phe, Trp, and Leu. We propose that the retroinverso d-peptide adopts a right-handed helical conformation to achieve functional mimicry of the p53 peptide.     

Relative strength of cation-pi vs salt-bridge interactions: the Gtalpha(340-350) peptide/rhodopsin system

Interactions between cationic and aromatic side chains of amino acid residues, the so-called cation-pi interaction, are thought to contribute to the overall stability of the folded structure of peptides and proteins. The transferred NOE NMR structure of the G(t)alpha(340-350) peptide bound to photoactivated rhodopsin (R*) geometrically suggests a cation-pi interaction stabilizing the structure between the epsilon-amine of Lys341 and the aromatic ring of the C-terminal residue, Phe350. This interaction has been explored by varying substituents on the phenyl ring to alter the electron density of the aromatic ring of Phe350 and observing the impact on binding of the peptide to R*. The results suggest that while a cation-pi interaction geometrically exists in the G(t)alpha(340-350) peptide when bound to R*, its energetic contribution to the stability of the receptor-bound structure is relatively insignificant, as it was not observed experimentally. The presence of an adjacent and competing salt-bridge interaction between the epsilon-amine of Lys341 and the C-terminal carboxylate of Phe350 effectively shields the charge of the ammonium group. Experimental data supporting a significant cation-pi interaction can be regained through a series of Phe350 analogues where the C-terminal carboxyl has been converted to the neutral carboxamide, thus eliminating the shielding salt-bridge. TrNOE NMR experiments confirmed the existence of the cation-pi interaction in the carboxamide analogues. Various literature estimates of the strength of cation-pi interactions, including some that estimate strengths in excess of salt-bridges, are compromised by omission of the relevant anion in the calculations.     

Cation-pi interaction in model alpha-helical peptides

Cation-pi interactions are increasingly recognized as important in chemistry and biology. Here we investigate the cation-pi interaction by determining its effect on the helicity of model peptides using a combination of CD and NMR spectroscopy. The data show that a single Trp/Arg interaction on the surface of a peptide can make a significant net favorable free energy contribution to helix stability if the two residues are positioned with appropriate spacing and orientation. The solvent-exposed Trp-->Arg (i, i + 4) interaction in helices can contribute -0.4 kcal/mol to the helix stability, while no free energy gain is detected if the two residues have the reversed orientation, Arg-->Trp (i, i + 4). The derived free energy is consistent with other experimental results studied in proteins or model peptides on cation-pi interactions. However in the same system the postulated Phe/Arg (i, i + 4) cation-pi interaction provides no net free energy to helix stability. Thus the Trp-->Arg interaction is stronger than Phe-->Arg. The cation-pi interactions are not sensitive to the screening effect by adding neutral salt as indicated by salt titration. Our results are in qualitative agreement with theoretical calculations emphasizing that cation-pi interactions can contribute significantly to protein stability with the order Trp > Phe. However, our and other experimental values are significantly smaller than estimates from theoretical calculations.     

Peptide/protein-polymer conjugates: synthetic strategies and design concepts

This feature article provides a compilation of tools available for preparing well-defined peptide/protein-polymer conjugates, which are defined as hybrid constructs combining (i) a defined number of peptide/protein segments with uniform chain lengths and defined monomer sequences (primary structure) with (ii) a defined number of synthetic polymer chains. The first section describes methods for post-translational, or direct, introduction of chemoselective handles onto natural or synthetic peptides/proteins. Addressed topics include the residue- and/or site-specific modification of peptides/proteins at Arg, Asp, Cys, Gln, Glu, Gly, His, Lys, Met, Phe, Ser, Thr, Trp, Tyr and Val residues and methods for producing peptides/proteins containing non-canonical amino acids by peptide synthesis and protein engineering. In the second section, methods for introducing chemoselective groups onto the side-chain or chain-end of synthetic polymers produced by radical, anionic, cationic, metathesis and ring-opening polymerization are described. The final section discusses convergent and divergent strategies for covalently assembling polymers and peptides/proteins. An overview of the use of chemoselective reactions such as Heck, Sonogashira and Suzuki coupling, Diels-Alder cycloaddition, Click chemistry, Staudinger ligation, Michael's addition, reductive alkylation and oxime/hydrazone chemistry for the convergent synthesis of peptide/protein-polymer conjugates is given. Divergent approaches for preparing peptide/protein-polymer conjugates which are discussed include peptide synthesis from synthetic polymer supports, polymerization from peptide/protein macroinitiators or chain transfer agents and the polymerization of peptide side-chain monomers.     

Interaction between N-terminal loop and beta-scaffold of photoactive yellow protein

During the photoreaction cycle of photoactive yellow protein (PYP), a physiologically active intermediate (PYP(M)) is formed as a consequence of global protein conformational change. Previous studies have demonstrated that the photocycle of PYP is regulated by the N-terminal loop region, which is located across the central beta-sheet from the p-coumaric acid chromophore. In this paper, the hydrophobic interaction between N-terminal loop and beta-sheet was studied by characterizing PYP mutants of the hydrophobic residues. The rate constants and structural changes of the photocycle of L15A and L23A possibly participating in such an interaction were more similar to wild-type than F6A, showing that the CH/pi interaction between Phe6 and Lys123 is the most essential as reported previously. To better understand the interactions between N-terminal tail and beta-sheet of PYP, Phe6 and Phe121 were replaced by Cys and linked by a disulfide bond. Since the photocycle kinetics, structural change and thermal stability of F6C/F121C were similar to F6A, the CH/pi interaction between Phe6 and Lys123 is not substitutable. It is likely that the detachment of position 6 from position 123 substantially alters the nature of PYP.     

Optimization of a biomimetic transamination reaction

For a range of protein substrates, N-terminal transamination offers a convenient way to install a reactive ketone or aldehyde functional group at a single location. We report herein the effects of the identity of N-terminal residues on the product distribution generated upon reaction with pyridoxal 5'-phosphate (PLP). This study was accomplished through the combination of solid-phase peptide synthesis with detailed liquid chromatography-mass spectrometry analysis. Many N-terminal amino acids provided high yields of the desired transaminated products, but some residues (His, Trp, Lys, and Pro) generated adducts with PLP itself. N-terminal Cys and Ser residues were observed to undergo beta-elimination in addition to transamination, and the transamination product of N-terminal Gln was resistant to subsequent oxime formation attempts. The information generated through the screening of peptide substrates was successfully applied to a protein target, changing an initially unreactive terminus into one that could be modified in high (70%) yield. Thus, these studies have increased our predictive power for the reaction, both in terms of improving conversion and suppressing reaction byproducts. An initial set of guidelines that may be used to increase the applicability of this reaction to specific proteins of interest is provided.     

Quantum mechanical studies of residue-specific hydrophobic interactions in p53-MDM2 binding

Quantum chemistry calculations at the levels of MP2/cc-pVDZ and MP2/cc-PVTZ have been carried out to study residue-specific interactions at the hydrophobic p53-MDM2 binding interface. The result of the calculation, based on structures from nanosecond molecular dynamics simulation, revealed that (19)Phe, (22)Leu, and (23)Trp of p53 have the strongest binding interaction with MDM2 followed by (26)Leu and (27)Pro. The specific residues of MDM2 that have dominant binding interactions with p53 are specifically identified to be (51)Lys, (54)Leu, (62)Met, (67)Tyr, (72)Gln, (94)Lys, (96)His, and (100)Tyr. The p53-MDM2 binding interaction is dominated by van der Waals interaction and to a lesser degree by electrostatic interaction. The MP2 results are in generally good agreement with those from the force field calculation while the DFT/B3LYP calculation failed to give attractive interaction energies for certain residue-residue interactions due to the lack of dispersion energy.     

Empirical formulation and parameterization of cation-π interactions for protein modeling

Cation-π interaction is comparable and as important as other main molecular interaction types, such as hydrogen bond, electrostatic interaction, van der Waals interaction, and hydrophobic interaction. Cation-π interactions frequently occur in protein structures, because six (Phe, Tyr, Trp, Arg, Lys, and His) of 20 natural amino acids and all metallic cations could be involved in cation-π interaction. Cation-π interactions arise from complex physicochemical nature and possess unique interaction behaviors, which cannot be modeled and evaluated by existing empirical equations and force field parameters that are widely used in the molecular dynamics. In this study, the authors present an empirical approach for cation-π interaction energy calculations in protein interactions. The accurate cation-π interaction energies of aromatic amino acids (Phe, Tyr, and Try) with protonated amino acids (Arg and Lys) and metallic cations (Li(+), Na(+), K(+), and Ca(2+)) are calculated using B3LYP/6-311+G(d,p) method as the benchmark for the empirical formulization and parameterization. Then, the empirical equations are built and the parameters are optimized based on the benchmark calculations. The cation-π interactions are distance and orientation dependent. Correspondingly, the empirical equations of cation-π interactions are functions of two variables, the distance r and the orientation angle θ. Two types of empirical equations of cation-π interactions are proposed. One is a modified distance and orientation dependent Lennard-Jones equation. The second is a polynomial function of two variables r and θ. The amino acid-based empirical equations and parameters provide simple and useful tools for evaluations of cation-π interaction energies in protein interactions.     

Mechanistic insights into mode of action of rice allene oxide synthase on hydroxyperoxides: An intermediate step in herbivory-induced jasmonate pathway

Various types of oxygenated fatty acids termed ‘oxylipins’ are involved in plant response to herbivory. Oxylipins like jasmonic acid (JA) and green leafy volatiles (GLVs) are formed by the action of enzymes like allene oxide synthase (AOS) and hydroxyperoxide lyase (HPL) respectively. In this study, we focus on AOS of Oryza sativa sb. Japonica, that interact with 9- and 13- hydroxyperoxides to produce intermediates of jasmonate pathway and compare it with rice HPL that yields GLVs. We attempt to elucidate the interaction pattern by computational docking protocols keeping the Arabidopsis AOS system as the reference model system. Both 9-hydroxyperoxide and 13-hydroxyperoxide fit into the active site of AOS completely with Phe347, Phe92, Ile463, Val345, and Asn278 being the common interacting residues. Phe347 and Phe92 were mutated with Leucine and docked again with the hydroxyperoxides. The Phe347 → Leu347 mutant showed a different mode of action than AOS-hydroxyperoxide complex with Trp413 in direct bonding with the OOH group of 9-hydroxyperoxide. The loss of Lys88-OOH interaction in 13-hydroxyperoxide and loss-of-interaction of Leu347 indicated the importance of Phe347 residue in hydroxyperoxide catalysis. The second mutant Phe92 → Leu92 also shows a very different interaction pattern with 13-hydroxyperoxide but not with 9-hydroxyperoxide.Therefore, it can be concluded that Phe347 is more crucial for AOS functionality than Phe92. The aromatic ring of a Phenylalanine residue is important for catalysis and its mutation affects the binding of the two ligands. Another important residue is Asn278 which is an important part of the AOS catalytic site for maintaining stability and can be compared with the Arabidopsis AOS residue Asn321. Lastly, the interaction of HPL with these two derivatives involves Leu363 residue instead of Phe347 and thus, validating the importance of Phe → Leu substitution to be the reason of different modes of action that result in completely different products from same substrates.