Parametrization of MARTINI for Modeling Hinging Motions in Membrane Proteins.

Helical hinges are common in transmembrane (TM) proteins. Hinging motions of TM helices are associated with the functional dynamics of many important membrane proteins such as various G-protein-coupled receptors and potassium channels. Although coarse-grained (CG) simulations are useful for studying membrane proteins, they are limited in describing accurately the hinge-related protein dynamics. In this study, we have overcome this limitation through a further development of the MARTINI CG model. The key improvement lies in implementation of additional local structural types and parameterization of the model against a structural library of TM helices obtained from the Protein Data Bank. Through simulations of 12 membrane proteins, we have demonstrated that the improved model not only accurately describes the local hinge structures of these proteins but also reproduces the overall structures of TM domains better than the original MARTINI is able to do. Furthermore, we show that the improved model, when combined with an elastic network, can now be used to explore the deactivation of the human β2 adrenergic receptor and the gating motions of the KcsA potassium channel. The important details of these hinge-associated transitions captured in our simulations agree well with previous experiments and all-atom simulations, and it was concluded that the improved model shows promise as a useful tool for the study of functional dynamics of membrane proteins.

[1]  Roland Faller,et al.  Coarse-Grained Modeling of Polystyrene in Various Environments by Iterative Boltzmann Inversion , 2012 .

[2]  M. Bansal,et al.  HELANAL: A Program to Characterize Helix Geometry in Proteins , 2000, Journal of biomolecular structure & dynamics.

[3]  Charles R. Sanders,et al.  The Amyloid Precursor Protein Has a Flexible Transmembrane Domain and Binds Cholesterol , 2012, Science.

[4]  Helgi I. Ingólfsson,et al.  Martini straight: Boosting performance using a shorter cutoff and GPUs , 2016, Comput. Phys. Commun..

[5]  Durba Sengupta,et al.  Identification of cholesterol binding sites in the serotonin1A receptor. , 2012, The journal of physical chemistry. B.

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

[7]  Hyeon Joo,et al.  OPM database and PPM web server: resources for positioning of proteins in membranes , 2011, Nucleic Acids Res..

[8]  R. Wade,et al.  On the application of the MARTINI coarse-grained model to immersion of a protein in a phospholipid bilayer. , 2015, The Journal of chemical physics.

[9]  Marc Baaden,et al.  Coarse-grain modelling of protein-protein interactions. , 2013, Current opinion in structural biology.

[10]  Tim Werner,et al.  Kink Characterization and Modeling in Transmembrane Protein Structures , 2013, J. Chem. Inf. Model..

[11]  James U Bowie,et al.  Shifting hydrogen bonds may produce flexible transmembrane helices , 2012, Proceedings of the National Academy of Sciences.

[12]  Youxing Jiang,et al.  Crystal structure and mechanism of a calcium-gated potassium channel , 2002, Nature.

[13]  S. Rasmussen,et al.  Structure of a nanobody-stabilized active state of the β2 adrenoceptor , 2010, Nature.

[14]  P. S. Kim,et al.  X-ray structure of the GCN4 leucine zipper, a two-stranded, parallel coiled coil. , 1991, Science.

[15]  Albert C. Pan,et al.  Activation mechanism of the β2-adrenergic receptor , 2011, Proceedings of the National Academy of Sciences.

[16]  Cinque S. Soto,et al.  Structural and dynamic mechanisms for the function and inhibition of the M2 proton channel from influenza A virus. , 2011, Current opinion in structural biology.

[17]  Benoît Roux,et al.  Ion conduction and selectivity in K(+) channels. , 2005, Annual review of biophysics and biomolecular structure.

[18]  Sebastian Kelm,et al.  Examining the Conservation of Kinks in Alpha Helices , 2016, PloS one.

[19]  H. Weinstein,et al.  Structural mimicry of proline kinks: tertiary packing interactions support local structural distortions. , 2002, Journal of molecular biology.

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

[21]  Eduardo Perozo,et al.  Structural mechanism of C-type inactivation in K+ channels , 2010, Nature.

[22]  Helgi I. Ingólfsson,et al.  Computational Lipidomics with insane: A Versatile Tool for Generating Custom Membranes for Molecular Simulations. , 2015, Journal of chemical theory and computation.

