Relating Structure and Dynamics in RNA Biology

Synthesis and Enzymatic Incorporation of a Dual-App Nucleotide Probe That Reports Antibiotics-Induced Conformational Change in the Bacterial Ribosomal Decoding Site RNA

Natural nucleosides are nonfluorescent and do not have intrinsic labels that can be readily utilized for analyzing nucleic acid structure and recognition. In this regard, researchers typically use the so-called "one-label, one-technique" approach to study nucleic acids. However, we envisioned that a responsive dual-app nucleoside system that harnesses the power of two complementing biophysical techniques namely, fluorescence and 19F NMR, will allow the investigation of nucleic acid conformations more comprehensively than before. We recently introduced a nucleoside analogue by tagging trifluoromethyl-benzofuran at the C5 position of 2'-deoxyuridine, which serves as an excellent fluorescent and 19F NMR probe to study G-quadruplex and i-motif structures. Taking forward, here, we report the development of a ribonucleotide version of the dual-app probe to monitor antibiotics-induced conformational changes in RNA. The ribonucleotide analog is derived by conjugating trifluoromethyl-benzofuran at the C5 position of uridine (TFBF-UTP). The analog is efficiently incorporated by T7 RNA polymerase to produce functionalized RNA transcripts. Detailed photophysical and 19F NMR of the nucleoside and nucleotide incorporated into RNA oligonucleotides revealed that the analog is structurally minimally invasive and can be used for probing RNA conformations by fluorescence and 19F NMR techniques. Using the probe, we monitored and estimated aminoglycoside antibiotics binding to the bacterial ribosomal decoding site RNA (A-site, a very important RNA target). While 2-aminopurine, a famous fluorescent nucleic acid probe, fails to detect structurally similar aminoglycoside antibiotics binding to the A-site, our probe reports the binding of different aminoglycosides to the A-site. Taken together, our results demonstrate that TFBF-UTP is a very useful addition to the nucleic acid analysis toolbox and could be used to devise discovery platforms to identify new RNA binders of therapeutic potential.

Capturing heterogeneous conformers of cobalamin riboswitch by cryo-EM

RNA conformational heterogeneity often hampers its high-resolution structure determination, especially for large and flexible RNAs devoid of stabilizing proteins or ligands. The adenosylcobalamin riboswitch exhibits heterogeneous conformations under 1 mM Mg2+ concentration and ligand binding reduces conformational flexibility. Among all conformers, we determined one apo (5.3 Å) and four holo cryo-electron microscopy structures (overall 3.0-3.5 Å, binding pocket 2.9-3.2 Å). The holo dimers exhibit global motions of helical twisting and bending around the dimer interface. A backbone comparison of the apo and holo states reveals a large structural difference in the P6 extension position. The central strand of the binding pocket, junction 6/3, changes from an 'S'- to a 'U'-shaped conformation to accommodate ligand. Furthermore, the binding pocket can partially form under 1 mM Mg2+ and fully form under 10 mM Mg2+ within the bound-like structure in the absence of ligand. Our results not only demonstrate the stabilizing ligand-induced conformational changes in and around the binding pocket but may also provide further insight into the role of the P6 extension in ligand binding and selectivity.

Chemical crosslinking and ligation methods for in vivo analysis of RNA structures and interactions

RNA structures and interactions in living cells drive a variety of biological processes and play critical roles in physiology and disease states. However, studies of RNA structures and interactions have been challenging due to limitations in available technologies. Direct determination of structures in vitro has been only possible to a small number of RNAs with limited sizes and conformations. We recently introduced two chemical crosslink-ligation techniques that enabled studies of transcriptome-wide secondary and tertiary structures and their dynamics. In a dramatically improved version of the psoralen analysis of RNA interactions and structures (PARIS2) method, we detailed the synthesis and use of amotosalen, a highly soluble psoralen analogue, and enhanced enzymology for higher efficiency duplex capture. We also introduced spatial 2'-hydroxyl acylation reversible crosslinking (SHARC) with exonuclease (exo) trimming, a method which utilizes a novel crosslinker class that targets the 2'-OH to capture three-dimensional (3D) structures. Both are powerful orthogonal approaches for solving in vivo RNA structure and interactions, integrating crosslinking, exo trimming, proximity ligation, and high throughput sequencing. In this chapter, we present a detailed protocol for the methods and highlight steps that outperform existing crosslink-ligation approaches.

