Synthesis and study of peptides with semirigid i and i + 7 side-chain bridges designed for alpha-helix stabilization

A search for conformational constraints on the peptide alpha-helical conformation indicated that para-substituted amino acid derivatives of a benzene ring might be suitable for linking pairs of side chains that are separated by two turns of the helix. A 14-residue synthetic, amphiphilic alpha-helical peptide model system has been used to study the helix stabilizing effects of a series of four such bridges having constitutionally isomeric structures. These bridges were used to link positions 3 and 10 of the model peptides. The peptides were synthesized in good yield by standard solid-phase methods, including cyclization on the solid support. They were then studied for their solution conformations and melting behavior by circular dichroism (CD) spectropolarimetry, and for their elution behavior on reversed-phase HPLC columns. In aqueous solution and in 50% (v/v) trifluoroethanol, the most effective bridge for helix stabilization consisted of a 4-(aminomethyl)phenylacetic acid residue (AMPA) linked by amide bonds to the side chain functional groups of a (S)-2,3-diaminopropionic acid residue (Dap) in position 3 of the model peptide and an aspartic acid residue in position 10. This Dap3(AMPA), Asp10 bridge was about as effective as two Lys(i), Asp(i+4) lactam bridges incorporated linking residues 3 and 7, and 10 and 14, in the same model peptide sequence. This suggests that it is worth about 1 kcal/mol of helix stabilization energy.     

Conformational preferences of amino acid side chains in collagen

Conformational energy computations were carried out on collagenlike triple-stranded conformations of several poly(tripeptide)s with the general structure CH₃CO? (Gly? X? Y)₃? NHCH₃. The sequences considered had various amino acid residues in position X or Y of the central tripeptide, with either Pro or Ala as a neighbor, i.e., Gly-X-Pro, Gly-X-Ala, Gly-Pro-Y, and Gly-Ala-Y. Minimum-energy conformations were computed for the side chains, and their distributions were compared for the four sequences. The residues used were Abu (= α-aminobutyric acid), Leu, Phe, Ser, Asp, Asn, Val, Ile, and Thr. The conformational energy of a ? Ch₂? CH₃ side chain in Abu was mapped as a function of the dihedral angle χ¹. Intrastrand interactions with neighboring residues do not affect the conformations of a side chain in position Y, and they have a minor effect on it in the X-Ala sequence, but they strongly restrict the conformational freedom of the side chain in the X-Pro sequence. Conversely, interstrand interactions do not affect side chains in position X, but they strongly restrict the conformational freedom of a side chain in position Y if there is a nearby Pro residue in a neighboring strand. Hydrogen bonds with the backbone can be formed in some conformations of long polar side chains, such as Asp, Asn, or Gln. All amino acid residues can be accommodated in collagen. Because of the interactions mentioned above, steric and energetic constraints can be correlated with observed preferences of certain amino acids for positions X or Y in collagen. Hence, these preferences may be explained, in part, in terms of differences in the conformational freedom of the side chains in the triple-stranded structure.     

Effect of lysine side chain length on intra-helical glutamate--lysine ion pairing interactions

Ion-pairing interactions are important for protein stabilization. Despite the apparent electrostatic nature of these interactions, natural positively charged amino acids Lys and Arg have multiple methylenes linking the charged functionality to the backbone. Interestingly, the amino acids Lys and Orn have positively charged side chains that differ by only one methylene. However, only Lys is encoded and incorporated into proteins. To investigate the effect of side chain length of Lys on ion-pairing interactions, a series of 12 monomeric alpha-helical peptides containing potential Glu-Xaa (i, i+3), (i, i+4) and (i, i+5) (Xaa = Lys, Orn, Dab, Dap) interactions were studied by circular dichroism (CD) spectroscopy at pH 7 and 2. At pH 7, no Glu-Xaa (i, i+5) interaction was observed, regardless of the Xaa side chain length. Furthermore, only Lys was capable of supporting Glu-Xaa (i, i+3) interactions, whereas any Xaa side chain length supported Glu-Xaa (i, i+4) interactions. Side chain conformational analysis by molecular mechanics calculations showed that the side chain length of Lys enables the Glu-Xaa (i, i+3) interaction with lower energy conformations compared to residues with side chain lengths shorter than that of Lys. Furthermore, these calculated low energy conformers were consistent with conformations of intra-helical Glu-Lys salt bridges in a non-redundant protein structure database. Importantly, the CD spectra for peptides with Glu-Lys interactions did not alter significantly upon changing the pH because of a greater contribution to these interactions by forces other than electrostatics. Incorporating side chains just one methylene shorter (Orn) resulted in significant pH dependence or lack of interaction, suggesting that nature has chosen Lys to form durable interactions with negatively charged functional groups.     

