Multidimensional umbrella sampling and replica‐exchange molecular dynamics simulations for structure prediction of transmembrane helix dimers

Structural information of a transmembrane (TM) helix dimer is useful in understanding molecular mechanisms of important biological phenomena such as signal transduction across the cell membrane. Here, we describe an umbrella sampling (US) scheme for predicting the structure of a TM helix dimer in implicit membrane using the interhelical crossing angle and the TM–TM relative rotation angles as the reaction coordinates. This scheme conducts an efficient conformational search on TM–TM contact interfaces, and its robustness is tested by predicting the structures of glycophorin A (GpA) and receptor tyrosine kinase EphA1 (EphA1) TM dimers. The nuclear magnetic resonance (NMR) structures of both proteins correspond to the global free‐energy minimum states in their free‐energy landscapes. In addition, using the landscape of GpA as a reference, we also examine the protocols of temperature replica‐exchange molecular dynamics (REMD) simulations for structure prediction of TM helix dimers in implicit membrane. A wide temperature range in REMD simulations, for example, 250–1000 K, is required to efficiently obtain a free‐energy landscape consistent with the US simulations. The interhelical crossing angle and the TM–TM relative rotation angles can be used as reaction coordinates in multidimensional US and be good measures for conformational sampling of REMD simulations. © 2013 Wiley Periodicals, Inc.

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

[2]  Joseph Schlessinger,et al.  Signal transduction by receptors with tyrosine kinase activity , 1990, Cell.

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

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

[5]  Taehoon Kim,et al.  Novel free energy calculations to explore mechanisms and energetics of membrane protein structure and function , 2009, J. Comput. Chem..

[6]  Alexander D. MacKerell,et al.  All-atom empirical potential for molecular modeling and dynamics studies of proteins. , 1998, The journal of physical chemistry. B.

[7]  M. Karplus,et al.  Simulation of activation free energies in molecular systems , 1996 .

[8]  Jianpeng Ma,et al.  CHARMM: The biomolecular simulation program , 2009, J. Comput. Chem..

[9]  Dieter Langosch,et al.  GxxxG motifs within the amyloid precursor protein transmembrane sequence are critical for the etiology of Aβ42 , 2007 .

[10]  Yuko Okamoto,et al.  Prediction of transmembrane helix configurations by replica-exchange simulations , 2003, cond-mat/0309338.

[11]  N. Go.  Theoretical studies of protein folding. , 1983, Annual review of biophysics and bioengineering.

[12]  Y. Sugita,et al.  Multidimensional replica-exchange method for free-energy calculations , 2000, cond-mat/0009120.

[13]  G. Torrie,et al.  Nonphysical sampling distributions in Monte Carlo free-energy estimation: Umbrella sampling , 1977 .

[14]  E. N. Tkach,et al.  Left-handed dimer of EphA2 transmembrane domain: Helix packing diversity among receptor tyrosine kinases. , 2010, Biophysical journal.

[15]  Martin A. Schwartz,et al.  Cell adhesion: integrating cytoskeletal dynamics and cellular tension , 2010, Nature Reviews Molecular Cell Biology.

[16]  Ying Xu,et al.  Energetics and stability of transmembrane helix packing: A replica‐exchange simulation with a knowledge‐based membrane potential , 2006, Proteins.

[17]  D Thirumalai,et al.  Transmembrane structures of amyloid precursor protein dimer predicted by replica-exchange molecular dynamics simulations. , 2009, Journal of the American Chemical Society.

[18]  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.

[19]  Chungho Kim,et al.  The structure of the integrin αIIbβ3 transmembrane complex explains integrin transmembrane signalling , 2009, The EMBO journal.

[20]  J. Onuchic,et al.  Theory of protein folding: the energy landscape perspective. , 1997, Annual review of physical chemistry.

[21]  Andrei L. Lomize,et al.  OPM: Orientations of Proteins in Membranes database , 2006, Bioinform..

