Michael Brenowitz

Monovalent cations mediate formation of native tertiary structure of the Tetrahymena thermophila ribozyme

The formation of individual tertiary contacts of the Tetrahymena L-21 Sca I ribozyme has been monitored by hydroxyl radical footprinting and its global conformation by analytical ultracentrifugation as a function of monovalent ion concentration in the absence of divalent ions. Advanced methods of data analysis, which allow the hydroxyl radical reactivity of every nucleotide to be quantified, permit monitoring of each and every structural element of the RNA. Monovalent ion-mediated global compaction of the ribozyme is accompanied by the formation of native tertiary contacts; most native tertiary contacts are evident except several that are located near where divalent ions are observed in crystallographic structures. Non-native tertiary contacts are also observed at low but not high concentrations of monovalent ions. In light of recent studies that have shown that the presence of monovalent ions greatly accelerates the Mg2+-dependent folding of the Tetrahymena ribozyme, the present studies suggest that Na+ concentration changes not only the starting position of the RNA on its folding funnel but also pushes it deep into the well by forming native tertiary contacts and, thus, favoring fast and correct folding pathways.

Semiautomated and rapid quantification of nucleic acid footprinting and structure mapping experiments

We have developed protocols for rapidly quantifying the band intensities from nucleic acid chemical mapping gels at single-nucleotide resolution. These protocols are implemented in the software SAFA (semi-automated footprinting analysis) that can be downloaded without charge from http://safa.stanford.edu. The protocols implemented in SAFA have five steps: (i) lane identification, (ii) gel rectification, (iii) band assignment, (iv) model fitting and (v) band-intensity normalization. SAFA enables the rapid quantitation of gel images containing thousands of discrete bands, thereby eliminating a bottleneck to the analysis of chemical mapping experiments. An experienced user of the software can quantify a gel image in ∼20 min. Although SAFA was developed to analyze hydroxyl radical (·OH) footprints, it effectively quantifies the gel images obtained with other types of chemical mapping probes. We also present a series of tutorial movies that illustrate the best practices and different steps in the SAFA analysis as a supplement to this protocol.

Distinct contribution of electrostatics, initial conformational ensemble, and macromolecular stability in RNA folding

We distinguish the contribution of the electrostatic environment, initial conformational ensemble, and macromolecular stability on the folding mechanism of a large RNA using a combination of time-resolved “Fast Fenton” hydroxyl radical footprinting and exhaustive kinetic modeling. This integrated approach allows us to define the folding landscape of the L-21 Tetrahymena thermophila group I intron structurally and kinetically from its earliest steps with unprecedented accuracy. Distinct parallel pathways leading the RNA to its native form upon its Mg2+-induced folding are observed. The structures of the intermediates populating the pathways are not affected by variation of the concentration and type of background monovalent ions (electrostatic environment) but are altered by a mutation that destabilizes one domain of the ribozyme. Experiments starting from different conformational ensembles but folding under identical conditions show that whereas the electrostatic environment modulates molecular flux through different pathways, the initial conformational ensemble determines the partitioning of the flux. This study showcases a robust approach for the development of kinetic models from collections of local structural probes.

Stability and cooperativity of individual tertiary contacts in RNA revealed through chemical denaturation

For proteins, understanding tertiary interactions involved in local versus global unfolding has become increasingly important for understanding the nature of the native state ensemble, the mechanisms of unfolding, and the stability of both the native and intermediate states in folding. In this work we have addressed related questions with respect to RNA structure by combining chemical denaturation and hydroxyl radical footprinting methods. We have determined unfolding isotherms for each of 26 discrete sites of protection located throughout the Tetrahymena thermophila group I ribozyme. The cooperativity of folding, m-value, and the free energy, ΔG°N−U, associated with formation of each tertiary contact was determined by analysis of the isotherms. The ΔG°N−U values measured in this study vary from 1.7 ± 0.2 to 7.6 ± 1.2 kcal mol−1. Thus, the stability of these discrete tertiary contacts vary by almost 104. In addition, an intradomain contact and three interdomain contacts show high cooperativity (m-values of 1.1 ± 0.2 to 1.7 ± 0.3 kcal mol−1 M−1) indicating that these contacts exhibit global cooperatively in their folding behavior. This new approach to examining RNA stability provides an exciting comparison to our understanding of protein structure and folding mechanisms.

Conformational States of Human Purine Nucleoside Phosphorylase at Rest, at Work, and with Transition State Analogues

Human purine nucleoside phosphorylase (PNP) is a homotrimer binding tightly to the transition state analogues Immucillin-H (ImmH; K(d) = 56 pM) and DATMe-ImmH-Immucillin-H (DATMe-ImmH; K(d) = 8.6 pM). ImmH binds with a larger entropic penalty than DATMe-ImmH, a chemically more flexible inhibitor. The testable hypothesis is that PNP conformational states are more relaxed (dynamic) with DATMe-ImmH, despite tighter binding than with ImmH. PNP conformations are probed by peptide amide deuterium exchange (HDX) using liquid chromatography high-resolution Fourier transform ion cyclotron resonance mass spectrometry and by sedimentation rates. Catalytically equilibrating Michaelis complexes (PNP.PO(4).inosine <--> PNP.Hx.R-1-P) and inhibited complexes (PNP.PO(4).DATMe-ImmH and PNP.PO(4).ImmH) show protection from HDX at 9, 13, and 15 sites per subunit relative to resting PNP (PNP.PO(4)) in extended incubations. The PNP.PO(4).ImmH complex is more compact (by sedimentation rate) than the other complexes. HDX kinetic analysis of ligand-protected sites corresponds to peptides near the catalytic sites. HDX and sedimentation results establish that PNP protein conformation (dynamic motion) correlates more closely with entropy of binding than with affinity. Catalytically active turnover with saturated substrate sites causes less change in HDX and sedimentation rates than binding of transition state analogues. DATMe-ImmH more closely mimics the transition of human PNP than does ImmH and achieves strong binding interactions at the catalytic site while causing relatively modest alterations of the protein dynamic motion. Transition state analogues causing the most rigid, closed protein conformation are therefore not necessarily the most tightly bound. Close mimics of the transition state are hypothesized to retain enzymatic dynamic motions related to transition state formation.

RNA molecules with conserved catalytic cores but variable peripheries fold along unique energetically optimized pathways

Functional and kinetic constraints must be efficiently balanced during the folding process of all biopolymers. To understand how homologous RNA molecules with different global architectures fold into a common core structure we determined, under identical conditions, the folding mechanisms of three phylogenetically divergent group I intron ribozymes. These ribozymes share a conserved functional core defined by topologically equivalent tertiary motifs but differ in their primary sequence, size, and structural complexity. Time-resolved hydroxyl radical probing of the backbone solvent accessible surface and catalytic activity measurements integrated with structural-kinetic modeling reveal that each ribozyme adopts a unique strategy to attain the conserved functional fold. The folding rates are not dictated by the size or the overall structural complexity, but rather by the strength of the constituent tertiary motifs which, in turn, govern the structure, stability, and lifetime of the folding intermediates. A fundamental general principle of RNA folding emerges from this study: The dominant folding flux always proceeds through an optimally structured kinetic intermediate that has sufficient stability to act as a nucleating scaffold while retaining enough conformational freedom to avoid kinetic trapping. Our results also suggest a potential role of naturally selected peripheral A-minor interactions in balancing RNA structural stability with folding efficiency.