Peptide inhibitors of E. collagenolyticum bacterial collagenase--effect of N-methylation. Consequences on biological activity and conformational properties.

Peptide inhibitors of E. collagenolyticum bacterial collagenase, HS-CH2-CH2-CO-Pro-Yaa (Yaa = Ala, Leu, Nle), have been N-methylated at the Yaa position. The N-methylation slightly increases the inhibitory potency of the modified peptides as compared to the parent compounds. The conformational effects of the N-methylation have been investigated by both 1H 2D-NMR and molecular mechanics energy minimization. Three low-energy conformers have been predicted for the unmethylated parent compounds (Yaa = Ala, Leu, Nle). They are characterized by the psi value of the central proline residue: psi Pro = 150 degrees (trans' conformation), psi Pro = 70 degrees (C7 conformation) and psi Pro = -50 degrees (cis' conformation). The N-methylation has been found to strongly increase the energy of the C7 conformer and to a less extent the energy of the cis' conformer. This leaves the trans' conformation as the only low-energy conformer. The ROESY experiments have established that both the N-methyl peptides and the parent compounds adopt the same preferred backbone conformation in water solution, i.e. the trans' conformation. Based on these results, the activities of the N-methyl peptides are discussed and a possible conformation of the inhibitor in the bound state is proposed.     

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.     

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.     

Influence of local interactions on protein structure. III. Conformational energy studies of N-acety-N′-methylamides of Gly-X and X-Gly dipeptides

Conformational energy calculations using an empirical conformational energy program for peptides (ECEPP) were carried out on 20 N-acetyl- N′-methylamides of Gly-X and X-Gly depeptides, where X = Ala, Asn, Asp, Gly, Phe, Ser, Thr, Tyr, Val, and Pro, and also of Leu-Gly. Each depeptde was found to have 25 or more low-energy minima, except Gly-Thr, which had only 11 low-energy minima because of the stable side chian-backbone hydrogen present in all low-energy conformation. As a group, the stble chain-backbone hydrogen bonds present in all low-energy conformations. As a group, the Gly-containing dipeptides were calculated in all low-energy prpensity for formation of bends than the Ala-containing depeptides. The X- Gly dipeptides were calculated to favor bends more than the Gly-X dipeptides, primarlly because of the high stability of the type II bend in X-Gly dipeptides. These results are in agreement with obseved occurrences of bends in the x-ray structures of globular proteins. The calculated conformation properties were found to be in good agreement with experimental results.     

The effect of the L-azetidine-2-carboxylic acid residue on protein conformation. I. Conformations of the residue and of dipeptides

The L-azetidine-2-carboxylic acid (Aze) residue can be incorporated into proteins in the place of L-proline, of which it is the lower homologue. This substitution alters the properties of proteins, especially of collagen. Conformational constraints in N-acetyl-Aze-N'-methylamide and in several dipeptides containing Aze have been analyzed by means of energy computations. They have been compared with peptides containing Pro. The overall conformational preferences of Aze and Pro are similar, but several significant differences occur between them. In general, peptides containing Aze are somewhat more flexible than corresponding peptides containing Pro, because of a decrease in constraints caused by repulsive nonconvalent interactions of the atoms of the ring with neighboring residues. This results in an entropic effect that lessens the stability of ordered polypeptide conformations with respect to the disordered statistical coil. The collagen-like near-extended conformation is energetically less favorable for Aze than for Pro in the single residue and in dipeptides. This effect also contributes to a destabilization of the collagen triple helix. The influence of Aze on the conformation of polypeptides is discussed in the accompanying papers.     

Conformational choice at alpha,alpha-di-n-propylglycine residues: helical or fully extended structures?

The conformational analysis of peptides containing a single alpha, alpha-di-n-propylglycine (Dpg) residue incorporated into valine-rich sequences has been undertaken in order to delineate the possible role of sequence effects in stabilizing fully extended (C(5)) or local helical conformations at this residue. The three peptides Boc-Val-Dpg-Val-OMe (3), Boc-Val-Val-Dpg-Val-OMe (4), Boc-Val-Val-Dpg-Val-Val-OMe (5), have been studied by (1)H-nmr methods in chloroform (CDCl(3)) and dimethylsulfoxide (DMSO) solutions. Even in a relatively poorly solvating medium like CDCl(3), all the valine NH groups appear to be solvent-exposed, suggesting an absence of folded beta-turn conformations. However, in both CDCl(3) and DMSO the Dpg NH groups in all the three peptides appear to behave like apparently solvent-inaccessible groups. In fully extended C(5) conformations, the proximity of the NH and CO groups of Dpg may preclude effective solvation due to a combination of stereoelectronic factors. Nuclear Overhauser effects provide support for the largely extended backbones. The crystal structure of peptide 3 reveals an extended conformation at Dpg (2) with straight phi = -176 degrees, psi = 180 degrees. A correlation between the crystallographically observed backbone conformation and solution nmr parameters in DMSO has been attempted using available data. Dpg residues placed in poor helix stabilizing environments may be expected to favor a local C(5) conformation.     

