Prevalence of syn nucleobases in the active sites of functional RNAs

Structure of HIV-1 RRE stem-loop II identifies two conformational states of the high-affinity Rev binding site

During HIV infection, specific RNA-protein interaction between the Rev response element (RRE) and viral Rev protein is required for nuclear export of intron-containing viral mRNA transcripts. Rev initially binds the high-affinity site in stem-loop II, which promotes oligomerization of additional Rev proteins on RRE. Here, we present the crystal structure of RRE stem-loop II in distinct closed and open conformations. The high-affinity Rev-binding site is located within the three-way junction rather than the predicted stem IIB. The closed and open conformers differ in their non-canonical interactions within the three-way junction, and only the open conformation has the widened major groove conducive to initial Rev interaction. Rev binding assays show that RRE stem-loop II has high- and low-affinity binding sites, each of which binds a Rev dimer. We propose a binding model, wherein Rev-binding sites on RRE are sequentially created through structural rearrangements induced by Rev-RRE interactions.

High-resolution RNA tertiary structures in Zika virus stem-loop A for the development of inhibitory small molecules

Flaviviruses such as Zika (ZIKV) and dengue virus (DENV) are positive-sense RNA viruses belonging to Flaviviridae The flavivirus genome contains a 5' end stem-loop promoter sequence known as stem-loop A (SLA) that is recognized by the flavivirus polymerase NS5 during viral RNA synthesis and 5' guanosine cap methylation. The crystal structures of ZIKV and DENV SLAs show a well-defined fold, consisting of a bottom stem, side loop, and top stem-loop, providing unique interaction sites for small molecule inhibitors to disrupt the promoter function. To facilitate the identification of small molecule binding sites in flavivirus SLA, we determined high-resolution structures of the bottom and top stems of ZIKV SLA, which contain a single U- or G-bulge, respectively. Both bulge nucleotides exhibit multiple orientations, from folded back on the adjacent nucleotide to flipped out of the helix, and are stabilized by stacking or base triple interactions. These structures suggest that even a single unpaired nucleotide can provide flexibility to RNA structures, and its conformation is mainly determined by the stabilizing chemical environment. To facilitate discovery of small molecule inhibitors that interfere with the functions of ZIKV SLA, we screened and identified compounds that bind to the bottom and top stems of ZIKV SLA.

Conformational dynamics of RNA G4C2 and G2C4 repeat expansions causing ALS/FTD using NMR and molecular dynamics studies

G4C2 and G2C4 repeat expansions in chromosome 9 open reading frame 72 (C9orf72) are the most common cause of genetically defined amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD), or c9ALS/FTD. The gene is bidirectionally transcribed, producing G4C2 repeats [r(G4C2)exp] and G2C4 repeats [r(G2C4)exp]. The c9ALS/FTD repeat expansions are highly structured, and structural studies showed that r(G4C2)exp predominantly folds into a hairpin with a periodic array of 1 × 1 G/G internal loops and a G-quadruplex. A small molecule probe revealed that r(G4C2)exp also adopts a hairpin structure with 2 × 2 GG/GG internal loops. We studied the conformational dynamics adopted by 2 × 2 GG/GG loops using temperature replica exchange molecular dynamics (T-REMD) and further characterized the structure and underlying dynamics using traditional 2D NMR techniques. These studies showed that the loop's closing base pairs influence both structure and dynamics, particularly the configuration adopted around the glycosidic bond. Interestingly, r(G2C4) repeats, which fold into an array of 2 × 2 CC/CC internal loops, are not as dynamic. Collectively, these studies emphasize the unique sensitivity of r(G4C2)exp to small changes in stacking interactions, which is not observed in r(G2C4)exp, providing important considerations for further principles in structure-based drug design.

Structural and functional diversity of bacterial cyclic nucleotide perception by CRP proteins

Cyclic AMP (cAMP) is a ubiquitous second messenger synthesized by most living organisms. In bacteria, it plays highly diverse roles in metabolism, host colonization, motility, and many other processes important for optimal fitness. The main route of cAMP perception is through transcription factors from the diverse and versatile CRP-FNR protein superfamily. Since the discovery of the very first CRP protein CAP in Escherichia coli more than four decades ago, its homologs have been characterized in both closely related and distant bacterial species. The cAMP-mediated gene activation for carbon catabolism by a CRP protein in the absence of glucose seems to be restricted to E. coli and its close relatives. In other phyla, the regulatory targets are more diverse. In addition to cAMP, cGMP has recently been identified as a ligand of certain CRP proteins. In a CRP dimer, each of the two cyclic nucleotide molecules makes contacts with both protein subunits and effectuates a conformational change that favors DNA binding. Here, we summarize the current knowledge on structural and physiological aspects of E. coli CAP compared with other cAMP- and cGMP-activated transcription factors, and point to emerging trends in metabolic regulation related to lysine modification and membrane association of CRP proteins.

G-QINDER Tool: Bioinformatically Predicted Formation of Different Four-Stranded DNA Motifs from (GT)n and (GA)n Repeats

The recently introduced semi-orthogonal system of nucleic acid imaging offers a greatly improved method of identifying DNA sequences that are capable of adopting noncanonical structures. This paper uses our newly developed G-QINDER tool to identify specific repeat sequences that adopt unique structural motifs in DNA: TG and AG repeats. The structures were found to adopt a left-handed G-quadruplex form under extreme crowding conditions and a unique tetrahelical motif under certain other conditions. The tetrahelical structure likely consists of stacked AGAG-tetrads but, unlike G-quadruplexes, their stability does not appear to be dependent on the type of monovalent cation present. The occurrence of TG and AG repeats in genomes is not rare, and they are also found frequently in the regulatory regions of nucleic acids, so it is reasonable to assume that putative structural motifs, like other noncanonical forms, could play an important regulatory role in cells. This hypothesis is supported by the structural stability of the AGAG motif; its unfolding can occur even at physiological temperatures since the melting temperature is primarily dependent on the number of AG repeats in the sequence.

Rare ribosomal RNA sequences from archaea stabilize the bacterial ribosome

The ribosome serves as the universally conserved translator of the genetic code into proteins and supports life across diverse temperatures ranging from below freezing to above 120°C. Ribosomes are capable of functioning across this wide range of temperatures even though the catalytic site for peptide bond formation, the peptidyl transferase center, is nearly universally conserved. Here we find that Thermoproteota, a phylum of thermophilic Archaea, substitute cytidine for uridine at large subunit rRNA positions 2554 and 2555 (Escherichia coli numbering) in the A loop, immediately adjacent to the binding site for the 3'-end of A-site tRNA. We show by cryo-EM that E. coli ribosomes with uridine to cytidine mutations at these positions retain the proper fold and post-transcriptional modification of the A loop. Additionally, these mutations do not affect cellular growth, protect the large ribosomal subunit from thermal denaturation, and increase the mutational robustness of nucleotides in the peptidyl transferase center. This work identifies sequence variation across archaeal ribosomes in the peptidyl transferase center that likely confers stabilization of the ribosome at high temperatures and develops a stable mutant bacterial ribosome that can act as a scaffold for future ribosome engineering efforts.

