Study of the conformational profile of the cyclohexane analogs of L-phenylalanine.

The conformational profile of the conformationally constrained cyclohexane analogs of phenylalanine (1-amino-2-phenylcyclohexanecarboxylic acids, c6Phe) was assessed using computational methods. For this purpose, the conformational space of the N-acetyl methylamide derivatives of the stereoisomers (2S,3R)c6Phe and (2S,3S)c6Phe was explored by computing their respective Ramachandran maps, and low-energy minima were characterized at molecular mechanics level by means of the AMBER program, using the parm94 force field set of parameters. In order to assess the performance of the molecular mechanics calculations, each of the low-energy conformations was also investigated further at the ab initio level. Accordingly, the molecular mechanics geometries were used as starting conformations to perform full geometry optimizations at the Hartree-Fock level, using a 6-31G(d) basis set. Analysis of the results revealed that the cyclohexane structure directly induces some restrictions on the backbone, and constrains the orientation of the aromatic side-chain to two narrow regions for each stereoisomer. The conformational profile of these amino acids is then explained on the grounds of the interaction between the rigidly held phenyl ring and the main chain NH and CO groups. The results obtained are in good accordance with the experimental observations.     

A molecular dynamics study of the pentacyclo-undecane (PCU) cage polypeptides of the type Ac-3Ala-Cage-3Ala-NHMe

An extensive exploration of the conformational space of the seven-residue peptide sequences, Ac-Ala-Ala-Ala-Cage-Ala-Ala-Ala-NHMe and the model peptide Ac-Ala-Ala-Ala-Ala-Ala-Ala-Ala-NHMe, was carried out using single trajectories of molecular dynamics (MD) in the solution phase using the periodic boundary conditions. Our MD studies revealed that the majority of the motifs of the PCU cage peptide exist as type I–III β-turns along with their mirror conformations, viz. type I′–III′ β-turns. This peptide sequence adopted a U-shaped backbone, with alpha-helical characteristics. The results reported here provide further evidence that the PCU cage amino acid exhibits C_7eq, C_7aq, α_R and α_L conformations in aqueous solution.     

Crystal-state conformation of Calpha,alpha-dialkylated peptides containing chiral beta-homo-residues.

Secondary structure formation and stability are essential features in the knowledge of complex folding topology of biomolecules. To better understand the relationships between preferred conformations and functional properties of beta-homo-amino acids, the synthesis and conformational characterization by X-ray diffraction analysis of peptides containing conformationally constrained Calpha,alpha-dialkylated amino acid residues, such as alpha-aminoisobutyric acid or 1-aminocyclohexane-1-carboxylic acid and a single beta-homoamino acid, differently displaced along the peptide sequence have been carried out. The peptides investigated are: Boc-betaHLeu-(Ac6c)2-OMe, Boc-Ac6c-betaHLeu-(Ac6c)2-OMe and Boc-betaHVal-(Aib)5-OtBu, together with the C-protected beta-homo-residue HCl.H-betaHVal-OMe. The results indicate that the insertion of a betaH-residue at position 1 or 2 of peptides containing strong helix-inducing, bulky Calpha,alpha-disubstituted amino acid residues does not induce any specific conformational preferences. In the crystal state, most of the NH groups of beta-homo residues of tri- and tetrapeptides are not involved in intramolecular hydrogen bonds, thus failing to achieve helical structures similar to those of peptides exclusively constituted of Calpha,alpha-disubstituted amino acid residues. However, by repeating the structural motifs observed in the molecules investigated, a beta-pleated sheet secondary structure, and a new helical structure, named (14/15)-helix, were generated, corresponding to calculated minimum-energy conformations. Our findings, as well as literature data, strongly indicate that conformations of betaH-residues, with the micro torsion angle equal to -60 degrees, are very unlikely.     

Azetidine-2-carboxylic acid derivatives from seeds of Fagus silvatica L. and a revised structure for nicotianamine

