Folding and unfolding single RNA molecules under tension

Inhalation delivery of nucleic acid gene therapies in preclinical drug development

INTRODUCTION

Inhaled gene therapy programs targeting diseases of the lung have seen increasing interest in recent years, though as of yet no product has successfully entered the market. Preclinical research to support such programs is critically important in maximizing the chances of developing successful candidates.

AREAS COVERED

Aspects of inhalation delivery of gene therapies are reviewed, with a focus on preclinical research in animal models. Various barriers to inhalation delivery of gene therapies are discussed, including aerosolization stresses, aerosol behavior in the respiratory tract, and disposition processes post-deposition. Important aspects of animal models are considered, including determinations of biologically relevant determinations of dose and issues related to translatability.

EXPERT OPINION

Development of clinically-efficacious inhaled gene therapies has proven difficult owing to numerous challenges. Fit-for-purpose experimental and analytical methods are necessary for determinations of biologically relevant doses in preclinical animal models. Further developments in disease-specific animal models may aid in improving the translatability of results in future work, and we expect to see accelerated interests in inhalation gene therapies for various diseases. Sponsors, researchers, and regulators are encouraged to engage in early and frequent discussion regarding candidate therapies, and additional dissemination of preclinical methodologies would be of immense value in avoiding common pitfalls.

Single-molecule kinetic locking allows fluorescence-free quantification of protein/nucleic-acid binding

Fluorescence-free micro-manipulation of nucleic acids (NA) allows the functional characterization of DNA/RNA processing proteins, without the interference of labels, but currently fails to detect and quantify their binding. To overcome this limitation, we developed a method based on single-molecule force spectroscopy, called kinetic locking, that allows a direct in vitro visualization of protein binding while avoiding any kind of chemical disturbance of the protein’s natural function. We validate kinetic locking by measuring accurately the hybridization energy of ultrashort nucleotides (5, 6, 7 bases) and use it to measure the dynamical interactions of Escherichia coli/E. coli RecQ helicase with its DNA substrate.

Rieu et al. present a magnetic tweezers based single-molecule manipulation method, called kinetic locking, for direct detection of biomolecular binding without use of fluorescent probes. By measuring dynamical interactions of E. coli RecQ helicase with its DNA substrate, authors show that this method holds promise for studying DNA-DNA and DNA-protein interactions while avoiding the need for labelling.

Mechanical characterization of base analogue modified nucleic acids by force spectroscopy

We use mechanical unfolding of single DNA hairpins with modified bases to accurately assess intra- and intermolecular forces in nucleic acids. As expected, the modification stabilizes the hybridized hairpin, but we also observe intriguing stacking interactions in the unfolded hairpin. Our study highlights the benefit of using base-modified nucleic acids in force-spectroscopy.

Dramatic Shape Changes Occur as Cytochrome c Folds

Extensive experimental studies on the folding of cytochrome c (Cyt c) make this small protein an ideal target for atomic detailed simulations for the purposes of quantitatively characterizing the structural transitions and the associated time scales for folding to the native state from an ensemble of unfolded states. We use previously generated atomically detailed folding trajectories by the stochastic difference equation in length to calculate the time-dependent changes in the small-angle X-ray scattering (SAXS) profiles. Excellent agreement is obtained between experiments and simulations for the time-dependent SAXS spectra, allowing us to identify the structures of the folding intermediates, which shows that Cyt c reaches the native state by a sequential folding mechanism. Using the ensembles of structures along the folding pathways, we show that compaction and the sphericity of Cyt c change dramatically from the prolate ellipsoid shape in the unfolded state to the spherical native state. Our data, which are in unprecedented quantitative agreement with all aspects of time-resolved SAXS experiments, show that hydrophobic collapse and amide group protection coincide on the 100 microseconds time scale, which is in accordance with ultrafast hydrogen/deuterium exchange studies. Based on these results, we propose that compaction of polypeptide chains, accompanied by dramatic shape changes, is a universal characteristic of globular proteins, regardless of the underlying folding mechanism.

Measuring Aptamer Folding Energy Using a Molecular Clamp

Folding energy (ΔG_fold) offers a useful metric for characterizing the stability and function of aptamers. However, experimentally measuring the folding energy is challenging, and there is currently no general technique to measure this parameter directly. In this work, we present a simple approach for measuring aptamer folding energy. First, the aptamer is stretched under equilibrium conditions with a double-stranded DNA "molecular clamp" that is coupled to the aptamer ends. We then measure the total internal energy of stressed DNA molecules using time-lapse gel electrophoresis and compare the folding and unfolding behavior of molecular clamp-stressed molecules that incorporate either the aptamer or unstructured random single-stranded DNA in order to derive the aptamer folding energy. Using this approach, we measured a folding energy of 10.40 kJ/mol for the HD22 thrombin aptamer, which is consistent with other predictions and estimates. We also analyzed a simple hairpin structure, generating a folding energy result of 9.05 kJ/mol, consistent with the value predicted by computational models (9.24 kJ/mol). We believe our strategy offers an accessible and generalizable approach for obtaining such measurements with virtually any aptamer.

The influence of device handles in single-molecule experiments

We deduce a fully analytical model to predict the artifacts of the device handles in single molecule force spectroscopy experiments. As we show, neglecting the handle stiffness can lead to crucial overestimation or underestimation of the stability properties and unfolding thresholds of multistable molecules.

SNARE machinery is optimized for ultrafast fusion

