Two Dimensional Window Exchange Umbrella Sampling for Transmembrane Helix Assembly.

The method of window exchange umbrella sampling molecular dynamics (WEUSMD) with a pre-optimized parameter set was recently used to obtain the most probable conformations and the energetics of transmembrane (TM) helix assembly of a generic TM sequence. When applied to glycophorin A TM domain (GpA-TM) using the restraint potentials along the helix-helix distance, however, tight interfacial packing of GpA-TM resulted in insufficient conformational sampling at short helix-helix separation. This sampling issue is addressed by extending the WEUSMD into two dimensions with the restraint potentials along the helix-helix distance and crossing angle. The two-dimensional WEUSMD results demonstrate that the incomplete sampling in the one-dimensional WEUSMD arises from high barriers along the crossing angle between the GpA-TM helices. Together with the faster convergence in both the assembled conformations and the potential of mean force, the 2D-WEUSMD can be a general and efficient approach in computational studies of TM helix assembly.

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

[2]  J. Bowie,et al.  Dramatic destabilization of transmembrane helix interactions by features of natural membrane environments. , 2011, Journal of the American Chemical Society.

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

[4]  Anton A Polyansky,et al.  Structural, dynamic, and functional aspects of helix association in membranes: a computational view. , 2011, Advances in protein chemistry and structural biology.

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

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

[7]  Sunhwan Jo,et al.  Cholesterol flip-flop: insights from free energy simulation studies. , 2010, The journal of physical chemistry. B.

[8]  T. Lazaridis,et al.  Calculating the free energy of association of transmembrane helices. , 2006, Biophysical journal.

[9]  Christopher G. Tate,et al.  Overcoming barriers to membrane protein structure determination , 2011, Nature Biotechnology.

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

[11]  Eric Darve,et al.  Adaptive biasing force method for scalar and vector free energy calculations. , 2008, The Journal of chemical physics.

[12]  Wang,et al.  Replica Monte Carlo simulation of spin glasses. , 1986, Physical review letters.

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

[14]  Ulrich H E Hansmann,et al.  Generalized ensemble and tempering simulations: a unified view. , 2006, Physical review. E, Statistical, nonlinear, and soft matter physics.

[15]  Robert Fredriksson,et al.  Mapping the human membrane proteome : a majority of the human membrane proteins can be classified according to function and evolutionary origin , 2015 .

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

[17]  Wei Yang,et al.  Random walk in orthogonal space to achieve efficient free-energy simulation of complex systems , 2008, Proceedings of the National Academy of Sciences.

[18]  A. Laio,et al.  Metadynamics: a method to simulate rare events and reconstruct the free energy in biophysics, chemistry and material science , 2008 .

[19]  Gerhard König,et al.  Non‐Boltzmann sampling and Bennett's acceptance ratio method: How to profit from bending the rules , 2011, J. Comput. Chem..

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

[21]  P. Tavan,et al.  Efficiency of exchange schemes in replica exchange , 2009 .

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

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

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

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

[26]  C. Chothia,et al.  Helix to helix packing in proteins. , 1981, Journal of molecular biology.

[27]  A. Pohorille,et al.  Insights into the recognition and association of transmembrane alpha-helices. The free energy of alpha-helix dimerization in glycophorin A. , 2005, Journal of the American Chemical Society.