Pith. sign in

REVIEW 3 major objections 6 minor 3 references

Enhanced Sampling of Protein Conformational Changes via True Reaction Coordinates from Energy Relaxation

T0 review · 3 major / 6 minor · reviewed 2026-08-11 · deepseek-v4-flash

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

Pith's one-line read Relaxation-derived coordinates speed protein transitions by 10^5- to 10^15-fold

desk verdict A genuinely new route to reaction coordinates from cheap relaxation runs, with one solid implicit-solvent check and an unverified leap to explicit solvent. read the letter →

arxiv 2412.04400 v1 pith:ZKTSKFZU submitted 2024-12-05 physics.chem-ph

Huiyu Li , Ao Ma This is my paper · ORCID
classification physics.chem-ph
keywords truereactioncoordinatesenhancedsamplinggeneralizedworkfunctionalenergyrelaxationHIV-1proteasePDZdomainallosterytransitionpathcollectivevariables
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 claims that the coordinates along which a protein dissipates excess energy are the same coordinates it must activate to change shape, so the true reaction coordinates (the few variables that predict which side of a barrier a trajectory will end on) can be obtained from a 5 ps energy-relaxation simulation of a single structure. If correct, this dissolves the paradox that optimal collective variables could only be found from unbiased reactive trajectories, which are exactly what good collective variables are needed to produce. The authors report that biasing these relaxation-derived coordinates opens HIV-1 protease flaps and dissociates DRV and MA/CA ligands within 200 ps, a claimed acceleration of $10^{5}$ to $10^{15}$ over unbiased simulation and experiment, and that the resulting trajectories pass through transition states and can seed natural reactive trajectories via transition path sampling. In a blind test on PDZ2, the same protocol predicts a transient cleft-opening mechanism for a decades-old allostery puzzle.

What carries the argument

The generalized work functional (GWF) is a tensor-valued generalization of mechanical work whose elements are products of force and displacement components; its singular value decomposition produces an orthonormal set of singular coordinates ranked by the potential energy flowing through each coordinate. The leading singular coordinates, carrying the largest potential energy flows during energy relaxation, are the paper's claimed true reaction coordinates. The physical bridge is Onsager's regression hypothesis: the same channels that dissipate deposited energy in a 5 ps relaxation are the channels that must absorb energy to cross a barrier during activation.

What would settle it

Apply the same protocol to a protein whose true reaction coordinates are already known from unbiased transition path sampling: if the leading singular coordinates from energy relaxation have inner products far from 1 with the known tRCs, or if bias along them fails to produce trajectories that pass TPS shooting, the claimed identity between relaxation and activation coordinates is refuted.

Watch

Extended reading notes

Core claim

The central discovery is that true reaction coordinates control both activation and energy relaxation. Using the generalized work functional, the authors compute singular coordinates ranked by potential energy flow through each coordinate during relaxation; the leading singular coordinates from a 5 ps ensemble match, to near unit inner product, the tRCs previously identified from natural reactive trajectories for flap opening in HIV-PR in implicit solvent. In explicit solvent, biases applied to these leading singular coordinates open the flaps and drive complete DRV and MA/CA unbinding in 200 ps, with the experimental DRV unbinding half-life as reference, while shooting from conformations on these trajectories produces unbiased natural reactive trajectories. The same procedure applied to PDZ2 yields natural trajectories showing transient opening of the alpha2-beta2 cleft and binding groove during ligand unbinding, which the authors propose as the physical basis of PDZ allostery.

Load-bearing premise

The load-bearing premise is that the channels a protein uses to shed a burst of excess energy in a few picoseconds are the same channels it must energize to cross a barrier in a rare conformational change, so ranking coordinates by relaxation energy flow reveals the coordinates that control the transition.

Editorial extensions

If this is right

  • Biasing relaxation-derived tRCs should reproducibly open HIV-PR flaps and dissociate DRV and MA/CA ligands within about 200 ps, a claimed 10^5- to 10^15-fold acceleration over unbiased MD and experiment.
  • RC-uncovered trajectories should provide transition-state conformations that make transition path sampling practical, yielding the first natural reactive trajectories for ligand dissociation from HIV-PR.
  • Because tRCs come from a single structure, the method should predict functional conformational changes and ligand release for proteins with no prior reactive trajectory data.
  • For PDZ2, the predicted tRCs imply that allosteric effectors such as Cdc42 and the alpha_A helix alter ligand affinity by sterically interfering with transient cleft and groove opening during unbinding.
  • For multi-step reactions, repeating the relaxation-and-bias protocol basin by basin should build a network of natural reactive trajectories and rate constants.