2(S),3′(S)-N-(3-Amino-3-carboxypropyl)azetidine-2-carboxylic acid and 2(S),3′(S),3″(S)-N-[N-(3-amino-3-carboxypropyl)-3-amino-3-carboxypropyl]azetidine-2-carboxylic acid have been isolated from seeds of Fagus silvatica L. (beechnuts). The structures have been established by PMR- and ¹³C-NMR-spectroscopy and by synthesis from l-azetidine-2-carboxylic acid. The second of the new amino acids is identical with nicotianamine. previously isolated from Nicotiana tabacum but assigned a different formula. The ring opening reactions of azetidine-2-carboxylic acid in neutral solution have been studied and the chemical and possibly biochemical precursor role of this amino acid for various amino acids including the two new ones described here, nicotianine [N-(3-amino-3-carboxypropyl)nicotinic acid] and methionine is discussed.     

Folding of peptide fragments comprising the complete sequence of proteins. Models for initiation of protein folding. II. Plastocyanin.

In an attempt to understand the earliest events in the protein folding pathway, the complete sequence of French bean plastocyanin has been synthesized as a series of short peptide fragments, and the conformational preferences of each peptide examined in aqueous solution using proton n.m.r. methods. Plastocyanin consists largely of beta-sheet, with reverse turns and loops between the strands of the sheet, and one short helix. The n.m.r. experiments indicate that most of the peptides derived from the plastocyanin sequence have remarkably little propensity to adopt folded conformations in aqueous solution, in marked contrast to the peptides derived from the helical protein, myohemerythrin (accompanying paper). For most plastocyanin peptides, the backbone dihedral angles are predominantly in the beta-region of conformational space. Some of the peptides show weak NOE connectivities between adjacent amide protons, indicative of small local populations of backbone conformations in the a region of (phi,psi) space. A conformational preference for a reverse turn is seen in the sequence Ala65-Pro-Gly-Glu68, where a turn structure is found in the folded protein. Significantly, the peptide sequences that populate the alpha-region of (phi,psi) space are mostly derived from turn and loop regions in the protein. The addition of trifluoroethanol does not drive the peptides into helical conformations. In one region of the sequence, the n.m.r. spectra provide evidence of the formation of a hydrophobic cluster involving aromatic and aliphatic side-chains. These results have significance for understanding the initiation of protein folding. From these studies of the fragments of plastocyanin (this paper) and myohemerythrin (accompanying paper), it appears that there is a pre-partitioning of the conformational space sampled by the polypeptide backbone that is related to the secondary structure in the final folded state.     

Conformational basis for the biological activity of TOAC-labeled angiotensin II and bradykinin: electron paramagnetic resonance, circular dichroism, and fluorescence studies

N-Terminally and internally labeled analogues of the hormones angiotensin (AII, DRVYIHPF) and bradykinin (BK, RPPGFSPFR) were synthesized containing the paramagnetic amino acid 2,2,6,6-tetramethylpiperidine-1-oxyl-4-amino-4-carboxylic acid (TOAC). TOAC replaced Asp1 (TOAC1-AII) and Val3 (TOAC3-AII) in AII and was inserted prior to Arg1 (TOAC0-BK) and replacing Pro3 (TOAC3-BK) in BK. The peptide conformational properties were examined as a function of trifluoroethanol (TFE) content and pH. Electron paramagnetic resonance spectra were sensitive to both variables and showed that internally labeled analogues yielded rotational correlation times (tauC) considerably larger than N-terminally labeled ones, evincing the greater freedom of motion of the N-terminus. In TFE, tauC increased due to viscosity effects. Calculation of tau(Cpeptide)/tau(CTOAC) ratios indicated that the peptides acquired more folded conformations. Circular dichroism spectra showed that, except for TOAC1-AII in TFE, the N-terminally labeled analogues displayed a conformational behavior similar to that of the parent peptides. In contrast, under all conditions, the TOAC3 derivatives acquired more restricted conformations. Fluorescence spectra of AII and its derivatives were especially sensitive to the ionization of Tyr4. Fluorescence quenching by the nitroxide moiety was much more pronounced for TOAC3-AII. The conformational behavior of the TOAC derivatives bears excellent correlation with their biological activity, since, while the N-terminally labeled peptides were partially active, their internally labeled counterparts were inactive [Nakaie, C. R., et al., Peptides 2002, 23, 65-70]. The data demonstrate that insertion of TOAC in the middle of the peptide chain induces conformational restrictions that lead to loss of backbone flexibility, not allowing the peptides to acquire their receptor-bound conformation.     

