From Enigma to Revelation: Unravelling Biological Functions of Ubiquitous Small Ribozymes

RNA, widely recognized as an information-carrier molecule, is capable of catalyzing essential biological processes through ribozymes. Despite their ubiquity, specific functions in a biological context and phenotypes based on the ribozymes' activity are often unknown. Here, we present the discovery of a subgroup of minimal HDV-like ribozymes, which reside 3' to viral tRNAs and appear to cleave the 3'-trailers of viral premature tRNA transcripts. This proposed tRNA-processing function is unprecedented for any ribozymes, thus, we designate this subgroup as theta ribozymes. Most theta ribozymes were identified in Caudoviricetes bacteriophages, the main constituent (>90%) of the mammalian gut virome. Intriguingly, our findings further suggest the involvement of theta ribozymes in the transition of certain bacteriophages between distinct genetic codes, thus possibly contributing to the phage lysis trigger. Our discovery expands the limited repertoire of biological functions attributed to HDV-like ribozymes and provides insights into the fascinating world of RNA catalysis.

Identification of HDV-like theta ribozymes involved in tRNA-based recoding of gut bacteriophages

Trillions of microorganisms, collectively known as the microbiome, inhabit our bodies with the gut microbiome being of particular interest in biomedical research. Bacteriophages, the dominant virome constituents, can utilize suppressor tRNAs to switch to alternative genetic codes (e.g., the UAG stop-codon is reassigned to glutamine) while infecting hosts with the standard bacterial code. However, what triggers this switch and how the bacteriophage manipulates its host is poorly understood. Here, we report the discovery of a subgroup of minimal hepatitis delta virus (HDV)-like ribozymes - theta ribozymes - potentially involved in the code switch leading to the expression of recoded lysis and structural phage genes. We demonstrate their HDV-like self-scission behavior in vitro and find them in an unreported context often located with their cleavage site adjacent to tRNAs, indicating a role in viral tRNA maturation and/or regulation. Every fifth associated tRNA is a suppressor tRNA, further strengthening our hypothesis. The vast abundance of tRNA-associated theta ribozymes - we provide 1753 unique examples - highlights the importance of small ribozymes as an alternative to large enzymes that usually process tRNA 3'-ends. Our discovery expands the short list of biological functions of small HDV-like ribozymes and introduces a previously unknown player likely involved in the code switch of certain recoded gut bacteriophages.

Organometallic Pillarplexes That Bind DNA 4-Way Holliday Junctions and Forks

Holliday 4-way junctions are key to important biological DNA processes (insertion, recombination, and repair) and are dynamic structures that adopt either open or closed conformations, the open conformation being the biologically active form. Tetracationic metallo-supramolecular pillarplexes display aryl faces about a cylindrical core, an ideal structure to interact with open DNA junction cavities. Combining experimental studies and MD simulations, we show that an Au pillarplex can bind DNA 4-way (Holliday) junctions in their open form, a binding mode not accessed by synthetic agents before. Pillarplexes can bind 3-way junctions too, but their large size leads them to open up and expand that junction, disrupting the base pairing, which manifests in an increased hydrodynamic size and lower junction thermal stability. At high loading, they rearrange both 4-way and 3-way junctions into Y-shaped forks to increase the available junction-like binding sites. Isostructural Ag pillarplexes show similar DNA junction binding behavior but lower solution stability. This pillarplex binding contrasts with (but complements) that of metallo-supramolecular cylinders, which prefer 3-way junctions and can rearrange 4-way junctions into 3-way junction structures. The pillarplexes' ability to bind open 4-way junctions creates exciting possibilities to modulate and switch such structures in biology, as well as in synthetic nucleic acid nanostructures. In human cells, the pillarplexes do reach the nucleus, with antiproliferative activity at levels similar to those of cisplatin. The findings provide a new roadmap for targeting higher-order junction structures using a metallo-supramolecular approach, as well as expanding the toolbox available to design bioactive junction binders into organometallic chemistry.

The structural features of the ligand-free moaA riboswitch and its ion-dependent folding

Riboswitches are structural elements of mRNA involved in the regulation of gene expression by responding to specific cellular metabolites. To fulfil their regulatory function, riboswitches prefold into an active state, the so-called binding competent form, that guarantees metabolite binding and allows a consecutive refolding of the RNA. Here, we describe the folding pathway to the binding competent form as well as the ligand free structure of the moaA riboswitch of E. coli. This RNA proposedly responds to the molybdenum cofactor (Moco), a highly oxygen-sensitive metabolite, essential in the carbon and sulfur cycles of eukaryotes. K⁺- and Mg²⁺-dependent footprinting assays and spectroscopic investigations show a high degree of structure formation of this RNA already at very low ion-concentrations. Mg²⁺ facilitates additionally a general compaction of the riboswitch towards its proposed active structure. We show that this fold agrees with the earlier suggested secondary structure which included also a long-range tetraloop/tetraloop-receptor like interaction. Metal ion cleavage assays revealed specific Mg²⁺-binding pockets within the moaA riboswitch. These Mg²⁺ binding pockets are good indicators for the potential Moco binding site, since in riboswitches, Mg²⁺ was shown to be necessary to bind phosphate-carrying metabolites. The importance of the phosphate and of other functional groups of Moco is highlighted by binding assays with tetrahydrobiopterin, the reduced and oxygen-sensitive core moiety of Moco. We demonstrate that the general molecular shape of pterin by its own is insufficient for the recognition by the riboswitch.

Group II Introns: Highly Structured yet Dynamic

