Backbone conformational preferences and pseudorotational ring puckering of 1-aminocyclopentane-1-carboxylic acid

Origin of Helical Screw Sense Selectivity Induced by Chiral Constrained Cα-Tetrasubstituted α-Amino Acids in Aib-based Peptides

The mechanisms behind the propensity of chiral constrained Cα-tetrasubstituted amino acids (cCTAAs) to induce one particular helical screw sense, when included in the Ac-Aib2-cCTAA-Aib2-NHMe peptide model, were studied through replica exchange molecular dynamics, potential of mean force, and quantum theory of atoms in molecules calculations. We observed that cCTAAs exert their effect on helical screw sense selectivity through the positioning of the side chain to generate steric hindrance in either the (-x, +y, +z) or (+x, +y, -z) sectors of a right-handed 3D Cartesian space, where the z axis corresponds to the axis of the helix and the Cα lies on the +y semiaxis (0, +y, 0). The different strengthening of the noncovalent interactions, also comprising C-H···O interactions, exerted by the cCTAA in the right-handed or left-handed helix was also found important to define the preference of a cCTAA for a particular helix screw sense.

Mechanism of stabilization of helix secondary structure by constrained Cα-tetrasubstituted α-amino acids

The theoretical basis behind the ability of constrained Cα-tetrasubstituted amino acids (CTAAs) to induce stable helical conformations has been studied through Replica Exchange Molecular Dynamics Potential of Mean Force Quantum Theory of Atoms In Molecules calculations on Ac-l-Ala-CTAA-l-Ala-Aib-l-Ala-NHMe peptide models. We found that the origin of helix stabilization by CTAAs can be ascribed to at least two complementary mechanisms limiting the backbone conformational freedom: steric hindrance predominantly in the (+x,+y,-z) sector of a right-handed 3D Cartesian space, where the z axis coincides with the helical axis and the Cα of the CTAA lies on the +y axis (0,+y,0), and the establishment of additional and relatively strong C-H···O interactions involving the CTAA.

Solvent-induced conformational flexibility of a bicyclic proline analogue: Octahydroindole-2-carboxylic acid

The conformational preferences of the N-acetyl-N'-methylamide derivatives of the four octahydroindole-2-carboxylic acid (Oic) stereoisomers have been investigated in the gas-phase and in aqueous solution using quantum mechanical calculations. In addition to the conformational effects provoked by the stereochemical diversity of Oic, which presents three chiral centers, results provide evidence of interesting and rather unusual features. The conformational preferences of the Oic stereoisomers in solution are only well described by applying a complete and systematic search process, results achieved by simple re-optimization of the gas-phase minima being very imprecise. This is because the conformational rigidity detected in the gas-phase, which is imposed by the chemical restrictions of the fused bicyclic skeleton, disappears in aqueous solution, the four stereoisomers behaving as flexible molecules in this environment. Thus, in general, the γ-turn is the only minimum energy conformation in the gas-phase while in aqueous solution the helical, polyproline-II and γ-turn motifs are energetically favored. Molecular dynamics simulations indicate that the conformational flexibility predicted by quantum mechanical calculations for the four Oic stereoisomers in solution is satisfactorily reproduced by classical force-fields.

Exploring the energy landscape of a molecular engineered analog of a tumor-homing peptide

Recently a new non-coded amino acid was designed as a replacement for Arg, to protect the tumor-homing pentapeptide CREKA (Cys-Arg-Glu-Lys-Ala) from proteases. This constrained Arg analog, denoted c(5)Arg, was engineered to also promote the stability of the CREKA bioactive conformation. The conformational profile of the CREKA analog obtained by replacing Arg by c(5)Arg has been extensively investigated in this work. Two molecular dynamics simulations-based strategies have been employed: a modified simulated annealing and replica exchange. Results obtained using both techniques show that the conformational features of the new analog fulfill the purpose of its design. The new CREKA analog not only preserves the main structural attributes found for the bioactive conformation of the parent peptide but also shows lower flexibility. Moreover, the conformational profile of the mutated peptide narrows towards the most stable structures previously observed for the parent CREKA peptide.

