REVIEW 4 major objections 5 minor 1 cited by
BioMD: All-atom Generative Model for Biomolecular Dynamics Simulation
T0 review · 4 major / 5 minor · reviewed 2026-08-05 · deepseek-v4-flash
Pith's one-line read BioMD claims to be the first all-atom generative model that simulates long-timescale protein-ligand dynamics, generating physically plausible all-atom trajectories and complete ligand unbinding pathways for up to 97.1% of protein-ligand sys
desk verdict A well-engineered generative model for all-atom trajectories whose headline 97.1% unbinding rate is likely inflated by a weak convex-hull criterion and AR drift. 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
The central object is a unified trajectory flow-matching model with a 'noising-as-masking' schedule: each frame is independently noised according to a time variable, so known conditioning frames are kept clean and frames to be generated are initialized from noise, with different masking schedules realizing coarse-grained forecasting or fine-grained interpolation within the same architecture. This reduces sequence length by decoupling long-term evolution from local dynamics and helps manage error accumulation for long trajectories. An SE(3)-equivariant graph transformer encodes the initial conformation as conditional embeddings, and the FlowTrajectoryTransformer uses AttentionPairBias for int
What would settle it
Re-run molecular dynamics or enhanced sampling seeded from BioMD's predicted exit structures and check whether they are metastable dissociated states; if the predicted exits collapse back to the binding site or sit in high free-energy barriers, the claimed unbinding paths fail. A faster check: apply the same convex-hull success metric to random-walk trajectories and compare the resulting Success@10 rate.
Extended reading notes
Core claim
BioMD's core discovery claim is that long-timescale biomolecular motion can be generated all-atom by one conditional flow-matching velocity network, with the distinction between forecasting and interpolation encoded only as different noise/masking schedules. The model takes the first frame as a clean condition, denoises coarse frames spaced every k steps, then refills the gaps from two clean anchors; at inference the auto-regressive mode reconditions on its own history. The paper reports that this scheme yields complete ligand unbinding trajectories for 92.9% of systems within five attempts and 97.1% within ten, while also reproducing two known metadynamics pathways on 6EY8 and discovering a
Load-bearing premise
The load-bearing premise is that a ligand counted as 'unbound' because its center moved outside the protein's starting shape has really dissociated, not just drifted or distorted as the model's own error-accumulation numbers suggest.
Editorial extensions
If this is right
- Long-timescale ligand dissociation, normally requiring enhanced sampling, becomes a direct generative sampling problem.
- The same architecture handles both coordinate accuracy (BioMD-abs) and exploratory sampling (BioMD-rel) by switching prediction target or masking schedule.
- Generated trajectories retain local chemical plausibility, with bond and angle errors below thermal fluctuation thresholds, even with autoregressive error accumulation.
- The method can run in seconds on a single GPU, making pathway exploration practical for drug-discovery screening.
- Protein flexibility, not just ligand motion, is captured when many baselines treat the receptor as static.
Reading between the lines
- If convex-hull exits are not validated energetically, the reported success rates may partly reflect accumulated drift; a fair test would compare free-energy profiles or committor probabilities of generated exits against reference dissociation paths.
- The forecasting-then-interpolation design suggests a general recipe for other rare-event trajectory tasks, such as folding or conformational transitions, where coarse milestones are known but fine dynamics are not.
- BioMD-abs versus BioMD-rel trade-off implies a Pareto frontier: the same model family can be pushed toward accurate reproduction or toward exploration, and a tuned interpolation may combine both.
- Because the model is trained on simulated trajectories, its 'discovered' novel pathways are only as good as the training distribution's coverage; genuinely new pathways would need experimental or enhanced-sampling confirmation.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper introduces BioMD, an all-atom generative model for protein-ligand dynamics, built on a unified conditional flow-matching framework with hierarchical coarse-grained forecasting and fine-grained interpolation. The model is evaluated on two datasets: MISATO for ligand dynamics in the binding pocket and DD-13M for ligand unbinding. The authors report high physical stability and competitive RMSF correlations on MISATO, and on DD-13M report that the auto-regressive variant BioMD-rel (AR-5) achieves 70.9%/92.9%/97.1% unbinding success at @1/@5/@10 attempts, with the claim that BioMD is the first all-atom generative model to simulate long-timescale ligand-unbinding pathways. The central scientific claim is the 97.1% unbinding success and the physical meaningfulness of the generated unbinding paths.
Significance. If the central claim is correct, the work would be a notable advance: a generative model that produces plausible all-atom trajectories, including ligand dissociation events, at orders-of-magnitude lower cost than metadynamics. The methodological core is sensible: using a single flow-matching model with different masking schedules for forecasting and interpolation is elegant, and the architecture is described in sufficient detail to be reproducible. The qualitative 6EY8 example, in which BioMD reproduces two known metadynamics pathways and finds a third, is compelling but anecdotal. The paper also honestly acknowledges the error-accumulation problem in Section 5.3. The main weakness is that the headline unbinding metric is not physically validated, and the paper's own path-RMSD and clash numbers are consistent with the concern that many 'successful' unbinding events are artifacts of cumulative AR drift or distorted geometries rather than genuine dissociation transitions. The result is potentially important, but the current empirical support is insufficient to sustain the strong conclusions.
major comments (4)
- [Sec. A.3.4 / Table 2] The unbinding success metric—defined as the ligand centroid leaving the convex hull of the initial protein heavy atoms—is necessary but not sufficient for a physical dissociation pathway. The paper's own numbers in Table 2 are consistent with the concern that convex-hull exits are produced by cumulative drift: BioMD-rel (AR-5), which has the highest success rate, has an unbinding path RMSD (0.7055 Å) worse than the Static baseline (0.6504 Å) and an protein-ligand clash score (0.6375) far above Static (0). Section 5.3 documents AR error accumulation. The authors should either replace this metric with a stricter one (e.g., sustained outward displacement of the ligand from the binding site, requiring the exit to be irreversible and not solely occurring in the last frame) or provide additional physical validation: contact maps at the exit event, ligand-pocket distance profiles, free-energy o
- [Table 2 / Sec. A.3.3] The Unbinding Path RMSD metric needs clarification and appears internally inconsistent. The DD-13M dataset consists of dissociation trajectories, so a static ligand centroid should be very far from a reference unbinding path. Yet the Static baseline reports an RMSD of only 0.6504 Å, comparable to the best BioMD variants. This suggests either that the metric's best-match/resampling procedure is not sensitive to unbinding, or that the reference trajectories used for comparison are not predominantly unbinding-like, or that an alignment step removes the signal. The authors should specify the alignment, how trajectories are resampled, whether only successful reference dissociation paths are used, and how the best-match search is performed. As reported, Table 2's path-accuracy column does not provide meaningful evidence for the physical fidelity of the generated unbinding paths.
- [Sec. 5.3 / Table 2] The auto-regressive variant that achieves the headline unbinding success also shows substantially degraded physical stability: bond MAE increases from 0.0308 (BioMD-rel) to 0.0580 (BioMD-rel AR-5), angle MAE from 0.0606 to 0.0918, and protein-ligand clashes from 0.0004 to 0.6375 (Table 2). Section 5.3 acknowledges this but argues the errors are correctable by 'a simple local refinement step' without demonstrating that such refinement preserves the unbinding pathways. If the reported 97.1% success is computed on unrefined trajectories where the ligand may be sterically clashing with the protein, the physical interpretation of 'complete unbinding paths' is questionable. The authors should report results after the proposed refinement, or at minimum show that the geometric distortions do not drive the hull-exit events.
- [Sec. A.1 / Algorithms 3-4] The Euler solver uses dt=0.1 with 10 integration steps (Algorithm 3: τ in {0, 0.1, ..., 0.9}), and the paper does not provide a sensitivity analysis of success/path metrics to the integrator step size. Since the headline results depend on long autoregressive generation, the numerical integration scheme is a potentially important source of error. A convergence check or a comparison with a smaller dt would strengthen the claim that the generated trajectories are stable and not numerical artifacts.
minor comments (5)
- [Abstract / Sec. 1] The claim 'first all-atom generative model to simulate long-timescale protein-ligand dynamics' is strong. Given that MDGen and other trajectory models exist, the 'first' claim should be qualified with respect to protein-ligand systems specifically and justified in the related-work section.
- [Sec. 4.2.1] Typo: 'nmasked' should be 'unmasked'.
- [Algorithm 2] Typo: 'Indepentent' should be 'Independent'.
- [A.3.3] The 'best-match search' over reference trajectories needs a precise statement of whether the same reference path can be selected for multiple generated trajectories and how the reported mean is aggregated over complexes; otherwise the RMSD values are hard to interpret.
- [Sec. 5.1 / Table 1] The RMSF correlation for BioMD-abs on protein atoms (0.685) is good, but the authors do not report confidence intervals or significance tests against NeuralMD; adding these would strengthen the comparison.
Circularity Check
No circular derivation chain found; minor self-citation via DD-13M dataset and training-aligned geometric metrics, but core predictions are on held-out data.
full rationale
BioMD's central derivation is a supervised generative model trained on trajectory data; evaluations on held-out complexes (DD-13M) and held-out systems (MISATO) are genuine predictions rather than fits. The 'unbinding success' criterion (convex-hull exit) is a self-defined, permissive proxy, but it is not an input to training and the model is not fitted to it, so it is a validity concern rather than a circular reduction. The DD-13M dataset originates from the same research group (ref [18]), which constitutes a self-citation, but it provides a fixed held-out test set and the model is also validated on the external MISATO dataset. The reported bond-length and steric-clash stability metrics are partially aligned with the auxiliary training losses (Lgeom in Appendix A.2), yet they are measured on test data and compared against baselines, so low values are not forced by construction. No step in the claimed derivation reduces to its own input or to a self-citation chain; therefore no significant circularity is identified.
Assumptions & free parameters
free parameters (6)
- coarse-grain interval k =
10
- AR block size j =
5
- data noise scale sigma_data =
16
- collision thresholds zeta_pl, zeta_ll =
3.0 Å, 2.0 Å
- auxiliary loss weights lambda_col, lambda_bond, lambda_ctr =
not reported
- Euler solver step dt =
0.1
assumptions (6)
- standard math Flow matching regresses a conditional velocity field that is a valid generative model for trajectory frames
- domain assumption Independent per-frame noise ('noising-as-masking') supports arbitrary conditioning patterns, so one model can do forecasting and interpolation with different mask schedules
- domain assumption Short-timescale conformational changes are small, so a trajectory can be decomposed into coarse forecasting at intervals of k=10 frames plus interpolation without losing the true dynamics
- domain assumption The reference trajectories (MISATO MD and DD-13M metadynamics) are an adequate ground truth for the dynamics to be learned and validated
- domain assumption The auxiliary geometric losses do not push the learned distribution away from the Boltzmann distribution
- domain assumption The 36 held-out DD-13M complexes and the 1,031 MISATO targets (sequence length <= 800) are representative and independent of the training set
Cite this review
Pith. "Pith review of BioMD: All-atom Generative Model for Biomolecular Dynamics Simulation." pith.science (2026). https://pith.science/paper/EUSOEM5I
@misc{pith2026250902642,
author = {Pith},
title = {Pith review of: BioMD: All-atom Generative Model for Biomolecular Dynamics Simulation},
year = {2026},
howpublished = {\url{https://pith.science/paper/EUSOEM5I}},
note = {Machine review of arXiv:2509.02642}
}
read the original abstract
Molecular dynamics (MD) simulations are essential tools in computational chemistry and drug discovery, offering crucial insights into dynamic molecular behavior. However, their utility is significantly limited by substantial computational costs, which severely restrict accessible timescales for many biologically relevant processes. Despite the encouraging performance of existing machine learning (ML) methods, they struggle to generate extended biomolecular system trajectories, primarily due to the lack of MD datasets and the large computational demands of modeling long historical trajectories. Here, we introduce BioMD, the first all-atom generative model to simulate long-timescale protein-ligand dynamics using a hierarchical framework of forecasting and interpolation. We demonstrate the effectiveness and versatility of BioMD on the DD-13M (ligand unbinding) and MISATO datasets. For both datasets, BioMD generates highly realistic conformations, showing high physical plausibility and low reconstruction errors. Besides, BioMD successfully generates ligand unbinding paths for 97.1% of the protein-ligand systems within ten attempts, demonstrating its ability to explore critical unbinding pathways. Collectively, these results establish BioMD as a tool for simulating complex biomolecular processes, offering broad applicability for computational chemistry and drug discovery.
Forward citations
Cited by 1 Pith paper
-
Spectral Diffusion for Protein Dynamics
Diffusion over DCT spectral volumes of Cα displacements yields fast, temperature-conditioned protein trajectories with RMSF Pearson r of 0.844 on held-out mdCATH.
Reviewed August 5, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.