Biophysical Insights from Temperature-Dependent Single-Molecule Förster Resonance Energy Transfer

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.

Kinetic and Thermodynamic Control of G-Quadruplex Polymorphism by Na+ and K+ Cations

G-Quadruplexes (G4s) are ubiquitous nucleic acid folding motifs that exhibit structural diversity that is dependent on cationic conditions. In this work, we exploit temperature-controlled single-molecule fluorescence resonance energy transfer (smFRET) to elucidate the kinetic and thermodynamic mechanisms by which monovalent cations (K+ and Na+) impact folding topologies for a simple G-quadruplex sequence (5'-GGG-(TAAGGG)3-3') with a three-state folding equilibrium. Kinetic measurements indicate that Na+ and K+ influence G4 formation in two distinctly different ways: the presence of Na+ modestly enhances an antiparallel G4 topology through an induced fit (IF) mechanism with a low affinity (Kd = 228 ± 26 mM), while K+ drives G4 into a parallel/hybrid topology via a conformational selection (CS) mechanism with much higher affinity (Kd = 1.9 ± 0.2 mM). Additionally, temperature-dependent studies of folding rate constants and equilibrium ratios reveal distinctly different thermodynamic driving forces behind G4 binding to K+ (ΔH°bind > 0, ΔS°bind > 0) versus Na+ (ΔH°bind < 0, ΔS°bind < 0), which further illuminates the diversity of the possible pathways for monovalent facilitation of G-quadruplex folding.

Entropic control of the free-energy landscape of an archetypal biomolecular machine

Biomolecular machines are complex macromolecular assemblies that utilize thermal and chemical energy to perform essential, multistep, cellular processes. Despite possessing different architectures and functions, an essential feature of the mechanisms of action of all such machines is that they require dynamic rearrangements of structural components. Surprisingly, biomolecular machines generally possess only a limited set of such motions, suggesting that these dynamics must be repurposed to drive different mechanistic steps. Although ligands that interact with these machines are known to drive such repurposing, the physical and structural mechanisms through which ligands achieve this remain unknown. Using temperature-dependent, single-molecule measurements analyzed with a time-resolution-enhancing algorithm, here, we dissect the free-energy landscape of an archetypal biomolecular machine, the bacterial ribosome, to reveal how its dynamics are repurposed to drive distinct steps during ribosome-catalyzed protein synthesis. Specifically, we show that the free-energy landscape of the ribosome encompasses a network of allosterically coupled structural elements that coordinates the motions of these elements. Moreover, we reveal that ribosomal ligands which participate in disparate steps of the protein synthesis pathway repurpose this network by differentially modulating the structural flexibility of the ribosomal complex (i.e., the entropic component of the free-energy landscape). We propose that such ligand-dependent entropic control of free-energy landscapes has evolved as a general strategy through which ligands may regulate the functions of all biomolecular machines. Such entropic control is therefore an important driver in the evolution of naturally occurring biomolecular machines and a critical consideration for the design of synthetic molecular machines.

Ionic Cooperativity between Lysine and Potassium in the Lysine Riboswitch: Single-Molecule Kinetic and Thermodynamic Studies

Functionality in many biological systems, including proteins and nucleic acid structures, including protein and nucleic acid riboswitch structures, can depend on cooperative kinetic behavior between multiple small molecule ligands. In this work, single-molecule FRET data on the Bacillus subtilis lysine riboswitch reveals that affinity for the cognate lysine ligand increases significantly with K⁺, providing evidence for synergism between lysine/K⁺ binding to the aptamer and successful folding of the riboswitch. To describe/interpret this more complex kinetic scenario, we explore the conventional 4-state ("square") model for aptamer binding as a function of K⁺. Extension into this additional dimension generates a novel "cube" model for riboswitch folding dynamics with respect to lysine/K⁺ binding, revealing that riboswitch folding (k_fold) and unfolding (k_unfold) rate constants increase and decrease dramatically with K⁺, respectively. Furthermore, temperature-dependent single-molecule kinetic studies indicate that the presence of K⁺ entropically enhances the transition state barrier to folding but partially compensates for this by increasing the overall exothermicity for lysine binding. We rationalize this behavior as evidence that K⁺ facilitates hydrogen bonding between the negatively charged carboxyl group of lysine and the RNA, increasing structural rigidity and lowering entropy in the binding pocket. Finally, we explore the effects of cation size with Na⁺ and Cs⁺ studies to demonstrate that K⁺ is optimally suited for bridging interactions between lysine and the riboswitch aptamer domain. Regulation of lysine production and transport, dictated by the riboswitch's ability to recognize and bind lysine, is therefore intimately tied to the presence of K⁺ in the binding pocket and is strongly modulated by local cation conditions. The results suggest an increase in lysine riboswitch functionality by sensitivity to additional species in the cellular riboswitch environment.