Conformational profile of a proline-arginine hybrid

The intrinsic conformational preferences of a new nonproteinogenic amino acid have been explored by computational methods. This tailored molecule, named ((β)Pro)Arg, is conceived as a replacement for arginine in bioactive peptides when the stabilization of folded turn-like conformations is required. The new residue features a proline skeleton that bears the guanidilated side chain of arginine at the C(β) position of the five-membered pyrrolidine ring, in either a cis or a trans orientation with respect to the carboxylic acid. The conformational profiles of the N-acetyl-N'-methylamide derivatives of the cis and trans isomers of ((β)Pro)Arg have been examined in the gas phase and in solution by B3LYP/6-31+G(d,p) calculations and molecular dynamics simulations. The main conformational features of both isomers represent a balance between geometric restrictions imposed by the five-membered pyrrolidine ring and the ability of the guanidilated side chain to interact with the backbone through hydrogen bonds. Thus, both cis- and trans-((β)Pro)Arg exhibit a preference for the α(L) conformation as a consequence of the interactions established between the guanidinium moiety and the main-chain amide groups.

NCAD, a database integrating the intrinsic conformational preferences of non-coded amino acids

Peptides and proteins find an ever-increasing number of applications in the biomedical and materials engineering fields. The use of non-proteinogenic amino acids endowed with diverse physicochemical and structural features opens the possibility to design proteins and peptides with novel properties and functions. Moreover, non-proteinogenic residues are particularly useful to control the three-dimensional arrangement of peptidic chains, which is a crucial issue for most applications. However, information regarding such amino acids--also called non-coded, non-canonical, or non-standard--is usually scattered among publications specialized in quite diverse fields as well as in patents. Making all these data useful to the scientific community requires new tools and a framework for their assembly and coherent organization. We have successfully compiled, organized, and built a database (NCAD, Non-Coded Amino acids Database) containing information about the intrinsic conformational preferences of non-proteinogenic residues determined by quantum mechanical calculations, as well as bibliographic information about their synthesis, physical and spectroscopic characterization, conformational propensities established experimentally, and applications. The architecture of the database is presented in this work together with the first family of non-coded residues included, namely, alpha-tetrasubstituted alpha-amino acids. Furthermore, the NCAD usefulness is demonstrated through a test-case application example.

Conformational preferences of 1-amino-2-phenylcyclohexanecarboxylic acid, a phenylalanine cyclohexane analogue

The intrinsic conformational preferences of the restricted phenylalanine analogue generated by including the alpha and beta carbon atoms into a cyclohexane ring (1-amino-2-phenylcyclohexanecarboxylic acid, c(6)Phe) have been determined using quantum mechanical calculations. Specifically, the conformational profile of the N-acetyl-N'-methylamide derivative of the c(6)Phe stereoisomers exhibiting either a cis or a trans relative orientation between the amino and phenyl substituents has been analyzed in different environments (gas phase, chloroform, and aqueous solutions). Calculations were performed using B3LYP, MP2, and HF methods combined with the 6-31+G(d,p) and 6-311++G(d,p) basis sets, and a self-consistent reaction-field (SCRF) method was applied to analyze the influence of the solvent. The amino acids investigated can be viewed as constrained phenylalanine analogues with a rigidly oriented aromatic side chain that may interact with the peptide backbone not only sterically but also electronically through the aromatic pi orbitals. Their conformational propensities have been found to be strongly influenced by the specific orientation of the aromatic substituent in each stereoisomer and the conformation adopted by the cyclohexane ring, as well as by the environment.

In silico molecular engineering for a targeted replacement in a tumor-homing peptide