Lone Pair…π Contacts and Structure Signatures of r(UNCG) Tetraloops, Z-Turns, and Z-Steps: A WebFR3D Survey

Z-DNA and Z-RNA have long appeared as oddities to nucleic acid scientists. However, their Z-step constituents are recurrently observed in all types of nucleic acid systems including ribosomes. Z-steps are NpN steps that are isostructural to Z-DNA CpG steps. Among their structural features, Z-steps are characterized by the presence of a lone pair…π contact that involves the stacking of the ribose O4′ atom of the first nucleotide with the 3′-face of the second nucleotide. Recently, it has been documented that the CpG step of the ubiquitous r(UNCG) tetraloops is a Z-step. Accordingly, such r(UNCG) conformations were called Z-turns. It has also been recognized that an r(GAAA) tetraloop in appropriate conditions can shapeshift to an unusual Z-turn conformation embedding an ApA Z-step. In this report, we explore the multiplicity of RNA motifs based on Z-steps by using the WebFR3D tool to which we added functionalities to be able to retrieve motifs containing lone pair…π contacts. Many examples that underscore the diversity and universality of these motifs are provided as well as tutorial guidance on using WebFR3D. In addition, this study provides an extensive survey of crystallographic, cryo-EM, NMR, and molecular dynamics studies on r(UNCG) tetraloops with a critical view on how to conduct database searches and exploit their results.

Evaluating Geometric Definitions of Stacking for RNA Dinucleoside Monophosphates Using Molecular Mechanics Calculations

RNA modulation via small molecules is a novel approach in pharmacotherapies, where the determination of the structural properties of RNA motifs is considered a promising way to develop drugs capable of targeting RNA structures to control diseases. However, due to the complexity and dynamic nature of RNA molecules, the determination of RNA structures using experimental approaches is not always feasible, and computational models employing force fields can provide important insight. The quality of the force field will determine how well the predictions are compared to experimental observables. Stacking in nucleic acids is one such structural property, originating mainly from London dispersion forces, which are quantum mechanical and are included in molecular mechanics force fields through nonbonded interactions. Geometric descriptions are utilized to decide if two residues are stacked and hence to calculate the stacking free energies for RNA dinucleoside monophosphates (DNMPs) through statistical mechanics for comparison with experimental thermodynamics data. Here, we benchmark four different stacking definitions using molecular dynamics (MD) trajectories for 16 RNA DNMPs produced by two different force fields (RNA-IL and ff99OL3) and show that our stacking definition better correlates with the experimental thermodynamics data. While predictions within an accuracy of 0.2 kcal/mol at 300 K were observed in RNA CC, CU, UC, AG, GA, and GG, stacked states of purine-pyrimidine and pyrimidine-purine DNMPs, respectively, were typically underpredicted and overpredicted. Additionally, population distributions of RNA UU DNMPs were poorly predicted by both force fields, implying a requirement for further force field revisions. We further discuss the differences predicted by each RNA force field. Finally, we show that discrete path sampling (DPS) calculations can provide valuable information and complement the MD simulations. We propose the use of experimental thermodynamics data for RNA DNMPs as benchmarks for testing RNA force fields.

Multiscale computational investigations of the translesion synthesis bypass of tobacco-derived DNA adducts: critical insights that complement experimental biochemical studies

Among the numerous agents that damage DNA, tobacco products remain one of the most lethal and result in the most diverse set of DNA lesions. This perspective aims to provide an overview of computational work conducted to complement experimental biochemical studies on the mutagenicity of adducts derived from the most potent tobacco carcinogen, namely 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone (nicotine-derived nitrosaminoketone or NNK). Lesions ranging from the smallest methylated thymine derivatives to the larger, flexible pyridyloxobutyl (POB) guanine adducts are considered. Insights are obtained from density functional theory (DFT) calculations and molecular dynamics (MD) simulations into the damaged nucleobase and nucleoside structures, the accommodation of the lesions in the active site of key human polymerases, the intrinsic base pairing potentials of the adducts, and dNTP incorporation opposite the lesions. Overall, the computational data provide atomic level information that can rationalize the differential mutagenic properties of tobacco-derived lesions and uncover important insights into the impact of adduct size, nucleobase, position, and chemical composition of the bulky moiety.

MA'AT Analysis of Aldofuranosyl Rings: Unbiased Modeling of Conformational Equilibria and Dynamics in Solution

MA'AT analysis has been applied to methyl β-d-ribofuranoside (3) and methyl 2-deoxy-β-d-erythro-pentofuranoside (4) to demonstrate the ability of this new experimental method to determine multi-state conformational equilibria in solution. Density functional theory (DFT) was used to obtain parameterized equations for >20 NMR spin-coupling constants sensitive to furanose ring conformation in 3 and 4, and these equations were used in conjunction with experimental spin-couplings to produce unbiased MA'AT models of ring pseudorotation. These models describe two-state north-south conformational exchange consistent with results obtained from traditional treatments of more limited sets of NMR spin-couplings (e.g., PSEUROT). While PSEUROT, MA'AT, and aqueous molecular dynamics models yielded similar two-state models, MA'AT analysis gives more reliable results since significantly more experimental observables are employed compared to PSEUROT, and no assumptions are needed to render the fitting tractable. MA'AT models indicate a roughly equal distribution of north and south ring conformers of 4 in aqueous (²H₂O) solution compared to ∼80% north forms for 3. Librational motion about the mean pseudorotation phase angles P of the preferred north and south conformers of 3 in solution is more constrained than that for 4. The greater rigidity of the β-ribo ring may be caused by synergistic stereoelectronic effects and/or noncovalent (e.g., hydrogen-bonding) interactions in solution that preferentially stabilize north forms of 3. MA'AT analysis of oligonucleotides and other furanose ring-containing biomolecules promises to improve current experimental models of sugar ring behavior in solution and help reveal context effects on ring conformation in more complex biologically important systems.

Nuclear Magnetic Resonance Spectra and AMBER OL3 and ROC-RNA Simulations of UCUCGU Reveal Force Field Strengths and Weaknesses for Single-Stranded RNA