Visualizing RNA conformational and architectural heterogeneity in solution

RNA flexibility is reflected in its heterogeneous conformation. Through direct visualization using atomic force microscopy (AFM) and the adenosylcobalamin riboswitch aptamer domain as an example, we show that a single RNA sequence folds into conformationally and architecturally heterogeneous structures under near-physiological solution conditions. Recapitulated 3D topological structures from AFM molecular surfaces reveal that all conformers share the same secondary structural elements. Only a population-weighted cohort, not any single conformer, including the crystal structure, can account for the ensemble behaviors observed by small-angle X-ray scattering (SAXS). All conformers except one are functionally active in terms of ligand binding. Our findings provide direct visual evidence that the sequence-structure relationship of RNA under physiologically relevant solution conditions is more complex than the one-to-one relationship for well-structured proteins. The direct visualization of conformational and architectural ensembles at the single-molecule level in solution may suggest new approaches to RNA structural analyses.

Biomotors, viral assembly, and RNA nanobiotechnology: Current achievements and future directions

The International Society of RNA Nanotechnology and Nanomedicine (ISRNN) serves to further the development of a wide variety of functional nucleic acids and other related nanotechnology platforms. To aid in the dissemination of the most recent advancements, a biennial discussion focused on biomotors, viral assembly, and RNA nanobiotechnology has been established where international experts in interdisciplinary fields such as structural biology, biophysical chemistry, nanotechnology, cell and cancer biology, and pharmacology share their latest accomplishments and future perspectives. The results summarized here highlight advancements in our understanding of viral biology and the structure-function relationship of frame-shifting elements in genomic viral RNA, improvements in the predictions of SHAPE analysis of 3D RNA structures, and the understanding of dynamic RNA structures through a variety of experimental and computational means. Additionally, recent advances in the drug delivery, vaccine design, nanopore technologies, biomotor and biomachine development, DNA packaging, RNA nanotechnology, and drug delivery are included in this critical review. We emphasize some of the novel accomplishments, major discussion topics, and present current challenges and perspectives of these emerging fields.

Lessons Learned and Yet-to-Be Learned on the Importance of RNA Structure in SARS-CoV-2 Replication

SARS-CoV-2, the etiological agent responsible for the COVID-19 pandemic, is a member of the virus family Coronaviridae, known for relatively extensive (~30-kb) RNA genomes that not only encode for numerous proteins but are also capable of forming elaborate structures. As highlighted in this review, these structures perform critical functions in various steps of the viral life cycle, ultimately impacting pathogenesis and transmissibility. We examine these elements in the context of coronavirus evolutionary history and future directions for curbing the spread of SARS-CoV-2 and other potential human coronaviruses. While we focus on structures supported by a variety of biochemical, biophysical, and/or computational methods, we also touch here on recent evidence for novel structures in both protein-coding and noncoding regions of the genome, including an assessment of the potential role for RNA structure in the controversial finding of SARS-CoV-2 integration in “long COVID” patients. This review aims to serve as a consolidation of previous works on coronavirus and more recent investigation of SARS-CoV-2, emphasizing the need for improved understanding of the role of RNA structure in the evolution and adaptation of these human viruses.

High-resolution view of HIV-1 reverse transcriptase initiation complexes and inhibition by NNRTI drugs