A new amino acid has been designed as a replacement for arginine (Arg, R) to protect the tumor-homing pentapeptide CREKA (Cys-Arg-Glu-Lys-Ala) from proteases. This amino acid, denoted (Pro)hArg, is characterized by a proline skeleton bearing a specifically oriented guanidinium side chain. This residue combines the ability of Pro to induce turn-like conformations with the Arg side-chain functionality. The conformational profile of the CREKA analogue incorporating this Arg substitute has been investigated by a combination of simulated annealing and molecular dynamics. Comparison of the results with those previously obtained for the natural CREKA shows that (Pro)hArg significantly reduces the conformational flexibility of the peptide. Although some changes are observed in the backbone...backbone and side-chain...side-chain interactions, the modified peptide exhibits a strong tendency to accommodate turn conformations centered at the (Pro)hArg residue and the overall shape of the molecule in the lowest energy conformations characterized for the natural and the modified peptides exhibit a high degree of similarity. In particular, the turn orients the backbone such that the Arg, Glu, and Lys side chains face the same side of the molecule, which is considered important for bioactivity. These results suggest that replacement of Arg by (Pro)hArg in CREKA may be useful in providing resistance against proteolytic enzymes while retaining conformational features which are essential for tumor-homing activity.

Ab initio investigation of the hydration of deprotonated amino acids

The complexation of five deprotonated anionic amino acids (glycine, L-alanine, L-valine, L-Aspartic acid, and L-glutamine) with one water molecule, has been investigated using a MP2/63-11++G(d,p) approach fully accounting for the basis set superposition errors. For each amino acid, several energetic minima have been identified, and we provide spectroscopic information allowing to discriminate them. Our results strongly suggest that two complexes should coexist under the experimental conditions for [Ala - H](-), [Val - H](-), and [Asp - H](-). Comparisons with the experimental enthalpies, entropies, and Gibbs free energies recently obtained by Wincel [J. Am. Soc. Mass Spectrom. 2008, 19, 1091-1097] show that our simulation reproduces the most significant structure/energy experimental trends, though the entropic changes induced by hydration are slightly overestimated.

Side-chain to backbone interactions dictate the conformational preferences of a cyclopentane arginine analogue

The intrinsic conformational preferences of the nonproteinogenic amino acids constructed by incorporating the arginine side chain in the beta position of 1-aminocyclopentane-1-carboxylic acid (either in a cis or a trans orientation relative to the amino group) have been investigated by using computational methods. These compounds may be considered as constrained analogues of arginine (denoted as c(5)Arg) in which the orientation of the side chain is fixed by the cyclopentane moiety. Specifically, the N-acetyl-N'-methylamide derivatives of cis- and trans-c(5)Arg have been examined in the gas phase and in solution by using B3LYP/6-311+G(d,p) calculations and Molecular Dynamics simulations. Results indicate that the conformational space available to these compounds is highly restricted, their conformational preferences being dictated by the ability of the guanidinium group in the side chain to establish hydrogen bond interactions with the backbone. A comparison with the behavior previously described for the analogous phenylalanine derivatives is presented.

Recent Progress on the Stereoselective Synthesis of Cyclic Quaternary alpha-Amino Acids

The most recent papers describing the stereoselective synthesis of cyclic quaternary alpha-amino acids are collected in this review. The diverse synthetic approaches are classified according to the size of the ring and taking into account the bond that is formed to complete the quaternary skeleton.

Structural analysis of a beta-helical protein motif stabilized by targeted replacements with conformationally constrained amino acids

