Integrated NMR/Molecular Dynamics Determination of the Ensemble Conformation of a Thermodynamically Stable CUUG RNA Tetraloop

Insight into the conformational ensembles formed by U-U and T-T mismatches in RNA and DNA duplexes from a structure-based survey, NMR, and molecular dynamics simulations

Nucleic acid base pairs interconvert between alternative conformations on a free energy landscape, and these dynamics play critical roles in recognition, folding, and catalysis. U-U and T-T mismatches can adopt two nearly isoenergetic wobble conformations, distinguished by their relative shearing displacements. Experimental NMR evidence suggests that these conformations dynamically interconvert in RNA motifs containing tandem U-U mismatches. However, whether such motions occur ubiquitously across U-U and T-T mismatches remains unknown, as high-resolution nucleic acid structures typically report only a single conformation. Here, we used NMR spectroscopy, a structure-based survey of the Protein Data Bank, and molecular dynamics (MD) simulations to investigate wobble dynamics in U-U and T-T mismatches when flanked by canonical Watson-Crick base pairs in RNA and DNA duplexes. The structure-based survey revealed that U-U mismatches have propensities to adopt alternative wobble conformations even when controlling for sequence and identified potential intermediates along the wobble transition. Off-resonance R1ρ relaxation dispersion experiments detected no micro- to millisecond dynamics for U-U mismatches in duplex RNA and T-T mismatches in duplex DNA. However, alternative conformer refinement of the electron density in X-ray structures, inter-proton NOEs, carbonyl carbon chemical shifts, an RDC-derived conformational ensemble, and MD simulations indicated that U-U and T-T mismatches exist in a dynamic equilibrium between two wobble conformations, with the minor state exceeding 30% and the transitions occurring on the nanosecond timescale. Our findings suggest that U-U and T-T ubiquitously undergo sub-microsecond wobble motions, contributing to the energetic landscape and dynamic plasticity of nucleic acids, with important implications for processes that generate and act on these mismatches.

Optimization of Structure-Guided Development of Chemical Probes for the Pseudoknot RNA of the Frameshift Element in SARS-CoV-2

Targeting the RNA genome of SARS-CoV-2 is a viable option for antiviral drug development. We explored three ligand binding sites of the core pseudoknot RNA of the SARS-CoV-2 frameshift element. We iteratively optimized ligands, based on improved affinities, targeting these binding sites and report on structural and dynamic properties of the three identified binding sites. Available experimental 3D structures of the pseudoknot element were compared to SAXS and NMR data to validate its dominant folding state in solution. In order to experimentally map in silico predicted binding sites, NMR assignments of the majority of nucleobases were achieved by segmental labeling of the pseudoknot RNA and isotope-filtered NMR experiments at 1.2 GHz, demonstrating the value of NMR spectroscopy to supplement modelling and docking data. Optimized ligands with enhanced affinity were shown to specifically inhibit frameshifting without affecting 0-frame translation in cell-free translation assays, establishing the frameshift element as target for drug-like ligands of low molecular weight.

Inferring phase transitions and critical exponents from limited observations with thermodynamic maps

Phase transitions are ubiquitous across life, yet hard to quantify and describe accurately. In this work, we develop an approach for characterizing generic attributes of phase transitions from very limited observations made deep within different phases' domains of stability. Our approach is called thermodynamic maps (TM), which combines statistical mechanics and molecular simulations with score-based generative models. TM enable learning the temperature dependence of arbitrary thermodynamic observables across a wide range of temperatures. We show its usefulness by calculating phase transition attributes such as melting temperature, temperature-dependent heat capacities, and critical exponents. For instance, we demonstrate the ability of TM to infer the ferromagnetic phase transition of the Ising model, including temperature-dependent heat capacity and critical exponents, despite never having seen samples from the transition region. In addition, we efficiently characterize the temperature-dependent conformational ensemble and compute melting curves of the two RNA systems: a GCAA tetraloop and the HIV-TAR RNA, which are notoriously hard to sample due to glassy-like energy landscapes.

Nexus between RNA conformational dynamics and functional versatility

