Efficient molecular mechanics simulations of the folding, orientation, and assembly of peptides in lipid bilayers using an implicit atomic solvation model

Membrane proteins comprise a significant fraction of the proteomes of sequenced organisms and are the targets of approximately half of marketed drugs. However, in spite of their prevalence and biomedical importance, relatively few experimental structures are available due to technical challenges. Computational simulations can potentially address this deficit by providing structural models of membrane proteins. Solvation within the spatially heterogeneous membrane/solvent environment provides a major component of the energetics driving protein folding and association within the membrane. We have developed an implicit solvation model for membranes that is both computationally efficient and accurate enough to enable molecular mechanics predictions for the folding and association of peptides within the membrane. We derived the new atomic solvation model parameters using an unbiased fitting procedure to experimental data and have applied it to diverse problems in order to test its accuracy and to gain insight into membrane protein folding. First, we predicted the positions and orientations of peptides and complexes within the lipid bilayer and compared the simulation results with solid-state NMR structures. Additionally, we performed folding simulations for a series of host–guest peptides with varying propensities to form alpha helices in a hydrophobic environment and compared the structures with experimental measurements. We were also able to successfully predict the structures of amphipathic peptides as well as the structures for dimeric complexes of short hexapeptides that have experimentally characterized propensities to form beta sheets within the membrane. Finally, we compared calculated relative transfer energies with data from experiments measuring the effects of mutations on the free energies of translocon-mediated insertion of proteins into lipid bilayers and of combined folding and membrane insertion of a beta barrel protein.

[1]  Trevor Hastie,et al.  The Elements of Statistical Learning , 2001 .

[2]  Eric R. Ziegel,et al.  The Elements of Statistical Learning , 2003, Technometrics.

[3]  H. Scheraga,et al.  Conformational studies of poly‐L‐alanine in water , 1968, Biopolymers.

[4]  R. Hodges,et al.  Relationship of sidechain hydrophobicity and alpha-helical propensity on the stability of the single-stranded amphipathic alpha-helix. , 1995, Journal of peptide science : an official publication of the European Peptide Society.

[5]  Paul A. Wiggins,et al.  Emerging roles for lipids in shaping membrane-protein function , 2009, Nature.

[6]  P. Y. Chou,et al.  Conformational parameters for amino acids in helical, beta-sheet, and random coil regions calculated from proteins. , 1974, Biochemistry.

[7]  K. H. Kalk,et al.  Structural evidence for dimerization-regulated activation of an integral membrane phospholipase. , 1999 .

[8]  P M Cullis,et al.  Affinities of amino acid side chains for solvent water. , 1981, Biochemistry.

[9]  Per Lincoln,et al.  Tryptophan orientation in model lipid membranes. , 2007, Biochemical and biophysical research communications.

[10]  W. Yau,et al.  Electrostatic contributions to indole-lipid interactions. , 2005, The journal of physical chemistry. B.

[11]  Alessandro Senes,et al.  E(z), a depth-dependent potential for assessing the energies of insertion of amino acid side-chains into membranes: derivation and applications to determining the orientation of transmembrane and interfacial helices. , 2007, Journal of molecular biology.

[12]  Velin Z. Spassov,et al.  Introducing an Implicit Membrane in Generalized Born/Solvent Accessibility Continuum Solvent Models , 2002 .

[13]  Nathaniel J. Traaseth,et al.  Structure and topology of monomeric phospholamban in lipid membranes determined by a hybrid solution and solid-state NMR approach , 2009, Proceedings of the National Academy of Sciences.

[14]  R. L. Baldwin,et al.  Helix propensities of the amino acids measured in alanine‐based peptides without helix‐stabilizing side‐chain interactions , 1994, Protein science : a publication of the Protein Society.

[15]  Iupaciubcommissiononbiochemic IUPAC-IUB Commission on Biochemical Nomenclature. Abbreviations and symbols for the description of the conformation of polypeptide chains. , 1971, Journal of molecular biology.

[16]  R G Efremov,et al.  A solvent model for simulations of peptides in bilayers. I. Membrane-promoting alpha-helix formation. , 1999, Biophysical journal.

[17]  Yves F Dufrêne,et al.  Antiparallel beta-sheet: a signature structure of the oligomeric amyloid beta-peptide. , 2009, The Biochemical journal.

[18]  R. Read,et al.  Atomic solvation parameters in the analysis of protein‐protein docking results , 1995, Protein science : a publication of the Protein Society.

[19]  H. Scheraga,et al.  Accessible surface areas as a measure of the thermodynamic parameters of hydration of peptides. , 1987, Proceedings of the National Academy of Sciences of the United States of America.

[20]  P S Kim,et al.  Context is a major determinant of beta-sheet propensity. , 1994, Nature.

