Optical and conformational properties of the 7-(beta-D-ribofuranosyl) purines

The optical and conformational properties of certain 7-ribosyl purines have been studied using several experimental and theoretical methods. The CD, MCD, and absorption spectra of the 7-ribosyl purines and their cations have been determined experimentally with some solvent effects being noted. Theoretically we have used the SCF-CI and CNDO molecular orbital calculations to determine the optical properties of the bases; Simpson's bond exciton theory to determine the optical properties of the bases; Simpson's bond exciton theory to determine the optical properties of the ribose moiety; the optical activity was determined by the dipole coupling and electricmagnetic coupling equations; and the conformational energy calculations are included as a basis for excluding highly improbable conformations. On the basis of these calculations, we concluded that the oxy derivatives and the amino derivatives are predominately in the “extreme anti” conformation and “standard anti” conformation, respectively. However, the conclusions are tentative since at the present time the band assignments are not unequivocal. At least one n-π* transition has been identified in the spectra which is not notably sensitive to pH and solvent effects. Several implications of this observation are discussed.     

Location of proline derivatives in conformational space. II. Theoretical optical activity

The coupling of theoretical optical calculations with experimental data provides a check of the validity of the theory or provides conformational information. The theory was validated by studies in which the approximate conformation was located independently. These studies have shown that a theory restricted to the two lowest energy transitions for each chromophore gives qualitative agreement with experiment. On the other hand, for some of the proline derivatives, the theoretical treatment allows detailed conformational assignments. Both successes and failures in correlating theory with experiment are discussed. The results presented provide a basis for assessing the prospects for relating protein and polypeptide optical activity to their conformations.     

The conformational and vibrational behavior of the inhibitory neuropeptide derived from beta-endorphin

In this study, conformational behavior, structural, and vibrational characterization of the carboxy terminal dipeptide of β-endorphin (glycy-l-glutamine, glycyl-glutamine, beta-endorphin_30-31), which is an inhibitory neuropeptide synthesized from beta-endorphin_1-31 in brain stem regions, has been investigated. The theoretically possible stable conformers were searched by means of molecular mechanics method to determine their energetically preferred conformations. The 360 different conformations were calculated with the φ, Ψ, χ dihedral angles using the Ramachandran maps. The most stable conformation of the title molecule is characterized by the extended backbone shape (e) in the BR conformational range with ?.78 kcal/mol energy. The cis- and trans-dimeric forms of the dipeptide were also formed and energetically preferred conformations of dimers were investigated. The experimental methods (FT-IR, micro-Raman spectroscopies) coupled with quantum chemical calculations based on density functional theory (DFT) have been used to identify the geometrical, energetic, and vibrational characteristics of the dipeptide. The assignment of the vibrational spectra was performed based on the potential energy distribution of the vibrational modes. To investigate the electronic properties, such as nonlinear optical properties, the electric dipole moment, the mean polarizability, the mean first hyperpolarizability, and HOMO–LUMO energy gaps were computed using the DFT with the B3LYP/6-31++G(d,p) basis set combination. The second-order interaction energies were derived from natural bonding orbital analysis. The focus of this study is to determine possible stable conformation on inhibitory neuropeptide and to investigate molecular geometry, molecular vibrations of monomeric and dimeric forms, and hydrogen bonding interactions of glycy-l-glutamine dipeptide.     

Prediction of the native conformation of a polypeptide by a statistical-mechanical procedure. III. Probable and average conformations of enkephalin