SNARE proteins zipper to form complexes (SNAREpins) that power vesicle fusion with target membranes in a variety of biological processes. A single SNAREpin takes about 1 s to fuse two bilayers, yet a handful can ensure release of neurotransmitters from synaptic vesicles much faster: in a 10th of a millisecond. We propose that, similar to the case of muscle myosins, the ultrafast fusion results from cooperative action of many SNAREpins. The coupling originates from mechanical interactions induced by confining scaffolds. Each SNAREpin is known to have enough energy to overcome the fusion barrier of 25-[Formula: see text]; however, the fusion barrier only becomes relevant when the SNAREpins are nearly completely zippered, and from this state, each SNAREpin can deliver only a small fraction of this energy as mechanical work. Therefore, they have to act cooperatively, and we show that at least three of them are needed to ensure fusion in less than a millisecond. However, to reach the prefusion state collectively, starting from the experimentally observed half-zippered metastable state, the SNAREpins have to mechanically synchronize, which takes more time as the number of SNAREpins increases. Incorporating this somewhat counterintuitive idea in a simple coarse-grained model results in the prediction that there should be an optimum number of SNAREpins for submillisecond fusion: three to six over a wide range of parameters. Interestingly, in situ cryoelectron microscope tomography has very recently shown that exactly six SNAREpins participate in the fusion of each synaptic vesicle. This number is in the range predicted by our theory.

Single-Molecule Mechanical Folding and Unfolding of RNA Hairpins: Effects of Single A-U to A·C Pair Substitutions and Single Proton Binding and Implications for mRNA Structure-Induced -1 Ribosomal Frameshifting

A wobble A·C pair can be protonated at near physiological pH to form a more stable wobble A⁺·C pair. Here, we constructed an RNA hairpin (rHP) and three mutants with one A-U base pair substituted with an A·C mismatch on the top (near the loop, U22C), middle (U25C), and bottom (U29C) positions of the stem, respectively. Our results on single-molecule mechanical (un)folding using optical tweezers reveal the destabilization effect of A-U to A·C pair substitution and protonation-dependent enhancement of mechanical stability facilitated through an increased folding rate, or decreased unfolding rate, or both. Our data show that protonation may occur rapidly upon the formation of an apparent mechanical folding transition state. Furthermore, we measured the bulk -1 ribosomal frameshifting efficiencies of the hairpins by a cell-free translation assay. For the mRNA hairpins studied, -1 frameshifting efficiency correlates with mechanical unfolding force at equilibrium and folding rate at around 15 pN. U29C has a frameshifting efficiency similar to that of rHP (∼2%). Accordingly, the bottom 2-4 base pairs of U29C may not form under a stretching force at pH 7.3, which is consistent with the fact that the bottom base pairs of the hairpins may be disrupted by ribosome at the slippery site. U22C and U25C have a similar frameshifting efficiency (∼1%), indicating that both unfolding and folding rates of an mRNA hairpin in a crowded environment may affect frameshifting. Our data indicate that mechanical (un)folding of RNA hairpins may mimic how mRNAs unfold and fold in the presence of translating ribosomes.

Rare Dissipative Transitions Punctuate the Initiation of Chemical Denaturation in Proteins

Protein unfolding dynamics are bound by their degree of entropy production, a quantity that relates the amount of heat dissipated by a nonequilibrium process to a system's forward and time-reversed trajectories. We here explore the statistics of heat dissipation that emerge in protein molecules subjected to a chemical denaturant. Coupling large molecular dynamics datasets and Markov state models with the theory of entropy production, we demonstrate that dissipative processes can be rigorously characterized over the course of the urea-induced unfolding of the protein chymotrypsin inhibitor 2. By enumerating full entropy production probability distributions as a function of time, we first illustrate that distinct passive and dissipative regimes are present in the denaturation dynamics. Within the dissipative dynamical region, we next find that chymotrypsin inhibitor 2 is strongly driven into unfolded states in which the protein's hydrophobic core has been penetrated by urea molecules and disintegrated. Detailed analyses reveal that urea's interruption of key hydrophobic contacts between core residues causes many of the protein's native structural features to dissolve.

Conformational dynamics of the frameshift stimulatory structure in HIV-1

Programmed ribosomal frameshifting (PRF) in HIV-1 is thought to be stimulated by a hairpin in the mRNA, although a pseudoknot-like triplex has also been proposed. Because the conformational dynamics of the stimulatory structure under tension applied by the ribosomal helicase during translation may play an important role in PRF, we used optical tweezers to apply tension to the HIV stimulatory structure and monitor its unfolding and refolding dynamics. The folding and unfolding kinetics and energy landscape of the hairpin were measured by ramping the force on the hairpin up and down, providing a detailed biophysical characterization. Unexpectedly, whereas unfolding reflected the simple two-state behavior typical of many hairpins, refolding was more complex, displaying significant heterogeneity. Evidence was found for multiple refolding pathways as well as previously unsuspected, partially folded intermediates. Measuring a variant mRNA containing only the sequence required to form the proposed triplex, it behaved largely in the same way. Nonetheless, very rarely, high-force unfolding events characteristic of pseudoknot-like structures were observed. The rare occurrence of the triplex suggests that the hairpin is the functional stimulatory structure. The unusual heterogeneity of the hairpin dynamics under tension suggests a possible functional role in PRF similar to the dynamics of other stimulatory structures.

Direct Observation of Folding Energy Landscape of RNA Hairpin at Mechanical Loading Rates

By applying a controlled mechanical load using optical tweezers, we measured the diffusive barrier crossing in a 49 nt long P5ab RNA hairpin. We find that in the free-energy landscape the barrier height (G^‡) and transition distance (x^‡) are dependent on the loading rate (r) along the pulling direction, x, as predicted by Bell. The barrier shifted toward the initial state, whereas ΔG^‡ reduced significantly from 50 to 5 kT, as r increased from 0 to 32 pN/s. However, the equilibrium work (ΔG) during strand separation, as estimated by Crook's fluctuation theorem, remained unchanged at different rates. Previously, helix formation and denaturation have been described as two-state (F ↔ U) transitions for P5ab. Herein, we report three intermediate states I₁, I, and I₂ located at 4, 11, and 16 nm respectively, from the folded conformation. The intermediates were observed only when the hairpin was subjected to an optimal r, 7.6 pN/s. The results indicate that the complementary strands in P5ab can zip and unzip through complex routes, whereby mismatches act as checkpoints and often impose barriers. The study highlights the significance of loading rates in force-spectroscopy experiments that are increasingly being used to measure the folding properties of biomolecules.

