Parameterization of PACE Force Field for Membrane Environment and Simulation of Helical Peptides and Helix-Helix Association.

The recently developed PACE force field was further parametrized so that it can be applied to the studies of membrane systems. Parameters for the interactions between united-atom protein particles and lipid hydrophobic tails were developed by reproducing the solvation free energies of small organic molecules in hexadecane. Interactions between protein particles and lipid heads were parametrized by fitting the potential of mean force of the corresponding all-atom simulation. The force field was applied to the study of five helical peptides in membrane environments. The calculated tilt angles of WALP and GWALP and their mutations are in good agreement with experimental data. The association of two glycophorin A (GpA) helices was simulated for 6 μs. Root-mean-square-deviation of the simulated dimer from the nuclear magnetic resonance structure was found to be 0.272 nm, better than all results obtained so far. These findings demonstrate the high accuracy and applicability of the PACE force field in studying membrane proteins.

[1]  Roman G. Efremov,et al.  Association of transmembrane helices: what determines assembling of a dimer? , 2006, J. Comput. Aided Mol. Des..

[2]  Siewert J. Marrink,et al.  Release of content through mechano-sensitive gates in pressurized liposomes , 2010, Proceedings of the National Academy of Sciences.

[3]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[4]  Maarten G. Wolf,et al.  g_membed: Efficient insertion of a membrane protein into an equilibrated lipid bilayer with minimal perturbation , 2010, J. Comput. Chem..

[5]  James H. Prestegard,et al.  A Transmembrane Helix Dimer: Structure and Implications , 1997, Science.

[6]  D. Engelman,et al.  Detergents modulate dimerization, but not helicity, of the glycophorin A transmembrane domain. , 1999, Journal of molecular biology.

[7]  J. Killian,et al.  Interfacial anchor properties of tryptophan residues in transmembrane peptides can dominate over hydrophobic matching effects in peptide-lipid interactions. , 2003, Biochemistry.

[8]  Luís M. S. Loura,et al.  Is there a preferential interaction between cholesterol and tryptophan residues in membrane proteins? , 2008, Biochemistry.

[9]  Yundong Wu,et al.  Coarse-Grained Protein Model Coupled with a Coarse-Grained Water Model:  Molecular Dynamics Study of Polyalanine-Based Peptides. , 2007, Journal of chemical theory and computation.

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

[11]  Michael L Klein,et al.  Probing Membrane Insertion Activity of Antimicrobial Polymers via Coarse-grain Molecular Dynamics. , 2006, Journal of chemical theory and computation.

[12]  Yun-Dong Wu,et al.  Toward a Coarse-Grained Protein Model Coupled with a Coarse-Grained Solvent Model: Solvation Free Energies of Amino Acid Side Chains. , 2008, Journal of chemical theory and computation.

[13]  G. Lorigan,et al.  Determining the helical tilt angle of a transmembrane helix in mechanically aligned lipid bilayers using EPR spectroscopy. , 2007, Journal of the American Chemical Society.

[14]  Paj Peter Hilbers,et al.  Lipid-based mechanisms for vesicle fission. , 2007, The journal of physical chemistry. B.

[15]  B. Maigret,et al.  Self assembly of peptides near or within membranes using coarse grained MD simulations , 2009 .

[16]  W. DeGrado,et al.  Association of a model transmembrane peptide containing gly in a heptad sequence motif. , 2004, Biophysical journal.

[17]  G. von Heijne,et al.  Membrane protein structure: prediction versus reality. , 2007, Annual review of biochemistry.

[18]  C. Brooks,et al.  Interfacial folding and membrane insertion of designed peptides studied by molecular dynamics simulations. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[19]  A. Pohorille,et al.  Insights into the recognition and association of transmembrane α-helices. The free energy of α-helix dimerization in glycophorin A , 2005 .

[20]  D. Engelman,et al.  Lipid bilayer thickness varies linearly with acyl chain length in fluid phosphatidylcholine vesicles. , 1983, Journal of molecular biology.

[21]  D. Tieleman,et al.  The MARTINI force field: coarse grained model for biomolecular simulations. , 2007, The journal of physical chemistry. B.

[22]  A. Kukol,et al.  Systematic molecular dynamics searching in a lipid bilayer: application to the glycophorin A and oncogenic ErbB-2 transmembrane domains. , 2006, Journal of molecular graphics & modelling.

[23]  G. Heijne.  Membrane-protein topology , 2006, Nature Reviews Molecular Cell Biology.

