Pith. sign in

REVIEW 2 major objections 2 minor 95 references

Integrated Alchemical and Conformational Enhanced Sampling for Solvation Free Energy Calculations

T0 review · 2 major / 2 minor · reviewed 2026-08-16 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read The paper claims flexible-solute logP accuracy requires simultaneously enhanced sampling of the alchemical coordinate and solute conformations, which λ-OPES-d-AFED delivers, cutting mean absolute error from 0.75 to 0.30 log units.

desk verdict Solid lambda-OPES protocol with a convincing rigid-solute benchmark; the flexible-solute logP claim needs the d-AFED reweighting quantified before it can be trusted. read the letter →

arxiv 2608.12691 v1 pith:6CYSBBKY submitted 2026-08-13 physics.chem-ph physics.comp-ph

Gabriela B. Correa , Charlles R. A. Abreu , Nisanth N. Nair , Mark E. Tuckerman This is my paper · ORCID
classification physics.chem-phphysics.comp-ph
keywords solvationfreeenergylambda-dynamicsalchemicalenhancedsamplingOPESd-AFEDoctanol-waterpartitioncoefficientconformationalmassscaling
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper addresses a known bottleneck in molecular simulation: solvation free energies converge slowly because three families of motions—reorganization of the surrounding solvent, exploration of the alchemical coordinate $\lambda$ that switches solute-solvent interactions on and off, and internal rotations of a flexible solute—relax on timescales far beyond ordinary simulation lengths. The authors claim all three can be handled in a single trajectory by combining solvent mass scaling, an adaptive OPES bias on $\lambda$, and d-AFED, which drives selected dihedral angles by coupling them to hot auxiliary variables. For rigid solutes they show that cutting octanol atomic masses tenfold accelerates Hamiltonian replica exchange convergence by more than fivefold without changing the equilibrium free energy, and that a dual-bias $\lambda$-OPES protocol matches replica exchange on comparable timescales without predefined $\lambda$ windows. For seven flexible $N$-acetyl amino-acid amides, adding d-AFED to $\lambda$-OPES lowers the mean absolute error of octanol-water partition coefficients from 0.75 to 0.30 log units relative to experiment, and shifts the predicted partitioning of Ac-Leu-NH$_2$ and Ac-Tyr-NH$_2$ into the experimentally observed water-favoring regime. The paper concludes that for flexible molecules, accurate partition coefficients require simultaneous alchemical and conformational sampling, and that the combined protocol provides a practical single-trajectory route to that goal.

What carries the argument

The machinery is an extended Hamiltonian in which the alchemical coordinate $\lambda$ is a dynamical variable propagated alongside the atomic coordinates, with a mirror-periodic mapping $\theta \in [0,2] \to \lambda \in [0,1]$ that lets the trajectory traverse the alchemical pathway in both directions. The alchemical coupling uses the linear basis function (LBF) formulation, where $\lambda$ enters only through switching functions $h_A(\lambda)$ and $h_B(\lambda)$ that smoothly turn on the full solute-solvent interaction and remove a softened Lennard-Jones capping potential, avoiding soft-core nonlinearities. OPES supplies the bias along $\lambda$ as an adaptive potential built from a real-time kernel density estimate; the dual-bias variant applies OPES-Standard and OPES-Explore simultaneously, the former flattening the free energy profile and the latter driving barrier crossing between the end states. d-AFED accelerates conformational sampling by coupling each selected dihedral $\xi_k$ to a harmonic auxiliary variable $s_k$ thermostatted at 1000 K with a large fictitious mass, chosen to preserve adiabatic separation from the physical degrees of freedom. The free energy is recovered by thermodynamic integration of the $\lambda$-conditional mean force, which the paper argues is unaffected by the marginal sampling changes introduced by the biases and the auxiliary thermostat.

What would settle it

A reader could settle the question by re-analyzing the stored d-AFED trajectories of the seven $N$-acetyl amino-acid amides with the explicit torsional reweighting procedure for d-AFED, and comparing the reweighted solvation free energies with the un-reweighted values reported here; if the differences exceed the reported statistical uncertainties, the 0.30 log-unit mean absolute error is contaminated by the auxiliary thermostat. A second decisive test is to repeat the $\lambda$-OPES-d-AFED simulations at auxiliary temperatures of 500 K and 2000 K: if the predicted partition coefficients shift by more than the statistical error, the results depend on the nonphysical thermostat setting rather than on a converged conformational equilibrium.

Watch

Extended reading notes

Core claim

The central claim is that the three coupled slow variables in solvation free energy calculations—solvent relaxation, the alchemical coordinate, and solute conformational modes—can be sampled simultaneously and accurately within one extended-dynamics simulation. The paper establishes this in stages: mass scaling of octanol by a factor of ten accelerates convergence of Hamiltonian replica exchange by more than fivefold while leaving equilibrium free energies unchanged; a dual-bias $\lambda$-OPES protocol, in which OPES-Standard flattens the free energy profile along $\lambda$ and a weaker OPES-Explore sustains transitions between the end states, reproduces HREX results in a single trajectory; and, for flexible $N$-acetyl amino-acid amides, coupling $\lambda$-OPES with d-AFED on backbone and side-chain dihedrals lowers the mean absolute error of predicted octanol-water partition coefficients from 0.75 to 0.30 log units relative to experiment, including correcting the sign of partitioning for Ac-Leu-NH$_2$ and Ac-Tyr-NH$_2$. The paper interprets these results as demonstrating that simultaneous treatment of alchemical and conformational sampling is required for accurate partition coefficients of flexible molecules, and that the $\lambda$-OPES-d-AFED protocol supplies that treatment in a single simulation.

