Stability issues of covalently and noncovalently bonded peptide subunits

The present study focuses on important questions associated with modeling of peptide and protein stability. Computing at different levels of theory (RHF, B3LYP) for a representative ensemble of conformers of di‐ and tripeptides of alanine, we found that the Gibbs Free Energy values correlate significantly with the total electronic energy of the molecules (0.922 ≤ R2). For noncovalently attached but interacting peptide subunits, such as [For‐NH2]2 or [For‐L‐Ala‐NH2]2, we have found, as expected, that the basis set superimposition error (BSSE) is large in magnitude for small basis set but significantly smaller when larger basis sets [e.g., B3LYP/6‐311++G(d,p)] are used. Stability of the two hydrogen bonds of antiparallel β‐pleated sheets were quantitatively determined as a function of the molecular structure, S10 and S14, computed as 4.0 ± 0.5 and 8.1 ± 1.1 kcal/mol, respectively. Finally, a suitable thermoneutral isodesmic reaction was introduced to scale both covalently and noncovalently attached peptide units onto a common stability scale. We found that a suitable isodesmic reaction can result in the total electronic energy as well as the Gibbs free energy of a molecule, from its “noninteracting” fragments, as accurate as a few tenths of a kcal per mol. The latter observation seems to hold for peptides regardless of their length (1 ≤ n ≤ 8) or the level of theory applied. © 2004 Wiley Periodicals, Inc. J Comput Chem 25: 1084–1100, 2004

[1]  ANDRÁS PERCZEL,et al.  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 , 2003, J. Comput. Chem..

[2]  T. Keiderling,et al.  Differentiation of β-Sheet-Forming Structures: Ab Initio-Based Simulations of IR Absorption and Vibrational CD for Model Peptide and Protein β-Sheets , 2001 .

[3]  S. F. Boys,et al.  The calculation of small molecular interactions by the differences of separate total energies. Some procedures with reduced errors , 1970 .

[4]  P. Karlson,et al.  Biochemical nomenclature , 1978, Nature.

[5]  James C. Weisshaar,et al.  Conformational Analysis of Alanine Dipeptide from Dipolar Couplings in a Water-Based Liquid Crystal , 2003 .

[6]  D. Beveridge,et al.  Exploratory studies of ab initio protein structure prediction: Multiple copy simulated annealing, AMBER energy functions, and a generalized born/solvent accessibility solvation model , 2002, Proteins.

[7]  Young Kee Kang,et al.  Ab initio MO and density functional studies on trans and cis conformers of N-methylacetamide , 2001 .

[8]  David A. Dixon,et al.  Conformational Study of the Alanine Dipeptide at the MP2 and DFT Levels , 2002 .

[9]  H. Hofmann,et al.  Stability Order of Basic Peptide Conformations Reflected by Density Functional Theory , 1998 .

[10]  D. Weaver,et al.  A data mining and ab initio study of the interaction between the aromatic and backbone amide groups in proteins , 2000 .

[11]  G R Marshall,et al.  Ab initio modeling of small, medium, and large loops in proteins. , 2001, Biopolymers.

[12]  Ching-Hsing Yu,et al.  Local Geometry Trends and Torsional Sensitivity in N-Formyl-l-alanyl-l-alanine Amide and the Limitations of the Dipeptide Approximation , 1999 .

[13]  O. Olsen,et al.  A simple and realistic model system for studying hydrogen bonds in β-sheets , 2003 .

[14]  A. Perczel,et al.  Peptide models 6. New .beta.-turn conformations from ab initio calculations confirmed by x-ray data of proteins , 1993 .

[15]  P. Pulay.  Convergence acceleration of iterative sequences. the case of scf iteration , 1980 .

[16]  Peter Pulay,et al.  Geometry optimization by direct inversion in the iterative subspace , 1984 .

[17]  W. Goddard,et al.  Stabilization of α-Helices by Dipole−Dipole Interactions within α-Helices , 2000 .

[18]  Ranbir Singh,et al.  J. Mol. Struct. (Theochem) , 1996 .

[19]  Organon Scientific.  Commission on Biochemical Nomenclature , 1987 .

[20]  Harry A. Stern,et al.  Strength of the N-HiiiO=C and C-HiiiO=C Bonds in Formamide and N-Methylacetamide Dimers , 2001 .

[21]  IUPAC-IUB commission on biochemical nomenclature , 2005, Molecular and Cellular Biochemistry.

[22]  H. Hofmann,et al.  Structural and energetic relations between β turns , 1997 .

[23]  I. Jákli,et al.  PEPTIDE MODELS XXII. A CONFORMATIONAL MODEL FOR AROMATIC AMINO ACID RESIDUES IN PROTEINS. A COMPREHENSIVE ANALYSIS OF ALL THE RHF/6-31+G CONFORMERS OF FOR-L-PHE-NH2 , 1998 .

[24]  Ching-Hsing Yu,et al.  Ab initio conformational analysis of N-formyl l-alanine amide including electron correlation☆ , 2001 .

[25]  M. Gerhards,et al.  Structures of the protected amino acid Ac–Phe–OMe and its dimer: A β-sheet model system in the gas phase , 2002 .

[26]  H. Schaefer,et al.  Extensive theoretical studies of the hydrogen‐bonded complexes (H2O)2, (H2O)2H+, (HF)2, (HF)2H+, F2H−, and (NH3)2 , 1986 .

[27]  G. N. Ramachandran,et al.  Stereochemical criteria for polypeptide and protein chain conformations. II. Allowed conformations for a pair of peptide units. , 1965, Biophysical journal.

[28]  David A. Dixon,et al.  Amide-Water and Amide-Amide Hydrogen Bond Strengths , 1994 .

[29]  P. Bouř,et al.  Ab initio quantum mechanical models of peptide helices and their vibrational spectra. , 2002, Biopolymers.

[30]  P. Pulay.  Improved SCF convergence acceleration , 1982 .

[31]  Walter Thiel,et al.  Description of peptide and protein secondary structures employing semiempirical methods , 2001 .

[32]  Imre G. Csizmadia,et al.  Peptide models. 1. Topology of selected peptide conformational potential energy surfaces (glycine and alanine derivatives) , 1991 .

[33]  I. Jákli,et al.  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 , 2000 .

[34]  Donald G. Truhlar,et al.  Systematic study of basis set superposition errors in the calculated interaction energy of two HF molecules , 1985 .

[35]  Miquel Duran,et al.  How does basis set superposition error change the potential surfaces for hydrogen-bonded dimers? , 1996 .

[36]  R. L. Baldwin.  Relation between peptide backbone solvation and the energetics of peptide hydrogen bonds. , 2002, Biophysical chemistry.

[37]  I. C. O. B. Nomenclature.  IUPAC-IUB Commission on Biochemical Nomenclature. Abbreviations and symbols for the description of the conformation of polypeptide chains. Tentative rules (1969). , 1970, Biochemistry.