[21]  J. Berg,et al.  Thermodynamic beta-sheet propensities measured using a zinc-finger host peptide. , 1993, Nature.

[22]  H. Scheraga,et al.  Energy parameters in polypeptides. 10. Improved geometrical parameters and nonbonded interactions for use in the ECEPP/3 algorithm, with application to proline-containing peptides , 1994 .

[23]  S. Lifson,et al.  On the Theory of Helix—Coil Transition in Polypeptides , 1961 .

[24]  C. Deber,et al.  Uncoupling Hydrophobicity and Helicity in Transmembrane Segments , 1998, The Journal of Biological Chemistry.

[25]  A. Waring,et al.  Impact of single-residue mutations on the structure and function of ovispirin/novispirin antimicrobial peptides. , 2002, Protein engineering.

[26]  P S Kim,et al.  Measurement of the beta-sheet-forming propensities of amino acids. , 1994, Nature.

[27]  Marco Punta,et al.  The New York Consortium on Membrane Protein Structure (NYCOMPS): a high-throughput platform for structural genomics of integral membrane proteins , 2010, Journal of Structural and Functional Genomics.

[28]  P. S. Kim,et al.  Measurement of the β-sheet-forming propensities of amino acids , 1994, Nature.

[29]  Stephen H. White,et al.  Experimentally determined hydrophobicity scale for proteins at membrane interfaces , 1996, Nature Structural Biology.

[30]  S H White,et al.  Folding of beta-sheet membrane proteins: a hydrophobic hexapeptide model. , 1998, Journal of molecular biology.

[31]  Ruben Abagyan,et al.  ICM—A new method for protein modeling and design: Applications to docking and structure prediction from the distorted native conformation , 1994, J. Comput. Chem..

[32]  A. Barabasi,et al.  Drug—target network , 2007, Nature Biotechnology.

[33]  G. Heijne,et al.  Molecular code for transmembrane-helix recognition by the Sec61 translocon , 2007, Nature.

[34]  A. D. McLachlan,et al.  Solvation energy in protein folding and binding , 1986, Nature.

[35]  C. Brooks,et al.  An implicit membrane generalized born theory for the study of structure, stability, and interactions of membrane proteins. , 2003, Biophysical journal.

[36]  Arieh Warshel,et al.  On the energetics of translocon-assisted insertion of charged transmembrane helices into membranes , 2010, Proceedings of the National Academy of Sciences.

[37]  G. Heijne,et al.  Genome‐wide analysis of integral membrane proteins from eubacterial, archaean, and eukaryotic organisms , 1998, Protein science : a publication of the Protein Society.

[38]  J. Gesell,et al.  Structures of the M2 channel-lining segments from nicotinic acetylcholine and NMDA receptors by NMR spectroscopy , 1999, Nature Structural Biology.

[39]  H. Scheraga,et al.  Energy parameters in polypeptides. 9. Updating of geometrical parameters, nonbonded interactions, and hydrogen bond interactions for the naturally occurring amino acids , 1983 .

[40]  B. Zimm,et al.  Theory of the Phase Transition between Helix and Random Coil in Polypeptide Chains , 1959 .

[41]  C. Deber,et al.  A measure of helical propensity for amino acids in membrane environments , 1994, Nature Structural Biology.

[42]  J. Gesell,et al.  Two-dimensional 1H NMR experiments show that the 23-residue magainin antibiotic peptide is an α-helix in dodecylphosphocholine micelles, sodium dodecylsulfate micelles, and trifluoroethanol/water solution , 1997, Journal of biomolecular NMR.

[43]  R. Hodges,et al.  Relationship of sidechain hydrophobicity and α‐helical propensity on the stability of the single‐stranded amphipathic α‐helix , 1995 .

[44]  E. Lindahl,et al.  Membrane proteins: molecular dynamics simulations. , 2008, Current opinion in structural biology.

[45]  Mei Hong Structure, topology, and dynamics of membrane peptides and proteins from solid-state NMR spectroscopy. , 2007, The journal of physical chemistry. B.

[46]  Andrei L Lomize,et al.  Positioning of proteins in membranes: A computational approach , 2006, Protein science : a publication of the Protein Society.

[47]  T. Lazaridis Effective energy function for proteins in lipid membranes , 2003, Proteins.

[48]  P. S. Kim,et al.  Context is a major determinant of β-sheet propensity , 1994, Nature.

[49]  Martin B Ulmschneider,et al.  Properties of integral membrane protein structures: Derivation of an implicit membrane potential , 2005, Proteins.

[50]  D. Eisenberg,et al.  Atomic solvation parameters applied to molecular dynamics of proteins in solution , 1992, Protein science : a publication of the Protein Society.

[51]  Jakob P Ulmschneider,et al.  A generalized born implicit-membrane representation compared to experimental insertion free energies. , 2007, Biophysical journal.

