Conformational Preferences of Proline Oligopeptides

Proline, a unique amino acid whose polymer, polyproline II helix, and its analogues are involved in many biological processes: a review

Proline is a unique amino acid in that its side-chain is cyclised to the backbone, thus giving proline an exceptional rigidity and a considerably restricted conformational space. Polyproline forms two well-characterized helical structures: a left-handed polyproline helix (PPII) and a right-handed polyproline helix (PPI). Usually, sequences made only of prolyl residues are in PPII conformation, but even sequences not rich in proline but which are rich in glycine, lysine, glutamate, or aspartate have also a tendency to form PPII helices. Currently, the only way to study unambiguously PPII structure in solution is to use spectroscopies based on optical activity such as circular dichroism, vibrational circular dichroism and Raman optical activity. The importance of the PPII structure is emphasized by its ubiquitous presence in different organisms from yeast to human beings where proline-rich motifs and their binding domains are believed to be involved in vital biological processes. Some of the domains that are bound by proline-rich motifs include SH3 domains, WW domains, GYF domains and UEV domains, etc. The PPII structure has been demonstrated to be essential to biological activities such as signal transduction, transcription, cell motility, and immune response.

Deducing the conformational space for an octa-proline helix

A molecular dyad, PY-P8-PER, comprising a proline octamer sandwiched between pyrene and perylene terminals has been synthesized in order to address the dynamics of electronic energy transfer (EET) along the oligo-proline chain. A simple pyrene-based control compound equipped with a bis-proline attachment serves as a reference for spectroscopic studies. The N-H NMR signal at the terminal pyrene allows distinction between cis and trans amides and, although the crystal structure for the control has the trans conformation, temperature-dependent NMR studies provide clear evidence for trans/cis isomerisation in D6-DMSO. Polar solvents tend to stabilise the trans structure for the pyrene amide group, even for longer oligo-proline units. Circular dichroism shows that the proline spacer for PY-P8-PER exists mainly in the all-trans geometry in methanol. Preferential excitation of the pyrene chromophore is possible at wavelengths in the 320-350 nm range and, for the dyad, is followed by efficacious EET to the perylene emitter. The probability for intramolecular EET, obtained from analysis of steady-state spectroscopic data, is ca. 80-90% in solvents of disparate polarity. Comparison with the Förster critical distance suggests the terminals are ca. 18 Å apart. Time-resolved fluorescence spectroscopy, in conjunction with DFT calculations, indicates the dyad exists as a handful of conformers displaying a narrow range of EET rates. Optimisation of a distributive model allows accurate simulation of the EET dynamics in terms of reasonable structures based on isomerisation of certain amide groups.

A Unique and Stable Polyproline I Helix Sorted out from Conformational Equilibrium by Solvent Polarity

Polyproline I helical structures are often considered as the hidden face of their most famous geminal sibling, Polyproline II, as PPI is generally spotted only within a conformational equilibrium. We designed and synthesized a stable Polyproline I structure exploiting the striking tendency of (S)-indoline-2-carboxylic acid to drive the peptide bond conformation toward the cis amide isomer, when dissolved in polar solvents. The cooperative effect of only four amino acidic units is sufficient to form a preferential structure in solution. We shed light on this rare secondary structure with a thorough analysis of the spectroscopic and chiroptical properties of the tetramer, supported by X-ray crystallography and computational studies.

Pseudoprolines as stereoelectronically tunable proline isosteres

The cyclic structure of proline (Pro) confers unique conformational properties on this natural amino acid that influences polypeptide structure and function. Pseudoprolines are a family of Pro isosteres that incorporate a heteroatom, most prominently oxygen or sulfur but also silicon and selenium, to replace the C_β or C_γ carbon atom of the pyrrolidine ring. These readily synthetically accessible structural motifs can facilitate facile molecular editing in a fashion that allows modulation of the amide bond topology of dipeptide elements and influence over ring pucker. While the properties of pseudoprolines have been exploited most prominently in the design of oligopeptide analogues, they have potential application in the design and optimization of small molecules. In this Digest, we summarize the physicochemical properties of pseudoprolines and illustrate their potential in drug discovery by surveying examples of applications in the design of bioactive molecules.

The Alanine World Model for the Development of the Amino Acid Repertoire in Protein Biosynthesis