Computational study of the conformational preferences of the (R)-8-amino-pentacyclo[5.4.0.0(2,6).0(3,10).0(5,9)]undecane-8-carboxylic acid monopeptide.

alpha-Amino acids are important building blocks for the synthesis of a large number of bioactive compounds and pharmaceutical drugs. However, a literature survey revealed that no theoretical conformational study of alpha-amino acids with cage carbon frameworks has been performed to date. This paper reports the results of a conformational study on the (R)-8-amino-pentacyclo[5.4.0.0(2,6).0(3,10).0(5,9)]undecane-8-carboxylic acid monopeptide (cage monopeptide), using molecular mechanics and ab initio methods. The in vacuo Ramachandran maps computed using the different parameterizations of the AMBER force field show the C7eq structure as the most favourable conformation, in contrast to the C7ax structure, that is the lowest energy conformation at the ab initio level. Analysis of these maps reveals the helical preference for the monopeptide and provides the potential for the cage residue to be incorporated into constrained peptide analogues.     

Characterization of the conformational domains of bradykinin by computational methods

The AMBER 4.0 force field was used to perform a characterization of the conformational profile of the nonapeptide bradykinin. A thorough conformational search was carried out using molecular dynamics as sampling technique, by computing cycles of high (900 K) and low (300 K) temperature trajectories. A total of 2400 minima were generated and subsequently clustered using the root-mean-square of the backbone dihedral angles as criterium. After the use of a tolerance value of 20deg;, the conformations were clustered in 233 unique conformations with energies up to 40 kcal/mol above the lowest minimum. The analysis of the low-energy conformations indicate that the peptide exhibits a high tendency to adopt a β-turn at the C-terminus and a propensity to adopt a bent structure at the N-terminus. These results are in agreement with the experimental evidence reported in the literature and provide detailed information necessary to understand the conformational preferences of the peptide.     

Classification of conformation for sugar amino acid systems using chiroptical spectroscopy

A number of structurally related sugar amino acid systems have been examined by chiroptical methods to aid interpretation of their conformational preference. The use of circular dichroism, in addition to NMR and solution IR, has enabled classification of the conformations adopted by sugar amino acid systems as hydrogen-bonded regular, non-hydrogen-bonded regular, and non-hydrogen-bonded irregular. A set of tetrameric SAAs are examined and the effect of change in primary structure related to conformation.     

Reduced peptide bond cyclic somatostatin based opioid octapeptides. Synthesis, conformational properties and pharmacological characterization.

The conformational and pharmacological properties that result from peptide bond reduction as well as the use of secondary amino acids in a series of cyclic peptides related to the mu opioid receptor selective antagonist D-Phe1-Cys2-Tyr3-D-Trp4-Orn5-Thr6-Pen7+ ++-Thr8-NH2 (IV), have been investigated. Peptide analogues that contain [CH2NH] and [CH2N] pseudo-peptide bonds (in primary and secondary amino acids, respectively) were synthesized on a solid support. Substitution of Tyr3 in IV by the cyclic, secondary amino acid 1,2,3,4-tetrahydroisoquinoline carboxylate (Tic) and of D-Trp4 with D-1,2,3,4-tetrahydro-beta-carboline(D-Tca4), gave peptides 4 and 1, respectively. Both analogues displayed reduced affinities for mu opioid receptors. Conformational analysis based on extensive NMR investigations demonstrated that the backbone conformations of 1 and 4 are similar to those of the potent and selective analogue D-Phe-Cys-Tyr-D-Trp-Lys-Thr-Pen-Thr-NH2 (I), while the conformational properties of the side chains of Tic3 (4) and D-Tca4 (1) resulted in topographical properties that were not well recognized by the mu opioid receptor. Peptide bond modifications were made including (Tyr3-psi[CH2NH]-D-Trp4), 3; (Tyr3-psi[CH2N]-D-Tca4), 2; and (Cys2-psi[CH2N]-Tic3), 6. These analogues showed decreases in their mu opioid receptor affinities relative to the parent compounds IV, 1, and 4, respectively. 1H NMR based conformational analysis in conjunction with receptor binding data led to the conclusion that the reduced peptide bonds in 2, 3, 5, and 6 do not contribute to the process of discrimination between mu and delta opioid receptors, and in spite of their different dynamic behaviors (relative to 1 and 4), they are still capable of attaining similar receptor bound conformations, possibly due to their increased flexibility.     