The program SMAPPS (Statistical-Mechanical Algorithm for Predicting Protein Structure) was originally designed to determine the probable and average backbone (?, ψ) conformations of a polypeptide by the application of equilibrium statistical mechanics in conjunction with an adaptive importance sampling Monte Carlo procedure. In the present paper, the algorithm has been extended to include the variation of all side-chain (χ) and peptide-bond (ω) dihedral angles of a polypeptide during the Monte Carlo search of the conformational space. To test the effectiveness of the generalized algorithm, SMAPPS was used to calculate the probable and average conformations of Met-enkephalin for which all dihedral angles of the pentapeptide were allowed to vary. The total conformational energy for each randomly generated structure of Met-enkephalin was obtained by summing over the interaction energies of all pairs of nonbonded atoms of the whole molecule. The interaction energies were computed by the program ECEPP /2 (Empirical Conformational Energy Program for Peptides). Solvent effects were not included in the computation. The results of the Monte Carlo calculation of the structure of Met-enkephalin indicate that the thermodynamically preferred conformation of the pentapeptide contains a γ-turn involving the three residues Gly²-Gly³-Phe⁴. The γ-turn conformation, however, does not correspond to the structure of lowest conformational energy. Rather, the global minimum-energy conformation, recently determined by a new optimization technique developed in this laboratory, contains a type II′ β-bend that is formed by the interaction of the four residues Gly²-Gly³-Phe⁴-Met⁵. A similar minimum-energy conformation is found by the SMAPPS procedure. The thermodynamically preferred γ-turn structure has a conformational energy of 4.93 kcal/mole higher than the β-bend structure of lowest energy but, because of the inclusion of entropy in the SMAPPS procedure, it is estimated to be ～ 9 kcal/mole lower in free energy. The calculation of the average conformation of Met-enkephalin was repeated until a total of ten independent average conformations were established. As far as the phenylalanine residue of the pentapeptide is concerned, the results of the ten independent average conformations were all found to lie in the region of conformational space corresponding to the γ-turn. These results further support the conclusion that the γturn conformation is thermodynamically favored.     

Optical activity of aromatic chromophores. I. O, m, and p-tyrosine

The optical activity and the conformational energy of the amino acids tyrosine, o-tyrosine, and m-tyrosine have been calculated as a function of molecular conformation. A new graphic technique, rotatory strength–conformational energy maps, was developed for the presentation of these calculations. Experimental circular dichroism spectra were determined by utilization of a new experimental technique that involves complete computer control of a Cary spectropolarimeter. This permitted repetitive scanning with signal averaging over extended time periods and resulted in a greatly enhanced signal to noise ratio. On the basis of these calculations, it, was concluded that the Kirkwood “coupled oscillator” mechanism is capable of explaining the observed optical activity of the aromatic chromophores of the molecules that were investigated.     

Prediction of the native conformation of a polypeptide by a statistical-mechanical procedure. II. Average backbone structure of enkephalin

The average conformation of Met-enkephalin was determined by using an adaptive, importance-sampling Monte Carlo algorithm (SMAPPS—Statistical Mechanical Algorithm for Predicting Protein Structure). In the calculation, only the backbone dihedral angles (? and ψ) were allowed to vary; i.e., all side-chain (χ) and peptide-bond (ω) dihedral angles were kept fixed at the values corresponding to a low-energy structure of the pentapeptide. The total conformational energy for each randomly generated structure of the polypeptide was obtained by summing over the interaction energies of all pairs of nonbonded atoms of the whole molecule. The interaction energies were computed by the program ECEPP/2 (Empirical Conformational Energy Program for Peptides). Solvent effects were not included in the computation. The calculation was repeated until a total of 10 independent average conformations were established. The regions of conformational space occupied by the average structures were compared with the regions of low conditional free energy obtained by SMAPPS in the first paper of this series. Such a comparison provides an analysis of the capacity of SMAPPS to adjust the Monte Carlo search to regions of highest probability. The results demonstrate that the ability of SMAPPS to focus the Monte Carlo search is excellent. Finally, the 10 independent average conformations and the mean of the 10 average structures were utilized as the initial conformations for a direct energy minimization of the pentapeptide. Of the 11 final energy-minimized structures, three of the conformations were found to be equivalent to the conformation of lowest energy determined previously. In addition, all but two of the remaining energy-minimized structures were found to correspond to one of the two other conformations of high probability obtained in the first paper of this series. These results indicate that a set of independent average conformations can provide a rational, unbiased choice for the initial conformation, to be used in a direct energy minimization of a polypeptide. The final energy-minimized structures consequently constitute a set of low-energy conformations, which include the global energy minimum.     

