SOP‐GPU: influence of solvent‐induced hydrodynamic interactions on dynamic structural transitions in protein assemblies

Hydrodynamic interactions (HI) are incorporated into Langevin dynamics of the Cα‐based protein model using the Truncated Expansion approximation (TEA) to the Rotne–Prager–Yamakawa diffusion tensor. Computational performance of the obtained GPU realization demonstrates the model's capability for describing protein systems of varying complexity (102–105 residues), including biological particles (filaments, virus shells). Comparison of numerical accuracy of the TEA versus exact description of HI reveals similar results for the kinetics and thermodynamics of protein unfolding. The HI speed up and couple biomolecular transitions through cross‐communication among protein domains, which result in more collective displacements of structure elements governed by more deterministic (less variable) dynamics. The force‐extension/deformation spectra from nanomanipulations in silico exhibit sharper force signals that match well the experimental profiles. Hence, biomolecular simulations without HI overestimate the role of tension/stress fluctuations. Our findings establish the importance of incorporating implicit water‐mediated many‐body effects into theoretical modeling of dynamic processes involving biomolecules. © 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.

[1]  A. Elcock,et al.  Striking Effects of Hydrodynamic Interactions on the Simulated Diffusion and Folding of Proteins. , 2009, Journal of chemical theory and computation.

[2]  A. Elcock,et al.  Absolute protein-protein association rate constants from flexible, coarse-grained Brownian dynamics simulations: the role of intermolecular hydrodynamic interactions in barnase-barstar association. , 2010, Biophysical journal.

[3]  A. Zhmurov,et al.  Generation of random numbers on graphics processors: forced indentation in silico of the bacteriophage HK97. , 2011, The journal of physical chemistry. B.

[4]  Jason C. Crane,et al.  The folding mechanism of a -sheet: the WW domain1 , 2001 .

[5]  Duncan Poole,et al.  Routine Microsecond Molecular Dynamics Simulations with AMBER on GPUs. 2. Explicit Solvent Particle Mesh Ewald. , 2013, Journal of chemical theory and computation.

[6]  F. Young Biochemistry , 1955, The Indian Medical Gazette.

[7]  Hiromi Yamakawa,et al.  Transport Properties of Polymer Chains in Dilute Solution: Hydrodynamic Interaction , 1970 .

[8]  R. Dima,et al.  Molecular investigations into the mechanics of actin in different nucleotide states. , 2011, The journal of physical chemistry. B.

[9]  A. Fersht,et al.  Ultrafast folding of WW domains without structured aromatic clusters in the denatured state , 2001, Proceedings of the National Academy of Sciences of the United States of America.

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

[12]  Marek Cieplak,et al.  Influence of hydrodynamic interactions on mechanical unfolding of proteins , 2007, 0706.1102.

[13]  E. R. Smith,et al.  Hydrodynamic interactions in Brownian dynamics I. Periodic boundary conditions for computer simulations , 1987 .

[14]  Olga Kononova,et al.  Fluctuating Nonlinear Spring Model of Mechanical Deformation of Biological Particles , 2015, PLoS Comput. Biol..

[15]  Tihamér Geyer,et al.  An O(N2) approximation for hydrodynamic interactions in Brownian dynamics simulations. , 2008, The Journal of chemical physics.

[16]  M J Harvey,et al.  ACEMD: Accelerating Biomolecular Dynamics in the Microsecond Time Scale. , 2009, Journal of chemical theory and computation.

[17]  Uwe Winter,et al.  Coarse grained simulations of a small peptide: Effects of finite damping and hydrodynamic interactions , 2009, 0908.1685.

[18]  S Kim,et al.  Hydrodynamic steering effects in protein association. , 1994, Proceedings of the National Academy of Sciences of the United States of America.

[19]  P. Casey,et al.  Crystal Structure of Protein Farnesyltransferase at 2.25 Angstrom Resolution , 1997, Science.

[20]  Klaus Schulten,et al.  Accelerating Molecular Modeling Applications with GPU Computing , 2009 .

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

[22]  R. Ranganathan,et al.  Structural and Functional Analysis of the Mitotic Rotamase Pin1 Suggests Substrate Recognition Is Phosphorylation Dependent , 1997, Cell.

[23]  H. Nguyen,et al.  Tuning the free-energy landscape of a WW domain by temperature, mutation, and truncation , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[24]  M. Cieplak,et al.  Hydrodynamic interactions in protein folding. , 2008, The Journal of chemical physics.

[25]  André E. X. Brown,et al.  Mechanism of fibrin(ogen) forced unfolding. , 2011, Structure.

[26]  D. Thirumalai,et al.  Relative stability of helices determines the folding landscape of adenine riboswitch aptamers. , 2008, Journal of the American Chemical Society.

[27]  O. Dym,et al.  Thermal stabilization of the protozoan Entamoeba histolytica alcohol dehydrogenase by a single proline substitution , 2008 .

[28]  Gerhard Hummer,et al.  Position-dependent diffusion coefficients and free energies from Bayesian analysis of equilibrium and replica molecular dynamics simulations , 2005 .

