Combined use of molecular dynamics simulations and NMR to explore peptide bond isomerization and multiple intramolecular hydrogen-bonding possibilities in a cyclic pentapeptide, cyclo(Gly-Pro-D-Phe-Gly-Val).

The conformational behavior of a model cyclic pentapeptide--cyclo(Gly-L-Pro-D-Phe-Gly-L-Val)--has been explored through the combined use of in vacuo molecular dynamics simulations and a range of nmr experiments (preceding paper). The molecular dynamics analysis suggests that, despite the conformational constraints imposed by formation of the pentapeptide cycle, this pentapeptide undergoes conformational transitions between various hydrogen-bonded conformations, characterized by low energy barriers. An inverse gamma turn with Pro in position i + 1 and a gamma turn with D-Phe in position i + 1 are two alternatives occurring frequently. Like other DLDDL cyclic pentapeptides, cyclo(Gly-Pro-D-Phe-Gly-Val) is also stabilized by an inverse gamma-turn structure with the beta-branched Val residue in position i + 1, and this hydrogen bond is retained in the different conformational families. The gamma-turn around D-Phe3 and the inverse gamma turn around Val5 are consistent with the nmr observations. 3JNH-CH alpha coupling constants of the all-trans forms were calculated from one of the molecular dynamics trajectories and are comparable to nmr experimental data, suggesting that the conformational states visited during the simulation are representative of the conformational distribution in solution. In addition to the equilibrium among various hydrogen-bonded all-trans conformers, the observation in nmr spectra of two sets of resonances for all peptide protons indicated a slow conformational interconversion of the Gly-Pro peptide bond between trans and cis isomers. The activation energy between these two conformers was determined experimentally by magnetization transfer and was calculated by high temperature constrained molecular dynamics simulation. Both methods yield a free energy of activation of ca. 20 kcal/mol. Furthermore, the free energy of activation is dependent on the direction of rotation of the Gly-Pro peptide bond.     

Conformational analysis of endomorphin-1 by molecular dynamics methods.

Endomorphin-1 (EM1, H-Tyr-Pro-Trp-Phe-NH2) is a highly potent and selective agonist for the mu-opioid receptor. A conformational analysis of this tetrapeptide was carried out by simulated annealing and molecular dynamics methods. EM1 was modeled in the neutral (NH2-) and cationic (NH-) forms of the N-terminal amino group. The results of NMR measurements were utilized to perform simulations with restrained cis and trans Tyr1-Pro2 peptide bonds. Preferred conformational regions in the Phi 2-Psi 2, Phi 3-Psi 3 and Phi 4-Psi 4 Ramachandran plots were identified. The g(+), g(-) and trans rotamer populations of the side-chains of the Tyr1, Trp3 and Phe4 residues were determined in chi 1 space. The distances between the N-terminal N atom and the other backbone N and O atoms, and the distances between the centers of the aromatic side-chain rings and the Pro2 ring were measured. The preferred secondary structures were determined as different types of beta-turns and gamma-turns. In the conformers of trans-EM1, an inverse gamma-turn can be formed in the N-terminal region, but in the conformers of cis-EM1 the N-terminal inverse gamma-turn is absent. Regular and inverse gamma-turns were observed in the C-terminal region in both isomers. These beta- and gamma-turns were stabilized by intramolecular H-bonds and bifurcated H-bonds.     

Preferred conformations of cyclic Ac-Cys-Pro-Xaa-Cys-NHMe peptides: a model for chain reversal and active site of disulfide oxidoreductase

The conformational study on cyclic Ac-Cys-Pro-Xaa-Cys-NHMe (Ac-CPXC-NHMe; X=Ala, Val, Leu, Aib, Gly, His, Phe, Tyr, Asn and Ser) peptides has been carried out using the Empirical Conformational Energy Program for Peptides, version 3 (ECEPP/3) force field and the hydration shell model in the unhydrated and hydrated states. This work has been undertaken to investigate structural implications of the CPXC sequence as the chain reversal for the initiation of protein folding and as the motif for active site of disulfide oxidoreductases. The backbone conformation DAAA is commonly the most feasible for cyclic CPXC peptides in the hydrated state, which has a type I beta-turn at the Pro-Xaa sequence. The proline residue and the hydrogen bond between backbones of two cystines as well as the formation of disulfide bond appear to play a role in stabilizing this preferred conformation of cyclic CPXC peptides. However, the distributions of backbone conformations and beta-turns may indicate that the cyclic CPXC peptide seems to exist as an ensemble of beta-turns and coiled conformations in aqueous solution. The intrinsic stability of the cyclic CPXC motif itself for the active conformation seems to play a role in determining electrochemical properties of disulfide oxidoreductases.     