Conformational analysis of the 20-residue membrane-bound portion of melittin by conformational space annealing

The conformational space of the 20-residue membrane-bound portion of melittin has been investigated extensively with the conformational space annealing (CSA) method and the ECEPP/3 (Empirical Conformational Energy Program for Peptides) algorithm. Starting from random conformations, the CSA method finds that there are at least five different classes of conformations, within 4 kcal/mol, which have distinct backbone structures. We find that the lowest energy conformation of this peptide from previous investigations is not the global minimum-energy conformation (GMEC); but it belongs to the second lowest energy class of the five classes found here. In four independent runs, one conformation is found repeatedly as the lowest energy conformation of the peptide (two of the four lowest energy conformations are identical; the other two have essentially identical backbone conformations but slightly different side-chain conformations). We propose this conformation, whose energy is lower than that found previously by 1.9 kcal/mol, as the GMEC of the ECEPP/3 force field. The structure of the proposed GMEC is less helical and more compact than the previous one. It appears that the CSA method can find several classes of conformations of a 20-residue peptide starting from random conformations utilizing only its amino acid sequence information. The proposed GMEC has also been found with a modified electrostatically driven Monte Carlo method [D. R. Ripoll, A. Liwo, and H.A. Scheraga (1998) “New Developments of the Electrostatically Driven Monte Carlo Method: Test on the Membrane-Bound Portion of Melittin,” Biopolymers, Vol. 46, pp. 117–126]. © 1998 John Wiley & Sons, Inc. Biopoly 46: 103–115, 1998     

Stereochemistry of nucleic acids and polynucleotides. VI. Minimum energy conformations of dimethyl phosphate

As part of a study on the conformation of polynucleotides and nucleic acids the preferred conformations of the model conpound dimethyl phosphate are worked out using potential energy functions. In calculating the total potential energy associated with the conformation, nonbonded, torsional, and electrostatic terms have been considered. The variation of the total conformational energy is represented as a function of two torsion angles ? and ψ which are the rotations about the two phosphoester bonds. The most stable conformations are found to be the gauche–gauche conformations about these bonds. The conformations observed for phosphodiesters in the solid state and in the proposed structures of polynucleotides and nucleic acids cluster around the minimum. Also, regions of minimum energy correspond well with the typical allowed regions of a representative dinucleotide.     

Conformational energy calculations of enzyme-substrate complexes of lysozyme. I. Energy minimization of monosaccharide and oligosaccharide inhibitors and substrates of lysozyme

Conformational energy calculations were performed on monosaccharide and oligosaccharide inhibitors and substrates of lysozyme to examine the preferred conformations of these molecules. A grid-search method was used to locate all of the low-energy conformational regions for N-acetyl-β-D -glycosamine (NAG), and energy minimization was then carried out in each of these regions. Three stable positions for the N-acetyl group have ben located, in two of which the plane of the amide unit is normal to the mean plane of the pyranosyl ring. Nine local energy minima were located for the —CH₂OH group. The positions of the two vicinal cis —OH groups are determined predominantly by interactions with either the —CH₂OH or the N-acetyl group. The most stable conformations of β-N-acetylmuramic acid (NAM) were determined from the study of the low-energy conformations of NAG. In the two stable orientations for the D -lactic acid side chain, the O—C—C′ plane (C′ being the carbon atom of the terminal carboxyl group) was found to be normal to the mean plane of the pyranosyl ring. The low-energy positions for the COOH group of NAM are determined mainly by interactions with neighboring groups. The conformational preferences of the α-anomers of NAG and NAM were also explored. The calculated conformation of the N-acetyl group for α-NAG was quite close to that determined by X-ray analysis. Two of the three lowest energy conformations of α-NAM are similar to the corresponding conformations of the β-anomer. A third low-energy structure, which has a hydrogen bond from the NH of the N-acetyl group to the C?O of the lactic acid group, corresponds very closely to the X-ray structure of this molecule. The preferred conformations of the disaccharides NAG–NAG, NAM–NAG and NAG–NAM were also investigated. Two preferred orientations of the reducing pyranosyl ring relative to the nonreducing ring were found for all of these disaccharides, both of which are close to the extended conformation. In one of these conformations, a hydrogen bond can form between the OH group attached to C₃ of the reducing sugar and the ring oxygen of the preceding residue. Each conformation can be stabilized further by a hydrogen bond between the CH₂OH (donor) of residue i + 1 and the C?O of residue i (acceptor). The interactions that determine conformations for all oligosaccharides containing both NAG and NAM are shown to be exclusively intraresidue and nearest neighbor interactions, so that it is possible to predict all stable conformations of oligosaccharides containing NAG and NAM in any sequence.     