The biophysics of disordered proteins from the point of view of single-molecule fluorescence spectroscopy

Intrinsically disordered proteins (IDPs) and regions (IDRs) have emerged as key players across many biological functions and diseases. Differently from structured proteins, disordered proteins lack stable structure and are particularly sensitive to changes in the surrounding environment. Investigation of disordered ensembles requires new approaches and concepts for quantifying conformations, dynamics, and interactions. Here, we provide a short description of the fundamental biophysical properties of disordered proteins as understood through the lens of single-molecule fluorescence observations. Single-molecule Förster resonance energy transfer (FRET) and fluorescence correlation spectroscopy (FCS) provides an extensive and versatile toolbox for quantifying the characteristics of conformational distributions and the dynamics of disordered proteins across many different solution conditions, both in vitro and in living cells.

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.

Measuring Excess Heat Capacities of Deoxyribonucleic Acid (DNA) Folding at the Single-Molecule Level

Measurements of the thermodynamic properties of biomolecular folding (ΔG°, ΔH°, ΔS°, etc.) provide a wealth of information on the folding process and have long played a central role in biophysical investigation. In particular, the excess heat capacity of folding (ΔC_P) is crucial, as typically measured in bulk ensemble studies by differential scanning calorimetry (DSC) and isothermal titration calorimetry (ITC). Here, we report the first measurements of ΔC_P at the single-molecule level using the single-molecule fluorescence resonance energy transfer (smFRET) as well as the very first measurements of the heat capacity change associated with achieving the transition state (ΔC^‡_P) for nucleic acid folding. The deoxyribonucleic acid (DNA) hairpin used in these studies exhibits an excess heat capacity for hybridization (ΔC_P = -340 ± 60 J/mol/K per base pair) consistent with the range of literature expectations (ΔC_P = -100 to -420 J/mol/K per base pair). Furthermore, the measured activation heat capacities (ΔC^‡_P) for such hairpin unfolding are consistent with a folding transition state containing few fully formed base pairs, in agreement with prevailing models of DNA hybridization.

Kinetic Microscale Thermophoresis for Simultaneous Measurement of Binding Affinity and Kinetics

Microscale thermophoresis (MST) is a versatile technique to measure binding affinities of binder-ligand systems, based on the directional movement of molecules in a temperature gradient. We extended MST to measure binding kinetics as well as binding affinity in a single experiment by increasing the thermal dissipation of the sample. The kinetic relaxation fingerprints were derived from the fluorescence changes during thermodynamic re-equilibration of the sample after local heating. Using this method, we measured DNA hybridization on-rates and off-rates in the range 10⁴ -10⁶  m^-1  s^-1 and 10^-4 -10^-1  s^-1 , respectively. We observed the expected exponential dependence of the DNA hybridization off-rates on salt concentration, strand length and inverse temperature. The measured on-rates showed a linear dependence on salt concentration and weak dependence on strand length and temperature. For biomolecular interactions with large enthalpic contributions, the kinetic MST technique offers a robust, cost-effective and immobilization-free determination of kinetic rates and binding affinity simultaneously, even in crowded solutions.

Chirality-Dependent Amino Acid Modulation of RNA Folding