[23]  M. Sansom,et al.  Organization and Dynamics of Receptor Proteins in a Plasma Membrane. , 2015, Journal of the American Chemical Society.

[24]  Nagarajan Vaidehi,et al.  Position of helical kinks in membrane protein crystal structures and the accuracy of computational prediction. , 2009, Journal of molecular graphics & modelling.

[25]  A. Laio,et al.  Escaping free-energy minima , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[26]  D. Tieleman,et al.  Perspective on the Martini model. , 2013, Chemical Society reviews.

[27]  Adolfo B Poma,et al.  Combining the MARTINI and Structure-Based Coarse-Grained Approaches for the Molecular Dynamics Studies of Conformational Transitions in Proteins. , 2017, Journal of chemical theory and computation.

[28]  D. Tieleman,et al.  Improving Internal Peptide Dynamics in the Coarse-Grained MARTINI Model: Toward Large-Scale Simulations of Amyloid- and Elastin-like Peptides , 2012, Journal of chemical theory and computation.

[29]  Hengameh Shams,et al.  α-Actinin Induces a Kink in the Transmembrane Domain of β3-Integrin and Impairs Activation via Talin. , 2017, Biophysical journal.

[30]  J. McCoy,et al.  Molecular interactions involved in proton-dependent gating in KcsA potassium channels , 2013, The Journal of general physiology.

[31]  A. Kolinski,et al.  Coarse-Grained Protein Models and Their Applications. , 2016, Chemical reviews.

[32]  A. Thompson,et al.  Molecular mechanism of pH sensing in KcsA potassium channels , 2008, Proceedings of the National Academy of Sciences.

[33]  James U Bowie,et al.  Membrane Protein Twists and Turns , 2013, Science.

[34]  Lea Thøgersen,et al.  Conformational flexibility of the leucine binding protein examined by protein domain coarse-grained molecular dynamics , 2013, Journal of Molecular Modeling.

[35]  Gwyndaf Evans,et al.  Membrane protein structure determination — The next generation , 2014, Biochimica et biophysica acta.

[36]  Thomas Huber,et al.  G protein-coupled receptors self-assemble in dynamics simulations of model bilayers. , 2007, Journal of the American Chemical Society.

[37]  W F Drew Bennett,et al.  Improved Parameters for the Martini Coarse-Grained Protein Force Field. , 2013, Journal of chemical theory and computation.

[38]  W. Han,et al.  Initial Substrate Binding of γ-Secretase: The Role of Substrate Flexibility. , 2017, ACS chemical neuroscience.

[39]  W. Catterall,et al.  Overview of Molecular Relationships in the Voltage-Gated Ion Channel Superfamily , 2005, Pharmacological Reviews.

[40]  B. Chait,et al.  The structure of the potassium channel: molecular basis of K+ conduction and selectivity. , 1998, Science.

[41]  Duan Yang,et al.  The evolution of transmembrane helix kinks and the structural diversity of G protein-coupled receptors. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[42]  S. White,et al.  Membrane protein folding and stability: physical principles. , 1999, Annual review of biophysics and biomolecular structure.

[43]  S. Rasmussen,et al.  Crystal Structure of the β2Adrenergic Receptor-Gs protein complex , 2011, Nature.

[44]  R. Stevens,et al.  GPCR Engineering Yields High-Resolution Structural Insights into β2-Adrenergic Receptor Function , 2007, Science.

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

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

[47]  Alexander D. MacKerell,et al.  CHARMM-GUI Input Generator for NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM Simulations Using the CHARMM36 Additive Force Field , 2015, Journal of chemical theory and computation.

[48]  David S. Goodsell,et al.  The RCSB protein data bank: integrative view of protein, gene and 3D structural information , 2016, Nucleic Acids Res..

[49]  C.-M. Chen,et al.  Statistical analyses and computational prediction of helical kinks in membrane proteins , 2012, Journal of Computer-Aided Molecular Design.

[50]  Berk Hess,et al.  GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers , 2015 .

[51]  F. Cordes,et al.  Conformational dynamics of helix S6 from Shaker potassium channel: simulation studies. , 2002, Biopolymers.