Reverse transcription of the HIV-1 viral RNA genome (vRNA) is an integral step in virus replication. Upon viral entry, HIV-1 reverse transcriptase (RT) initiates from a host tRNA^Lys₃ primer bound to the vRNA genome and is the target of key antivirals, such as non-nucleoside reverse transcriptase inhibitors (NNRTIs). Initiation proceeds slowly with discrete pausing events along the vRNA template. Despite prior medium-resolution structural characterization of reverse transcriptase initiation complexes (RTICs), higher-resolution structures of the RTIC are needed to understand the molecular mechanisms that underlie initiation. Here we report cryo-EM structures of the core RTIC, RTIC-nevirapine, and RTIC-efavirenz complexes at 2.8, 3.1, and 2.9 Å, respectively. In combination with biochemical studies, these data suggest a basis for rapid dissociation kinetics of RT from the vRNA-tRNA^Lys₃ initiation complex and reveal a specific structural mechanism of nucleic acid conformational stabilization during initiation. Finally, our results show that NNRTIs inhibit the RTIC and exacerbate discrete pausing during early reverse transcription.

Drug Development for Target Ribosomal Protein rpL35/uL29 for Repair of LAMB3R635X in Rare Skin Disease Epidermolysis Bullosa

INTRODUCTION

Epidermolysis bullosa (EB) describes a family of rare genetic blistering skin disorders. Various subtypes are clinically and genetically heterogeneous, and a lethal postpartum form of EB is the generalized severe junctional EB (gs-JEB). gs-JEB is mainly caused by premature termination codon (PTC) mutations in the skin anchor protein LAMB3 (laminin subunit beta-3) gene. The ribosome in majority of translational reads of LAMB3PTC mRNA aborts protein synthesis at the PTC signal, with production of a truncated, nonfunctional protein. This leaves an endogenous readthrough mechanism needed for production of functional full-length Lamb3 protein albeit at insufficient levels. Here, we report on the development of drugs targeting ribosomal protein L35 (rpL35), a ribosomal modifier for customized increase in production of full-length Lamb3 protein from a LAMB3PTC mRNA.

METHODS

Molecular docking studies were employed to identify small molecules binding to human rpL35. Molecular determinants of small molecule binding to rpL35 were further characterized by titration of the protein with these ligands as monitored by nuclear magnetic resonance (NMR) spectroscopy in solution. Changes in NMR chemical shifts were used to map the docking sites for small molecules onto the 3D structure of the rpL35.

RESULTS

Molecular docking studies identified 2 FDA-approved drugs, atazanavir and artesunate, as candidate small-molecule binders of rpL35. Molecular interaction studies predicted several binding clusters for both compounds scattered throughout the rpL35 structure. NMR titration studies identified the amino acids participating in the ligand interaction. Combining docking predictions for atazanavir and artesunate with rpL35 and NMR analysis of rpL35 ligand interaction, one binding cluster located near the N-terminus of rpL35 was identified. In this region, the nonidentical binding sites for atazanavir and artesunate overlap and are accessible when rpL35 is integrated in its natural ribosomal environment.

CONCLUSION

Atazanavir and artesunate were identified as candidate compounds binding to ribosomal protein rpL35 and may now be tested for their potential to trigger a rpL35 ribosomal switch to increase production of full-length Lamb3 protein from a LAMB3PTC mRNA for targeted systemic therapy in treating gs-JEB.

FRET-based dynamic structural biology: Challenges, perspectives and an appeal for open-science practices

Single-molecule FRET (smFRET) has become a mainstream technique for studying biomolecular structural dynamics. The rapid and wide adoption of smFRET experiments by an ever-increasing number of groups has generated significant progress in sample preparation, measurement procedures, data analysis, algorithms and documentation. Several labs that employ smFRET approaches have joined forces to inform the smFRET community about streamlining how to perform experiments and analyze results for obtaining quantitative information on biomolecular structure and dynamics. The recent efforts include blind tests to assess the accuracy and the precision of smFRET experiments among different labs using various procedures. These multi-lab studies have led to the development of smFRET procedures and documentation, which are important when submitting entries into the archiving system for integrative structure models, PDB-Dev. This position paper describes the current 'state of the art' from different perspectives, points to unresolved methodological issues for quantitative structural studies, provides a set of 'soft recommendations' about which an emerging consensus exists, and lists openly available resources for newcomers and seasoned practitioners. To make further progress, we strongly encourage 'open science' practices.

