Crystal structure of cyclo(Gly-L-Pro-D-Phe-Gly-Val): an example of a new type of three-residue turn

Elucidating Solution Structures of Cyclic Peptides Using Molecular Dynamics Simulations

Protein-protein interactions are vital to biological processes, but the shape and size of their interfaces make them hard to target using small molecules. Cyclic peptides have shown promise as protein-protein interaction modulators, as they can bind protein surfaces with high affinity and specificity. Dozens of cyclic peptides are already FDA approved, and many more are in various stages of development as immunosuppressants, antibiotics, antivirals, or anticancer drugs. However, most cyclic peptide drugs so far have been natural products or derivatives thereof, with de novo design having proven challenging. A key obstacle is structural characterization: cyclic peptides frequently adopt multiple conformations in solution, which are difficult to resolve using techniques like NMR spectroscopy. The lack of solution structural information prevents a thorough understanding of cyclic peptides' sequence-structure-function relationship. Here we review recent development and application of molecular dynamics simulations with enhanced sampling to studying the solution structures of cyclic peptides. We describe novel computational methods capable of sampling cyclic peptides' conformational space and provide examples of computational studies that relate peptides' sequence and structure to biological activity. We demonstrate that molecular dynamics simulations have grown from an explanatory technique to a full-fledged tool for systematic studies at the forefront of cyclic peptide therapeutic design.

Understanding and designing head-to-tail cyclic peptides

Cyclic peptides (CPs) are an exciting class of molecules with a variety of applications. However, design strategies for CP therapeutics, for example, are generally limited by a poor understanding of their sequence-structure relationships. This knowledge gap often leads to a trial-and-error approach for designing CPs for a specific purpose, which is both costly and time-consuming. Herein, we describe the current experimental and computational efforts in understanding and designing head-to-tail CPs along with their respective challenges. In addition, we provide several future directions in the field of computational CP design to improve its accuracy, efficiency and applicability. These advances, combined with experimental techniques, shall ultimately provide a better understanding of these interesting molecules and a reliable working platform to rationally design CPs with desired characteristics.

Small head-to-tail macrocyclic α-peptoids

A convenient and efficient methodology for the head-to-tail macrocyclization of small 3-mer, 4-mer, and 5-mer α-peptoid acids (9-, 12-, and 15-atom N-substituted glycine oligomers) is described. The cyclic trimer has a ccc amide sequence in the crystal structure, whereas the tetramer has ctct and the pentamer has ttccc stereochemistry. NMR analysis reveals rigid structures in solution. These synthetic macrocycles may prove useful in medicinal and materials applications.

An extensible and systematic force field, ESFF, for molecular modeling of organic, inorganic, and organometallic systems

ESFF is a rule-based force field designed for modeling organic, inorganic, and organometallic systems. To cover this broad range of molecular systems, ESFF was developed in an extensible and systematic manner. Several unique features were introduced including pseudoangle and a dot product function representing torsion energy terms. The partial atomic charges that are topology-dependent are determined from ab initio (DFT) calculated electronegativity and hardness for valence orbitals. The van der Waals parameters are charge-dependent, and correlated with the ionization potential for atoms in various valence states. To obtain a set of well-defined and physically meaningful parameters, ESFF employs semiempirical rules to translate atomic-based parameters to parameters typically associated with a covalent valence force field. The atomic parameters depend not only on atom type, but also on internal type, thus resulting in a more accurate force field. This article presents the theory and the method used to develop the force field. The force field has been applied to molecular simulations of a wide variety of systems including nucleic acids, peptides, hydrocarbons, porphyrins, transition metal complexes, zeolites, and organometallic compounds. Agreement with the experimental results indicates that ESFF is a valuable tool in molecular simulations for understanding and predicting both crystal and gas phase molecular structures.

Structure of cyclic peptides: the crystal and solution conformation of cyclo(Phe-Phe-Aib-Leu-Pro)