Reading between the lines

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

  • Not in the paper: if relaxation and activation share tRCs in general, protein-level fluctuation-dissipation reasoning could extend to far-from-equilibrium energy injections, a generalization the authors only gesture at.
  • A reader could test the method against machine-learned slow-mode collective variables on the same proteins; where the two agree, energy-flow and statistical descriptions point to the same physics, and where they disagree the discrepancy would reveal which motions relaxation actually resolves.
  • The method's input is a single conformation, so a natural blind test is a protein with a known cryptic pocket or experimentally characterized allostery; success would mean tRCs can be predicted before any transition is observed.
  • The paper's tRCs are linear combinations of backbone dihedrals; whether linearity is generic remains open, and curved tRCs would require the piecewise linearization the authors note.
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

3 major / 6 minor

Summary. The paper proposes that true reaction coordinates (tRCs) for protein conformational changes can be computed from short (5 ps) energy-relaxation simulations of a single protein structure, using the generalized work functional (GWF) method. The authors validate this for ligand-free HIV-PR in implicit solvent by showing that the leading singular coordinates (SCs) from energy relaxation have inner products close to 1 with previously identified tRCs. They then apply bias potentials along these SCs to accelerate flap opening and ligand dissociation in HIV-PR in explicit solvent, and to propose a mechanism for PDZ2 allostery. The reported accelerations are 10^5- to 10^15-fold, and shooting moves from the biased trajectories are used to argue that the biased trajectories pass through transition states and follow natural transition pathways.

Significance. If the central claim is correct, the method would break the long-standing circularity that tRCs are needed for enhanced sampling but are themselves identified only from natural reactive trajectories, which require enhanced sampling. The implicit-solvent validation in Fig. 2a is a quantitative, direct check and is the strongest evidence in the paper. The GWF formalism is an original theoretical contribution, and the paper makes a falsifiable prediction about PDZ2 allostery that is mechanistically interesting. However, the validation in explicit solvent and for PDZ2 is indirect, the acceleration factors are rough single-trajectory comparisons without error bars, and the PDZ mechanism is not a quantitatively falsifiable prediction. With additional validation the work could be highly impactful, but in its present form the evidence is not commensurate with the strength of the claims.

major comments (3)
  1. [RC-uncovered trajectories follow natural transition pathways (p.16)] The shooting-move validation with only five pairs of trajectories per conformation is statistically underpowered. A single successful reactive pair out of five attempts is consistent with a wide range of committor values (e.g., p_B between roughly 0.1 and 0.9 at the 95% confidence level in a binomial model), so it does not establish that the chosen conformations are near p_B=0.5 or that the SCs are tRCs. The paper does not report the number of conformations attempted, the success rate, or the distribution of p_B. To support the claim that explicit-solvent RC-uncovered trajectories follow natural pathways, the authors should provide a committor histogram or at least a statistically meaningful count of shooting outcomes for at least one explicit-solvent system.
  2. [A hypothesis on activation and energy relaxation (p.9)] The foundational assumption that the leading SCs of energy relaxation are identical to the tRCs for the conformational transition is tested quantitatively only for one system (ligand-free HIV-PR in implicit solvent, Fig. 2a). For DRV-bound and MA/CA-bound HIV-PR in explicit solvent and for PDZ2, the tRC assignment is based on the PEF gap in Fig. 2b and Fig. S5a, but a gap in potential-energy-flow magnitudes does not establish that these coordinates determine the committor. The paper would be substantially strengthened by a committor test for at least one explicit-solvent system, even a small one, or by a systematic sensitivity analysis showing that the resulting SCs are robust to the choice of injection site, cluster selection, and the number of leading SCs.
  3. [Efficiency of enhanced sampling by tRCs (p.18)] The acceleration factors are not rigorously defined. Comparing the duration of a single biased RC-uncovered trajectory (200 ps) with an experimental half-life (8.9×10^? s) conflates a biased trajectory time with an unbiased mean first-passage time, and no error bars or replicate statistics are provided. The 10^15-fold factor appears to be obtained by multiplying a waiting-time reduction with a diffusive-motion reduction, but these are not measured on the same footing and are based on single trajectories (e.g., one NRT for DRV, one NRT for MA/CA). Please provide precise definitions of the acceleration factor, an ensemble average with uncertainty, and a discussion of the different physical meanings of biased trajectory time, NRT duration, and experimental lifetime.
