Novel Heat-Promoted Folding Dynamics of the yybP-ykoY Manganese Riboswitch: Kinetic and Thermodynamic Studies at the Single-Molecule Level

Two-dimensional (P/T) studies of secondary/tertiary conformational dynamics in nucleic acids: pressure induced melting and Maxwell relations at the single molecule level

A predictive understanding of conformational folding of nucleic acids depends crucially on the underlying competition between enthalpic and entropic contributions to overall free energy changes. In extreme environments (e.g. deep ocean vents), such free energy changes are in turn impacted by both pressure and temperature as strongly coupled intensive variables, emphasizing the importance of a detailed molecular level understanding of the underlying thermodynamics. In this work, single-molecule fluorescence energy resonance transfer (smFRET) microscopy methods are implemented for quantitative kinetic study of secondary structure (i.e., DNA hybridization) and tertiary structure (i.e., RNA Mn2+ riboswitch folding) equilibria as a function of both pressure (P = 1 to 1000 bar) and temperature (T = 21 to 27 °C). Temperature dependent studies at a series of fixed pressures reveal the single molecule DNA and RNA constructs to be stabilized and destabilized, respectively, with increasing T. Interestingly, results for the Mn2+ riboswitch indicate a positive entropy change (ΔS > 0) for achieving the native tertiary conformation at all external pressures. This is contrary to more common physical expectations of increased order and thus an entropically penalty for tertiary folding of the RNA. On the other hand, pressure dependent scans for a series of constant temperature conditions confirm that both the DNA hairpin and RNA riboswitch constructs destabilize ("melt") under increasing pressure, which by van't Hoff analysis implies a positive volume change (ΔV0 > 0) for both secondary and tertiary folding into the native state. Furthermore, slices through these two-dimensional free energy surfaces permit parameters for isobaric thermal expansion in both secondary and tertiary conformational folding coordinates to be extracted. Finally, experimental control of multiple variables allows determination of the folding free energies as a two-dimensional function of pressure and temperature (ΔG(P, T)). To the best of our knowledge, this facilitates a first experimental confirmation of the underlying Maxwell relation (∂ΔS/∂P)T = -(∂ΔV/∂T)P for the exact free energy differential dG(P, V, S, T) at the single molecule level.

Thermodynamic compensation to temperature extremes in B. subtilis vs T. maritima lysine riboswitches

T. maritima and B. subtilis are bacteria that inhabit significantly different thermal environments, ∼80 vs. ∼40°C, yet employ similar lysine riboswitches to aid in the transcriptional regulation of the genes involved in the synthesis and transport of amino acids. Despite notable differences in G-C basepair frequency and primary sequence, the aptamer moieties of each riboswitch have striking similarities in tertiary structure, with several conserved motifs and long-range interactions. To explore genetic adaptation in extreme thermal environments, we compare the kinetic and thermodynamic behaviors in T. maritima and B. subtilis lysine riboswitches via single-molecule fluorescence resonance energy transfer analysis. Kinetic studies reveal that riboswitch folding rates increase with lysine concentration while the unfolding rates are independent of lysine. This indicates that both riboswitches bind lysine through an induced-fit ("bind-then-fold") mechanism, with lysine binding necessarily preceding conformational changes. Temperature-dependent van't Hoff studies reveal qualitative similarities in the thermodynamic landscapes for both riboswitches in which progression from the open, lysine-unbound state to both transition states (‡) and closed, lysine-bound conformations is enthalpically favored yet entropically penalized, with comparisons of enthalpic and entropic contributions extrapolated to a common [K+] = 100 mM in quantitative agreement. Finally, temperature-dependent Eyring analysis reveals the TMA and BSU riboswitches to have remarkably similar folding/unfolding rate constants when extrapolated to their respective (40 and 80°C) environmental temperatures. Such behavior suggests a shared strategy for ligand binding and aptamer conformational change in the two riboswitches, based on thermodynamic adaptations in number of G-C basepairs and/or modifications in tertiary structure that stabilize the ligand-unbound conformation to achieve biocompetence under both hyperthermophilic and mesothermophilic conditions.

The effect of pseudoknot base pairing on cotranscriptional structural switching of the fluoride riboswitch