Peptide conformations--49(1): synthesis and structure-activity relationships of side chain modified peptides of cyclo(-D-Pro-Phe-Thr-Lys-Trp-Phe.)

Cyclic hexapeptide analogues representing the modified retro sequence of the amino acid residues 7-11 of natural somatostatin are known to protect liver cells from phalloidin poisoning. To determine the influence of steric, lipophilic, and charge effects on (a) the conformation of the backbone and the aromatic side chains and (b) the biological response, the side chains of Phe2, Lys4, and Phe6 of cyclo(-D-Pro1-Phe2-Thr3-Lys(Z)4-Trp5-Phe6-), 1a, one of the most active peptides found so far, were modified by various residues. The discussion of conformationally relevant parameters proves that neither backbone conformations nor populations of aromatic side chain rotamers were altered by these substitutions. The potency of these derivatives in a cytoprotection assay varies by at most one order of magnitude (more or less active than the parent peptide 1a). A qualitative evaluation of lipophilic, steric, and charge effects reveals the dominance of lipophilic effects of aromatic residues; the most potent compounds contain aromatic substructures in the side chain of Lys4.     

Electrostatic effects on the alpha-helix and beta-strand formation of BPTI(16-36) studied by Monte Carlo simulated annealing.

The electrostatic effects on the secondary structure forming tendencies of a peptide fragment with residues 16-36 of bovine pancreatic trypsin inhibitor, BPTI(16-36), are studied using Monte Carlo simulated annealing simulations. We consider three dielectric functions epsilon(r) of distance r: constant dielectric function (epsilon = 2; strong electrostatic interactions) and sigmoidal functions varying from epsilon(0) = 2 to epsilon(infinity) = 47 (intermediate) and to epsilon(infinity) = 78 (weak). Simulations with epsilon = 2 suggest that this peptide exhibits a significant propensity for beta-strand formations in accordance with a beta-sheet structure of the relevant segment in native BPTI. The tendency for alpha-helix formations becomes almost comparable with that of beta-strands in the simulation with epsilon(infinity) = 47, and there appears no appreciable conformational propensity for this case. Finally, the results with epsilon(infinity) = 78 generate low-energy conformations with conspicuous alpha-helices. These findings suggest the possibility that the change in electrostatic interactions can be the key factor for the conformational transitions of peptides between alpha-helix and beta-sheet that have recently been observed in experiments. These changes in electrostatic interactions can arise from those in various environmental factors such as conformations of the rest of the protein molecule and solvent conditions.     

Highly potent side-chain to side-chain cyclized enkephalin analogues containing a carbonyl bridge: synthesis, biology and conformation.

Six novel cyclic enkephalin analogues have been synthesized. Cyclization of the linear peptides containing basic amino acid residues in position 2 and 5 was achieved by treatment with bis(4-nitrophenyl)carbonate. It was found that some of the compounds exibit unusually high mu-opioid activity in the guinea pig ileum (GPI) assay. The 18-membered analogue cyclo(N(epsilon),N(beta)-carbonyl-D-Lys2,Dap5)-enkephalinamide turned out to be one of the most potent mu-agonists reported so far. NMR spectra of the peptides were recorded and structural parameters were determined. The conformational space was exhaustively examined for each of them using the electrostatically driven Monte Carlo method. Each peptide was finally described as an ensemble of conformations. A model of the bioactive conformation of this class of opioid peptides was proposed.     

Potent side-chain to side-chain cyclized dermorphin analogues containing a carbonyl bridge.