Single-stranded regions of RNA are important for folding of sequences into 3D structures and for design of therapeutics targeting RNA. Prediction of ensembles of 3D structures for single-stranded regions often involves classical mechanical approximations of interactions defined by quantum mechanical calculations on small model systems. Nuclear magnetic resonance (NMR) spectra and molecular dynamics (MD) simulations of short single strands provide tests for how well the approximations model many of the interactions. Here, the NMR spectra for UCUCGU at 2, 15, and 30 °C are compared to simulations with the AMBER force fields, OL3 and ROC-RNA. This is the first such comparison to an oligoribonucleotide containing an internal guanosine nucleotide (G). G is particularly interesting because of its many H-bonding groups, large dipole moment, and proclivity for both syn and anti conformations. Results reveal formation of a G amino to phosphate non-bridging oxygen H-bond. The results also demonstrate dramatic differences in details of the predicted structures. The variations emphasize the dependence of predictions on individual parameters and their balance with the rest of the force field. The NMR data can serve as a benchmark for future force fields.

Metabolomic Profiling of Human Urine Samples Using LC-TIMS-QTOF Mass Spectrometry

The identification of metabolites in biological samples is challenging due to their chemical and structural diversity. Ion mobility spectrometry (IMS) separates ionized molecules based on their mobility in a carrier buffer gas giving information about the ionic shape by measuring the rotationally averaged collision cross-section (CCS) value. This orthogonal descriptor, in combination with the m/z, isotopic pattern distribution, and MS/MS spectrum, has the potential to improve the identification of molecular molecules in complex mixtures. Urine metabolomics can reveal metabolic differences, which arise as a result of a specific disease or in response to therapeutic intervention. It is, however, complicated by the presence of metabolic breakdown products derived from a wide range of lifestyle and diet-related byproducts, many of which are poorly characterized. In this study, we explore the use of trapped ion mobility spectrometry (TIMS) via LC parallel accumulation with serial fragmentation (PASEF) for urine metabolomics. A total of 362 urine metabolites were characterized from 80 urine samples collected from healthy volunteers using untargeted metabolomics employing HILIC and RP chromatography. Additionally, three analytes (Trp, Phe, and Tyr) were selected for targeted quantification. Both the untargeted and targeted data was highly reproducible and reported CCS measurements for identified metabolites were robust in the presence of the urine matrix. A comparison of CCS values among different laboratories was also conducted, showing less than 1.3% ΔCCS values across different platforms. This is the first report of a human urine metabolite database compiled with CCS values experimentally acquired using an LC-PASEF TIMS-qTOF platform.

Conformational energies and equilibria of cyclic dinucleotides in vacuo and in solution: computational chemistry vs. NMR experiments

Performance of computational methods in modelling cyclic dinucleotides - an important and challenging class of compounds - has been evaluated by two different benchmarks: (1) gas-phase conformational energies and (2) qualitative agreement with NMR observations of the orientation of the χ-dihedral angle in solvent. In gas-phase benchmarks, where CCSD(T) and DLPNO-CCSD(T) methods have been used as the reference, most of the (dispersion corrected) density functional approximations are accurate enough to justify prioritizing computational cost and compatibility with other modelling options as the criterion of choice. NMR experiments of 3'3'-c-di-AMP, 3'3'-c-GAMP, and 3'3'-c-di-GMP show the overall prevalence of the anti-conformation of purine bases, but some population of syn-conformations is observed for guanines. Implicit solvation models combined with quantum-chemical methods struggle to reproduce this behaviour, probably due to a lack of dynamics and explicitly modelled solvent, leading to structures that are too compact. Molecular dynamics simulations overrepresent the syn-conformation of guanine due to the overestimation of an intramolecular hydrogen bond. Our combination of experimental and computational benchmarks provides "error bars" for modelling cyclic dinucleotides in solvent, where such information is generally difficult to obtain, and should help gauge the interpretability of studies dealing with binding of cyclic dinucleotides to important pharmaceutical targets. At the same time, the presented analysis calls for improvement in both implicit solvation models and force-field parameters.

The model structure of the hammerhead ribozyme formed by RNAs of reciprocal chirality

RNA-based tools are frequently used to modulate gene expression in living cells. However, the stability and effectiveness of such RNA-based tools is limited by cellular nuclease activity. One way to increase RNA’s resistance to nucleases is to replace its D-ribose backbone with L-ribose isomers. This modification changes chirality of an entire RNA molecule to L-form giving it more chance of survival when introduced into cells. Recently, we have described the activity of left-handed hammerhead ribozyme (L-Rz, L-HH) that can specifically hydrolyse RNA with the opposite chirality at a predetermined location. To understand the structural background of the RNA specific cleavage in a heterochiral complex, we used circular dichroism (CD) and nuclear magnetic resonance (NMR) spectroscopy as well as performed molecular modelling and dynamics simulations of homo- and heterochiral RNA complexes. The active ribozyme-target heterochiral complex showed a mixed chirality as well as low field imino proton NMR signals. We modelled the 3D structures of the oligoribonucleotides with their ribozyme counterparts of reciprocal chirality. L- or D-ribozyme formed a stable, homochiral helix 2, and two short double heterochiral helixes 1 and 3 of D- or L-RNA strand thorough irregular Watson–Crick base pairs. The formation of the heterochiral complexes is supported by the result of simulation molecular dynamics. These new observations suggest that L-catalytic nucleic acids can be used as tools in translational biology and diagnostics.

Structural basis of KdpD histidine kinase binding to the second messenger c-di-AMP

The KdpDE two-component system regulates potassium homeostasis and virulence in various bacterial species. The KdpD histidine kinases (HK) of this system contain a universal stress protein (USP) domain which binds to the second messenger cyclic-di-adenosine monophosphate (c-di-AMP) for regulating transcriptional output from this two-component system in Firmicutes such as Staphylococcus aureus. However, the structural basis of c-di-AMP specificity within the KdpD-USP domain is not well understood. Here, we resolved a 2.3 Å crystal structure of the S. aureus KdpD-USP domain (USPSa) complexed with c-di-AMP. Binding affinity analyses of USPSa mutants targeting the observed USPSa:c-di-AMP structural interface enabled the identification of the sequence residues that are required for c-di-AMP specificity. Based on the conservation of these residues in other Firmicutes, we identified the binding motif, (A/G/C)XSXSX2N(Y/F), which allowed us to predict c-di-AMP binding in other KdpD HKs. Furthermore, we found that the USPSa domain contains structural features distinct from the canonical standalone USPs that bind ATP as a preferred ligand. These features include inward-facing conformations of its β1-α1 and β4-α4 loops, a short α2 helix, the absence of a triphosphate-binding Walker A motif, and a unique dual phospho-ligand binding mode. It is therefore likely that USPSa-like domains in KdpD HKs represent a novel subfamily of the USPs.

Experimental and theoretical rationalization for the base pairing abilities of inosine, guanosine, adenosine, and their corresponding 8-oxo-7,8-dihydropurine, and 8-bromopurine analogues within A-form duplexes of RNA