A central question in the evolution of the modern translation machinery is the origin and chemical ethology of the amino acids prescribed by the genetic code. The RNA World hypothesis postulates that templated protein synthesis has emerged in the transition from RNA to the Protein World. The sequence of these events and principles behind the acquisition of amino acids to this process remain elusive. Here we describe a model for this process by following the scheme previously proposed by Hartman and Smith, which suggests gradual expansion of the coding space as GC–GCA–GCAU genetic code. We point out a correlation of this scheme with the hierarchy of the protein folding. The model follows the sequence of steps in the process of the amino acid recruitment and fits well with the co-evolution and coenzyme handle theories. While the starting set (GC-phase) was responsible for the nucleotide biosynthesis processes, in the second phase alanine-based amino acids (GCA-phase) were recruited from the core metabolism, thereby providing a standard secondary structure, the α-helix. In the final phase (GCAU-phase), the amino acids were appended to the already existing architecture, enabling tertiary fold and membrane interactions. The whole scheme indicates strongly that the choice for the alanine core was done at the GCA-phase, while glycine and proline remained rudiments from the GC-phase. We suggest that the Protein World should rather be considered the Alanine World, as it predominantly relies on the alanine as the core chemical scaffold.

Which DFT levels of theory are appropriate in predicting the prolyl cis–trans isomerization in solution?

DFTs were assessed for the conformational preferences of the peptides containing Pro and its derivatives in chloroform and water.

Exploring hydrophobicity limits of polyproline helix with oligomeric octahydroindole-2-carboxylic acid

The polyproline-II helix is the most extended naturally occurring helical structure and is widely present in polar, exposed stretches and "unstructured" denatured regions of polypeptides. Can it be hydrophobic? In this study, we address this question using oligomeric peptides formed by a hydrophobic proline analogue, (2S,3aS,7aS)-octahydroindole-2-carboxylic acid (Oic). Previously, we found the molecular principles underlying the structural stability of the polyproline-II conformation in these oligomers, whereas the hydrophobicity of the peptide constructs remains to be examined. Therefore, we investigated the octan-1-ol/water partitioning and inclusion in detergent micelles of the oligo-Oic peptides. The results showed that the hydrophobicity is remarkably enhanced in longer oligomeric sequences, and the oligo-Oic peptides with 3 to 4 residues and higher are specific towards hydrophobic environments. This contrasts significantly to the parent oligoproline peptides, which were moderately hydrophilic. With these findings, we have demonstrated that the polyproline-II structure is compatible with nonpolar media, whereas additional manipulations of the terminal functionalities feature solubility in extremely nonpolar solvents such as hexane.

Response speed control of helicity inversion based on a "regulatory enzyme"-like strategy

In biological systems, there are many signal transduction cascades in which a chemical signal is transferred as a series of chemical events. Such successive reaction systems are advantageous because the efficiency of the functions can be finely controlled by regulatory enzymes at an earlier stage. However, most of artificial responsive molecules developed so far rely on single-step conversion, whose response speeds have been difficult to be controlled by external stimuli. In this context, developing artificial conversion systems that have a regulation step similar to the regulatory enzymes has been anticipated. Here we report a novel artificial two-step structural conversion system in which the response speed can be controlled based on a regulatory enzyme-like strategy. In this system, addition of fluoride ion caused desilylation of the siloxycarboxylate ion attached to a helical complex, resulting in the subsequent helicity inversion. The response speeds of the helicity inversion depended on the reactivity of the siloxycarboxylate ions; when a less-reactive siloxycarboxylate ion was used, the helicity inversion rate was governed by the desilylation rate. This is the first artificial responsive molecule in which the overall response speed can be controlled at the regulation step separated from the function step.

Design, synthesis and biological evaluation of novel peptide MC2 analogues from Momordica charantia as potential anti-diabetic agents