The preponderance of a specific d- or l-chirality in fats, sugars, amino acids, nucleic acids, and so on is ubiquitous in nature, yet the biological origin of such chiral dominance (i.e., with one enantiomer overwhelmingly present) remains an open question. One plausible proposal for the predominance of l-chirality in amino acids could be through evolutionary templating of chiral RNA-folding via chaperone activity. To help evaluate this possibility, single molecule fluorescence experiments have been performed that measure the chiral dependence of chaperone folding dynamics for the simple tetraloop-tetraloop receptor (TL-TLR) tertiary binding motif in the presence of a series of chiral amino acids. Specifically, d- vs l-arginine is found to accelerate the unfolding of this RNA motif in a chirally selective fashion, with temperature-dependent studies of the kinetics performed to extract free energy, enthalpy, and entropy landscapes for the underlying thermodynamics. Furthermore, all-atom molecular dynamics (MD) simulations are pursued to provide additional physical insight into this chiral sensitivity, which reveal enantiomer-specific sampling of nucleic acid surfaces by d- vs l-arginine and support a putative mechanism for chirally specific denaturation of RNA tertiary structure by arginine but not other amino acids.

Advances in RNA 3D Structure Modeling Using Experimental Data

RNA is a unique bio-macromolecule that can both record genetic information and perform biological functions in a variety of molecular processes, including transcription, splicing, translation, and even regulating protein function. RNAs adopt specific three-dimensional conformations to enable their functions. Experimental determination of high-resolution RNA structures using x-ray crystallography is both laborious and demands expertise, thus, hindering our comprehension of RNA structural biology. The computational modeling of RNA structure was a milestone in the birth of bioinformatics. Although computational modeling has been greatly improved over the last decade showing many successful cases, the accuracy of such computational modeling is not only length-dependent but also varies according to the complexity of the structure. To increase credibility, various experimental data were integrated into computational modeling. In this review, we summarize the experiments that can be integrated into RNA structure modeling as well as the computational methods based on these experimental data. We also demonstrate how computational modeling can help the experimental determination of RNA structure. We highlight the recent advances in computational modeling which can offer reliable structure models using high-throughput experimental data.

Exploring the effects of cosolutes and crowding on the volumetric and kinetic profile of the conformational dynamics of a poly dA loop DNA hairpin: a single-molecule FRET study

We investigated the volumetric and kinetic profile of the conformational landscape of a poly dA loop DNA hairpin (Hp) in the presence of salts, osmolytes and crowding media, mimicking the intracellular milieu, using single-molecule FRET methodology. Pressure modulation was applied to explore the volumetric and hydrational characteristics of the free-energy landscape of the DNA Hp, but also because pressure is a stress factor many organisms have to cope with, e.g. in the deep sea where pressures even up to the kbar level are encountered. Urea and pressure synergistically destabilize the closed conformation of the DNA Hp due to a lower molar partial volume in the unfolded state. Conversely, multivalent salts, trimethylamine-N-oxide and Ficoll strongly populate the closed state and counteract deteriorating effects of pressure. Complementary smFRET measurements under immobilized conditions at ambient pressure allowed us to dissect the equilibrium data in terms of folding and unfolding rate constants of the conformational transitions, leading to a deeper understanding of the stabilization mechanisms of the cosolutes. Our results show that the free-energy landscape of the DNA Hp is a rugged one, which is markedly affected by the ionic strength of the solution, by preferential interaction and exclusion of cosolvents as well as by pressure.

Raman and Quantum Yield Studies of Trp48-d5 in Azurin: Closed-Shell and Neutral Radical Species

Isotopologues are valuable vibrational probes that shift features in a vibrational spectrum while preserving the electronic structure of the molecule. We report the vibrational and electronic spectra of perdeuterated tryptophan in solution (l-Trp-d₅), as Trp48-d₅ in azurin, and as the photogenerated neutral tryptophan radical, Trp48-d₅^•, in azurin. The UV resonance Raman bands of the perdeuterated closed-shell tryptophan in solution and in azurin are lower in frequency relative to the protiated counterpart. The observed decrease in frequencies of l-Trp-d₅ bands relative to l-Trp-h₅ enables the analysis of vibrational markers of other amino acids, e.g., phenylalanine, that overlap with some modes of l-Trp-h₅. The Raman intensities vary between l-Trp-d₅ and l-Trp-h₅; these differences likely reflect modifications in normal mode composition upon perdeuteration. Analysis of the W3, W6, and W17 modes suggests that the W3 mode retains its utility as a conformational marker; however, the H-bond markers W6 and W17 appear to be less sensitive upon perdeuteration. The neutral tryptophan radical, Trp48-d₅^•, was generated in azurin with a slightly lower radical quantum yield than for Trp48-h₅^•. The visible resonance Raman spectrum of Trp48-d₅^• is different from that of Trp48-h₅^•, especially in terms of relative intensities, and all assignable peaks decreased in frequency upon perdeuteration. The absorption and emission spectra of the perdeuterated closed-shell and radical species exhibited hypsochromic shifts of less than 1 nm relative to the protiated species. The data presented here indicate that l-Trp-d₅ is a valuable probe of vibrational structure, with minimal modification of photoreactivity and photophysics compared to l-Trp-h₅.