Structural compromise of disallowed conformations in peptide and protein structures

Using a data set of 454 crystal structures of peptides and 80 crystal structures of non-homologous proteins solved at ultra high resolution of 1.2 A or better we have analyzed the occurrence of disallowed Ramachandran (phi, psi) angles. Out of 1492 and 13508 non-glycyl residues in peptides and proteins respectively 12 and 76 residues in the two datasets adopt clearly disallowed combinations of Ramachandran angles. These examples include a number of conformational points which are far away from any of the allowed regions in the Ramachandran map. According to the Ramachandran map a given (phi, psi) combination is considered disallowed when two non-bonded atoms in a system of two-linked peptide units with ideal geometry are prohibitively proximal in space. However, analysis of the disallowed conformations in peptide and protein structures reveals that none of the observations of disallowed conformations in the crystal structures correspond to a short contact between non-bonded atoms. A further analysis of deviations of bond lengths and angles, from the ideal peptide geometry, at the residue positions of disallowed conformations in the crystal structures suggest that individual bond lengths and angles are all within acceptable limits. Thus, it appears that the rare tolerance of disallowed conformations is possible by gentle and acceptable deviations in a number of bond lengths and angles, from ideal geometry, over a series of bonds resulting in a net gross effect of acceptable non-bonded inter-atomic distances.     

Interpretation of some Conformational Data in Myoglobin and Lysozyme

THE conformations of the polypeptide chains of myoglobin¹ and lysozyme² have been successfully simulated with the aid of computed Van der Waals contact and energy maps of the theoretical independent peptide unit (IPU)^3–5. The non-glycyl experimental points plotted on an alanyl IPU are rather scattered on the allowed conformational regions of the map⁶, especially in the case of lysozyme. By contrast, well defined clusters of points can be observed when only the amino-acid residues in segments of the helical secondary structure (mainly α and β chains) are plotted. In addition, clusters of points, albeit less well defined, can be observed by plotting the points relative to the experimental conformations of the first non-helical amino-acid residue next to a more or less folded segment of that α-helical type so frequently present in globular proteins (Fig. 1).     

Conformation dynamics of cyclic disulfides and selenosulfides in CXXC(U) (X = Gly,Ala) tetrapeptide redox motifs

Thioredoxin fold proteins often contain a Cys‐(Xxx)_n‐Cys(Sec) or CX_nC(U) motif, where the active cysteine (C) or selenocysteine (U) is bridged by X residues, which vary with protein function. The effect of the X residues on the conformation space of the oxidized disulfide and selenosulfide forms of the CXXC(U) motif has been investigated using molecular dynamics (MD) and density functional theory. Multi‐microsecond‐length MD simulations of the CGGC, CGAC, and CAGC cyclic peptides show that CGGC rings readily exchange between several conformations over the course of the simulation, but steric interactions with the methyl group of Ala limit the conformation space available to the cyclic peptide, especially for CGAC. The potential for the motif to be reduced, as measured by the energy of the lowest unoccupied molecular orbitals, is dependent upon the ring conformation. These results suggest that control of available conformations by the bridging residues and the protein tertiary structure may be important for defining the function of the CXXC motif. Theoretical ⁷⁷Se chemical shifts of the selenosulfide moiety are dependent upon the conformation and/or intramolecular Se···O interactions with the backbone carbonyl group of the C‐terminal U residue.     

Amino acid determinants of beta-hairpin conformation in erythropoeitin receptor agonist peptides derived from a phage display library