minor comments (6)
  1. [Throughout] Many superscripted numbers and exponents are garbled (e.g., '10! to 10"!-fold', '8.9×10! s', '10>-fold' in Sections 2, 13, and 18), making quantitative claims unreadable. The typesetting must be fixed.
  2. [Fig. 2a (p.10)] The inner products are described only as 'close to 1'; the numerical values and any uncertainty estimates should be reported in the text or a table so that the reader can judge the agreement.
  3. [p.9, before Fig. 1] The stray text 'dfdf' appears immediately before Fig. 1 and should be removed.
  4. [Methods (p.25)] The statement that testing ligand temperature increases of 150K and 800K gives 'the same results' is not substantiated with data. Please include a supporting figure or table in the Supplementary Information.
  5. [Methods (p.26)] Candidate TS conformations for shooting are selected 'based on visual inspection and intuition.' This is not a reproducible criterion; please provide a quantitative selection procedure or at least describe the visual criteria in detail.
  6. [Comparison with empirical CVs (p.18, Fig. 6)] The metadynamics comparison uses Gaussian heights up to 10,000-fold the recommended value, which is an extreme bias regime. The non-physical features observed in Fig. 6e may be aggravated by this regime; the paper should acknowledge this limitation and ideally also show a comparison at a more moderate bias that still achieves dissociation.

Circularity Check

1 steps flagged · score 5.0 of 10

Partial circularity: 'tRCs control energy relaxation' restates the GWF definition of tRCs as maximal-PEF SCs, so computing tRCs from energy relaxation is in part a construction, not a prediction.

  1. self definitional [Results, 'Potential energy flow'/'Generalized work functional'; 'A hypothesis on activation and energy relaxation'; 'Leading SCs of energy relaxation are the same as tRCs...' (Fig. 2a)]
    "Consequently, tRCs are identified as the SCs with the highest PEFs 34,40,43. ... we propose that the coordinates essential for energy relaxation (i.e. its leading SCs) are identical to those governing activation (i.e. the tRCs). ... Figure 2a shows the inner product between the six leading SCs for the energy relaxation of HIV-PR in implicit solvent and the six tRCs of its flap opening we identified in ref. 34. All inner products are close to 1—confirming our hypothesis."

    Once tRCs are operationally defined as the SCs with the highest PEFs, 'tRCs control energy relaxation' is true by construction, and computing 'tRCs' from an ER trajectory is exactly computing its leading PEF SCs via GWF. The inner-product check compares this ER output with the authors' previously GWF-derived tRCs; it does not recompute the committor for the ER-derived coordinates. For explicit-solvent HIV-PR and PDZ2, 'tRC' is assigned from the PEF gap (Fig. 2b, Fig. S5a) without any committor test, so the subsequent claims that biasing 'tRCs' accelerates transitions and follows natural pathways inherit this definitional identification. The independent residue is ref. 34's prior committor validation and the TPS shooting tests, which do not repair the tautological energy-relaxation half.

full rationale

The GWF/SVD construction itself is internally consistent, and the committor-based definition of tRCs is not in question. The circularity enters at the moment the paper equates tRCs with the SCs having the largest PEFs and then presents 'tRCs control energy relaxation' as a discovery. With that equation, the energy-relaxation half of the central claim is a restatement of the method, and the Fig. 2a confirmation is a comparison between two GWF-derived sets rather than a fresh committor test. The explicit-solvent and PDZ applications additionally assign 'tRC' from the PEF gap without direct committor validation. These are genuine circularity/evidentiary concerns. Other weaknesses are evidentiary rather than circular: the shooting validation uses only five pairs per candidate selected 'based on visual inspection and intuition' (Results, 'RC-uncovered trajectories follow natural transition pathways'; Methods), and the reported 10^5-10^15-fold accelerations depend on hand-tuned spring constants ('we empirically selected k=2000 kJ mol^-1 ... k=500 kJ mol^-1'; supplementary videos use k=30,000 kJ/mol). Those parameters make the acceleration numbers contingent, but they are not the core definitional reduction. Score 5 reflects the partial circularity in the ER-to-tRC identification while acknowledging the independent TPS shooting checks and the external experimental time-scale comparisons.

Assumptions & free parameters 6 free parameters · 5 assumptions · 0 invented entities

The central method depends on several hand-picked parameters (bias spring constants, inclusion threshold, number of SCs, cluster choice) and on the untested assumption that relaxation and activation share the same coordinates.