A new family of cyclic opioid peptide analogues related to the 1-4 sequence of dermorphin/deltorphin (Tyr-D-Aaa2-Phe-Aaa4-NH2) has been synthesized. The synthesis of the linear precursor peptides was accomplished by the solid-phase method and ring formation was achieved via a ureido group incorporating the side chain amino functions of D-Aaa2 (D-Lys, D-Orn) and Aaa4 (Lys, Orn, Dab, Dap). The peptides were tested in the guinea-pig ileum (GPI) and mouse vas deferens (MVD) assays. Most showed very high agonist potency in the GPI assay. The peptide containing D-Lys in position 2 and Dab in position 4 was 210 times more active than enkephalin, and that containing Orn and Dab, respectively, was 150 times more active than enkephalin. The latter peptide was also very active in the MVD assay, and showed an IC50 MVD/GPI ratio of 0.816. NMR spectra of selected peptides were recorded, and structural parameters were determined. The conformational space of the peptides was examined using the electrostatically driven Monte Carlo method. With the help of the NMR spectra each peptide was described as an ensemble of conformations. The conformations have been interpreted with regard to the opioid activities, and comparisons have been made with a model proposed earlier for enkephalin analogues.     

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.     

Calculations of the conformations of Iturin A in relation with NMR studies

Iturin A is an antifungal antibiotic which was isolated from a strain of Bacillus subtilis, and contains a lipophilic beta amino acid closing an heptapeptide cycle with polar L and D residues. Iturin A belongs to a lipopeptide family of which the LDDLLDL sequence is kept constant. NMR spectroscopy and semi-empirical energy calculations are combined to design the conformations of Iturin A in pyridine solution. J coupling constants and nOes (nuclear Overhauser enhancements) are used as guiding line for energy calculations. This preliminary study shows that Iturin A in pyridine appears as rather rigid, especially in the L Pro 5-D Asn 6 region, probably involved in a beta turn. The polar side chains can form different networks of intramolecular hydrogen bonds. The Tyr side chain, relatively mobile, could be involved in interactions with an hydrophobic environment as the beta amino acid side chain found away from the peptide cycle.     

Synthesis of novel adamantylacetyl derivative of peptidoglycan monomer--biological evaluation of immunomodulatory peptidoglycan monomer and respective derivatives with lipophilic substituents on amino group

Novel synthetic analogue of immunomodulatory peptidoglycan monomer 1 (PGM), (adamant-1-yl)-CH2CO-PGM (2), was prepared by acylation of epsilon-amino group of diaminopimelic acid with symmetrical (adamant-1-yl)-acetic acid anhydride in the presence of triethylamine. The product was isolated by gel filtration on Sephadex G-25, followed by ion exchange chromatography on SP-Sephadex C-25. The susceptibility of (adamant-1-yl)-CH2CO-PGM to hydrolysis with N-acetylmuramyl-L-alanine amidase was demonstrated, and the product of hydrolysis, (adamant-1-yl) CH2CO-pentapeptide 3, was characterized. Both 2 and 3 are water soluble and non-pyrogenic compounds. Immunomodulatory activity of PGM (adamant-1-yl)-CH2CO-PGM and structurally related derivative Boc-Tyr-PGM was compared in experiments in vivo, in mice, using ovalbumin (OVA) as an antigen. All three tested compounds exhibited comparable immunostimulating effects with respect to the induction of anti-ovalbumin immunoglobulin G. The results of evaluation of biological activity show that the substitution of free amino group in the parent peptidoglycan molecule with bulky lipophilic substituents did not affect the susceptibility to hydrolysis with N-acetylmuramyl-L-alanine amidase and did not alter markedly the immunostimulating activity. The results also indicate that the free amino group in the peptide chain is not a necessary requirement in the mechanism of immunostimulation of tested immunomodulators.     

Analysis of side-chain conformational distributions in neutrophil peptide-5 NMR structures

The side-chain conformations have been analyzed in the antimicrobial peptide, Neutrophil Peptide-5 (NP-5), whose structure was independently generated from nmr-derived distance constraints using a distance geometry algorithm. The side-chain and peptide dihedral angle distributions in the nmr structures were compared with those constructed from a data base of high-resolution protein crystal structures. The side-chain conformational preferences for NP-5 in solution are significantly different from those observed in the crystal structure data base. These results indicate that the side-chain conformations are quite disordered for many of the residues of NP-5. The absence of a correlation between the width of the conformational distribution and surface accessibility suggests that the disorder may be due to limitations in the structural information extracted from the nmr data rather than to molecular motion. However, it is also observed that the degree of conformational disorder is only weakly correlated with the number of nuclear Overhauser enhancements to a given side chain. Possible reasons for this are discussed. Molecular mechanics refinement of these structures did not significantly change the side-chain populations. Anomolously wide distributions are observed for rotations about the peptide bonds and the disulfide bonds in the NP-5 distance geometry structures, which are improved by the refinement. The very high degree of order observed for the central dihedral angle of the disulfide bond in the high-resolution crystal data base suggests that the rotation about this bond in proteins is determined by the local potential.     