Display of peptide libraries on filamentous phage has led to the identification of peptides of the form X(2-5)CX(2)GPXTWXCX(2-5) (where X is a variable residue) that bind to the extra-cellular portion of the erythropoietin receptor (EPO-R). These peptides adopt beta-hairpin conformations when co-crystallized with EPO-R. Solution NMR studies reveal that the peptide is conformationally heterogeneous in the absence of receptor due to cis-trans isomerization about the Gly-Pro peptide bond. Replacement of the conserved threonine residue with glycine at the turn i+3 position produces a stable beta-hairpin conformation in solution, although this peptide no longer has activity in an EPO-R-dependent cell proliferation assay. A truncated form of the EPO-R-binding peptide (containing the i+3 glycine residue) also forms a highly populated, monomeric beta-hairpin. In contrast, phage-derived peptide antagonists of insulin-like growth factor binding protein 1 (IGFBP-1) have a high level of sequence identity with the truncated EPO-R peptide (eight of 12 residues) yet adopt a turn-alpha-helix conformation in solution. Peptides containing all possible pairwise amino acid substitutions between the EPO-R and IGFBP-1 peptides have been analyzed to assess the degree to which the non-conserved residues stabilize the hairpin or helix conformation. All four residues present in the original sequence are required for maximum population of either the beta-hairpin or alpha-helix conformation, although some substitutions have a more dominant effect. The results demonstrate that, within a given sequence, the observed conformation can be dictated by a small subset of the residues (in this case four out of 12).     

Conformational preferences and prolyl cis–trans isomerization of phosphorylated Ser/Thr‐Pro motifs

The conformational study on Ac‐pSer‐Pro‐NHMe and Ac‐pThr‐Pro‐NHMe peptides has been carried out using hybrid density functional methods with the implicit solvation reaction field theory at the B3LYP/ 6‐311++G(d,p)//B3LYP/6‐31+G(d) level of theory in the gas phase and in solution (chloroform and water). For both pSer‐Pro and pThr‐Pro peptides in the gas phase and in chloroform, the most preferred conformation has the α‐helical structure for the pSer/pThr residue, the down‐puckered polyproline I structure for the Pro residue, and the cis prolyl peptide bond between the two residues, in which two hydrogen bonds between the phosphate oxygens with the backbone N? H groups seem to play a role. However, the trans conformations that have a single hydrogen bond of the phosphate oxygen with either of two backbone N? H groups become most preferred for both peptides in water. This is because the hydration free energy of the anionic oxygen of the phosphate group is expected to dramatically decrease for the cis conformation upon formation of the hydrogen bond with the backbone N? H groups. These calculated results are consistent with the observations by NMR and IR experiments, suggesting the existence of hydrogen bonds between the charged phosphoryl group and the backbone amide protons in solution. The calculated cis populations of 14.7 and 14.2% and rotational barriers of 19.87 and 20.57 kcal/mol to the cis‐to‐trans isomerization for pSer‐Pro and pThr‐Pro peptides in water, respectively, are consistent with the observed values for pSer‐Pro and pThr‐Pro containing peptides from NMR experiments. However, the hydrogen bond between the prolyl nitrogen and the following amide N? H group, which was suggested to be capable of catalyzing the prolyl isomerization, does not play a role in stabilizing the preferred transition state for the pSer/pThr‐Pro peptides in water. Instead, the amide hydrogen of the NHMe group is involved in a bifurcated hydrogen bond with the anionic oxygen and phosphoester oxygen of the phosphate group. © 2009 Wiley Periodicals, Inc. Biopolymers 93: 330–339, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Context-dependent conformation of diethylglycine residues in peptides.

Diethylglycine (Deg) residues incorporated into peptides can stabilize fully extended (C5) or helical conformations. The conformations of three tetrapeptides Boc-Xxx-Deg-Xxx-Deg-OMe (Xxx=Gly, GD4; Leu, LD4 and Pro, PD4) have been investigated by NMR. In the Gly and Leu peptides, NOE data suggest that the local conformations at the Deg residues are fully extended. Low temperature coefficients for the Deg(2) and Deg(4) NH groups are consistent with their inaccessibility to solvent, in a C5 conformation. NMR evidence supports a folded beta-turn conformation involving Deg(2)-Gly(3), stabilized by a 4-->1 intramolecular hydrogen bond between Pro(1) CO and Deg(4) NH in the proline containing peptide (PD4). The crystal structure of GD4 reveals a hydrated multiple turn conformation with Gly(1)-Deg(2) adopting a distorted type II/II' conformation, while the Deg(2)-Pro(3) segment adopts a type III/III' structure. A lone water molecule is inserted into the potential 4-->1 hydrogen bond of the Gly(1)-Deg(2) beta-turn.     

