Search papers, labs, and topics across Lattice.
This paper introduces the Bayes-optimal RMSD (BRMSD), a novel metric for structural comparison in molecular dynamics that optimally weights per-atom contributions based on their mobility. By maximizing a Bayesian posterior, BRMSD addresses the limitations of traditional RMSD, which can be skewed by flexible regions while neglecting the rigid functional core. The framework includes various analysis modules for structural alignment and clustering, demonstrating significant improvements in performance on benchmarked molecular systems.
BRMSD redefines structural comparison by allowing researchers to focus on rigid domains while minimizing the influence of flexible regions, enhancing the accuracy of biomolecular analyses.
Root-mean-square deviation (RMSD) is the standard metric of structural comparison in molecular dynamics (MD) simulations. In its conventional form, RMSD assigns equal weight to all atoms regardless of mobility. Hence, flexible loops and disordered regions can dominate a global RMSD, while the rigid functional core contributes negligibly to the overall metric. To address this issue, we introduce the Bayes-optimal RMSD (BRMSD), which optimizes per-atom weights jointly with structural averages by maximizing a Bayesian posterior. In a trade-off between low RMSD and weight uniformity, a position-fluctuation parameter $\sigma$ controls the transition from classical RMSD ($\sigma \to \infty$) to a progressive focus on a rigid core ($\sigma \to 0$). The BRMSD framework supports analysis modules for structural alignment, focused alignment onto a user-specified domain, trajectory smoothing, soft $K$-means conformational clustering, and rigid-domain identification. These modules are implemented in the open-source Python package BRMSD (https://github.com/bio-phys/BRMSD) and benchmarked on two MD systems, the endoplasmic reticulum translocon-associated protein SND3 and the phosphotransferase adenylate kinase.
