Insights into the structural stability of major groove RNA triplexes by WAXS-guided MD simulations
Explore the source record for details and available documents.
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
The ribose 2'-hydroxyl is the key chemical difference between RNA and DNA and primary source of their divergent structural and functional characteristics. Macromolecular X-ray diffraction experiments typically do not reveal the positions of hydrogen atoms. Thus, standard crystallography cannot determine 2'-OH orientation (H2'-C2'-O2'-HO2' torsion angle) and its potential roles in sculpting the RNA backbone and the expansive fold space. Here, we report the first neutron crystal structure of an RNA, the Escherichia coli rRNA Sarcin-Ricin Loop (SRL). 2'-OD orientations were established for all 27 residues and revealed O-D bonds pointing toward backbone (O3', 13 observations), nucleobase (11) or sugar (3). Most riboses in the SRL stem region show a 2'-OD backbone-orientation. GAGA-tetraloop riboses display a 2'-OD base-orientation. An atypical C2'-endo sugar pucker is strictly correlated with a 2'-OD sugar-orientation. Neutrons reveal the strong preference of the 2'-OH to donate in H-bonds and that 2'-OH orientation affects both backbone geometry and ribose pucker. We discuss 2'-OH and water molecule orientations in the SRL neutron structure and compare with results from a solution phase 10 μs MD simulation. We demonstrate that joint cryo-neutron/X-ray crystallography offers an all-in-one approach to determine the complete structural properties of RNA, i.e. geometry, conformation, protonation state and hydration structure.
Coronavirus genomes sequester their start codons within stem-loop 5 (SL5), a structured, 5' genomic RNA element. In most alpha- and betacoronaviruses, the secondary structure of SL5 is predicted to contain a four-way junction of helical stems, some of which are capped with UUYYGU hexaloops. Here, using cryogenic electron microscopy (cryo-EM) and computational modeling with biochemically determined secondary structures, we present three-dimensional structures of SL5 from six coronaviruses. The SL5 domain of betacoronavirus severe-acute-respiratory-syndrome-related coronavirus 2 (SARS-CoV-2), resolved at 4.7 Å resolution, exhibits a T-shaped structure, with its UUYYGU hexaloops at opposing ends of a coaxial stack, the T’s “arms.” Further analysis of SL5 domains from SARS-CoV-1 and MERS (7.1 and 6.4 to 6.9 Å resolution, respectively) indicate that the junction geometry and inter-hexaloop distances are conserved features across these human-infecting betacoronaviruses. The MERS SL5 domain displays an additional tertiary interaction, which is also observed in the non-human-infecting betacoronavirus BtCoV-HKU5 (5.9 to 8.0 Å resolution). SL5s from human-infecting alphacoronaviruses, HCoV-229E and HCoV-NL63 (6.5 and 8.4 to 9.0 Å resolution, respectively), exhibit the same coaxial stacks, including the UUYYGU-capped arms, but with a phylogenetically distinct crossing angle, an X-shape. As such, all SL5 domains studied herein fold into stable tertiary structures with cross-genus similarities and notable differences, with implications for potential protein-binding modes and therapeutic targets.
La-related protein 7 (LARP7) are a family of RNA chaperones that protect the 3'-end of RNA and are components of specific ribonucleoprotein complexes (RNP). Here, in Tetrahymena thermophila telomerase, LARP7 protein p65 together with telomerase reverse transcriptase (TERT) and telomerase RNA (TER) form the core RNP. p65 has four known domains—N-terminal domain (NTD), La motif (LaM), RNA recognition motif 1 (RRM1), and C-terminal xRRM2. To date, only the xRRM2 and LaM and their interactions with TER have been structurally characterized. Conformational dynamics leading to low resolution in cryo-EM density maps have limited our understanding of how full-length p65 specifically recognizes and remodels TER for telomerase assembly. Here, we combined focused classification of Tetrahymena telomerase cryo-EM maps with NMR spectroscopy to determine the structure of p65–TER. Three previously unknown helices are identified, one in the otherwise intrinsically disordered NTD that binds the La module, one that extends RRM1, and another preceding xRRM2, that stabilize p65–TER interactions. The extended La module (αN, LaM and RRM1) interacts with the four 3' terminal U nucleotides, while LaM and αN additionally interact with TER pseudoknot, and LaM with stem 1 and 5' end. Our results reveal the extensive p65–TER interactions that promote TER 3'-end protection, TER folding, and core RNP assembly and stabilization. The structure of full-length p65 with TER also sheds light on the biological roles of genuine La and LARP7 proteins as RNA chaperones and core RNP components.
Structure-based virtual high-throughput screening is used in early-stage drug discovery. Over the years, docking protocols and scoring functions for protein–ligand complexes have evolved to improve the accuracy in the computation of binding strengths and poses. In the past decade, RNA has also emerged as a target class for new small-molecule drugs. However, most ligand docking programs have been validated and tested for proteins and not RNA. Here, we test the docking power (pose prediction accuracy) of three state-of-the-art docking protocols on 173 RNA–small molecule crystal structures. The programs are AutoDock4 (AD4) and AutoDock Vina (Vina), which were designed for protein targets, and rDock, which was designed for both protein and nucleic acid targets. AD4 performed relatively poorly. For RNA targets for which a crystal structure of a bound ligand used to limit the docking search space is available and for which the goal is to identify new molecules for the same pocket, rDock performs slightly better than Vina, with success rates of 48% and 63%, respectively. However, in the more common type of early-stage drug discovery setting, in which no structure of a ligand–target complex is known and for which a larger search space is defined, rDock performed similarly to Vina, with a low success rate of ~27%. Further, Vina was found to have bias for ligands with certain physicochemical properties, whereas rDock performs similarly for all ligand properties. Thus, for projects where no ligand–protein structure already exists, Vina and rDock are both applicable. However, the relatively poor performance of all methods relative to protein–target docking illustrates a need for further methods refinement.
Single-particle cryogenic electron microscopy (cryo-EM) has become a standard technique for determining protein structures at atomic resolution1,2,3. However, cryo-EM studies of protein-free RNA are in their early days. The Tetrahymena thermophila group I self-splicing intron was the first ribozyme to be discovered and has been a prominent model system for the study of RNA catalysis and structure–function relationships4, but its full structure remains unknown. Here we report cryo-EM structures of the full-length Tetrahymena ribozyme in substrate-free and bound states at a resolution of 3.1 Å. Additionally, newly resolved peripheral regions form two coaxially stacked helices; these are interconnected by two kissing loop pseudoknots that wrap around the catalytic core and include two previously unforeseen (to our knowledge) tertiary interactions. The global architecture is nearly identical in both states; only the internal guide sequence and guanosine binding site undergo a large conformational change and a localized shift, respectively, upon binding of RNA substrates. These results provide a long-sought structural view of a paradigmatic RNA enzyme and signal a new era for the cryo-EM-based study of structure–function relationships in ribozymes.
Fungal Trl1 is an essential tRNA splicing enzyme composed of C-terminal cyclic phosphodiesterase and central polynucleotide kinase end-healing domains that convert the 2′,3′-cyclic-PO 4 and 5′-OH ends of tRNA exons into the 3′-OH,2′-PO 4 and 5′-PO 4 termini required for sealing by an N-terminal ATP-dependent ligase domain. Trifunctional Trl1 enzymes are present in most human fungal pathogens and are untapped targets for antifungal drug discovery. Mucorales species, deemed high-priority human pathogens by WHO, elaborate a noncanonical tRNA splicing apparatus in which a stand-alone monofunctional RNA ligase enzyme joins 3′-OH,2′-PO 4 and 5′-PO 4 termini. Here we identify a stand-aloneMucor circinelloidespolynucleotide kinase (MciKIN) and affirm its biological activity in tRNA splicing by genetic complementation in yeast. Recombinant MciKIN catalyzes magnesium-dependent phosphorylation of 5′-OH RNA and DNA ends in vitro. MciKIN displays a strong preference for GTP as the phosphate donor in the kinase reaction, a trait shared with the stand-alone RNA kinase homologs from Mucorales speciesRhizopus azygosporus(RazKIN) andLichtheimia corymbifera(LcoKIN) and with the kinase domains of fungal Trl1 enzymes. We report a 1.65 Å crystal structure of RazKIN in complex with GDP•Mg 2+ that illuminates the basis for guanosine nucleotide specificity.
Elongating RNA polymerase II (Pol II) can be paused or arrested by a variety of obstacles. These obstacles include DNA lesions, DNA-binding proteins, and small molecules. Hairpin pyrrole-imidazole (Py-Im) polyamides bind to the minor groove of DNA in a sequence-specific manner and induce strong transcriptional arrest. Remarkably, this Py-Im–induced Pol II transcriptional arrest is persistent and cannot be rescued by transcription factor TFIIS. In contrast, TFIIS can effectively rescue the transcriptional arrest induced by a nucleosome barrier. The structural basis of Py-Im–induced transcriptional arrest and why TFIIS cannot rescue this arrest remain elusive. Here we determined the X-ray crystal structures of four distinct Pol II elongation complexes (Pol II ECs) in complex with hairpin Py-Im polyamides as well as of the hairpin Py-Im polyamides–dsDNA complex. We observed that the Py-Im oligomer directly interacts with RNA Pol II residues, introduces compression of the downstream DNA duplex, prevents Pol II forward translocation, and induces Pol II backtracking. These results, together with biochemical studies, provide structural insight into the molecular mechanism by which Py-Im blocks transcription. Our structural study reveals why TFIIS fails to promote Pol II bypass of Py-Im–induced transcriptional arrest.
Abstract RNA modifications shape gene expression through a variety of chemical changes to canonical RNA bases. Although numbering in the hundreds, only a few RNA modifications are well characterized, in part due to the absence of methods to identify modification sites. Antibodies remain a common tool to identify modified RNA and infer modification sites through straightforward applications. However, specificity issues can result in off-target binding and confound conclusions. This work utilizes in silico λ-dynamics to efficiently estimate binding free energy differences of modification-targeting antibodies between a variety of naturally occurring RNA modifications. Crystal structures of inosine and N6-methyladenosine (m6A) targeting antibodies bound to their modified ribonucleosides were determined and served as structural starting points. λ-Dynamics was utilized to predict RNA modifications that permit or inhibit binding to these antibodies. In vitro RNA-antibody binding assays supported the accuracy of these in silico results. High agreement between experimental and computed binding propensities demonstrated that λ-dynamics can serve as a predictive screen for antibody specificity against libraries of RNA modifications. More importantly, this strategy is an innovative way to elucidate how hundreds of known RNA modifications interact with biological molecules without the limitations imposed by in vitro or in vivo methodologies.
Abstract The challenge of targeting RNA with small molecules necessitates a better understanding of RNA–ligand interaction mechanisms. However, the dynamic nature of nucleic acids, their ligand‐induced stabilization, and how conformational changes influence gene expression pose significant difficulties for experimental investigation. This work employs a combination of computational and experimental methods to address these challenges. By integrating structure‐informed design, crystallography, and machine learning‐augmented all‐atom molecular dynamics simulations (MD), we synthesized, biophysically and biochemically characterized, and studied the dissociation of a library of small molecule activators of the 5‐aminoimidazole–4–carboxamide ribonucleotide triphosphate (ZTP) riboswitch, a ligand‐binding RNA motif that regulates bacterial gene expression. We uncovered key interaction mechanisms, revealing valuable insights into the role of ligand binding kinetics on riboswitch activation. Further, we established that ligand on‐rates determine activation potency as opposed to binding affinity and elucidated RNA structural differences, which provide mechanistic insights into the interplay of RNA structure on riboswitch activation.
Not Available
Transcription-coupled repair is essential for the removal of DNA lesions from the transcribed genome. The pathway is initiated by CSB protein binding to stalled RNA polymerase II. Mutations impairing CSB function cause severe genetic disease. Yet, the ATP-dependent mechanism by which CSB powers RNA polymerase to bypass certain lesions while triggering excision of others is incompletely understood. Here we build structural models of RNA polymerase II bound to the yeast CSB ortholog Rad26 in nucleotide-free and bound states. This enables simulations and graph-theoretical analyses to define partitioning of this complex into dynamic communities and delineate how its structural elements function together to remodel DNA. We identify an allosteric pathway coupling motions of the Rad26 ATPase modules to changes in RNA polymerase and DNA to unveil a structural mechanism for CSB-assisted progression past less bulky lesions. Our models allow functional interpretation of the effects of Cockayne syndrome disease mutations.
In bacterial defense and genome editing applications, the CRISPR-associated protein Cas9 searches millions of DNA base pairs to locate a 20-nucleotide, guide RNA-complementary target sequence that abuts a protospacer-adjacent motif (PAM). Target capture requires Cas9 to unwind DNA at candidate sequences using an unknown ATP-independent mechanism. Here we show that Cas9 sharply bends and undertwists DNA on PAM binding, thereby flipping DNA nucleotides out of the duplex and toward the guide RNA for sequence interrogation. Cryogenic-electron microscopy (cryo-EM) structures of Cas9–RNA–DNA complexes trapped at different states of the interrogation pathway, together with solution conformational probing, reveal that global protein rearrangement accompanies formation of an unstacked DNA hinge. Bend-induced base flipping explains how Cas9 ‘reads’ snippets of DNA to locate target sites within a vast excess of nontarget DNA, a process crucial to both bacterial antiviral immunity and genome editing. This mechanism establishes a physical solution to the problem of complementarity-guided DNA search and shows how interrogation speed and local DNA geometry may influence genome editing efficiency.
T-box riboswitches are unique riboregulators where gene regulation is mediated through interactions between two highly structured RNAs. Despite extensive structural insights, how RNA-RNA interactions drive the folding and structural transitions of T-box to achieve functional conformations remains unclear. Here, by combining SAXS, single-molecule FRET and computational modeling, we elaborate the folding energy landscape of a translational T-box aptamer consisting of stems I, II and IIA/B, which Mg 2+ -induced global folding and tRNA binding are cooperatively coupled. smFRET measurements reveal that high Mg 2+ stabilizes IIA/B and its stacking on II, which drives the pre-docking of I and II into a competent conformation, subsequent tRNA binding promotes docking of I and II to form a high-affinity tRNA binding groove, of which the essentiality of IIA/B and S-turn in II is substantiated with mutational analysis. We highlight a delicate balance among Mg 2+ , the intra- and intermolecular RNA-RNA interactions in modulating RNA folding and function.
RNA polymerase III (Pol III) is responsible for transcribing 5S ribosomal RNA (5S rRNA), tRNAs, and other short non-coding RNAs. Its recruitment to the 5S rRNA promoter requires transcription factors TFIIIA, TFIIIC, and TFIIIB. Here, we use cryoelectron microscopy (cryo-EM) to visualize the S. cerevisiae complex of TFIIIA and TFIIIC bound to the promoter. Gene-specific factor TFIIIA interacts with DNA and acts as an adaptor for TFIIIC-promoter interactions. We also visualize DNA binding of TFIIIB subunits, Brf1 and TBP (TATA-box binding protein), which results in the full-length 5S rRNA gene wrapping around the complex. Our smFRET study reveals that the DNA within the complex undergoes both sharp bending and partial dissociation on a slow timescale, consistent with the model predicted from our cryo-EM results. Our findings provide new insights into the transcription initiation complex assembly on the 5S rRNA promoter and allow us to directly compare Pol III and Pol II transcription adaptations.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
RtcB enzymes comprise a widely distributed family of manganese- and GTP-dependent RNA repair enzymes that join 2',3'-cyclic phosphate ends to 5'-OH ends via RtcB-(histidinyl-N)–GMP, RNA 3'-phosphate, and RNA 3 'pp 5 'G intermediates. RtcB can ligate either 5'-OH RNA or 5'-OH DNA strands in vitro. The nucleic acid contacts of RtcB are uncharted. In this work, we report a 2.7 Å crystal structure of Pyrococcus horikoshii RtcB in complex with a 6-mer 5'-OH DNA oligonucleotide HO A 1 pT 2 pG 3 pT 4 pC 5 pC 6 , which reveals enzymic contacts of Asn202 to the terminal 5'-OH nucleophile; Arg238 to the A1pT2 and T2pG3 phosphates; Arg190 and Gln194 to the T 2 pG 3 phosphate; and an Arg190 π-cation interaction with the G 3 nucleobase. The structural insights affirm functional studies of E. coli RtcB that implicated the conserved counterpart of Arg238 in engagement of the 5'-OH strand for ligation. The essential active site Cys98 that coordinates two manganese ions is oxidized to cysteine sulfonic acid in our structure, raising the prospect that RtcB activity might be sensitive to modulation during oxidative stress.