Here we study conformational stabilization induced in a beta-helical nanostructure by position-specific mutations. The nanostructure is constructed through the self-assembly of the beta-helical building block excised from E. coli galactoside acetyltransferase (PDB code 1krr , chain A; residues 131-165). The mutations involve substitutions by cyclic, conformationally constrained amino acids. Specifically, a complete structural analysis of the Pro-Xaa-Val sequence [with Xaa being Gly, Ac 3c (1-aminocyclopropane-1-carboxylic acid) and Ac 5c (1-aminocyclopentane-1-carboxylic acid)], corresponding to the 148-150 loop region in the wild-type (Gly) and mutated (Ac 3c and Ac 5c) 1krr , has been performed using Molecular Dynamics simulations and X-ray crystallography. Simulations have been performed for the wild-type and mutants of three different systems, namely the building block, the nanoconstruct and the isolated Pro-Xaa-Val tripeptide. Furthermore, the crystalline structures of five peptides of Pro-Xaa-Val or Xaa-Val sequences have been solved by X-ray diffraction analysis and compared with theoretical predictions. Both the theoretical and crystallographic studies indicate that the Pro-Ac n c-Val sequences exhibit a high propensity to adopt turn-like conformations, and this propensity is little affected by the chemical environment. Overall, the results indicate that replacement of Gly149 by Ac 3c or Ac 5c significantly reduce the conformational flexibility of the target site enhancing the structural specificity of the building block and the nanoconstruct derived from the 1krr beta-helical motif.

Intrinsic conformational preferences of C(alpha,alpha)-dibenzylglycine

The intrinsic conformational preferences of C (alpha,alpha)-dibenzylglycine, a symmetric alpha,alpha-dialkylated amino acid bearing two benzyl substituents on the alpha-carbon atom, have been determined using quantum chemical calculations at the B3LYP/6-31+G(d,p) level. A total of 46 minimum energy conformations were found for the N-acetyl- N'-methylamide derivative, even though only nine of them showed a relative energy lower than 5.0 kcal/mol. The latter involves C 7, C 5, and alpha' backbone conformations stabilized by intramolecular hydrogen bonds and/or N-H...pi interactions. Calculation of the conformational free energies in different environments (gas-phase, carbon tetrachloride, chloroform, methanol, and water solutions) indicates that four different minima (two C 5 and two C 7) are energetically accessible at room temperature in the gas phase, while in methanol and aqueous solutions one such minimum (C 5) becomes the only significant conformation. Comparison with results recently reported for C (alpha,alpha)-diphenylglycine indicates that substitution of phenyl side groups by benzyl enhances the conformational flexibility leading to (i) a reduction of the strain of the peptide backbone and (ii) alleviating the repulsive interactions between the pi electron density of the phenyl groups and the lone pairs of the carbonyl oxygen atoms.

Application of 1-aminocyclohexane carboxylic acid to protein nanostructure computer design

Conformationally restricted amino acids are promising candidates to serve as basic pieces in redesigned protein motifs which constitute the basic modules in synthetic nanoconstructs. Here we study the ability of constrained cyclic amino acid 1-aminocyclohexane-1-carboxylic acid (Ac6c) to stabilize highly regular beta-helical motifs excised from naturally occurring proteins. Calculations indicate that the conformational flexibility observed in both the ring and the main chain is significantly higher than that detected for other 1-aminocycloalkane-1-carboxylic acids (Acnc, where n refers to the size of the ring) with smaller cycles. Incorporation of Ac6c into the flexible loops of beta-helical motifs indicates that the stability of such excised building blocks as well as the nanoassemblies derived from them is significantly enhanced. Thus, the intrinsic Ac6c tendency to adopt folded conformations combined with the low structural strain of the cyclohexane ring confers the ability to both self-adapt to the beta-helix motif and to stabilize the overall structure by absorbing part of its conformational fluctuations. Comparison with other Acnc residues indicates that the ability to adapt to the targeted position improves considerably with the ring size, i.e., when the rigidity introduced by the strain of the ring decreases.

1-amino-2-phenylcyclopentane-1-carboxylic acid: a conformationally restricted phenylalanine analogue