Folding of polypeptide chains induced by the amino acid side-chains

Conformational calculations with the use of semi-empirical potential functions have been applied to the analysis of the folding of peptide chains. In particular, the part played by the amino acid side-chains in the adoption of folded conformations has been investigated.The results show that the preferred conformations of short peptides are mostly extended ones. However, from a given peptide chain-length, the side-chain to backbone and side-chain to side-chain interactions become strong enough so that definite sequences of amino acids can induce a transition from extended to folded conformations. We propose to call these folded structures “conformational nuclei”. The type of “nucleus” formed is dependent on both the amino acid composition and the sequence.Our results strongly support the hypothesis that folding of polypeptide chains can occur through a nucleation process that could be induced by the side-chains.     

Backbone conformational flexibility of the lipid modified membrane anchor of the human N-Ras protein investigated by solid-state NMR and molecular dynamics simulation

The lipid modified human N-Ras protein, implicated in human cancer development, is of particular interest due to its membrane anchor that determines the activity and subcellular location of the protein. Previous solid-state NMR investigations indicated that this membrane anchor is highly dynamic, which may be indicative of backbone conformational flexibility. This article aims to address if a dynamic exchange between three structural models exist that had been determined previously. We applied a combination of solid-state nuclear magnetic resonance (NMR) methods and replica exchange molecular dynamics (MD) simulations using a Ras peptide that represents the terminal seven amino acids of the human N-Ras protein. Analysis of correlations between the conformations of individual amino acids revealed that Cys 181 and Met 182 undergo collective conformational exchange. Two major structures constituting about 60% of all conformations could be identified. The two conformations found in the simulation are in rapid exchange, which gives rise to low backbone order parameters and nuclear spin relaxation as measured by experimental NMR methods. These parameters were also determined from two 300 ns conventional MD simulations, providing very good agreement with the experimental data.     

Folding of peptide fragments comprising the complete sequence of proteins. Models for initiation of protein folding. I. Myohemerythrin.

In an attempt to delineate potential folding initiation sites for different protein structural motifs, we have synthesized series of peptides that span the entire length of the polypeptide chain of two proteins, and examined their conformational preferences in aqueous solution using proton nuclear magnetic resonance and circular dichroism spectroscopy. We describe here the behavior of peptides derived from a simple four-helix bundle protein, myohemerythrin. The peptides correspond to the sequences of the four long helices (the A, B, C and D helices), the N- and C-terminal loops and the connecting sequences between the helices. The peptides corresponding to the helices of the folded protein all exhibit preferences for helix-like conformations in solution. The conformational ensembles of the A- and D-helix peptides contain ordered helical forms, as shown by extensive series of medium-range nuclear Overhauser effect connectivities, while the B- and C-helix peptides exhibit conformational preferences for nascent helix. All four peptides adopt ordered helical conformations in mixtures of trifluoroethanol and water. The terminal and interconnecting loop peptides also appear to contain appreciable populations of conformers with backbone phi and psi angles in the alpha-region and include highly populated hydrophobic cluster and/or turn conformations in some cases. Trifluoroethanol is unable to drive these peptides towards helical conformations. Overall, the peptide fragments of myohemerythrin have a marked preference towards secondary structure formation in aqueous solution. In contrast, peptide fragments derived from the beta-sandwich protein plastocyanin are relatively devoid of secondary structure in aqueous solution (see accompanying paper). These results suggest that the two different protein structural motifs may require different propensities for formation of local elements of secondary structure to initiate folding, and that there is a prepartitioning of conformational space determined by the local amino acid sequence that is different for the helical and beta-sandwich structural motifs.     

Tryptophan rich peptides: influence of indole rings on backbone conformation

