Structural Ensembles of Ribonucleic Acids From Solvent Accessibility Data: Application to the S-Adenosylmethionine (SAM)-Responsive Riboswitch

Abstract

Riboswitches are regulatory ribonucleic acid (RNA) elements that act as ligand-dependent conformational switches. In the apo form, the aptamer domain, the region of a riboswitch that binds to its cognate ligand, is dynamic, thus requiring an ensemble-representation of its structure. Analysis of such ensembles can provide molecular insights into the sensing mechanism and capabilities of riboswitches. Here, as a proof-of-concept, we constructed a pair of atomistic ensembles of the well-studied S-adenosylmethionine (SAM)-responsive riboswitch in the absence (-SAM) and presence (+SAM) of SAM. To achieve this, we first generated a large conformational pool and then reweighted conformers in the pool using solvent accessible surface area (SASA) data derived from recently reported light-activated structural examination of RNA (LASER) reactivities, measured in the -SAM and +SAM states of the riboswitch. The differences in the resulting -SAM and +SAM ensembles are consistent with a SAM-dependent reshaping of the free landscape of the riboswitch. Interestingly, within the -SAM ensemble, we identified a conformer that harbors a hidden binding pocket, which was discovered using ensemble docking. The method we have applied to the SAM riboswitch is general, and could, therefore, be used to construct atomistic ensembles for other riboswitches, and more broadly, other classes of structured RNAs.