A central question in biology is how RNA sequence changes influence dynamic conformational changes during cotranscriptional folding. Here we investigated this question through the study of transcriptional fluoride riboswitches, non-coding RNAs that sense the fluoride anion through the coordinated folding and rearrangement of a pseudoknotted aptamer domain and a downstream intrinsic terminator expression platform. Using a combination of Escherichia coli RNA polymerase in vitro transcription and cellular gene expression assays, we characterized the function of mesophilic and thermophilic fluoride riboswitch variants. We showed that only variants containing the mesophilic pseudoknot function at 37°C. We next systematically varied the pseudoknot sequence and found that a single wobble base pair is critical for function. Characterizing thermophilic variants at 65°C through Thermus aquaticus RNA polymerase in vitro transcription showed the importance of this wobble pair for function even at elevated temperatures. Finally, we performed all-atom molecular dynamics simulations which supported the experimental findings, visualized the RNA structure switching process, and provided insight into the important role of magnesium ions. Together these studies provide deeper insights into the role of riboswitch sequence in influencing folding and function that will be important for understanding of RNA-based gene regulation and for synthetic biology applications.

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.

Multidimensional Applications and Challenges of Riboswitches in Biosensing and Biotherapy

Riboswitches have received significant attention over the last two decades for their multiple functionalities and great potential for applications in various fields. This article highlights and reviews the recent advances in biosensing and biotherapy. These fields involve a wide range of applications, such as food safety detection, environmental monitoring, metabolic engineering, live cell imaging, wearable biosensors, antibacterial drug targets, and gene therapy. The discovery, origin, and optimization of riboswitches are summarized to help readers better understand their multidimensional applications. Finally, this review discusses the multidimensional challenges and development of riboswitches in order to further expand their potential for novel applications.

The Effect of Pseudoknot Base Pairing on Cotranscriptional Structural Switching of the Fluoride Riboswitch

A central question in biology is how RNA sequence changes influence dynamic conformational changes during cotranscriptional folding. Here we investigated this question through the study of transcriptional fluoride riboswitches, non-coding RNAs that sense the fluoride anion through the coordinated folding and rearrangement of a pseudoknotted aptamer domain and a downstream intrinsic terminator expression platform. Using a combination of E. coli RNA polymerase in vitro transcription and cellular gene expression assays, we characterized the function of mesophilic and thermophilic fluoride riboswitch variants. We showed that only variants containing the mesophilic pseudoknot function at 37 °C. We next systematically varied the pseudoknot sequence and found that a single wobble base pair is critical for function. Characterizing thermophilic variants at 65 °C through Thermus aquaticus RNA polymerase in vitro transcription showed the importance of this wobble pair for function even at elevated temperatures. Finally, we performed all-atom molecular dynamics simulations which supported the experimental findings, visualized the RNA structure switching process, and provided insight into the important role of magnesium ions. Together these studies provide deeper insights into the role of riboswitch sequence in influencing folding and function that will be important for understanding of RNA-based gene regulation and for synthetic biology applications.

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.

Entropy Driving the Mg2+-Induced Folding of TPP Riboswitch RNA

Mg²⁺ is well known to facilitate the structural folding of RNA. However, the thermodynamic and dynamic roles of Mg²⁺ in RNA folding remain elusive. Here, we exploit single-molecule fluorescence resonance energy transfer (smFRET) and isothermal titration calorimetry (ITC) to study the mechanism of Mg²⁺ in facilitating the folding of thiamine pyrophosphate (TPP) riboswitch RNA. The results of smFRET identify that the presence of Mg²⁺ compacts the RNA and enlarges the conformational dispersity among individual RNA molecules, resulting in a large gain of entropy. The compact yet flexible conformations triggered by Mg²⁺ may help the riboswitch recognize its specific ligand and further fold. This is supported by the ITC experiments, in which the Mg²⁺-induced RNA folding is driven by entropy (ΔS) instead of enthalpy (ΔH). Our results complement the understanding of the Mg²⁺-induced RNA folding. The strategy developed in this work can be used to model other RNAs' folding under different conditions.

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.

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.

Applications of isothermal titration calorimetry in pure and applied research from 2016 to 2020

