Path corrected functionals of stochastic trajectories: towards relative free energy and reaction coordinate calculations

Abstract A numerical method for estimating the functional of a stochastic differential equation is evaluated for possible application to biomolecular systems. The approach uses variance reduction for stochastic trajectories to estimate relative free energy differences and define paths between initial and end points. The method shows promise as evidenced by application to simple one and two dimensional trial surfaces. This suggests that the approach may be applicable to more complex biomolecular systems.

[1]  Shankar Kumar,et al.  Multidimensional free‐energy calculations using the weighted histogram analysis method , 1995, J. Comput. Chem..

[2]  Peter A. Kollman,et al.  FREE ENERGY CALCULATIONS : APPLICATIONS TO CHEMICAL AND BIOCHEMICAL PHENOMENA , 1993 .

[3]  Ron Elber,et al.  Calculation of classical trajectories with a very large time step: Formalism and numerical examples , 1996 .

[4]  R. Elber,et al.  Reaction path study of conformational transitions in flexible systems: Applications to peptides , 1990 .

[5]  D. L. Freeman,et al.  The construction of double-ended classical trajectories , 1994 .

[6]  K. Müller,et al.  Reaction Paths on Multidimensional Energy Hypersurfaces , 1980 .

[7]  B. Roux.  The calculation of the potential of mean force using computer simulations , 1995 .

[8]  Hans Christian Öttinger,et al.  Stochastic Processes in Polymeric Fluids , 1996 .

[9]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .

[10]  W. Wagner,et al.  Undiased monte carlo estimators for functionals of weak solutions of stochastic diffretial equations , 1989 .

[11]  R. Berry,et al.  Finding saddles on multidimensional potential surfaces , 1988 .

[12]  Lawrence R. Pratt,et al.  A statistical method for identifying transition states in high dimensional problems , 1986 .

[13]  H. C. Öttinger.  Variance Reduced Brownian Dynamics Simulations , 1994 .

[14]  J. Straub,et al.  The MaxFlux algorithm for calculating variationally optimized reaction paths for conformational transitions in many body systems at finite temperature , 1997 .

[15]  W. Wagner.  Unbiased Monte Carlo evaluation of certain functional integrals , 1987 .

[16]  W. Wagner.  Monte carlo evaluation of functionals of solutions of stochastic differential equations. variance reduction and numerical examples , 1988 .

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

[18]  C. Gardiner.  Handbook of Stochastic Methods , 1983 .

[19]  Kent R. Wilson,et al.  Shadowing, rare events, and rubber bands. A variational Verlet algorithm for molecular dynamics , 1992 .