[24]  D. Engelman,et al.  Membrane protein folding and oligomerization: the two-stage model. , 1990, Biochemistry.

[25]  M. Sansom,et al.  Transmembrane helix-helix interactions: comparative simulations of the glycophorin a dimer. , 2006, Biochemistry.

[26]  Michael R. Shirts,et al.  Extremely precise free energy calculations of amino acid side chain analogs: Comparison of common molecular mechanics force fields for proteins , 2003 .

[27]  J. Bowie.  Solving the membrane protein folding problem , 2005, Nature.

[28]  D. Engelman,et al.  Helical membrane protein folding, stability, and evolution. , 2000, Annual review of biochemistry.

[29]  Yun-Dong Wu,et al.  PACE Force Field for Protein Simulations. 1. Full Parameterization of Version 1 and Verification. , 2010, Journal of chemical theory and computation.

[30]  Wonpil Im,et al.  Membrane assembly of simple helix homo-oligomers studied via molecular dynamics simulations. , 2007, Biophysical journal.

[31]  R. Larson,et al.  Molecular dynamics simulations of model trans-membrane peptides in lipid bilayers: a systematic investigation of hydrophobic mismatch. , 2006, Biophysical journal.

[32]  James U Bowie,et al.  A simple method for modeling transmembrane helix oligomers. , 2003, Journal of molecular biology.

[33]  Berk Hess,et al.  LINCS: A linear constraint solver for molecular simulations , 1997 .

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

[35]  Peter L. Freddolino,et al.  Ten-microsecond molecular dynamics simulation of a fast-folding WW domain. , 2008, Biophysical journal.

[36]  Graham R. Smith,et al.  Setting up and optimization of membrane protein simulations , 2002, European Biophysics Journal.

[37]  H. Fritz,et al.  Dimerisation of the glycophorin A transmembrane segment in membranes probed with the ToxR transcription activator. , 1996, Journal of molecular biology.

[38]  J. Salgado,et al.  The dynamic orientation of membrane-bound peptides: bridging simulations and experiments. , 2007, Biophysical journal.

[39]  D. Engelman,et al.  Effect of detergents on the association of the glycophorin a transmembrane helix. , 2003, Biophysical journal.

[40]  Mutagenesis data in the automated prediction of transmembrane helix dimers , 2007, Proteins.

[41]  B. L. de Groot,et al.  Progress in the analysis of membrane protein structure and function , 2002, FEBS letters.

[42]  D. Engelman,et al.  The glycophorin A transmembrane domain dimer: sequence-specific propensity for a right-handed supercoil of helices. , 1992, Biochemistry.

[43]  Joseph A. Bank,et al.  Supporting Online Material Materials and Methods Figs. S1 to S10 Table S1 References Movies S1 to S3 Atomic-level Characterization of the Structural Dynamics of Proteins , 2022 .

[44]  K. Fleming.  Standardizing the free energy change of transmembrane helix-helix interactions. , 2002, Journal of molecular biology.

[45]  R. Larson,et al.  The MARTINI Coarse-Grained Force Field: Extension to Proteins. , 2008, Journal of chemical theory and computation.

[46]  A. Holt,et al.  Order parameters of a transmembrane helix in a fluid bilayer: case study of a WALP peptide. , 2010, Biophysical journal.

[47]  Wonpil Im,et al.  Revisiting hydrophobic mismatch with free energy simulation studies of transmembrane helix tilt and rotation. , 2010, Biophysical journal.

[48]  M. Sansom,et al.  Changes in transmembrane helix alignment by arginine residues revealed by solid-state NMR experiments and coarse-grained MD simulations. , 2010, Journal of the American Chemical Society.

[49]  Christian Kandt,et al.  Setting up and running molecular dynamics simulations of membrane proteins. , 2007, Methods.

[50]  Berend Smit,et al.  Simulation studies of protein-induced bilayer deformations, and lipid-induced protein tilting, on a mesoscopic model for lipid bilayers with embedded proteins. , 2005, Biophysical journal.

[51]  Justin L. MacCallum,et al.  Partitioning of Amino Acid Side Chains into Lipid Bilayers: Results from Computer Simulations and Comparison to Experiment , 2007, The Journal of general physiology.

[52]  N. Kučerka,et al.  Influence of cholesterol on the bilayer properties of monounsaturated phosphatidylcholine unilamellar vesicles , 2007, The European physical journal. E, Soft matter.

[53]  Peter J Bond,et al.  Insertion and assembly of membrane proteins via simulation. , 2006, Journal of the American Chemical Society.