Monte Carlo calculations on the conformations of models for the glycopeptide linkage of glycoproteins

In order to search for probable conformations of the peptide, the amino acid side chain, and the carbohydrate linkage in glycoproteins, conformational energy surfaces of glycopeptide model compounds were studied by Monte Carlo methods using the Metropolis algorithm. The potential energies were composed of empirical energy functions which include nonbonded interactions, electrostatics, hydrogen bonding, and torsional energies specified by parameters which have been used for peptides. Calculations were performed on 1-N-acetyl-2-acetamido-beta-D-glucopyranosyl amine and the glycosylated dipeptide N-acetyl-delta-N-(2-acetamido-beta-D-glucopyranosyl)-L-asparaginyl-N'-methyl amide as models for N-glycosylated peptides and on methyl-2-acetamido-alpha-D-galactopyranoside as well as the glycosylated dipeptides N-acetyl-gamma-O-(2-acetamido-alpha-D-galactopyranosyl)-L-threonyl-N'-methyl amide and its seryl analog as models for O-glycosylated glycoproteins. The probable conformations of these compounds were analyzed by single-angle probability tables and by two-dimensional conformation density maps projected from the Markov chains which contained up to six independently varied conformational dihedral angles. The presence of high barriers to rotation required the use of search strategies which resulted in a rather low acceptance rate for new conformations in the Metropolis algorithm in order to avoid trapping of the Markov chain in local energy minima. This problem contributed to the failure of these calculations to attain complete convergence to the thermodynamic limit for the glycosylated dipeptide models in which six dihedral angles were independently varied. Analysis of the results shows that the conformational space available to the highly crowded axial glycosides of the alpha-O-GalNAc type is much more restricted than that for the N-asparaginyl glycopeptides. The most probable conformation for the O-glycosylated peptides is is a beta-turn while N-glycosylated peptides may be either in a beta-turn or an extended conformation.     

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.     

Applications of simulated annealing to peptides

We report the application of a new conformation searching algorithm called simulated annealing to the location of the global minimum energy conformation of peptides. Simulated annealing is a Metropolis Monte Carlo approach to conformation generation in which both the energy and temperature dependence of the Boltzmann distribution guides the search for the global minimum. Both uphill and downhill moves are possible, which allows the molecule to escape from local minima. Applications to the 20 natural amino acid "dipeptide models" as well as to polyalanines up to Ala80 are very successful in finding the lowest energy conformation. A history file of the simulated annealing process allows reconstruction and examination of the random walk around conformation space. A separate program, Conf-Gen, reads the history file and extracts all low-energy conformations visited during the run.     

CD-n.m.r. study of the solution conformation of bradykinin analogs containing alpha-aminoisobutyric acid

The conformation in aqueous solution of several alpha-aminoisobutyric acid (AIB)-containing analogs of bradykinin (BK) has been probed by complementary CD and 1H n.m.r. measurements. The conclusion reached is that substitution of AIB for Pro2 and/or Pro3 in BK stabilizes a degree of beta-turn conformation in the N-terminal tetrapeptide moiety of the resulting analogs. Changing the solvent from water to DMSO or TFE further enhances the contribution of particular hydrogen bonded structures to the time-averaged conformation of these peptides. Bradykinin and [AIB7]-BK adopt similar hydrogen bonded conformations in TFE, apparently with a contribution from a beta-turn involving their common Arg1-Pro2-Pro3-Gly4 moiety. The contrasting biological activities of BK and its AIB-analogs are considered in terms of the conformational analogy between the AIB-residue and cis' Pro and the propensity for a beta-turn at the N-terminus of the peptide.     

Use of molecular mechanics for secondary structure prediction. Is it possible to reveal alpha-helix?

A new approach to predicting protein standard conformations is suggested. The idea consists in modeling by molecular mechanics tools a continuous alpha-helical conformation for the whole protein. The profile of energy along the model alpha-helix reveals minima corresponding to real alpha-helical segments in the native protein. The 3/10-helices and beta-turns including a local alpha-helical conformation may be detected as well. All alpha-helical segments in the test sample are delineated; mean residue by residue accuracy Q(3alpha) is 79%. This non-statistical approach can shed light on the physical grounds of alpha-helix formation.     

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.     

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.