Three series of Momordica charantia (MC)2 analogues were designed, synthesized and evaluated for their anti-hyperglycaemic effects. Alanine scanning focusing on the peptide MC2 indicated the importance of the residues proline (Pro)(3), serine (Ser)(6), isoleucine (Ile)(7) and Ser(10) for anti-hyperglycaemic effects. Among the first series of MC2 analogues, peptide I-4 exhibited a better anti-hyperglycaemic effect and was chosen for further modification. A further two series of conformationally constrained analogues were designed by scanning the residues Pro(3), Ser(6), Ile(7), and Ser(10) with an i - (i + 2) lactam bridge consisting of a glutamic acid-Xaa-lysine (Glu-Xaa-Lys) scaffold and a diproline fragment. By screening in normal mice and mice with diabetes mellitus, peptides II-1, II-2 and III-3 showed a significant improvement in anti-hyperglycaemic and anti-oxidative activities compared with I-4. These data suggest that II-1, II-2 and III-3 could be candidates for future treatment of diabetes mellitus.

Synthesis and helical structures of poly(ω-alkynamide)s having chiral side chains: effect of solvent on their screw-sense inversion

New ω-alkynamides, (S)-HC≡CCH2CONHCH2CH(CH3)CH2CH3 (1) and (S)-HC≡CCH2CH2CONHCH(CH3)CH2CH2CH2CH2CH3 (2) were synthesized and polymerized with a rhodium catalyst in CHCl3 to obtain cis-stereoregular poly(ω-alkynamide)s (poly(1) and poly(2)). Polarimetric, CD, and IR spectroscopic studies revealed that in solution the polymers adopted predominantly one-handed helical structures stabilized by intramolecular hydrogen bonds between the pendent amide groups. This behavior was similar to that of the corresponding poly(N-alkynylamide) counterparts (poly(3) and poly(4)) reported previously, whereas the helical senses were opposite to each other. The helical structures of the poly(ω-alkynamide)s were stable upon heating similar to those of the poly(N-alkynylamide)s, but the solvent response was completely different. An increase in MeOH content in CHCl3/MeOH resulted in inversion of the predominant screw-sense for poly(1) and poly(2). Conversely, poly(3) was transformed into a random coil, and poly(4) maintained the predominant screw-sense irrespective of MeOH content. The solvent dependence of predominant screw-sense for poly(1) and poly(2) was reasonably explained by molecular orbital studies using the conductor-like screening model (COSMO).

Molecular dynamics of the proline switch and its role in Crk signaling

The Crk adaptor proteins play a central role as a molecular timer for the formation of protein complexes including various growth and differentiation factors. The loss of regulation of Crk results in many kinds of cancers. A self-regulatory mechanism for Crk was recently proposed, which involves domain–domain rearrangement. It is initiated by a cis–trans isomerization of a specific proline residue (Pro238 in chicken Crk II) and can be accelerated by Cyclophilin A. To understand how the proline switch controls the autoinhibition at the molecular level, we performed large-scale molecular dynamics and metadynamics simulations in the context of short peptides and multidomain constructs of chicken Crk II. We found that the equilibrium and kinetic properties of the macrostates are regulated not only by the local environments of specified prolines but also by the global organization of multiple domains. We observe the two macrostates (cis closed/autoinhibited and trans open/uninhibited) consistent with NMR experiments and predict barriers. We also propose an intermediate state, the trans closed state, which interestingly was reported to be a prevalent state in human Crk II. The existence of this macrostate suggests that the rate of switching off the autoinhibition by Cyp A may be limited by the relaxation rate of this intermediate state.

Solvent interactions stabilize the polyproline II conformation of glycosylated oligoprolines

In nature, proline residues carry several post-translational modifications (PTMs), including 4R hydroxylation and glycosylation. A recent study synthesized contiguously hydroxylated and glycosylated nonaproline peptides and revealed that both PTMs lead to a significant increase in the thermal stability of PPII relative to the unmodified oligoproline. The increased stability of the hydroxylated peptide can be explained by increased stability of the trans isomer due to stereoelectronic effects. However, the effects of glycosylation cannot be completely explained by stereoelectronics since previous experimental results indicate that 4R-glycosylation does not produce observable changes in the trans preference compared to 4R-hydroxylation. We therefore used sophisticated molecular modeling techniques to determine the reason for the further increase in thermal stability upon glycosylation. Free energy estimates obtained from adaptively biased molecular dynamics calculations in implicit (explicit) solvent are -9 kcal mol(-1) (-20 kcal mol(-1)) for the hydroxylated compound and -9 kcal mol(-1) (-46 kcal mol(-1)) for the glycosylated compound, indicating that direct solvent-peptide interactions are vital for explaining the glycosylation effects on PPII stability. Our data reveals for the first time that interactions between the hydroxyl groups in the glycosylated compound and water act in a complementary fashion with stereoelectronic effects to stabilize the PPII conformation in these substituted oligoproline peptides.