Load-bearing premise

The load-bearing premise is that the auxiliary torsional variables, thermostatted at 1000 K with large fictitious masses, leave the $\lambda$-conditional mean force unbiased, so that solvation free energies can be read directly from the biased trajectory; the paper states that reweighting produced negligible changes without showing a quantitative comparison, and if that adiabatic-separation assumption silently fails, the improved mean absolute error of 0.30 log units could be an artifact of distorted torsional populations rather than a genuine conformational-sampling correction.

Editorial extensions

If this is right

  • If the central claim is right, solvation free energy and partition coefficient calculations on flexible solutes should treat conformational sampling as a first-class requirement; otherwise errors on the order of half a log unit may be misattributed to force-field inaccuracy.
  • The dual-bias $\lambda$-OPES protocol matches Hamiltonian replica exchange results on comparable timescales while running as a single trajectory with no predefined $\lambda$ windows and no replica coordination, which the paper notes is advantageous when replica-based approaches are computationally demanding.
  • Tenfold mass scaling of octanol speeds convergence more than fivefold without altering the equilibrium free energy, providing a cheap acceleration that can be combined with any sampling strategy in viscous organic solvents.
  • The protocol's validation on seven amino-acid amides establishes a foundation the paper proposes for applying alchemical free energy methods to larger and more conformationally complex solutes.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the pattern in this paper generalizes, a portion of the systematic error in computed octanol-water partition coefficients of flexible drug-like molecules that is currently blamed on force fields is really incomplete torsional sampling; comparing $\lambda$-OPES-d-AFED against plain $\lambda$-OPES on a diverse drug-like panel would test this directly.
  • The paper reports the d-AFED reweighting correction as negligible without quantitative support; a natural extension is to report reweighted and un-reweighted free energies explicitly for solutes with higher torsional barriers, where the 1000 K auxiliary population can deviate more strongly from the 298 K physical one.
  • The tenfold solvent mass-scaling result was established in octanol; applying the same rule to other viscous solvents (for example DMSO or ethylene glycol) would test whether the acceleration is robust beyond this one solvent.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

2 major / 2 minor

Summary. The paper proposes an integrated enhanced sampling protocol for solvation free energy calculations. It combines (i) solvent mass scaling to accelerate relaxation in viscous phases such as octanol, (ii) a dual-bias lambda-dynamics scheme (“λ-OPES”) that applies OPES-Standard and OPES-Explore simultaneously on the alchemical coordinate, and (iii) driven adiabatic free energy dynamics (d-AFED) on selected solute dihedrals to accelerate conformational sampling. The manuscript reports that for seven rigid organic solutes, λ-OPES agrees with Hamiltonian replica exchange (HREX) and literature values (MAE below 0.06 kcal/mol in water and 0.17 kcal/mol in octanol), and that a tenfold reduction of octanol atomic masses accelerates convergence by more than fivefold without changing equilibrium free energies. For seven flexible N-acetyl amino-acid amides, the paper claims that coupling λ-OPES with d-AFED reduces the mean absolute error in octanol-water logP from 0.75 to 0.30 log units relative to experiment. The central methodological claims are the single-trajectory capability of λ-OPES and the necessity of simultaneous alchemical and conformational enhanced sampling for flexible solutes.

Significance. If the flexible-solute claim is correct, this is a practically valuable contribution: it offers a single-trajectory route to solvation free energies with coupled alchemical and conformational sampling, avoiding the need for predefined λ windows or replica coordination. The rigid-solute results are convincing and well controlled: λ-OPES reproduces HREX reference values and external literature values, and the mass-scaling acceleration in octanol is a broadly useful finding that is made credible by the consistent convergence data. The paper also provides open-source simulation scripts and a clear description of the LBF alchemical route, which strengthens reproducibility. However, the headline flexible-solute claim rests on a single, hand-picked set of d-AFED and OPES parameters for seven compounds, and on an unquantified statement that reweighting the d-AFED torsional ensemble changes solvation free energies negligibly. These points need to be substantiated before the central claim is established.