Inosine is an important RNA modification, furthermore RNA oxidation has gained interest due, in part, to its potential role in the development/progression of disease as well as on its impact on RNA structure and function. In this report we established the base pairing abilities of purine nucleobases G, I, A, as well as their corresponding, 8-oxo-7,8-dihydropurine (common products of oxidation at the C8-position of purines), and 8-bromopurine (as probes to explore conformational changes), derivatives, namely 8-oxoG, 8-oxoI, 8-oxoA, 8-BrG, and 8-BrI. Dodecamers of RNA were obtained using standard phosphoramidite chemistry via solid-phase synthesis, and used as models to establish the impact that each of these nucleobases have on the thermal stability of duplexes, when base pairing to canonical and noncanonical nucleobases. Thermal stabilities were obtained from thermal denaturation transition (Tm ) measurements, via circular dichroism (CD). The results were then rationalized using models of base pairs between two monomers, via density functional theory (DFT), that allowed us to better understand potential contributions from H-bonding patterns arising from distinct conformations. Overall, some of the important results indicate that: (a) an anti-I:syn-A base pair provides thermal stability, due to the absence of the exocyclic amine; (b) 8-oxoG base pairs like U, and does not induce destabilization within the duplex when compared to the pyrimidine ring; (c) a U:G wobble-pair is only stabilized by G; and (d) 8-oxoA displays an inherited base pairing promiscuity in this sequence context. Gaining a better understanding of how this oxidatively generated lesions potentially base pair with other nucleobases will be useful to predict various biological outcomes, as well as in the design of biomaterials and/or nucleotide derivatives with biological potential.

Assessment of the nucleotide modifications in the high-resolution cryo-electron microscopy structure of the Escherichia coli 50S subunit

Post-transcriptional ribosomal RNA (rRNA) modifications are present in all organisms, but their exact functional roles and positions are yet to be fully characterized. Modified nucleotides have been implicated in the stabilization of RNA structure and regulation of ribosome biogenesis and protein synthesis. In some instances, rRNA modifications can confer antibiotic resistance. High-resolution ribosome structures are thus necessary for precise determination of modified nucleotides' positions, a task that has previously been accomplished by X-ray crystallography. Here, we present a cryo-electron microscopy (cryo-EM) structure of the Escherichia coli 50S subunit at an average resolution of 2.2 Å as an additional approach for mapping modification sites. Our structure confirms known modifications present in 23S rRNA and additionally allows for localization of Mg2+ ions and their coordinated water molecules. Using our cryo-EM structure as a testbed, we developed a program for assessment of cryo-EM map quality. This program can be easily used on any RNA-containing cryo-EM structure, and an associated Coot plugin allows for visualization of validated modifications, making it highly accessible.

Conformational ensembles of an RNA hairpin using molecular dynamics and sparse NMR data

Solution nuclear magnetic resonance (NMR) experiments allow RNA dynamics to be determined in an aqueous environment. However, when a limited number of peaks are assigned, it is difficult to obtain structural information. We here show a protocol based on the combination of experimental data (Nuclear Overhauser Effect, NOE) and molecular dynamics simulations with enhanced sampling methods. This protocol allows to (a) obtain a maximum entropy ensemble compatible with NMR restraints and (b) obtain a minimal set of metastable conformations compatible with the experimental data (maximum parsimony). The method is applied to a hairpin of 29 nt from an inverted SINEB2, which is part of the SINEUP family and has been shown to enhance protein translation. A clustering procedure is introduced where the annotation of base-base interactions and glycosidic bond angles is used as a metric. By reweighting the contributions of the clusters, minimal sets of four conformations could be found which are compatible with the experimental data. A motif search on the structural database showed that some identified low-population states are present in experimental structures of other RNA transcripts. The introduced method can be applied to characterize RNA dynamics in systems where a limited amount of NMR information is available.

Small Molecule Rescue and Glycosidic Conformational Analysis of the Twister Ribozyme

The number of self-cleaving ribozymes has increased sharply in recent years, giving rise to elaborations of the four known ribozyme catalytic strategies, α, β, γ, and δ. One such extension is utilized by the twister ribozyme, which is hypothesized to conduct δ, or general acid catalysis, via N3 of the syn adenine +1 nucleobase indirectly via buffer catalysis at biological pH and directly at lower pH. Herein, we test the δ catalysis role of A1 via chemical rescue and the catalytic relevance of the syn orientation of the nucleobase by conformational analysis. Using inhibited twister ribozyme variants with A1(N3) deaza or A1 abasic modifications, we observe >100-fold chemical rescue effects in the presence of protonatable biological small molecules such as imidazole and histidine, similar to observed rescue values previously reported for C75U/C76Δ in the HDV ribozyme. Brønsted plots for the twister variants support a model in which small molecules rescue catalytic activity via a proton transfer mechanism, suggesting that A1 in the wild type is involved in proton transfer, most likely general acid catalysis. Additionally, through glycosidic conformational analysis in an appropriate background that accommodates the bromine atom, we observe that an 8BrA1-modified twister ribozyme is up to 10-fold faster than a nonmodified A1 ribozyme, supporting crystallographic data that show that A1 is syn when conducting proton transfer. Overall, this study provides functional evidence that the nucleotide immediately downstream of the cleavage site participates directly or indirectly in general acid-base catalysis in the twister ribozyme while occupying the syn conformation.

Structural basis of amino acid surveillance by higher-order tRNA-mRNA interactions

Amino acid availability in Gram-positive bacteria is monitored by T-box riboswitches. T-boxes directly bind tRNAs, assess their aminoacylation state, and regulate the transcription or translation of downstream genes to maintain nutritional homeostasis. Here, we report cocrystal and cryo-EM structures of Geobacillus kaustophilus and Bacillus subtilis T-box-tRNA complexes, detailing their multivalent, exquisitely selective interactions. The T-box forms a U-shaped molecular vise that clamps the tRNA, captures its 3' end using an elaborate 'discriminator' structure, and interrogates its aminoacylation state using a steric filter fashioned from a wobble base pair. In the absence of aminoacylation, T-boxes clutch tRNAs and form a continuously stacked central spine, permitting transcriptional readthrough or translation initiation. A modeled aminoacyl disrupts tRNA-T-box stacking, severing the central spine and blocking gene expression. Our data establish a universal mechanism of amino acid sensing on tRNAs and gene regulation by T-box riboswitches and exemplify how higher-order RNA-RNA interactions achieve multivalency and specificity.

In silico survey of the central conserved regions in viroids of the Pospiviroidae family for conserved asymmetric loop structures