Recipes for free energy calculations in biomolecular systems

During the last decade, several methods for sampling phase space and calculating various free energies in biomolecular systems have been devised or refined for molecular dynamics (MD) simulations. Thus, state-of-the-art methodology and the ever increasing computer power allow calculations that were forbidden a decade ago. These calculations, however, are not trivial as they require knowledge of the methods, insight into the system under study, and, quite often, an artful combination of different methodologies in order to avoid the various traps inherent in an unknown free energy landscape. In this chapter, we illustrate some of these concepts with two relatively simple systems, a sugar ring and proline oligopeptides, whose free energy landscapes still offer considerable challenges. In order to explore the configurational space of these systems, and to surmount the various free energy barriers, we combine three complementary methods: a nonequilibrium umbrella sampling method (adaptively biased MD, or ABMD), replica-exchange molecular dynamics (REMD), and steered molecular dynamics (SMD). In particular, ABMD is used to compute the free energy surface of a set of collective variables; REMD is used to improve the performance of ABMD, to carry out sampling in space complementary to the collective variables, and to sample equilibrium configurations directly; and SMD is used to study different transition mechanisms.

Parameterization of the proline analogue Aze (azetidine-2-carboxylic acid) for molecular dynamics simulations and evaluation of its effect on homo-pentapeptide conformations

We have parameterized and evaluated the proline homologue Aze (azetidine-2-carboxylic acid) for the gromos56a3 force-field for use in molecular dynamics simulations using GROMACS. Using bi-phasic cyclohexane/water simulation systems and homo-pentapeptides, we measured the Aze solute interaction potential energies, ability to hydrogen bond with water, and overall compaction, for comparison to Pro, Gly, and Lys. Compared to Pro, Aze has a slightly higher H-bonding potential, and stronger electrostatic but weaker non-electrostatic interactions with water. The 20-ns simulations revealed the preferential positioning of Aze and Pro at the interface of the water and cyclohexane layers, with Aze spending more time in the aqueous layer. We also demonstrated through simulations of the homo-pentapeptides that Aze has a greater propensity than Pro to undergo trans→cis peptide bond isomerization, which results in a severe 180° bend in the polypeptide chain. The results provide evidence for the hypothesis that the misincorporation of Aze within proline-rich regions of proteins could disrupt the formation of poly-proline type II structures and compromise events such as recognition and binding by SH3-domains.

Conformational preference and cis-trans isomerization of 4-methylproline residues

Conformational preferences and prolyl cis-trans isomerizations of the (2S,4S)-4-methylproline (4S-MePro) and (2S,4R)-4-methylproline (4R-MePro) residues are explored at the M06-2X/cc-pVTZ//M06-2X/6-31+G(d) level of theory in the gas phase and in water, where solvation free energies were calculated using the implicit SMD model. In the gas phase, the down-puckered γ-turn structure with the trans prolyl peptide bond is most preferred for both Ac-4S-MePro-NHMe and Ac-4R-MePro-NHMe, in which the C(7) hydrogen bond between two terminal groups seems to play a role, as found for Ac-Pro-NHMe. Because of the C(7) hydrogen bonds weakened by the favorable direct interactions between the backbone C==O and H--N groups and water molecules, the 4S-MePro residue has a strong preference of the up-puckered polyproline II (PP(II)) structure over the down-puckered PP(II) structure in water, whereas the latter somewhat prevails over the former for the 4R-MePro residue. However, these two structures are nearly equally populated for Ac-Pro-NHMe. The calculated populations for the backbone structures of Ac-4S-MePro-NHMe and Ac-4R-MePro-NHMe in water are reasonably consistent with CD and NMR experiments. In particular, our calculated results on the puckering preference of the 4S-MePro and 4R-MePro residues with the PP(II) structures are in accord with the observed results for the stability of the (X-Y-Gly)(7) triple helix with X = 4R-MePro or Pro and Y = 4S-MePro or Pro. The calculated rotational barriers indicate that the cis-trans isomerization may in common proceed through the anticlockwise rotation for Ac-4S-MePro-NHMe, Ac-4R-MePro-NHMe, and Ac-Pro-NHMe in water. The lowest rotational barriers become higher by 0.24-1.43 kcal/mol for Ac-4S-MePro-NHMe and Ac-4R-MePro-NHMe than those for Ac-Pro-NHMe in water.

