Single-molecule mechanical unfolding and folding of a pseudoknot in human telomerase RNA

Deciphering the sequence-dependent unfolding pathways of an RNA pseudoknot with steered molecular dynamics

Programmed ribosomal frameshifting in Simian retrovirus-1 (SRV-1) is sensitive to the mechanical properties of an RNA pseudoknot. Unravelling these mechanical intricacies via unfolding reveals fundamental insights into their structural dynamics. Using constant velocity steered molecular dynamics (CV-SMD) simulations, we explored the unfolding dynamics and the impact of mutations on the unfolding pathway of the pseudoknot. Except for A28C, A/U to C mutations that disrupt base triples between the loop 2 and stem 1 significantly weaken the pseudoknot and make it more susceptible to unfolding. Complementary mutations in 3 base pairs of the stem region (S1) enhanced its susceptibility to disruption except for Mut5 (S2). We quantitatively assessed the variations in unfolding pathways by analysing the opening of distinct Canonical (WC) and non-canonical (NWC) interactions, force-extension curves, and potential mean force profiles (as a guiding decision for planning mutations). These findings offer a quantified perspective, showcasing the potential of utilizing the unfolding pathways of RNA pseudoknots to explore the programmability of RNA structures. This insight proves valuable for designing RNA-PROTACS and RNA-aptamers, allowing for the assessment and manipulation of their biological folding/unfolding processes.

FAIR data for optical tweezers experiments

The single-molecule biophysics community has delivered significant impacts to our understanding of fundamental biological processes, yet the field is also siloed and has fragmented data structures, which impede data sharing and limit the ability to conduct comprehensive meta-analyses. To advance the field of optical tweezers in single-molecule biophysics, it is important that the field adopts open and collaborative data sharing that facilitate meta-analyses that combine diverse resources and supports more advanced analyses, akin to those seen in projects such as the Protein Data Bank and the 1000 Genomes Project. Here, we assess the state of data findability, accessibility, interoperability, and reusability (the FAIR principles) within the single-molecule optical tweezers field. By combining a qualitative review with quantitative tools from bibliometrics, our analysis suggests that the field has significant room for improvement in terms of FAIR adherence. Finally, we discuss the potential of compulsory data deposition and a minimal set of metadata standards to ensure reproducibility and interoperability between systems. While implementing these measures may not be straightforward, they are key steps that will enhance the integration of optical tweezers biophysics with the broader biomedical literature.

Exploring the conformational dynamics of the SARS-CoV-2 SL4 hairpin by combining optical tweezers and base analogues

The parasitic nature of the SARS-CoV-2 virus demands selective packaging of its RNA genome (gRNA) from the abundance of other nucleic acids present in infected cells. Despite increasing evidence that stem-loop 4 (SL4) of the gRNA 5' UTR is involved in the initiation of this process by binding the nucleocapsid (N) protein, little is known about its conformational dynamics. Here, we unravel the stability, dynamics and (un)folding pathways of SL4 using optical tweezers and a base analogue, tCO, that provides a local and subtle increase in base stacking without perturbing hydrogen bonding. We find that SL4 (un)folds mainly in a single step or through an intermediate, encompassing nucleotides from the central U bulge to the hairpin loop. Due to an upper-stem CU mismatch, SL4 is prone to misfold, the extent of which can be tuned by incorporating tCO at different positions. Our study contributes to a better understanding of SARS-CoV-2 packaging and the design of drugs targeting SL4. We also highlight the generalizability of using base analogues in optical tweezers experiments for probing intramolecular states and conformational transitions of various nucleic acids at the level of single molecules and with base-pair resolution.

A versatile and high-throughput flow-cell system combined with fluorescence imaging for simultaneous single-molecule force measurement and visualization

A flow-cell offers many advantages for single-molecule studies. But, its merit as a quantitative single-molecule tool has long been underestimated. In this work, we developed a gas-pumped fully calibrated flow-cell system combined with fluorescence imaging for simultaneous single-molecule force measurement and visualization. Such a flow-cell system has considered the hydrodynamic drags on biomolecules and hence can apply and measure force up to more than 100 pN in sub-pN precision with an ultra-high force stability (force drift <0.01 pN in 10 minutes) and tuning accuracy (∼0.04 pN). Meanwhile, it also allows acquiring force signals and fluorescence images at the same time, parallelly tracking hundreds of protein motors in real time as well as monitoring the conformational changes of biomolecules under a well-controlled force, as demonstrated by a series of single-molecule experiments in this work, including the studies of DNA overstretching dynamics, transcription under force and DNA folding/unfolding dynamics. Interesting findings, such as the very tight association of single-stranded binding (SSB) proteins with ssDNA and the reversed transcription, have also been made. These results together lay down an essential foundation for a flow-cell to be used as a versatile, quantitative and high-throughput tool for single-molecule manipulation and visualization.

Stability of DNA and RNA hairpins: a comparative study based on ox-DNA

