Exploring Conformational Transitions of Adenine RNA Dimer via Machine Learning Potentials

Abstract

RNA is a flexible biopolymer that adopts diverse conformations while forming structural motifs essential for its function. Classical RNA force fields often show limited transferability and inefficient sampling of transitions between stable states, particularly in moderately large RNA. To address these limitations, quantum-informed machine learning (ML) potentials have recently emerged as a promising alternative, offering improved accuracy and transferability relative to classical force fields. Here, we assess ML potentials for exploring RNA conformations using the adenine–adenine dinucleoside monophosphate (ApA) dimer, a fundamental RNA building block. We generated an extensive quantum-mechanical (QM) dataset of physicochemical properties for ApA conformations obtained from temperature replica exchange molecular dynamics (TREMD) simulations. Despite its small size, the ApA dimer exhibits a complex energetic landscape with six well-defined conformational clusters in which quantum effects and solvent-mediated interactions play a crucial role. Using this dataset, we parametrized ML potentials based on the equivariant MACE architecture and informed by both ab initio and semiempirical property data. The resulting potentials reproduce key conformational features of the ApA system, including base stacking, sugar geometry, and backbone flexibility, and provide broader coverage of structural transitions than the general-purpose SO3LR and MACE-POLAR-1 models. These findings underscore the importance of comprehensive QM datasets for RNA building blocks to support the structural and energetic characterization of RNA complexes and emphasize the need for robust and efficient validation metrics for ML potentials.