Prediction of the native conformation of a polypeptide by a statistical-mechanical procedure. I. Backbone structure of enkephalin

A new methodology for theoretically predicting the native, three-dimensional structure of a polypeptide is presented. Based on equilibrium statistical mechanics, an algorithm has been designed to determine the probable conformation of a polypeptide by calculating conditional free-energy maps for each residue of the macromolecule. The conditional free-energy map of each residue is computed from a set of probability integrals, obtained by summing over the interaction energies of all pairs of nonbonded atoms of the whole molecule. By locating the region(s) of lowest free energy for each map, the probable conformation for each residue can be identified. The native structure of the polypeptide is assumed to be the combination of the probable conformations of the individual residues. All multidimensional probability integrals are evaluated by an adaptive Monte Carlo algorithm (SMAPPS —Statistical-Mechanical Algorithm for Predicting Protein Structure). The Monte Carlo algorithm searches the entire conformational space, adjusting itself automatically to concentrate its sampling in regions where the magnitude of the integrand is largest (“importance sampling”). No assumptions are made about the native conformation. The only prior knowledge necessary for the prediction of the native conformation is the amino acid sequence of the polypeptide. To test the effectiveness of the algorithm, SMAPPS was applied to the prediction of the native conformation of the backbone of Met-enkephalin, a pentapeptide. In the calculations, only the backbone dihedral angles (? and ψ) were allowed to vary; all side-chain (χ) and peptide-bond (ω) dihedral angles were kept fixed at the values corresponding to the alleged global minimum energy previously determined by direct energy minimization. For each conformation generated randomly by the Monte Carlo algorithm, the total conformational energy of the polypeptide was obtained from established empirical potential energy functions. Solvent effects were not included in the computations. With this initial application of SMAPPS , three distinct low-free-energy β-bend structures of Met-enkephalin were found. In particular, one of the structures has a conformation remarkably similar to the one associated with the previously alleged global minimum energy. The two additional structures of the pentapeptide have conformational energies lower than the previously computed low-energy structure. However, the Monte Carlo results are in agreement with an improved energy-minimization procedure. These initial results on the backbone structure of Met-enkephalin indicate that an equilibrium statistical-mechanical procedure, coupled with an adaptive Monte Carlo algorithm, can overcome many of the problems associated with the standard methods of direct energy minimization.     

Proline-induced constraints in alpha-helices

