Sequential Folding of the Nickel/Cobalt Riboswitch Is Facilitated by a Conformational Intermediate: Insights from Single-Molecule Kinetics and Thermodynamics

Effects of mRNA conformational switching on translational noise in gene circuits

Intragenic translational heterogeneity describes the variation in translation at the level of transcripts for an individual gene. A factor that contributes to this source of variation is the mRNA structure. Both the composition of the thermodynamic ensemble, i.e., the stationary distribution of mRNA structures, and the switching dynamics between those play a role. The effect of the switching dynamics on intragenic translational heterogeneity remains poorly understood. We present a stochastic translation model that accounts for mRNA structure switching and is derived from a Markov model via approximate stochastic filtering. We assess the approximation on various timescales and provide a method to quantify how mRNA structure dynamics contributes to translational heterogeneity. With our approach, we allow quantitative information on mRNA switching from biophysical experiments or coarse-grain molecular dynamics simulations of mRNA structures to be included in gene regulatory chemical reaction network models without an increase in the number of species. Thereby, our model bridges a gap between mRNA structure kinetics and gene expression models, which we hope will further improve our understanding of gene regulatory networks and facilitate genetic circuit design.

Mechanism of Fluoride Ion Encapsulation by Magnesium Ions in a Bacterial Riboswitch

Riboswitches sense various ions in bacteria and activate gene expression to synthesize proteins that help maintain ion homeostasis. The crystal structure of the aptamer domain (AD) of the fluoride riboswitch shows that the F- ion is encapsulated by three Mg2+ ions bound to the ligand-binding domain (LBD) located at the core of the AD. The assembly mechanism of this intricate structure is unknown. To this end, we performed computer simulations using coarse-grained and all-atom RNA models to bridge multiple time scales involved in riboswitch folding and ion binding. We show that F- encapsulation by the Mg2+ ions bound to the riboswitch involves multiple sequential steps. Broadly, two Mg2+ ions initially interact with the phosphate groups of the LBD using water-mediated outer-shell coordination and transition to a direct inner-shell interaction through dehydration to strengthen their interaction with the LBD. We propose that the efficient binding mode of the third Mg2+ and F- is that they form a water-mediated ion pair and bind to the LBD simultaneously to minimize the electrostatic repulsion between three Mg2+ bound to the LBD. The tertiary stacking interactions among the LBD nucleobases alone are insufficient to stabilize the alignment of the phosphate groups to facilitate Mg2+ binding. We show that the stability of the whole assembly is an intricate balance of the interactions among the five phosphate groups, three Mg2+, and the encapsulated F- ion aided by the Mg2+ solvated water. These insights are helpful in the rational design of RNA-based ion sensors and fast-switching logic gates.

Ligand-Dependent Volumetric Characterization of Manganese Riboswitch Folding: A High-Pressure Single-Molecule Kinetic Study

Nanoscopic differences in free volume result in pressure-dependent changes in free energies which can therefore impact folding/unfolding stability of biomolecules. Although such effects are typically insignificant under ambient pressure conditions, they are crucially important for deep ocean marine life, where the hydraulic pressure can be on the kilobar scale. In this work, single molecule FRET spectroscopy is used to study the effects of pressure on both the kinetics and overall thermodynamics for folding/unfolding of the manganese riboswitch. Detailed pressure-dependent analysis of the conformational kinetics allows one to extract precision changes (σ ≲ 4-8 ų) in free volumes not only between the fully folded/unfolded conformations but also with respect to the folding transition state of the manganese riboswitch. This permits first extraction of a novel "reversible work" free energy (PΔV) landscape, which reveals a monotonic increase in manganese riboswitch volume along the folding coordinate. Furthermore, such a tool permits exploration of pressure-dependent effects on both Mn²⁺ binding and riboswitch folding, which demonstrate that ligand attachment stabilizes the riboswitch under pressure by decreasing the volume increase upon folding (ΔΔV < 0). Such competition between ligand binding and pressure-induced denaturation dynamics could be of significant evolutionary advantage, compensating for a weakening in riboswitch tertiary structure with pressure-mediated ligand binding and promotion of folding response.

Synergism in the Molecular Crowding of Ligand-Induced Riboswitch Folding: Kinetic/Thermodynamic Insights from Single-Molecule Spectroscopy