[54]  S. O. Smith,et al.  Structure of the transmembrane dimer interface of glycophorin A in membrane bilayers. , 2001, Biochemistry.

[55]  R. C. Reeder,et al.  A Coarse Grain Model for Phospholipid Simulations , 2001 .

[56]  Michael H. Abraham,et al.  Thermodynamics of solute transfer from water to hexadecane , 1990 .

[57]  D. Engelman,et al.  Modulation of glycophorin A transmembrane helix interactions by lipid bilayers: molecular dynamics calculations. , 2000, Journal of molecular biology.

[58]  D. van der Spoel,et al.  GROMACS: A message-passing parallel molecular dynamics implementation , 1995 .

[59]  Durba Sengupta,et al.  Polarizable Water Model for the Coarse-Grained MARTINI Force Field , 2010, PLoS Comput. Biol..

[60]  F. Jiang,et al.  Influence of side chain conformations on local conformational features of amino acids and implication for force field development. , 2010, The journal of physical chemistry. B.

[61]  S. Vishveshwara,et al.  Geometry of proline-containing alpha-helices in proteins. , 2009, International journal of peptide and protein research.

[62]  Reinhard Lipowsky,et al.  Solvent-exposed tails as prestalk transition states for membrane fusion at low hydration. , 2010, Journal of the American Chemical Society.

[63]  Philip W Fowler,et al.  Helix-helix interactions in membrane proteins: coarse-grained simulations of glycophorin a helix dimerization. , 2008, Biochemistry.

[64]  J. Killian,et al.  Phospholipid flop induced by transmembrane peptides in model membranes is modulated by lipid composition. , 2003, Biochemistry.

[65]  Alessandro Senes,et al.  Folding of helical membrane proteins: the role of polar, GxxxG-like and proline motifs. , 2004, Current opinion in structural biology.

[66]  Gregory D. Hawkins,et al.  Extension of the platform of applicability of the SM5.42R universal solvation model , 1999 .

[67]  Peter J Bond,et al.  Bilayer deformation by the Kv channel voltage sensor domain revealed by self-assembly simulations , 2007, Proceedings of the National Academy of Sciences.

[68]  Berend Smit,et al.  Molecular simulations of lipid-mediated protein-protein interactions. , 2008, Biophysical journal.

[69]  Siewert J Marrink,et al.  Hybrid simulations: combining atomistic and coarse-grained force fields using virtual sites. , 2011, Physical chemistry chemical physics : PCCP.

[70]  Berk Hess,et al.  Improving efficiency of large time‐scale molecular dynamics simulations of hydrogen‐rich systems , 1999, Journal of computational chemistry.

[71]  A F Smeijers,et al.  Coarse-grained transmembrane proteins: hydrophobic matching, aggregation, and their effect on fusion. , 2006, The journal of physical chemistry. B.

[72]  M. Sansom,et al.  Molecular dynamics simulations of the dimerization of transmembrane alpha-helices. , 2010, Accounts of chemical research.

[73]  J. Killian,et al.  Hydrophobic mismatch between helices and lipid bilayers. , 2003, Biophysical journal.

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

[75]  K. Schulten,et al.  Molecular Dynamics Simulations of Micelle Formation around Dimeric Glycophorin A Transmembrane Helices. , 2004, Biophysical journal.

[76]  J. Killian,et al.  On the orientation of a designed transmembrane peptide: toward the right tilt angle? , 2007, Journal of the American Chemical Society.

[77]  Richard Wolfenden,et al.  Comparing the polarities of the amino acids: side-chain distribution coefficients between the vapor phase, cyclohexane, 1-octanol, and neutral aqueous solution , 1988 .

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

[79]  F. Cordes,et al.  Proline-induced distortions of transmembrane helices. , 2002, Journal of molecular biology.

[80]  A. Holt,et al.  Orientation and dynamics of transmembrane peptides: the power of simple models , 2009, European Biophysics Journal.

[81]  W F Drew Bennett,et al.  Distribution of amino acids in a lipid bilayer from computer simulations. , 2008, Biophysical journal.

[82]  J. Killian,et al.  Influence of flanking residues on tilt and rotation angles of transmembrane peptides in lipid bilayers. A solid-state 2H NMR study. , 2005, Biochemistry.

[83]  Yun-Dong Wu,et al.  PACE Force Field for Protein Simulations. 2. Folding Simulations of Peptides. , 2010, Journal of chemical theory and computation.

[84]  D. Engelman,et al.  Introduction to the membrane protein reviews: the interplay of structure, dynamics, and environment in membrane protein function. , 2006, Annual review of biochemistry.