The preferred conformation of the tripeptide Ala-Phe-Ala in water is an inverse gamma-turn: implications for protein folding and drug design

Recent studies have provided evidence that peptides as short as tripeptides do adopt preferred conformations. Here we report that the tripeptide Ala-Phe-Ala (AFA) in aqueous solution preferentially forms an inverse gamma-turn. Circular dichroism (CD) indicated the presence of a predominant turn structure, and Fourier transform infrared (FTIR) bands suggested the presence of a gamma-turn forming a bifurcated H-bond with the solvent molecules. The high-resolution structure was obtained by a combined use of NMR spectroscopy and calculations. On the basis of 30 unambiguous ROESY-derived distance restraints (including the Halpha-NH NOE between Ala(1) and Ala(3) and a hydrogen bond between the CO group of Ala(1) and the NH group of Ala(3)), calculations clearly demonstrated the presence of an inverse gamma-turn centered on Phe(2). From NOE data, we estimated a mole fraction for the gamma-turn of 0.65. Since for AFA an extended beta-strand was also reported [Eker, F., Griebenow, K., Cao, X., Nafie, L. A., and Schweitzer-Stenner, R. (2004) Proc. Natl. Acad. Sci. U.S.A. 101, 10054-10059], we investigated the possibility that gamma-turn and beta-strand may represent two major conformations. By using a best-fit procedure that calculated experimental NOEs as weighted averages of the effects originating from both structures, we were able to calculate with good accuracy the backbone NOEs at 280 K in terms of the two limiting conformers, yielding a mole fraction for the gamma-turn and beta-strand conformations of 0.60 and 0.40, respectively, in good agreement with those found by NOE data. The implication of the existence of a preferred conformation by a small structural element is discussed in the context of the nucleation of protein folding events and the design of small peptide and peptidomimetic drugs.     

Backbone modifications in cyclic peptides. Conformational analysis of a cyclic pseudopentapeptide containing a thiomethylene ether amide bond replacement

NMR and X-ray crystallographic studies have shown that cyclic pentapeptides of the general structure cyclo(D-Xxx-Pro-Gly-Pro-Gly) possess beta- and gamma-turn intramolecular hydrogen bonds. As part of our continuing series surveying the compatibility of various amide bond replacements on peptide structure, we have synthesized cyclo(D-Phe-Pro psi[CH2S]Gly-Pro-Gly). The pseudopeptide was prepared by solid phase methods and cleaved from the resin by a new procedure involving phase transfer catalysis using K2CO3 and tetrabutylammonium hydrogen sulfate. Cyclization was carried out with the use of DPPA, HOBt, and DMAP to afford the product in 69% yield. The conformational behavior of the pseudopeptide was analyzed by 1H and 13C (1D and 2D) NMR techniques. The backbone modification replaced the amide bond that is involved in a gamma-turn intramolecular hydrogen bond in the all-amide structure. In CDCl3, the pseudopeptide adopted the same all-trans conformation as its parent, although the remaining beta-turn hydrogen bond was weaker according to delta delta/delta TNH measurements. In DMSO-d6, the all-trans conformer and a second conformer were observed in a ratio of 55:45. These conformers, which slowly interconverted on the NMR time scale, could be separately assigned; peaks due to chemical exchange were readily distinguishable by the ROESY technique as reported earlier by others. 13C and ROESY experiments suggested the minor conformer contained one cis amide bond at the Gly1-Pro2 position. Thus, both the location and type of amide surrogate are important determinants affecting the compatibility of the replacement with a particular conformational feature.     

Side-chain conformational entropy in protein folding.