The effect of active fluctuations on the dynamics of particles, motors and DNA-hairpins

Inspired by recent experiments on the dynamics of particles and polymers in artificial cytoskeletons and in cells, we introduce a modified Langevin equation for a particle in an environment that is a viscoelastic medium and that is brought out of equilibrium by the action of active fluctuations caused by molecular motors. We show that within such a model, the motion of a free particle crosses over from superdiffusive to subdiffusive as observed for tracer particles in an in vitro cytoskeleton or in a cell. We investigate the dynamics of a particle confined by a harmonic potential as a simple model for the motion of the tethered head of kinesin-1. We find that the probability that the head is close to its binding site on the microtubule can be enhanced by a factor of two due to active forces. Finally, we study the dynamics of a particle in a double well potential as a model for the dynamics of DNA-hairpins. We show that the active forces effectively lower the potential barrier between the two minima and study the impact of this phenomenon on the zipping/unzipping rate.

Enhanced stochastic fluctuations to measure steep adhesive energy landscapes

Free-energy landscapes govern the behavior of all interactions in the presence of thermal fluctuations in the fields of physical chemistry, materials sciences, and the biological sciences. From the energy landscape, critical information about an interaction, such as the reaction kinetic rates, bond lifetimes, and the presence of intermediate states, can be determined. Despite the importance of energy landscapes to understanding reaction mechanisms, most experiments do not directly measure energy landscapes, particularly for interactions with steep force gradients that lead to premature jump to contact of the probe and insufficient sampling of transition regions. Here we present an atomic force microscopy (AFM) approach for measuring energy landscapes that increases sampling of strongly adhesive interactions by using white-noise excitation to enhance the cantilever's thermal fluctuations. The enhanced fluctuations enable the recording of subtle deviations from a harmonic potential to accurately reconstruct interfacial energy landscapes with steep gradients. Comparing the measured energy landscape with adhesive force measurements reveals the existence of an optimal excitation voltage that enables the cantilever fluctuations to fully sample the shape and depth of the energy surface.

High Spatiotemporal-Resolution Magnetic Tweezers: Calibration and Applications for DNA Dynamics

The observation of biological processes at the molecular scale in real time requires high spatial and temporal resolution. Magnetic tweezers are straightforward to implement, free of radiation or photodamage, and provide ample multiplexing capability, but their spatiotemporal resolution has lagged behind that of other single-molecule manipulation techniques, notably optical tweezers and AFM. Here, we present, to our knowledge, a new high-resolution magnetic tweezers apparatus. We systematically characterize the achievable spatiotemporal resolution for both incoherent and coherent light sources, different types and sizes of beads, and different types and lengths of tethered molecules. Using a bright coherent laser source for illumination and tracking at 6 kHz, we resolve 3 Å steps with a 1 s period for surface-melted beads and 5 Å steps with a 0.5 s period for double-stranded-dsDNA-tethered beads, in good agreement with a model of stochastic bead motion in the magnetic tweezers. We demonstrate how this instrument can be used to monitor the opening and closing of a DNA hairpin on millisecond timescales in real time, together with attendant changes in the hairpin dynamics upon the addition of deoxythymidine triphosphate. Our approach opens up the possibility of observing biological events at submillisecond timescales with subnanometer resolution using camera-based detection.

Statistical mechanics of the Huxley-Simmons model

The chemomechanical model of Huxley and Simmons (HS) [A. F. Huxley and R. M. Simmons, Nature 233, 533 (1971)NATUAS0028-083610.1038/233533a0] provides a paradigmatic description of mechanically induced collective conformational changes relevant in a variety of biological contexts, from muscles power stroke and hair cell gating to integrin binding and hairpin unzipping. We develop a statistical mechanical perspective on the HS model by exploiting a formal analogy with a paramagnetic Ising model. We first study the equilibrium HS model with a finite number of elements and compute explicitly its mechanical and thermal properties. To model kinetics, we derive a master equation and solve it for several loading protocols. The developed formalism is applicable to a broad range of allosteric systems with mean-field interactions.

Connecting thermal and mechanical protein (un)folding landscapes

Molecular dynamics simulations supplement single-molecule pulling experiments by providing the possibility of examining the full free energy landscape using many coordinates. Here, we use an all-atom structure-based model to study the force and temperature dependence of the unfolding of the protein filamin by applying force at both termini. The unfolding time-force relation τ(F) indicates that the force-induced unfolding behavior of filamin can be characterized into three regimes: barrier-limited low- and intermediate-force regimes, and a barrierless high-force regime. Slope changes of τ(F) separate the three regimes. We show that the behavior of τ(F) can be understood from a two-dimensional free energy landscape projected onto the extension X and the fraction of native contacts Q. In the low-force regime, the unfolding rate is roughly force-independent due to the small (even negative) separation in X between the native ensemble and transition state ensemble (TSE). In the intermediate-force regime, force sufficiently separates the TSE from the native ensemble such that τ(F) roughly follows an exponential relation. This regime is typically explored by pulling experiments. While X may fail to resolve the TSE due to overlap with the unfolded ensemble just below the folding temperature, the overlap is minimal at lower temperatures where experiments are likely to be conducted. The TSE becomes increasingly structured with force, whereas the average order of structural events during unfolding remains roughly unchanged. The high-force regime is characterized by barrierless unfolding, and the unfolding time approaches a limit of ∼10 μs for the highest forces we studied. Finally, a combination of X and Q is shown to be a good reaction coordinate for almost the entire force range.