[29]  J. Kästner Umbrella sampling , 2011 .

[30]  Joshua A. Anderson,et al.  General purpose molecular dynamics simulations fully implemented on graphics processing units , 2008, J. Comput. Phys..

[31]  David P. Anderson,et al.  High-Throughput All-Atom Molecular Dynamics Simulations Using Distributed Computing , 2010, J. Chem. Inf. Model..

[32]  J. Straub,et al.  Propensity to form amyloid fibrils is encoded as excitations in the free energy landscape of monomeric proteins. , 2014, Journal of molecular biology.

[33]  D. Ermak,et al.  Brownian dynamics with hydrodynamic interactions , 1978 .

[34]  Y. Kholodov,et al.  Evaluation of CDK6 and p16/INK4a-Derived Peptides Interaction , 2013 .

[35]  Changbong Hyeon,et al.  Revealing the bifurcation in the unfolding pathways of GFP by using single-molecule experiments and simulations , 2007, Proceedings of the National Academy of Sciences.

[36]  D. Thirumalai,et al.  From mechanical folding trajectories to intrinsic energy landscapes of biopolymers , 2013, Proceedings of the National Academy of Sciences.

[37]  J. Onuchic,et al.  Ligand-induced global transitions in the catalytic domain of protein kinase A , 2009, Proceedings of the National Academy of Sciences.

[38]  K. Marx,et al.  Molecular Mechanisms, Thermodynamics, and Dissociation Kinetics of Knob-Hole Interactions in Fibrin* , 2013, The Journal of Biological Chemistry.

[39]  R. Doolittle,et al.  Crystal structure of human fibrinogen. , 2009, Biochemistry.

[40]  R. Rosenfeld Nature , 2009, Otolaryngology--head and neck surgery : official journal of American Academy of Otolaryngology-Head and Neck Surgery.

[41]  J. Skolnick,et al.  Krylov subspace methods for computing hydrodynamic interactions in brownian dynamics simulations. , 2012, The Journal of chemical physics.

[42]  Olga Kononova,et al.  Structural transitions and energy landscape for Cowpea Chlorotic Mottle Virus capsid mechanics from nanomanipulation in vitro and in silico. , 2013, Biophysical journal.

[43]  Olga Kononova,et al.  Tubulin Bond Energies and Microtubule Biomechanics Determined from Nanoindentation in Silico , 2014, Journal of the American Chemical Society.

[44]  D. Thirumalai,et al.  Kinetics of allosteric transitions in S-adenosylmethionine riboswitch are accurately predicted from the folding landscape. , 2013, Journal of the American Chemical Society.

[45]  C. Tanford Macromolecules , 1994, Nature.

[46]  Changbong Hyeon,et al.  Pathways and kinetic barriers in mechanical unfolding and refolding of RNA and proteins. , 2006, Structure.

[47]  M. Fixman Inclusion of hydrodynamic interaction in polymer dynamical simulations , 1981 .

[48]  V Barsegov,et al.  Order statistics theory of unfolding of multimeric proteins. , 2010, Biophysical journal.

[49]  R. Tehver,et al.  Rigor to post-rigor transition in myosin V: link between the dynamics and the supporting architecture. , 2010, Structure.

[50]  Olga Kononova,et al.  Mechanical transition from α-helical coiled coils to β-sheets in fibrin(ogen). , 2012, Journal of the American Chemical Society.

[51]  L. Christophorou Science , 2018, Emerging Dynamics: Science, Energy, Society and Values.

[52]  V. Barsegov,et al.  Order statistics inference for describing topological coupling and mechanical symmetry breaking in multidomain proteins. , 2015, The Journal of chemical physics.

[53]  R. R. Schmidt,et al.  Comparison of Brownian dynamics algorithms with hydrodynamic interaction. , 2011, The Journal of chemical physics.

[54]  Zydrunas Gimbutas,et al.  A fast multipole method for the Rotne-Prager-Yamakawa tensor and its applications , 2013, J. Comput. Phys..

[55]  V Barsegov,et al.  Sop‐GPU: Accelerating biomolecular simulations in the centisecond timescale using graphics processors , 2010, Proteins.

[56]  S. Prager,et al.  Variational Treatment of Hydrodynamic Interaction in Polymers , 1969 .

[57]  S. Günther,et al.  Structural basis of enzyme encapsulation into a bacterial nanocompartment , 2008, Nature Structural &Molecular Biology.

[58]  D. Thirumalai,et al.  Protein folding kinetics: timescales, pathways and energy landscapes in terms of sequence-dependent properties. , 1996, Folding & design.

[59]  Weitao Yang,et al.  Mechanical anisotropy of ankyrin repeats. , 2012, Biophysical journal.

[60]  D. Klimov,et al.  Analyzing forced unfolding of protein tandems by ordered variates, 2: dependent unfolding times. , 2007, Biophysical journal.

[61]  Duncan Poole,et al.  Routine Microsecond Molecular Dynamics Simulations with AMBER on GPUs. 1. Generalized Born , 2012, Journal of chemical theory and computation.

[62]  E. Nogales,et al.  High-Resolution Model of the Microtubule , 1999, Cell.