Spatial structure of Thr-Pro-Ala-Glu-Asp-Phe-Met-Arg-Phe-NH2 molecule

The spatial structure of cardioactive Thr-Pro-Ala-Glu-Asp-Phe-Met-Arg-Phe-NH₂ molecule has been investigated using a theoretical conformational analysis. The low-energy conformations of the molecule were found, the values of the backbone and side T-T chain dihedral angles of amino acid residues constituting the peptide were determined, and the energies of intra- and interresidual interactions were estimated. It is revealed that the spatial structure of this molecule can exist in 11 stable backbone forms.     

Tertiary structure of myelopeptides. II. Conformational analysis of Phe-Arg-Pro-Arg-Ile-Met-Thr-Pro, Val-Val-Tyr-Pro-Asp, and Val-Asp-Pro-Pro

Theoretical conformational analysis was used to study the spatial structure and conformational properties of myelopeptides, bone marrow peptide mediators. The low-energy conformations of myelopeptides MP-4 (Phe-Arg-Pro-Arg-Ile-Met-Thr-Pro), MP-5 (Val-Val-Tyr-Pro-Asp), and MP-6 (Val-Asp-Pro-Pro) were found; the values of dihedral angles of backbone and side chains of the amino acid residues were determined; and the energies of intra- and interresidual interactions were estimated.     

Spatial structure of myelopeptides: I. Conformational analysis of MP-1, MP-2, and MP-3

Theoretical conformational analysis was used to study the spatial structure and conformational properties of myelopeptides, bone-marrow peptide mediators. The low-energy conformations of three hexapeptides MP-1 (Phe-Leu-Gly-Phe-Pro-Thr), MP-2 (Leu-Val-Val-Tyr-Pro-Trp), and MP-3 (Leu-Val-Cys-Tyr-Pro-Gln) were found, the values of dihedral angles of the backbone and side chains of the amino acid residues constituting these peptides were determined, and the energies of intra- and interresidual interactions were estimated.     

A tree-based algorithm for determining the effects of solvation on the structure of salivary gland tripeptide NH3+-D-PHE-D-GLU-GLY-COO-

A D-enantiomeric analog of the submandibular gland rat-1 tripeptide FEG (Seq: NH(3)(+)-Phe-Glu-Gly-COO(-)) called feG (Seq: NH(3)(+)-D-Phe-D-Glu-Gly-COO(-)) was examined by molecular dynamics simulations in water. Previous in vacuo simulations suggested a conformation consisting predominantly of interactions between the Phe side chain and glutamyl-carboxyl group and a carboxyl/amino termini interaction. The solvated peptide was simulated using two approaches which were compared-a single 400-ns simulation and a "simulation tree." The "tree" approach utilized 45 10-ns simulations with different conformations used as initial structures for given trajectories. We demonstrate that multiple short duration simulations are able to describe the same conformational space as that described by longer simulations. Furthermore, previously described in vacuo interactions were confirmed with amendments: the previously described head-to-tail arrangement of the amino and carboxyl termini, was not observed; the interaction between the glutamyl carboxyl and Phe side chain describes only one of a continuum of conformations present wherein the aromatic residue remains in close proximity to the glutamyl carbonyl group, and also interacts with either of the two available carboxyl groups. Finally, utilizing only two separate 10-ns trajectories, we were able to better describe the conformational space than a single 60-ns trajectory, realizing a threefold decrease in the computational complexity of the problem.     

Influence of hydrogen bonding,main chain–side chain interactions,and protecting groups on the excimer formation of bis(pyrenylalanine) peptides

