# A multi-grained symmetric differential equation model for learning protein-ligand binding dynamics

**Authors:** Shengchao Liu, Weitao Du, Hannan Xu, Yanjing Li, Zhuoxinran Li, Vignesh Bhethanabotla, Divin Yan, Christian Borgs, Anima Anandkumar, Hongyu Guo, Jennifer Chayes

PMC · DOI: 10.1038/s41467-025-67808-z · Nature Communications · 2025-12-30

## TL;DR

NeuralMD is a new AI framework that improves the accuracy and efficiency of simulating how proteins and ligands bind, which could speed up drug discovery.

## Contribution

NeuralMD introduces a physics-informed, multi-grained, group-symmetric framework for modeling protein-ligand binding dynamics.

## Key findings

- NeuralMD achieves up to 15× lower reconstruction error compared to existing ML baselines.
- The model shows 70% higher validity in predicting protein-ligand binding dynamics.
- Predicted oscillations closely align with ground-truth dynamics in simulations.

## Abstract

Molecular dynamics (MD) simulation is a key tool in drug discovery for predicting protein-ligand binding affinities, transport properties, and pocket dynamics. While advances in numerical and machine learning (ML) methods have improved MD efficiency, accurately modeling long-timescale dynamics remains challenging. We introduce NeuralMD, an ML surrogate that accelerates and enhances MD simulations of protein-ligand binding. NeuralMD employs a physics-informed, multi-grained, group-symmetric framework comprising (1) BindingNet, which enforces symmetry via vector frames and captures multi-level protein-ligand interactions, and (2) an augmented neural differential equation solver that learns trajectories under Newtonian mechanics. Across ten single-trajectory and three multi-trajectory tasks, NeuralMD achieves up to 15 × lower reconstruction error and 70% higher validity than existing ML baselines. The predicted oscillations closely align with ground-truth dynamics, establishing NeuralMD as a foundation for next-generation protein-ligand simulation research.

Artificial intelligence is advancing molecular simulation. Here, the authors introduce NeuralMD, a physics-informed AI framework that efficiently models and predicts protein-ligand binding dynamics, providing opportunities for accelerating drug discovery.

## Full-text entities

- **Genes:** TTC41P (tetratricopeptide repeat domain 41, pseudogene) [NCBI Gene 253724] {aka GNN, GNNP}, NFKB2 (nuclear factor kappa B subunit 2) [NCBI Gene 4791] {aka CVID10, H2TF1, LYT-10, LYT10, NF-kB2, p100}, CHUK (component of inhibitor of nuclear factor kappa B kinase complex) [NCBI Gene 1147] {aka BPS2, IKBKA, IKK-1, IKK-alpha, IKK1, IKKA}, MAP3K14 (mitogen-activated protein kinase kinase kinase 14) [NCBI Gene 9020] {aka FTDCR1B, HS, HSNIK, IMD112, NIK}, RELB (RELB proto-oncogene, NF-kB subunit) [NCBI Gene 5971] {aka I-REL, IMD53, IREL, REL-B}, NFKB1 (nuclear factor kappa B subunit 1) [NCBI Gene 4790] {aka CVID12, EBP-1, KBF1, NF-kB, NF-kB1, NF-kappa-B1}, Map3k14 (mitogen-activated protein kinase kinase kinase 14) [NCBI Gene 53859] {aka Nik, aly}
- **Diseases:** type 2 diabetes (MESH:D003924), solid cancers (MESH:D009369), GNN-MD (MESH:D000092242), autoimmune diseases (MESH:D001327), obesity (MESH:D009765), ML (MESH:D007859), inflammation (MESH:D007249), hematologic malignancies (MESH:D019337), cardiovascular disease (MESH:D002318)
- **Chemicals:** 6-alkynylindoline (-), amino acid (MESH:D000596), water (MESH:D014867), polymer (MESH:D011108), carbon (MESH:D002244)
- **Species:** Homo sapiens (human, species) [taxon 9606], Mus musculus (house mouse, species) [taxon 10090]

## Full text

_Full body text omitted from this summary view._ Fetch the complete paper as Markdown: https://tomesphere.com/paper/PMC12848303/full.md

## Figures

3 figures with captions in the complete paper: https://tomesphere.com/paper/PMC12848303/full.md

## References

15 references — full list in the complete paper: https://tomesphere.com/paper/PMC12848303/full.md

---
Source: https://tomesphere.com/paper/PMC12848303