RNA conformational dynamics is pivotal for functional regulations in biology. RNA can function as versatile as protein but adopts multiple distinct structures. In this review, we provide a focused review of the recent advances in studies of RNA conformational dynamics and address some of the misconceptions about RNA structure and its conformational dynamics. We discuss why the traditional methods for structure determination come up short in describing RNA conformational space. The examples discussed provide illustrations of the structure-based mechanisms of RNAs with diverse roles, including viral, long noncoding, and catalytic RNAs, one of which focuses on the debated area of conformational heterogeneity of an RNA structural element in the HIV-1 genome.

Conformational penalties: New insights into nucleic acid recognition

The energy cost accompanying changes in the structures of nucleic acids when they bind partner molecules is a significant but underappreciated thermodynamic contribution to binding affinity and specificity. This review highlights recent advances in measuring conformational penalties and determining their contribution to the recognition, folding, and regulatory activities of nucleic acids. Notable progress includes methods for measuring and structurally characterizing lowly populated conformational states, obtaining ensemble information in high throughput, for large macromolecular assemblies, and in complex cellular environments. Additionally, quantitative and predictive thermodynamic models have been developed that relate conformational penalties to nucleic acid-protein association and cellular activity. These studies underscore the crucial role of conformational penalties in nucleic acid recognition.

Dissecting the Conformational Heterogeneity of Stem-Loop Substructures of the Fifth Element in the 5'-Untranslated Region of SARS-CoV-2

Throughout the family of coronaviruses, structured RNA elements within the 5' region of the genome are highly conserved. The fifth stem-loop element from SARS-CoV-2 (5_SL5) represents an example of an RNA structural element, repeatedly occurring in coronaviruses. It contains a conserved, repetitive fold within its substructures SL5a and SL5b. We herein report the detailed characterization of the structure and dynamics of elements SL5a and SL5b that are located immediately upstream of the SARS-CoV-2 ORF1a/b start codon. Exploiting the unique ability of solution NMR methods, we show that the structures of both apical loops are modulated by structural differences in the remote parts located in their stem regions. We further integrated our high-resolution models of SL5a/b into the context of full-length 5_SL5 structures by combining different structural biology methods. Finally, we evaluated the impact of the two most common VoC mutations within 5_SL5 with respect to individual base-pair stability.

Synthesis and Application of D- and 13C-Labelled tert-Butyl Hoechst Dye

Herein, the successful syntheses of D3- and 13C-N-methyl and D9-tert-butyl Hoechst dyes are presented. This includes the preparation of the labelled D3- and 13C-N-methyl piperazines and D9-tert-butylated hydroxytoluene precursors. The tert-butyl Hoechst dye is known to bind a specific RNA aptamer. Spectroscopic NMR studies of the labelled Hoechst dye-aptamer complexes allowed for the unambiguous assignment of chemical shifts, as well as the dynamics of the bound dye.

Insights into the A-C Mismatch Conformational Ensemble in Duplex DNA and its Role in Genetic Processes through a Structure-based Review

Knowing the conformational ensembles formed by mismatches is crucial for understanding how they are generated and repaired and how they contribute to genomic instability. Here, we review structural and energetic studies of the A-C mismatch in duplex DNA and use the information to identify critical conformational states in its ensemble and their significance in genetic processes. In the 1970s, Topal and Fresco proposed the A-C wobble stabilized by two hydrogen bonds, one requiring protonation of adenine-N1. Subsequent NMR and X-ray crystallography studies showed that the protonated A-C wobble was in dynamic equilibrium with a neutral inverted wobble. The mismatch was shown to destabilize duplex DNA in a sequence- and pH-dependent manner by 2.4-3.8 kcal/mol and to have an apparent pKa ranging between 7.2 and 7.7. The A-C mismatch conformational repertoire expanded as structures were determined for damaged and protein-bound DNA. These structures included Watson-Crick-like conformations forming through tautomerization of the bases that drive replication errors, the reverse wobble forming through rotation of the entire nucleotide proposed to increase the fidelity of DNA replication, and the Hoogsteen base-pair forming through the flipping of the adenine base which explained the unusual specificity of DNA polymerases that bypass DNA damage. Thus, the A-C mismatch ensemble encompasses various conformational states that can be selectively stabilized in response to environmental changes such as pH shifts, intermolecular interactions, and chemical modifications, and these adaptations facilitate critical biological processes. This review also highlights the utility of existing 3D structures to build ensemble models for nucleic acid motifs.

RNA dynamics from experimental and computational approaches

