Binding Free Energy Calculations and Molecular Dynamics Studies ...
Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier...
Transcript of Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier...
![Page 1: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/1.jpg)
Free Energy Calculation
in MD Simulation
![Page 2: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/2.jpg)
Basic Thermodynamics
● Helmoholtz free energy
A = U – TS + Σ μi Ni
dA = wrev (reversible, const N V T) eq (22.9) McQuarrie & Simon
● Gibbs free energy
G = U + PV – TS + Σ μi Ni
= H – TS + Σ μi Ni
dG = wnonPV
(reversible, const N P T) eq (22.16) McQuarrie & Simon
U internal energy P Pressure μi Chemical potential
V Volume T Temperature
S Entropy dS = dq/T S = kB ln W
![Page 3: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/3.jpg)
Implication of Free Energy
● A ↔ B Keq = [A]/[B]
Keq = exp (-ΔG0/RT)
ΔG0 = -RT ln Keq
ΔG = ΔG0 + RT ln Q
ΔG > 0 Unfavorable
ΔG = 0
ΔG < 0 Favorable
![Page 4: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/4.jpg)
Statistical Mechanics
● System can be described using Hamiltonian
H(p1,p2,.........pN, r1, r2,.......r
N)
● Different ensemble (fixed system quantities)
- Canonical ensemble (N,V,T)
- NPT ensemble
![Page 5: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/5.jpg)
Statistical Mechanics
![Page 6: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/6.jpg)
Common Free Energy Type
● Solvation Free Energy / Transfer Free Energy
● Binding Free energy
● Confomational Free Energy
![Page 7: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/7.jpg)
Calculating Free Energy
● Experimentally
- probabilities of finding the system at given states
ΔG = -RT ln (P1/P
0)
- reversible work of moving the system between two states
● Computationally
- pobability
- reversible work → more efficient
![Page 8: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/8.jpg)
Probability Method ( Brute Force )
● Ergodic hypothesis
time average = ensemble average
● ΔG = -RT ln (P1/P
0)
● Is this true? In general, NO!
- sampling time is too short computationally
- hard to sample every state of the system
For system state of zero sampling, PA = 0 → ΔG
0 → A = ∞
A=limt ∞
1t∫A t dt=
1M∑
MA i=⟨A⟩
![Page 9: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/9.jpg)
Alanine Dipeptide
![Page 10: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/10.jpg)
![Page 11: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/11.jpg)
![Page 12: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/12.jpg)
![Page 13: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/13.jpg)
![Page 14: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/14.jpg)
Challenge
● Accurate calculations of absolute free energy is nearly impossible due to insufficient sampling in a finite length and time scale simulation.
● Need different methods to estimate free energy.● Common method
- Thermodynamical integration
- Free energy perturbation
- Umbrella sampling
- Potential of mean force
![Page 15: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/15.jpg)
Thermodynamic Cycle
![Page 16: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/16.jpg)
H x , px ,=H 0 x , p xH b x , px 1−H a x , p x
![Page 17: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/17.jpg)
Thermodynamic Integration
![Page 18: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/18.jpg)
Thermodynamic Integration
● the value of <dA/dλ> is accurately determined for a number of intermediate values of λ, the total free energy is determined with numerical integration methods based on these values
● A(λ) is smooth enough and converged● Intermediate values:
fluctuation of ∂H/∂λ for each value of dA/dλ
![Page 19: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/19.jpg)
Free Energy Perturbation
![Page 20: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/20.jpg)
![Page 21: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/21.jpg)
Error Estimation
Chris Chipot @ NAMD mailing list
● In FEP, convergence may be probed by monitoring the time-evolution of the ensemble average. This is, however, a necessary, but not sufficient condition for convergence, because apparent plateaus of the ensemble average often conceal anomalous overlap of the density of states characterizing the initial and the final states. The latter should be the key-criterion to ascertain the local convergence of the simulation for those degrees of freedom that are effectively sampled.
● Statistical errors in FEP calculations may be estimated by means of a first-order expansion of the free energy, which involves an estimation of the sampling ratio of the latter of the calculation (Straatsma, 1986).
http://www.cirm.univ-rs.fr/videos/2006/exposes/02_LeBris/Chipot.pdf
![Page 22: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/22.jpg)
Example
● Mutate PHE to HIS
![Page 23: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/23.jpg)
Phenylalanine Histidine
![Page 24: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/24.jpg)
F to H SWNT dU
![Page 25: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/25.jpg)
Phenylalanine vs. Histidine
ΔGPHE = -2.734 kcal/mol
ΔGHIS = -1.743 kcal/mol H2N
CH C
CH2 OH
O
![Page 26: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/26.jpg)
● From amino acid – SWNT free energy profile, we can estimate ΔG1 and ΔG2 :
ΔG1 = -2.734 kcal/mol ΔG2 = -1.743 kcal/mol
ΔΔG = ΔG2 - ΔG1 = +0.991 kcal/mol+0.991 kcal/mol
● From alchemical transformation: ΔG1
alch = -28.58 kcal/mol ΔG2alch = -27.88 kcal/mol
ΔΔG’ = ΔG2alch -ΔG1
alch = +0.7 kcal/mol+0.7 kcal/mol
1alch
2alch12(aq)
ΔG(aq)
ΔGΔG
1alch2
2alch1(aq)
ΔG(aq)
ΔGΔGΔGΔG SWNT-HIS HIS
ΔGΔGΔGΔG SWNT-PHE PHE
2
2
alch
1
alch
1
−=− →
↓↓
+=+→
![Page 27: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/27.jpg)
● To be continued............
![Page 28: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/28.jpg)
Thermodynamic Integration
![Page 29: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/29.jpg)
H p1 , p2 , ....... pN , r1 , r 2 , .......... rN
Aa b=Ab−Aa =∫dAd
d
dAd
=
∫ ∂H r , p
∂exp−H r , pdr dp
∫exp−H r , pdr dp
Aa b=∫⟨∂H r , p
∂⟩d
⟨∂ H r , p
∂ ⟩=−⟨F ⟩=constraint force
Reaction Coordinate ξ
![Page 30: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/30.jpg)
● Constraint force : force required to constrain the system at a fix reaction coordinate ξ
![Page 31: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/31.jpg)
Example
● Translocation of C60 into DPPC bilayer● Reaction coordinate : Distance between the
COM of C60 and COM of DPPC membrane● 33 constrained points : 0 – 32 Å, 1 Å spacing● Harmonic constraint force ● k = 100 kcal/mol/Å2
![Page 32: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/32.jpg)
![Page 33: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/33.jpg)
Umbrella Sampling
● Sample with umbrella potential U'(x)
● Compute biased probability P'(x)
● Estimate unbiased free energy
● A x = −k BT ln P' x − U' x F
● F is undetermined but irrelevant
![Page 34: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/34.jpg)
● Complex surfaces have multiple barriers
● Need to know the free energy surface to know an efficient bias
● Harmonic biasing function
● Multiple simulations
● Put the minimum of the bias in a different place for each simulation (sampling windows)
● Estimate P'(x) for each simulation
● Combine results from all simulations
● From one simulation
A x = −k BT ln P' x − U' x F
● F depends on U'(x)
● Different simulations have different offsets
● Not obvious how to combine different windows
![Page 35: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/35.jpg)
![Page 36: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/36.jpg)
WHAM
● Weighted Histogram Analysis Method● Determines optimal F values for combining
simulations● Kumar, et al. J Comput Chem, 13, 1011-1021,
1992● Generalizations
– Multidimension reaction coordinates
– Multiple temperatures
![Page 37: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/37.jpg)
● Nsims
= number of simulations
● ni(x)= number of counts in histogram bin associated with x
● Ubias,i
, Fi = biasing potential and free energy shift from simulation i
● P(x) = best estimate of unbiased probability distribution
● Fi and P(x) are unknowns
● Solve by iteration to self consistency
![Page 38: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/38.jpg)
Running US simulation
● Choose the reaction coordinate● Choose the number of windows and the biasing
potential● Run the simulations● Compute time series for the value of the
reaction coordinate (histograms)● Apply the WHAM equations
![Page 39: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/39.jpg)
![Page 40: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/40.jpg)
Example
● Translocation of C60 into DPPC bilayer● Reaction coordinate : Distance between the
COM of C60 and COM of DPPC membrane● 9 windows : 0 – 32 Å, 4 Å spacing● Harmonic constraint force ● k = 0.5 kcal/mol/Å2
![Page 41: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/41.jpg)
![Page 42: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/42.jpg)
![Page 43: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/43.jpg)
![Page 44: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/44.jpg)
Other Methods
● Steered MD
● Adaptive Biased Force Method (ABF)
● Parallel Tempering (Replica exchange)
- overcome free energy barriers
● etc.
![Page 45: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/45.jpg)
Replica ExchangeOvercoming Free Energy Barrier
● Non-directed method
(no reaction coordinate)● How to sample unfavorable
states?● At high T, barriers are
easier to overcome.● Heat and cool the system
to push it over barriers to sample new configurations
![Page 46: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/46.jpg)
Replica ExchangeOvercoming Free Energy Barrier
● Launch simulations at different temperatures
● Swap configurations based on a monte carlo criterion
● This criterion guarantees that the lowest (target) temperature “trajectory” samples from the Boltzmann distribution.
● Swapping configurations effectively improves sampling
![Page 47: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/47.jpg)
Overcoming Free Energy Barrier
● When doing parallel tempering molecular dynamics, one must take care in the interpretation of the results. A parallel tempering exchange is an ‘unphysical’ move, and so one cannot draw conclusions about dynamics. That is, when using parallel tempering molecular dynamics, one is only really doing a form of sampling and not ‘true’ molecular dynamics.
Deem et al. Phys . Chem. Chem. Phys . 2 0 0 5 , 7 , 3 9 1 0 – 3 9 1 6
![Page 48: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/48.jpg)
Garcia et al. PROTEINS: Structure, Function, and Bioinformatics 59:783–790 (2005)
![Page 49: Free Energy Calculation in MD Simulationson051000/comp/FreeE.pdf · Overcoming Free Energy Barrier When doing parallel tempering molecular dynamics, one must take care in the interpretation](https://reader034.fdocument.pub/reader034/viewer/2022042707/5a76b12e7f8b9a93088d54ab/html5/thumbnails/49.jpg)
Reference
● Anna Johansson xray.bmc.uu.se/~calle/md_phd/free_energ.pdf
● David Mathews
http://rna.urmc.rochester.edu/teaching.html● Chris Chipot
http://www.cirm.univ-rs.fr/videos/2006/exposes/02_LeBris/Chipot.pdf