An important, but often neglected, contribution to the thermodynamics of protein folding is the loss of entropy that results from restricting the number of accessible side-chain conformers in the native structure. Conformational entropy changes can be found by comparing the number of accessible rotamers in the unfolded and folded states, or by estimating fusion entropies. Comparison of several sets of results using different techniques shows that the mean conformational free energy change (T delta S) is 1 kcal.mol-1 per side chain or 0.5 kcal.mol-1 per bond. Changes in vibrational entropy appear to be negligible compared to the entropy change resulting from the loss of accessible rotamers. Side-chain entropies can help rationalize alpha-helix propensities, predict protein/inhibitor complex structures, and account for the distribution of side chains on the protein surface or interior.     

Single residue modifications of the delta opioid receptor selective peptide, [D-Pen2,D-Pen5]-enkephalin (DPDPE). Correlation of pharmacological effects with structural and conformational features

Six analogs of the highly delta opioid receptor selective, conformationally restricted, cyclic peptide [D-Pen2,D-Pen5]enkephalin, Tyr-D-Pen-Gly-Phe-D-PenOH (DPDPE), were synthesized and evaluated for opioid activity in rat brain receptor binding and mouse vas deferens (MVD) smooth muscle assays. All analogs were single amino acid modifications of DPDPE and employed amino acid substitutions of known effects in linear enkephalin analogs. The effect on binding affinity and MVD potency of each modification within the DPDPE structural framework was consistent with the previous reports on similarly substituted linear analogs. Conformational features of four of the modified DPDPE analogs were examined by 1H NMR spectroscopy and compared with DPDPE. From these studies it was concluded that the observed pharmacological differences with DPDPE displayed by diallyltyrosine1-DPDPE ([DAT1]DPDPE) and phenylglycine4-DPDPE ([Pgl4]DPDPE) are due to structural and/or conformational differences localized near the substituted amino acid. The observed enhanced mu receptor binding affinity of the carboxamide terminal DPDPE-NH2 appears to be founded solely upon electronic differences, the NMR data suggesting indistinguishable conformations. The observation that the alpha-aminoisobutyric acid substituted analog [Aib3]DPDPE displays similar in vitro opioid behavior as DPDPE while apparently assuming a significantly different solution conformation suggests that further detailed conformational analysis of this analog will aid the elucidation of the key structural and conformational features required for action at the delta opioid receptor.     

Equilibrium of the Cis-Trans Isomerisation of the Peptide Bond with N-alkyl Amino Acids Measured by 2D NMR

The conformational cis-trans equilibrium around the peptide bond in model tripeptides has been determined by 2D NMR methods (HOHAHA, ROESY). The study was limited to three different N-substituted amino acids in position 2, namely Pro (proline), Tic (1,2,3,4-tetrahydroisoquinoline-3-carboxylic acid), and N-MePhe (N-methylphenylalanine). In all cases the amino acid in position 1 was tyrosine and in position 3, phenylalanine. The results of our studies show that the cis-trans ratio depends mostly on the configuration of the amino acids forming the peptide bond undergoing the cis-trans isomerisation. The amino acid following the sequence (in position 3) does not have much influence on the cis-trans isomerisation, indicating that there is no interaction of the side chains between these amino acids. The model peptides with the L-Tyr-L-AA-(L- or D-)Phe (where AA is N-substituted amino acid) chiralities give 80–100% more of the cis form in comparison to the corresponding peptides with the D-Tyr-L-AA-(L-or D-)Phe chiralities. These results indicate that the incorporation of N-substituted amino acids in small peptides with the same chirality as the precedent amino acid involved in the peptide bond undergoing the cis/trans isomerisation moves the equilibrium to a significant amount of the cis form.     

Empirical rules predicting conformation of cyclic tetrapeptides from primary structure