Viroids are the smallest replicative pathogens, consisting of RNA circles (∼300 nucleotides) that require host machinery to replicate. Structural RNA elements recruit these host factors. Currently, many of these structural elements and the nature of their interactions are unknown. All Pospiviroidae have homology in the central conserved region (CCR). The CCR of potato spindle tuber viroid (PSTVd) contains a sarcin/ricin domain (SRD), the only viroid structural element with an unequivocal replication role. We assumed that every member of this family uses this region to recruit host factors, and that each CCR has an SRD-like asymmetric loop within it. Potential SRD or SRD-like motifs were sought in the CCR of each Pospiviroidae member as follows. Motif location in each CCR was predicted with MUSCLE alignment and Vienna RNAfold. Viroid-specific models of SRD-like motifs were built by superimposing noncanonical base pairs and nucleotides on a model of an SRD. The RNA geometry search engine FR3D was then used to find nucleotide groups close to the geometry suggested by this superimposition. Atomic resolution structures were assembled using the molecular visualization program Chimera, and the stability of each motif was assessed with molecular dynamics (MD). Some models required a protonated cytosine. To be stable within a cell, the pK_a of that cytosine must be shifted up. Constant pH-replica exchange MD analysis showed such a shift in the proposed structures. These data show that every Pospiviroidae member could form a motif that resembles an SRD in its CCR, and imply there could be undiscovered mimics of other RNA domains.

Uncovering a unique approach for damaged DNA replication: A computational investigation of a mutagenic tobacco-derived thymine lesion

4-(Methylnitrosamino)-1-(3-pyridyl)-1-butanone is a potent nicotine carcinogen that leads to many DNA lesions, the most persistent being the O2-[4-oxo-4-(3-pyridyl)butyl]thymine adduct (POB-T). Although the experimental mutagenic profile for the minor groove POB-T lesion has been previously reported, the findings are puzzling in terms of the human polymerases involved. Specifically, while pol κ typically replicates minor groove adducts, in vivo studies indicate pol η replicates POB-T despite being known for processing major groove adducts. Our multiscale modeling approach reveals that the canonical (anti) glycosidic orientation of POB-T can fit in the pol κ active site, but only a unique (syn) POB-T conformation is accommodated by pol η. These distinct binding orientations rationalize the differential in vitro mutagenic spectra based on the preferential stabilization of dGTP and dTTP opposite the lesion for pol κ and η, respectively. Overall, by uncovering the first evidence for the replication of a damaged pyrimidine in the syn glycosidic orientation, the current work provides the insight necessary to clarify a discrepancy in the DNA replication literature, expand the biological role of the critical human pol η, and understand the mutational signature in human cancers associated with tobacco exposure.

Assembly and cryo-EM structures of RNA-specific measles virus nucleocapsids provide mechanistic insight into paramyxoviral replication