Advances in single-molecule experiments on macromolecular crowding urgently need an efficient simulation method to resolve their discrepancies quantitatively. Ox-DNA model has been since reworked to treat the thermodynamics and mechanical properties of DNA/RNA hairpin at a stretching force. In hopping experiments, the critical forces of RNA hairpins at different temperatures are greater than those of DNA hairpins, in addition, the Gibbs free energy at a fixed temperature required to convert an RNA hairpin into a single-stranded molecule at zero force is obviously greater than that of DNA hairpin and gradually decreases by increasing the temperature. As far as force-ramping experiments are concerned, the first-rupture forces of RNA/DNA hairpins corresponding to the maximum probability density linearly pertain to the force-loading rate, with those of RNA hairpins being greater. The extended ox-DNA model could potentially identify the interaction between biologically inert polymer and RNA/DNA hairpins in crowded environments.

An overview of single-molecule techniques and applications in the study of nucleic acid structure and function

Nucleic acids are an indispensable component in all known life forms. The biological processes are regulated by Nucleic acids, which associate to form special high-order structures. since the high-level structures of nucleic acids are related to gene expression in cancer cells or viruses, it is very likely to become a potential drug target. Traditional biochemical methods are limited to distinguish the conformational distribution and dynamic transition process of single nucleic acid structure. The ligands based on the intermediate and transition states between different conformations are not designed by traditional biochemical methods. The single-molecule techniques enable real-time observation of the individual nucleic acid behavior due to its high resolution. Here, we introduce the application of single-molecule techniques in the study of small molecules to recognize nucleic acid structures, such as single-molecule FRET, magnetic tweezers, optical tweezers and atomic force microscopy. At the same time, we also introduce the specific advantages of single-molecule technology compared with traditional biochemical methods and some problems arisen in current research.

Optical tweezers for drug discovery

The time taken and the cost of producing novel therapeutic drugs presents a significant burden - a typical target-based drug discovery process involves computational screening of drug libraries, compound assays and expensive clinical trials. This review summarises the value of dynamic conformational information obtained by optical tweezers and how this information can target 'undruggable' proteins. Optical tweezers provide insights into the link between biological mechanisms and structural conformations, which can be used in drug discovery. Developing workflows including software and sample preparation will improve throughput, enabling adoption of optical tweezers in biopharma. As a complementary tool, optical tweezers increase the number of drug candidates, improve the understanding of a target's complex structural dynamics and elucidate interactions between compounds and their targets.

Single-molecule fluorescence imaging techniques reveal molecular mechanisms underlying deoxyribonucleic acid damage repair

Advances in single-molecule techniques have uncovered numerous biological secrets that cannot be disclosed by traditional methods. Among a variety of single-molecule methods, single-molecule fluorescence imaging techniques enable real-time visualization of biomolecular interactions and have allowed the accumulation of convincing evidence. These techniques have been broadly utilized for studying DNA metabolic events such as replication, transcription, and DNA repair, which are fundamental biological reactions. In particular, DNA repair has received much attention because it maintains genomic integrity and is associated with diverse human diseases. In this review, we introduce representative single-molecule fluorescence imaging techniques and survey how each technique has been employed for investigating the detailed mechanisms underlying DNA repair pathways. In addition, we briefly show how live-cell imaging at the single-molecule level contributes to understanding DNA repair processes inside cells.

Applicability of the nearest-neighbour model for pseudoknot RNAs

The validity of the nearest-neighbour (NN) model was verified in an RNA pseudoknot (PK) structure. The thermodynamic parameters of the second hairpin stem (S2) region, which separates the PK from a hairpin structure, were monitored using CD and UV melting. Different PKs with identical NN base pairs in the S2 region exhibited similar thermodynamic parameters, highlighting the validity of the NN model in this RNA tertiary structure motif.

Force-Dependent Folding Kinetics of Single Molecules with Multiple Intermediates and Pathways

Most single-molecule studies derive the kinetic rates of native, intermediate, and unfolded states from equilibrium hopping experiments. Here, we apply the Kramers kinetic diffusive model to derive the force-dependent kinetic rates of intermediate states from nonequilibrium pulling experiments. From the kinetic rates, we also extract the force-dependent kinetic barriers and the equilibrium folding energies. We apply our method to DNA hairpins with multiple folding pathways and intermediates. The experimental results agree with theoretical predictions. Furthermore, the proposed nonequilibrium single-molecule approach permits us to characterize kinetic and thermodynamic properties of native, unfolded, and intermediate states that cannot be derived from equilibrium hopping experiments.

Structural and molecular basis for Cardiovirus 2A protein as a viral gene expression switch

Programmed –1 ribosomal frameshifting (PRF) in cardioviruses is activated by the 2A protein, a multi-functional virulence factor that also inhibits cap-dependent translational initiation. Here we present the X-ray crystal structure of 2A and show that it selectively binds to a pseudoknot-like conformation of the PRF stimulatory RNA element in the viral genome. Using optical tweezers, we demonstrate that 2A stabilises this RNA element, likely explaining the increase in PRF efficiency in the presence of 2A. Next, we demonstrate a strong interaction between 2A and the small ribosomal subunit and present a cryo-EM structure of 2A bound to initiated 70S ribosomes. Multiple copies of 2A bind to the 16S rRNA where they may compete for binding with initiation and elongation factors. Together, these results define the structural basis for RNA recognition by 2A, show how 2A-mediated stabilisation of an RNA pseudoknot promotes PRF, and reveal how 2A accumulation may shut down translation during virus infection.