Rigidity generation by nonthermal fluctuations

Active stabilization in systems with zero or negative stiffness is an essential element of a wide variety of biological processes. We study a prototypical example of this phenomenon and show how active rigidity, interpreted as a formation of a pseudowell in the effective energy landscape, can be generated in an overdamped stochastic system. We link the transition from negative to positive rigidity with time correlations in the additive noise, and we show that subtle differences in the out-of-equilibrium driving may compromise the emergence of a pseudowell.

Ligand-induced changes of the apparent transition-state position in mechanical protein unfolding

Force-spectroscopic measurements of ligand-receptor systems and the unfolding/folding of nucleic acids or proteins reveal information on the underlying energy landscape along the pulling coordinate. The slope Δx(‡) of the force-dependent unfolding/unbinding rates is interpreted as the distance from the folded/bound state to the transition state for unfolding/unbinding and, hence, often related to the mechanical compliance of the sample molecule. Here we show that in ligand-binding proteins, the experimentally inferred Δx(‡) can depend on the ligand concentration, unrelated to changes in mechanical compliance. We describe the effect in single-molecule, force-spectroscopy experiments of the calcium-binding protein calmodulin and explain it in a simple model where mechanical unfolding and ligand binding occur on orthogonal reaction coordinates. This model predicts changes in the experimentally inferred Δx(‡), depending on ligand concentration and the associated shift of the dominant barrier between the two reaction coordinates. We demonstrate quantitative agreement between experiments and simulations using a realistic six-state kinetic scheme using literature values for calcium-binding kinetics and affinities. Our results have important consequences for the interpretation of force-spectroscopic data of ligand-binding proteins.

Probing the structural dynamics of proteins and nucleic acids with optical tweezers

Conformational changes are an essential feature of most molecular processes in biology. Optical tweezers have emerged as a powerful tool for probing conformational dynamics at the single-molecule level because of their high resolution and sensitivity, opening new windows on phenomena ranging from folding and ligand binding to enzyme function, molecular machines, and protein aggregation. By measuring conformational changes induced in a molecule by forces applied by optical tweezers, new insight has been gained into the relationship between dynamics and function. We discuss recent advances from studies of how structure forms in proteins and RNA, including non-native structures, fluctuations in disordered proteins, and interactions with chaperones assisting native folding. We also review the development of assays probing the dynamics of complex protein-nucleic acid and protein-protein assemblies that reveal the dynamic interactions between biomolecular machines and their substrates.

Determining intrachain diffusion coefficients for biopolymer dynamics from single-molecule force spectroscopy measurements

The conformational diffusion coefficient for intrachain motions in biopolymers, D, sets the timescale for structural dynamics. Recently, force spectroscopy has been applied to determine D both for unfolded proteins and for the folding transitions in proteins and nucleic acids. However, interpretation of the results remains unsettled. We investigated how instrumental effects arising from the force probes used in the measurement can affect the value of D recovered via force spectroscopy. We compared estimates of D for the folding of DNA hairpins found from measurements of rates and energy landscapes made using optical tweezers with estimates obtained from the same single-molecule trajectories via the transition path time. The apparent D obtained from the rates was much lower than the result found from the same data using transition time analysis, reflecting the effects of the mechanical properties of the force probe. Deconvolution of the finite compliance effects on the measurement allowed the intrinsic value to be recovered. These results were supported by Brownian dynamics simulations of the effects of force-probe compliance and bead size.

Effects of Water on the Single-Chain Elasticity of Poly(U) RNA

Water, the dominant component under the physiological condition, is a complicated solvent which greatly affects the properties of solute molecules. Here, we utilize atomic force microscope-based single-molecule force spectroscopy to study the influence of water on the single-molecule elasticity of an unstructured single-stranded RNA (poly(U)). In nonpolar solvents, RNA presents its inherent elasticity, which is consistent with the theoretical single-chain elasticity calculated by quantum mechanics calculations. In aqueous buffers, however, an additional energy of 1.88 kJ/mol·base is needed for the stretching of the ssRNA chain. This energy is consumed by the bound water rearrangement (Ew) during chain elongation. Further experimental results indicate that the Ew value is uncorrelated to the salt concentrations and stretching velocity. The results obtained in an 8 M guanidine·HCl solution provide more evidence that the bound water molecules around RNA give rise to the observed deviation between aqueous and nonaqueous environments. Compared to synthetic water-soluble polymers, the value of Ew of RNA is much lower. The weak interference of water is supposed to be the precondition for the RNA secondary structure to exist in aqueous solution.

Reconstructing folding energy landscapes from splitting probability analysis of single-molecule trajectories

Structural self-assembly in biopolymers, such as proteins and nucleic acids, involves a diffusive search for the minimum-energy state in a conformational free-energy landscape. The likelihood of folding proceeding to completion, as a function of the reaction coordinate used to monitor the transition, can be described by the splitting probability, p(fold)(x). P(fold) encodes information about the underlying energy landscape, and it is often used to judge the quality of the reaction coordinate. Here, we show how p(fold) can be used to reconstruct energy landscapes from single-molecule folding trajectories, using force spectroscopy measurements of single DNA hairpins. Calculating p(fold)(x) directly from trajectories of the molecular extension measured for hairpins fluctuating in equilibrium between folded and unfolded states, we inverted the result expected from diffusion over a 1D energy landscape to obtain the implied landscape profile. The results agreed well with the landscapes reconstructed by established methods, but, remarkably, without the need to deconvolve instrumental effects on the landscape, such as tether compliance. The same approach was also applied to hairpins with multistate folding pathways. The relative insensitivity of the method to the instrumental compliance was confirmed by simulations of folding measured with different tether stiffnesses. This work confirms that the molecular extension is a good reaction coordinate for these measurements, and validates a powerful yet simple method for reconstructing landscapes from single-molecule trajectories.