The disrupting effect of a prolyl residue on an α-helix has been analyzed by means of conformational energy computations. In the preferred, nearly α-helical conformations of Ac-Ala₄-Pro-NHMe and of Ac-Ala₇-Pro-Ala₇-NHMe, only the residue preceding Pro is not α-helical, while all other residues can occur in the α-helical A conformation; i.e., it is sufficient to introduce a conformational change of only one residue in order to accommodate proline in a distorted α-helix. Other low-energy conformations exist in which the conformational state of three residues preceding proline is altered considerably; on the other hand, another conformation in which these three residues retain the near-α-helical A-conformational state (with up to 26° changes of their dihedral angles ? and ψ, and a 48° change in one ω from those of the ideal α-helix) has a considerably higher energy. These conclusions are not altered by the substitution of other residues in the place of the Ala preceding Pro. The conformations of the peptide chain next to prolyl residues in or near an α-helix have been analyzed in 58 proteins of known structure, based on published atomic coordinates. Of 331 α-helices, 61 have a Pro at or next to their N-terminus, 21 have a Pro next to their C-terminus, and 30 contain a Pro inside the helix. Of the latter, 16 correspond to a break in the helix, 9 are located inside distorted first turns of the helix, and 5 are parts of irregular helices. Thus, the reported occurrence of prolyl residues next to or inside observed α-helices in proteins is consistent with the computed steric and energetic requirements of prolyl peptides.     

On the Bioactive Conformation of a Small Peptide and its Set of Thermodynamically Accessible Conformations

Abstract

In order to investigate the relationship between the bioactive conformation of a peptide and its set of thermodynamically accessible structures in solution, the conformational profile of the tetrapeptide Ac-Pro-Ala-Pro-Tyr-OH was characterized by computational methods. Search of the conformational space was performed within the molecular mechanics framework using the AMBER4.0 force field with an effective dielectric constant of 80. Unique structures of the peptide were compared with its bioactive conformation for the protein Streptomyces griseus Protease A, as taken from the crystal structure of the enzyme-peptide complex. The results show that the bound conformation is close to one of the unique conformations characterized in the conformational search of the isolated peptide. Moreover, the lowest energy minimum characterized in the conformational search exhibits large deviations when compared to the bound conformation of the crystal structure.     

Flexible-geometry conformational energy maps for the amino acid residue preceding a proline.

Previously calculated conformational energy maps suggest that the alpha-helical conformation for the residue preceding a proline is disfavored relative to the extended conformation by more than 7 kcal/mol. In known protein structures this conformation is observed, however, to occur for about 9% of all prolines. In addition, introduction or removal of prolines at theoretically unfavorable positions in proteins and peptides can have modest effects on stability and structure. To investigate the discrepancy between calculation and experiment, we have determined how the conformation of the proline affects the calculated energy. We have also explored the effect of bond length and bond angle relaxation on the conformational energy map. The conformational energy of the preceding residue is found to be unaffected by the conformation of the proline, but the effect of allowing covalent bond relaxation is dramatic. If bond lengths and angles, and dihedral angles within the pyrrolidine ring, are allowed to relax, a calculated energy difference between the alpha and beta conformations of 1.1 kcal/mol is obtained, in reasonable agreement with experiment. The detailed shape of the calculated energy surface is also in excellent agreement with the observed conformational distributions in known protein structures.     

Prediction of protein--ligand interactions: the complex of porcine pancreatic elastase with a valine-derived benzoxazinone

Prior to availability of the crystal structure of the complex, we evaluated models of the complex between porcine pancreatic elastase and a t-Boc–Val-derived benzoxazinone inhibitor. Models of the noncovalent and covalent complex were generated using computer graphics and each model was subjected to energy minimization using molecular mechanics. After the crystal structure became available, we found that the model with the lowest energy was in good agreement with the crystal structure, except for the position of the His57 side chain. Permissible conformations of the inhibitor were based on information from x-ray crystal structures and an earlier conformational energy investigation of t-Boc–amino acids. We did not, however, limit ourselves to these conformations. The conformation of the inhibitor in the lowest energy model and crystal structure, was not similar to any of the minimum-energy conformations of t-Boc–amino acids. This suggests that limiting proposed binding modes only to the lowest energy conformations of a ligand (prior to binding) may sometimes unfairly bias the procedure.     

The multiple-minima problem in the conformational analysis of polypeptides. III. An Electrostatically Driven Monte Carlo Method: Tests on enkephalin

