Dynamic Energy Landscapes of Riboswitches Help Interpret Conformational Rearrangements and Function

Riboswitches are RNAs that modulate gene expression by ligand-induced conformational changes. However, the way in which sequence dictates alternative folding pathways of gene regulation remains unclear. In this study, we compute energy landscapes, which describe the accessible secondary structures for a range of sequence lengths, to analyze the transcriptional process as a given sequence elongates to full length. In line with experimental evidence, we find that most riboswitch landscapes can be characterized by three broad classes as a function of sequence length in terms of the distribution and barrier type of the conformational clusters: low-barrier landscape with an ensemble of different conformations in equilibrium before encountering a substrate; barrier-free landscape in which a direct, dominant “downhill” pathway to the minimum free energy structure is apparent; and a barrier-dominated landscape with two isolated conformational states, each associated with a different biological function. Sharing concepts with the “new view” of protein folding energy landscapes, we term the three sequence ranges above as the sensing, downhill folding, and functional windows, respectively. We find that these energy landscape patterns are conserved in various riboswitch classes, though the order of the windows may vary. In fact, the order of the three windows suggests either kinetic or thermodynamic control of ligand binding. These findings help understand riboswitch structure/function relationships and open new avenues to riboswitch design

Three-dimensional RNA structure refinement by hydroxyl radical probing

Molecular modeling guided by experimentally-derived structural information is an attractive approach for three-dimensional structure determination of complex RNAs that are not amenable to study by high-resolution methods. Hydroxyl radical probing (HRP), performed routinely in many laboratories, provides a measure of solvent accessibility at individual nucleotides. HRP measurements have, to date, only been used to evaluate RNA models qualitatively. Here, we report development of a quantitative structure refinement approach using HRP measurements to drive discrete molecular dynamics simulations for RNAs ranging in size from 80 to 230 nucleotides. HRP reactivities were first used to identify RNAs that form extensive helical packing interactions. For these RNAs, we achieved highly significant structure predictions, given inputs of RNA sequence and base pairing. This HRP-directed tertiary structure refinement approach generates robust structural hypotheses useful for guiding explorations of structure-function interrelationships in RNA

The structural basis of cephalosporin formation in a mononuclear ferrous enzyme

Deacetoxycephalosporin-C synthase (DAOCS) is a mononuclear ferrous enzyme that transforms penicillins into cephalosporins by inserting a carbon atom into the penicillin nucleus. In the first half-reaction, dioxygen and 2-oxoglutarate produce a reactive iron-oxygen species, succinate and CO2. The oxidizing iron species subsequently reacts with penicillin to give cephalosporin and water. Here we describe high-resolution structures for ferrous DAOCS in complex with penicillins, the cephalosporin product, the cosubstrate and the coproduct. Steady-state kinetic data, quantum-chemical calculations and the new structures indicate a reaction sequence in which a ‘booby-trapped’ oxidizing species is formed. This species is stabilized by the negative charge of succinate on the iron. The binding sites of succinate and penicillin overlap, and when penicillin replaces succinate, it removes the stabilizing charge, eliciting oxidative attack on itself. Requisite groups of penicillin are within 1 Å of the expected position of a ferryl oxygen in the enzyme–penicillin complex.