Many RNA viruses employ programmed –1 ribosomal frameshifting (PRF) to expand their coding capacity and optimize production of viral proteins. Here, the authors report structural and biophysical analysis of protein 2A from a cardiovirus, with insights into the mechanism of its PRF-stimulatory function.

Complex Conformational Dynamics of the Heart Failure-Associated Pre-miRNA-377 Hairpin Revealed by Single-Molecule Optical Tweezers

Pre-miRNA-377 is a hairpin-shaped regulatory RNA associated with heart failure. Here, we use single-molecule optical tweezers to unzip pre-miRNA-377 and study its stability and dynamics. We show that magnesium ions have a strong stabilizing effect, and that sodium ions stabilize the hairpin more than potassium ions. The hairpin unfolds in a single step, regardless of buffer composition. Interestingly, hairpin folding occurs either in a single step (type 1) or through the formation of intermediates, in multiple steps (type 2) or gradually (type 3). Type 3 occurs only in the presence of both sodium and magnesium, while type 1 and 2 take place in all buffers, with type 1 being the most prevalent. By reducing the size of the native hairpin loop from fourteen to four nucleotides, we demonstrate that the folding heterogeneity originates from the large size of the hairpin loop. Further, while efficient pre-miRNA-377 binders are lacking, we demonstrate that the recently developed C2 ligand displays bimodal activity: it enhances the mechanical stability of the pre-miRNA-377 hairpin and perturbs its folding. The knowledge regarding pre-miRNA stability and dynamics that we provide is important in understanding its regulatory function and how it can be modulated to achieve a therapeutic effect, e.g., in heart failure treatment.

Structural dynamics of single SARS-CoV-2 pseudoknot molecules reveal topologically distinct conformers

The RNA pseudoknot that stimulates programmed ribosomal frameshifting in SARS-CoV-2 is a possible drug target. To understand how it responds to mechanical tension applied by ribosomes, thought to play a key role during frameshifting, we probe its structural dynamics using optical tweezers. We find that it forms multiple structures: two pseudoknotted conformers with different stability and barriers, and alternative stem-loop structures. The pseudoknotted conformers have distinct topologies, one threading the 5′ end through a 3-helix junction to create a knot-like fold, the other with unthreaded 5′ end, consistent with structures observed via cryo-EM and simulations. Refolding of the pseudoknotted conformers starts with stem 1, followed by stem 3 and lastly stem 2; Mg²⁺ ions are not required, but increase pseudoknot mechanical rigidity and favor formation of the knot-like conformer. These results resolve the SARS-CoV-2 frameshift signal folding mechanism and highlight its conformational heterogeneity, with important implications for structure-based drug-discovery efforts.

The RNA pseudoknot of SARS-CoV-2 promotes -1 programmed ribosomal frameshifting. Here the authors use single molecule force spectroscopy to study the folding of this pseudoknot, showing that it forms at least two different pseudoknot conformers with distinct fold topologies.

Modeling the structure of the frameshift-stimulatory pseudoknot in SARS-CoV-2 reveals multiple possible conformers

The coronavirus causing the COVID-19 pandemic, SARS-CoV-2, uses -1 programmed ribosomal frameshifting (-1 PRF) to control the relative expression of viral proteins. As modulating -1 PRF can inhibit viral replication, the RNA pseudoknot stimulating -1 PRF may be a fruitful target for therapeutics treating COVID-19. We modeled the unusual 3-stem structure of the stimulatory pseudoknot of SARS-CoV-2 computationally, using multiple blind structural prediction tools followed by μs-long molecular dynamics simulations. The results were compared for consistency with nuclease-protection assays and single-molecule force spectroscopy measurements of the SARS-CoV-1 pseudoknot, to determine the most likely conformations. We found several possible conformations for the SARS-CoV-2 pseudoknot, all having an extended stem 3 but with different packing of stems 1 and 2. Several conformations featured rarely-seen threading of a single strand through junctions formed between two helices. These structural models may help interpret future experiments and support efforts to discover ligands inhibiting -1 PRF in SARS-CoV-2.

Tertiary Base Triple Formation in the SRV-1 Frameshifting Pseudoknot Stabilizes Secondary Structure Components

Minor-groove base triples formed between stem 1 and loop 2 of the simian retrovirus type 1 (SRV-1) mRNA frameshifting pseudoknot are essential in stimulating -1 ribosomal frameshifting. How tertiary base triple formation affects the local stabilities of secondary structures (stem 1 and stem 2) and thus ribosomal frameshifting efficiency is not well understood. We made a short peptide nucleic acid (PNA) that is expected to invade stem 1 of the SRV-1 pseudoknot by PNA-RNA duplex formation to mimic the stem 1 unwinding process by a translating ribosome. In addition, we used a PNA for invading stem 2 in the SRV-1 pseudoknot. Our nondenaturing polyacrylamide gel electrophoresis data for the binding of PNA to the SRV-1 pseudoknot and mutants reveal that mutations in loop 2 disrupting base triple formation between loop 2 and stem 1 in the SRV-1 pseudoknot result in enhanced invasion by both PNAs. Our data suggest that tertiary stem 1-loop 2 base triple interactions in the SRV-1 pseudoknot can stabilize both of the secondary structural components, stem 1 and stem 2. Stem 2 stability is thus coupled to the structural stability of stem 1-loop 2 base triples, mediated through a long-range effect. The apparent dissociation constants of both PNAs are positively correlated with the pseudoknot mechanical stabilities and frameshifting efficiencies. The relatively simple PNA local invasion experiment may be used to characterize the energetic contribution of tertiary interactions and ligand binding in many other RNA and DNA structures.