A Bayes-inspired theory for optimally building an efficient coarse-grained folding force field

Because of their potential utility in predicting conformational changes and assessing folding dynamics, coarse-grained (CG) RNA folding models are appealing for rapid characterization of RNA molecules. Previously, we reported the iterative simulated RNA reference state (IsRNA) method for parameterizing a CG force field for RNA folding, which consecutively updates the simulation force field to reflect marginal distributions of folding coordinates in the structure database and extract various energy terms. While the IsRNA model was validated by showing close agreement between the IsRNA-simulated and experimentally observed distributions, here, we expand our theoretical understanding of the model and, in doing so, improve the parameterization process to optimize the subset of included folding coordinates, which leads to accelerated simulations. Using statistical mechanical theory, we analyze the underlying, Bayesian concept that drives parameterization of the energy function, providing a general method for developing predictive, knowledge-based, polymer force fields on the basis of limited data. Furthermore, we propose an optimal parameterization procedure, based on the principal of maximum entropy.

Structural Heterogeneities of the Ribosome: New Frontiers and Opportunities for Cryo-EM

The extent of ribosomal heterogeneity has caught increasing interest over the past few years, as recent studies have highlighted the presence of structural variations of the ribosome. More precisely, the heterogeneity of the ribosome covers multiple scales, including the dynamical aspects of ribosomal motion at the single particle level, specialization at the cellular and subcellular scale, or evolutionary differences across species. Upon solving the ribosome atomic structure at medium to high resolution, cryogenic electron microscopy (cryo-EM) has enabled investigating all these forms of heterogeneity. In this review, we present some recent advances in quantifying ribosome heterogeneity, with a focus on the conformational and evolutionary variations of the ribosome and their functional implications. These efforts highlight the need for new computational methods and comparative tools, to comprehensively model the continuous conformational transition pathways of the ribosome, as well as its evolution. While developing these methods presents some important challenges, it also provides an opportunity to extend our interpretation and usage of cryo-EM data, which would more generally benefit the study of molecular dynamics and evolution of proteins and other complexes.

Deleterious effects of carbon-carbon dipolar coupling on RNA NMR dynamics

Many regulatory RNAs undergo dynamic exchanges that are crucial for their biological functions and NMR spectroscopy is a versatile tool for monitoring dynamic motions of biomolecules. Meaningful information on biomolecular dynamics requires an accurate measurement of relaxation parameters such as longitudinal (R₁) rates, transverse (R₂) rates and heteronuclear Overhauser effect (hNOE). However, earlier studies have shown that the large ¹³C-¹³C interactions complicate analysis of the carbon relaxation parameters. To investigate the effect of ¹³C-¹³C interactions on RNA dynamic studies, we performed relaxation measurements on various RNA samples with different labeling patterns and compared these measurements with the computational simulations. For uniformly labeled samples, contributions of the neighboring carbon to R₁ measurements were observed. These contributions increased with increasing magnetic field and overall correlation time ([Formula: see text]) for R₁ rates, necessitating more careful analysis for uniformly labeled large RNAs. In addition, the hNOE measurements were also affected by the adjacent carbon nuclei. Unlike R₁ rates, R_1ρ rates showed relatively good agreement between uniformly- and site-selectively labeled samples, suggesting no dramatic effect from their attached carbon, in agreement with previous observations. Overall, having more accurate rate measurements avoids complex analysis and will be a key for interpreting ¹³C relaxation rates for molecular motion that can provide valuable insights into cellular molecular recognition events.

Distinct Conformational States Underlie Pausing during Initiation of HIV-1 Reverse Transcription