The three-dimensional conformation of Met-enkephalin, corresponding to the lowest minimum of the empirical potential energy function ECEPP/2 (empirical conformational energy program for peptides), has been determined using a new algorithm, viz. the Electrostatically Driven Monte Carlo Method. This methodology assumes that a polypeptide or protein molecule is driven toward the native structure by the combined action of electrostatic interactions and stochastic conformational changes associated with thermal movements. These features are included in the algorithm that produces a Monte Carlo search in the conformational hyperspace of the polypeptide, using electrostatic predictions and a random sampling technique to locate low-energy conformations. In addition, we have incorporated an alternative mechanism that allows the structure to escape from some conformational regions representing metastable local energy minima and even from regions of the conformational space with great stability. In 33 test calculations on Met-enkephalin, starting from arbitrary or completely random conformations, the structure corresponding to the global energy minimum was found inall the cases analyzed, with a relatively small search of the conformational space. Some of these starting conformations wereright orleft-handed -helices, characterized by good electrostatic interactions involving their backbone peptide dipoles; nevertheless, the procedure was able to convert such locally stable structures to the global-minimum conformation.On leave from the National University of San Luis, Faculty of Sciences and Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Instituto de Matemática Aplicada, San Luis, Ejército de los Andes 950, 5700 San Luis, Argentina.     

Molecular simulation uncovers the conformational space of the λ Cro dimer in solution

The significant variation among solved structures of the λ Cro dimer suggests its flexibility. However, contacts in the crystal lattice could have stabilized a conformation which is unrepresentative of its dominant solution form. Here we report on the conformational space of the Cro dimer in solution using replica exchange molecular dynamics in explicit solvent. The simulated ensemble shows remarkable correlation with available x-ray structures. Network analysis and a free energy surface reveal the predominance of closed and semi-open dimers, with a modest barrier separating these two states. The fully open conformation lies higher in free energy, indicating that it requires stabilization by DNA or crystal contacts. Most NMR models are found to be unstable conformations in solution. Intersubunit salt bridging between Arg⁴ and Glu⁵³ during simulation stabilizes closed conformations. Because a semi-open state is among the low-energy conformations sampled in simulation, we propose that Cro-DNA binding may not entail a large conformational change relative to the dominant dimer forms in solution.     

Conformational analysis of some phenethylamine molecules

A set of empirical potential functions (EPF), previously used in conformational energy calculations of polymers, was employed in the study of the conformational properties of a number of methyl-substituted phenethylmines, as well as phenylmethylamine, phenyl-n-propylamine, and 3,4,5-trimethoxyamphetamine. The conformational free energy was computed for each of these molecular species in four states: neutral charge-vacuo (I), neutral charge-aqueous solution (II), positive charge-vacuo (III), positive charge-aqueous solution (IV). The molecules generally adopt one of two stable conformations: a folded conformation with the amine chain perpendicular to the ring, and the amine group nearest to the ring; and an extended conformation with the amine chain perpendicular to the ring, and the amine group far from the ring. The folded conformation is usually preferred for states I, II and III, while the extended form is adopted for state IV. By using empirical potential functions it was also possible to calculate the conformational entropies associated with the minimum energy conformations, thereby allowing the Boltzmann probabilities to be determined. These probabilities are a measure of the population density of each of the various low energy regions. Some of the molecules studied have a steric “bulge” below the plane of the benzene ring. All of the compounds studied which possess this “bulge” are psychotropically inactive, and, in most cases, also pharmacologically inactive. All active compounds studied do not possess this “bulge”.     

A competition smFRET assay to study ligand-induced conformational changes of the dengue virus protease

Ligand binding to proteins often is accompanied by conformational transitions. Here, we describe a competition assay based on single molecule Förster resonance energy transfer (smFRET) to investigate the ligand-induced conformational changes of the dengue virus (DENV) NS2B-NS3 protease, which can adopt at least two different conformations. First, a competitive ligand was used to stabilize the closed conformation of the protease. Subsequent addition of the allosteric inhibitor reduced the fraction of the closed conformation and simultaneously increased the fraction of the open conformation, demonstrating that the allosteric inhibitor stabilizes the open conformation. In addition, the proportions of open and closed conformations at different concentrations of the allosteric inhibitor were used to determine its binding affinity to the protease. The K_D value observed is in accordance with the IC₅₀ determined in the fluorometric assay. Our novel approach appears to be a valuable tool to study conformational transitions of other proteases and enzymes.     