Conformational dynamics in riboswitches involves ligand binding and folding of RNA, each of which can be influenced by excluded volume effects under "crowded" in vivo cellular conditions and thus incompletely characterized by in vitro studies under dilute buffer conditions. In this work, temperature-dependent single-molecule fluorescence resonance energy transfer (FRET) spectroscopy is used to characterize the thermodynamics of (i) cognate ligand and (ii) molecular crowders (PEG, polyethylene glycol) on folding of the B. subtilis LysC lysine riboswitch. With the help of detailed kinetic analysis, we isolate and study the effects of PEG on lysine binding and riboswitch folding steps individually, from which we find that PEG crowding facilitates riboswitch folding primarily via a surprising increase in affinity for the cognate ligand. This is furthermore confirmed by temperature-dependent studies, which reveal that PEG crowding is not purely entropic and instead significantly impacts both enthalpic and entropic contributions to the free energy landscape for folding. The results indicate that PEG molecular crowding/stabilization of the lysine riboswitch is more mechanistically complex and requires extension beyond the conventional picture of purely repulsive solvent-solute steric interactions arising from excluded volume and entropy. Instead, the current experimental FRET data support an alternative multistep mechanism, whereby PEG first entropically crowds the unfolded riboswitch into a "pre-folded" conformation, which in turn greatly increases the ligand binding affinity and thereby enhances the overall equilibrium for riboswitch folding.

Reconsidering the czcD (NiCo) Riboswitch as an Iron Riboswitch

Recent work has proposed a new mechanism of bacterial iron regulation: riboswitches that undergo a conformational change in response to Fe^II. The czcD (NiCo) riboswitch was initially proposed to be specific for Ni^II and Co^II, but we recently showed via a czcD-based fluorescent sensor that Fe^II is also a plausible physiological ligand for this riboswitch class. Here, we provide direct evidence that this riboswitch class responds to Fe^II. Isothermal titration calorimetry studies of the native czcD riboswitches from three organisms show no response to Mn^II, a weak response to Zn^II, and similar dissociation constants (∼1 μM) and conformational responses for Fe^II, Co^II, and Ni^II. Only the iron response is in the physiological concentration regime; the riboswitches’ responses to Co^II, Ni^II, and Zn^II require 10³-, 10⁵-, and 10⁶-fold higher “free” metal ion concentrations, respectively, than the typical availability of those metal ions in cells. By contrast, the “Sensei” RNA, recently claimed to be an iron-specific riboswitch, exhibits no response to Fe^II. Our results demonstrate that iron responsiveness is a conserved property of czcD riboswitches and clarify that this is the only family of iron-responsive riboswitch identified to date, setting the stage for characterization of their physiological function.

TPP Riboswitch Populates Holo-Form-like Structure Even in the Absence of Cognate Ligand at High Mg2+ Concentration

Riboswitches are noncoding RNA that regulate gene expression by folding into specific three-dimensional structures (holo-form) upon binding by their cognate ligand in the presence of Mg²⁺. Riboswitch functioning is also hypothesized to be under kinetic control requiring large cognate ligand concentrations. We ask the question under thermodynamic conditions, can the riboswitches populate structures similar to the holo-form only in the presence of Mg²⁺ and absence of cognate ligand binding. We addressed this question using thiamine pyrophosphate (TPP) riboswitch as a model system and computer simulations using a coarse-grained model for RNA. The folding free energy surface (FES) shows that with the initial increase in Mg²⁺ concentration ([Mg²⁺]), the aptamer domain (AD) of TPP riboswitch undergoes a barrierless collapse in its dimensions. On further increase in [Mg²⁺], intermediates separated by barriers appear on the FES, and one of the intermediates has a TPP ligand-binding competent structure. We show that site-specific binding of the Mg²⁺ aids in the formation of tertiary contacts. For [Mg²⁺] greater than physiological concentration, AD folds into a structure similar to the crystal structure of the TPP holo-form even in the absence of the TPP ligand. The folding kinetics shows that TPP AD populates an intermediate due to the misalignment of two arms present in the structure, which acts as a kinetic trap, leading to larger folding timescales. The predictions of the intermediate structures from the simulations are amenable for experimental verification.

Lysine-Dependent Entropy Effects in the B. subtilis Lysine Riboswitch: Insights from Single-Molecule Thermodynamic Studies

Riboswitches play an important role in RNA-based sensing/gene regulation control for many bacteria. In particular, the accessibility of multiple conformational states at physiological temperatures allows riboswitches to selectively bind a cognate ligand in the aptamer domain, which triggers secondary structural changes in the expression platform, and thereby "switching" between on or off transcriptional or translational states for the downstream RNA. The present work exploits temperature-controlled, single-molecule total internal reflection fluorescence (TIRF) microscopy to study the thermodynamic landscape of such ligand binding/folding processes, specifically for the Bacillus subtilis lysine riboswitch. The results confirm that the riboswitch folds via an induced-fit (IF) mechanism, in which cognate lysine ligand first binds to the riboswitch before structural rearrangement takes place. The transition state to folding is found to be enthalpically favored (ΔH_fold^‡ < 0), yet with a free-energy barrier that is predominantly entropic (-TΔS_fold^‡ > 0), which results in folding (unfolding) rate constants strongly dependent (independent) of lysine concentration. Analysis of the single-molecule kinetic "trajectories" reveals this rate constant dependence of k_fold on lysine to be predominantly entropic in nature, with the additional lysine conferring preferential advantage to the folding process by the presence of ligands correctly oriented with respect to the riboswitch platform. By way of contrast, van't Hoff analysis reveals enthalpic contributions to the overall folding thermodynamics (ΔH⁰) to be surprisingly constant and robustly independent of lysine concentration. The results demonstrate the crucial role of hydrogen bonding between the ligand and riboswitch platform but with only a relatively modest fraction (45%) of the overall enthalpy change needed to access the transition state and initiate transcriptional switching.