A hallmark of the initiation step of HIV-1 reverse transcription, in which viral RNA genome is converted into double-stranded DNA, is that it is slow and non-processive. Biochemical studies have identified specific sites along the viral RNA genomic template in which reverse transcriptase (RT) stalls. These stalling points, which occur after the addition of three and five template dNTPs, may serve as checkpoints to regulate the precise timing of HIV-1 reverse transcription following viral entry. Structural studies of reverse transcriptase initiation complexes (RTICs) have revealed unique conformations that may explain the slow rate of incorporation; however, questions remain about the temporal evolution of the complex and features that contribute to strong pausing during initiation. Here we present cryo-electron microscopy and single-molecule characterization of an RTIC after three rounds of dNTP incorporation (+3), the first major pausing point during reverse transcription initiation. Cryo-electron microscopy structures of a +3 extended RTIC reveal conformational heterogeneity within the RTIC core. Three distinct conformations were identified, two of which adopt unique, likely off-pathway, intermediates in the canonical polymerization cycle. Single-molecule Förster resonance energy transfer experiments confirm that the +3 RTIC is more structurally dynamic than earlier-stage RTICs. These alternative conformations were selectively disrupted through structure-guided point mutations to shift single-molecule Förster resonance energy transfer populations back toward the on-pathway conformation. Our results support the hypothesis that conformational heterogeneity within the HIV-1 RTIC during pausing serves as an additional means of regulating HIV-1 replication.

NMR Chemical Exchange Measurements Reveal That N6-Methyladenosine Slows RNA Annealing

N⁶-Methyladenosine (m⁶A) is an abundant epitranscriptomic modification that plays important roles in many aspects of RNA metabolism. While m⁶A is thought to mainly function by recruiting reader proteins to specific RNA sites, the modification can also reshape RNA-protein and RNA-RNA interactions by altering RNA structure mainly by destabilizing base pairing. Little is known about how m⁶A and other epitranscriptomic modifications might affect the kinetic rates of RNA folding and other conformational transitions that are also important for cellular activity. Here, we used NMR R_1ρ relaxation dispersion and chemical exchange saturation transfer to noninvasively and site-specifically measure nucleic acid hybridization kinetics. The methodology was validated on two DNA duplexes and then applied to examine how a single m⁶A alters the hybridization kinetics in two RNA duplexes. The results show that m⁶A minimally impacts the rate constant for duplex dissociation, changing k_off by ∼1-fold but significantly slows the rate of duplex annealing, decreasing k_on by ∼7-fold. A reduction in the annealing rate was observed robustly for two different sequence contexts at different temperatures, both in the presence and absence of Mg²⁺. We propose that rotation of the N⁶-methyl group from the preferred syn conformation in the unpaired nucleotide to the energetically disfavored anti conformation required for Watson-Crick pairing is responsible for the reduced annealing rate. The results help explain why in mRNA m⁶A slows down tRNA selection and more generally suggest that m⁶A may exert cellular functions by reshaping the kinetics of RNA conformational transitions.

Refining RNA solution structures with the integrative use of label-free paramagnetic relaxation enhancement NMR

NMR structure calculation is inherently integrative, and can incorporate new experimental data as restraints. As RNAs have lower proton densities and are more conformational heterogenous than proteins, the refinement of RNA structures can benefit from additional types of restraints. Paramagnetic relaxation enhancement (PRE) provides distance information between a paramagnetic probe and protein or RNA nuclei. However, covalent conjugation of a paramagnetic probe is difficult for RNAs, thus limiting the use of PRE NMR for RNA structure characterization. Here, we show that the solvent PRE can be accurately measured for RNA labile imino protons, simply with the addition of an inert paramagnetic cosolute. Demonstrated on three RNAs that have increasingly complex topologies, we show that the incorporation of the solvent PRE restraints can significantly improve the precision and accuracy of RNA structures. Importantly, the solvent PRE data can be collected for RNAs without isotope enrichment. Thus, the solvent PRE method can work integratively with other biophysical techniques for better characterization of RNA structures.