[22]  Y. Sugita,et al.  Replica-exchange molecular dynamics method for protein folding , 1999 .

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

[24]  Wonpil Im,et al.  Two Dimensional Window Exchange Umbrella Sampling for Transmembrane Helix Assembly. , 2013, Journal of chemical theory and computation.

[25]  Michael Feig,et al.  MMTSB Tool Set: enhanced sampling and multiscale modeling methods for applications in structural biology. , 2004, Journal of molecular graphics & modelling.

[26]  M Scott Shell,et al.  Two-dimensional replica exchange approach for peptide-peptide interactions. , 2011, The Journal of chemical physics.

[27]  Klaus Schulten,et al.  A glycophorin A-like framework for the dimerization of photosynthetic core complexes. , 2009, Journal of the American Chemical Society.

[28]  M. V. Goncharuk,et al.  Dimeric Structure of the Transmembrane Domain of Glycophorin A in Lipidic and Detergent Environments , 2010, Acta naturae.

[29]  Sarel J Fleishman,et al.  A novel scoring function for predicting the conformations of tightly packed pairs of transmembrane alpha-helices. , 2002, Journal of molecular biology.

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

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

[32]  Gevorg Grigoryan,et al.  Transmembrane communication: general principles and lessons from the structure and function of the M2 proton channel, K⁺ channels, and integrin receptors. , 2011, Annual review of biochemistry.

[33]  Roman G. Efremov,et al.  Spatial Structure and pH-dependent Conformational Diversity of Dimeric Transmembrane Domain of the Receptor Tyrosine Kinase EphA1* , 2008, Journal of Biological Chemistry.

[34]  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.

[35]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[36]  Anton Arkhipov,et al.  Architecture and Membrane Interactions of the EGF Receptor , 2013, Cell.

[37]  S. Takada,et al.  On the Hamiltonian replica exchange method for efficient sampling of biomolecular systems: Application to protein structure prediction , 2002 .

[38]  A. Laio,et al.  Free-energy landscape for beta hairpin folding from combined parallel tempering and metadynamics. , 2006, Journal of the American Chemical Society.

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

[40]  Alexander D. MacKerell,et al.  Extending the treatment of backbone energetics in protein force fields: Limitations of gas‐phase quantum mechanics in reproducing protein conformational distributions in molecular dynamics simulations , 2004, J. Comput. Chem..

[41]  Jun Qin,et al.  Structure of an integrin αIIbβ3 transmembrane-cytoplasmic heterocomplex provides insight into integrin activation , 2009, Proceedings of the National Academy of Sciences.

[42]  Wonpil Im,et al.  Implementation and application of helix–helix distance and crossing angle restraint potentials , 2007, J. Comput. Chem..

[43]  S. Chandrasekhar.  Stochastic problems in Physics and Astronomy , 1943 .

[44]  Zsuzsanna Dosztányi,et al.  Transmembrane proteins in the Protein Data Bank: identification and classification , 2004, Bioinform..

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

[46]  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.

[47]  Wonpil Im,et al.  Transmembrane helix assembly by window exchange umbrella sampling. , 2012, Physical review letters.

[48]  Themis Lazaridis,et al.  Membrane protein native state discrimination by implicit membrane models , 2013, J. Comput. Chem..

[49]  K. Garcia,et al.  Structural basis of T cell recognition. , 1999, Annual review of immunology.

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

[51]  S. Rasmussen,et al.  The structure and function of G-protein-coupled receptors , 2009, Nature.

[52]  David Baker,et al.  The structure of a receptor with two associating transmembrane domains on the cell surface: integrin alphaIIbbeta3. , 2009, Molecular cell.

[53]  Siewert J Marrink,et al.  Structural determinants of the supramolecular organization of G protein-coupled receptors in bilayers. , 2012, Journal of the American Chemical Society.

[54]  M. Karplus,et al.  Stochastic boundary conditions for molecular dynamics simulations of ST2 water , 1984 .

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