Abstract:Molecular dynamics simulations of biomolecular processes are often discussed in terms of diffusive motion on a low-dimensional free energy landscape F(𝒙). To provide a theoretical basis for this interpretation, one may invoke the system-bath ansatz á la Zwanzig. That is, by assuming a time scale separation between the slow motion along the system coordinate x and the fast fluctuations of the bath, a memory-free Langevin equation can be derived that describes the system's motion on the free energy landscape F(𝒙), which is damped by a friction field and driven by a stochastic force that is related to the friction via the fluctuation-dissipation theorem. While the theoretical formulation of Zwanzig typically assumes a highly idealized form of the bath Hamiltonian and the system-bath coupling, one would like to extend the approach to realistic data-based biomolecular systems. Here a practical method is proposed to construct an analytically defined global model of structural dynamics. Given a molecular dynamics simulation and adequate collective coordinates, the approach employs an "empirical valence bond"-type model which is suitable to represent multidimensional free energy landscapes as well as an approximate description of the friction field. Adopting alanine dipeptide and a three-dimensional model of heptaalanine as simple examples, the resulting Langevin model is shown to reproduce the results of the underlying all-atom simulations. Because the Langevin equation can also be shown to satisfy the underlying assumptions of the theory (such as a delta-correlated Gaussian-distributed noise), the global model provides a correct, albeit empirical, realization of Zwanzig's formulation. As an application, the model can be used to investigate the dependence of the system on parameter changes and to predict the effect of site-selective mutations on the dynamics.

Protein conformational change based on a two-dimensional generalized langevin equation

Conformation Studies of Two-Dimensional Model Molecules of Proteins in the Process of Folding

On the Stochastic Dynamics of Molecular Conformation

Langevin Equation and Its Applications in the Simulation of Protein Folding Process

Generalized Langevin Equation With Fractional Gaussian Noise: Subdiffusion Within A Single Protein Molecule

Multiple Time Scale Dynamics of Distance Fluctuations in a Semiflexible Polymer: a One-Dimensional Generalized Langevin Equation Treatment.

Distance Fluctuation of a Single Molecule in Lennard-Jones Liquid Based on Generalized Langevin Equation and Mode Coupling Theory

Structural Fluctuation of Protein in Water around Its Native State: A New Statistical Mechanics Formulation

Equations for Stochastic Macromolecular Mechanics of Single Proteins: Equilibrium Fluctuations, Transient Kinetics and Nonequilibrium Steady-State

Non-Markovian effects of conformational fluctuations on the global diffusivity in Langevin equation with fluctuating diffusivity

A quantum mechanical computational method for modeling electrostatic and solvation effects of protein

Global Langevin model of multidimensional biomolecular dynamics

A generalized Langevin equation approach for barrier crossing dynamics in conformational transitions of proteins

Molecular Dynamics Simulation of Protein Conformation Transition in Solution

Dimension Conversion and Scaling of Disordered Protein Chains

Biomolecular conformational changes and ligand binding: from kinetics to thermodynamics.

On the derivation of the generalized Langevin equation and the fluctuation-dissipation theorem

A simple way to compute protein dynamics without a mechanical model.

Langevin Dynamics Of Correlated Subdiffusion And Normal Diffusion

Experimental and Theoretical Studies of Protein Folding‐Unfolding

Langevin analogy between particle trajectories and polymer configurations