Dipeptides of the aromatic fluorescent amino acid, pyrenylalanine, are studied using both stationary and transient fluorescence techniques. Since the conformational transitions of the peptide chain are slow compared to the decay of the pyrene excited state, both ground state conformations, adopted by the peptide, i.e., C₅ and C₇, can be monitored separately. Kinetic models are proposed to describe the molecular dynamics of the peptide chain as probed by the intramolecular excimer formation between both pyrene chromophores. These kinetic schemes explain the influence of solvent, chain chirality, main chain–side chain interactions, and nature of the protecting groups on the emission spectrum and the fluorescence decay profile of these model peptides. These schemes also provide a tool to calculate rate constants of conformational transitions and excimer formation. By comparing the kinetic and thermodynamic parameters of the various compounds, the influence of a structural modification on the molecular dynamics of the peptide chain is determined.     

Statistical and energetic analysis of side-chain conformations in oligopeptides

The distributions of side-chain conformations in 258 crystal structures of oligopeptides have been analyzed. The sample contains 321 residues having side chains that extend beyond the C beta atom. Statistically observed preferences of side-chain dihedral angles are summarized and correlated with stereochemical and energetic constraints. The distributions are compared with observed distributions in proteins of known X-ray structures and with computed minimum-energy conformations of amino acid derivatives. The distributions are similar in all three sets of data, and they appear to be governed primarily by intraresidue interactions. In side chains with no beta-branching, the most important interactions that determine chi 1 are those between the C gamma H2 group and atoms of the neighboring peptide groups. As a result, the g- conformation (chi 1 congruent to -60 degrees) occurs most frequently for rotation around the C alpha-C beta bond in oligopeptides, followed by the t conformation (chi 1 congruent to 180 degrees), while the g+ conformation (chi 1 congruent to 60 degrees) is least favored. In residues with beta-branching, steric repulsions between the C gamma H2 or C gamma H3 groups and backbone atoms govern the distribution of chi 1. The extended (t) conformation is highly favored for rotation around the C beta-C gamma and C gamma-C delta bonds in unbranched side chains, because the t conformer has a lower energy than the g+ and g- conformers in hydrocarbon chains. This study of the observed side-chain conformations has led to a refinement of one of the energy parameters used in empirical conformational energy computations.     

Spatial structure of the cardioactive octapeptide

The spatial structure of the cardioactive octapeptide Pro1-Gln2-Asp3-Pro4-Phe5-Leu6-Arg7-Ile8-NH2 was investigated using the theoretical conformational analysis. The low-energy conformations of the octapeptide molecule were found, the values of dihedral angles of the backbone and side chains of the amino acid residues constituting the peptide were determined, and the energies of intra-and interresidual interactions were estimated. It was shown that the spatial structure of this molecule represent six stable low-energy forms of the main chain.     

Conformational studies of antimetastatic laminin-1 derived peptides in different solvent systems, using solution NMR spectroscopy.

Due to its critical role in cancer progression, interactions between laminin-1 and the 67 kDa Laminin-Binding Protein (the 67 kDa LBP) have been the focus of a number of structural and biological studies. As laminin-1 is such a large and complex molecule, research interests have turned to the investigation of bioactive peptides derived from binding domains of laminin-1. Two peptides of interest, CDPGYIGSR (peptide 11) and YIGSR, both derived from the beta1 chain of laminin-1, have been shown to block invasion of basement membranes by tumor cells. Substituting the C-terminal arginine to lysine, a conservative substitution, results in a loss of peptide antimetastatic activity. This difference in bioactivity has been attributed, based on numerous modeling studies of free peptide conformations, to structural differences between YIGSR and YIGSK. Yet the nature of the 'active' free peptide backbone conformation has been a matter of debate and controversy. In order to test the validity of the structural modeling claims, we have undertaken detailed conformational studies of the two laminin-1 derived peptides YIGSR and CDPGYIGSR along with the biologically inactive YIGSK analog by two-dimensional solution 1H NMR spectroscopy in three different solvent systems. Herein we report that although both the active (YIGSR, CDPGYIGSR) and the inactive (YIGSK) peptides can adopt several closely related conformations in solution, the two peptides share similar conformational preferences, and there are no significant structural differences between the active and inactive peptides, contrary to previously reported modeling data. We conclude that the basis of the peptide biological activity, in contrast to published models, cannot be attributed to well-defined structural preferences of the free peptides. We infer that the difference in bioactivity observed between YIGSR and YIGSK originates primarily from the chemical nature of the arginine versus lysine sidechain substitution, rather than being due to a structural change in the free peptide conformations.