Assembly of paramyxoviral nucleocapsids on the RNA genome is an essential step in the viral cycle. The structural basis of this process has remained obscure due to the inability to control encapsidation. We used a recently developed approach to assemble measles virus nucleocapsid-like particles on specific sequences of RNA hexamers (poly-Adenine and viral genomic 5') in vitro, and determined their cryoelectron microscopy maps to 3.3-Å resolution. The structures unambiguously determine 5' and 3' binding sites and thereby the binding-register of viral genomic RNA within nucleocapsids. This observation reveals that the 3' end of the genome is largely exposed in fully assembled measles nucleocapsids. In particular, the final three nucleotides of the genome are rendered accessible to the RNA-dependent RNA polymerase complex, possibly enabling efficient RNA processing. The structures also reveal local and global conformational changes in the nucleoprotein upon assembly, in particular involving helix α6 and helix α13 that form edges of the RNA binding groove. Disorder is observed in the bound RNA, localized at one of the two backbone conformational switch sites. The high-resolution structure allowed us to identify putative nucleobase interaction sites in the RNA-binding groove, whose impact on assembly kinetics was measured using real-time NMR. Mutation of one of these sites, R195, whose sidechain stabilizes both backbone and base of a bound nucleic acid, is thereby shown to be essential for nucleocapsid-like particle assembly.

Synthesis of 1',2'-methano-2',3'-dideoxynucleosides as potential antivirals

The synthesis of constrained nucleosides has become an important tool to understand the SAR in the interaction between biological and synthetic nucleotides in the context of antisense oligonucleotide therapy. The incorporation of a cyclopropane into a furanose ring of a nucleoside induces some degree of constrain without affecting significantly the steric environment of a nucleoside. Here, we report a new, short and stereocontrolled synthesis of two constrained nucleosides analogues, 1',2'- methano-2',3'-dideoxyuridine 9, and the corresponding cytidine analog 12. X-ray crystallography revealed that the furanose ring in the constrained uridine and cytidine analogues was flattened with virtual loss of pseudorotation. The phosphoramidate esters of the novel constrained uridine and cytidine nucleosides, intended as prodrugs, were tested in cell-based assays for viral replication across the herpes virus family and HIV inhibition courtesy of Merck laboratories, Rahway. They were also tested in antiproliferative assays against colorectal and melanoma cell lines. Unfortunately, none of the compounds showed activity in these assays.

ykkC riboswitches employ an add-on helix to adjust specificity for polyanionic ligands

The ykkC family of bacterial riboswitches combines several widespread classes that have similar secondary structures and consensus motifs but control different genes in response to different cellular metabolites. Here we report the crystal structures of two distinct ykkC riboswitches specifically bound to their cognate ligands ppGpp, a second messenger involved in stress response, or PRPP, a precursor in purine biosynthesis. Both RNAs adopt similar structures and contain a conserved core previously observed in the guanidine-specific ykkC riboswitch. However, ppGpp and PRPP riboswitches uniquely employ an additional helical element that joins the ends of the ligand-sensing domains and creates a tunnel for direct and Mg²⁺-mediated binding of ligands. Mutational and footprinting experiments highlight the importance of conserved nucleotides forming the tunnel and long-distance contacts for ligand binding and genetic response. Our work provides new insights into the specificity of riboswitches and gives a unique opportunity for future studies of RNA evolution.

Crystal structures of a group II intron lariat primed for reverse splicing

The 2'-5' branch of nuclear premessenger introns is believed to have been inherited from self-splicing group II introns, which are retrotransposons of bacterial origin. Our crystal structures at 3.4 and 3.5 angstrom of an excised group II intron in branched ("lariat") form show that the 2'-5' branch organizes a network of active-site tertiary interactions that position the intron terminal 3'-hydroxyl group into a configuration poised to initiate reverse splicing, the first step in retrotransposition. Moreover, the branchpoint and flanking helices must undergo a base-pairing switch after branch formation. A group II-based model of the active site of the nuclear splicing machinery (the spliceosome) is proposed. The crucial role of the lariat conformation in active-site assembly and catalysis explains its prevalence in modern splicing.

Understanding in-line probing experiments by modeling cleavage of nonreactive RNA nucleotides

Ribonucleic acid (RNA) is involved in many regulatory and catalytic processes in the cell. The function of any RNA molecule is intimately related with its structure. In-line probing experiments provide valuable structural data sets for a variety of RNAs and are used to characterize conformational changes in riboswitches. However, the structural determinants that lead to differential reactivities in unpaired nucleotides have not been investigated yet. In this work, we used a combination of theoretical approaches, i.e., classical molecular dynamics simulations, multiscale quantum mechanical/molecular mechanical calculations, and enhanced sampling techniques in order to compute and interpret the differential reactivity of individual residues in several RNA motifs, including members of the most important GNRA and UNCG tetraloop families. Simulations on the multinanosecond timescale are required to converge the related free-energy landscapes. The results for uGAAAg and cUUCGg tetraloops and double helices are compared with available data from in-line probing experiments and show that the introduced technique is able to distinguish between nucleotides of the uGAAAg tetraloop based on their structural predispositions toward phosphodiester backbone cleavage. For the cUUCGg tetraloop, more advanced ab initio calculations would be required. This study is the first attempt to computationally classify chemical probing experiments and paves the way for an identification of tertiary structures based on the measured reactivity of nonreactive nucleotides.

Predicting the Kinetics of RNA Oligonucleotides Using Markov State Models

Nowadays different experimental techniques, such as single molecule or relaxation experiments, can provide dynamic properties of biomolecular systems, but the amount of detail obtainable with these methods is often limited in terms of time or spatial resolution. Here we use state-of-the-art computational techniques, namely, atomistic molecular dynamics and Markov state models, to provide insight into the rapid dynamics of short RNA oligonucleotides, to elucidate the kinetics of stacking interactions. Analysis of multiple microsecond-long simulations indicates that the main relaxation modes of such molecules can consist of transitions between alternative folded states, rather than between random coils and native structures. After properly removing structures that are artificially stabilized by known inaccuracies of the current RNA AMBER force field, the kinetic properties predicted are consistent with the time scales of previously reported relaxation experiments.

Purine 3':5'-cyclic nucleotides with the nucleobase in a syn orientation: cAMP, cGMP and cIMP

Purine 3':5'-cyclic nucleotides are very well known for their role as the secondary messengers in hormone action and cellular signal transduction. Nonetheless, their solid-state conformational details still require investigation. Five crystals containing purine 3':5'-cyclic nucleotides have been obtained and structurally characterized, namely adenosine 3':5'-cyclic phosphate dihydrate, C10H12N5O6P·2H2O or cAMP·2H2O, (I), adenosine 3':5'-cyclic phosphate 0.3-hydrate, C10H12N5O6P·0.3H2O or cAMP·0.3H2O, (II), guanosine 3':5'-cyclic phosphate pentahydrate, C10H12N5O7P·5H2O or cGMP·5H2O, (III), sodium guanosine 3':5'-cyclic phosphate tetrahydrate, Na(+)·C10H11N5O7P(-)·4H2O or Na(cGMP)·4H2O, (IV), and sodium inosine 3':5'-cyclic phosphate tetrahydrate, Na(+)·C10H10N4O7P(-)·4H2O or Na(cIMP)·4H2O, (V). Most of the cyclic nucleotide zwitterions/anions [two from four cAMP present in total in (I) and (II), cGMP in (III), cGMP(-) in (IV) and cIMP(-) in (V)] are syn conformers about the N-glycosidic bond, and this nucleobase arrangement is accompanied by Crib-H...Npur hydrogen bonds (rib = ribose and pur = purine). The base orientation is tuned by the ribose pucker. An analysis of data obtained from the Cambridge Structural Database made in the context of syn-anti conformational preferences has revealed that among the syn conformers of various purine nucleotides, cyclic nucleotides and dinucleotides predominate significantly. The interactions stabilizing the syn conformation have been indicated. The inter-nucleotide contacts in (I)-(V) have been systematized in terms of the chemical groups involved. All five structures display three-dimensional hydrogen-bonded networks.

'Z-DNA like' fragments in RNA: a recurring structural motif with implications for folding, RNA/protein recognition and immune response

Since the work of Alexander Rich, who solved the first Z-DNA crystal structure, we have known that d(CpG) steps can adopt a particular structure that leads to forming left-handed helices. However, it is still largely unrecognized that other sequences can adopt ‘left-handed’ conformations in DNA and RNA, in double as well as single stranded contexts. These ‘Z-like’ steps involve the coexistence of several rare structural features: a C2’-endo puckering, a syn nucleotide and a lone pair–π stacking between a ribose O4’ atom and a nucleobase. This particular arrangement induces a conformational stress in the RNA backbone, which limits the occurrence of Z-like steps to ≈0.1% of all dinucleotide steps in the PDB. Here, we report over 600 instances of Z-like steps, which are located within r(UNCG) tetraloops but also in small and large RNAs including riboswitches, ribozymes and ribosomes. Given their complexity, Z-like steps are probably associated with slow folding kinetics and once formed could lock a fold through the formation of unique long-range contacts. Proteins involved in immunologic response also specifically recognize/induce these peculiar folds. Thus, characterizing the conformational features of these motifs could be a key to understanding the immune response at a structural level.

The role of RNA conformation in RNA-protein recognition

Interactions between protein and RNA play a key role in many biological processes in the gene expression pathway. Those interactions are mediated through a variety of RNA-binding protein domains, among them the highly abundant RNA recognition motif (RRM). Here we studied protein-RNA complexes from different RNA binding domain families solved by NMR and x-ray crystallography. Characterizing the structural properties of the RNA at the binding interfaces revealed an unexpected number of nucleotides with unusual RNA conformations, specifically found in RNA-RRM complexes. Moreover, we observed that the RNA nucleotides that are directly involved in interactions with the RRM domains, via hydrogen bonds and hydrophobic contacts, are significantly enriched with unique RNA conformations. Further examination of the sequences binding the RRM domain showed a preference for G nucleotides in syn conformation to precede or to follow U nucleotides in the anti-conformation, and U nucleotides in C2' endo conformation to precede U and G nucleotides possessing the more common C3' endo conformation. These findings imply a possible mode of RNA recognition by the RRM domains which enables the recognition of a wide variety of different RNA sequences and shapes. Overall, this study suggests an additional way by which the RRM domain recognizes its RNA target, involving a conformational readout.

Structural Flexibility and Conformation Features of Cyclic Dinucleotides in Aqueous Solutions

Cyclic dinucleotides are able to trigger the innate immune system by activating STING. It was found that the binding affinity of asymmetric 2'3'-cGAMP to symmetric dimer of STING is 3 orders of magnitude higher than that of the symmetric 3'3'-cyclic dinucleotides. Such a phenomenon has not been understood yet. Here we show that the subtle changes in phosphodiester linkage of CDNs lead to their distinct structural properties which correspond to the varied binding affinities. 2'-5' and/or 3'-5' linked CDNs adopt specific while different types of ribose puckers and backbone conformations. That ribose conformations and base types have different propensities for anti or syn glycosidic conformations further affects the overall flexibility of CDNs. The counterbalance between backbone ring tension and electrostatic repulsion, both affected by the ring size, also contributes to the different flexibility of CDNs. Our calculations reveal that the free energy cost for 2'3'-cGAMP to adopt the STING-bound structure is smaller than that for 3'3'-cGAMP and cyclic-di-GMP. These findings may serve as a reference for design of CDN-analogues as vaccine adjuvants. Moreover, the cyclization pattern of CDNs closely related to their physiological roles suggests the importance of understanding structural properties in the study of protein-ligand interactions.

Embracing proteins: structural themes in aptamer-protein complexes

Understanding the structural rules that govern specific, high-affinity binding characteristic of aptamer-protein interactions is important in view of the increasing use of aptamers across many applications. From the modest number of 16 aptamer-protein structures currently available, trends are emerging. The flexible phosphodiester backbone allows folding into precise three-dimensional structures using known nucleic acid motifs as scaffolds that orient specific functional groups for target recognition. Still, completely novel motifs essential for structure and function are found in modified aptamers with diversity-enhancing side chains. Aptamers and antibodies, two classes of macromolecules used as affinity reagents with entirely different backbones and composition, recognize protein epitopes of similar size and with comparably high shape complementarity.

A Two-Metal-Ion-Mediated Conformational Switching Pathway for HDV Ribozyme Activation

RNA enzymes serve as a potentially powerful platform from which to design catalysts and engineer new biotechnology. A fundamental understanding of these systems provides insight to guide design. The hepatitis delta virus ribozyme (HDVr) is a small, self-cleaving RNA motif widely distributed in nature, that has served as a paradigm for understanding basic principles of RNA catalysis. Nevertheless, questions remain regarding the precise roles of divalent metal ions and key nucleotides in catalysis. In an effort to establish a reaction mechanism model consistent with available experimental data, we utilize molecular dynamics simulations to explore different conformations and metal ion binding modes along the HDVr reaction path. Building upon recent crystallographic data, our results provide a dynamic model of the HDVr reaction mechanism involving a conformational switch between multiple non-canonical G25:U20 base pair conformations in the active site. These local nucleobase dynamics play an important role in catalysis by modulating the metal binding environments of two Mg²⁺ ions that support catalysis at different steps of the reaction pathway. The first ion plays a structural role by inducing a base pair flip necessary to obtain the catalytic fold in which C75 moves towards to the scissile phosphate in the active site. Ejection of this ion then permits a second ion to bind elsewhere in the active site and facilitate nucleophile activation. The simulations collectively describe a mechanistic scenario that is consistent with currently available experimental data from crystallography, phosphorothioate substitutions, and chemical probing studies. Avenues for further experimental verification are suggested.

Stacking in RNA: NMR of Four Tetramers Benchmark Molecular Dynamics

Molecular dynamics (MD) simulations for RNA tetramers r(AAAA), r(CAAU), r(GACC), and r(UUUU) are benchmarked against 1H-1H NOESY distances and 3J scalar couplings to test effects of RNA torsion parametrizations. Four different starting structures were used for r(AAAA), r(CAAU), and r(GACC), while five starting structures were used for r(UUUU). On the basis of X-ray structures, criteria are reported for quantifying stacking. The force fields, AMBER ff99, parmbsc0, parm99χ_Yil, ff10, and parmTor, all predict experimentally unobserved stacks and intercalations, e.g., base 1 stacked between bases 3 and 4, and incorrect χ, ϵ, and sugar pucker populations. The intercalated structures are particularly stable, often lasting several microseconds. Parmbsc0, parm99χ_Yil, and ff10 give similar agreement with NMR, but the best agreement is only 46%. Experimentally unobserved intercalations typically are associated with reduced solvent accessible surface area along with amino and hydroxyl hydrogen bonds to phosphate nonbridging oxygens. Results from an extensive set of MD simulations suggest that recent force field parametrizations improve predictions, but further improvements are necessary to provide reasonable agreement with NMR. In particular, intramolecular stacking and hydrogen bonding interactions may not be well balanced with the TIP3P water model. NMR data and the scoring method presented here provide rigorous benchmarks for future changes in force fields and MD methods.

High-resolution structure of the Escherichia coli ribosome

Protein synthesis by the ribosome is highly dependent on the ionic conditions in the cellular environment, but the roles of ribosome solvation remain poorly understood. Moreover, the function of modifications to ribosomal RNA and ribosomal proteins are unclear. Here we present the structure of the Escherichia coli 70S ribosome to 2.4 Å resolution. The structure reveals details of the ribosomal subunit interface that are conserved in all domains of life, and suggest how solvation contributes to ribosome integrity and function. The structure also suggests how the conformation of ribosomal protein uS12 likely impacts its contribution to messenger RNA decoding. This structure helps to explain the phylogenetic conservation of key elements of the ribosome, including posttranscriptional and posttranslational modifications and should serve as a basis for future antibiotic development.

Structural basis for the fast self-cleavage reaction catalyzed by the twister ribozyme

Twister is a recently discovered RNA motif that is estimated to have one of the fastest known catalytic rates of any naturally occurring small self-cleaving ribozyme. We determined the 4.1-Å resolution crystal structure of a twister sequence from an organism that has not been cultured in isolation, and it shows an ordered scissile phosphate and nucleotide 5' to the cleavage site. A second crystal structure of twister from Orzyza sativa determined at 3.1-Å resolution exhibits a disordered scissile phosphate and nucleotide 5' to the cleavage site. The core of twister is stabilized by base pairing, a large network of stacking interactions, and two pseudoknots. We observe three nucleotides that appear to mediate catalysis: a guanosine that we propose deprotonates the 2'-hydroxyl of the nucleotide 5' to the cleavage site and a conserved adenosine. We suggest the adenosine neutralizes the negative charge on a nonbridging phosphate oxygen atom at the cleavage site. The active site also positions the labile linkage for in-line nucleophilic attack, and thus twister appears to simultaneously use three strategies proposed for small self-cleaving ribozymes. The twister crystal structures (i) show its global structure, (ii) demonstrate the significance of the double pseudoknot fold, (iii) provide a possible hypothesis for enhanced catalysis, and (iv) illuminate the roles of all 10 highly conserved nucleotides of twister that participate in the formation of its small and stable catalytic pocket.

Optimization of an AMBER force field for the artificial nucleic acid, LNA, and benchmarking with NMR of L(CAAU)

Locked Nucleic Acids (LNAs) are RNA analogues with an O2'-C4' methylene bridge which locks the sugar into a C3'-endo conformation. This enhances hybridization to DNA and RNA, making LNAs useful in microarrays and potential therapeutics. Here, the LNA, L(CAAU), provides a simplified benchmark for testing the ability of molecular dynamics (MD) to approximate nucleic acid properties. LNA χ torsions and partial charges were parametrized to create AMBER parm99_LNA. The revisions were tested by comparing MD predictions with AMBER parm99 and parm99_LNA against a 200 ms NOESY NMR spectrum of L(CAAU). NMR indicates an A-Form equilibrium ensemble. In 3000 ns simulations starting with an A-form structure, parm99_LNA and parm99 provide 66% and 35% agreement, respectively, with NMR NOE volumes and (3)J-couplings. In simulations of L(CAAU) starting with all χ torsions in a syn conformation, only parm99_LNA is able to repair the structure. This implies methods for parametrizing force fields for nucleic acid mimics can reasonably approximate key interactions and that parm99_LNA will improve reliability of MD studies for systems with LNA. A method for approximating χ population distribution on the basis of base to sugar NOEs is also introduced.

Ribonucleotides as nucleotide excision repair substrates

The incorporation of ribonucleotides in DNA has attracted considerable notice in recent years, since the pool of ribonucleotides can exceed that of the deoxyribonucleotides by at least 10-20-fold, and single ribonucleotide incorporation by DNA polymerases appears to be a common event. Moreover ribonucleotides are potentially mutagenic and lead to genome instability. As a consequence, errantly incorporated ribonucleotides are rapidly repaired in a process dependent upon RNase H enzymes. On the other hand, global genomic nucleotide excision repair (NER) in prokaryotes and eukaryotes removes damage caused by covalent modifications that typically distort and destabilize DNA through the production of lesions derived from bulky chemical carcinogens, such as polycyclic aromatic hydrocarbon metabolites, or via crosslinking. However, a recent study challenges this lesion-recognition paradigm. The work of Vaisman et al. (2013) [34] reveals that even a single ribonucleotide embedded in a deoxyribonucleotide duplex is recognized by the bacterial NER machinery in vitro. In their report, the authors show that spontaneous mutagenesis promoted by a steric-gate pol V mutant increases in uvrA, uvrB, or uvrC strains lacking rnhB (encoding RNase HII) and to a greater extent in an NER-deficient strain lacking both RNase HI and RNase HII. Using purified UvrA, UvrB, and UvrC proteins in in vitro assays they show that despite causing little distortion, a single ribonucleotide embedded in a DNA duplex is recognized and doubly-incised by the NER complex. We present the hypothesis to explain the recognition and/or verification of this small lesion, that the critical 2'-OH of the ribonucleotide - with its unique electrostatic and hydrogen bonding properties - may act as a signal through interactions with amino acid residues of the prokaryotic NER complex that are not possible with DNA. Such a mechanism might also be relevant if it were demonstrated that the eukaryotic NER machinery likewise incises an embedded ribonucleotide in DNA.

Conformational readout of RNA by small ligands

RNA molecules have highly versatile structures that can fold into myriad conformations, providing many potential pockets for binding small molecules. The increasing number of available RNA structures, in complex with proteins, small ligands and in free form, enables the design of new therapeutically useful RNA-binding ligands. Here we studied RNA ligand complexes from 10 RNA groups extracted from the protein data bank (PDB), including adaptive and non-adaptive complexes. We analyzed the chemical, physical, structural and conformational properties of binding pockets around the ligand. Comparing the properties of ligand-binding pockets to the properties of computed pockets extracted from all available RNA structures and RNA-protein interfaces, revealed that ligand-binding pockets, mainly the adaptive pockets, are characterized by unique properties, specifically enriched in rare conformations of the nucleobase and the sugar pucker. Further, we demonstrate that nucleotides possessing the rare conformations are preferentially involved in direct interactions with the ligand. Overall, based on our comprehensive analysis of RNA-ligand complexes, we suggest that the unique conformations adopted by RNA nucleotides play an important role in RNA recognition by small ligands. We term the recognition of a binding site by a ligand via the unique RNA conformations "RNA conformational readout." We propose that "conformational readout" is a general way by which RNA binding pockets are recognized and selected from an ensemble of different RNA states.

Fingerprinting noncanonical and tertiary RNA structures by differential SHAPE reactivity

Many RNA structures are composed of simple secondary structure elements linked by a few critical tertiary interactions. SHAPE chemistry has made interrogation of RNA dynamics at single-nucleotide resolution straightforward. However, de novo identification of nucleotides involved in tertiary interactions remains a challenge. Here we show that nucleotides that form noncanonical or tertiary contacts can be detected by comparing information obtained using two SHAPE reagents, N-methylisatoic anhydride (NMIA) and 1-methyl-6-nitroisatoic anhydride (1M6). Nucleotides that react preferentially with NMIA exhibit slow local nucleotide dynamics and usually adopt the less common C2'-endo ribose conformation. Experiments and first-principles calculations show that 1M6 reacts preferentially with nucleotides in which one face of the nucleobase allows an unhindered stacking interaction with the reagent. Differential SHAPE reactivities were used to detect noncanonical and tertiary interactions in four RNAs with diverse structures and to identify preformed noncanonical interactions in partially folded RNAs. Differential SHAPE reactivity analysis will enable experimentally concise, large-scale identification of tertiary structure elements and ligand binding sites in complex RNAs and in diverse biological environments.

Thermodynamics of ligand binding to a heterogeneous RNA population in the malachite green aptamer

The malachite green aptamer binds two closely related ligands, malachite green (MG) and tetramethylrosamine (TMR), with nearly equal affinity. The MG ligand consists of three phenyl rings emanating from a central carbon, while TMR has two of the three rings connected by an ether linkage. The binding pockets for MG and TMR in the aptamer, known from high-resolution structures, differ only in the conformation of a few nucleotides. Herein, we applied isothermal titration calorimetry (ITC) to compare the thermodynamics of binding of MG and TMR to the aptamer. Binding heat capacities were obtained from ITC titrations over the temperature range of 15-60 °C. Two temperature regimes were found for MG binding: one from 15 to 45 °C where MG bound with a large negative heat capacity and an apparent stoichiometry (n) of ~0.4 and another from 50 to 60 °C where MG bound with a positive heat capacity and an n of ~1.1. The binding of TMR, on the other hand, revealed only one temperature regime for binding, with a more modest negative heat capacity and an n of ~1.2. The large difference in heat capacity between the two ligands suggests that significantly more conformational rearrangement occurs upon the binding of MG than that of TMR, which is consistent with differences in solvent accessible surface area calculated for available ligand-bound structures. Lastly, we note that the binding stoichiometry of MG was improved not only by an increase in the temperature but also by a decrease in the concentration of Mg(2+) or an increase in the time between ITC injections. These studies suggest that binding of a dynamical ligand to a functional RNA requires the RNA itself to have significant dynamics.