Cooperative Analysis of Structural Dynamics in RNA-Protein Complexes by Single-Molecule Förster Resonance Energy Transfer Spectroscopy

RNA-protein complexes (RNPs) are essential components in a variety of cellular processes, and oftentimes exhibit complex structures and show mechanisms that are highly dynamic in conformation and structure. However, biochemical and structural biology approaches are mostly not able to fully elucidate the structurally and especially conformationally dynamic and heterogeneous nature of these RNPs, to which end single molecule Förster resonance energy transfer (smFRET) spectroscopy can be harnessed to fill this gap. Here we summarize the advantages of strategic smFRET studies to investigate RNP dynamics, complemented by structural and biochemical data. Focusing on recent smFRET studies of three essential biological systems, we demonstrate that investigation of RNPs on a single molecule level can answer important functional questions that remained elusive with structural or biochemical approaches alone: The complex structural rearrangements throughout the splicing cycle, unwinding dynamics of the G-quadruplex (G4) helicase RHAU, and aspects in telomere maintenance regulation and synthesis.

The Role of Proton Transfer on Mutations

Hydrogen bonds play a critical role in nucleobase studies as they encode genes, map protein structures, provide stability to the base pairs, and are involved in spontaneous and induced mutations. Proton transfer mechanism is a critical phenomenon that is related to the acid-base characteristics of the nucleobases in Watson-Crick base pairs. The energetic and dynamical behavior of the proton can be depicted from these characteristics and their adjustment to the water molecules or the surrounding ions. Further, new pathways open up in which protonated nucleobases are generated by proton transfer from the ionized water molecules and elimination of a hydroxyl radical in this review, the analysis will be focused on understanding the mechanism of untargeted mutations in canonical, wobble, Hoogsteen pairs, and mutagenic tautomers through the non-covalent interactions. Further, rare tautomer formation through the single proton transfer (SPT) and the double proton transfer (DPT), quantum tunneling in nucleobases, radiation-induced bystander effects, role of water in proton transfer (PT) reactions, PT in anticancer drugs-DNA interaction, displacement and oriental polarization, possible models for mutations in DNA, genome instability, and role of proton transfer using kinetic parameters for RNA will be discussed.

Complex dynamics under tension in a high-efficiency frameshift stimulatory structure

Specific structures in mRNA can stimulate programmed ribosomal frameshifting (PRF). PRF efficiency can vary enormously between different stimulatory structures, but the features that lead to efficient PRF stimulation remain uncertain. To address this question, we studied the structural dynamics of the frameshift signal from West Nile virus (WNV), which stimulates -1 PRF at very high levels and has been proposed to form several different structures, including mutually incompatible pseudoknots and a double hairpin. Using optical tweezers to apply tension to single mRNA molecules, mimicking the tension applied by the ribosome during PRF, we found that the WNV frameshift signal formed an unusually large number of different metastable structures, including all of those previously proposed. From force-extension curve measurements, we mapped 2 mutually exclusive pathways for the folding, each encompassing multiple intermediates. We identified the intermediates in each pathway from length changes and the effects of antisense oligomers blocking formation of specific contacts. Intriguingly, the number of transitions between the different conformers of the WNV frameshift signal was maximal in the range of forces applied by the ribosome during -1 PRF. Furthermore, the occupancy of the pseudoknotted conformations was far too low for static pseudoknots to account for the high levels of -1 PRF. These results support the hypothesis that conformational heterogeneity plays a key role in frameshifting and suggest that transitions between different conformers under tension are linked to efficient PRF stimulation.

Monovalent ions modulate the flux through multiple folding pathways of an RNA pseudoknot

The assembly mechanism of RNA, vital to describing its functions, depends on both the sequence and the metal ion concentration. How the latter influences the folding trajectories remains an important unsolved problem. Here, we examine the folding pathways of an RNA pseudoknot (PK) with key functional roles in transcription and translation, using a combination of experiments and simulations. We demonstrate that the PK, consisting of two hairpins with differing stabilities, folds by parallel pathways. Surprisingly, the flux between them is modulated by monovalent salt concentration. Our work shows that the order of assembly of PKs is determined by the relative stability of the hairpins, implying that the folding landscape can be controlled by sequence and ion concentration.