Synthetic peptides with defined secondary structure scaffolds, namely hairpins and helices, containing tryptophan residues, have been investigated in this study to probe the influence of a large number of aromatic amino acids on backbone conformations. Solution NMR investigations of Boc-W-L-W-(D)P-G-W-L-W-OMe (peptide 1), designed to form a well-folded hairpin, clearly indicates the influence of flanking aromatic residues at the (D)Pro-Gly region on both turn nucleation and strand propagation. Indole-pyrrolidine interactions in this peptide lead to the formation of the less-frequent type I' turn at the (D)Pro-Gly segment and frayed strand regions, with the strand residues adopting local helical conformations. An analog of peptide 1 with an Aib-Gly turn-nucleated hairpin (Boc-W-L-W-U-G-W-L-W-OMe (peptide 2)) shows a preference for helical structures in solution, in both chloroform and methanol. Peptides with either one (Boc-W-L-W-U-W-L-W-OMe (peptide 3)) or two (Boc-U-W-L-W-U-W-L-W-OMe (peptide 4)) helix-nucleating Aib residues give rise to the well-folded helical conformations in the chloroform solution. The results are indicative of a preference for helical folding in peptides containing a large number of Trp residues. Investigation of a tetrapeptide analog of peptide 2, Boc-W-U-G-W-OMe (peptide 5), in solution and in the crystal state (by X-ray diffraction), also indicates a preference for a helical fold. Additionally, peptide 5 is stabilized in crystals by both aromatic interactions and an array of weak interactions. Examination of Trp-rich sequences in protein structures, however, reveals no secondary structure preference, suggesting that other stabilizing interactions in a well-folded protein may offset the influence of indole rings on backbone conformations.     

Conformational features responsible for the binding of cyclic analogues of enkephalin to opioid receptors. III. Probable binding conformations of mu-agonists with phenylalanine in position 3.

Theoretical conformational analysis was carried out for several tetrapeptide analogues of beta-casomorphin and dermorphin containing a Phe residue in position 3. Sets of low-energy backbone structures of the mu-selective peptides [N-Me-Phe3, D-Pro4]-morphiceptin and Tyr-D-Orn-Phe-Asp-NH2 were obtained. These sets of structures were compared for geometrical similarity between themselves and with the low-energy conformations found for the delta-selective peptide Tyr-D-Cys-Phe-D-Pen-OH and nonactive peptide Tyr-Orn-Phe-Asp-NH2. Two pairs of geometrically similar conformations of mu-selective peptides, sharing no similarity with the conformations of peptides showing low affinity to the mu-receptor, were selected as two alternative models of probable mu-receptor-bound backbone conformations. Both models share geometrical similarity with the low-energy structures of the linear mu-selective peptide Tyr-D-Ala-Phe-Phe-NH2. Putative binding conformations of Tyr1 and Phe3 side chains are also discussed.     

Structural and biological effects of a beta2- or beta3-amino acid insertion in a peptide.

Molecular mechanics calculations on conformers of Ac-HGly-NHMe, Ac-beta2-HAla-NHMe and Ac-beta3-HAla-NHMe indicate that low-energy conformations of the beta-amino acids backbone, corresponding to gauche rotamers around the Calpha-Cbeta bond, may overlap canonical backbone conformers observed for alpha-amino acids. Therefore, Substance P (SP) was used as a model peptide to analyse the structural and biological consequences of the substitution of Phe7 and Phe8 by (R)-beta2-HPhe and of Gly9 by HGly (R)-beta2-HAla or (S)-beta3-HAla. [(R)-beta2-HAla9]SP has pharmacological potency similar to that of SP while [HGly9]SP and [(S)-beta3-HAla9]SP show a 30- to 50-fold decrease in biological activities. The three analogues modified at position 9 are more resistant to degradation by angiotensin converting enzyme than SP and [Ala9]SP. NMR analysis of these SP analogues suggest that a beta-amino acid insertion in position 9 does not affect the overall backbone conformation. Altogether these data suggest that [HGly9]SP, [(S)-beta3-HAla9]SP and [(R)-beta2-HAla9]SP could adopt backbone conformations similar to that of SP, [Ala9]SP and [Pro9]SP. In contrast, incorporation of beta2-HPhe in position 7 and 8 of SP led to peptides that are almost devoid of biological activity. Thus, a beta-amino acid could replace an alpha-amino acid within the sequence of a bioactive peptide provided that the additional methylene group does not cause steric hindrance and does not confine orientations of the side chain to regions of space different from those permitted in the alpha-amino acid.     

Configurationally driven folding of model tetrapeptides containing L- or D-morpholine-3-carboxylic acids as beta-turn nucleators