DFT calculations at the B3LYP/6-311G(d,p) level have been used to investigate the intrinsic conformational preferences of 1-amino-2-phenylcyclopentane-1-carboxylic acid (c5Phe), a constrained analogue of phenylalanine in which the alpha and beta carbons are included in a cyclopentane ring. Specifically, the N-acetyl-N'-methylamide derivatives of the cis and trans stereoisomers, where cis and trans refer to the relative position between the amino group and the phenyl ring, have been calculated. Solvent effects have been examined using a self-consistent reaction field (SCRF) method. Results indicate that the conformational space of the cis stereoisomer is much more restricted than that of the trans derivative both in the gas phase and in solution.

Stability of tubular structures based on beta-helical proteins: self-assembled versus polymerized nanoconstructs and wild-type versus mutated sequences

In this work we used atomistic molecular dynamics simulations to examine different aspects of tubular nanostructures constructed using protein building blocks with a beta-helical conformation. Initially, we considered two different natural protein building blocks, which were extracted from the protein data base, to compare the relative stabilities of the nanotubes obtained made of self-assembled and covalently linked repeats. Results show nanotubes constructed by linking building blocks through covalent bonds are very stable suggesting that the basic principles of polymer physics are valid when the repeating units are made of large fragments of proteins. In contrast, the stability of self-assembled nanostructures strongly depends on the attractive nonbonding interactions associated to building blocks aligned in a complementary manner. On the other hand, we investigated the ability of a conformationally constrained synthetic amino acid to enhance the stability of both self-assembled and polymerized nanotubes when it is used to substitute natural residues. Specifically, we considered 1-aminocyclopentane-1-caboxylic acid, which involves strong stereochemical constraints produced by the cyclopentane side chain. We found that the incorporation of this amino acid within the more flexible regions of the beta-helical building blocks is an excellent strategy to enhance the stability of the nanotubes. Thus, when a single mutation is performed in the loop region of the beta-helix, the bend architecture of the whole loop is stabilized since the conformational mobility is reduced not only at the mutated position but also at the adjacent positions.

Conformations of proline analogues having double bonds in the ring

The intrinsic conformational preferences of proline analogues having double bonds between carbon atoms in their rings have been investigated using quantum mechanical calculations at the B3LYP/6-31+G(d,p) level. For this purpose, the potential energy surface of the N-acety-N'-methylamide derivatives of three dehydroprolines (proline analogues unsaturated at alpha,beta; beta,gamma; and gamma,delta) and pyrrole (proline analogue with unsaturations at both alpha,beta and gamma,delta) have been explored, and the results are compared with those obtained for the derivative of the nonmodified proline. We found that the double bonds affect the ring puckering and the geometric internal parameters, even though the backbone conformation was influenced the most. Results indicate that the formation of double bonds between carbon atoms in the pyrrolidine ring should be considered as an effective procedure to restrict the conformational flexibility of prolines. Interestingly, we also found that the N-acetyl-N'-methylamide derivative of pyrrole shows a high probability of having a cis peptide bond preceding the proline analogue.

Intrinsic conformational characteristics of alpha,alpha-diphenylglycine

Quantum mechanical calculations at the B3LYP/6-31+G(d,p) level have been used to investigate the intrinsic conformational preferences of alpha,alpha-diphenylglycine, a simple alpha,alpha-dialkylated amino acid bearing two phenyl substituents on the alpha-carbon, in both the gas phase and aqueous solution. Nine minimum energy conformations have been characterized for the N-acetyl-N'-methylamide derivative within a relative energy range of about 9 kcal/mol. The relative stability of these structures is largely influenced by specific backbone...side chain and side chain...side chain interactions that can be attractive (N-H...pi and C-H...pi) or repulsive (C=O...pi). On the other hand, comparison with the minimum energy conformations calculated for alpha-aminoisobutyric acid, in which the two phenyl substituents are replaced by methyl groups, revealed that the bulky aromatic rings of alpha,alpha-diphenylglycine induce strain in the internal geometry of the peptide. Finally, a set of force-field parameters for classical Molecular Mechanics calculations was developed for the investigated amino acid. Molecular Dynamics simulations in aqueous solutions have been carried out to validate the parameters obtained.