The functions of RNA pseudoknots (PKs), which are minimal tertiary structural motifs and an integral part of several ribozymes and ribonucleoprotein complexes, are determined by their structure, stability, and dynamics. Therefore, it is important to elucidate the general principles governing their thermodynamics/folding mechanisms. Here, we combine laser temperature-jump experiments and coarse-grained simulations to determine the folding/unfolding pathways of VPK, a variant of the mouse mammary tumor virus (MMTV) PK involved in ribosomal frameshifting. Fluorescent nucleotide analogs (2-aminopurine and pyrrolocytidine) placed at different stem/loop positions in the PK serve as local probes allowing us to monitor the order of assembly of VPK that has two constituent hairpins with different intrinsic stabilities. We show that at 50 mM KCl, the dominant folding pathway populates only the more stable hairpin intermediate; as the salt concentration is increased, a parallel folding pathway emerges involving the less stable hairpin as an alternate intermediate. Notably, the flux between the pathways is modulated by the ionic strength. Our findings support the principle that the order of PK structure formation is determined by the relative stabilities of the hairpins, which can be altered by sequence variations or salt concentrations. The experimental results of salt effects on the partitioning between the two folding pathways are in remarkable agreement with simulations that were performed with no adjustable parameters. Our study not only unambiguously demonstrates that VPK folds by parallel pathways but also showcases the power of combining experiments and simulations for a more enriched description of RNA self-assembly.

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.

Combined single molecule experimental and computational approaches for understanding the unfolding pathway of a viral translation enhancer that participates in a conformational switch

How plus-strand [+]RNA virus genomes transition from translation templates to replication templates is a matter of much speculation. We have previously proposed that, for Turnip crinkle virus, binding of the encoded RNA-dependent RNA polymerase (RdRp) to the 3'UTR of the [+]RNA template promotes a regional wide-spread conformational switch to an alternative structure that disassembles the cap-independent translation enhancer (CITE) in the 3'UTR. The active 3'CITE folds into a tRNA-like T-shaped structure (TSS) that binds to 80S ribosomes and 60S subunits in the P-site. In this Point-of-View, we discuss the history of our research on the TSS and our recent report combining coarse level single molecule force spectroscopy (optical tweezers) with fine-grain computer simulations of this experimental process and biochemical approaches to obtain a detailed understanding of how RdRp binding in the TSS vicinity might lead to an extensive rearrangement of the RNA structure.

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.

Single-Molecule Studies of Telomeres and Telomerase

Telomeres are specialized chromatin structures that protect chromosome ends from dangerous processing events. In most tissues, telomeres shorten with each round of cell division, placing a finite limit on cell growth. In rapidly dividing cells, including the majority of human cancers, cells bypass this growth limit through telomerase-catalyzed maintenance of telomere length. The dynamic properties of telomeres and telomerase render them difficult to study using ensemble biochemical and structural techniques. This review describes single-molecule approaches to studying how individual components of telomeres and telomerase contribute to function. Single-molecule methods provide a window into the complex nature of telomeres and telomerase by permitting researchers to directly visualize and manipulate the individual protein, DNA, and RNA molecules required for telomere function. The work reviewed in this article highlights how single-molecule techniques have been utilized to investigate the function of telomeres and telomerase.

Folding behavior of a T-shaped, ribosome-binding translation enhancer implicated in a wide-spread conformational switch

Turnip crinkle virus contains a T-shaped, ribosome-binding, translation enhancer (TSS) in its 3'UTR that serves as a hub for interactions throughout the region. The viral RNA-dependent RNA polymerase (RdRp) causes the TSS/surrounding region to undergo a conformational shift postulated to inhibit translation. Using optical tweezers (OT) and steered molecular dynamic simulations (SMD), we found that the unusual stability of pseudoknotted element H4a/Ψ3 required five upstream adenylates, and H4a/Ψ3 was necessary for cooperative association of two other hairpins (H5/H4b) in Mg2+. SMD recapitulated the TSS unfolding order in the absence of Mg2+, showed dependence of the resistance to pulling on the 3D orientation and gave structural insights into the measured contour lengths of the TSS structure elements. Adenylate mutations eliminated one-site RdRp binding to the 3'UTR, suggesting that RdRp binding to the adenylates disrupts H4a/Ψ3, leading to loss of H5/H4b interaction and promoting a conformational switch interrupting translation and promoting replication.

Mechanical unfolding kinetics of the SRV-1 gag-pro mRNA pseudoknot: possible implications for -1 ribosomal frameshifting stimulation

Minus-one ribosomal frameshifting is a translational recoding mechanism widely utilized by many RNA viruses to generate accurate ratios of structural and catalytic proteins. An RNA pseudoknot structure located in the overlapping region of the gag and pro genes of Simian Retrovirus type 1 (SRV-1) stimulates frameshifting. However, the experimental characterization of SRV-1 pseudoknot (un)folding dynamics and the effect of the base triple formation is lacking. Here, we report the results of our single-molecule nanomanipulation using optical tweezers and theoretical simulation by steered molecular dynamics. Our results directly reveal that the energetic coupling between loop 2 and stem 1 via minor-groove base triple formation enhances the mechanical stability. The terminal base pair in stem 1 (directly in contact with a translating ribosome at the slippery site) also affects the mechanical stability of the pseudoknot. The -1 frameshifting efficiency is positively correlated with the cooperative one-step unfolding force and inversely correlated with the one-step mechanical unfolding rate at zero force. A significantly improved correlation was observed between -1 frameshifting efficiency and unfolding rate at forces of 15-35 pN, consistent with the fact that the ribosome is a force-generating molecular motor with helicase activity. No correlation was observed between thermal stability and -1 frameshifting efficiency.