A useful set of empirical rules is put forward to predict the conformations of cyclic tetrapeptides and cyclic tetradepsipeptides on the basis of primary structure, briefly presented as follows: A conformation allowing an intramolecular hydrogen bond (IMHB) of gamma-turn is preferred, and an ester bond always adopts a trans form. On a right-handed peptide ring, the carbonyl group acylating a D residue is oriented to the upper side of the main ring. The carbonyl group acylating a D proline or an N-methyl-D-amino acid residue is oriented to the lower side of the ring, forming a cis bond. The LDDL configurational sequence adopts a cis-trans-cis-trans backbone with Ci symmetry. A glycine residue behaves as a D residue in an L-peptide. Conformations of cyclotetrapeptides containing two glycine residues at diametric positions or containing an N-methyl-dehydroamino acid residue are predicted by use of appendices of rule 5. Almost all conformations of cyclic tetrapeptides are predicted by these rules. Energetical rationalization of the rules and prediction of possible new conformations are described. Conformations of cyclo(-L-Pro-L-Leu-D-Tyr(Me)-L-Ile-)(1) and cyclo (-L-Pro-D-Leu-D-Tyr(Me)-L-Ile)(2) are compared. Results of n.m.r. experiments showed that compound 1 adopts a unique cis-trans-trans-trans backbone with a gamma-turn IMHB, and 2 has a cis-trans-cis-trans backbone with Ci symmetry. These observations confirmed the rules described above. Peptides 1 and 2 are the first diastereomeric peptides with trans (LD) and cis (DD) secondary amide bonds.     

Synthesis, receptor binding affinities and conformational properties of cyclic methylenedithioether analogues of angiotensin II

Cyclic 12-, 13- and 14-membered ring angiotensin II analogues related to disulfides but encompassing methylene-dithioether bridges have been prepared. The affinity data from these derivatives were compared to those from the disulfides. The methylenedithioether analogues displayed good binding affinities to rat liver AT1 receptors although in most cases somewhat lower than their disulfide counterparts. One of the methylenedithioethers with a 13-membered ring system demonstrated the highest binding affinity among the thioethers. Theoretical conformational analysis of model compounds of the two series were performed suggesting a similarity between the disulfide and the corresponding methylenedithioether analogues and also between the ring size homologues. This analysis also suggested that some of the model compounds were prone to adopt inverse gamma-turn conformations, which was further supported by use of NMR spectroscopy of the 12-membered ring analogue in the series. The easily executed methylenedithioether cyclization should constitute a valuable complement to the common disulfide methodology for fine-tuning and for probing the bioactive conformation of peptides.     

Synthesis and properties of L-cysteinyl-L-cysteine disulfides

Oligomeric cyclic disulfides, obtained by mild oxidation of the fully protected dipeptide L -cysteinyl-L -cysteine, have been isolated by gel and thin-layer chromatography. Polymeric material was recycled by a thiol–disulfide exchange-reaction performed at basic pH. Spectroscopic investigations of the monomer and the two dimers indicate that conformers characterized by dihedral angles about the S? S bond close to ±90° are preferred. Moreover, chiroptical and ¹H-nmr data for these compounds suggest higher mobility for the two dimers. The antiparallel dimeric disulfide can be considered a model compound for the hinge region formed at the subunit interface of the bovine seminal ribonuclease, a dimeric enzyme showing a complex kinetic behavior.     

Solution conformation of tuftsin.

Tuftsin, a natural linear tetrapeptide (Thr-Lys-Pro-Arg) of potential antitumor activity, has been studied in DMSO-d6 solution by 2D NMR spectroscopy. 1H and 13C spectra show the presence of two families of conformations characterized by a trans or cis Lys-Pro bond, respectively. The family of conformers containing the cis peptide bond is a mixture of extended structures as expected for a short linear peptide. On the contrary, the trans isomer appears to be a rigid, folded conformer, as indicated by crucial NOEs and by the exceptionally low temperature coefficient of Arg NH. Analysis of the solution data by means of energy calculations leads to a unique structure, characterized by a Lys-Pro inverse gamma-turn.     

Vibrational analysis of peptides, polypeptides, and proteins. XXX. Normal mode analyses of gamma-turns

The normal modes have been calculated for three kinds of low energy gamma-turn structures resulting from recent conformational energy calculations by Némethy. Frequencies have been computed for a gamma-turn, a mirror-related gamma-turn, and an inverse gamma-turn of CH3-CO-(L-Ala)n-NH-CH3, with n = 3 and n = 5, and for certain 14C and 15N derivatives of the n = 3 molecule. Correlations are evident between amide frequencies and gamma-turn structures, and it is found that only amide I modes of peptide groups in the turn are relatively insensitive to the lengths of attached chains.     

