Alternative approaches to potential of mean force calculations: Free energy perturbation versus thermodynamic integration. The case study of hydrophobic interactions

1995 ◽  
Author(s):  
Christophe Chipot ◽  
Bernard Maigret ◽  
David A. Pearlman
Keyword(s):  
Free Energy ◽  
Hydrophobic Interactions ◽  
Potential Of Mean Force ◽  
Thermodynamic Integration ◽  
Free Energy Perturbation ◽  
Perturbation Versus ◽  
Alternative Approaches ◽  
Energy Perturbation

scholarly journals Alchemical Free Energy Estimators and Molecular Dynamics Engines: Accuracy, Precision and Reproducibility

2021 ◽  
Author(s):  
Alexander Wade ◽  
Agastya Bhati ◽  
Shunzhou Wan ◽  
Peter Coveney
Keyword(s):  
Molecular Dynamics ◽  
Free Energy ◽  
Binding Free Energy ◽  
Thermodynamic Integration ◽  
Free Energy Perturbation ◽  
Free Energies ◽  
Ensemble Averaging ◽  
Relative Binding ◽  
Energy Perturbation ◽  
Binding Free Energies

The binding free energy between a ligand and its target protein is an essential quantity to know at all stages of the drug discovery pipeline. Assessing this value computationally can offer insight into where efforts should be focused in the pursuit of effective therapeutics to treat myriad diseases. In this work we examine the computation of alchemical relative binding free energies with an eye to assessing reproducibility across popular molecular dynamics packages and free energy estimators. The focus of this work is on 54 ligand transformations from a diverse set of protein targets: MCL1, PTP1B, TYK2, CDK2 and thrombin. These targets are studied with three popular molecular dynamics packages: OpenMM, NAMD2 and NAMD3. Trajectories collected with these packages are used to compare relative binding free energies calculated with thermodynamic integration and free energy perturbation methods. The resulting binding free energies show good agreement between molecular dynamics packages with an average mean unsigned error between packages of 0.5 $kcal/mol$ The correlation between packages is very good with the lowest Spearman's, Pearson's and Kendall's tau correlation coefficient between two packages being 0.91, 0.89 and 0.74 respectively. Agreement between thermodynamic integration and free energy perturbation is shown to be very good when using ensemble averaging.

scholarly journals A computational study of barium blockades in the KcsA potassium channel based on multi-ion potential of mean force calculations and free energy perturbation

2013 ◽  
Vol 142 (4) ◽  
pp. 451-463 ◽  
Author(s):  
Christopher N. Rowley ◽  
Benoît Roux
Keyword(s):  
Free Energy ◽  
Computational Study ◽  
Potential Of Mean Force ◽  
Free Energy Perturbation ◽  
K Channel ◽  
Rate Theory ◽  
Crystallographic Structure ◽  
Energy Perturbation ◽  
Lock In Effect ◽  
Lock In

Electrophysiological studies have established that the permeation of Ba2+ ions through the KcsA K+-channel is impeded by the presence of K+ ions in the external solution, while no effect is observed for external Na+ ions. This Ba2+ “lock-in” effect suggests that at least one of the external binding sites of the KcsA channel is thermodynamically selective for K+. We used molecular dynamics simulations to interpret these lock-in experiments in the context of the crystallographic structure of KcsA. Assuming that the Ba2+ is bound in site S2 in the dominant blocked state, we examine the conditions that could impede its translocation and cause the observed “lock-in” effect. Although the binding of a K+ ion to site S1 when site S2 is occupied by Ba2+ is prohibitively high in energy (>10 kcal/mol), binding to site S0 appears to be more plausible (ΔG > 4 kcal/mol). The 2D potential of mean force (PMF) for the simultaneous translocation of Ba2+ from site S2 to site S1 and of a K+ ion on the extracellular side shows a barrier that is consistent with the concept of external lock-in. The barrier opposing the movement of Ba2+ is very high when a cation is in site S0, and considerably smaller when the site is unoccupied. Furthermore, free energy perturbation calculations show that site S0 is selective for K+ by 1.8 kcal/mol when S2 is occupied by Ba2+. However, the same site S0 is nonselective when site S2 is occupied by K+, which shows that the presence of Ba2+ affects the selectivity of the pore. A theoretical framework within classical rate theory is presented to incorporate the concentration dependence of the external ions on the lock-in effect.

Sign in / Sign up

Export Citation Format

Share Document