major comments (2)
  1. [Section II.A and Section III.D] The assertion that reweighting the d-AFED torsional ensemble produces “only negligible changes in the solvation free energies” is unquantified and load-bearing. With κ = 1000 kJ mol^-1 rad^-2 and T_s = 1000 K (Section II.F), the conditional width of (ξ_k − s_k) at 298 K is about sqrt(k_B T/κ) ≈ 0.05 rad, so each physical dihedral is effectively slaved to its auxiliary variable, whose distribution is thermostatted at 1000 K. The sampled torsional population is therefore far from the 298 K physical population. Because Figures 9A/B show that ΔG changes by 0.5–1.2 kcal/mol when the torsional sampling protocol is changed, the mean force in Eq. (2) is demonstrably sensitive to the torsional ensemble. The statement in Section II.A that the conditional mean force is unaffected by the auxiliary temperature conflates the OPES bias (which depends only on λ and cancels in the conditional average) with the d-AFED coupling (which depends on r through ξ(r) and does not cancel). The authors must report a quantitative comparison of ΔG obtained with and without the Cuendet–Tuckerman reweighting for each of the seven solutes in water and octanol. If reweighting is not actually negligible, the MAE reduction from 0.75 to 0.30 log units (Figure 10) may reflect a high-temperature conformational ensemble rather than unbiased 298 K conformational sampling.
  2. [Section II.F and Section III.D] The flexible-solute result (MAE 0.30) is established for a single hand-picked parameter set: d-AFED masses of 50 Da nm^2 rad^-2, coupling constants of 1000 kJ mol^-1 rad^-2, auxiliary temperature 1000 K, and OPES barriers of 50 and 5 kcal/mol. No sensitivity analysis is presented for these parameters in the peptide-like systems. This matters because the authors themselves show in Section III.B (Figures S5 and S6) that the closely related λ-AFED results are sensitive to the extended-variable mass. A demonstration that the 0.75→0.30 improvement is robust to, e.g., varying T_s between 800 and 1200 K or κ by a factor of two, is needed to exclude the possibility that the improvement is a result of parameter tuning on a seven-compound set.
minor comments (2)
  1. [Section II.A] The sentence describing the mirror-periodic representation of λ cites Wu, Hu, and Yang; a brief explanation of why the generalized force changes sign in the interval (1,2] would improve accessibility.
  2. [Table S3] The superscripts in the tabulated ΔG values (presumably standard errors over the three runs) are not defined in the table caption or in the text; please add an explicit legend.

Circularity Check

0 steps flagged · score 1.0 of 10

No significant circularity: the d-AFED and LBF components come from independently published derivations, and the headline MAE reduction is benchmarked against external experimental and HREX data rather than fitted inputs.

full rationale

The paper's central derivation is a standard thermodynamic-integration route in which the solvation free energy is obtained by integrating the conditional mean force along lambda, and the enhancement methods (OPES on lambda, d-AFED on torsions, solvent mass scaling) are applied as biasing or dynamics-acceleration devices. The load-bearing assertion that the d-AFED auxiliary temperature and the OPES bias leave the conditional mean force unbiased is supported by the authors' earlier d-AFED theory (refs. 71 and 73), but those are parameter-free published derivations whose stated assumptions do not include the target logP values, so under the review rules they count as independent support rather than circular self-citation. The LBF alchemical model is likewise taken from a prior independent publication (ref. 21). The headline result, a reduction in mean absolute error from 0.75 to 0.30 log units against experimental octanol-water partition coefficients, is an external benchmark comparison and is not obtained by fitting any parameter to those experimental values; the reported OPES barriers, d-AFED masses, coupling constants, and auxiliary temperatures are set by stated physical criteria or taken from earlier method papers. The only notable weakness is the unquantified statement that torsional reweighting produced 'negligible changes' in the solvation free energies; this is an omitted numerical support detail and a possible correctness risk, but it is not a circular step because the reweighting correction is an empirical check on the reported free energies rather than an input reused as the output. Accordingly, no specific circular reduction can be exhibited, and the appropriate finding is essentially non-circular.

Assumptions & free parameters 7 free parameters · 4 assumptions · 0 invented entities

The protocol rests on standard statistical-mechanics results and on domain assumptions about force field accuracy and d-AFED adiabatic separation. The hand-selected OPES and d-AFED parameters are the main free inputs; none are fitted to the experimental partition coefficients, keeping circularity low but leaving the quantitative MAE claims sensitive to parameter choices.

free parameters (7)
  • OPES barrier DeltaE_standard = 30 kcal/mol (rigid solutes), 50 kcal/mol (peptide solutes)
    Hand-chosen in Section II.F; controls the bias factor gamma; no sensitivity analysis reported.
  • OPES barrier DeltaE_explore = 5 kcal/mol (dual-bias protocol)
    Hand-chosen as the weaker exploratory component; no sensitivity analysis reported.
  • d-AFED auxiliary mass m_s = 50 Da nm2 rad-2
    Chosen to be about 300 times the effective torsional mass following Cuendet et al.; no explicit check shown.
  • d-AFED coupling constant kappa = 1000 kJ mol-1 rad-2
    Targets an approximate 40 fs oscillation period for the coupled torsional motion.
  • d-AFED auxiliary temperature T_s = 1000 K
    Taken from previous d-AFED applications to promote barrier crossing while preserving adiabatic separation.
  • octanol mass scaling factor = 10
    Largest stable factor at a 1 fs time step; scaling beyond 10 led to numerical instabilities.
  • lambda mass m_lambda = 1 Da nm2
    Shown insensitive for lambda-OPES on quinone in water over 1 to 1e3 Da nm2, but still hand-selected for all other runs.
