Abstract
The Conformational Free Energy Thermodynamic Integration (CFTI) method, a new multidimensional approach for conformational free energy simulations, is presented. The method is applied to two problems of biochemical interest: exploration of the free energy surfaces of helical alanine (Ala) and α-methylalanine (Aib) homopeptides in vacuum and the cost of pre-organization of the opioid peptide Tyr-D-Pen-Gly-Phe-D-Pen (DPDPE) peptide for disulfide bond formation. In the CFTI approach a single molecular dynamics simulation with all ø and ψ dihedrals kept fixed yields the complete conformational free energy gradient for the studied peptides. For regular structures of model peptides (Ala)n and (Aib)n where n=6,8,10 and Aib is a-methylalanine in vacuum, free energy maps in the helical region of ø- ψ space are calculated, and used to roughly locate stable states. The locations of the free energy minima are further refined by the novel procedure of free energy optimization by steepest descent down the gradient, leading to structures in excellent agreement with experimental data. The stability of the minima with respect to deformations is studied by analysis of second derivatives of the free energy surface. Analysis of free energy components and molecular structures uncovers the molecular mechanism for the propensity of Aib peptides for the 310-helix structure in the interplay between the quality and quantity of hydrogen bonds. For the linear form of the DPDPE peptide in solution, free energy differences are calculated between four conformers: Cyc, representing the structure adopted by the linear peptide prior to disulfide bond formation, β c and β E, two slightly different β-turns previously identified as representative, stable structures of the peptide, and Ext, an extended structure. The simulations indicate that β E is the most stable of the studied conformers of linear DPDPE in aqueous solution, with β c , Cyc and Ext having free energies higher by 2.3, 6.3, and 28.2 kcal/mol, respectively. The free energy differences of 4.0 kcal/mol between β c and Cyc, and 6.3 kcal/mol between β E and Cyc, reflect the cost of pre-organizing the linear peptide into a conformation conducive for disulfide bond formation. Such a conformational change is a pre-requisite for the chemical reaction of S-S bond formation to proceed.
Access this chapter
Tax calculation will be finalised at checkout
Purchases are for personal use only
Preview
Unable to display preview. Download preview PDF.
Similar content being viewed by others
References
T. P. Straatsma and J. A. McCammon. Computational alchemy. Annu. Rev. Phys. Chem., 43: 407–435, 1992.
K. Kuczera. One-and multidimensional conformational free energy simulations. J. Comp. Chem., 17: 1726–1749, 1996.
T.P. Straatsma and J.A. McCammon. Multiconfiguration thermodynamic integration. J. Chem. Phys., 95: 1175–1188, 1991.
P. E. Smith and B. M. Pettitt. Amino acid side-chain populations in aqueous and saline solution: Bis(penicillamine) enkephalin. Biopolymers, 32: 1623–1629, 1992.
D. E. Smith and A. D. J. Haymet. Free energy, entropy and internal energy of hydrophobic interactions: computer simulations. J. Chem. Phys., 98: 6445–6454, 1993.
E. M. Boczko and C. L. Brooks III. First principles calculation of the folding free energy of a three helix bundle protein. Science, 269: 393–396, 1995.
M. Mezei and D. L. Beveridge. Free energy simulations. Ann. N. Y. Acad. Sci., 482: 1–23, 1986.
Y. Wang and K. Kuczera. Exploration of conformational free energy surface of helical Alan and Aibn peptides. J. Phys. Chem. B, 101: 5205–5213, 1997.
Y. Wang. Computational Approach to the Influence of the Bisulfide Bond on Peptide Properties. PhD thesis, University of Kansas, 1997.
Y. Wang and K. Kuczera. Conformational free energy surface of the linear DPDPE peptide: Cost of pre-organization for disulfide bond formation. J. Am. Chem. Soc, submitted, 1997.
H. I. Mosberg, R. Hurst, V. J. Hruby, K. Gee, H. I. Yamamura, J. J. Galligan, and T. F. Burks. Bis-penicillamine enkephalins possess highly improved specificity toward δ opioid receptors. Proc. Natl. Acad. Sci. U.S.A., 80: 5871–5874, 1983.
E. A. Carter, G. Ciccotti, J. T. Hynes, and R. Kapral. Constrained reaction coordinate dynamics for the simulation of rare events. Chem. Phys. Lett, 156: 472–477, 1989.
C. L. Brooks III, M. Karplus, and B. M. Petitt. Proteins: A Theoretical Perspective of Dynamics, Structure, and Thermodynamics. John Wiley and Sons, New York, 1988.
Y. Wang and K. Kuczera. Molecular dynamics simulations of cyclic and linear DPDPE: Influence of the disulfide bond on peptide flexibility. J. Phys. Chem., 100: 2555–2563, 1996.
J. Gao, K. Kuczera, B. Tidor, and M. Karplus. Hidden thermodynamics of mutant proteins: A molecular dynamics analysis. Science, 244: 1069–1072, 1989.
K. Kuczera, J. Gao, B. Tidor, and M. Karplus. Free energy of sickling: A simulation analysis. Proc. Natl. Acad. Sci. USA, 87: 8481–8485, 1990.
Author information
Authors and Affiliations
Editor information
Editors and Affiliations
Rights and permissions
Copyright information
© 1999 Springer-Verlag Berlin Heidelberg
About this paper
Cite this paper
Kuczera, K. (1999). Exploration of Peptide Free Energy Surfaces. In: Deuflhard, P., Hermans, J., Leimkuhler, B., Mark, A.E., Reich, S., Skeel, R.D. (eds) Computational Molecular Dynamics: Challenges, Methods, Ideas. Lecture Notes in Computational Science and Engineering, vol 4. Springer, Berlin, Heidelberg. https://doi.org/10.1007/978-3-642-58360-5_9
Download citation
DOI: https://doi.org/10.1007/978-3-642-58360-5_9
Publisher Name: Springer, Berlin, Heidelberg
Print ISBN: 978-3-540-63242-9
Online ISBN: 978-3-642-58360-5
eBook Packages: Springer Book Archive