Salt Effects on the Thermodynamics of a Frameshifting RNA Pseudoknot under Tension

Because of the potential link between -1 programmed ribosomal frameshifting and response of a pseudoknot (PK) RNA to force, a number of single-molecule pulling experiments have been performed on PKs to decipher the mechanism of programmed ribosomal frameshifting. Motivated in part by these experiments, we performed simulations using a coarse-grained model of RNA to describe the response of a PK over a range of mechanical forces (fs) and monovalent salt concentrations (Cs). The coarse-grained simulations quantitatively reproduce the multistep thermal melting observed in experiments, thus validating our model. The free energy changes obtained in simulations are in excellent agreement with experiments. By varying f and C, we calculated the phase diagram that shows a sequence of structural transitions, populating distinct intermediate states. As f and C are changed, the stem-loop tertiary interactions rupture first, followed by unfolding of the 3'-end hairpin (I⇌F). Finally, the 5'-end hairpin unravels, producing an extended state (E⇌I). A theoretical analysis of the phase boundaries shows that the critical force for rupture scales as (logCm)(α) with α=1(0.5) for E⇌I (I⇌F) transition. This relation is used to obtain the preferential ion-RNA interaction coefficient, which can be quantitatively measured in single-molecule experiments, as done previously for DNA hairpins. A by-product of our work is the suggestion that the frameshift efficiency is likely determined by the stability of the 5'-end hairpin that the ribosome first encounters during translation.

Coarse-grained modeling of RNA 3D structure

Functional RNA molecules depend on three-dimensional (3D) structures to carry out their tasks within the cell. Understanding how these molecules interact to carry out their biological roles requires a detailed knowledge of RNA 3D structure and dynamics as well as thermodynamics, which strongly governs the folding of RNA and RNA-RNA interactions as well as a host of other interactions within the cellular environment. Experimental determination of these properties is difficult, and various computational methods have been developed to model the folding of RNA 3D structures and their interactions with other molecules. However, computational methods also have their limitations, especially when the biological effects demand computation of the dynamics beyond a few hundred nanoseconds. For the researcher confronted with such challenges, a more amenable approach is to resort to coarse-grained modeling to reduce the number of data points and computational demand to a more tractable size, while sacrificing as little critical information as possible. This review presents an introduction to the topic of coarse-grained modeling of RNA 3D structures and dynamics, covering both high- and low-resolution strategies. We discuss how physics-based approaches compare with knowledge based methods that rely on databases of information. In the course of this review, we discuss important aspects in the reasoning process behind building different models and the goals and pitfalls that can result.

Mimicking Ribosomal Unfolding of RNA Pseudoknot in a Protein Channel

Pseudoknots are a fundamental RNA tertiary structure with important roles in regulation of mRNA translation. Molecular force spectroscopic approaches such as optical tweezers can track the pseudoknot's unfolding intermediate states by pulling the RNA chain from both ends, but the kinetic unfolding pathway induced by this method may be different from that in vivo, which occurs during translation and proceeds from the 5' to 3' end. Here we developed a ribosome-mimicking, nanopore pulling assay for dissecting the vectorial unfolding mechanism of pseudoknots. The pseudoknot unfolding pathway in the nanopore, either from the 5' to 3' end or in the reverse direction, can be controlled by a DNA leader that is attached to the pseudoknot at the 5' or 3' ends. The different nanopore conductance between DNA and RNA translocation serves as a marker for the position and structure of the unfolding RNA in the pore. With this design, we provided evidence that the pseudoknot unfolding is a two-step, multistate, metal ion-regulated process depending on the pulling direction. Most notably, unfolding in both directions is rate-limited by the unzipping of the first helix domain (first step), which is Helix-1 in the 5' → 3' direction and Helix-2 in the 3' → 5' direction, suggesting that the initial unfolding step in either pulling direction needs to overcome an energy barrier contributed by the noncanonical triplex base-pairs and coaxial stacking interactions for the tertiary structure stabilization. These findings provide new insights into RNA vectorial unfolding mechanisms, which play an important role in biological functions including frameshifting.

Free energy landscape and multiple folding pathways of an H-type RNA pseudoknot

How RNA sequences fold to specific tertiary structures is one of the key problems for understanding their dynamics and functions. Here, we study the folding process of an H-type RNA pseudoknot by performing a large-scale all-atom MD simulation and bias-exchange metadynamics. The folding free energy landscapes are obtained and several folding intermediates are identified. It is suggested that the folding occurs via multiple mechanisms, including a step-wise mechanism starting either from the first helix or the second, and a cooperative mechanism with both helices forming simultaneously. Despite of the multiple mechanism nature, the ensemble folding kinetics estimated from a Markov state model is single-exponential. It is also found that the correlation between folding and binding of metal ions is significant, and the bound ions mediate long-range interactions in the intermediate structures. Non-native interactions are found to be dominant in the unfolded state and also present in some intermediates, possibly hinder the folding process of the RNA.

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).

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.

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 fluorescence resonance energy transfer studies of the human telomerase RNA pseudoknot: temperature-/urea-dependent folding kinetics and thermodynamics

