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

The determination of the three-dimensional structure of large RNA macromolecules in solution is a challenging task that often requires the use of several experimental and computational techniques. Small-angle X-ray scattering can provide insight into some geometrical properties of the probed molecule, but this data must be properly interpreted in order to generate a three-dimensional model. Here, we propose a multiscale pipeline which introduces SAXS data into modelling the global shape of RNA in solution, which can be hierarchically refined until reaching atomistic precision in explicit solvent. The low-resolution helix model (Ernwin) deals with the exploration of the huge conformational space making use of the SAXS data, while a nucleotide-level model (SPQR) removes clashes and disentangles the proposed structures, leading the structure to an all-atom representation in explicit water. We apply the procedure on four different known pdb structures up to 159 nucleotides with promising results. Additionally, we predict an all-atom structure for the Plasmodium falceparum signal recognition particle ALU RNA based on SAXS data deposited in the SASBDB, which has an alternate conformation and better fit to the SAXS data than the previously published structure based on the same data but other modelling methods.