The last 5 years have seen a series of advances in the application of isothermal titration microcalorimetry (ITC) and interpretation of ITC data. ITC has played an invaluable role in understanding multiprotein complex formation including proteolysis-targeting chimeras (PROTACS), and mitochondrial autophagy receptor Nix interaction with LC3 and GABARAP. It has also helped elucidate complex allosteric communication in protein complexes like trp RNA-binding attenuation protein (TRAP) complex. Advances in kinetics analysis have enabled the calculation of kinetic rate constants from pre-existing ITC data sets. Diverse strategies have also been developed to study enzyme kinetics and enzyme-inhibitor interactions. ITC has also been applied to study small molecule solvent and solute interactions involved in extraction, separation, and purification applications including liquid-liquid separation and extractive distillation. Diverse applications of ITC have been developed from the analysis of protein instability at different temperatures, determination of enzyme kinetics in suspensions of living cells to the adsorption of uremic toxins from aqueous streams.

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.

High pressure single-molecule FRET studies of the lysine riboswitch: cationic and osmolytic effects on pressure induced denaturation

Deep sea biology is known to thrive at pressures up to ≈1 kbar, which motivates fundamental biophysical studies of biomolecules under such extreme environments. In this work, the conformational equilibrium of the lysine riboswitch has been systematically investigated by single molecule FRET (smFRET) microscopy at pressures up to 1500 bar. The lysine riboswitch preferentially unfolds with increasing pressure, which signals an increase in free volume (ΔV0 > 0) upon folding of the biopolymer. Indeed, the effective lysine binding constant increases quasi-exponentially with pressure rise, which implies a significant weakening of the riboswitch-ligand interaction in a high-pressure environment. The effects of monovalent/divalent cations and osmolytes on folding are also explored to acquire additional insights into cellular mechanisms for adapting to high pressures. For example, we find that although Mg2+ greatly stabilizes folding of the lysine riboswitch (ΔΔG0 < 0), there is negligible impact on changes in free volume (ΔΔV0 ≈ 0) and thus any pressure induced denaturation effects. Conversely, osmolytes (commonly at high concentrations in deep sea marine species) such as the trimethylamine N-oxide (TMAO) significantly reduce free volumes (ΔΔV0 < 0) and thereby diminish pressure-induced denaturation. We speculate that, besides stabilizing RNA structure, enhanced levels of TMAO in cells might increase the dynamic range for competent riboswitch folding by suppressing the pressure-induced denaturation response. This in turn could offer biological advantage for vertical migration of deep-sea species, with impacts on food searching in a resource limited environment.

DNA Hairpin Hybridization under Extreme Pressures: A Single-Molecule FRET Study

Organisms have evolved to live in a variety of complex environments, which clearly has required cellular biology to accommodate to extreme conditions of hydraulic pressure and elevated temperature. In this work, we exploit single-molecule Forster resonance energy transfer (FRET) spectroscopy to probe structural changes in DNA hairpins as a function of pressure and temperature, which allows us to extract detailed thermodynamic information on changes in free energy (ΔG°), free volume (ΔV°), enthalpy (ΔH°), and entropy (ΔS°) associated with DNA loop formation and sequence-dependent stem hybridization. Specifically, time-correlated single-photon counting experiments on freely diffusing 40A DNA hairpin FRET constructs are performed in a 50 μm × 50 μm square quartz capillary cell pressurized from ambient pressure up to 3 kbar. By pressure-dependent van't Hoff analysis of the equilibrium constants, ΔV° for hybridization of the DNA hairpin can be determined as a function of stem length (n_stem = 7-10) with single base-pair resolution, which further motivates a simple linear deconstruction into additive stem (ΔV°_stem = ΔV°_bp x n_stem) and loop (ΔV°_loop) contributions. We find that increasing pressure destabilizes the DNA hairpin stem region [ΔV°_bp = +1.98(16) cm³/(mol bp)], with additional positive free volume changes [ΔV°_loop = +7.0(14) cm³/mol] we ascribe to bending and base stacking disruption of the 40-dA loop. From a van't Hoff temperature-dependent analysis of the DNA 40A hairpin equilibria, the data support a similar additive loop/stem deconstruction of enthalpic (ΔH° = ΔH°_loop + ΔH°_stem) and entropic (ΔS° = ΔS°_loop + ΔS°_stem) contributions, which permits insightful comparison with predictions from nearest-neighbor thermodynamic models for DNA duplex formation. In particular, the stem thermodynamics is consistent with exothermically favored (ΔH°_stem < 0) and entropically penalized (ΔS°_stem < 0) hydrogen bonding but with additional enthalpic (ΔH°_loop > 0) and entropic (ΔS°_loop > 0) contributions due to loop bending effects consistent with distortion of dA base stacking in the 40-dA linker.