New 21- and 24-atom Aib-containing cyclopeptides.

The synthesis of two Aib-containing cyclic octapeptides, cyclo(Leu-Aib-Phe-Aib-Gly-Aib-Phe-Aib) (1) and cyclo(Gly-Aib-Aib-Val-Aib-Leu-Aib-Phe) (2), and a cyclic heptapeptide cyclo(Gly-Aib-Val-Aib-Leu-Aib-Phe) (3), is described. The linear precursors of 1-3 were prepared using solution-phase techniques, and the cyclisation was also accomplished in solution. Among the coupling reagents examined in the final macrolactamisation step, PyAOP, HATU and DEPC/DEPBT efficiently yielded cyclised products. However, the success of the cyclisation was found to be dependent on the coupling reagent used. The two octapeptides 1 and 2 were obtained in much better yields (up to 63%) than the cycloheptapeptide 3 (30-37%). In addition, crystal-state conformational analysis of 2 was performed by X-ray crystallography. Six intramolecular hydrogen bonds stabilise the structure characterised by two consecutive type II'/I beta-turns, two consecutive type II/III' beta-turns, one gamma-turn, and one inverse gamma-turn. Copyright (c) 2008 European Peptide Society and John Wiley & Sons, Ltd.     

Effect of peptide conformation on membrane permeability.

The effect of peptide conformational constraint on the peptide permeation across the model membranes was examined by determining the permeability of pairs of cyclic and acyclic peptides related to c[d-Pen2, d-Pen5] enkephalin (DPDPE). The peptides were cyclized by formation of an intramolecular disulfide bridge between the second and fifth residues composed of either d-penicillamine or cysteine. In each case the acyclic peptide was three to seven times more permeable than corresponding cyclic peptide. The possibility that the differences in permeability of cyclic and acyclic peptides is based on the greater conformational freedom of the acyclic peptides in the presence of membrane was examined in more detail by isothermal titration calorimetric studies of Trp6-DPDPE and its acyclic analog. The membrane binding of the acyclic peptide is a more exothermic process than binding of its cyclic Trp6-DPDPE. The transfer of acyclic peptide from water to membrane is an enthalpy driven process, whereas the transfer of the cyclic peptide is driven by entropy.     

Depsipeptide analogues of elastin repeating sequences: conformational analysis

In this work the effect of elimination of a specific hydrogen bond on the conformation of the repeating peptides of elastin was studied. These repeating sequences are the pentapeptide Val-Pro-Gly-Val-Gly and the hexapeptide Val-Ala-Pro-Gly-Val-Gly. These sequences have been proposed to occur in a beta-turn conformation with a hydrogen bond involving the amide NH of the internal valine residue and the carbonyl oxygen of the residue preceding proline. In the depsipeptide analogues studied in this work, this 4-1 beta-turn hydrogen bond cannot occur. We studied the depsipeptide sequences Val-Pro-Gly-Hiv-Gly and Val-Ala-Pro-Gly-Hiv-Gly (Hiv denotes S-alpha-hydroxyisovaleric acid, the hydroxy acid analogue of valine), as well as the peptide sequences Val-Pro-Gly-Val-Gly and Val-Ala-Pro-Gly-Val-Gly. Compounds studied included sequences with the Boc and benzyl ester protecting groups, derivatives with the acetyl and N-methylamide end groups and polymers of the above sequences. Our conclusions are based on a comparison of depsipeptides with analogous peptides. Conformational analysis was carried out by nmr, CD, and ir spectroscopy. We propose that in the repeating sequences of elastin an equilibrium exists between a gamma-turn structure and a beta-turn structure in the Pro-Gly segment resulting in a structure that combines flexibility with strong conformational preferences. The C7 involves the amide NH of the internal glycine and the carbonyl oxygen of the residue preceding proline. In the N-methylamide derivatives a similar equilibrium exists in the Gly-Val-Gly segment. In the depsipeptides the beta-turn cannot occur and only the gamma-turn is seen. In the polydepsipeptides the major conformational feature is a type I beta-turn involving Gly5 NH and Pro CO.     