assumptions (4)
  • standard math Thermodynamic integration of the conditional mean force gives the solvation free energy under OPES and d-AFED biases (Eq. 2).
    Invoked in Section II.A with citations 71 and 73; relies on adiabatic separation between physical and auxiliary variables.
  • standard math Changing atomic masses does not change equilibrium configurational properties.
    Used in Section III.A to justify mass scaling; standard classical statistical mechanics.
  • domain assumption Force fields GAFF/TIP3P, AMBER14SB/OPC/GAFF2 describe water and octanol solvation well enough for the accuracy claims.
    Stated in Section II.F and used to interpret MAE against experimental logP values.
  • ad hoc to paper The d-AFED torsional bias changes solvation free energies negligibly, so no reweighting is needed.
    Self-reported in Section II.A as a negligible correction; no quantitative comparison is shown.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Integrated Alchemical and Conformational Enhanced Sampling for Solvation Free Energy Calculations." pith.science (2026). https://pith.science/paper/6CYSBBKY

@misc{pith2026260812691,
  author       = {Pith},
  title        = {Pith review of: Integrated Alchemical and Conformational Enhanced Sampling for Solvation Free Energy Calculations},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/6CYSBBKY}},
  note         = {Machine review of arXiv:2608.12691}
}
abstract

Accurate solvation free energies from molecular dynamics simulations require efficient sampling of coupled slow variables, including solvent coordinates, solute conformational modes, and the alchemical coordinate $\lambda$. Here, we develop a $\lambda$-dynamics framework that combines mass scaling, on-the-fly probability enhanced sampling (OPES), and driven adiabatic free energy dynamics (d-AFED) to address these sampling challenges within a unified protocol. For rigid organic solutes, Hamiltonian replica exchange with mass scaling is first used to quantify the effect of octanol solvent relaxation. Reducing all octanol atomic masses by a factor of ten accelerates convergence by more than fivefold while preserving equilibrium solvation free energies. These calculations then provide reference benchmarks for $\lambda$-OPES, a dual-bias $\lambda$-dynamics strategy that combines the "standard" and "explore" variants of OPES to promote transitions along the alchemical coordinate. This approach reaches convergence on timescales comparable to replica exchange, but without predefined $\lambda$ windows or multiple parallel simulations. For flexible $N$-acetyl amino-acid amide solutes, $\lambda$-OPES is coupled with d-AFED on selected backbone and side-chain dihedrals to enable simultaneous alchemical and conformational enhanced sampling. This combined strategy improves agreement with experimental octanol-water partition coefficients and reduces the mean absolute error from 0.75 log units with $\lambda$-OPES alone to 0.30 log units with $\lambda$-OPES-d-AFED. Overall, this work establishes an integrated enhanced sampling protocol for solvation free energy calculations across rigid organic solutes and flexible peptide-like solutes, and provides a foundation for the application of alchemical free energy methods to larger and more conformationally complex systems.

Figures

Figures reproduced from arXiv: 2608.12691 by the authors.