RNA splicing, the removal of introns and ligation of exons, is a crucial process during mRNA maturation. Group II introns are large ribozymes that self-catalyze their splicing, as well as their transposition. They are living fossils of spliceosomal introns and eukaryotic retroelements. The yeast mitochondrial Sc.ai5γ is the first identified and best-studied self-splicing group II intron. A combination of biochemical, biophysical, and computational tools enables studying its catalytic properties, structure, and dynamics, while also serving to develop new therapeutic and biotechnological tools. We survey the history of group II intron studies paralleling the trends in RNA methodology with Sc.ai5γ in the spotlight.

Deconvolution of 1D NMR spectra: A deep learning-based approach

The analysis of nuclear magnetic resonance (NMR) spectra to detect peaks and characterize their parameters, often referred to as deconvolution, is a crucial step in the quantification, elucidation, and verification of the structure of molecular systems. However, deconvolution of 1D NMR spectra is a challenge for both experts and machines. We propose a robust, expert-level quality deep learning-based deconvolution algorithm for 1D experimental NMR spectra. The algorithm is based on a neural network trained on synthetic spectra. Our customized pre-processing and labeling of the synthetic spectra enable the estimation of critical peak parameters. Furthermore, the neural network model transfers well to the experimental spectra and demonstrates low fitting errors and sparse peak lists in challenging scenarios such as crowded, high dynamic range, shoulder peak regions as well as broad peaks. We demonstrate in challenging spectra that the proposed algorithm is superior to expert results.

A blind benchmark of analysis tools to infer kinetic rate constants from single-molecule FRET trajectories

Single-molecule FRET (smFRET) is a versatile technique to study the dynamics and function of biomolecules since it makes nanoscale movements detectable as fluorescence signals. The powerful ability to infer quantitative kinetic information from smFRET data is, however, complicated by experimental limitations. Diverse analysis tools have been developed to overcome these hurdles but a systematic comparison is lacking. Here, we report the results of a blind benchmark study assessing eleven analysis tools used to infer kinetic rate constants from smFRET trajectories. We test them against simulated and experimental data containing the most prominent difficulties encountered in analyzing smFRET experiments: different noise levels, varied model complexity, non-equilibrium dynamics, and kinetic heterogeneity. Our results highlight the current strengths and limitations in inferring kinetic information from smFRET trajectories. In addition, we formulate concrete recommendations and identify key targets for future developments, aimed to advance our understanding of biomolecular dynamics through quantitative experiment-derived models.

A blind benchmark of analysis tools to infer kinetic rate constants from single-molecule FRET trajectories

Single-molecule FRET (smFRET) is a versatile technique to study the dynamics and function of biomolecules since it makes nanoscale movements detectable as fluorescence signals. The powerful ability to infer quantitative kinetic information from smFRET data is, however, complicated by experimental limitations. Diverse analysis tools have been developed to overcome these hurdles but a systematic comparison is lacking. Here, we report the results of a blind benchmark study assessing eleven analysis tools used to infer kinetic rate constants from smFRET trajectories. We test them against simulated and experimental data containing the most prominent difficulties encountered in analyzing smFRET experiments: different noise levels, varied model complexity, non-equilibrium dynamics, and kinetic heterogeneity. Our results highlight the current strengths and limitations in inferring kinetic information from smFRET trajectories. In addition, we formulate concrete recommendations and identify key targets for future developments, aimed to advance our understanding of biomolecular dynamics through quantitative experiment-derived models.

Magnesium(II)‐ATP Complexes in 1‐Ethyl‐3‐Methylimidazolium Acetate Solutions Characterized by 31 Mg β‐Radiation‐Detected NMR Spectroscopy

The complexation of Mg^II with adenosine 5′‐triphosphate (ATP) is omnipresent in biochemical energy conversion, but is difficult to interrogate directly. Here we use the spin‐1/2 β‐emitter ³¹Mg to study Mg^II‐ATP complexation in 1‐ethyl‐3‐methylimidazolium acetate (EMIM‐Ac) solutions using β‐radiation‐detected nuclear magnetic resonance (β‐NMR). We demonstrate that (nuclear) spin‐polarized ³¹Mg, following ion‐implantation from an accelerator beamline into EMIM‐Ac, binds to ATP within its radioactive lifetime before depolarizing. The evolution of the spectra with solute concentration indicates that the implanted ³¹Mg initially bind to the solvent acetate anions, whereafter they undergo dynamic exchange and form either a mono‐ (³¹Mg‐ATP) or di‐nuclear (³¹MgMg‐ATP) complex. The chemical shift of ³¹Mg‐ATP is observed up‐field of ³¹MgMg‐ATP, in accord with quantum chemical calculations. These observations constitute a crucial advance towards using β‐NMR to probe chemistry and biochemistry in solution.

Coordination Chemistry of Nucleotides and Antivirally Active Acyclic Nucleoside Phosphonates, including Mechanistic Considerations