[85]  A. Holt,et al.  Tilt and rotation angles of a transmembrane model peptide as studied by fluorescence spectroscopy. , 2009, Biophysical journal.

[86]  R. Koeppe,et al.  Charged or Aromatic Anchor Residue Dependence of Transmembrane Peptide Tilt* , 2010, The Journal of Biological Chemistry.

[87]  J. Bowie,et al.  Method to measure strong protein–protein interactions in lipid bilayers using a steric trap , 2010, Proceedings of the National Academy of Sciences.

[88]  M. Gerstein,et al.  Statistical analysis of amino acid patterns in transmembrane helices: the GxxxG motif occurs frequently and in association with beta-branched residues at neighboring positions. , 2000, Journal of molecular biology.

[89]  Huan-Xiang Zhou,et al.  A method to determine dielectric constants in nonhomogeneous systems: application to biological membranes. , 2008, Biophysical journal.

[90]  R. Friesner,et al.  Evaluation and Reparametrization of the OPLS-AA Force Field for Proteins via Comparison with Accurate Quantum Chemical Calculations on Peptides† , 2001 .

[91]  Justin L. MacCallum,et al.  Calculation of the water–cyclohexane transfer free energies of neutral amino acid side‐chain analogs using the OPLS all‐atom force field , 2003, J. Comput. Chem..

[92]  K. Schulten,et al.  Computational studies of membrane channels. , 2004, Structure.

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

[94]  D. Tieleman,et al.  Interpretation of 2H-NMR experiments on the orientation of the transmembrane helix WALP23 by computer simulations. , 2010, Biophysical journal.

[95]  Klaus Schulten,et al.  Molecular dynamics simulations of proteins in lipid bilayers. , 2005, Current opinion in structural biology.

[96]  Klaus Schulten,et al.  Coarse grained protein-lipid model with application to lipoprotein particles. , 2006, The journal of physical chemistry. B.

[97]  Y. Okamoto,et al.  Prediction of membrane protein structures by replica-exchange Monte Carlo simulations: case of two helices. , 2004, The Journal of chemical physics.

[98]  Erik Strandberg,et al.  Tilt angles of transmembrane model peptides in oriented and non-oriented lipid bilayers as determined by 2H solid-state NMR. , 2004, Biophysical journal.

[99]  M Mezei,et al.  Free Energy Simulations a , 1986, Annals of the New York Academy of Sciences.

[100]  Gregory A Voth,et al.  Solvent-free lipid bilayer model using multiscale coarse-graining. , 2009, The journal of physical chemistry. B.

[101]  R. Koeppe,et al.  Influence of proline upon the folding and geometry of the WALP19 transmembrane peptide. , 2009, Biochemistry.

[102]  Qiang Shi,et al.  Mixed atomistic and coarse-grained molecular dynamics: simulation of a membrane-bound ion channel. , 2006, The journal of physical chemistry. B.

[103]  J. Salgado,et al.  Orientational landscapes of peptides in membranes: prediction of (2)H NMR couplings in a dynamic context. , 2009, Biochemistry.

[104]  V. Cherezov,et al.  Transmembrane peptides stabilize inverted cubic phases in a biphasic length-dependent manner: implications for protein-induced membrane fusion. , 2006, Biophysical journal.

[105]  J Wang,et al.  Structure of the transmembrane region of the M2 protein H+ channel , 2001, Protein science : a publication of the Protein Society.

[106]  A. Ulrich,et al.  Orientation and dynamics of peptides in membranes calculated from 2H-NMR data. , 2009, Biophysical journal.

[107]  J. Killian,et al.  Self-association of transmembrane alpha-helices in model membranes: importance of helix orientation and role of hydrophobic mismatch. , 2005, The Journal of biological chemistry.

[108]  Durba Sengupta,et al.  Lipid-mediated interactions tune the association of glycophorin A helix and its disruptive mutants in membranes. , 2010, Physical chemistry chemical physics : PCCP.

[109]  J. Killian,et al.  Sensitivity of single membrane-spanning alpha-helical peptides to hydrophobic mismatch with a lipid bilayer: effects on backbone structure, orientation, and extent of membrane incorporation. , 2001, Biochemistry.

[110]  D. Engelman,et al.  The GxxxG motif: a framework for transmembrane helix-helix association. , 2000, Journal of molecular biology.

[111]  D. Langosch,et al.  The dimerization motif of the glycophorin A transmembrane segment in membranes: Importance of glycine residues , 1998, Protein science : a publication of the Protein Society.