A U⋅U Pair-to-U⋅C Pair Mutation-Induced RNA Native Structure Destabilisation and Stretching-Force-Induced RNA Misfolding

Little is known about how a non-Watson-Crick pair affects the RNA folding dynamics. We studied the effects of a U⋅U-to-U⋅C pair mutation on the folding of a hairpin in human telomerase RNA. The ensemble thermal melting of the hairpins shows an on-pathway intermediate with the disruption of the internal loop structure containing the U⋅U/U⋅C pairs. By using optical tweezers, we applied a stretching force on the terminal ends of the hairpins to probe directly the non-nearest-neighbour effects upon the mutations. The single U⋅U to U⋅C mutations are observed to 1) lower the mechanical unfolding force by approximately 1 picoNewton (pN) per mutation without affecting the unfolding reaction transition-state position (thus suggesting that removing a single hydrogen bond affects the structural dynamics at least two base pairs away), 2) result in more frequent misfolding into a small hairpin at approximately 10 pN and 3) shift the folding reaction transition-state position towards the native hairpin structure and slightly increase the mechanical folding kinetics (thus suggesting that untrapping from the misfolded state is not the rate-limiting step).

Transition-Path Probability as a Test of Reaction-Coordinate Quality Reveals DNA Hairpin Folding Is a One-Dimensional Diffusive Process

Chemical reactions are typically described in terms of progress along a reaction coordinate. However, the quality of reaction coordinates for describing reaction dynamics is seldom tested experimentally. We applied a framework for gauging reaction-coordinate quality based on transition-path analysis to experimental data for the first time, looking at folding trajectories of single DNA hairpin molecules measured under tension applied by optical tweezers. The conditional probability for being on a reactive transition path was compared with the probability expected for ideal diffusion over a 1D energy landscape based on the committor function. Analyzing measurements and simulations of hairpin folding where end-to-end extension is the reaction coordinate, after accounting for instrumental effects on the analysis, we found good agreement between transition-path and committor analyses for model two-state hairpins, demonstrating that folding is well-described by 1D diffusion. This work establishes transition-path analysis as a powerful new tool for testing experimental reaction-coordinate quality.

Facile and Cost-Effective Detection of Saxitoxin Exploiting Aptamer Structural Switching

A simple method to detect saxitoxin (STX), one of the main components of the paralytic shellfish poison from red tide, has been developed. By using a next generation dye for double-stranded DNA we were able to differentiate fluorescence from STX-binding aptamers when exposed to different concentrations of STX, suggesting a change in aptamer folding upon target binding. The developed method is extremely rapid, only requiring small sample volumes, with quantitative results in the concentration range of 15 ng/mL to 3 µg/mL of STX, with a detection limit of 7.5 ng/mL.

Reconstructing folding energy landscape profiles from nonequilibrium pulling curves with an inverse Weierstrass integral transform

The energy landscapes that drive structure formation in biopolymers are difficult to measure. Here we validate experimentally a novel method to reconstruct landscape profiles from single-molecule pulling curves using an inverse Weierstrass transform (IWT) of the Jarzysnki free-energy integral. The method was applied to unfolding measurements of a DNA hairpin, replicating the results found by the more-established weighted histogram (WHAM) and inverse Boltzmann methods. Applying both WHAM and IWT methods to reconstruct the folding landscape for a RNA pseudoknot having a stiff energy barrier, we found that landscape features with sharper curvature than the force probe stiffness could not be recovered with the IWT method. The IWT method is thus best for analyzing data from stiff force probes such as atomic force microscopes.

The purine riboswitch as a model system for exploring RNA biology and chemistry

Over the past decade the purine riboswitch, and in particular its nucleobase-binding aptamer domain, has emerged as an important model system for exploring various aspects of RNA structure and function. Its relatively small size, structural simplicity and readily observable activity enable application of a wide variety of experimental approaches towards the study of this RNA. These analyses have yielded important insights into small molecule recognition, co-transcriptional folding and secondary structural switching, and conformational dynamics that serve as a paradigm for other RNAs. In this article, the current state of understanding of the purine riboswitch family and how this growing knowledge base is starting to be exploited in the creation of novel RNA devices are examined. This article is part of a Special Issue entitled: Riboswitches.

Reconstructing folding energy landscapes by single-molecule force spectroscopy

Folding may be described conceptually in terms of trajectories over a landscape of free energies corresponding to different molecular configurations. In practice, energy landscapes can be difficult to measure. Single-molecule force spectroscopy (SMFS), whereby structural changes are monitored in molecules subjected to controlled forces, has emerged as a powerful tool for probing energy landscapes. We summarize methods for reconstructing landscapes from force spectroscopy measurements under both equilibrium and nonequilibrium conditions. Other complementary, but technically less demanding, methods provide a model-dependent characterization of key features of the landscape. Once reconstructed, energy landscapes can be used to study critical folding parameters, such as the characteristic transition times required for structural changes and the effective diffusion coefficient setting the timescale for motions over the landscape. We also discuss issues that complicate measurement and interpretation, including the possibility of multiple states or pathways and the effects of projecting multiple dimensions onto a single coordinate.

Nanomanipulation of single RNA molecules by optical tweezers