A preliminary survey of the peptoid folding landscape

We present an analysis of the conformational preferences of N-substituted glycine peptoid oligomers. We survey the backbone conformations observed in experimentally determined peptoid structures and provide a comparison with high-level quantum mechanics calculations of short peptoid oligomers. The dominant sources of structural variation derive from: side-chain dependent cis/trans isomerization of backbone amide bonds, side chain stereochemistry, and flexibility in the psi dihedral angle. We find good agreement between the clustering of experimentally determined peptoid torsion angles and local torsional minima predicted by theory for a disarcosine model. The calculations describe a well-defined conformational map featuring distinct energy minima. The general features of the peptoid backbone conformational landscape are consistent across a range of N-alkyl glycine side chains. Alteration of side chain types, however, creates subtle but potentially significant variations in local folding propensities. We identify a limited number of low energy local conformations, which may be preferentially favored by incorporation of particular monomer units. Greater variation in backbone dihedral angles are accessible in peptoids featuring trans amide bond geometries. These results confirm that computational approaches can play a valuable role in guiding the design of complex peptoid architectures and may lead to strategies for introducing constraints that select among a limited number of low energy local conformations.

Conformational preferences and cis-trans isomerization of L-3,4-dehydroproline residue

The conformational study of N-acetyl-N'-methylamide of L-3,4-dehydroproline (Ac-Dhp-NHMe, the Dhp dipeptide) is carried out using hybrid density functional methods with the self-consistent reaction field method in the gas phase and in solution (chloroform and water). The incorporation of a double bond between C(beta) and C(gamma) into the prolyl ring results in the puckering, backbone population, and barriers to prolyl cis-trans isomerization different from those of the Pro dipeptide. For local minima of the Dhp dipeptide in the gas phase and in water, the C(beta)-C(gamma) bonds become shorter by approximately 0.2 A and the bond angles C(alpha)-C(beta)-C(gamma) and C(beta)-C(gamma)-C(delta) are widened by approximately 8 degrees than those of the Pro dipeptide, and the puckering amplitude is computed to be 0.01-0.07 A, indicating that the 3,4-dehydroprolyl ring is quite less puckered. However, polyproline-like conformations become preferred and the relative stability of the conformation tC with a C(7) intramolecular hydrogen bond decreases as the solvent polarity increases, as found for the Pro dipeptide. The barriers to cis-trans isomerization of the Ac-Dhp peptide bond increase with the increase of solvent polarity and the isomerization is likely to proceed through the clockwise rotation in water, as found for the prolyl peptide bond. The hydrogen bond between the prolyl nitrogen and the following amide N-H group seems to contribute in stabilizing the transition state structures.

Conformations and free energy landscapes of polyproline peptides

The structure of the proline amino acid allows folded polyproline peptides to exist as both left- (PPII) and right-handed (PPI) helices. We have characterized the free energy landscapes of hexamer, nanomer, and tridecamer polyproline peptides in gas phase and implicit water as well as explicit hexane and 1-propanol for the nanomer. To enhance the sampling provided by regular molecular dynamics, we used the recently developed adaptively biased molecular dynamics method, which describes Landau free energy maps in terms of relevant collective variables. These maps, as a function of the collective variables of handedness, radius of gyration, and three others based on the peptide torsion angle omega, were used to determine the relative stability of the different structures, along with an estimate of the transition pathways connecting the different minima. Results show the existence of several metastable isomers and therefore provide a complementary view to experimental conclusions based on photo-induced electron transfer experiments with regard to the existence of stable heterogeneous subpopulations in PPII polyproline.

Conformational preferences of beta- and gamma-aminated proline analogues

