Computer simulation of the conformational properties of retro-inverso peptides. I. Empirical force field calculations of rigid and flexible geometries of N-acetylglycine-N'-methylamide, bis(acetamido) methane, and N,N'-dimethylmalonamide and their corresponding C alpha-methylated analogs

Oxalyl retro-peptide gelators. Synthesis, gelation properties and stereochemical effects

In this work we report on gelation properties, self-assembly motifs, chirality effects and morphological characteristics of gels formed by chiral retro-dipeptidic gelators in the form of terminal diacids (1a-5a) and their dimethyl ester (1b-5b) and dicarboxamide (1c-5c) derivatives. Terminal free acid retro-dipeptides (S,S)-bis(LeuLeu) 1a, (S,S)-bis(PhgPhg) 3a and (S,S)-bis(PhePhe) 5a showed moderate to excellent gelation of highly polar water/DMSO and water/DMF solvent mixtures. Retro-peptides incorporating different amino acids (S,S)-(LeuPhg) 2a and (S,S)-(PhgLeu) 4a showed no or very weak gelation. Different gelation effectiveness was found for racemic and single enantiomer gelators. The heterochiral (S,R)-1c diastereoisomer is capable of immobilizing up to 10 and 4 times larger volumes of dichloromethane/DMSO and toluene/DMSO solvent mixtures compared to homochiral (S,S)-1c. Based on the results of (1)H NMR, FTIR, CD investigations, molecular modeling and XRPD studies of diasteroisomeric diesters (S,S)-1b/(S,R)-1b and diacids (S,S)-1b/(S,R)-1a, a basic packing model in their gel aggregates is proposed. The intermolecular hydrogen bonding between extended gelator molecules utilizing both, the oxalamide and peptidic units and layered organization were identified as the most likely motifs appearing in the gel aggregates. Molecular modeling studies of (S,S)-1a/(S,R)-1a and (S,S)-1b/(S,R)-1b diasteroisomeric pairs revealed a decisive stereochemical influence yielding distinctly different low energy conformations: those of (S,R)-diastereoisomers with lipophilic i-Bu groups and polar carboxylic acid or ester groups located on the opposite sides of the oxalamide plane resembling bola amphiphilic structures and those of (S,S)-diasteroisomers possessing the same groups located at both sides of the oxalamide plane. Such conformational characteristics were found to strongly influence both, gelator effectiveness and morphological characteristics of gel aggregates.

Force-field parametrization of retro-inverso modified residues: Development of torsional and electrostatic parameters

Torsional and the electrostatic parameters for molecular mechanics studies of retro-inverso modified peptides have been developed using quantum mechanical calculations. The resulting parameters have been compared with those calculated for conventional peptides. Rotational profiles, which were obtained spanning the corresponding dihedral angle, were corrected by removing the energy contributions associated to changes in interactions different from torsion under study. For this purpose, the torsional energy associated to each point of the profiles was estimated as the corresponding quantum mechanical energy minus the bonding and nonbonding energy contributions produced by the perturbations that the variation of the spanned dihedral angle causes in the bond distances, bond angles and the other dihedral angles. These energies were calculated using force-field expressions. The corrected profiles were fitted to a three-term Fourier expansion to derive the torsional parameters. Atomic charges for retro-inverso modified residues were derived from the rigorously calculated quantum mechanical electrostatic potential. Furthermore, the reliability of electrostatic models based on geometry-dependent charges and fixed charges has been examined.

Interactions of 5-deazapteridine derivatives with Mycobacterium tuberculosis and with human dihydrofolate reductases

There are major differences between the structures of human dihydrofolate reductase (hDHFR) and Mycobacterium tuberculosis dihydrofolate reductase (mtDHFR). These differences may allow us to design more selective mtDHFR inhibitors. In this paper we study the reactions of six different compounds derived from 5-deazapteridine with human and bacterial enzymes. Results suggest that the addition of hydrophobic groups to the aminophenyl ring would increase mtDHFR-inhibitor affinity and selectivity.

Extended form of a retro-inverso peptide stabilized by beta-sheet unidirectional H-bonds: Crystallographic and NMR evidence