Considering that practically all reactions that involve nucleotides also involve metal ions, it is evident that the coordination chemistry of nucleotides and their derivatives is an essential corner stone of biological inorganic chemistry. Nucleotides are either directly or indirectly involved in all processes occurring in Nature. It is therefore no surprise that the constituents of nucleotides have been chemically altered-that is, at the nucleobase residue, the sugar moiety, and also at the phosphate group, often with the aim of discovering medically useful compounds. Among such derivatives are acyclic nucleoside phosphonates (ANPs), where the sugar moiety has been replaced by an aliphatic chain (often also containing an ether oxygen atom) and the phosphate group has been replaced by a phosphonate carrying a carbon-phosphorus bond to make the compounds less hydrolysis-sensitive. Several of these ANPs show antiviral activity, and some of them are nowadays used as drugs. The antiviral activity results from the incorporation of the ANPs into the growing nucleic acid chain-i.e., polymerases accept the ANPs as substrates, leading to chain termination because of the missing 3'-hydroxyl group. We have tried in this review to describe the coordination chemistry (mainly) of the adenine nucleotides AMP and ATP and whenever possible to compare it with that of the dianion of 9-[2-(phosphonomethoxy)ethyl]adenine (PMEA2- = adenine(N9)-CH2-CH2-O-CH2-PO32) [or its diphosphate (PMEApp4-)] as a representative of the ANPs. Why is PMEApp4- a better substrate for polymerases than ATP4-? There are three reasons: (i) PMEA2- with its anti-like conformation (like AMP2-) fits well into the active site of the enzyme. (ii) The phosphonate group has an enhanced metal ion affinity because of its increased basicity. (iii) The ether oxygen forms a 5-membered chelate with the neighboring phosphonate and favors thus coordination at the Pα group. Research on ANPs containing a purine residue revealed that the kind and position of the substituent at C2 or C6 has a significant influence on the biological activity. For example, the shift of the (C6)NH2 group in PMEA to the C2 position leads to 9-[2-(phosphonomethoxy)ethyl]-2-aminopurine (PME2AP), an isomer with only a moderate antiviral activity. Removal of (C6)NH2 favors N7 coordination, e.g., of Cu2+, whereas the ether O atom binding of Cu2+ in PMEA facilitates N3 coordination via adjacent 5- and 7-membered chelates, giving rise to a Cu(PMEA)cl/O/N3 isomer. If the metal ions (M2+) are M(α,β)-M(γ)-coordinated at a triphosphate chain, transphosphorylation occurs (kinases, etc.), whereas metal ion binding in a M(α)-M(β,γ)-type fashion is relevant for polymerases. It may be noted that with diphosphorylated PMEA, (PMEApp4-), the M(α)-M(β,γ) binding is favored because of the formation of the 5-membered chelate involving the ether O atom (see above). The self-association tendency of purines leads to the formation of dimeric [M2(ATP)]2(OH)- stacks, which occur in low concentration and where one half of the molecule undergoes the dephosphorylation reaction and the other half stabilizes the structure-i.e., acts as the "enzyme" by bridging the two ATPs. In accord herewith, one may enhance the reaction rate by adding AMP2- to the [Cu2(ATP)]2(OH)- solution, as this leads to the formation of mixed stacked Cu3(ATP)(AMP)(OH)- species, in which AMP2- takes over the structuring role, while the other "half" of the molecule undergoes dephosphorylation. It may be added that Cu3(ATP)(PMEA) or better Cu3(ATP)(PMEA)(OH)- is even a more reactive species than Cu3(ATP)(AMP)(OH)-. - The matrix-assisted self-association and its significance for cell organelles with high ATP concentrations is summarized and discussed, as is, e.g., the effect of tryptophanate (Trp-), which leads to the formation of intramolecular stacks in M(ATP)(Trp)3- complexes (formation degree about 75%). Furthermore, it is well-known that in the active-site cavities of enzymes the dielectric constant, compared with bulk water, is reduced; therefore, we have summarized and discussed the effect of a change in solvent polarity on the stability and structure of binary and ternary complexes: Opposite effects on charged O sites and neutral N sites are observed, and this leads to interesting insights.

Single-Molecule Kinetic Studies of Nucleic Acids by Förster Resonance Energy Transfer

Single-molecule microscopy is often used to observe and characterize the conformational dynamics of nucleic acids (NA). Due to the large variety of NA structures and the challenges specific to single-molecule observation techniques, the data recorded in such experiments must be processed via multiple statistical treatments to finally yield a reliable mechanistic view of the NA dynamics. In this chapter, we propose a comprehensive protocol to analyze single-molecule trajectories in the scope of single-molecule Förster resonance energy transfer (FRET) microscopy. The suggested protocol yields the conformational states common to all molecules in the investigated sample, together with the associated conformational transition kinetics. The given model resolves states that are indistinguishable by their observed FRET signals and is estimated with 95% confidence using error calculations on FRET states and transition rate constants. In the end, a step-by-step user guide is given to reproduce the protocol with the Multifunctional Analysis Software to Handle single-molecule FRET data (MASH-FRET).

Chemical Dual End-Labeling of Large Ribozymes

Fast and efficient site-specific labeling of long RNAs is one of the main bottlenecks limiting distance measurements by means of Förster resonance energy transfer (FRET) or electron paramagnetic resonance (EPR) spectroscopy. Here, we present an optimized protocol for dual end-labeling with different fluorophores at the same time meeting the restrictions of highly labile and degradation-sensitive RNAs. We describe in detail the dual-labeling of a catalytically active wild-type group II intron as a typical representative of long functional RNAs. The modular procedure chemically activates the 5'-phosphate and the 3'-ribose for bioconjugation with a pair of fluorophores, as shown herein, or with spin labels. The mild reaction conditions preserve the structural and functional integrity of the biomacromolecule and results in covalent, dual-labeled RNA in its pre-catalytic state in yields suitable for both ensemble and single-molecule FRET experiments.

FRETraj: Integrating single-molecule spectroscopy with molecular dynamics

SUMMARY

Quantitative interpretation of single-molecule FRET experiments requires a model of the dye dynamics to link experimental energy transfer efficiencies to distances between atom positions. We have developed FRETraj, a Python module to predict FRET distributions based on accessible-contact volumes (ACV) and simulated photon statistics. FRETraj helps to identify optimal fluorophore positions on a biomolecule of interest by rapidly evaluating donor-acceptor distances. FRETraj is scalable and fully integrated into PyMOL and the Jupyter ecosystem. Here we describe the conformational dynamics of a DNA hairpin by computing multiple ACVs along a molecular dynamics trajectory and compare the predicted FRET distribution with single-molecule experiments. FRET-assisted modeling will accelerate the analysis of structural ensembles in particular dynamic, non-coding RNAs and transient protein-nucleic acid complexes.