[52]  Robert M. Graham,et al.  Non-α-helical elements modulate polytopic membrane protein architecture , 2001 .

[53]  K. Hahm,et al.  Role of the hinge region and the tryptophan residue in the synthetic antimicrobial peptides, cecropin A(1-8)-magainin 2(1-12) and its analogues, on their antibiotic activities and structures. , 2000, Biochemistry.

[54]  H. Scheraga,et al.  Energy parameters in polypeptides. VII. Geometric parameters, partial atomic charges, nonbonded interactions, hydrogen bond interactions, and intrinsic torsional potentials for the naturally occurring amino acids , 1975 .

[55]  R. Abagyan,et al.  Biased probability Monte Carlo conformational searches and electrostatic calculations for peptides and proteins. , 1994, Journal of molecular biology.

[56]  S. White,et al.  The preference of tryptophan for membrane interfaces. , 1998, Biochemistry.

[57]  K. G. Fleming,et al.  Side-chain hydrophobicity scale derived from transmembrane protein folding into lipid bilayers , 2011, Proceedings of the National Academy of Sciences.

[58]  M. Feig,et al.  A generalized Born formalism for heterogeneous dielectric environments: application to the implicit modeling of biological membranes. , 2005, The Journal of chemical physics.

[59]  A. Mark,et al.  Coarse grained model for semiquantitative lipid simulations , 2004 .

[60]  Syma Khalid,et al.  Coarse-grained molecular dynamics simulations of membrane proteins and peptides. , 2007, Journal of structural biology.

[61]  S H White,et al.  Structure of a fluid dioleoylphosphatidylcholine bilayer determined by joint refinement of x-ray and neutron diffraction data. III. Complete structure. , 1992, Biophysical journal.

[62]  R. Abagyan,et al.  Soft protein–protein docking in internal coordinates , 2002, Protein science : a publication of the Protein Society.

[63]  S. Opella,et al.  Structure determination of a membrane protein with two trans-membrane helices in aligned phospholipid bicelles by solid-state NMR spectroscopy. , 2006, Journal of the American Chemical Society.

[64]  Portland Press Ltd IUPAC-IUB Commission on Biochemical Nomenclature. Abbreviations and symbols for the description of the conformation of polypeptide chains. Tentative rules (1969). , 1971, Biochemistry.

[65]  W. Kabsch,et al.  Dictionary of protein secondary structure: Pattern recognition of hydrogen‐bonded and geometrical features , 1983, Biopolymers.

[66]  Sara M. Butterfield,et al.  Amyloidogenic Protein—Membrane Interactions: Mechanistic Insight from Model Systems , 2010 .

[67]  Catherine Etchebest,et al.  Determining membrane protein structures: still a challenge! , 2007, Trends in biochemical sciences.

[68]  Wilfred F. van Gunsteren,et al.  Computer Simulation of Biomolecular Systems: Theoretical and Experimental Applications , 1989 .

[69]  H. Scheraga,et al.  Effect of amino acid composition on the twist and the relative stability of parallel and antiparallel .beta.-sheets , 1983 .

[70]  Roman G. Efremov,et al.  A Solvent Model for Simulations of Peptides in Bilayers. I. Membrane-Promoting α-Helix Formation , 1999 .

[71]  S. Opella,et al.  Simultaneous assignment and structure determination of a membrane protein from NMR orientational restraints , 2003, Protein science : a publication of the Protein Society.

[72]  P. Kollman,et al.  Protein structure prediction with a combined solvation free energy-molecular mechanics force field , 1993 .

[73]  W. Wimley,et al.  Folding of beta-sheets in membranes: specificity and promiscuity in peptide model systems. , 2001, Journal of molecular biology.

[74]  J. Thornton,et al.  Helix geometry in proteins. , 1988, Journal of molecular biology.

[75]  Ruben Abagyan,et al.  Protein structure prediction by global energy optimization , 1997 .

[76]  M. Karplus,et al.  Effective energy function for proteins in solution , 1999, Proteins.

[77]  B Honig,et al.  Extracting hydrophobic free energies from experimental data: relationship to protein folding and theoretical models. , 1991, Biochemistry.

[78]  J Novotny,et al.  Non-alpha-helical elements modulate polytopic membrane protein architecture. , 2001, Journal of molecular biology.

[79]  Cinque S. Soto,et al.  Structure of the Amantadine Binding Site of Influenza M2 Proton Channels In Lipid Bilayers , 2010, Nature.

[80]  Jeremy M. Berg,et al.  Thermodynamic β -sheet propensities measured using a zinc-finger host peptide , 1993, Nature.

[81]  Erik Lindahl,et al.  Protein contents in biological membranes can explain abnormal solvation of charged and polar residues , 2009, Proceedings of the National Academy of Sciences.