The crystallographic investigation of the retro-inverso peptide Bz-S-gAla-R-mAla-NHPh reveals an extended backbone conformation where the NH groups of the gem-diamino alkyl moiety and the CO groups of the malonyl residue face side by side. This extended conformation, presenting all carbonyls on opposite sides of the NH groups, is stabilized by interstrand H-bonds running in a single direction of the parallel beta-sheets that characterize the crystal packing. These sheets differ from the beta-sheets formed by native amino acids only. (1)H-NMR nuclear Overhauser effect spectroscopy (NOESY) experiments suggest that a conformation similar to that found in the crystal also prevails in dimethylsulfoxide solution. Previous potential energy calculations of gem-diamino alkyl (g) and malonyl (m) Ala residues predicted that extended forms were less stable than the helical ones because of strong electrostatic repulsions between the parallel polar groups. Similar arguments were invoked to give more weight to helical forms of the retro-peptide units in the proposal of packing models of some nylons in their crystalline polar regions. The present findings show that both g and m Ala residues can experience the extended conformation in the beta-sheet aggregation. The energy increase occurring in one strand, due to the parallel orientation of consecutive peptide dipoles, is more than compensated by favorable cooperative interactions among head-to-tail aligned peptide dipoles of facing strands, resulting in the formation of two C==O...H==N H-bonds per residue.

Crystal structures of peptides and modified peptides

The X-ray diffraction experiments on peptides and related molecules which have been carried out in Western Europe, except Italy, in the last eight years are reviewed. The crystal structures of some bioactive peptides such as Leu-enkephalin (a neurotransmitter), cyclosporin A (an immunomodulator in both the free and protein-bound state), balhimycin (an antibiotic) and octreotide (a somatostatin analogue) are briefly presented. Crystallized N- and C-protected model peptides have given an insight into the folding tendency and folding modes depending on the peptide sequences. The crystal structures of various pseudopeptide molecules reveal how the three-dimensional structure of peptide analogues can be modulated by substituting non-peptide groups for the peptide bond. A few examples of structural mimetics of the beta- and gamma-turns, and of templates for alpha-helix induction are also presented.

Conformational studies of retro-inverso peptides: the crystal and molecular structure of the hydantoin from H-Ala-g-Ala-mGly-OBzl

