Abstract:Small-angle X-ray scattering (SAXS) experiments are increasingly used to probe RNA structure. A number of \emph{forward models} that relate measured SAXS intensities and structural features, and that are suitable to model either explicit-solvent effects or solute dynamics, have been proposed in the past years. Here we introduce an approach that integrates atomistic molecular dynamics simulations and SAXS experiments to reconstruct RNA structural ensembles while simultaneously accounting for both RNA conformational dynamics and explicit-solvent effects. Our protocol exploits SAXS pure-solute forward models and enhanced sampling methods to sample an heterogenous ensemble of structures, with no information towards the experiments provided on-the-fly. The generated structural ensemble is then reweighted through the maximum entropy principle so as to match reference SAXS experimental data at multiple ionic conditions. Importantly, accurate explicit-solvent forward models are used at this reweighting stage. We apply this framework to the GTPase-associated center, a relevant RNA molecule involved in protein translation, in order to elucidate its ion-dependent conformational ensembles. We show that (a) both solvent and dynamics are crucial to reproduce experimental SAXS data and (b) the resulting dynamical ensembles contain an ion-dependent fraction of extended structures.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to reconstruct the structural ensemble of RNA under different ionic conditions by combining atom - molecular dynamics simulations and small - angle X - ray scattering (SAXS) experiments, while taking into account RNA conformational dynamics and explicit solvent effects. Specifically, the paper aims to: 1. **Develop a new method**: Integrate atom - molecular dynamics simulations and SAXS experimental data to reconstruct the structural ensemble of RNA. This method takes into account not only the conformational dynamics of RNA but also the explicit solvent effects. 2. **Verify the effectiveness of the method**: By applying this method to study the GTPase - associated center (GAC), a 57 - nucleotide - long RNA molecule related to protein translation, verify the effectiveness and accuracy of the method under different ionic conditions. 3. **Explore the influence of ions on RNA structure**: In particular, study how magnesium ions (Mg²⁺) and potassium ions (K⁺) affect the structural ensemble of GAC, and how these ions cause the transition of RNA molecules from a compact structure to an extended structure. ### Main contributions of the paper - **Method innovation**: Propose a new method that combines the pure solute forward model and enhanced sampling methods to generate heterogeneous structural ensembles, and then re - weight them through the maximum entropy principle to match the reference SAXS experimental data. - **Experimental verification**: Verify the effectiveness of this method through the SAXS experimental data of GAC RNA molecules under different ionic conditions. - **Combination of theory and experiment**: Demonstrate how to more accurately understand the dynamic behavior and structural changes of RNA molecules in solution through the combination of theoretical simulations and experimental data. ### Specific implementation - **Molecular dynamics simulation**: Use the AMBER force field and 4 - point optimal point charge (OPC) water model for simulation, taking into account different ionic conditions. - **Enhanced sampling**: Adopt the Hamiltonian replica exchange (HREX) method and the metadynamics method to enhance the sampling of extended structures. - **Inverse calculation of SAXS spectra**: Use different software (such as PLUMED, Crysol, WAXSiS, and Capriqorn) to calculate SAXS spectra, among which Capriqorn explicitly takes into account the solvent effect. - **Re - weighting**: Re - weight the generated structural ensemble through the maximum entropy principle to match the experimental SAXS data. ### Results - **Differences in SAXS spectra**: The SAXS spectra calculated by different methods show significant differences in the high - q value region, especially when comparing the methods that explicitly consider the solvent effect (such as WAXSiS and Capriqorn) with those that only consider pure solutes (such as PLUMED and Crysol). - **Reconstruction of structural ensembles**: Through re - weighting, the structural ensembles of GAC under different ionic conditions are successfully reconstructed. In particular, under K⁺ conditions, a significant increase in the extended structure is observed. - **Statistical analysis**: Through the calculation of Kish effective sample size, it is found that only a few structures make significant contributions to the final SAXS spectra, which reflects the insufficient sampling of extended structures. In conclusion, this paper provides a new method to study the structural dynamics of RNA under different ionic conditions by combining molecular dynamics simulations and SAXS experimental data, providing an important tool for understanding and predicting the behavior of RNA molecules.

Reweighting of molecular simulations with explicit-solvent SAXS restraints elucidates ion-dependent RNA ensembles

Learning Effective Molecular Models from Experimental Observables.

Combining simulations and solution experiments as a paradigm for RNA force field refinement

Sampling globally and locally correct RNA 3D structures using Ernwin, SPQR and experimental SAXS data

Selective deuteration of an RNA:RNA complex for structural analysis using small-angle scattering

Structure and ensemble refinement against SAXS data: combining MD simulations with Bayesian inference or with the maximum entropy principle

Ensemble Refinement of mismodeled cryo-EM RNA Structures Using All-Atom Simulations

Small-Angle X-ray Scattering (SAXS) Combined with SAXS-Driven Molecular Dynamics for Structural Analysis of Multistranded RNA Assemblies

Molecular Dynamics Simulations with Grand-Canonical Reweighting Suggest Cooperativity Effects in RNA Structure Probing Experiments

Refinement of molecular dynamics ensembles using experimental data and flexible forward models

Bayesian refinement of protein structures and ensembles against SAXS data using molecular dynamics

Conformational ensembles of an RNA hairpin using molecular dynamics and sparse NMR data

Deriving RNA topological structure from SAXS

Characterizing RNA ensembles from NMR data with kinematic models

Predicting solution scattering patterns with explicit-solvent molecular simulations

Modeling Rna Topological Structures Using Small Angle X-Ray Scattering

A novel implicit solvent model for simulating the molecular dynamics of RNA.

SAXS ensemble refinement of ESCRT-III CHMP3 conformational transitions

Combining molecular dynamics simulations with small-angle X-ray and neutron scattering data to study multi-domain proteins in solution.

Predicting RNA Structure and Dynamics with Deep Learning and Solution Scattering

Saxs-Oriented Ensemble Refinement of Flexible Biomolecules