Effects of Molecular Crowders on Single-Molecule Nucleic Acid Folding: Temperature-Dependent Studies Reveal True Crowding vs Enthalpic Interactions

Biomolecular folding in cells can be strongly influenced by spatial overlap/excluded volume interactions (i.e., "crowding") with intracellular solutes. As a result, traditional in vitro experiments with dilute buffers may not accurately recapitulate biomolecule folding behavior in vivo. In order to account for such ubiquitous excluded volume effects, biologically inert polyethylene glycol (PEG) and polysaccharides (dextran and Ficoll) are often used as in vitro crowding agents to mimic in vivo crowding conditions, with a common observation that high concentrations of these polymers stabilize the more compact biomolecule conformation. However, such an analysis can be distorted by differences in polymer interactions with the folded vs unfolded conformers, requiring temperature-dependent analysis of the thermodynamics to reliably assess competing enthalpic vs entropic contributions and thus the explicit role of excluded volume. In this work, temperature-controlled single-molecule fluorescence resonance energy transfer (smFRET) is used to characterize the thermodynamic interaction between nucleic acids and common polymer crowders PEG, dextran, and Ficoll. The results reveal that PEG promotes secondary and tertiary nucleic acid folding by simultaneously increasing the folding rate while decreasing the unfolding rate, with temperature-dependent studies confirming that the source of PEG stabilization is predominantly entropic and consistent with a true excluded volume crowding mechanism. By way of contrast, neither dextran nor Ficoll induces any significant concentration-dependent change in nucleic acid folding stability at room temperature, but instead, stabilization effects gradually appear with a temperature increase. Such a thermal response indicates that both folding enthalpies and entropies are impacted by dextran and Ficoll. A detailed thermodynamic analysis of the kinetics suggests that, instead of true entropic molecular crowding, dextran and Ficoll associate preferentially with the unfolded vs folded nucleic acid conformer as a result of larger solvent accessible surface area, thereby skewing the free energy landscapes through both significant entropic/enthalpic contributions that compete and fortuitously cancel near room temperature.

Smaller molecules crowd better: Crowder size dependence revealed by single-molecule FRET studies and depletion force modeling analysis

The cell is an extremely crowded environment, which is known to have a profound impact on the thermodynamics, functionality, and conformational stability of biomolecules. Speculations from recent theoretical molecular dynamics studies suggest an intriguing size dependence to such purely entropic crowding effects, whereby small molecular weight crowders under constant enthalpy conditions are more effective than larger crowders on a per volume basis. If experimentally confirmed, this would be profoundly significant, as the cellular cytoplasm is also quite concentrated in smaller molecular weight solutes such as inorganic ions, amino acids, and various metabolites. The challenge is to perform such studies isolating entropic effects under isoenthalpic conditions. In this work, we first present results from single-molecule FRET spectroscopy (smFRET) on the molecular size-dependent crowding stabilization of a simple RNA tertiary motif (the GAAA tetraloop-tetraloop receptor), indeed providing evidence in support of the surprising notion in the crowding literature that "smaller is better." Specifically, systematic smFRET studies as a function of crowder solute size reveal that smaller molecules both significantly increase the RNA tertiary folding rate and, yet, simultaneously decrease the unfolding rate, predicting strongly size-dependent stabilization of RNA tertiary structures under crowded cellular conditions. The size dependence of these effects has been explored via systematic variation of crowder size over a broad range of molecular weights (90-3000 amu). Furthermore, corresponding temperature dependent studies indicate the systematic changes in the folding equilibrium to be predominantly entropic in origin, i.e., consistent with a fundamental picture of entropic molecular crowding without additional enthalpic interactions. Most importantly, all trends in the single-molecule crowding data can be quantitatively recapitulated by a simple analytic depletion force model, whereby excluded volume interactions represent the major thermodynamic driving force toward folding. Our study, thus, not only provides experimental evidence and theoretical support for small molecule crowding but also predicts further enhancement of crowding effects for even smaller molecules on a per volume basis.