AVAILABILITY

FRETraj is implemented as a cross-platform Python package available under the GPL-3.0 on Github (https://github.com/RNA-FRETools/fretraj) and is documented at https://RNA-FRETools.github.io/fretraj.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Site-Specific Dual-Color Labeling of Long RNAs

Labeling of large RNAs with reporting entities, e.g., fluorophores, has significant impact on RNA studies in vitro and in vivo. Here, we describe a minimally invasive RNA labeling method featuring nucleotide and position selectivity, which solves the long-standing challenge of how to achieve accurate site-specific labeling of large RNAs with a least possible influence on folding and/or function. We use a custom-designed reactive DNA strand to hybridize to the RNA and transfer the alkyne group onto the targeted adenine or cytosine. Simultaneously, the 3'-terminus of RNA is converted to a dialdehyde moiety under the experimental condition applied. The incorporated functionalities at the internal and the 3'-terminal sites can then be conjugated with reporting entities via bioorthogonal chemistry. This method is particularly valuable for, but not limited to, single-molecule fluorescence applications. We demonstrate the method on an RNA construct of 275 nucleotides, the btuB riboswitch of Escherichia coli.

Screening for potential interaction partners with surface plasmon resonance imaging coupled to MALDI mass spectrometry

We coupled SPR imaging (SPRi) with matrix-assisted laser desorption/ionization mass spectrometry (MALDI MS) to identify new potential RNA binders. Here, we improve this powerful method, especially by optimizing the proteolytic digestion (type of reducing agent, its concentration, and incubation time), to work with complex mixtures, specifically a lysate of the rough mitochondrial fraction from yeast. The advantages of this hyphenated method compared to column-based or separate analyses are (i) rapid and direct visual readout from the SPRi array, (ii) possibility of high-throughput analysis of different interactions in parallel, (iii) high sensitivity, and (iv) no sample loss or contamination due to elution or micro-recovery procedures. The model system used is a catalytically active RNA (group IIB intron from Saccharomyces cerevisiae, Sc.ai5γ) and its cofactor Mss116. The protein supports the RNA folding process and thereby the subsequent excision of the intronic RNA from the coding part. Using the novel approach of coupling SPR with MALDI MS, we report the identification of potential RNA-binding proteins from a crude yeast mitochondrial lysate in a non-targeted approach. Our results show that proteins other than the well-known cofactor Mss116 interact with Sc.ai5γ (Dbp8, Prp8, Mrp13, and Cullin-3), suggesting that the intron folding and splicing are regulated by more than one cofactor in vivo.

Metal ions and sugar puckering balance single-molecule kinetic heterogeneity in RNA and DNA tertiary contacts

The fidelity of group II intron self-splicing and retrohoming relies on long-range tertiary interactions between the intron and its flanking exons. By single-molecule FRET, we explore the binding kinetics of the most important, structurally conserved contact, the exon and intron binding site 1 (EBS1/IBS1). A comparison of RNA-RNA and RNA-DNA hybrid contacts identifies transient metal ion binding as a major source of kinetic heterogeneity which typically appears in the form of degenerate FRET states. Molecular dynamics simulations suggest a structural link between heterogeneity and the sugar conformation at the exon-intron binding interface. While Mg2+ ions lock the exon in place and give rise to long dwell times in the exon bound FRET state, sugar puckering alleviates this structural rigidity and likely promotes exon release. The interplay of sugar puckering and metal ion coordination may be an important mechanism to balance binding affinities of RNA and DNA interactions in general.

Encapsulation of Fluorescently Labeled RNAs into Surface-Tethered Vesicles for Single-Molecule FRET Studies in TIRF Microscopy

Imaging fluorescently labeled biomolecules on a single-molecule level is a well-established technique to follow intra- and intermolecular processes in time, usually hidden in the ensemble average. The classical approach comprises surface immobilization of the molecule of interest, which increases the risk of restricting the natural behavior due to surface interactions. Encapsulation of such biomolecules into surface-tethered phospholipid vesicles enables to follow one molecule at a time, freely diffusing and without disturbing surface interactions. Further, the encapsulation allows to keep reaction partners (reactants and products) in close proximity and enables higher temperatures otherwise leading to desorption of the direct immobilized biomolecules.Here, we describe a detailed protocol for the encapsulation of a catalytically active RNA starting from surface passivation over RNA encapsulation to data evaluation of single-molecule FRET experiments in TIRF microscopy. We present an optimized procedure that preserves RNA functionality and applies to investigations of, e.g., large ribozymes and RNAs, where direct immobilization is structurally not possible.

NMR solution structure of tricyclo-DNA containing duplexes: insight into enhanced thermal stability and nuclease resistance

Tc-DNA is a conformationally constrained oligonucleotide analogue which shows significant increase in thermal stability when hybridized with RNA, DNA or tc-DNA. Remarkably, recent studies revealed that tc-DNA antisense oligonucleotides (AO) hold great promise for the treatment of Duchenne muscular dystrophy and spinal muscular atrophy. To date, no high-resolution structural data is available for fully modified tc-DNA duplexes and little is known about the origins of their enhanced thermal stability. Here, we report the structures of a fully modified tc-DNA oligonucleotide paired with either complementary RNA, DNA or tc-DNA. All three investigated duplexes maintain a right-handed helical structure with Watson-Crick base pairing and overall geometry intermediate between A- and B-type, but closer to A-type structures. All sugars of the tc-DNA and RNA residues adopt a North conformation whereas the DNA deoxyribose are found in a South-East-North conformation equilibrium. The conformation of the tc-DNA strand in the three determined structures is nearly identical and despite the different nature and local geometry of the complementary strand, the overall structures of the examined duplexes are very similar suggesting that the tc-DNA strand dominates the duplex structure.

Concerted dynamics of metallo-base pairs in an A/B-form helical transition

Metal-mediated base pairs expand the repertoire of nucleic acid structures and dynamics. Here we report solution structures and dynamics of duplex DNA containing two all-natural C-Hg^II-T metallo base pairs separated by six canonical base pairs. NMR experiments reveal a 3:1 ratio of well-resolved structures in dynamic equilibrium. The major species contains two (N3)T-Hg^II-(N3)C base pairs in a predominantly B-form helix. The minor species contains (N3)T-Hg^II-(N4)C base pairs and greater A-form characteristics. Ten-fold different ¹J coupling constants (¹⁵N,¹⁹⁹Hg) are observed for (N3)C-Hg^II (114 Hz) versus (N4)C-Hg^II (1052 Hz) connectivities, reflecting differences in cytosine ionization and metal-bonding strengths. Dynamic interconversion between the two types of C-Hg^II-T base pairs are coupled to a global conformational exchange between the helices. These observations inspired the design of a repetitive DNA sequence capable of undergoing a global B-to-A-form helical transition upon adding Hg^II, demonstrating that C-Hg^II-T has unique switching potential in DNA-based materials and devices.

Site-specific dual-color labeling of long RNAs for single-molecule spectroscopy

Labeling of long RNA molecules in a site-specific yet generally applicable manner is integral to many spectroscopic applications. Here we present a novel covalent labeling approach that is site-specific and scalable to long intricately folded RNAs. In this approach, a custom-designed DNA strand that hybridizes to the RNA guides a reactive group to target a preselected adenine residue. The functionalized nucleotide along with the concomitantly oxidized 3'-terminus can subsequently be conjugated to two different fluorophores via bio-orthogonal chemistry. We validate this modular labeling platform using a regulatory RNA of 275 nucleotides, the btuB riboswitch of Escherichia coli, demonstrate its general applicability by modifying a base within a duplex, and show its site-selectivity in targeting a pair of adjacent adenines. Native folding and function of the RNA is confirmed on the single-molecule level by using FRET as a sensor to visualize and characterize the conformational equilibrium of the riboswitch upon binding of its cofactor adenosylcobalamin. The presented labeling strategy overcomes size and site constraints that have hampered routine production of labeled RNA that are beyond 200 nt in length.

Optimal molecular crowding accelerates group II intron folding and maximizes catalysis

Unlike in vivo conditions, group II intron ribozymes are known to require high magnesium(II) concentrations ([Mg2+]) and high temperatures (42 °C) for folding and catalysis in vitro. A possible explanation for this difference is the highly crowded cellular environment, which can be mimicked in vitro by macromolecular crowding agents. Here, we combined bulk activity assays and single-molecule Förster Resonance Energy Transfer (smFRET) to study the influence of polyethylene glycol (PEG) on catalysis and folding of the ribozyme. Our activity studies reveal that PEG reduces the [Mg2+] required, and we found an "optimum" [PEG] that yields maximum activity. smFRET experiments show that the most compact state population, the putative active state, increases with increasing [PEG]. Dynamic transitions between folded states also increase. Therefore, this study shows that optimal molecular crowding concentrations help the ribozyme not only to reach the native fold but also to increase its in vitro activity to approach that in physiological conditions.

Influence of pH and Mg(ii) on the catalytic core domain 5 of a bacterial group II intron

RNA molecules fold into complex structures that allow them to perform specific functions. To compensate the relative lack of diversity of functional groups within nucleotides, metal ions work as crucial co-factors. In addition, shifted pK_as are observed in RNA, enabling acid-base reactions at ambient pH. The central catalytic domain 5 (D5) hairpin of the Azotobacter vinelandii group II intron undergoes both metal ion binding and pH dependence, presumably playing an important functional role in the ribozyme's reaction. By NMR spectroscopy we have here characterized the metal ion binding sites and affinities for the hairpin's internal G-A mismatch, bulge, and pentaloop. The influence of Mg(ii) and pH on the local conformation of the catalytically crucial region is also explored by fluorescence spectroscopy.

Reliable State Identification and State Transition Detection in Fluorescence Intensity-Based Single-Molecule Förster Resonance Energy-Transfer Data

Single-molecule Förster resonance energy transfer (smFRET) is a powerful technique to probe biomolecular structure and dynamics. A popular implementation of smFRET consists of recording fluorescence intensity time traces of surface-immobilized, chromophore-tagged molecules. This approach generates large and complex data sets, the analysis of which is to date not standardized. Here, we address a key challenge in smFRET data analysis: the generation of thermodynamic and kinetic models that describe with statistical rigor the behavior of FRET trajectories recorded from surface-tethered biomolecules in terms of the number of FRET states, the corresponding mean FRET values, and the kinetic rates at which they interconvert. For this purpose, we first perform Monte Carlo simulations to generate smFRET trajectories, in which a relevant space of experimental parameters is explored. Then, we provide an account on current strategies to achieve such model selection, as well as a quantitative assessment of their performances. Specifically, we evaluate the performance of each algorithm (change-point analysis, STaSI, HaMMy, vbFRET, and ebFRET) with respect to accuracy, reproducibility, and computing time, which yields a range of algorithm-specific referential benchmarks for various data qualities. Data simulation and analysis were performed with our MATLAB-based multifunctional analysis software for handling smFRET data (MASH-FRET).

Simulations of camera-based single-molecule fluorescence experiments

Single-molecule microscopy has become a widely used technique in (bio)physics and (bio)chemistry. A popular implementation is single-molecule Förster Resonance Energy Transfer (smFRET), for which total internal reflection fluorescence microscopy is frequently combined with camera-based detection of surface-immobilized molecules. Camera-based smFRET experiments generate large and complex datasets and several methods for video processing and analysis have been reported. As these algorithms often address similar aspects in video analysis, there is a growing need for standardized comparison. Here, we present a Matlab-based software (MASH-FRET) that allows for the simulation of camera-based smFRET videos, yielding standardized data sets suitable for benchmarking video processing algorithms. The software permits to vary parameters that are relevant in cameras-based smFRET, such as video quality, and the properties of the system under study. Experimental noise is modeled taking into account photon statistics and camera noise. Finally, we survey how video test sets should be designed to evaluate currently available data analysis strategies in camera-based sm fluorescence experiments. We complement our study by pre-optimizing and evaluating spot detection algorithms using our simulated video test sets.

An atomistic view on carbocyanine photophysics in the realm of RNA

Carbocyanine dyes have a long-standing tradition in fluorescence imaging and spectroscopy, due to their photostability and large spectral separation between individual dye species. Herein, we explore the versatility of cyanine dyes to probe the dynamics of nucleic acids and we report on the interrelation of fluorophores, RNA, and metal ions, namely K⁺ and Mg²⁺. Photophysical parameters including the fluorescence lifetime, quantum yield and dynamic anisotropy are monitored as a function of the nucleic acid composition, conformation, and metal ion abundance. Occasional excursions to a non-fluorescent cis-state hint at the remarkable sensitivity of carbocyanines to their local environment. Comparison of time-correlated single photon experiments with all-atom molecular dynamics simulations demonstrate that the propensity of photoisomerization is dictated by sterical constraints imposed on the fluorophore. Structural features in the vicinity of the dye play a crucial role in RNA recognition and have far-reaching implications on the mobility of the fluorescent probe. An atomic level description of the mutual interactions will ultimately benefit the quantitative interpretation of single-molecule FRET measurements on large RNA systems.

Metal ion complexes of nucleoside phosphorothioates reflecting the ambivalent properties of lead(ii)

The lead(ii)-lone pair leads to ambivalency: hemidirected (distorted, non-spherical) coordination spheres result from electronegative O-coordination and holodirected (symmetric, spherical) ones from less electronegative S-coordination.

The Role of Lead(II) in Nucleic Acids

Although lead(II) is naturally not associated with nucleic acids, this metal ions has been applied with DNA and RNA in various contexts. Pb2+ is an excellent hydrolytic metal ion for nucleic acids, which is why it is mainly used as probing agent for secondary structure and to determine metal ion binding sites both in vitro and in vivo. A further application of lead(II) is in structural studies, i.e., NMR, but also in X-ray crystallography, mostly using this heavy metal to solve the phase problem in the latter method. The structures of tRNAPhe, RNase P, HIV-1 DIS, and the leadzyme are discussed here in detail. A major part of this review is devoted to the cleavage properties of lead(II) with RNA because of its excellence in catalyzing phosphodiester cleavage. Metal ion binding sites in large naturally occurring ribozymes are regularly determined by Pb2+ cleavage, and also in the in vitro selected socalled leadzyme, this metal ion is the decisive key to backbone cleavage at a specific site. Lead(II) was used in the first in vitro selection that yielded a catalytic DNA, i.e., the DNAzyme named GR5. Next to the GR5, the so-called 8-17E is the second most prominent DNAzyme today. Derivatives of these two lead(II)-dependent DNAzymes, as well as the G-quadruplex forming PS2.M have been applied to detect lead(II) in the lower nanomolar range not only in the test tube but also in body fluids. Due to the toxicity of lead(II) for living beings, this is a highly active research field. Finally, further applications of lead(II)-dependent DNAzymes, e.g., in the construction of nanocomputers, are also discussed.

Tb3+-Cleavage Assays Reveal Specific Mg2+ Binding Sites Necessary to Pre-fold the btuB Riboswitch for AdoCbl Binding

Riboswitches are RNA elements that bind specific metabolites in order to regulate the gene expression involved in controlling the cellular concentration of the respective molecule or ion. Ligand recognition is mostly facilitated by Mg²⁺ mediated pre-organization of the riboswitch to an active tertiary fold. To predict these specific Mg²⁺ induced tertiary interactions of the btuB riboswitch from E. coli, we here report Mg²⁺ binding pockets in its aptameric part in both, the ligand-free and the ligand-bound form. An ensemble of weak and strong metal ion binding sites distributed over the entire aptamer was detected by terbium(III) cleavage assays, Tb³⁺ being an established Mg²⁺ mimic. Interestingly many of the Mⁿ⁺ (n = 2 or 3) binding sites involve conserved bases within the class of coenzyme B₁₂-binding riboswitches. Comparison with the published crystal structure of the coenzyme B₁₂ riboswitch of S. thermophilum aided in identifying a common set of Mⁿ⁺ binding sites that might be crucial for tertiary interactions involved in the organization of the aptamer. Our results suggest that Mⁿ⁺ binding at strategic locations of the btuB riboswitch indeed facilitates the assembly of the binding pocket needed for ligand recognition. Binding of the specific ligand, coenzyme B₁₂ (AdoCbl), to the btuB aptamer does however not lead to drastic alterations of these Mⁿ⁺ binding cores, indicating the lack of a major rearrangement within the three-dimensional structure of the RNA. This finding is strengthened by Tb³⁺ mediated footprints of the riboswitch's structure in its ligand-free and ligand-bound state indicating that AdoCbl indeed induces local changes rather than a global structural rearrangement.

G-quadruplex DNA targeted metal complexes acting as potential anticancer drugs

This review summarizes the recent development of G4 DNA targeted metal complexes and discusses their potential as anticancer drugs.

Sequence-Specific Post-Synthetic Oligonucleotide Labeling for Single-Molecule Fluorescence Applications

The sequence-specific fluorescence labeling of nucleic acids is a prerequisite for various methods including single-molecule Förster resonance energy transfer (smFRET) for the detailed study of nucleic acid folding and function. Such nucleic acid derivatives are commonly obtained by solid-phase methods; however, yields decrease rapidly with increasing length and restrict the practicability of this approach for long strands. Here, we report a new labeling strategy for the postsynthetic incorporation of a bioorthogonal group into single stranded regions of both DNA and RNA of unrestricted length. A 12-alkyne-etheno-adenine modification is sequence-selectively formed using DNA-templated synthesis, followed by conjugation of the fluorophore Cy3 via a copper-catalyzed azide-alkyne cycloaddition (CuAAC). Evaluation of the labeled strands in smFRET measurements shows that the strategy developed here has the potential to be used for the study of long functional nucleic acids by (single-molecule) fluorescence or other methods. To prove the universal use of the method, its application was successfully extended to the labeling of a short RNA single strand. As a proof-of-concept, also the labeling of a large RNA molecule in form of a 633 nucleotide long construct derived from the Saccharomyces cerevisiae group II intron Sc.ai5γ was performed, and covalent attachment of the Cy3 fluorophore was shown with gel electrophoresis.

Secondary structure confirmation and localization of Mg2+ ions in the mammalian CPEB3 ribozyme

Most of today's knowledge of the CPEB3 ribozyme, one of the few small self-cleaving ribozymes known to occur in humans, is based on comparative studies with the hepatitis delta virus (HDV) ribozyme, which is highly similar in cleavage mechanism and probably also in structure. Here we present detailed NMR studies of the CPEB3 ribozyme in order to verify the formation of the predicted nested double pseudoknot in solution. In particular, the influence of Mg(2+), the ribozyme's crucial cofactor, on the CPEB3 structure is investigated. NMR titrations, Tb(3+)-induced cleavage, as well as stoichiometry determination by hydroxyquinoline sulfonic acid fluorescence and equilibrium dialysis, are used to evaluate the number, location, and binding mode of Mg(2+)ions. Up to eight Mg(2+)ions interact site-specifically with the ribozyme, four of which are bound with high affinity. The global fold of the CPEB3 ribozyme, encompassing 80%-90% of the predicted base pairs, is formed in the presence of monovalent ions alone. Low millimolar concentrations of Mg(2+)promote a more compact fold and lead to the formation of additional structures in the core of the ribozyme, which contains the inner small pseudoknot and the active site. Several Mg(2+)binding sites, which are important for the functional fold, appear to be located in corresponding locations in the HDV and CPEB3 ribozyme, demonstrating the particular relevance of Mg(2+)for the nested double pseudoknot structure.

Studying metal ion binding properties of a three-way junction RNA by heteronuclear NMR

Self-splicing group II introns are highly structured RNA molecules, containing a characteristic secondary and catalytically active tertiary structure, which is formed only in the presence of Mg(II). Mg(II) initiates the first folding step governed by the κζ element within domain 1 (D1κζ). We recently solved the NMR structure of D1κζ derived from the mitochondrial group II intron ribozyme Sc.ai5γ and demonstrated that Mg(II) is essential for its stabilization. Here, we performed a detailed multinuclear NMR study of metal ion interactions with D1κζ, using Cd(II) and cobalt(III)hexammine to probe inner- and outer-sphere coordination of Mg(II) and thus to better characterize its binding sites. Accordingly, we mapped (1)H, (15)N, (13)C, and (31)P spectral changes upon addition of different amounts of the metal ions. Our NMR data reveal a Cd(II)-assisted macrochelate formation at the 5'-end triphosphate, a preferential Cd(II) binding to guanines in a helical context, an electrostatic interaction in the ζ tetraloop receptor and various metal ion interactions in the GAAA tetraloop and κ element. These results together with our recently published data on Mg(II) interaction provide a much better understanding of Mg(II) binding to D1κζ, and reveal how intricate and complex metal ion interactions can be.

Mimicking the in vivo Environment--The Effect of Crowding on RNA and Biomacromolecular Folding and Activity

In vitro studies on macromolecules, like proteins and nucleic acids, are mostly carried out in highly diluted systems where the molecules are studied under artificial conditions. These experimental conditions are optimized for both the system under investigation and the technique used. However, these conditions often do not reflect the in vivo situation and are therefore inappropriate for a reliable prediction of the native behavior of the molecules and their interactions under in vivo conditions. The intracellular environment is packed with cosolutes (macromolecules, metabolites, etc.) that create 'macromolecular crowding'. The addition of natural or synthetic macromolecules to the sample solution enables crowding to be mimicked. In this surrounding most of the studied biomolecules show a more compact structure, an increased activity, and a decrease of salt requirement for structure formation and function. Herein, we refer to a collection of examples for proteins and nucleic acids and their interactions in crowding environments and present in detail the effect of cosolutes on RNA folding and activity using a group II intron ribozyme as an example.

Protonation-Dependent Base Flipping at Neutral pH in the Catalytic Triad of a Self-Splicing Bacterial Group II Intron

NMR spectroscopy has revealed pH-dependent structural changes in the highly conserved catalytic domain 5 of a bacterial group II intron. Two adenines with pK(a) values close to neutral pH were identified in the catalytic triad and the bulge. Protonation of the adenine opposite to the catalytic triad is stabilized within a G(syn)-AH(+) (anti) base pair. The pH-dependent anti-to-syn flipping of this G in the catalytic triad modulates the known interaction with the linker region between domains 2 and 3 (J23) and simultaneously the binding of the catalytic Mg(2+) ion to its backbone. Hence, this here identified shifted pK(a) value controls the conformational change between the two steps of splicing.

Strategy for Internal Labeling of Large RNAs with Minimal Perturbation by Using Fluorescent PNA

Fluorescence techniques for the investigation of biomolecules and their folding pathways require an efficient labeling strategy. A common method to internally label large RNAs involves the introduction of long loops for hybridization of fluorophore-carrying DNA strands. Such loops often disturb the structure, and thus the functionality, of the RNA. Here we show, in a proof of concept study with a >600 nucleotide group II intron ribozyme, that the usage of the nucleic acid analogue peptide nucleic acid (PNA) is more efficient in several aspects, minimizing the required structural modifications of the RNA. We demonstrate by various methods, including smFRET, that much smaller concentrations and shorter PNAs can be applied, compared to DNA, for rapid and specific internal RNA labeling. The folding pathway and catalytic activity of this large ribozyme is only minimally affected by the PNA, but the background signal is significantly reduced.

Explicit analytic equations for multimolecular thermal melting curves

The analysis of thermal melting curves requires the knowledge of equations for the temperature dependence of the relative fraction of folded and unfolded components. To implement these equations as standard tools for curve fitting, they should be as explicit as possible. From the van't Hoff formalism it is known that the equilibrium constant and hence the folded fraction is a function of the absolute temperature, the van't Hoff transition enthalpy, and the melting temperature. The work presented here is devoted to the mathematically self-contained derivation and the listing of explicit equations for the folded fraction as a function of the thermodynamic parameters in the case of arbitrary molecularities. Part of the results are known, others are new. It is in particular shown for the first time that the folded fraction is the composition of a universal function which depends solely on the molecularity and a dimensionless function which is governed by the concrete thermodynamic regime but is independent of the molecularity. The results will prove useful for extracting the thermodynamic parameters from experimental data on the basis of regression analysis. As supporting information, open-source Matlab scripts for the computer implementation of the equations are provided.

The role of Mg(II) in DNA cleavage site recognition in group II intron ribozymes: solution structure and metal ion binding sites of the RNA-DNA complex

Group II intron ribozymes catalyze the cleavage of (and their reinsertion into) DNA and RNA targets using a Mg2(+)-dependent reaction. The target is cleaved 3' to the last nucleotide of intron binding site 1 (IBS1), one of three regions that form base pairs with the intron's exon binding sites (EBS1 to -3).We solved the NMR solution structure of the d3' hairpin of the Sc.ai5γ intron containing EBS1 in its 11-nucleotide loop in complex with the dIBS1 DNA 7-mer and compare it with the analogous RNA-RNA contact. The EBS1-dIBS1 helix is slightly flexible and non-symmetric. NMR data reveal two major groove binding sites for divalent metal ions at the EBS1-dIBS1 helix, and surface plasmon resonance experiments show that low concentrations of Mg2(+) considerably enhance the affinity of dIBS1 for EBS1. Our results indicate that identification of both RNA and DNA IBS1 targets, presentation of the scissile bond, and stabilization of the structure by metal ions are governed by the overall structure of EBS1-dIBS1 and the surrounding loop nucleotides but are irrespective of different EBS1-(d)IBS1 geometries and interstrand affinities.

Structure determination of noncanonical RNA motifs guided by ¹H NMR chemical shifts

Structured non-coding RNAs underline fundamental cellular processes, but determining their 3D structures remains challenging. We demonstrate herein that integrating NMR ¹H chemical shift data with Rosetta de novo modeling can consistently return high-resolution RNA structures. On a benchmark set of 23 noncanonical RNA motifs, including 11 blind targets, Chemical-Shift-ROSETTA for RNA (CS-ROSETTA-RNA) recovered the experimental structures with high accuracy (0.6 to 2.0 Å all-heavy-atom rmsd) in 18 cases.

Monitoring global structural changes and specific metal-ion-binding sites in RNA by in-line probing and Tb(III) cleavage

In this chapter we describe the use of two methods, in-line probing as well as terbium(III) cleavage. Both methods can be applied to RNAs of any size, structure, and function. Aside from revealing directly metal ion-binding sites these techniques also provide structural information for longer RNA sequences that are out of range to be analyzed with other techniques such as NMR. The cleavage pattern derived from in-line probing experiments reflects local and overall conformational changes in RNA upon the addition of metal ions, metal complexes, or other ligands. On the other side, terbium(III) cleavage experiments are applied to locate specific metal ion-binding sites in RNA molecules.

BOBA FRET: bootstrap-based analysis of single-molecule FRET data

Time-binned single-molecule Förster resonance energy transfer (smFRET) experiments with surface-tethered nucleic acids or proteins permit to follow folding and catalysis of single molecules in real-time. Due to the intrinsically low signal-to-noise ratio (SNR) in smFRET time traces, research over the past years has focused on the development of new methods to extract discrete states (conformations) from noisy data. However, limited observation time typically leads to pronounced cross-sample variability, i.e., single molecules display differences in the relative population of states and the corresponding conversion rates. Quantification of cross-sample variability is necessary to perform statistical testing in order to assess whether changes observed in response to an experimental parameter (metal ion concentration, the presence of a ligand, etc.) are significant. However, such hypothesis testing has been disregarded to date, precluding robust biological interpretation. Here, we address this problem by a bootstrap-based approach to estimate the experimental variability. Simulated time traces are presented to assess the robustness of the algorithm in conjunction with approaches commonly used in thermodynamic and kinetic analysis of time-binned smFRET data. Furthermore, a pair of functionally important sequences derived from the self-cleaving group II intron Sc.ai5γ (d3'EBS1*/IBS1*) is used as a model system. Through statistical hypothesis testing, divalent metal ions are shown to have a statistically significant effect on both thermodynamic and kinetic aspects of their interaction. The Matlab source code used for analysis (bootstrap-based analysis of smFRET data, BOBA FRET), as well as a graphical user interface, is available via http://www.aci.uzh.ch/rna/.