A large portion of the human genome is transcribed but not translated. In this post genomic era, regulatory functions of RNA have been shown to be increasingly important. As RNA function often depends on its ability to adopt alternative structures, it is difficult to predict RNA three-dimensional structures directly from sequence. Single-molecule approaches show potentials to solve the problem of RNA structural polymorphism by monitoring molecular structures one molecule at a time. This work presents a method to precisely manipulate the folding and structure of single RNA molecules using optical tweezers. First, methods to synthesize molecules suitable for single-molecule mechanical work are described. Next, various calibration procedures to ensure the proper operations of the optical tweezers are discussed. Next, various experiments are explained. To demonstrate the utility of the technique, results of mechanically unfolding RNA hairpins and a single RNA kissing complex are used as evidence. In these examples, the nanomanipulation technique was used to study folding of each structural domain, including secondary and tertiary, independently. Lastly, the limitations and future applications of the method are discussed.

Single-molecule studies of riboswitch folding

The folding dynamics of riboswitches are central to their ability to modulate gene expression in response to environmental cues. In most cases, a structural competition between the formation of a ligand-binding aptamer and an expression platform (or some other competing off-state) determines the regulatory outcome. Here, we review single-molecule studies of riboswitch folding and function, predominantly carried out using single-molecule FRET or optical trapping approaches. Recent results have supplied new insights into riboswitch folding energy landscapes, the mechanisms of ligand binding, the roles played by divalent ions, the applicability of hierarchical folding models, and kinetic vs. thermodynamic control schemes. We anticipate that future work, based on improved data sets and potentially combining multiple experimental techniques, will enable the development of more complete models for complex RNA folding processes. This article is part of a Special Issue entitled: Riboswitches.

RNA folding pathways and kinetics using 2D energy landscapes

RNA folding pathways play an important role in various biological processes, such as (i) the hok/sok (host-killing/suppression of killing) system in E. coli to check for sufficient plasmid copy number, (ii) the conformational switch in spliced leader (SL) RNA from Leptomonas collosoma, which controls trans splicing of a portion of the '5 exon, and (iii) riboswitches--portions of the 5' untranslated region of messenger RNA that regulate genes by allostery. Since RNA folding pathways are determined by the energy landscape, we describe a novel algorithm, FFTbor2D, which computes the 2D projection of the energy landscape for a given RNA sequence. Given two metastable secondary structures A, B for a given RNA sequence, FFTbor2D computes the Boltzmann probability p(x, y) = Z(x,y)/Z that a secondary structure has base pair distance x from A and distance y from B. Using polynomial interpolationwith the fast Fourier transform,we compute p(x, y) in O(n(5)) time and O(n(2)) space, which is an improvement over an earlier method, which runs in O(n(7)) time and O(n(4)) space. FFTbor2D has potential applications in synthetic biology, where one might wish to design bistable switches having target metastable structures A, B with favorable pathway kinetics. By inverting the transition probability matrix determined from FFTbor2D output, we show that L. collosoma spliced leader RNA has larger mean first passage time from A to B on the 2D energy landscape, than 97.145% of 20,000 sequences, each having metastable structures A, B. Source code and binaries are freely available for download at http://bioinformatics.bc.edu/clotelab/FFTbor2D. The program FFTbor2D is implemented in C++, with optional OpenMP parallelization primitives.

Anti-frameshifting ligand reduces the conformational plasticity of the SARS virus pseudoknot

Programmed -1 ribosomal frameshifting (-1 PRF) stimulated by mRNA pseudoknots regulates gene expression in many viruses, making pseudoknots potential targets for anti-viral drugs. The mechanism by which pseudoknots trigger -1 PRF, however, remains controversial, with several competing models. Recent work showed that high -1 PRF efficiency was linked to high pseudoknot conformational plasticity via the formation of alternate conformers. We tested whether pseudoknots bound with an anti-frameshifting ligand exhibited a similar correlation between conformational plasticity and -1 PRF efficiency by measuring the effects of a ligand that was found to inhibit -1 PRF in the SARS coronavirus on the conformational dynamics of the SARS pseudoknot. Using single-molecule force spectroscopy to unfold pseudoknots mechanically, we found that the ligand binding effectively abolished the formation of alternate conformers. This result extends the connection between -1 PRF and conformational dynamics and, moreover, suggests that targeting the conformational dynamics of pseudoknots may be an effective strategy for anti-viral drug design.

The essential role of stacking adenines in a two-base-pair RNA kissing complex

In minimal RNA kissing complexes formed between hairpins with cognate GACG tetraloops, the two tertiary GC pairs are likely stabilized by the stacking of 5'-unpaired adenines at each end of the short helix. To test this hypothesis, we mutated the flanking adenines to various nucleosides and examined their effects on the kissing interaction. Electrospray ionization mass spectrometry was used to detect kissing dimers in a multiequilibria mixture, whereas optical tweezers were applied to monitor the (un)folding trajectories of single RNA molecules. The experimental findings were rationalized by molecular dynamics simulations. Together, the results showed that the stacked adenines are indispensable for the tertiary interaction. By shielding the tertiary base pairs from solvent and reducing their fraying, the stacked adenines made terminal pairs act more like interior base pairs. The purine double-ring of adenine was essential for effective stacking, whereas additional functional groups modulated the stabilizing effects through varying hydrophobic and electrostatic forces. Furthermore, formation of the kissing complex was dominated by base pairing, whereas its dissociation was significantly influenced by the flanking bases. Together, these findings indicate that unpaired flanking nucleotides play essential roles in the formation of otherwise unstable two-base-pair RNA tertiary interactions.

RNA under tension: Folding Landscapes, Kinetic partitioning Mechanism, and Molecular Tensegrity