Synthetic and conformational studies on dehydroalanine-containing model peptides

Three model dipeptides containing a dehydroalanine residue (delta Ala) at the C-terminal, Boc-X-delta Ala-NHCH3 [X = Ala, Val, and Phe,] have been synthesized and their solution conformations investigated by 1H-NMR, IR, and CD spectroscopy. NMR studies on these peptides in CDCl3 clearly indicate that the NH group of dehydroalanine is involved in an intramolecular hydrogen bond. This conclusion is supported by IR studies also. Nuclear Overhauser effect (NOE) studies are also accommodative of an inverse gamma-turn-type of conformation that is characterized by conformational angles of phi approximately -70 degrees and psi approximately +70 degrees around the X residue, and a C alpha i + 1 H-Ni + 2H interproton distance of 2.5 A. It appears that unlike dehydrophenylalanine or dehydroleucine, which tend to stabilize beta-turn type of structures occupying the i + 2 position of the turn, dehydroalanine favors the formation of an inverse gamma-turn, centered at the preceding L-residue in such solvents as CDCl3 and (CD3)2SO. A comparison of solution conformation of Boc Val-delta Ala-NHCH3 with the corresponding saturated analogue, Boc-Val-Ala-NHCH3, is also presented and shows that dehydroalanine is responsible for inducing the turn structure. It may be possible to design peptides with different preferred conformations using the suitable dehydroamino acid.     

Dermorphin tetra- and heptapeptide analogues containing a [3,4] amide bond replacement by a carbon-carbon double and single bond.

Seven dermorphin hepta- and tetrapeptide analogues containing [3,4] amide bond replacement by a carbon-carbon double and single bond were prepared. 1H NMR studies of the pseudoheptapeptide in DMSO indicate the presence of extended conformations with stacking of the side chains in the N-terminal part and an inverse gamma-turn around Ser7 in the conformational equilibrium. The binding data show that the affinity of the analogues for the mu-receptor is only slightly diminished in the D-Ala2 series and is more affected in the D-Arg2 series. Since the Gly4NH is not present in these compounds we conclude that this NH is not required to stabilize the bioactive conformation nor is it directly involved in binding to the receptor.     

Dynamical search for bis-penicillamine enkephalin conformations.

Quenched molecular dynamics is used as a conformational search technique for the constrained cyclic analog [D-Pen2,D-Pen5]enkephalin (DPDPE) in a continuum solvent. The results show a Gaussianlike distribution of conformations as a function of energy, unlike the distributions found for simple liquids which have sharp bands for different crystal forms and broad glasslike states are found. The lowest energy conformers have structural features in common with those obtained from constrained searches based on energy minimization. (Hruby, V. J., L-.F. Kao, B. M. Pettitt, and M. Karplus. 1988. J. Am. Chem. Soc. 110:3351-3359). Many of the low energy configurations are amphiphilic with the carbonyl groups on one surface and the hydrophobic groups on the other. This supports the conclusions from the previous modeling study, which yielded amphiphilic structures as the most probable conformations of DPDPE when NOE data were included.     

Beta-alanine containing peptides: a novel molecular tool for the design of gamma-turns.

In the present paper we describe the synthesis, purification, single crystal x-ray analysis, and solution conformational characterization of the cyclic tetrapeptide cyclo-(L-Pro-beta-Ala-L-Pro-beta-Ala). This peptide was synthesized by classical solution methods and the cyclization of the free tetrapeptide was accomplished in good yields in diluted methylene chloride solution using N,N-dicyclohexyl-carbodiimide (DCCI). The compound crystallizes in the orthorombic space group P2(1)2(1)2(1) from ethyl acetate. All peptide bonds are trans. The molecular conformation is stabilized by two intramolecular hydrogen bonds between the CO and NH groups of the two beta-alanine residues. These hydrogen bonds take place in a C7 structure in which both proline residues occupy the 2 position of an inverse gamma-turn. The two beta-alanine residues have a typical folded conformation (around the C alpha-C beta bond) observed in other cyclic peptides containing this residue. A detailed 1H-nmr analysis in CD3CN solution has been carried out. The molecule assumes a twofold symmetry in solution with a molecular conformation consistent with that observed in the solid state.