# Efficient sampling of large-scale transition pathways and intermediate conformations in sub-mesoscopic protein complexes

**Authors:** Domenico Scaramozzino, Byung Ho Lee, Laura Orellana

PMC · DOI: 10.1038/s41467-026-69809-y · 2026-03-02

## TL;DR

The paper introduces eBDIMS2, a new algorithm for simulating protein conformational changes in large complexes, enabling efficient exploration of transition pathways.

## Contribution

eBDIMS2 is a quasi-linear algorithm for simulating large-scale protein transitions, enabling desktop-level simulation of megadalton complexes.

## Key findings

- eBDIMS2 enables simulation of complex transitions in megadalton protein assemblies like ATP synthases on a desktop.
- The algorithm's pathways align with experimental intermediates and enhanced molecular dynamics simulations.
- eBDIMS2 integrates elastic networks with Brownian dynamics for sub-mesoscopic systems.

## Abstract

Protein conformational changes are the cornerstone of biological function. While conformers captured experimentally represent metastable states, the pathways connecting them have been elusive for experiments and simulations alike. Nowadays, cryogenic Electron Microscopy is providing rich structural data on proteins trapped in different states for increasingly large systems, but these are out of scope for most computational methods which exhibit an N2 dependence on size. Based on our previous eBDIMS algorithm, here we present eBDIMS2, an optimized version with quasi-linear size dependence, able to simulate on a desktop computer particularly complex transitions for megadalton protein assemblies, like the rotary motion of ATP synthases. Not only eBDIMS2 pathways spontaneously visit experimental intermediates but also overlap with enhanced and microsecond Molecular Dynamics simulations requiring extensive supercomputing resources. By integrating Elastic Networks with Brownian Dynamics, eBDIMS2 allows an unprecedented exploration of conformational changes of sub-mesoscopic systems previously inaccessible.

Protein conformational changes underlie many biological functions but remain challenging to model, particularly in large complexes. Here authors introduce eBDIMS2, an improved coarse-grained algorithm for simulating these complex transition pathways in large assemblies

## Full-text entities

- **Genes:** RYR2 (ryanodine receptor 2) [NCBI Gene 6262] {aka ARVC2, ARVD2, RYR-2, RyR, VACRDS, VTSIP}, RBP4 (retinol binding protein 4) [NCBI Gene 5950] {aka MCOPCB10, RDCCAS}, PTPA (protein phosphatase 2 phosphatase activator) [NCBI Gene 5524] {aka PARK25, PP2A, PPP2R4, PR53}, CBX4 (chromobox 4) [NCBI Gene 8535] {aka NBP16, PC2}, NF1 (neurofibromin 1) [NCBI Gene 4763] {aka NFNS, VRNF, WSS}, FAS (Fas cell surface death receptor) [NCBI Gene 355] {aka ALPS1A, APO-1, APT1, CD95, FAS1, FASTM}, PRKDC (protein kinase, DNA-activated, catalytic subunit) [NCBI Gene 5591] {aka DNA-PKC, DNA-PKcs, DNAPK, DNAPKc, DNPK1, HYRC}, MARVELD2 (MARVEL domain containing 2) [NCBI Gene 153562] {aka DFNB49, MARVD2, MRVLDC2, Tric}, MPEG1 (macrophage expressed 1) [NCBI Gene 219972] {aka IMD77, MPG1, MPS-1, MPS1, Mpg-1, P-2}, LRP2 (LDL receptor related protein 2) [NCBI Gene 4036] {aka DBS, GP330, LRP-2}, PAH (phenylalanine hydroxylase) [NCBI Gene 5053] {aka PH, PKU, PKU1}, ITPR3 (inositol 1,4,5-trisphosphate receptor type 3) [NCBI Gene 3710] {aka CMT1J, IMD132, IMD133, IP3R, IP3R-3, IP3R3}, HSPD1 (heat shock protein family D (Hsp60) member 1) [NCBI Gene 3329] {aka CPN60, GROEL, HLD4, HSP-60, HSP60, HSP65}, S (surface glycoprotein) [NCBI Gene 43740568] {aka spike glycoprotein}, RYR1 (ryanodine receptor 1) [NCBI Gene 6261] {aka CCO, CMYO1A, CMYO1B, CMYP1A, CMYP1B, KDS}, DROSHA (drosha ribonuclease III) [NCBI Gene 29102] {aka ETOHI2, HSA242976, RANSE3L, RN3, RNASE3L, RNASEN}, ACE2 (angiotensin converting enzyme 2) [NCBI Gene 59272] {aka ACEH}, ACLY (ATP citrate lyase) [NCBI Gene 47] {aka ACL, ATPCL, CLATP}
- **Diseases:** tuberculosis (MESH:D014376), BD (MESH:D000092242), cancer (MESH:D009369), skeletal muscle disorders (MESH:D005207)
- **Chemicals:** F0 (-), Na+ (MESH:D012964), amino acid (MESH:D000596), Cl- (MESH:D002713), ATP (MESH:D000255), AMPPCP (MESH:C005147), salt (MESH:D012492), E2 (MESH:D004958), N2 (MESH:D009584), water (MESH:D014867), chloride (MESH:D002712)
- **Species:** Escherichia coli (E. coli, species) [taxon 562], Mus musculus (house mouse, species) [taxon 10090], Severe acute respiratory syndrome coronavirus 2 (no rank) [taxon 2697049], Severe acute respiratory syndrome-related coronavirus (no rank) [taxon 694009]
- **Mutations:** N234A, N165A

## Figures

7 figures with captions in the complete paper: https://tomesphere.com/paper/PMC12960823/full.md

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