Non-coding RNA sequences play a great role in controlling a number of cellular functions, thus raising the need to understand their complex conformational dynamics in quantitative detail. In this perspective, we first show that single molecule pulling when combined with with theory and simulations can be used to quantitatively explore the folding landscape of nucleic acid hairpins, and riboswitches with tertiary interactions. Applications to riboswitches, which are non-coding RNA elements that control gene expression by undergoing dynamical conformational changes in response to binding of metabolites, lead to an organization principle that assembly of RNA is determined by the stability of isolated helices. We also point out the limitations of single molecule pulling experiments, with molecular extension as the only accessible parameter, in extracting key parameters of the folding landscapes of RNA molecules.

Mesoscopic models for DNA stretching under force: New results and comparison with experiments

Single-molecule experiments on double-stranded B-DNA stretching have revealed one or two structural transitions, when increasing the external force. They are characterized by a sudden increase of DNA contour length and a decrease of the bending rigidity. The nature and the critical forces of these transitions depend on DNA base sequence, loading rate, salt conditions and temperature. It has been proposed that the first transition, at forces of 60-80 pN, is a transition from B to S-DNA, viewed as a stretched duplex DNA, while the second one, at stronger forces, is a strand peeling resulting in single-stranded DNAs (ssDNA), similar to thermal denaturation. But due to experimental conditions these two transitions can overlap, for instance for poly(dA-dT). In an attempt to propose a coherent picture compatible with this variety of experimental observations, we derive an analytical formula using a coupled discrete worm-like chain-Ising model. Our model takes into account bending rigidity, discreteness of the chain, linear and non-linear (for ssDNA) bond stretching. In the limit of zero force, this model simplifies into a coupled model already developed by us for studying thermal DNA melting, establishing a connection with previous fitting parameter values for denaturation profiles. Our results are summarized as follows: i) ssDNA is fitted, using an analytical formula, over a nano-Newton range with only three free parameters, the contour length, the bending modulus and the monomer size; ii) a surprisingly good fit on this force range is possible only by choosing a monomer size of 0.2 nm, almost 4 times smaller than the ssDNA nucleobase length; iii) mesoscopic models are not able to fit B to ssDNA (or S to ss) transitions; iv) an analytical formula for fitting B to S transitions is derived in the strong force approximation and for long DNAs, which is in excellent agreement with exact transfer matrix calculations; v) this formula fits perfectly well poly(dG-dC) and λ-DNA force-extension curves with consistent parameter values; vi) a coherent picture, where S to ssDNA transitions are much more sensitive to base-pair sequence than the B to S one, emerges. This relatively simple model might allow one to further study quantitatively the influence of salt concentration and base-pairing interactions on DNA force-induced transitions.

Programmed -1 frameshifting efficiency correlates with RNA pseudoknot conformational plasticity, not resistance to mechanical unfolding

Programmed -1 frameshifting, whereby the reading frame of a ribosome on messenger RNA is shifted in order to generate an alternate gene product, is often triggered by a pseudoknot structure in the mRNA in combination with an upstream slippery sequence. The efficiency of frameshifting varies widely for different sites, but the factors that determine frameshifting efficiency are not yet fully understood. Previous work has suggested that frameshifting efficiency is related to the resistance of the pseudoknot against mechanical unfolding. We tested this hypothesis by studying the mechanical properties of a panel of pseudoknots with frameshifting efficiencies ranging from 2% to 30%: four pseudoknots from retroviruses, two from luteoviruses, one from a coronavirus, and a nonframeshifting bacteriophage pseudoknot. Using optical tweezers to apply tension across the RNA, we measured the distribution of forces required to unfold each pseudoknot. We found that neither the average unfolding force, nor the unfolding kinetics, nor the parameters describing the energy landscape for mechanical unfolding of the pseudoknot (energy barrier height and distance to the transition state) could be correlated to frameshifting efficiency. These results indicate that the resistance of pseudoknots to mechanical unfolding is not a primary determinant of frameshifting efficiency. However, increased frameshifting efficiency was correlated with an increased tendency to form alternate, incompletely folded structures, suggesting a more complex picture of the role of the pseudoknot involving the conformational dynamics.

Transition path times for nucleic Acid folding determined from energy-landscape analysis of single-molecule trajectories

The duration of structural transitions in biopolymers is only a fraction of the time spent searching diffusively over the configurational energy landscape. We found the transition time, τ(TP), and the diffusion constant, D, for DNA and RNA folding using energy landscapes obtained from single-molecule trajectories under tension in optical traps. DNA hairpins, RNA pseudoknots, and a riboswitch all had τ(TP)~10 μs and D~10(-13-14) m(2)/s, despite widely differing unfolding rates. These results show how energy-landscape analysis can be harnessed to characterize brief but critical events during folding reactions.

Structure and mechanism of purine-binding riboswitches

A riboswitch is a non-protein coding sequence capable of directly binding a small molecule effector without the assistance of accessory proteins to regulate expression of the mRNA in which it is embedded. Currently, over 20 different classes of riboswitches have been validated in bacteria with the promise of many more to come, making them an important means of regulating the genome in the bacterial kingdom. Strikingly, half of the known riboswitches recognize effector compounds that contain a purine or related moiety. In the last decade, significant progress has been made to determine how riboswitches specifically recognize these compounds against the background of many other similar cellular metabolites and transduce this signal into a regulatory response. Of the known riboswitches, the purine family containing guanine, adenine and 2'-deoxyguanosine-binding classes are the most extensively studied, serving as a simple and useful paradigm for understanding how these regulatory RNAs function. This review provides a comprehensive summary of the current state of knowledge regarding the structure and mechanism of these riboswitches, as well as insights into how they might be exploited as therapeutic targets and novel biosensors.

Enzymatic oxygen scavenging for photostability without pH drop in single-molecule experiments