Quantum mechanical calculations have been used to investigate how the incorporation of an amino group to the Cbeta- or Cgamma-positions of the pyrrolidine ring affects the intrinsic conformational properties of the proline. Specifically, a conformational study of the N-acetyl-N'-methylamide derivatives of four isomers of aminoproline, which differ not only in the beta- or gamma-position of the substituent but also in its cis or trans relative disposition, has been performed. To further understand the role of the intramolecular hydrogen bonds between the backbone carbonyl groups and the amino side group, a conformational study was also performed on the corresponding four analogues of (dimethylamino)proline. In addition, the effects of solvation on aminoproline and (dimethylamino)proline dipeptides have been evaluated using a self-consistent reaction field model, and considering four different solvents (carbon tetrachloride, chloroform, methanol and water). Results indicate that the incorporation of the amino substituent into the pyrrolidine ring affects the conformational properties, with backbone...side chain intramolecular hydrogen bonds detected when it is incorporated in a cis relative disposition. In general, the incorporation of the amino side group tends to stabilize those structures where the peptide bond involving the pyrrolidine nitrogen is arranged in cis. The aminoproline isomer with the substituent attached to the Cgamma-position with a cis relative disposition is the most stable in the gas phase and in chloroform, methanol and water solutions. Replacement of the amino side group by the dimethylamino substituent produces significant changes in the potential energy surfaces of the four investigated (dimethylamino)proline-containing dipeptides. Thus, these changes affect not only the number of minima, which increases considerably, but also the backbone and pseudorotational preferences. In spite of these effects, comparison of the conformational preferences, i.e., the more favored conformers, calculated for different isomers of aminoproline and (dimethylamino)proline dipeptides showed a high degree of consistency for the two families of compounds.

Conformations of proline analogues having double bonds in the ring

The intrinsic conformational preferences of proline analogues having double bonds between carbon atoms in their rings have been investigated using quantum mechanical calculations at the B3LYP/6-31+G(d,p) level. For this purpose, the potential energy surface of the N-acety-N'-methylamide derivatives of three dehydroprolines (proline analogues unsaturated at alpha,beta; beta,gamma; and gamma,delta) and pyrrole (proline analogue with unsaturations at both alpha,beta and gamma,delta) have been explored, and the results are compared with those obtained for the derivative of the nonmodified proline. We found that the double bonds affect the ring puckering and the geometric internal parameters, even though the backbone conformation was influenced the most. Results indicate that the formation of double bonds between carbon atoms in the pyrrolidine ring should be considered as an effective procedure to restrict the conformational flexibility of prolines. Interestingly, we also found that the N-acetyl-N'-methylamide derivative of pyrrole shows a high probability of having a cis peptide bond preceding the proline analogue.

Conformational Preferences and cis−trans Isomerization of Azaproline Residue

The conformational study of N-acetyl-N'-methylamide of azaproline (Ac-azPro-NHMe, the azPro dipeptide) is carried out using ab initio HF and density functional methods with the self-consistent reaction field method to explore the effects of the replacement of the backbone CHalpha group by the nitrogen atom on the conformational preferences and prolyl cis-trans isomerization in the gas phase and in solution (chloroform and water). The incorporation of the Nalpha atom into the prolyl ring results in the different puckering, backbone population, and barriers to prolyl cis-trans isomerization from those of Ac-Pro-NHMe (the Pro dipeptide). In particular, the azPro dipeptide has a dominant backbone conformation D (beta2) with the cis peptide bond preceding the azPro residue in both the gas phase and solution. This may be ascribed to the favorable electrostatic interaction or intramolecular hydrogen bond between the prolyl nitrogen and the amide hydrogen following the azPro residue and to the absence of the unfavorable interactions between electron lone pairs of the acetyl carbonyl oxygen and the prolyl Nalpha. This calculated higher population of the cis peptide bond is consistent with the results from X-ray and NMR experiments. As the solvent polarity increases, the conformations B and B* with the trans peptide bond become more populated and the cis population decreases more, which is opposite to the results for the Pro dipeptide. The conformation B lies between conformations D and A (alpha) and conformation B* is a mirror image of the conformation B on the phi-psi map. The barriers to prolyl cis-trans isomerization for the azPro dipeptide increase with the increase of solvent polarity, and the cis-trans isomerization proceeds through only the clockwise rotation with omega' approximately +120 degrees about the prolyl peptide bond for the azPro dipeptide in the gas phase and in solution, as seen for the Pro dipeptide. The pertinent distance d(N...H-NNHMe) and the pyramidality of imide nitrogen can describe the role of this hydrogen bond in stabilizing the transition state structure and the lower rotational barriers for the azPro dipeptide than those for the Pro dipeptide in the gas phase and in solution.