Abstract
Riboswitches are ligand dependent cis-acting regulatory elements found in diverse species, in some of their mRNA transcripts in the 5’-untranslated regions (UTRs). They are composed of two distinct, yet mutually interacting, domains: a ‘ligand-binding’ or aptamer domain and an expression platform involved in gene regulation. Conformational changes due to ligand binding affect the folding patterns in the expression platform, enabling riboswitches to exercise a variety of ligand dependent controls on downstream gene expressions. Purine riboswitches have been most widely investigated computationally as well as experimentally, for their remarkable ligand (adenine/guanine) specificity and variation in functional mechanism, despite their small size and their sequence similarities. Though the crystal structures of the ligand bound CLOSED states of the aptamer domains of two of the adenine and guanine riboswitches are known, the characterization of the ligand-free OPEN state is limited primarily to results derived from chemical probing, NMR, fluorescence, single labeled FRET and single molecule force spectroscopy experiments. Consequently, a complete molecular level understanding of mechanism of ligand binding, ligand discrimination and switching is still limited. We have addressed this issue for the aptamer domain of the add A-riboswitch by carrying out ab initio quantum chemical studies on the geometries and interaction energies of ligand-bound binding pocket1 and by studying explicit-solvent molecular dynamics simulations of the adenine riboswitch structure (1Y26), both in the adenine-bound (CLOSED) state and in the adenine-free (OPEN) state.2 Dynamic cross correlation (DCC), root mean square fluctuation (RMSF) and backbone torsion angle analyses of the two trajectories and analysis of the average OPEN and CLOSED state structures reveal, alternative hydrogen bonding schemes involving bases that are not directly involved in ligand binding, but, which constitute the lining of the binding pocket and possibly participate in the communication with the expression platform. Our model provides molecular level explanations for existing experimental data, including thermodynamic and kinetic data, related to the ligand binding and switching process.