Nmr studies of the molecular conformations in the linear oligopeptides H-(L-Ala)n-L-Pro-OH

The molecular conformations of the linear oligopeptides H-(L -Ala)_n-L -Pro-OH, with n = 1,2 and 3, have been investigated. ¹³C nmr observation of the equilibrium between the cis and trans forms of the Ala-Pro peptide bond indicated the occurrence of nonrandom conformations in solutions of these flexible peptides. The formation of the nonrandom species containing the cis form of the Ala-Pro bond was found to depend on the deprotonation of the carboxylic acid group of proline, the solvent, and the ionic strength in aqueous solution. The influence of intramolecular hydrogen bonding on the relative conformational energies of the species containing the cis and trans Ala-Pro peptide bond was studied by comparison of the peptides H-(Ala)_n-Pro-OH with analogous molecules where hydrogen bond formation was excluded by the covalent structure. In earlier work a hydrogen bond between the protonated terminal carboxylic acid group and the carbonyl oxygen of the penultimate amino acid residue had been suggested to stabilize conformations including trans proline. For the systems described here this hypothesis can be ruled out, since the cis:trans ratio is identical for molecules with methyl ester protected and free protonated terminal carboxylic acid groups of proline. Direct evidence for hydrogen bond formation between the deprotonated terminal carboxylic acid group and the amide proton of the penultimate amino acid residue in the molecular species containing cis proline was obtained from ¹H nmr studies. However, the cis:trans ratio of the Ala-Pro bond was not affected by N-methylation of the penultimate amino acid residue, which prevents formation of this hydrogen bond. Overall the experimental observations lead to the conclusion that the relative energies of the peptide conformations including cis or trans proline are mainly determined by intramolecular electrostatic interactions, whereas in the molecules considered, intramolecular hydrogen bonding is a consequence of specific peptide backbone conformations rather than a cause for the occurrence of energetically favored species. Independent support for this conclusion was obtained from model consideration which indicated that electrostatic interactions between the terminal carboxylic acid group and the carbonyl oxygen of the penultimate amino acid residue could indeed account for the observed relative conformational energies of the species containing cis and trans proline, respectively.     

Influence of local interactions on protein structure. I. Conformational energy studies of N-acetyl-N′-methylamides of pro-X and X-pro dipeptides

Conformational energy calculations using an Empirical Conformational Energy Program for Peptides (ECEPP) were carried out on the N-acetyl-N′-methylamides of Pro-X, where X = Ala, Asn, Asp, Gly, Leu, Phe, Ser, and Val, and of X-Pro, where X = Ala, Asn, Gly, and Pro. The conformational energy was minimized from starting conformations which included all combinations of low-energy single-residue minima and several standard bend structures. It was found that almost all resulting minima are combinations of low-energy single-residue minima, suggesting that intra residue interactions predominate in determining conformation. The calculations also indicate, however, that inter residue interactions can be important. In addition, librational entropy was found to influence the relative stabilities of some minima. Because of the existence of 10–100 low-energy minima for each dipeptide, the normalized statistical weight of an individual minimum rarely exceeds 0.3, suggesting that these dipeptides have considerable conformational flexibility and exist as statistical ensembles of low-energy structures. The propensity of each dipeptide to form bend conformations was calculated, and the results were compared with available experimental data. It was found that bends are favored in Pro-X dipeptides because ?_Pro is fixed by the pyrrolidine ring in a conformation which is frequently found in bends, but that bends are not favored in X-Pro dipeptides because interactions between the X residue and the pyrrolidine ring restrict the X residue to conformations which are not usually found in bends.