Figure 1
Figure 1. FIG. 1. Chemical structures of the rigid organic solutes. [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. Chemical structures of the flexible peptide-like solutes. [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3. Convergence of solvation free energies for (A) quinone [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figures from the paper (5 more)
Figure 5
Figure 5. Figure 5: FIG. 5. Cumulative free energy, ∆ [PITH_FULL_IMAGE:figures/full_fig_p008_5.png]
Figure 7
Figure 7. Figure 7: FIG. 7. Octanol-water partition coefficients for the seven rigid organic solutes obtained from HREX and [PITH_FULL_IMAGE:figures/full_fig_p009_7.png]
Figure 8
Figure 8. Figure 8: FIG. 8. Time evolution of the [PITH_FULL_IMAGE:figures/full_fig_p009_8.png]
Figure 9
Figure 9. Figure 9: FIG. 9. Comparison of solvation free energies obtained from [PITH_FULL_IMAGE:figures/full_fig_p010_9.png]
Figure 10
Figure 10. Figure 10: FIG. 10. Octanol-water partition coefficients for acetyl amino-acid amides obtained from [PITH_FULL_IMAGE:figures/full_fig_p011_10.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

95 extracted references · 57 canonical work pages

  1. [1]

    2016 , journal =

    Calculating Partition Coefficients of Small Molecules in Octanol/Water and Cyclohexane/Water , author =. 2016 , journal =. doi:10.1021/acs.jctc.6b00449 , language =

    work page doi:10.1021/acs.jctc.6b00449 2016

  2. [2]

    1995 , publisher =

    Exploring QSAR: Hydrophobic, Electronic, and Steric Constants , author =. 1995 , publisher =. doi:, language =

    1995

  3. [3]

    1989 , journal =

    Octanol-Water Partition Coefficients of Simple Organic Compounds , author =. 1989 , journal =. doi:10.1063/1.555833 , language =

    work page doi:10.1063/1.555833 1989

  4. [4]

    1971 , journal =

    Partition coefficients and their uses , author =. 1971 , journal =. doi:10.1021/cr60274a001 , language =

    work page doi:10.1021/cr60274a001 1971

  5. [5]

    1983 , journal =

    Hydrophobic parameters pi of amino-acid side chains from the partitioning of N-acetyl-amino-acid amides , author =. 1983 , journal =. doi:, language =

    1983

  6. [6]

    1992 , journal =

    Amino Acid Side-Chain Contributions to Free Energy of Transfer of Tripeptides from Water to Octanol , author =. 1992 , journal =. doi:10.1023/a:1015892313856 , language =

    work page doi:10.1023/a:1015892313856 1992

  7. [7]

    2004 , journal =

    Additive Transfer Free Energies of the Peptide Backbone Unit That Are Independent of the Model Compound and the Choice of Concentration Scale , author =. 2004 , journal =. doi:10.1021/bi035908r , language =

    work page doi:10.1021/bi035908r 2004

  8. [8]

    2013 , journal =

    Solvation Free Energy of the Peptide Group: Its Model Dependence and Implications for the Additive-Transfer Free-Energy Model of Protein Stability , author =. 2013 , journal =. doi:10.1016/j.bpj.2013.08.011 , language =

    work page doi:10.1016/j.bpj.2013.08.011 2013

Show all 95 references
  1. [9]

    2025 , journal =

    Predicting Distribution Coefficients (LogD) of Cyclic Peptides Using Molecular Dynamics Simulations , author =. 2025 , journal =. doi:10.1007/s11095-025-03850-2 , language =

    2025 doi

  2. [10]

    2026 , journal =

    Modeling Method Based on Traditional Machine Learning and Graph Neural Network for Interpretable Prediction of Short Peptide logD , author =. 2026 , journal =. doi:10.1021/acs.jpcb.6c00299 , language =

    2026 doi

  3. [11]

    2012 , journal =

    The influence of lipophilicity in drug discovery and design , author =. 2012 , journal =. doi:10.1517/17460441.2012.714363 , language =

    2012

  4. [12]

    2018 , journal =

    Lipophilic Efficiency as an Important Metric in Drug Design , author =. 2018 , journal =. doi:10.1021/acs.jmedchem.8b00077 , language =

    2018 doi

  5. [13]

    2016 , journal =

    Quality guidelines for oral drug candidates: dose, solubility and lipophilicity , author =. 2016 , journal =. doi:10.1016/j.drudis.2016.07.007 , language =

    2016 doi

  6. [14]

    2025 , journal =

    Therapeutic Peptides: Recent Advances in Discovery, Synthesis, and Clinical Translation , author =. 2025 , journal =. doi:10.3390/ijms26115131 , language =

    2025 doi

  7. [15]

    2022 , journal =

    Therapeutic peptides: current applications and future directions , author =. 2022 , journal =. doi:10.1038/s41392-022-00904-4 , language =

    2022 doi

  8. [16]

    2023 , journal =

    Recent advances in the development of therapeutic peptides , author =. 2023 , journal =. doi:10.1016/j.tips.2023.04.003 , language =

    2023 doi

  9. [17]

    2018 , journal =

    Calculation of the water-octanol partition coefficient of cholesterol for SPC, TIP3P, and TIP4P water , author =. 2018 , journal =. doi:10.1063/1.5054056 , language =

    2018 doi

  10. [18]

    2019 , journal =

    Calculating the Partition Coefficients of Organic Solvents in Octanol/Water and Octanol/Air , author =. 2019 , journal =. doi:10.1021/acs.jcim.9b00212 , language =

    2019 doi

  11. [19]

    2020 , journal =

    Prediction of octanol‑water partition coefficients for the SAMPL6‑logP molecules using molecular dynamics simulations with OPLS‑AA, AMBER and CHARMM force fields , author =. 2020 , journal =. doi:10.1007/s10822-019-00267-z , language =

    2020 doi

  12. [20]

    2021 , journal =

    Predicting octanol/water partition coefficients using molecular simulation for the SAMPL7 challenge: comparing the use of neat and water saturated 1-octanol , author =. 2021 , journal =. doi:10.1007/s10822-021-00415-4 , language =

    2021 doi

  13. [21]

    2015 , journal =

    Guidelines for the analysis of free energy calculations , author =. 2015 , journal =. doi:10.1007/s10822-015-9840-9 , language =

    2015 doi

  14. [22]

    2014 , journal =

    Practical Aspects of Free-Energy Calculations: A Review , author =. 2014 , journal =. doi:10.1021/ct500161f , language =

    2014 doi

  15. [23]

    Abreu, Charlles R. A. and Tuckerman, Mark E. , year =. Molecular. Journal of Chemical Theory and Computation , volume =. doi:10.1021/acs.jctc.0c00698 , language =

    doi

  16. [24]

    and Pham, Tri T

    Naden, Levi N. and Pham, Tri T. and Shirts, Michael R. , year =. Linear. Journal of Chemical Theory and Computation , volume =. doi:10.1021/ct4009188 , language =

    doi

  17. [25]

    and Shirts, Michael R

    Naden, Levi N. and Shirts, Michael R. , year =. Linear. Journal of Chemical Theory and Computation , volume =. doi:10.1021/ct501047e , language =

    doi

  18. [26]

    Journal of Chemical Theory and Computation , volume=

    A new formulation for the concerted alchemical calculation of van der Waals and Coulomb components of solvation free energies , author=. Journal of Chemical Theory and Computation , volume=. 2022 , doi=

    2022

  19. [27]

    Journal of Chemical Theory and Computation , volume=

    Computing Alchemical Free Energy Differences with Hamiltonian Replica Exchange Molecular Dynamics (H-REMD) Simulations , author=. Journal of Chemical Theory and Computation , volume=. 2011 , doi=

    2011

  20. [28]

    Journal of Chemical Physics , volume=

    Replica exchange and expanded ensemble simulations as Gibbs sampling: Simple improvements for enhanced mixing , author=. Journal of Chemical Physics , volume=. 2011 , doi=

    2011

  21. [29]

    1996 , journal =

    -dynamics: A new approach to free energy calculations , author =. 1996 , journal =. doi:10.1063/1.472109 , language =

    1996 doi

  22. [30]

    2017 , journal =

    Adaptive Landscape Flattening Accelerates Sampling of Alchemical Space in Multisite Dynamics , author =. 2017 , journal =. doi:10.1021/acs.jpcb.6b09656 , language =

    2017 doi

  23. [31]

    2024 , journal =

    How to Sample Dozens of Substitutions per Site with Dynamics , author =. 2024 , journal =. doi:10.1021/acs.jctc.4c00514 , language =

    2024 doi

  24. [32]

    2021 , journal =

    BLaDE: A Basic Lambda Dynamics Engine for GPU-Accelerated Molecular Dynamics Free Energy Calculations , author =. 2021 , journal =. doi:10.1021/acs.jctc.1c00833 , language =

    2021 doi

  25. [33]

    2006 , journal =

    Efficient and precise solvation free energies via alchemical adiabatic molecular dynamics , author =. 2006 , journal =. doi:10.1063/1.2232082 , language =

    2006 doi

  26. [34]

    2011 , journal =

    ‑Metadynamics Approach To Compute Absolute Solvation Free Energy , author =. 2011 , journal =. doi:10.1021/jz200808x , language =

    2011 doi

  27. [35]

    2014 , journal =

    Local Elevation Umbrella Sampling Applied to the Calculation of Alchemical Free-Energy Changes via ‑Dynamics: The ‑LEUS Scheme , author =. 2014 , journal =. doi:10.1021/ct5002686 , language =

    2014 doi

  28. [36]

    2024 , journal =

    Lambda-ABF: Simplified, Portable, Accurate, and Cost-Effective Alchemical Free-Energy Computation , author =. 2024 , journal =. doi:10.1021/acs.jctc.3c01249 , language =

    2024 doi

  29. [37]

    Rethinking

    Invernizzi, Michele and Parrinello, Michele , year =. Rethinking. Journal of Physical Chemistry Letters , volume =

  30. [38]

    and Parrinello, Michele , year =

    Invernizzi, Michele and Piaggi, Pablo M. and Parrinello, Michele , year =. Unified. Physical Review X , volume =

  31. [39]

    Exploration vs

    Invernizzi, Michele and Parrinello, Michele , year =. Exploration vs. Journal of Chemical Theory and Computation , volume =

  32. [40]

    2025 , journal =

    Lambda-ABF-OPES: Faster Convergence with High Accuracy in Alchemical Free Energy Calculations , author =. 2025 , journal =

    2025

  33. [41]

    2026 , journal =

    Dual-LAO for calculating fast and robust relative binding free energies of simple and complex transformations , author =. 2026 , journal =

    2026

  34. [42]

    2024 , journal =

    Eastman, Peter and Galvelis, Raimondas and Pel. 2024 , journal =

    2024

  35. [43]

    2023 , journal =

    Alchemical Metadynamics: Adding Alchemical Variables to Metadynamics to Enhance Sampling in Free Energy Calculations , author =. 2023 , journal =. doi:10.1021/acs.jctc.2c01258 , language =

    2023 doi

  36. [44]

    2013 , journal =

    Absolute Hydration Free Energies of Blocked Amino Acids: Implications for Protein Solvation and Stability , author =. 2013 , journal =. doi:10.1016/j.bpj.2012.12.008 , language =

    2013 doi

  37. [45]

    2010 , journal =

    Backbone additivity in the transfer model of protein solvation , author =. 2010 , journal =. doi:10.1002/pro.378 , language =

    2010 doi

  38. [46]

    2023 , journal =

    Solvent quality and solvent polarity in polypeptides , author =. 2023 , journal =. doi:10.1039/D2CP05214H , language =

    2023 doi

  39. [47]

    2024 , journal =

    Solvent quality and nonbiological oligomer folding: revisiting conventional paradigms , author =. 2024 , journal =. doi:10.1039/D4SM00727A , language =

    2024 doi

  40. [48]

    2015 , journal =

    Solvation thermodynamics of amino acid side chains on a short peptide backbone , author =. 2015 , journal =. doi:10.1063/1.4917076 , language =

    2015 doi

  41. [49]

    2019 , journal =

    Conformational perturbation of peptides in presence of polar organic solvents , author =. 2019 , journal =. doi:10.1016/j.jmgm.2019.02.009 , language =

    2019 doi

  42. [50]

    2016 , journal =

    Solvation Thermodynamics of Oligoglycine with Respect to Chain Length and Flexibility , author =. 2016 , journal =. doi:10.1016/j.bpj.2016.07.013 , language =

    2016 doi

  43. [51]

    2014 , journal =

    Hydration free energies calculated using the AMBER ff03 charge model for natural and unnatural amino acids and multiple water models , author =. 2014 , journal =. doi:10.1016/j.compchemeng.2014.07.017 , language =

    2014 doi

  44. [52]

    2016 , journal =

    Enthalpic and Entropic Contributions to Hydrophobicity , author =. 2016 , journal =. doi:10.1021/acs.jctc.6b00422 , language =

    2016 doi

  45. [53]

    2017 , journal =

    Solvent effects on ligand binding to a serine protease , author =. 2017 , journal =. doi:10.1039/C6CP07899K , language =

    2017 doi

  46. [54]

    2021 , journal =

    Therapeutic peptides for the treatment of cystic fibrosis: Challenges and perspectives , author =. 2021 , journal =. doi:10.1016/j.ejmech.2021.113191 , language =

    2021

  47. [55]

    2026 , journal =

    One for All, All for One: A Unified Framework for Free-Energy Calculations , author =. 2026 , journal =. doi:10.1021/acs.accounts.5c00666 , language =

    2026 doi

  48. [56]

    2024 , journal =

    Enhanced and Efficient Predictions of Dynamic Ionization through Constant-pH Adiabatic Free Energy Dynamics , author =. 2024 , journal =. doi:10.1021/acs.jctc.4c00704 , language =

    2024 doi

  49. [57]

    2025 , journal =

    Solute Tempered Adiabatic Free Energy Dynamics for Enhancing Conformational Space Sampling , author =. 2025 , journal =. doi:10.1021/acs.jctc.5c00717 , language =

    2025 doi

  50. [58]

    2023 , journal =

    An interoperable implementation of collective-variable based enhanced sampling methods in extended phase space within the OpenMM package , author =. 2023 , journal =. doi:10.1002/jcc.27182 , language =

    2023 doi

  51. [59]

    2002 , journal =

    On the use of the adiabatic molecular dynamics technique in the calculation of free energy profiles , author =. 2002 , journal =. doi:10.1063/1.1448491 , language =

    2002 doi

  52. [60]

    2024 , journal =

    Faster Sampling in Molecular Dynamics Simulations with TIP3P-F Water , author =. 2024 , journal =. doi:10.1021/acs.jctc.4c00990 , language =

    2024 doi

  53. [61]

    2010 , journal =

    Enhanced Conformational Sampling of Peptides via Reduced Side-Chain and Solvent Masses , author =. 2010 , journal =. doi:10.1021/jp109865y , language =

    2010 doi

  54. [62]

    2026 , journal =

    Accelerating Replica Exchange Molecular Dynamics: A Comparison of Hydrogen Mass Repartitioning and Light Water Models , author =. 2026 , journal =. doi:10.1021/acs.jctc.5c01929 , language =

    2026 doi

  55. [63]

    2015 , journal =

    Long-Time-Step Molecular Dynamics through Hydrogen Mass Repartitioning , author =. 2015 , journal =. doi:10.1021/ct5010406 , language =

    2015 doi

  56. [64]

    2015 , journal =

    Low-mass molecular dynamics simulation for configurational sampling enhancement: More evidence and theoretical explanation , author =. 2015 , journal =. doi:10.1016/j.bbrep.2015.08.023 , language =

    2015 doi

  57. [65]

    2024 , journal =

    Densities, Viscosities, and Self-Diffusion Coefficients of Octan-1-ol and Related Ether-Alcohols , author =. 2024 , journal =. doi:10.1021/acs.jced.4c00195 , language =

    2024 doi

  58. [66]

    2008 , journal =

    Statistically optimal analysis of samples from multiple equilibrium states , author =. 2008 , journal =

    2008

  59. [67]

    1935 , journal =

    Statistical mechanics of fluid mixtures , author =. 1935 , journal =. doi:10.1063/1.1749657 , language =

    1935 doi

  60. [68]

    2020 , journal =

    Improved Alchemical Free Energy Calculations with Optimized Smoothstep Softcore Potentials , author =. 2020 , journal =. doi:10.1021/acs.jctc.0c00237 , language =

    2020 doi

  61. [69]

    2023 , journal =

    OneOPES, a Combined Enhanced Sampling Method to Rule Them All , author =. 2023 , journal =. doi:10.1021/acs.jctc.3c00254 , language =

    2023 doi

  62. [70]

    2017 , journal =

    Exploring high dimensional free energy landscapes: Temperature accelerated sliced sampling , author =. 2017 , journal =. doi:10.1063/1.4977704 , language =

    2017 doi

  63. [71]

    2015 , journal =

    Efficient Sampling of High-Dimensional Free-Energy Landscapes with Parallel Bias Metadynamics , author =. 2015 , journal =. doi:10.1021/acs.jctc.5b00846 , language =

    2015 doi

  64. [72]

    2021 , journal =

    Hamiltonian based resonance-free approach for enabling very large time steps in multiple time-scale molecular dynamics , author =. 2021 , journal =. doi:10.1080/00268976.2021.1923848 , language =

    2021

  65. [73]

    2004 , journal =

    Molecular dynamics simulations of water and biomolecules with a Monte Carlo constant pressure algorithm , author =. 2004 , journal =. doi:10.1016/j.cplett.2003.12.039 , language =

    2004 doi

  66. [74]

    2009 , journal =

    Packmol: A package for building initial configurations for molecular dynamics simulations , author =. 2009 , journal =. doi:10.1002/jcc.21224 , language =

    2009 doi

  67. [75]

    1983 , month = jul, journal =

    Comparison of Simple Potential Functions for Simulating Liquid Water , author =. 1983 , month = jul, journal =. doi:10.1063/1.445869 , language =

    1983 doi

  68. [76]

    2004 , journal =

    Development and testing of a general AMBER force field , author =. 2004 , journal =. doi:10.1002/jcc.20035 , language =

    2004 doi

  69. [77]

    and Jack, David B

    Jakalian, Araz and Bush, Bruce L. and Jack, David B. and Bayly, Christopher I. , year =. Fast, Efficient Generation of High-Quality Atomic Charges. Journal of Computational Chemistry , volume =

  70. [78]

    Particle Mesh

    Darden, Tom and York, Darrin and Pedersen, Lee , year =. Particle Mesh. The Journal of Chemical Physics , volume =. doi:10.1063/1.464397 , language =

    doi

  71. [79]

    1995 , journal =

    A smooth particle mesh Ewald method , author =. 1995 , journal =. doi:10.1063/1.470117 , language =

    1995 doi

  72. [80]

    2014 , journal =

    Building Water Models: A Different Approach , author =. 2014 , journal =. doi:10.1021/jz501780a , language =

    2014 doi

  73. [81]

    2015 , journal =

    ff14SB: Improving the Accuracy of Protein Side Chain and Backbone Parameters from ff99SB , author =. 2015 , journal =. doi:10.1021/acs.jctc.5b00255 , language =

    2015 doi

  74. [82]

    2003 , journal =

    Generalized ensembles serve to improve the convergence of free energy simulations , author =. 2003 , journal =. doi:10.1016/S0009-2614(03)01057-1 , language =

    2003 doi

  75. [83]

    2014 , journal =

    Free Energy Reconstruction from Metadynamics or Adiabatic Free Energy Dynamics Simulations , author =. 2014 , journal =. doi:10.1021/ct500012b , language =

    2014 doi

  76. [84]

    2012 , journal =

    Heating and flooding: A unified approach for rapid generation of free energy surfaces , author =. 2012 , journal =. doi:10.1063/1.4733389 , language =

    2012 doi

  77. [85]

    2020 , journal =

    Best Practices for Alchemical Free Energy Calculations , author =. 2020 , journal =. doi:10.33011/livecoms.2.1.18378 , language =

    2020 doi

  78. [86]

    2022 , journal =

    Enhanced Sampling Methods for Molecular Dynamics Simulations , author =. 2022 , journal =. doi:10.33011/livecoms.4.1.1583 , language =

    2022 doi

  79. [87]

    2009 , journal =

    -dynamics free energy simulation methods , author =. 2009 , journal =. doi:10.1002/jcc.21295 , language =

    2009 doi

  80. [88]

    2023 , journal =

    Fast free energy estimates from -dynamics with bias-updated Gibbs sampling , author =. 2023 , journal =. doi:10.1038/s41467-023-44208-9 , language =

    2023 doi

  81. [89]

    1994 , journal =

    Avoiding singularities and numerical instabilities in free energy calculations based on molecular simulations , author =. 1994 , journal =. doi:10.1016/0009-2614(94)00397-1 , language =

    1994 doi

  82. [90]

    2023 , journal =

    AmberTools , author =. 2023 , journal =. doi:10.1021/acs.jcim.3c01153 , language =

    2023 doi

  83. [91]

    2019 , journal =

    Unified Efficient Thermostat Scheme for the Canonical Ensemble with Holonomic or Isokinetic Constraints via Molecular Dynamics , author =. 2019 , journal =. doi:10.1021/acs.jpca.9b02771 , language =

    2019 doi

  84. [92]

    2026 , journal =

    Convergence is not correctness: context-dependent performance of enhanced-sampling methods across biological complexity , author =. 2026 , journal =. doi:10.1038/s41467-026-74165-y , language =

    2026 doi

  85. [93]

    2012 , journal =

    Alchemical Free Energy Differences in Flexible Molecules from Thermodynamic Integration or Free Energy Perturbation Combined with Driven Adiabatic Dynamics , author =. 2012 , journal =. doi:10.1021/ct300090z , language =

    2012 doi

  86. [94]

    2026 , journal =

    Protonation states of class-C -lactamases determined by a combination of thermodynamic integration and unified free energy dynamics , author =. 2026 , journal =. doi:10.1007/s12039-025-02450-2 , language =

    2026 doi

  87. [95]

    2024 , journal =

    A Comprehensive Study of Factors Affecting the Prediction of the pKa Shift of Asp26 in Thioredoxin Protein , author =. 2024 , journal =. doi:10.1021/acs.jpcb.4c01516 , language =

    2024 doi

Pith tools

Reviewed August 16, 2026 · model on record in the stance chip above.