Sequence-Dependent Effects of Monovalent Cations on the Structural Dynamics of Trinucleotide-Repeat DNA Hairpins

Repetitive trinucleotide DNA sequences at specific genetic loci are associated with numerous hereditary, neurodegenerative diseases. The propensity of single-stranded domains containing these sequences to form secondary structure via extensive self-complementarity disrupts normal DNA processing to create genetic instabilities. To investigate these intrastrand structural dynamics, a DNA hairpin system was devised for single-molecule fluorescence study of the folding kinetics and energetics for secondary structure formation between two interacting, repetitive domains with specific numbers of the same trinucleotide motif (CXG), where X = T or A. Single-molecule fluorescence resonance energy transfer (smFRET) data show discrete conformational transitions between unstructured and closed hairpin states. The lifetimes of the closed hairpin states correlate with the number of repeats, with (CTG) _N/(CTG) _N domains maintaining longer-lived, closed states than equivalent-sized (CAG) _N/(CAG) _N domains. NaCl promotes similar degree of stabilization for the closed hairpin states of both repeat sequences. Temperature-based, smFRET experiments reveal that NaCl favors hairpin closing for (CAG) _N/(CAG) _N by preordering single-stranded repeat domains to accelerate the closing transition. In contrast, NaCl slows the opening transition of CTG hairpins; however, it promotes misfolded conformations that require unfolding. Energy diagrams illustrate the distinct folding pathways of (CTG) _N and (CAG) _N repeat domains and identify features that may contribute to their gene-destabilizing effects.

Amino Acid Stabilization of Nucleic Acid Secondary Structure: Kinetic Insights from Single-Molecule Studies

Amino acid and nucleic acid interactions are central in biology and may have played a role in the evolutionary development of protein-based life from an early "RNA Universe." To explore the possible role of single amino acids in promoting nucleic acid folding, single-molecule Förster resonance energy transfer experiments have been implemented with a DNA hairpin construct (7 nucleotide double strand with a 40A loop) as a simple model for secondary structure formation. Exposure to positively charged amino acids (arginine and lysine) is found to clearly stabilize the secondary structure. Kinetically, each amino acid promotes folding by generating a large increase in the folding rate with little change in the unfolding rate. From analysis as a function of temperature, arginine and lysine are found to significantly increase the overall exothermicity of folding while imposing only a small entropic penalty on the folding process. Detailed investigations into the kinetics and thermodynamics of this amino acid-induced folding stability reveal arginine and lysine to interact with nucleic acids in a manner reminiscent of monovalent cations. Specifically, these observations are interpreted in the context of an ion atmosphere surrounding the nucleic acid, in which amino acid salts stabilize folding qualitatively like small monovalent cations but also exhibit differences because of the composition of their side chains.

Antagonistic effects of natural osmolyte mixtures and hydrostatic pressure on the conformational dynamics of a DNA hairpin probed at the single-molecule level

Osmolyte mixtures from deep sea organisms are able to rescue nucleic acids from pressure-induced unfolding.

Protein folding transition path times from single molecule FRET

The transition path is the tiny segment of a single molecule trajectory when the free energy barrier between states is crossed and for protein folding contains all of the information about the self-assembly mechanism. As a first step toward obtaining structural information during the transition path from experiments, single molecule FRET spectroscopy has been used to determine average transition path times from a photon-by-photon analysis of fluorescence trajectories. These results, obtained for several different proteins, have already provided new and demanding tests that support both the accuracy of all-atom molecular dynamics simulations and the basic postulates of energy landscape theory of protein folding.