The ribonucleoprotein telomerase is an RNA-dependent DNA polymerase that catalyzes the repetitive addition of a short, species-specific, DNA sequence to the ends of linear eukaryotic chromosomes. The single RNA component of telomerase contains both the template sequence for DNA synthesis and a functionally critical pseudoknot motif, which can also exist as a less stable hairpin. Here we use a minimal version of the human telomerase RNA pseudoknot to study this hairpin-pseudoknot structural equilibrium using temperature-controlled single-molecule fluorescence resonance energy transfer (smFRET) experiments. The urea dependence of these experiments aids in determination of the folding kinetics and thermodynamics. The wild-type pseudoknot behavior is compared and contrasted to a mutant pseudoknot sequence implicated in a genetic disorder-dyskeratosis congenita. These findings clearly identify that this 2nt noncomplementary mutation destabilizes the folding of the wild-type pseudoknot by substantially reducing the folding rate constant (≈ 400-fold) while only nominally increasing the unfolding rate constant (≈ 5-fold). Furthermore, the urea dependence of the equilibrium and rate constants is used to develop a free energy landscape for this unimolecular equilibrium and propose details about the structure of the transition state. Finally, the urea-dependent folding experiments provide valuable physical insights into the mechanism for destabilization of RNA pseudoknots by such chemical denaturants.

Unraveling secrets of telomeres: one molecule at a time

Telomeres play important roles in maintaining the stability of linear chromosomes. Telomere maintenance involves dynamic actions of multiple proteins interacting with long repetitive sequences and complex dynamic DNA structures, such as G-quadruplexes, T-loops and t-circles. Given the heterogeneity and complexity of telomeres, single-molecule approaches are essential to fully understand the structure-function relationships that govern telomere maintenance. In this review, we present a brief overview of the principles of single-molecule imaging and manipulation techniques. We then highlight results obtained from applying these single-molecule techniques for studying structure, dynamics and functions of G-quadruplexes, telomerase, and shelterin proteins.

Folding a stable RNA pseudoknot through rearrangement of two hairpin structures

Folding messenger RNA into specific structures is a common regulatory mechanism involved in translation. In Escherichia coli, the operator of the rpsO gene transcript folds into a pseudoknot or double-hairpin conformation. S15, the gene product, binds only to the pseudoknot, thereby repressing its own synthesis when it is present in excess in the cell. The two RNA conformations have been proposed to exist in equilibrium. However, it remained unclear how structural changes can be achieved between these two topologically distinct conformations. We used optical tweezers to study the structural dynamics and rearrangements of the rpsO operator RNA at the single-molecule level. We discovered that the two RNA structures can be interchanged spontaneously and the pseudoknot can exist in conformations that exhibit various levels of stability. Conversion from the double hairpin to a pseudoknot through potential hairpin–hairpin interactions favoured the high-stability conformation. By contrast, mutations that blocked the formation of a hairpin typically resulted in alternative low-stability pseudoknots. These results demonstrate that specific tertiary interactions of RNA can be established and modulated based on the interactions and rearrangements between secondary structural components. Our findings provide new insight into the RNA folding pathway that leads to a regulatory conformation for target protein binding.

Molecular dynamics simulation of RNA pseudoknot unfolding pathway

Many biological functions of RNA molecules are related to their pseudoknot structures. It is significant for predicting the structure and function of RNA that learning about the stability and the process of RNA pseudoknot folding and unfolding. The structural features of mouse mammary tumor virus (MMTV) RNA pseudoknot in different ion concentration, the unfolding process of the RNA pseudoknot, and the two hairpin helices that constitute the RNA pseudoknot were studied with all atom molecule dynamics simulation method in this paper. We found that the higher cation concentration can cause structure of the RNA molecules more stable, and ions played an indispensable role in keeping the structure of RNA molecules stable; the unfolding process of hairpin structure was corresponding to the antiprocess of its folding process. The main pathway of pseudoknot unfolding was that the inner base pair opened first, and then, the two helices, which formed the RNA pseudoknot opened decussately, while the folding pathway of the RNA pseudoknot was a helix folding after formation of the other helix. Therefore, the unfolding process of RNA pseudoknot is different from the antiprocess of its folding process, and the unfolding process of each helix in the RNA pseudoknot is similar to the hairpin structure’s unfolding process, which means that both are the unzipping process.

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.

Revealing -1 programmed ribosomal frameshifting mechanisms by single-molecule techniques and computational methods