The crystal structure of the hydantoin 1-[(S)-1'-aminoethylmalonyl benzyl ester]-(S)-4-methylimidazolidin-2,5-dione (1) derived from the peptide H-Ala-gAla-mGly-OBzl, having the retro-inverso modification of the Ala-Gly bond, has been determined by x-ray diffraction analysis. The crystals are orthorhombic, space group P2(1)2(1)2(1) with a = 6.539, b = 14.721, c = 17.101 A, z = 4. The structure was solved by direct methods and refined with anisotropic thermal factors to a final R value of 0.067 for the 947 observed reflections. Reversal of the Ala-Gly amide bond perturbs the folding tendency of the backbone shown by the parent peptide t-BuCO-Ala-Gly-NHiPr. The gem-diamino residue, gAla, and the malonyl moieties are found in the helical and the extended conformations, respectively. Intramolecular hydrogen bonding is not observed. The molecules in the crystal are held together by the formation of two intermolecular hydrogen bonds of the N-H ... O=C type with N ... O distances of 2.86 and 3.17 A, respectively.

Recent developments in retro peptides and proteins--an ongoing topochemical exploration

Main-chain peptidomimetics based on peptide-bond reversal and inversion of chirality represent important structural alterations for peptides and proteins, and are highly significant for biotechnology; these modifications have been widely applied: the D-HIV-protease dimer cleaves only all-D substrate; an all-D-hexapeptide opioid is able to produce analgesia following intraperitoneal administration. Antigenicity and immunogenicity can be achieved by metabolically stable antigens such as all-D- and retro-inverso-isomers of natural antigenic peptides. Isomers, including the retro- and retro-inverso- forms, of hybrid peptides derived from cercropin A and melittin, maintain antimicrobial activity. Therefore, an insight is provided into structure-activity relationships and the rational design of biologically important isomeric peptides.

Coulombic attractions between partially charged main-chain atoms stabilise the right-handed twist found in most beta-strands

The use of Lennard-Jones potentials gives rise to an expected energy distribution for main-chain polypeptide conformations in the Ramachandran plot that matches well the observed distribution of phi, psi values in high-resolution proteins. The position of the energy minimum in the beta-strand conformation region is situated where there is a substantial contribution from the electrostatic attraction between the partial charge of the carbonyl carbon atom of one amino acid residue and that of the carbonyl oxygen atom of an adjacent residue. This attraction gives rise to a preference for the right-twisted beta-strand conformation compared with the left-twisted conformation. The majority of beta-sheets are twisted, almost always in one direction. Looking along a single strand, the twist is to the right. This twist also helps provide a rationale for the characteristic topology of the strand-helix-strand unit often observed in alpha/beta proteins. The electrostatic explanation for the twist we propose has not, to our knowledge, been explicitly suggested previously. The factor that has been most widely proposed to explain the twist is steric hindrance involving side-chain atoms. We provide evidence that the electrostatic effect is of comparable significance. Right-twisted beta-strands are geometrically closely related to polyproline II helices and to collagen helices, both of which are left-handed. Short regions of polyproline II type helices, which are sometimes, but not always, rich in proline residues, are common at protein surfaces. We point out that these helices are stabilised by the same carbonyl-carbonyl interactions as in right-twisted beta-strands.

Conformational study of alanine and alpha-aminoisobutyric psi [NHCO]-retroamide peptide analogues

The minimum-energy conformations of N,N'-dimethyl-2-methylmalondiamide (2mMA) and N,N'-dimethyl-2-dimethylmalondiamide (2dmMA) have been computed using the AM1 method. The results show that molecules tend to form an intramolecular six-membered ring hydrogen-bond system. The same trend is observed in crystals of malonamide derivatives substituted in the central carbon atom. The results are also compared with previous molecular-mechanics calculations. The conformational trends of 2mMA and 2dmMA are very similar to those showed by N,N'-dimethylmalondiamide. Details of the computed conformational energy maps for 2mMA and 2dmMA compounds are given.

SCF-MO study of the polyhydration of N,N'-dimethylmalonamide

The effect of hydration on the conformation of the N,3 nylons has been studied using AM1 SCF-MO on the model molecule N,N'-dimethylmalonamide (NNMA). The polyhydration model has been used, so that only the water molecules directly bonded to solute molecule are considered. The study reveals that, although hydration induces dramatic changes in the torsional angles of the molecule, the new conformation still remains in the helical region of the Ramachandran map.

14-membered cyclic opioids related to dermorphin and their partially retro-inverso modified analogues. II. Preferred conformations in solution as studied by 1H-NMR spectroscopy

The 1H-NMR studies were extensively carried out to elucidate preferred conformations of a series of 14-membered cyclic dermorphin analogues containing two phenylalanines at both the third and fourth positions, e.g., Tyr-c[D-A2bu-Phe-Phe-(L and D)-Leu], Tyr-c[D-A2bu-Phe-gPhe-(S and R)-mLeu], and Tyr-c[D-Glu-Phe-gPhe-(L and D)-rLeu]. The temperature coefficients of the amide proton chemical shifts, vicinal 1H-1H coupling constants for the NH-CH groupings, and nuclear Overhauser effects provided information regarding the preferred conformations of the backbones. The conformational preferences and flexibility of the side chains were also estimated from the vicinal 1H-1H coupling constants around the C-C beta and C beta-C bonds in the articulated side chains. A comparison of the results obtained was made with the results previously obtained for the corresponding enkephalin analogues containing a glycine at the third position. It was found that the replacement of the glycine with the phenylalanine at the third position increases the conformational flexibility of the molecules with an L-, or S-, residue at the fifth position but reduces the flexibility of the molecules with D-, or R-, residue at the same position. The rotating frame nuclear Overhauser experiments gave direct evidence for compact conformations, with the Tyr side chain folding back over the 14-membered ring in Tyr-c[D-Glu-Phe-gPhe-rLeu], which displays relatively high selectivity for the delta-receptor over the mu-receptor. This observation is in agreement with our model proposed for the cyclic enkephalin analogues: folded forms with close aromatic ring placement are required for the activity at the delta-receptor.

Cyclic retro-inverso dipeptides with two aromatic side chains. II. Conformational analysis

The conformations of cis and trans cyclic retro-inverso dipeptides--2-[(4-hydroxy)benzyl]-5-benzyl-4,6(1H,2H,3H,5H)-pyrimidinedi one (c[mTyr-gPhe]), and 2-benzyl-5-amino-5-[(4-hydroxy)benzyl]-4,6(1H,2H,3H,5H)-pyrimidinedione (c[mTyr-gPhe]), and 2-benzyl-5-amino-5-[(4-hydroxy)benzyl]-4,6(1H,2H,3H,5H)-pyrimidinedione (c[(alpha-amino)mTyr-gPhe])--and the parent cyclic dipeptides--c[tyrosyl-phenylalanine] (cis-c[L-Tyr-L-Phe]) and c[tyrosyl-D-phenylalanine] (trans-c[L-Tyr-D-Phe])--were studied by using 1H-nmr spectroscopy and semiempirical energy calculations. In the cis compounds of all the cyclic retro-inverso and parent dipeptides, the most stable conformer has both aromatic side chains sharing the space over the backbone ring in a "face-to-face" fashion. All the trans compounds predominantly assume a "sandwich" conformation in which the two aromatic rings are folded back over the backbone ring on opposite sides. However, different conformational preferences were observed for the backbones between the retro-inverso and parent cyclic dipeptides. The parent cyclic dipeptide trans-c[L-Tyr-D-Phe] adopts two types of boat structures with different side-chain orientations in almost equal amounts: one with the Tyr side chain in a pseudoaxial position and the Phe side chain in a pseudoequatorial position, the other with the Tyr side chain in a pseudoequatorial position and the Phe side chain in a pseudoaxial position. On the other hand, the cyclic retro-inverso dipeptides trans-c[mPhe-gTyr] and trans c[mTyr-gPhe] assume only one type of boat structure in which the malonyl side chain is in a pseudoequatorial and the gem-diamino side chain is in a pseudoaxial position. In addition to the preferred conformations, the conformational energies of the C alpha--C beta bonds in the malonyl and gem-diamino residues were estimated from the temperature variation of vicinal 1H--1H coupling constants for the H--C alpha--C beta--H groupings observed for the trans isomers of cyclic retro-inverso dipeptides. The energies were evaluated to be 1.1 and 1.8 kcal mol-1 for the malonyl and gem-diamino residues, respectively. Applying these energies to the parent cyclic dipeptide trans-c[L-Tyr-D-Phe], the observed fractions of three side-chain conformations are reasonably reproduced. The conformational energies as well as conformational properties of the molecules estimated in this investigation may be useful to refine force constants for both parent and retro-inverso peptides with aromatic side chains.

Molecular dynamics study of the structure and dynamics of a protein molecule in a crystalline ionic environment, Streptomyces griseus protease A

A large-scale molecular dynamics simulation of the behavior of a serine protease (Streptomyces griseus protease A) in a crystalline environment has been performed. All atoms (including hydrogens) of two protein molecules and the surrounding solvent of crystallization, consisting of both water and salt ions, were explicitly represented, and a relatively long range of interactions (up to 15 A) were included. The simulation is the longest so far reported for a protein in such an environment (60 ps). The use of the full crystalline environment allows a direct comparison of the structure and dynamic properties of the protein and surrounding solvent to be made with the experimental X-ray structure. Here we report the comparison of the protein structures and analyze the energetics of the system, including interaction with the aqueous environment. Subsequent papers will deal with other aspects of the simulation. The overall root mean square differences between the time-averaged molecular dynamics structure and that from crystallography, for all well-ordered, non-hydrogen atoms, are 1.67 and 1.25 A for the two molecules taken as the asymmetric unit. An extensive analysis of the conformation of substructural elements and individual residues and their deviation from experiment has revealed a strong influence of the ionic medium on their behavior. Implications of the results for free energy calculations and for future directions are also discussed.

Accessible conformations of melanin-concentrating hormone: a molecular dynamics approach

Molecular dynamics simulations have been used to search for the accessible conformations of the melanin-concentrating hormone (MCH). The studies have been performed on native MCH and two of its peptide fragments, a cyclic MCH(5-14) fragment and a linear MCH(5-14) fragment. An analysis of the molecular dynamics trajectories of the three peptides indicates that two regions of the peptide have characteristic conformational properties that may be important for the biological activity. One is a region around Gly8, which is conformationally mobile, and the other is around Pro13, which shows unusual rigidity. The molecular dynamics simulation results are discussed in terms of backbone structural features like beta turns, side-chain interactions, and orientations of the disulfide bridge. The results of this analysis are used to suggest new analogues that will modify the conformational features of the peptide and further define the conformational requirements for activity. Finally, the results are related to nmr studies of the peptide and reveal agreements between the experimental nuclear Overhauser effect constraints and some of the accessible conformations obtained from the simulation.

Conformational perturbations in retro-analogs of the tBuCO-Ala-Gly-NHiPr dipeptide. Crystal structure of the retro-dipeptide with a reversed Ala-Gly amide bond

The three retro-analogs of the tBuCO-Ala-Gly-NHiPr dipeptide, in which each amide bond had been successively reversed, were studied in solution by 1H-n.m.r. and i.r. spectroscopy with reference to the conformational properties of their parent dipeptide. Reversal of the Ala-Gly amide bond proved to perturb the folding tendency of the backbone less than the inversion of either of the terminal amide bonds. The crystal structure of the retro-peptide containing a reversed Ala-Gly amide bond was also solved by X-ray diffraction and constitutes the first available data for this retro-peptide series. In contrast to the beta II-folded structure of the parent dipeptide, the retro-peptide molecule adopts an open conformation in the crystal.

Pseudopeptides and beta folding: x-ray structures compared with structures in solution

In order to restrain the flexibility of the peptide molecules and reduce their biodegradation, modifications of the main chain are now introduced in pseudopeptide analogues. Surprisingly, there is very little data on the conformational properties of these derivatives. We have examined pseudopeptide analogues of RCO-X-Y-NHR' model dipeptides in the depsi, N-methylated, reduced, retro, alpha, beta-dehydro, beta-amino acid, and hydrazino series, in the solid state by x-ray diffraction, and in solution by ir and 1H-nmr spectroscopy. This study provides us with accurate dimensions of the peptide surrogates, and gives some information on the conformational tendencies induced by these substitutions, with reference to those of the related dipeptide sequences.

Partial retro-inverso analogues of somatostatin: pairwise modifications at residues 7 and 8 and at residues 8 and 9

Peptide bonds between residues 7 and 8 residues 8 and 9, postulated internal cleavage sites of the peptide hormone somatostatin, were subjected to "pairwise" retro-inverso modification, where atoms of these peptide bonds were interchanged to give the analogues [gPhe7-m-(RS)-Trp8]somatostatin (I) and [gTrp8-m-(RS)-Lys9]somatostatin (II). Key fragments containing the modifications were synthesized by using [bis(trifluoroacetoxy)iodo]benzene for the generation of gem-diaminoalkyl-containing precursors from peptide amides. The versatility of solution synthetic methods was utilized to allow the incorporation of the modified segments. Protecting groups, removable selectively and under mild conditions, included tert-butyl-based groups for the side chains and the tert-butylmercapto group for the cysteine thiols. The excellent results obtained in the syntheses of analogues I and II, and previously of somatostatin on a larger scale [Moroder, L., Gemeiner, M., Goehring, W., Faeger, E., Thamm, P., & Wunsch, E. (1981) Biopolymers 20, 17-31], suggest the general feasibility of this route for the synthesis of centrally modified analogues. The purification of the products by Sephadex LH-20 chromatography afforded the separation of diastereomers of both analogues. The two isomers of I showed significant but different activities while those of analogue II were marginally active.