Intrinsically stable secondary structure elements of proteins: a comprehensive study of folding units of proteins by computation and by analysis of data determined by X-ray crystallography

Different protein architectures show strong similarities regardless of their amino acid composition: the backbone folds of the different secondary structural elements exhibit nearly identical geometries. To investigate the principles of folding and stability properties, oligopeptide models (that is, HCO-(NH-L-CHR-CO)(n)-NH(2)) have been studied. Previously, ab initio structure determinations have provided a small amount of information on the conformational building units of di- and tripeptides. A maximum of nine differently folded backbone types is available for any natural alpha-amino acid residue, with the exception of proline. All of these conformers have different relative energies. The present study compiles an ab inito database of optimized HCO-(L-Xxx)(n)-NH(2) structures, where 1相似文献   

A theoretical case study of type I and type II beta-turns

NMR chemical shielding anisotropy tensors have been computed by employing a medium size basis set and the GIAO-DFT(B3LYP) formalism of electronic structure theory for all of the atoms of type I and type II beta-turn models. The models contain all possible combinations of the amino acid residues Gly, Ala, Val, and Ser, with all possible side-chain orientations where applicable in a dipeptide. The several hundred structures investigated contain either constrained or optimized phi, psi, and chi dihedral angles. A statistical analysis of the resulting large database was performed and multidimensional (2D and 3D) chemical-shift/chemical-shift plots were generated. The (1)H(alpha-13)C(alpha), (13)C(alpha-1)H(alpha-13)C(beta), and (13)C(alpha-1)H(alpha-13)C' 2D and 3D plots have the notable feature that the conformers clearly cluster in distinct regions. This allows straightforward identification of the backbone and side-chain conformations of the residues forming beta-turns. Chemical shift calculations on larger For-(L-Ala)(n)-NH(2) (n=4, 6, 8) models, containing a single type I or type II beta-turn, prove that the simple models employed are adequate. A limited number of chemical shift calculations performed at the highly correlated CCSD(T) level prove the adequacy of the computational method chosen. For all nuclei, statistically averaged theoretical and experimental shifts taken from the BioMagnetic Resonance Bank (BMRB) exhibit good correlation. These results confirm and extend our previous findings that chemical shift information from selected multiple-pulse NMR experiments could be employed directly to extract folding information for polypeptides and proteins.     

Mid-infrared spectroscopy of protected peptides in the gas phase: a probe of the backbone conformation

Infrared/UV hole-burning spectroscopy is performed on individual conformers of the protected dipeptide Z-Aib-Pro-NHMe. The extended IR range probed in this study allows one to elucidate both the H-bonding motif (5-7 microm) as well as the backbone structure (7-10 microm). Comparison with DFT calculations shows that the backbone is locally constrained to an alpha-conformation by Aib and to a gamma-turn by Pro. The gamma-turn motif observed here is intriguing since the condensed phase structure is known to be a beta-turn. This is the first actual observation of such a discrepancy, and it emphasizes the subtle balance between intra- and intermolecular forces, which is responsible for the relative stability of the different secondary structure motifs.     

Toward a rational design of beta-peptide structures

Intrinsic conformational characteristics of beta-peptides built up from simple achiral and chiral beta-amino acid residues (i.e., HCO-beta-Ala-NH2, HCO-beta-Abu-NH2) were studied using quantum chemical calculations and 1H-NMR spectroscopy. A conformer-based systematic and uniform nomenclature was introduced to differentiate conformers. Geometry optimizations were performed on all homoconformers of both HCO-(beta-Ala)(k)-NH2 and HCO-(beta-Abu)(k)-NH2 (1 < or = k < or = 6) model systems at the RHF/3-21G and RHF/6-311++G(d, p) levels of theory. To test for accuracy and precision, additional computations were carried out at several levels of theory [e.g., RHF/6-31G(d), and B3LYP/6-311++G(d, p)]. To display the folding preference, the relative stability of selected conformers as function of the length of the polypeptide chain was determined. Ab initio population distribution of hexapeptides and the conformational ensemble of synthetic models composed of beta-Ala and beta-Abu studied using 1H-NMR in different solvents were compared at a range of temperatures. Helical preference induced by various steric effects of nonpolar side chains was tested using higher level ab initio methods for well-known model systems such as: HCO-(beta-HVal-beta-HAla-beta-HLeu)2-NH2, HCO-(ACHC)6-NH2, HCO-(trans-ACPC)6-NH2, and HCO-(cis-ACPC)6-NH2. The relative stabilities determined by theoretical methods agreed well with most experimental data, supporting the theory that the local conformational preference influenced by steric effects is a key determining factor of the global fold both in solution and in the gas phase.     