Programmed ribosomal frameshifting (PRF) serves as an intrinsic translational regulation mechanism employed by some viruses to control the ratio between structural and enzymatic proteins. Most viral mRNAs which use PRF adapt an H-type pseudoknot to stimulate −1 PRF. The relationship between the thermodynamic stability and the frameshifting efficiency of pseudoknots has not been fully understood. Recently, single-molecule force spectroscopy has revealed that the frequency of −1 PRF correlates with the unwinding forces required for disrupting pseudoknots, and that some of the unwinding work dissipates irreversibly due to the torsional restraint of pseudoknots. Complementary to single-molecule techniques, computational modeling provides insights into global motions of the ribosome, whose structural transitions during frameshifting have not yet been elucidated in atomic detail. Taken together, recent advances in biophysical tools may help to develop antiviral therapies that target the ubiquitous −1 PRF mechanism among viruses.

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.

Fast folding of RNA pseudoknots initiated by laser temperature-jump

RNA pseudoknots are examples of minimal structural motifs in RNA with tertiary interactions that stabilize the structures of many ribozymes. They also play an essential role in a variety of biological functions that are modulated by their structure, stability, and dynamics. Therefore, understanding the global principles that determine the thermodynamics and folding pathways of RNA pseudoknots is an important problem in biology, both for elucidating the folding mechanisms of larger ribozymes as well as addressing issues of possible kinetic control of the biological functions of pseudoknots. We report on the folding/unfolding kinetics of a hairpin-type pseudoknot obtained with microsecond time-resolution in response to a laser temperature-jump perturbation. The kinetics are monitored using UV absorbance as well as fluorescence of extrinsically attached labels as spectroscopic probes of the transiently populated RNA conformations. We measure folding times of 1-6 ms at 37 °C, which are at least 100-fold faster than previous observations of very slow folding pseudoknots that were trapped in misfolded conformations. The measured relaxation times are remarkably similar to predictions of a computational study by Thirumalai and co-workers (Cho, S. S.; Pincus, D.L.; Thirumalai, D. Proc. Natl. Acad. Sci. U. S. A. 2009, 106, 17349-17354). Thus, these studies provide the first observation of a fast-folding pseudoknot and present a benchmark against which computational models can be refined.

Functional importance of telomerase pseudoknot revealed by single-molecule analysis

Telomerase ribonucleoprotein (RNP) employs an RNA subunit to template the addition of telomeric repeats onto chromosome ends. Previous studies have suggested that a region of the RNA downstream of the template may be important for telomerase activity and that the region could fold into a pseudoknot. Whether the pseudoknot motif is formed in the active telomerase RNP and what its functional role is have not yet been conclusively established. Using single-molecule FRET, we show that the isolated pseudoknot sequence stably folds into a pseudoknot. However, in the context of the full-length telomerase RNA, interference by other parts of the RNA prevents the formation of the pseudoknot. The protein subunits of the telomerase holoenzyme counteract RNA-induced misfolding and allow a significant fraction of the RNPs to form the pseudoknot structure. Only those RNP complexes containing a properly folded pseudoknot are catalytically active. These results not only demonstrate the functional importance of the pseudoknot but also reveal the critical role played by telomerase proteins in pseudoknot folding.

Atomistic analysis of pseudoknotted RNA unfolding

As important functional structures, RNA pseudoknots provide excellent models for studying the interplay between secondary and tertiary structures and the roles of triplexes, noncanonical interactions, and coaxial stacking in the folding/unfolding process. Here we report a first atomistic and statistical analysis of the unfolding of the pseudoknot within gene 32 mRNA of bacteriophage T2. Multiple unfolding pathways, diverse transition states, and various intermediate structures were observed. Water molecules were found to be coupled with the unfolding process via the expulsion or concurrent mechanism.

Cooperative and directional folding of the preQ1 riboswitch aptamer domain

Riboswitches are cis-acting RNA fragments that regulate gene expression by sensing cellular levels of the associated small metabolites. In bacteria, the class I preQ(1) riboswitch allows the fine-tuning of queuosine biosynthesis in response to the intracellular concentration of the queuosine anabolic intermediate preQ(1). When binding preQ(1), the aptamer domain undergoes a significant degree of secondary and tertiary structural rearrangement and folds into an H-type pseudoknot. Conformational "switching" of the riboswitch aptamer domain upon recognizing its cognate metabolite plays a key role in the regulatory mechanism of the preQ(1) riboswitch. We investigate the folding mechanism of the preQ(1) riboswitch aptamer domain using all-atom Go̅-model simulations. The folding pathway of such a single domain is found to be cooperative and sequentially coordinated, as the folding proceeds in the 5' → 3' direction. This kinetically efficient folding mechanism suggests a fast ligand-binding response in competition with RNA elongation.

Current trends in finite-time thermodynamics

The cornerstone of finite-time thermodynamics is all about the price of haste and how to minimize it. Reversible processes may be ultimately efficient, but they are unrealistically slow. In all situations-chemical, mechanical, economical-we pay extra to get the job done quickly. Finite-time thermodynamics can be used to develop methods to limit that extra expenditure, be it in energy, entropy production, money, or something entirely different. Finite-time thermodynamics also includes methods to calculate the optimal path or mode of operation to achieve this minimal expenditure. The concept is to place the system of interest in contact with a time-varying environment which will coax the system along the desired path, much like guiding a horse along by waving a carrot in front of it.