Time-resolved force distribution analysis

BackgroundBiomolecules or other complex macromolecules undergo conformational transitions upon exposure to an external perturbation such as ligand binding or mechanical force. To follow fluctuations in pairwise forces between atoms or residues during such conformational changes as observed in Molecular Dynamics (MD) simulations, we developed Time-Resolved Force Distribution Analysis (TRFDA).ResultsThe implementation focuses on computational efficiency and low-memory usage and, along with the wide range of output options, makes possible time series analysis of pairwise forces variation in long MD simulations and for large molecular systems. It also provides an exact decomposition of pairwise forces resulting from 3- and 4-body potentials and a unified treatment of pairwise forces between atoms or residues. As a proof of concept, we present a stress analysis during unfolding of ubiquitin in a force-clamp MD simulation.ConclusionsTRFDA can be used, among others, in tracking signal propagation at atomic level, for characterizing dynamical intermolecular interactions (e.g. protein-ligand during flexible docking), in development of force fields and for following stress distribution during conformational changes.

[1]  T. Straatsma,et al.  THE MISSING TERM IN EFFECTIVE PAIR POTENTIALS , 1987 .

[2]  Fei Xia,et al.  Dynamic Allostery in the Methionine Repressor Revealed by Force Distribution Analysis , 2009, PLoS Comput. Biol..

[3]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[4]  Carsten Kutzner,et al.  GROMACS 4:  Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. , 2008, Journal of chemical theory and computation.

[5]  Aidan P Thompson,et al.  General formulation of pressure and stress tensor for arbitrary many-body interaction potentials under periodic boundary conditions. , 2009, The Journal of chemical physics.

[6]  Berk Hess,et al.  P-LINCS:  A Parallel Linear Constraint Solver for Molecular Simulation. , 2008, Journal of chemical theory and computation.

[7]  Frauke Gräter,et al.  Implementation of force distribution analysis for molecular dynamics simulations , 2011, BMC Bioinformatics.

[8]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[9]  M. Parrinello,et al.  Polymorphic transitions in single crystals: A new molecular dynamics method , 1981 .

[10]  M. Parrinello,et al.  Canonical sampling through velocity rescaling. , 2007, The Journal of chemical physics.

[11]  C. Bugg,et al.  Structure of ubiquitin refined at 1.8 A resolution. , 1987, Journal of molecular biology.

[12]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[13]  Frauke Gräter,et al.  Mechanical Network in Titin Immunoglobulin from Force Distribution Analysis , 2009, PLoS Comput. Biol..

[14]  Berk Hess,et al.  LINCS: A linear constraint solver for molecular simulations , 1997, J. Comput. Chem..

[15]  Julio M Fernandez,et al.  Force-Clamp Spectroscopy Monitors the Folding Trajectory of a Single Protein , 2004, Science.

[16]  Gregory A Voth,et al.  A multiscale coarse-graining method for biomolecular systems. , 2005, The journal of physical chemistry. B.

[17]  W. L. Jorgensen,et al.  Development and Testing of the OPLS All-Atom Force Field on Conformational Energetics and Properties of Organic Liquids , 1996 .

[18]  Frauke Gräter,et al.  Force distribution reveals signal transduction in E. coli Hsp90. , 2012, Biophysical journal.

[19]  Frauke Gräter,et al.  Dynamic Prestress in a Globular Protein , 2012, PLoS Comput. Biol..

[20]  S. Nosé,et al.  Constant pressure molecular dynamics for molecular systems , 1983 .