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

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.

arXiv commentSupporting information included in ancillary files. This version includes corrections implemented after receiving feedbacks from the community

References in corpus (14)

• SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
• Canonical sampling through velocity-rescaling
• Escaping free-energy minima
• PLUMED 2: New feathers for an old bird
• Well-Tempered Metadynamics: A Smoothly Converging and Tunable Free-Energy Method
• Building Water Models, A Different Approach
• Metadynamics with adaptive Gaussians
• Replica-exchange multicanonical algorithm and multicanonical replica-exchange method for simulating systems with rough energy landscape
• Hamiltonian replica-exchange in GROMACS: a flexible implementation
• The Role of Nucleobase Interactions in RNA Structure and Dynamics
• Using the Maximum Entropy Principle to Combine Simulations and Solution Experiments
• Free Energy Landscape of GAGA and UUCG RNA Tetraloops
• Exploring RNA structure and dynamics through enhanced sampling simulations
• Accurate multiple time step in biased molecular simulations

Cited by in corpus (2)

• Integrating experimental data with molecular simulations to investigate RNA structural dynamics
• Refinement of molecular dynamics ensembles using experimental data and flexible forward models