A systematic study of minima in alanine dipeptide

The alanine dipeptide is a standard system to model dihedral angles in proteins. It is shown that obtaining the Ramachandran plot accurately is a hard problem because of many local minima; depending on the details of geometry optimizations, different Ramachandran plots can be obtained. To locate all energy minima, starting from geometries from MD simulations, 250,000 geometry optimizations were performed at the level of RHF/6-31G*, followed by re-optimizations of the located 827 minima at the level of MP2/6–311++G**, yielding 30 unique minima, most of which were not previously reported in literature. Both in vacuo and solvated structures are discussed. The minima are systematically categorized based on four backbone dihedral angles. The Gibbs energies are evaluated and the structural factors determining the relative stabilities of conformers are discussed. © 2018 Wiley Periodicals, Inc.     

Peptide models. XXXIII. Extrapolation of low-level Hartree-Fock data of peptide conformation to large basis set SCF,MP2, DFT,and CCSD(T) results. The Ramachandran surface of alanine dipeptide computed at various levels of theory

At the dawn of the new millenium, new concepts are required for a more profound understanding of protein structures. Together with NMR and X-ray-based 3D-structure determinations in silico methods are now widely accepted. Homology-based modeling studies, molecular dynamics methods, and quantum mechanical approaches are more commonly used. Despite the steady and exponential increase in computational power, high level ab initio methods will not be in common use for studying the structure and dynamics of large peptides and proteins in the near future. We are presenting here a novel approach, in which low- and medium-level ab initio energy results are scaled, thus extrapolating to a higher level of information. This scaling is of special significance, because we observed previously on molecular properties such as energy, chemical shielding data, etc., determined at a higher theoretical level, do correlate better with experimental data, than those originating from lower theoretical treatments. The Ramachandran surface of an alanine dipeptide now determined at six different levels of theory [RHF and B3LYP 3-21G, 6-31+G(d) and 6-311++G(d,p)] serves as a suitable test. Minima, first-order critical points and partially optimized structures, determined at different levels of theory (SCF, DFT), were completed with high level energy calculations such as MP2, MP4D, and CCSD(T). For the first time three different CCSD(T) sets of energies were determined for all stable B3LYP/6-311++G(d,p) minima of an alanine dipeptide. From the simplest ab initio data (e.g., RHF/3-21G) to more complex results [CCSD(T)/6-311+G(d,p)//B3LYP/6-311++G(d,p)] all data sets were compared, analyzed in a comprehensive manner, and evaluated by means of statistics.     

Secondary structures of short peptide chains in the gas phase: double resonance spectroscopy of protected dipeptides

The conformational structure of short peptide chains in the gas phase is studied by laser spectroscopy of a series of protected dipeptides, Ac-Xxx-Phe-NH(2), Xxx=Gly, Ala, and Val. The combination of laser desorption with supersonic expansion enables us to vaporize the peptide molecules and cool them internally; IR/UV double resonance spectroscopy in comparison to density functional theory calculations on Ac-Gly-Phe-NH(2) permits us to identify and characterize the conformers populated in the supersonic expansion. Two main conformations, corresponding to secondary structures of proteins, are found to compete in the present experiments. One is composed of a doubly gamma-fold corresponding to the 2(7) ribbon structure. Topologically, this motif is very close to a beta-strand backbone conformation. The second conformation observed is the beta-turn, responsible for the chain reversal in proteins. It is characterized by a relatively weak hydrogen bond linking remote NH and CO groups of the molecule and leading to a ten-membered ring. The present gas phase experiment illustrates the intrinsic folding properties of the peptide chain and the robustness of the beta-turn structure, even in the absence of a solvent. The beta-turn population is found to vary significantly with the residues within the sequence; the Ac-Val-Phe-NH(2) peptide, with its two bulky side chains, exhibits the largest beta-turn population. This suggests that the intrinsic stabilities of the 2(7) ribbon and the beta-turn are very similar and that weakly polar interactions occurring between side chains can be a decisive factor capable of controlling the secondary structure.     

Peptide models XXIII. Conformational model for polar side‐chain containing amino acid residues: A comprehensive analysis of RHF,DFT, and MP2 properties of HCO‐L‐SER‐NH2

Geometric and energetic properties of a diamide of serine, HCO‐NH‐L ‐CH(CH₂OH)CO‐NH₂, are investigated by standard methods of computational quantum chemistry. Similarly to other amino acid residues, conformational properties of HCO‐L ‐Ser‐NH₂ can be derived from the analysis of its E=E(ϕ,ψ;χ₁,χ₂) hypersurface. Reoptimization of 44 RHF/3‐21G conformers at the RHF/6‐311++G** level resulted in 36 minima. For all conformers, geometrical properties, including variation of H‐bond parameters and structural shifts in the torsional space, are thoroughly investigated. Results from further single‐point energy calculations at the RHF, DFT, and MP2 levels, performed on the entire conformational data set, form a database of 224 energy values, perhaps the largest set calculated so far for any single amino acid diamide. A comprehensive analysis of this database reveals significant correlation among energies obtained at six levels of ab initio theory. Regression parameters provide an opportunity for extrapolation in order to predict the energy of a conformer at a high level by doing explicit ab initio computations only for a few selected conformers. The computed conformational and relative energy data are compared with structural and occurrence results derived from a nonhomologous protein database incorporating 1135 proteins. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 626–655, 2000     

Peptide models XXXI. Conformational properties of hydrophobic residues shaping the core of proteins. An ab initio study of N‐formyl‐L‐valinamide and N‐formyl‐L‐phenylalaninamide

Employing introductory (3‐21G RHF) and medium‐size (6‐311++G** B3LYP) ab initio calculations, complete conformational libraries, containing as many as 27 conformers, have been determined for diamide model systems incorporating the amino acids valine (Val) and phenylalanine (Phe). Conformational and energetic properties of these libraries were analyzed. For example, significant correlation was found between relative energies from 6‐311++G** B3LYP and single‐point B3LYP/6‐311++G**//RHF/3‐21G calculations. Comparison of populations of molecular conformations of hydrophobic aromatic and nonaromatic residues, based on their ab initiorelative energies, with their natural abundance indicates that, at least for the hydrophobic core of proteins, the conformations of Val (Ile, Leu) and Phe (Tyr, Trp) are controlled by the local energetic preferences of the respective amino acids. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 732–751, 2001     

Probing the structural determinants of type II' beta-turn formation in peptides and proteins

The structural determinants of type II' beta-turns were probed through a comprehensive CD, NMR, and molecular dynamics analysis of 10 specially designed beta-hairpin peptides. The peptide model used in this study is a synthetic, water-soluble, 14-residue cyclic analogue of gramicidin S which contains two well-defined type II' beta-turns connected by a highly stable, amphipathic, antiparallel beta-sheet. A variety of coded and noncoded amino acids were systematically substituted in one of the two type II' turns to analyze the effects of backbone chirality, side-chain steric restriction, and side-chain/side-chain interactions. beta-Sheet content (as measured through a variety of experimental methods), molecular dynamics, and 3D structural analysis of the turn regions were used to assess the effects of each amino acid substitution on type II' beta-turn stabilization. Our results demonstrate that backbone heterochirality, which determines equatorial and axial side-chain orientation at the i+1 and i+2 residues of type II' turns, may account for up to 60% of type II' beta-turn stabilization. Steric restriction through side-chain N-alkylation appears to enhance type II' beta-turn propensity and may account for up to 20% of type II' beta-turn stabilization. Finally, aromatic/proline side-chain interactions appear to account for approximately 10% of type II' beta-turn stabilization. We believe this information could be particularly useful for the prediction of beta-turn propensity, the development of peptide-based drugs, and the de novo design of peptides, proteins, and peptidyl mimetics.     

Spiro beta-lactams as beta-turn mimetics. Design, synthesis, and NMR conformational analysis

Molecular modeling calculations using high-level ab initio methods (MP2/6-31+G) of a new type of spiro beta-lactams predict that these systems could adopt a beta-turn secondary structure in solution. Strong intramolecular hydrogen bonds stabilize the beta-turn conformation with a geometry that is very close to the ideal type II beta-turns. The synthesis of the spiro beta-lactams is achieved by Staudinger reaction of a cyclic ketene derived from N-bencyloxycarbonyl-L-proline acid chloride with an imine. This reaction allows the formation of the spiranic backbone in a single-step with high diastereoselectivity and good yields. The new spiro beta-lactams obtained are the core for the preparation of different types of peptidomimetics using well-established peptide chemistry. The NMR conformational analysis shows that these compounds adopt beta-turn conformation as predicted by the theoretical studies.     

Electronic structure and conformational properties of (carboxy-alkenyl)-phosphonic acids

The general conformational properties and electronic structure of (carboxy-alkenyl)-phosphonic derivatives were determined at RHF/STO-3G^* level. In all the series, low rotation barriers were found for the two C=C/P=O conformers. In the compounds in which the interactions between the carboxylic and phosphonic moieties are smaller, the most stable conformers are the C=C/P=O s-cis ones. In most of the conformers, the C=C/C=O system presents the disposition s-cis. The Z-(2-carboxy-vinyl) and Z-(2-carboxy-propenyl) phosphonic acids present intramolecular hydrogen bonds, existing in at least four conformer with internal hydrogen bonds. These last compounds were more rigorously studied at RHF/3-21G^* and RHF/6-31G^** levels. The most stable conformer shows a trans structure for the C=C/P=O angle, with an intramolecular hydrogen bond located between the hydroxylic hydrogen of phosphonic group and the carbonyl oxygen of carboxylic moiety. A secondary conformer is found with a double intramolecular hydrogen bond between two hydroxylic hydrogens of the phosphonic moiety and the oxygen of carboxylic bond. Another secondary conformer appears with an intramolecular hydrogen bond between the oxygen of the phosphoryl bond and the hydroxylic hydrogen of the carboxylic group. A study of the topology of charge densities is carried out. This analysis reveals bonds with an ionic participation. A very weak π conjugation, variable with the conformers, is found in the C=C/P=O system, as well as a strongly polarized P=O partial triple bond. The intramolecular hydrogen bonds give rise to cyclic structures.     

Extensive conformational searches of 13 representative dipeptides and an efficient method for dipeptide structure determinations based on amino acid conformers

Conformations of peptides are the basis for their property studies and the predictions of peptide structures are highly important in life science but very complex in practice. Here, thorough searches on the potential energy surfaces of 13 representative dipeptides by considering all possible combinations of the bond rotational degrees of freedom are performed using the density functional theory based methods. Careful analyses of the conformers of the 13 dipeptides and the corresponding amino acids reveal the connections between the structures of dipeptide and amino acids. A method for finding all important dipeptide conformers by optimizing a small number of trial structures generated by suitable superposition of the parent amino acid conformations is thus proposed. Applying the method to another eight dipeptides carefully examined by others shows that the new approach is both highly efficient and reliable by providing the most complete ensembles of dipeptide conformers and much improved agreements between the theoretical and experimental IR spectra. The method opens the door for the determination of the stable structures of all dipeptides with a manageable amount of effort. Preliminary result on the applicability of the method to the tripeptide structure determination is also presented. The results are the first step towards proving Anfinsen's hypothesis by revealing the relationships between the structures of the simplest peptide and its constituting amino acids. It implies that the structures of peptides are not only determined by their amino acid sequences, but also closely linked with the amino acid conformations. © 2009 Wiley Periodicals, Inc. J Comput Chem 2009     

Impact of azaproline on amide cis-trans isomerism: conformational analyses and NMR studies of model peptides including TRH analogues

The beta-turn is a well-studied motif in both proteins and peptides. Four residues, making almost a complete 180 degree-turn in the direction of the peptide chain, define the beta-turn. Several types of the beta-turn are defined according to Phi and Psi torsional angles of the backbone for residues i + 1 and i + 2. One special type of beta-turn, the type VI-turn, usually contains a proline with a cis-amide bond at residue i + 2. In an aza-amino acid, the alpha-carbon of the amino acid is changed to nitrogen. Peptides containing azaproline (azPro) have been shown to prefer the type VI beta-turn both in crystals and in organic solvents by NMR studies. MC/MD simulations using the GB/SA solvation model for water explored the conformational preferences of azPro-containing peptides in aqueous systems. An increase in the conformational preference for the cis-amide conformer of azPro was clearly seen, but the increased stability was relatively minor with respect to the trans-conformer as compared to previous suggestions. To test the validity of the calculations in view of the experimental data from crystal structures and NMR in organic solvents, [azPro(3)]-TRH and [Phe(2), azPro(3)]-TRH were synthesized, and their conformational preferences were determined by NMR in polar solvents as well as the impact of the azPro substitution on their biological activities.     

Local protein backbone folds determined by calculated NMR chemical shifts

NMR chemical shifts (CSs: δN(NH), δC(α), δC(β), δC', δH(NH), and δH(α)) were computed for the amino acid backbone conformers (α(L), β(L), γ(L), δ(L), ε(L), α(D), γ(D), δ(D), and ε(D) [Perczel et al., J Am Chem Soc 1991, 113, 6256]) modeled by oligoalanine structures. Topological differences of the extended fold were investigated on single β-strands, hairpins with type I and II β-turns, as well as double- and triple-stranded β-sheet models. The so-called "capping effect" was analyzed: residues at the termini of a homoconformer sequence unit usually have different CSs than the central residues of an adequately long homoconformer model. In heteroconformer sequences capping effect ruins the direct applicability of several chemical shift types (δH(NH), δC', and δN(NH)) for backbone structure determination of the parent residue. Experimental δH(α), δC(α), and δC(β) values retrieved from protein database are in good agreement with the relevant computed data in the case of the common backbone conformers (α(L), β(L), γ(L), and ε(L)), even though neighboring residue effects were not accounted for. Experimental and computed ΔδH(α)-ΔδC(α), ΔδH(α)-ΔδC(β), and ΔδC(α)-ΔδC(β) maps give qualitatively the same picture, that is, the positions of the backbone conformers relative to each other are very similar. This indicates that the H(α), C(α), and C(β) chemical shifts of alanine depend considerably on the backbone fold of the parent residue also in proteins. We provide tabulated CSs of the chiral amino acids that may predict the various structures of the residues.     

On the flexibility of beta-peptides

The full conformational space was explored for an achiral and two chiral beta-peptide models: namely For-beta-Ala-NH2, For-beta-Abu-NH2, and For-beta-Aib-NH2. Stability and conformational properties of all three model systems were computed at different levels of theory: RHF/3-21G, B3LYP/6-311++G(d,p)//RHF/3-21G, B3LYP/6-311++G(d,p), MP2//B3LYP/6-311++G(d,p), CCSD//B3LYP/6-311++G(d,p), and CCSD(T)//B3LYP/6-311++G(d,p). In addition, ab initio E = E(phi, micro, psi) potential energy hypersurfaces of all three models were determined, and their topologies were analyzed to determine the inherent flexibility properties of these beta-peptide models. Fewer points were found and assigned than expected on the basis of Multidimensional Conformational Analysis (MDCA). Furthermore, it has been demonstrated, that the four-dimensional surface, E = E(phi, mu, psi), can be reduced into a three-dimensional one: E = E[phi, f(phi), psi]. This reduction of dimensionality of freedom of motion suggests that beta-peptides are less flexible than one would have thought. This agrees with experimental data published on the conformational properties of peptides composed of beta-amino acid residues.     

Toward direct determination of conformations of protein building units from multidimensional NMR experiments I. A theoretical case study of For‐Gly‐NH2 and For‐L‐Ala‐NH2

NMR chemical shielding anisotropy tensors have been computed, employing several basis sets and the GIAO‐RHF and GIAO‐MP2 formalisms of electronic structure theory, for all the atoms of the five and nine typical backbone conformers of For‐Gly‐NH₂ and For‐L ‐Ala‐NH₂, respectively. Multidimensional chemical shift plots, as a function of the respective backbone fold, have been generated for both peptide models. On the 2D ¹H^NH‐¹⁵N^NH and ¹⁵N^NH‐¹³C^α plots the most notable feature is that at all levels of theory studied the backbone conformers cluster in different regions. Computed chemical shifts, as well as their averages, have been compared to relevant experimental values taken from the BioMagnetic Resonance Bank (BMRB). At the highest levels of theory, for all nuclei but the amide protons, deviations between statistically averaged theoretical and experimental shifts are as low as 1–3%. These results indicate that chemical shift information from selected multiple‐pulse NMR experiments (e.g., 2D‐HSQC and 3D‐HNCA) could directly be employed to extract folding information for polypeptides and proteins. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 882–900, 2000     

Delta-peptides and delta-amino acids as tools for peptide structure design--a theoretical study

An overview on all possible helix types in oligomers of delta-amino acids (delta-peptides) and their stabilities is given on the basis of a systematic conformational analysis employing various methods of ab initio MO theory (HF/6-31G*, B3LYP/6-31G*, PCM//HF/6-31G*). A wide variety of novel helical structures with hydrogen-bonded pseudocycles of different size are predicted. Since a delta-amino acid constituent may replace a dipeptide unit in alpha-peptides, there are close relationships between the secondary structures of peptides with delta-amino acid residues and typical secondary structures of alpha-peptides. However, the preference of gauche conformations at the central C(beta)-C(gamma) bonds of delta-amino acids, which correspond to the peptide linkages in alpha-peptides, over staggered ones makes completely novel structure alternatives for helices and turns more probable. The peculiarities of beta-turn formation by sugar amino acids derived from delta-amino acids are compared with the turn formation in delta-amino acid residues and in alpha-peptides. The considerable potential of secondary structure formation in delta-peptides and single delta-amino acid constituents predicted by ab initio MO theory may stimulate experimental work in the field of peptide and foldamer design.     

A photocontrolled beta-hairpin peptide

Beta-hairpins constitute the smallest beta-type structures in peptides and proteins. The development of highly stable, yet monomeric beta-hairpins based on the tryptophan zipper motif was therefore a remarkable success [A. G. Cochran, N. J. Skelton, M. A. Starovasnik, Proc. Natl. Acad. Sci USA 2001, 98, 5578-5583]. We have been able to design, synthesize and characterize a hairpin based on this motif which incorporates an azobenzene-based photoswitch, that allows for time-resolved folding studies of beta-structures with unprecedented time resolution. At room temperature the trans-azo isomer exhibits a mostly disordered structure; however, light-induced isomerization to the cis-azo form leads to a predominantly extended and parallel conformation of the two peptide parts, which are linked by the novel photoswitch, [3-(3-aminomethyl)phenylazo]phenylacetic acid (AMPP). While in the original sequence the dipeptide Asn-Gly forms a type I' beta-turn which connects the two strands of the hairpin, this role is adopted by the AMPP chromophore in our photoresponsive beta-hairpin that can apparently act as a beta I'-turn mimetic. The beta-hairpin structure was determined and confirmed by NMR spectroscopy, but the folding process can be monitored by pronounced changes in the CD, IR and fluorescence spectra. Finally, incorporation of the structurally and functionally important beta-hairpin motif into proteins by chemical ligation might allow for the photocontrol of protein structures and/or functions.