REVIEW 3 major objections 2 minor
Improving sampling of binding free energy differences between covalently bound ligands in alternate binding pockets using MT-REXEE
T0 review · 3 major / 2 minor · reviewed 2026-08-05 · deepseek-v4-flash
Pith's one-line read A multiple-topology replica-exchange method uses ligand chain length as a route between alternate binding pockets, yielding faster convergence of relative binding free energies.
desk verdict Plausible extension of MT-REXEE with a useful case study, but the abstract can't support the unbiased-convergence claim; peer review should demand a detailed-balance check. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
MT-REXEE (multiple topology replica exchange of expanded ensemble): a replica-exchange scheme in which each replica's Hamiltonian is an expanded ensemble over the ligand's alchemical coupling parameter, here the acyl chain length, and configurations are exchanged between replicas at different chain lengths. The key work it does is converting the physically slow transition between pockets into a fast reversible move along the alchemical dimension, so that growing or shrinking the chain visits the chain lengths where the ligand spontaneously switches pockets.
What would settle it
Run a long unbiased simulation of the FabB/ACP complex and measure the occupancy of the two binding pockets as a function of acyl chain length; if there is no chain length at which the ligand substantially occupies both pockets, the MT-REXEE swapping mechanism cannot work. Alternatively, apply MT-REXEE to a pair of binding pockets whose separation is insensitive to chain length and observe no convergence improvement.
Extended reading notes
Core claim
The central claim is that MT-REXEE can overcome high free-energy barriers for flexible ligand–protein systems by using the alchemical progress parameter itself, the ligand's chain length, as a route between otherwise disconnected configurational states. In the FabB/ACP complex, the acyl chain grows while covalently attached to ACP, and at intermediate lengths the chain has nonzero probability of being found in either of two FabB pockets. By running replicas at different chain lengths and exchanging configurations among them, each simulation can access both pockets even when the physical barrier between them is high. This removes the requirement, inherited by collective-variable methods, of s
Load-bearing premise
The enhanced sampling only works if, at some intermediate chain lengths, the ligand has a non-negligible probability of being found in both pockets, so that transitions between the pockets can be visited; if chain length does not open a route between the pockets, the speedup disappears.
Editorial extensions
If this is right
- For flexible complexes with high free-energy barriers, MT-REXEE can recover unbiased sampling of alternate binding modes without user-specified collective variables.
- Relative binding free energies between kinetically separated pockets can be computed with faster convergence, reducing simulation time.
- The open-source implementation makes the approach straightforward to apply to other ligand–protein systems with similar hidden barriers.
- The method extends alchemical free-energy machinery to covalently attached ligands whose chain length can be grown or shrunk, a setting where conventional alchemical mutations are awkward.
Reading between the lines
- If ligand chain length is a sufficient reaction route between any pair of binding sites whose preference changes with ligand size, MT-REXEE could serve as a general proxy for an unknown collective variable, not just for this system.
- A testable diagnostic follows: for a ligand whose two binding pockets are not connected by any physically accessible chain length, the method will not accelerate transitions; this could be checked by mutating the pocket or switching the ligand chemistry.
- The alchemical dimension might also be coupled to protein conformational changes, potentially sampling coupled pocket-opening events, though the paper only demonstrates ligand-pose exchange.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript describes an application and implementation of MT-REXEE (multiple-topology replica exchange of expanded ensemble) to compute relative binding free energies between two alternative binding pockets of the FabB/ACP complex. The central claim is that a new swapping approach within MT-REXEE yields significantly enhanced sampling, demonstrated by faster convergence of free energy estimates for kinetically separated binding pockets. The abstract further claims that the method achieves unbiased sampling of alternate configurational states without prior pocket definitions by letting each simulation visit chain lengths where transitions between pockets occur. The abstract contains only a qualitative statement of enhanced sampling; no quantitative results, error bars, or convergence criteria are provided.
Significance. If the claimed behavior is correct, this is a valuable methodological advance: the approach would address a recognized bottleneck in alchemical free energy calculations for flexible complexes, avoid hand-crafted collective variables, and provide an open-source implementation. These are real strengths. However, as presented, the central claim rests on a speed metric only, with no demonstration that the accelerated estimates are unbiased. The significance of the contribution cannot be assessed from the abstract alone; a quantitative comparison with an independent reference, or at least a statement of a rigorous correctness check, is required.
major comments (3)
- [Abstract] The core claim, 'significantly enhanced sampling ... demonstrated by faster convergence of free energy estimates,' is stated without quantitative support. No numerical free energy differences, error bars, replicate analyses, or convergence thresholds are reported. A reader cannot distinguish fast convergence to the correct value from fast convergence to a biased value. The abstract should report a specific convergence metric (e.g., time to a given uncertainty, agreement with an independent calculation, or variance across repeats) for both the new and baseline approaches.
- [Abstract] Faster convergence does not by itself establish unbiasedness. The novel swapping move must be shown to preserve the target expanded-ensemble distribution, for example by proving detailed balance or by comparing the final free energy estimates to an independent reference calculation. If the acceptance criterion is approximate, or if the adaptive weight updates are not equilibrated, the observed speed-up could be an artifact of systematic bias. The abstract currently offers no statement addressing this load-bearing correctness concern.
- [Abstract] The mechanistic premise that 'allowing each simulation to visit chain lengths where transitions between the pockets occur' is presented as the reason for the enhanced sampling. This is an assumption about the system's free-energy landscape; it should be supported by simulation data, e.g., observed pocket-occupancy transitions as a function of chain length. Without such evidence, the explanation for the claimed improvement is not established.
minor comments (2)
- [Abstract] The phrase 'unbiased sampling' is a strong claim. Consider qualifying it as 'unbiased in the tested cases' or adding a reference to a rigorous test, since no statistical or formal evidence is visible in the abstract.
- [Abstract] Some terms, such as 'MT-REXEE' and 'the new swapping approach,' are defined only by reference to prior work. A brief operational definition in the abstract would help readers assess the novelty and scope of the method.
Circularity Check
No circularity detected in the abstract-only text.
full rationale
The review is limited to the abstract; no equations, derivations, or cited prior work are available to inspect. The central claim is that the new swapping approach in MT-REXEE yields faster convergence of free energy estimates for relative binding affinities in kinetically separated pockets. This is an empirical performance claim, not a derivation that reduces to its inputs. No parameter is fitted and then renamed as a prediction, no target quantity is defined in terms of the method's own outputs, and no load-bearing self-citation appears in the abstract. The abstract's assertion of 'unbiased sampling' is not substantiated with details, but lack of proof is a correctness/validation concern, not circularity. Therefore the appropriate finding is no significant circularity.
Assumptions & free parameters
assumptions (4)
- domain assumption Alchemical free energy estimators are unbiased when sampling is ergodic
- domain assumption The molecular mechanics force field used for FabB/ACP represents the relevant binding free energies accurately enough for the comparison
- domain assumption Growing/shrinking acyl chain length connects the two binding pockets and enables transitions
- domain assumption Previously published experimental evidence of two binding pockets is correct
Cite this review
Pith. "Pith review of Improving sampling of binding free energy differences between covalently bound ligands in alternate binding pockets using MT-REXEE." pith.science (2026). https://pith.science/paper/76THIISJ
@misc{pith2026250806720,
author = {Pith},
title = {Pith review of: Improving sampling of binding free energy differences between covalently bound ligands in alternate binding pockets using MT-REXEE},
year = {2026},
howpublished = {\url{https://pith.science/paper/76THIISJ}},
note = {Machine review of arXiv:2508.06720}
}
read the original abstract
The primary limitation for the application of alchemical free energy methods to a wider variety of complex molecular systems is achieving reasonable sampling. Flexible binding complexes often have high free energy barriers, which require prohibitively long simulations or carefully tuned enhanced sampling methods in order to gather sufficient uncorrelated samples to obtain reliable free energy estimates. An example of such a flexible system is the complex formed between FabB, an elongating \b{eta}-ketoacyl-acyl carrier protein (ACP) synthase (KS) from Escherichia coli, and ACP, which carries acyl chains of varying lengths. Previous experimental evidence suggests that growing acyl chains can bind to at least two pockets in FabB. With the multiple topology replica exchange of expanded ensemble (MT-REXEE) enhanced sampling approach, we can obtain highly efficient sampling of both pockets by adaptively growing and shrinking the chains in the simulation ensemble, allowing each simulation to visit chain lengths where transitions between the pockets occur. This enables unbiased sampling of alternate configurational states for large complex systems without prior pocket definitions, as collective-variable based enhanced sampling methods would require. Using the new swapping approach gives significantly enhanced sampling even for this simpler problem, as demonstrated by faster convergence of free energy estimates of relative binding affinity between kinetically separated binding pockets. This case study demonstrates the utility of MT-REXEE and its open-source implementation for systems that feature high free energy barriers for a subset of ligands of interest, demonstrating a valuable addition to the existing stable of enhanced sampling methods.
Reviewed August 5, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.