A solid-state and solution conformation analyses of the cyclopentapeptide cyclo(Phe-Phe-Aib-Leu-Pro) has been carried out by X-ray diffraction and nuclear magnetic resonance techniques. The structure of the hexagonal crystals, grown from a methanol solution [a = b = 16.530(4) A, c = 21.356(9) A, space group P6(5), Z = 6], shows the presence of one intramolecular N-H ..O=C hydrogen bond with the formation of a gamma-turn (C7). The Aib3 residue, at the center of the gamma-turn, presents unexpected values of the torsion angles [phi = 70.5 degrees and psi = -73.8 degrees], which have been observed only once before for this helicogenic residue. A cis peptide bond occurs between Leu4 and Pro5; all other peptide bonds are trans. The overall conformation for the cyclopentapeptide with one cis-peptide bond on one side and an intramolecular gamma-turn on the opposite side results in an equatorial topology of the side-chains of the Phe1, Phe2 and Leu4 residues. Indeed, the Calpha-Cbeta and Cbeta-Cgamma bonds of these residues lie approximately in the mean plane of the cyclic ring system. The structure is compared with data in the literature on cyclic pentapeptides. In addition the Pro-Phe-Phe moiety shows a conformation similar to that observed in other larger cyclic bioactive peptides, which indicates a reduced number of conformations for this sequence. The solution study was carried out in three different solvent systems: chloroform, acetonitrile and methanol in the temperature interval 220-300 K. In all three solvents the room temperature spectra show that the peptide is conformationally nonhomogeneous. In acetonitrile at low temperatures it is possible to reduce the conformational equilibrium to two predominant conformers which differ for the cis-trans isomerism of the Leu4-Pro5 peptide bond.

Beta-alanine containing cyclic peptides with predetermined turned structure. V

In the present paper we describe the synthesis, purification, single crystal x-ray analysis, and solution structural characterization by nmr spectroscopy, combined with restrained molecular dynamic simulations, of the cyclic hexapeptide cyclo-(Pro-Phe-beta-Ala-Phe-Phe-beta-Ala). The peptide was synthesized by classical solution methods and the cyclization of the free hexapeptide was accomplished in good yields in diluted methylenechloride solution using N,N-dicyclohexyl-carbodiimide. The compound crystallizes in the monoclinic space group P2(1) from methanol/ethyl acetate. The molecule adopts in the solid state a conformation characterized by cis beta-Ala6-Pro1 peptide bond. The alpha-amino acid residues are at the corner positions of turned structures. The Pro1-Phe2 segment is incorporated in a pseudo type I beta-turn, while Phe4-Phe5 is in a typical type I beta-turn. Assignment of all 1H and 13C resonances was achieved by homo- and heteronuclear two-dimensional techniques in dimethylsulfoxide (DMSO) solutions. The conformational analysis was based on interproton distances derived from rotating frame nuclear Overhauser effect spectroscopy spectra and homonuclear coupling constants. Restrained molecular dynamic simulation in vacuo was also performed to built refined molecular models. The molecule is present in DMSO solution as two slowly interconverting conformers, characterized by a cis-trans isomerism around the beta-Ala6-Pro1 peptide bond. This work confirms our expectations on the low propensity of beta-alanyl residues to be positioned at the corners of turned structure.

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

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

Impact of a micellar environment on the conformations of two cyclic pentapeptides

In an effort to explore the influence of interfacial environments on reverse turns, we have performed a detailed analysis by nmr of the solution conformations of two cyclic pentapeptides in sodium dodecyl sulfate (SDS) micelles. The first peptide, cyclo (D-Phe1-Pro2-Gly3-D-Ala4-Pro5), adopts a single rigid conformation in solution (either chloroform or dimethylsulfoxide) and in crystals, whereas the second, cyclo (Gly1-Pro2-D-Phe3-Gly4-Val5), is much more flexible and adopts different conformations in the crystal and in solution. Both of these peptides are solubilized by SDS micelles, and nmr relaxation rates indicate that they are both partially immobilized by interaction with the micelles. Furthermore, some amide protons in both peptides participate in hydrogen bonds with water. In the presence of micelles, the former peptide retains a conformation essentially the same as that found in crystals and in solution, which consists of a beta turn and an inverse gamma turn. However, the micellar environment has a significant effect on the latter peptide. In particular, the population of a conformer containing a cis Gly-Pro peptide bond is increased significantly. The most likely conformation of the cis isomer, determined by a combination of nmr and restrained molecular dynamics, contains a Gly1-Pro2 delta turn and a gamma turn about D-Phe3. The nmr data on the trans isomer indicate that this isomer is averaging between two conformations that differ mainly in the orientation of the D-Phe3-Gly4 peptide bond.