The conformational analysis by NMR, IR, and molecular modeling of tetrapeptides containing morpholine-3-carboxylic acid (Mor) as a proline surrogate is presented. The relationship between the chirality of the cyclic amino acid at position i+1 and the turn propensity is maintained with respect to the reference proline-containing peptides, although marked differences in the type of folded structures were observed. The conformational profile of morpholine-containing turn peptides as a function of the chirality of the cyclic amino acid indicated that the heterochiral tetrapeptide containing the D-isomer of the cyclic amino acid is more prone to nucleate compact folded structures, although with no resemblance to the beta-turn structures of D-proline-containing peptides. Also, the solvation system proved to influence the organization of folded structures, as in the more interactive CD(3)CN the model peptides showed more compact conformations. The L-Mor-containing peptide displayed two rotamers at the Val-Mor amide bond. The trans isomer did not experience any turn structures, nor any intramolecular hydrogen-bonds, whereas the cis isomer showed a strong preference for a type VI beta-turn structure, thus providing a different conformational asset with respect to the beta-turn structure as reported for the reference L-proline model peptide.     

Dihedral angle preferences of amino acid residues forming various non-local interactions in proteins

In theory, a polypeptide chain can adopt a vast number of conformations, each corresponding to a set of backbone rotation angles. Many of these conformations are excluded due to steric overlaps. Ramachandran and coworkers were the first to look into this problem by plotting backbone dihedral angles in a two-dimensional plot. The conformational space in the Ramachandran map is further refined by considering the energetic contributions of various non-bonded interactions. Alternatively, the conformation adopted by a polypeptide chain may also be examined by investigating interactions between the residues. Since the Ramachandran map essentially focuses on local interactions (residues closer in sequence), out of interest, we have analyzed the dihedral angle preferences of residues that make non-local interactions (residues far away in sequence and closer in space) in the folded structures of proteins. The non-local interactions have been grouped into different types such as hydrogen bond, van der Waals interactions between hydrophobic groups, ion pairs (salt bridges), and ππ-stacking interactions. The results show the propensity of amino acid residues in proteins forming local and non-local interactions. Our results point to the vital role of different types of non-local interactions and their effect on dihedral angles in forming secondary and tertiary structural elements to adopt their native fold.     

Carbonyl-carbonyl interactions stabilize the partially allowed Ramachandran conformations of asparagine and aspartic acid

Asparagine and aspartate are known to adopt conformations in the left-handed alpha-helical region and other partially allowed regions of the Ramachandran plot more readily than any other non-glycyl amino acids. The reason for this preference has not been established. An examination of the local environments of asparagine and aspartic acid in protein structures with a resolution better than 1.5 A revealed that their side-chain carbonyls are frequently within 4 A of their own backbone carbonyl or the backbone carbonyl of the previous residue. Calculations using protein structures with a resolution better than 1.8 A reveal that this close contact occurs in more than 80% of cases. This carbonyl-carbonyl interaction offers an energetic sabilization for the partially allowed conformations of asparagine and aspartic acid with respect to all other non-glycyl amino acids. The non-covalent attractive interactions between the dipoles of two carbonyls has recently been calculated to have an energy comparable to that of a hydrogen bond. The preponderance of asparagine in the left-handed alpha-helical region, and in general of aspartic acid and asparagine in the partially allowed regions of the Ramachandran plot, may be a consequence of this carbonyl-carbonyl stacking interaction.     

AMBER force-field parameters for phosphorylated amino acids in different protonation states: phosphoserine, phosphothreonine, phosphotyrosine, and phosphohistidine

We report a consistent set of AMBER force-field parameters for the most common phosphorylated amino acids, phosphoserine, phosphothreonine, phosphotyrosine, and phosphohistidine in different protonation states. The calculation of atomic charges followed the original restrained electrostatic potential fitting procedure used to determine the charges for the parm94/99 parameter set, taking α-helical and β-strand conformations of the corresponding ACE-/NME-capped model peptide backbone into account. Missing force-field parameters were taken directly from the general AMBER force field (gaff) and the parm99 data set with minor modifications, or were newly generated based on ab initio calculations for model systems. Final parameters were validated by geometry optimizations and molecular-dynamics simulations. Template libraries for the phosphorylated amino acids in Leap format and corresponding frcmod parameter files are made available. Figure Schematic illustration of the systems used for parameter generation. Acid hydrogens are shown in red Electronic Supplementary Material Supplementary material is available for this article at