Over the past years, bottom-up bionanotechnology has been developed as a promising tool for future technological applications. Many of these biomolecule-based assemblies are characterized using various single-molecule techniques that require strict anaerobic conditions. The most common oxygen scavengers for single-molecule experiments are glucose oxidase and catalase (GOC) or protocatechuate dioxygenase (PCD). One of the pitfalls of these systems, however, is the production of carboxylic acids. These acids can result in a significant pH drop over the course of experiments and must thus be compensated by an increased buffer strength. Here, we present pyranose oxidase and catalase (POC) as a novel enzymatic system to perform single-molecule experiments in pH-stable conditions at arbitrary buffer strength. We show that POC keeps the pH stable over hours, while GOC and PCD cause an increasing acidity of the buffer system. We further verify in single-molecule fluorescence experiments that POC performs as good as the common oxygen-scavenging systems, but offers long-term pH stability and more freedom in buffer conditions. This enhanced stability allows the observation of bionanotechnological assemblies in aqueous environments under well-defined conditions for an extended time.

Single-molecule measurements of the binding between small molecules and DNA aptamers

Aptamers that bind small molecules can serve as basic biosensing platforms. Evaluation of the binding constant between an aptamer and a small molecule helps to determine the effectiveness of the aptamer-based sensors. Binding constants are often measured by a series of experiments with varying ligand or aptamer concentrations. Such experiments are time-consuming, material nonprudent, and prone to low reproducibility. Here, we use laser tweezers to determine the dissociation constant for aptamer-ligand interactions at the single-molecule level from only one ligand concentration. Using an adenosine 5'-triphosphate disodium salt (ATP) binding aptamer as an example, we have observed that the mechanical stabilities of aptamers bound with ATP are higher than those without a ligand. Comparison of the change in free energy of unfolding (ΔG(unfold)) between these two aptamers yields a ΔG of 33 ± 4 kJ/mol for the binding. By applying a Hess-like cycle at room temperature, we obtained a dissociation constant (K(d)) of 2.0 ± 0.2 μM, a value consistent with the K(d) obtained from our equilibrated capillary electrophoresis (CE) (2.4 ± 0.4 μM) and close to that determined by affinity chromatography in the literature (6 ± 3 μM). We anticipate that our laser tweezers and CE methodologies may be used to more conveniently evaluate the binding between receptors and ligands and also serve as analytical tools for force-based biosensing.

Single-molecule analysis of telomerase structure and function

The telomerase ribonucleoprotein is a specialized reverse transcriptase required to maintain protective chromosome end-capping structures called telomeres. In most cells, telomerase is not active and the natural shortening of telomeres with each round of DNA replication ultimately triggers cell growth arrest. In contrast, the presence of telomerase confers a high level of renewal capacity upon rapidly dividing cells. Telomerase is aberrantly activated in 90% of human cancers and thus represents an important target for anticancer therapeutics. However, the naturally low abundance of telomerase has hampered efforts to obtain high-resolution models for telomerase structure and function. To circumvent these challenges, single-molecule techniques have recently been employed to investigate telomerase assembly, structure, and catalysis.

Single-molecule force spectroscopy of the add adenine riboswitch relates folding to regulatory mechanism

Riboswitches regulate gene expression via ligand binding to an aptamer domain which induces conformational changes in a regulatory expression platform. By unfolding and refolding single add adenine riboswitch molecules in an optical trap, an integrated picture of the folding was developed and related to the regulatory mechanism. Force-extension curves (FECs) and constant-force folding trajectories measured on the aptamer alone revealed multiple partially-folded states, including several misfolded states not on the native folding pathway. All states were correlated to key structural components and interactions within hierarchical folding pathways. FECs of the full-length riboswitch revealed that the thermodynamically stable conformation switches upon ligand binding from a structure repressing translation to one permitting it. Along with rapid equilibration of the two structures in the absence of adenine, these results support a thermodynamically-controlled regulatory mechanism, in contrast with the kinetic control of the closely-related pbuE adenine riboswitch. Comparison of the folding of these riboswitches revealed many similarities arising from shared structural features but also essential differences related to their different regulatory mechanisms.

Free-energy landscapes and thermodynamic parameters of complex molecules from nonequilibrium simulation trajectories

Thermodynamic parameters such as free energies and heat capacities are important quantities for understanding processes involving structural transitions in complex molecules such as proteins. Computational investigations provide simulated data that can be used for calculating thermodynamic parameters. However, calculations give accurate results only if the simulations sample all of configuration space with the appropriate temperature-dependent Boltzmann equilibrium probabilities. For many systems, truly comprehensive sampling of configuration space is not computationally feasible. We present an approximation technique for the calculations that will give accurate values for thermodynamic parameters when the data is incomplete. Our work is applicable to systems in which there are two distinct, important regions of configuration space that must be sampled. Importantly, the results are also valid when the system is more complex than two-state systems. Transition pathways that involve intermediate configurations between two stable regions are allowed in this treatment, and therefore the results are valid for multistate systems.

HU protein induces incoherent DNA persistence length

HU is a highly conserved protein that is believed to play an important role in the architecture and dynamic compaction of bacterial DNA. Its ability to control DNA bending is crucial for functions such as transcription and replication. The effects of HU on the DNA structure have been studied so far mainly by single molecule methods that require us to apply stretching forces on the DNA and therefore may perturb the DNA-protein interaction. To overcome this hurdle, we study the effect of HU on the DNA structure without applying external forces by using an improved tethered particle motion method. By combining the results with DNA curvature analysis from atomic force microscopy measurements we find that the DNA consists of two different curvature distributions and the measured persistence length is determined by their interplay. As a result, the effective persistence length adopts a bimodal property that depends primarily on the HU concentration. The results can be explained according to a recently suggested model that distinguishes single protein binding from cooperative protein binding.