Bayesian Learning of Distance Metrics Beyond RMSD for Biomolecule Alignment, Clustering, and Domain Identification

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 (<i>BRMSD</i>), 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 σ controls the transition from classical RMSD (σ → ∞) to a progressive focus on a rigid core (σ → 0). The BRMSD framework supports analysis modules for structural alignment, focused alignment onto a user-specified domain, trajectory smoothing, soft <i>K</i>-means conformational clustering, and rigid-domain identification.