Inspired by thermodynamic integration, we propose a method for the calculation of time-independent free energy profiles from history-dependent biased simulations via Mean Force Integration (MFI). MFI circumvents the need for computing the ensemble average of the bias acting on the system c(t) and can be applied to different variants of metadynamics. Moreover, MFI naturally extends to aggregate information obtained from independent metadynamics simulations, allowing to converge free energy surfaces without the need to sample recrossing events in a single continuous trajectory. We validate MFI against one- and two-dimensional analytical potentials and by computing the conformational free energy landscape of ibuprofen in the bulk of its most common crystal phase.
REFERENCES
1.
C.
Chipot
and A.
Pohorille
, Free Energy Calculations
(Springer
, 2007
).2.
G. M.
Torrie
and J. P.
Valleau
, “Nonphysical sampling distributions in Monte Carlo free-energy estimation: Umbrella sampling
,” J. Comput. Phys.
23
(2
), 187
–199
(1977
).3.
S.
Marsili
, A.
Barducci
, R.
Chelli
, P.
Procacci
, and V.
Schettino
, “Self-healing umbrella sampling: A non-equilibrium approach for quantitative free energy calculations
,” J. Phys. Chem. B
110
(29
), 14011
–14013
(2006
).4.
H.
Grubmüller
, “Predicting slow structural transitions in macromolecular systems: Conformational flooding
,” Phys. Rev. E
52
(3
), 2893
(1995
).5.
A. F.
Voter
, “Hyperdynamics: Accelerated molecular dynamics of infrequent events
,” Phys. Rev. Lett.
78
(20
), 3908
(1997
).6.
E.
Darve
and A.
Pohorille
, “Calculating free energies using average force
,” J. Chem. Phys.
115
(20
), 9169
–9183
(2001
).7.
A.
Laio
and M.
Parrinello
, “Escaping free-energy minima
,” Proc. Natl. Acad. Sci. U. S. A.
99
(20
), 12562
–12566
(2002
).8.
A.
Barducci
, G.
Bussi
, and M.
Parrinello
, “Well-tempered metadynamics: A smoothly converging and tunable free-energy method
,” Phys. Rev. Lett.
100
(2
), 020603
(2008
).9.
A.
Barducci
, M.
Bonomi
, and M.
Parrinello
, “Metadynamics
,” Wiley Interdiscip. Rev.: Comput. Mol. Sci.
1
(5
), 826
–843
(2011
).10.
O.
Valsson
, P.
Tiwary
, and M.
Parrinello
, “Enhancing important fluctuations: Rare events and metadynamics from a conceptual viewpoint
,” Annu. Rev. Phys. Chem.
67
, 159
–184
(2016
).11.
A.
Laio
and F. L.
Gervasio
, “Metadynamics: A method to simulate rare events and reconstruct the free energy in biophysics, chemistry and material science
,” Rep. Prog. Phys.
71
(12
), 126601
(2008
).12.
P.
Tiwary
and M.
Parrinello
, “A time-independent free energy estimator for metadynamics
,” J. Phys. Chem. B
119
(3
), 736
–742
(2015
).13.
M.
Bonomi
, A.
Barducci
, and M.
Parrinello
, “Reconstructing the equilibrium Boltzmann distribution from well-tempered metadynamics
,” J. Comput. Chem.
30
(11
), 1615
–1621
(2009
).14.
R. W.
Zwanzig
, “High-temperature equation of state by a perturbation method. I. Nonpolar gases
,” J. Chem. Phys.
22
(8
), 1420
–1426
(1954
).15.
H.
Christoper
, J. C.
Sharp
, and F. G.
Prendergast
, “Tryptophan-47 rotational isomerization in variant-3 scorpion neurotoxin. A combination thermodynamic perturbation and umbrella sampling study
,” Biophys. J.
57
(6
), 1269
–1279
(1990
).16.
B.
Roux
, “The calculation of the potential of mean force using computer simulations
,” Comput. Phys. Commun.
91
(1-3
), 275
–282
(1995
).17.
S.
Kumar
, J. M.
Rosenberg
, D.
Bouzida
, R. H.
Swendsen
, and P. A.
Kollman
, “The weighted histogram analysis method for free-energy calculations on biomolecules. I. The method
,” J. Comput. Chem.
13
(8
), 1011
–1021
(1992
).18.
J.
Kästner
and W.
Thiel
, “Bridging the gap between thermodynamic integration and umbrella sampling provides a novel analysis method: “Umbrella integration”
,” J. Chem. Phys.
123
(14
), 144104
(2005
).19.
J. F.
Dama
, G.
Rotskoff
, M.
Parrinello
, and G. A.
Voth
, “Transition-tempered metadynamics: Robust, convergent metadynamics via on-the-fly transition barrier estimation
,” J. Chem. Theory Comput.
10
(9
), 3626
–3633
(2014
).20.
D.
Branduardi
, G.
Bussi
, and M.
Parrinello
, “Metadynamics with adaptive Gaussians
,” J. Chem. Theory Comput.
8
(7
), 2247
–2254
(2012
).21.
J.
Kästner
, “Umbrella integration in two or more reaction coordinates
,” J. Chem. Phys.
131
(3
), 034109
(2009
).22.
I.
Gimondi
, G. A.
Tribello
, and M.
Salvalaglio
, “Building maps in collective variable space
,” J. Chem. Phys.
149
(10
), 104104
(2018
).23.
G. A.
Tribello
, M.
Bonomi
, D.
Branduardi
, C.
Camilloni
, and G.
Bussi
, “Plumed 2: New feathers for an old bird
,” Comput. Phys. Commun.
185
(2
), 604
–613
(2014
).24.
D.
Van Der Spoel
, E.
Lindahl
, B.
Hess
, G.
Groenhof
, A. E.
Mark
, and H. J. C.
Berendsen
, “GROMACS: Fast, flexible, and free
,” J. Comput. Chem.
26
(16
), 1701
–1718
(2005
).25.
J.
Wang
, R. M.
Wolf
, J. W.
Caldwell
, P. A.
Kollman
, and D. A.
Case
, “Development and testing of a general amber force field
,” J. Comput. Chem.
25
(9
), 1157
–1174
(2004
).26.
V.
Marinova
, G. P. F.
Wood
, I.
Marziano
, and M.
Salvalaglio
, “Dynamics and thermodynamics of ibuprofen conformational isomerism at the crystal/solution interface
,” J. Chem. Theory Comput.
14
(12
), 6484
–6494
(2018
).27.
The PLUMED Consortium
, “Promoting transparency and reproducibility in enhanced molecular simulations
,” Nat. Methods
16
(8
), 670
–673
(2019
).28.
M.
Salvalaglio
, P.
Tiwary
, and M.
Parrinello
, “Assessing the reliability of the dynamics reconstructed from metadynamics
,” J. Chem. Theory Comput.
10
(4
), 1420
–1425
(2014
).29.
P.
Tiwary
and M.
Parrinello
, “From metadynamics to dynamics
,” Phys. Rev. Lett.
111
(23
), 230602
(2013
).© 2019 Author(s).
2019
Author(s)
You do not currently have access to this content.