Conformational dynamics is crucial for the biological function of RNA molecules and for their potential as therapeutic targets. This meeting report outlines key "take-home" messages that emerged from the presentations and discussions during the CECAM workshop "RNA dynamics from experimental and computational approaches" in Paris, June 26-28, 2023.

Characterization of SLA RNA promoter from dengue virus and its interaction with the viral non-structural NS5 protein

The Dengue virus (DENV) is the most significant arthropod-borne viral pathogen in humans with 400 million infections annually. DENV comprises four distinct serotypes (DENV-1 to -4) which complicates vaccine development. Any of the four serotypes can cause clinical illness but with distinctive infection dynamics. Variations in sequences identified within the four genomes induce structural differences in crucial RNA motifs that were suggested to be correlated to the degree of pathogenicity among DENV-1 to -4. In particular, the RNA Stem-loop A (SLA) at the 5'-end of the genome, acts as a key regulator of the viral replication cycle by interacting with the viral NS5 polymerase to initiate the minus-strand viral RNA synthesis and later to methylate and cap the synthesized RNA. The molecular details of this interaction remain not fully described. Here, we report the solution secondary structures of SLA from DENV-1 to -4. Our results highlight that the four SLA exhibit structural and dynamic differences. Secondly, to determine whether SLA RNA contains serotype-specific determinants for the recognition by the viral NS5 protein, we investigated interactions between SLA from DENV -1 to -4 and DENV2 NS5 using combined biophysical approaches. Our results show that NS5 from DENV2 is able to bind SLA from other serotypes, but that other viral or host factors may be necessary to stabilize the complex and promote the catalytically active state of the NS5. By contrast, we show that a serotype-specific binding is driven by specific interactions involving conformational changes within the SLA RNA.

NMR characterization and ligand binding site of the stem-loop 2 motif from the Delta variant of SARS-CoV-2

The stem-loop 2 motif (s2m) in SARS-CoV-2 (SCoV-2) is located in the 3'-UTR. Although s2m has been reported to display characteristics of a mobile genomic element that might lead to an evolutionary advantage, its function has remained unknown. The secondary structure of the original SCoV-2 RNA sequence (Wuhan-Hu-1) was determined by NMR in late 2020, delineating the base-pairing pattern and revealing substantial differences in secondary structure compared to SARS-CoV-1 (SCoV-1). The existence of a single G29742-A29756 mismatch in the upper stem of s2m leads to its destabilization and impedes a complete NMR analysis. With Delta, a variant of concern has evolved with one mutation compared to the original sequence that replaces G29742 by U29742. We show here that this mutation results in a more defined structure at ambient temperature accompanied by a rise in melting temperature. Consequently, we were able to identify >90% of the relevant NMR resonances using a combination of selective RNA labeling and filtered 2D NOESY as well as 4D NMR experiments. We present a comprehensive NMR analysis of the secondary structure, (sub)nanosecond dynamics, and ribose conformation of s2m Delta based on heteronuclear 13C NOE and T 1 measurements and ribose carbon chemical shift-derived canonical coordinates. We further show that the G29742U mutation in Delta has no influence on the druggability of s2m compared to the Wuhan-Hu-1 sequence. With the assignment at hand, we identify the flexible regions of s2m as the primary site for small molecule binding.

Currently Used Methods to Evaluate the Efficacy of Therapeutic Drugs and Kidney Safety

Kidney diseases with kidney failure or damage, such as chronic kidney disease (CKD) and acute kidney injury (AKI), are common clinical problems worldwide and have rapidly increased in prevalence, affecting millions of people in recent decades. A series of novel diagnostic or predictive biomarkers have been discovered over the past decade, enhancing the investigation of renal dysfunction in preclinical studies and clinical risk assessment for humans. Since multiple causes lead to renal failure, animal studies have been extensively used to identify specific disease biomarkers for understanding the potential targets and nephropathy events in therapeutic insights into disease progression. Mice are the most commonly used model to investigate the mechanism of human nephropathy, and the current alternative methods, including in vitro and in silico models, can offer quicker, cheaper, and more effective methods to avoid or reduce the unethical procedures of animal usage. This review provides modern approaches, including animal and nonanimal assays, that can be applied to study chronic nonclinical safety. These specific situations could be utilized in nonclinical or clinical drug development to provide information on kidney disease.