[52]  A. Labro,et al.  Being Flexible: The Voltage-Controllable Activation Gate of Kv Channels , 2012, Front. Pharmacol..

[53]  Rafael C. Bernardi,et al.  Molecular dynamics simulations of large macromolecular complexes. , 2015, Current opinion in structural biology.

[54]  Christian Blouin,et al.  Improved Helix and Kink Characterization in Membrane Proteins Allows Evaluation of Kink Sequence Predictors , 2010, J. Chem. Inf. Model..

[55]  Jesús Giraldo,et al.  Helix 3 acts as a conformational hinge in Class A GPCR activation: An analysis of interhelical interaction energies in crystal structures. , 2015, Journal of structural biology.

[56]  Helgi I Ingólfsson,et al.  CHARMM-GUI Martini Maker for Coarse-Grained Simulations with the Martini Force Field. , 2015, Journal of chemical theory and computation.

[57]  R. Stevens,et al.  Structure-function of the G protein-coupled receptor superfamily. , 2013, Annual review of pharmacology and toxicology.

[58]  J. Doye,et al.  Global Optimization by Basin-Hopping and the Lowest Energy Structures of Lennard-Jones Clusters Containing up to 110 Atoms , 1997, cond-mat/9803344.

[59]  Mark S. P. Sansom,et al.  The pore of voltage-gated potassium ion channels is strained when closed , 2013, Nature Communications.

[60]  Xavier Periole,et al.  Combining an Elastic Network With a Coarse-Grained Molecular Force Field: Structure, Dynamics, and Intermolecular Recognition. , 2009, Journal of chemical theory and computation.

[61]  C. Deane,et al.  Helix kinks are equally prevalent in soluble and membrane proteins , 2014, Proteins.

[62]  E. Campbell,et al.  Crystal Structure of a Mammalian Voltage-Dependent Shaker Family K+ Channel , 2005, Science.

[63]  Helgi I. Ingólfsson,et al.  Lipid–Protein Interactions Are Unique Fingerprints for Membrane Proteins , 2017, bioRxiv.

[64]  R. Swendsen,et al.  THE weighted histogram analysis method for free‐energy calculations on biomolecules. I. The method , 1992 .

[65]  Clare McCabe,et al.  Derivation of coarse-grained potentials via multistate iterative Boltzmann inversion. , 2014, The Journal of chemical physics.

[66]  Naomi R. Latorraca,et al.  GPCR Dynamics: Structures in Motion. , 2017, Chemical reviews.

[67]  Mark Johnson,et al.  NCBI BLAST: a better web interface , 2008, Nucleic Acids Res..

[68]  Xavier Periole,et al.  Interplay of G Protein-Coupled Receptors with the Membrane: Insights from Supra-Atomic Coarse Grain Molecular Dynamics Simulations. , 2017, Chemical reviews.

[69]  Dirk Reith,et al.  Deriving effective mesoscale potentials from atomistic simulations , 2002, J. Comput. Chem..

[70]  Youxing Jiang,et al.  High Resolution Structure of the Open NaK Channel , 2008, Nature Structural &Molecular Biology.

[71]  N. Vaidehi,et al.  Differences in allosteric communication pipelines in the inactive and active states of a GPCR. , 2014, Biophysical journal.

[72]  T. Hori,et al.  Internally bridging water molecule in transmembrane alpha-helical kink. , 2010, Current opinion in structural biology.

[73]  Jean-Pierre Mazat,et al.  Evidence for cardiolipin binding sites on the membrane-exposed surface of the cytochrome bc1. , 2013, Journal of the American Chemical Society.

[74]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

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

[76]  E. Perozo,et al.  A molecular mechanism for proton‐dependent gating in KcsA , 2010, FEBS letters.

[77]  Charles R. Sanders,et al.  Backbone Hydrogen Bond Strengths Can Vary Widely in Transmembrane Helices , 2017, Journal of the American Chemical Society.

[78]  Albert C. Pan,et al.  The Dynamic Process of β2-Adrenergic Receptor Activation , 2013, Cell.

[79]  Min Goo Lee,et al.  A protein sequence that can encode native structure by disfavoring alternate conformations , 2002, Nature Structural Biology.