Abstract:In 2012, Kim and Hirata derived two generalized Langevin equations (GLEs) for a biomolecule in water, one for the structural fluctuation of the biomolecule and the other for the density fluctuation of water, by projecting all the mechanical variables in phase space onto the two dynamic variables: the structural fluctuation defined by the displacement of atoms from their equilibrium positions, and the solvent density fluctuation. The equation has an expression similar to the classical Langevin equation (CLE) for a harmonic oscillator, possessing terms corresponding to the restoring force proportional to the structural fluctuation, as well as the frictional and random forces. However, there is a distinct difference between the two expressions that touches on the essential physics of the structural fluctuation, that is, the force constant, or Hessian, in the restoring force. In the CLE, this is given by the second derivative of the potential energy among atoms in a protein. So, the quadratic nature or the harmonicity is only valid at the minimum of the potential surface. On the contrary, the linearity of the restoring force in the GLE originates from the projection of the water’s degrees of freedom onto the protein’s degrees of freedom. Taking this into consideration, Kim and Hirata proposed an ansatz for the Hessian matrix. The ansatz is used to equate the Hessian matrix with the second derivative of the free-energy surface or the potential of the mean force of a protein in water, defined by the sum of the potential energy among atoms in a protein and the solvation free energy. Since the free energy can be calculated from the molecular mechanics and the RISM/3D-RISM theory, one can perform an analysis similar to the normal mode analysis (NMA) just by diagonalizing the Hessian matrix of the free energy. This method is referred to as the Generalized Langevin Mode Analysis (GLMA). This theory may be realized to explore a variety of biophysical processes, including protein folding, spectroscopy, and chemical reactions. The present article is devoted to reviewing the development of this theory, and to providing perspective in exploring life phenomena.

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

Protein Folding: A Perspective From Statistical Physics

Molecular Dynamics Simulation of Protein Conformation Transition in Solution

Cooperativity, Local-Nonlocal Coupling, and Nonnative Interactions: Principles of Protein Folding from Coarse-Grained Models

Protein Folding Kinetics: Time Scales, Pathways, and Energy Landscapes in Terms of Sequence Dependent Properties

Matching Simulation and Experiment: A New Simplified Model for Simulating Protein Folding

Modeling Coil–Globule–Helix Transition in Polymers by Self-Interacting Random Walks

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

Global Langevin model of multidimensional biomolecular dynamics

Statistical Mechanics Model for Protein Folding

Experimental and Theoretical Studies of Protein Folding‐Unfolding

Folding Simulations for A Three-Helix Bundle Protein

Dominant Pathways in Protein Folding

The effect of local thermal fluctuations on the folding kinetics: a study from the perspective of the nonextensive statistical mechanics

Structural Fluctuation, Relaxation, and Folding of Protein: An Approach Based on the Combined Generalized Langevin and RISM/3D-RISM Theories

Thermodynamics of Force-Dependent Folding and Unfolding of Small Protein and Nucleic Acid Structures.

Kinetic Analysis of Protein Folding Lattice Models

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

Discrete molecular dynamics studies of the folding of a protein-like model

Oscillatory molecular driving force for protein folding at high concentration: a molecular simulation.

Protein Folding Simulations Using the Activation-Relaxation Technique