Solution conformations of two flexible cyclic pentapeptides: cyclo(Gly-Pro-D-Phe-Gly-Ala) and cyclo(Gly-Pro-D-Phe-Gly-Val)

In an effort to explore the residue preferences in three-residue reverse turns (so-called gamma-turns), two cyclic pentapeptides--cyclo(Gly1-Pro2-D-Phe3-Gly4-Ala5) (I) and cyclo(Gly1-Pro2-D-Phe3-Gly4-Val5) (II)--have been synthesized and analyzed by nmr. It was anticipated that the Gly-Pro-D-Phe-Gly portions of these molecules would favor a beta-turn conformation, leaving the remainder of the molecule to adopt a gamma turn, as seen in several previously studied model cyclic pentapeptides. The nmr data for both peptides in CDCl3 (5% DMSO-d6) and in neat DMSO-d6 indicate that the most populated conformation contains a distorted beta turn around Pro2-D-Phe3, which includes a gamma turn around D-Phe3. The distortion in the beta turn does not impede the formation of an inverse gamma turn around residue 5, and indeed, this conformation is observed in both peptides. Both the alanine and the bulkier valine residues are therefore found to be compatible with an inverse gamma turn. Molecular dynamics simulations on the title peptides are reported in the following paper. These simulations indicate that there is conformational flexibility around the D-Phe3-Gly4 peptide bond, which enables the formation of the gamma turn around D-Phe3. The third paper in this series explores the impact of a micellar environment on conformational equilibria in II.

Cyclic pentapeptides as models for reverse turns: determination of the equilibrium distribution between type I and type II conformations of Pro-Asn and Pro-Ala beta-turns

Cyclic pentapeptides are excellent models for reverse turns and have been used extensively in our laboratory to explore the influence of different amino acid sequences on turn preference. This paper is divided into two parts: In the first, we review our previous studies of cyclic pentapeptides. We summarize work that demonstrates the range of conformations possible within the cyclic pentapeptide backbone, the importance of sequence chirality in determining the backbone fold, and the utility of these cyclic pentapeptides as models for various turns. In the second, we present new results on two cyclic pentapeptides that contain beta-turns with Pro-Ala or Pro-Asn sequences in the i + 1 and i + 2 positions. By stereochemical criteria, a type I beta-turn is expected to be preferred by such L-L sequences. On the other hand, in proteins Asn occurs frequently in the i + 2 position of type II turns. We asked whether the same propensity would be manifest in an isolated model peptide, and if so, what the interactions were that influenced the relative stability of the type I and type II turns. To address these questions we have compared the conformational behavior of two peptides: cyclo(Gly-Pro-Ala-D-Phe-Pro) and cyclo(D-Ala-Pro-Asn-Gly-Pro). From previous studies, we anticipated that both peptides would contain an inverse gamma-turn and a beta-turn which consisted of either Gly-Pro-Ala-D-Phe or D-Ala-Pro-Asn-Gly in positions i to i + 3, respectively. Nuclear magnetic resonance analysis confirms this overall backbone conformation. Furthermore, quantitative nuclear Overhauser effect measurements in combination with molecular dynamics simulations and torsionally-forced energy minimizations have enabled us to determine that both type I and type II beta-turns are present in equilibrium in these peptides. The introduction of Asn in position i + 2 shifts this equilibrium significantly towards type II. We have done preliminary assessment of the possible side-chain/backbone conformations that contribute to the shift in populations.