free parameters (6)
  • epsilon (epsilon) for dihedral inclusion = 0.03
    Dihedrals with coefficient magnitude above epsilon are included in the biased collective variable Q_i; chosen by hand to filter noise.
  • number of bins for adaptive bias = 100
    The target range of Q_i is divided into 100 bins over which the bias center moves; controls the bias schedule.
  • spring constant k for flap opening = 2000 kJ/mol
    Empirically selected to balance flap opening duration; see Methods.
  • spring constant k for flap closing = 500 kJ/mol
    Empirically selected for closing direction in reversible transitions.
  • spring constant k for ligand dissociation trajectories = 30000 kJ/mol
    Used for RC-uncovered trajectories in HIV-PR systems (Supplementary videos).
  • number of leading SCs selected as tRCs = 6 (HIV-PR implicit, PDZ2); 7 (DRV-bound HIV-PR)
    Determined by visual gap in potential energy flow spectra, which sets the dimensionality of the biased subspace.
assumptions (5)
  • domain assumption The GWF method identifies tRCs as the singular coordinates with the highest potential energy flows.
    Adopted from the authors' prior work (ref. 34); validated for HIV-PR flap opening in implicit solvent but assumed for other systems.
  • ad hoc to paper The leading singular coordinates of energy relaxation are identical to the tRCs for conformational activation.
    Core hypothesis introduced in 'A hypothesis on activation and energy relaxation'; based on Onsager's regression hypothesis and supported only by the implicit-solvent HIV-PR inner-product comparison.
  • domain assumption The tRCs are linear combinations of backbone dihedrals within numerical error.
    The method uses a fixed linear form Q_i = sum U_ij cos(chi_j - chi_j*); the authors note in the Discussion that linearity is an open question.
  • ad hoc to paper Energy relaxation from a single structure with deposited kinetic energy samples the same essential coordinates as the spontaneous transition.
    This is the basis for predictive sampling and is not independently validated.
  • ad hoc to paper The clustering of energy relaxation trajectories and the choice of cluster E1 for GWF analysis do not affect the resulting tRCs.
    The authors state results for E1 and E2 are identical, but provide no comparison for explicit-solvent systems or PDZ.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Enhanced Sampling of Protein Conformational Changes via True Reaction Coordinates from Energy Relaxation." pith.science (2026). https://pith.science/paper/ZKTSKFZU

@misc{pith2026241204400,
  author       = {Pith},
  title        = {Pith review of: Enhanced Sampling of Protein Conformational Changes via True Reaction Coordinates from Energy Relaxation},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/ZKTSKFZU}},
  note         = {Machine review of arXiv:2412.04400}
}
read the original abstract

The bottleneck in enhanced sampling lies in finding collective variables (CVs) that can effectively accelerate protein conformational changes. True reaction coordinates (tRCs) that can predict the committor are considered the optimal CVs, but identifying them requires unbiased natural reactive trajectories, which, paradoxically, depend on effective enhanced sampling. Using the generalized work functional method, we found that tRCs control both conformational changes and energy relaxation, enabling us to compute tRCs from energy relaxation simulations. Applying bias to tRCs accelerated conformational changes and ligand dissociation in HIV-1 protease and the PDZ2 domain by 10^5 to 10^15-fold. The resulting trajectories follow natural transition pathways, enabling efficient generation of natural reactive trajectories. In contrast, biased trajectories from empirical CVs often display non-physical features. Furthermore, by computing tRCs from a single protein structure, our method enables predictive sampling of conformational changes. These findings significantly broaden the range of protein functional processes accessible to molecular dynamics simulations.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

3 extracted references · 2 canonical work pages

  1. [29]

    Q. At any instant t, when the system configuration is at 𝑄=(𝑡)∈z𝑄=,C,𝑄=,CD

    instead. This is to mimic the process of dissipating the excess energy at the active site after ligand binding or enzymatic reaction. To test if the SCs depend on the specific amount of excess energy deposited into the ligand, we also tried increasing the ligand temperature by 150K and 800K. The results are the same. For the GWF analysis, the 8,000 ER tra...

    work page 2000

  2. [43]

    It remains to be seen whether this linearity is a general feature of tRCs in proteins and if there is a fundamental physical reason for it

    In practice, for all processes we have studied, the tRCs have been linear within numerical error 34,43. It remains to be seen whether this linearity is a general feature of tRCs in proteins and if there is a fundamental physical reason for it. Methods All simulations are constant NVE and use the CHARMM36m force field and TIP3P water model 61,62. For HIV-P...

    work page 2019

  3. [67]

    2∑zΔ𝑊|$(0→𝑛𝛿𝑡;𝛼)Δ𝑊|$(0→𝑛𝛿𝑡;𝛽) 2C4

    Additional Gaussian heights of 1, 100, 1000, and 2000 kJ/mol are also tested (Fig. S4) to achieve MA/CA dissociation. Data Availability All data generated or analyzed during this study are included in this published article and its Supplementary Information file. Code Availability All the custom codes used in this study are deposited in Code Ocean. Author...

    arXiv 1991

Pith tools

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