RNA extension drives a stepwise displacement of an initiation-factor structural module in initial transcription
Explore the source record for details and available documents.
SEARCH · Search NASA
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.
Mechanical anisotropy is an essential property for many biomolecules to assume their structures, functions and applications, however, the mechanisms for their direction-dependent mechanical responses remain elusive. Herein, by using a single-molecule nanopore sensing technique, we explore the mechanisms of directional mechanical stability of the xrRNA1 RNA from ZIKA virus (ZIKV), which forms a complex ring-like architecture. We reveal extreme mechanical anisotropy in ZIKV xrRNA1 which highly depends on Mg 2+ and the key tertiary interactions. The absence of Mg 2+ and disruption of the key tertiary interactions strongly affect the structural integrity and attenuate mechanical anisotropy. The significance of ring structures in RNA mechanical anisotropy is further supported by steered molecular dynamics simulations in combination with force distribution analysis. We anticipate the ring structures can be used as key elements to build RNA-based nanostructures with controllable mechanical anisotropy for biomaterial and biomedical applications.
Microbes commonly organize into communities consisting of hundreds of species involved in complex interactions with each other. 16S ribosomal RNA (16S rRNA) amplicon profiling provides snapshots that reveal the phylogenies and abundance profiles of these microbial communities. These snapshots, when collected from multiple samples, can reveal the co-occurrence of microbes, providing a glimpse into the network of associations in these communities. However, the inference of networks from 16S data involves numerous steps, each requiring specific tools and parameter choices. Moreover, the extent to which these steps affect the final network is still unclear. In this study, we perform a meticulous analysis of each step of a pipeline that can convert 16S sequencing data into a network of microbial associations. Through this process, we map how different choices of algorithms and parameters affect the co-occurrence network and identify the steps that contribute substantially to the variance. We further determine the tools and parameters that generate robust co-occurrence networks and develop consensus network algorithms based on benchmarks with mock and synthetic data sets. The Microbial Co-occurrence Network Explorer, or MiCoNE (available at https://github.com/segrelab/MiCoNE) follows these default tools and parameters and can help explore the outcome of these combinations of choices on the inferred networks. We envisage that this pipeline could be used for integrating multiple data sets and generating comparative analyses and consensus networks that can guide our understanding of microbial community assembly in different biomes.
Paramyxoviruses are enveloped, nonsegmented, negative-strand RNA viruses that cause a wide spectrum of human and animal diseases. The viral genome, packaged by the nucleoprotein (N), serves as a template for the polymerase complex, composed of the large protein (L) and the homo-tetrameric phosphoprotein (P). The ~250-kDa L possesses all enzymatic activities necessary for its function but requires P in vivo. Structural information is available for individual P domains from different paramyxoviruses, but how P interacts with L and how that affects the activity of L is largely unknown due to the lack of high-resolution structures of this complex in this viral family. In this study we determined the structure of the L–P complex from parainfluenza virus 5 (PIV5) at 4.3-Å resolution using cryoelectron microscopy, as well as the oligomerization domain (OD) of P at 1.4-Å resolution using X-ray crystallography. P-OD associates with the RNA-dependent RNA polymerase domain of L and protrudes away from it, while the X domain of one chain of P is bound near the L nucleotide entry site. The methyltransferase (MTase) domain and the C-terminal domain (CTD) of L adopt a unique conformation, positioning the MTase active site immediately above the poly-ribonucleotidyltransferase domain and near the likely exit site for the product RNA 5' end. Our study reveals a potential mechanism that mononegavirus polymerases may employ to switch between transcription and genome replication. This knowledge will assist in the design and development of antivirals against paramyxoviruses.
Prokaryotes and viruses have fought a long battle against each other. Prokaryotes use CRISPR–Cas-mediated adaptive immunity, while conversely, viruses evolve multiple anti-CRISPR (Acr) proteins to defeat these CRISPR–Cas systems. The type I-F CRISPR–Cas system in Pseudomonas aeruginosa requires the crRNA-guided surveillance complex (Csy complex) to recognize the invading DNA. Although some Acr proteins against the Csy complex have been reported, other relevant Acr proteins still need studies to understand their mechanisms. As such, here, we obtain three structures of previously unresolved Acr proteins (AcrF9, AcrF8, and AcrF6) bound to the Csy complex using electron cryo-microscopy (cryo-EM), with resolution at 2.57 Å, 3.42 Å, and 3.15 Å, respectively. The 2.57-Å structure reveals fine details for each molecular component within the Csy complex as well as the direct and water-mediated interactions between proteins and CRISPR RNA (crRNA). Our structures also show unambiguously how these Acr proteins bind differently to the Csy complex. AcrF9 binds to key DNA-binding sites on the Csy spiral backbone. AcrF6 binds at the junction between Cas7.6f and Cas8f, which is critical for DNA duplex splitting. AcrF8 binds to a distinct position on the Csy spiral backbone and forms interactions with crRNA, which has not been seen in other Acr proteins against the Csy complex. Our structure-guided mutagenesis and biochemistry experiments further support the anti-CRISPR mechanisms of these Acr proteins. Our findings support the convergent consequence of inhibiting degradation of invading DNA by these Acr proteins, albeit with different modes of interactions with the type I-F CRISPR–Cas system.
Tpt1 is a widely distributed enzyme that removes an internal RNA 2′-phosphate by transfer to NAD + , via a two-step reaction in which: (i) the RNA 2′-PO 4 attacks NAD + to form an RNA-2′-phospho-(ADP-ribose) intermediate and expel nicotinamide; and (ii) the ADP-ribose O2″ attacks the RNA 2′-phosphodiester to form 2′-OH RNA and ADP-ribose-1″,2″-cyclic phosphate products. Tpt1 can also execute a single-step ADP-ribosyltransferase reaction at a 5′-monophosphate nucleic acid terminus that installs a 5′-phospho-ADP-ribose cap structure. Here we present crystal structures of Tpt1 bound to an RNA containing an internal 2′-PO 4 mark (the substrate for the canonical Tpt1 pathway) and in a complex with 5′-AMP. We find that Tpt1 has distinct binding modes, whereby the RNA 2′-PO 4 and the AMP 5′-PO 4 are engaged by the same set of active site amino acids, but the 2′-PO 4 nucleoside and the 5′-nucleoside occupy different sites on the enzyme.
Abstract The Mediator complex coordinates regulatory input for transcription driven by RNA polymerase II in eukaryotes. reduced epidermal fluorescence4-3 (ref4-3) is a semidominant mutation that results in a single amino acid substitution in the Mediator tail subunit Med5b. Previous characterization of ref4-3 revealed altered expression of a variety of loci in Arabidopsis, including those contributing to phenylpropanoid biosynthesis. Examination of existing RNA-seq data indicated that loci enriched for the transcriptionally repressive chromatin modification H3K27me3 are overrepresented among genes that are misregulated in ref4-3. We used ChIP-seq and RNA-seq to examine the possibility that perturbation of H3K27me3 homeostasis in ref4-3 plants contributed to altered transcript levels. We observed that ref4-3 results in a modest global reduction of H3K27me3 at enriched loci and that this reduction is not dependent on gene expression; however, altered H3K27me3 was not strongly predictive of altered expression in ref4-3 plants. Instead, our analyses revealed a substantial enrichment of targets of the MYC2 transcriptional regulator among genes that exhibit decreased expression in ref4-3. Consistent with previous characterization of ref4-3, we observed that ref4-3-dependent decreased expression of MYC2 targets can be suppressed by loss of another Mediator tail subunit, MED25. This observation is consistent with previous biochemical characterization of MYC2. Our data highlight the diverse and distinct impacts that a single amino acid change in the tail subunit of Mediator can have on transcriptional circuits and raise the prospect that Mediator directly contributes to H3K27me3 homeostasis in plants.
Recent advancements in RNA structural biology have focused on unraveling the complexities of non-coding mRNA elements like riboswitches. These cis-acting regulatory regions undergo structural changes in response to specific cellular metabolites, leading to up or downregulation of downstream genes. The purine riboswitch family regulates many prokaryotic genes involved in purine degradation and biosynthesis. They feature an aptamer domain organized around a 3-way helical junction, where ligand encapsulation occurs at the junctional core. In our study, we chemically probed the aptamer domain of the 2’-dG-sensing purine riboswitch from Mesoplasma florum (dGsw) under various solution conditions to understand how Mg²⁺ and 2’-dG influence riboswitch folding. Here, we find that efficient 2’-dG binding strongly depends on Mg²⁺, indicating that Mg²⁺ is essential for priming dGsw for ligand interactions. We identified a previously undescribed sequence in the 5’ tail of dGsw that is complementary to a conserved helix. The inclusion of this region in a construct led to intramolecular competition between the alternate helix, Palt, and P1. Mutational analysis confirmed that 5’ flanking end of the aptamer domain forms an alternate helix in the absence of ligand. Molecular dynamics simulations revealed that this alternative conformation is stable. This helix may, therefore, facilitate the formation of an anti-terminator helix by opening the 3-way junction surrounding the 2’-dG binding site. Our study further establishes the importance of a closed terminal P1 helix conformation for metabolite binding and suggests that the delicate interplay between P1 and Palt may fine-tune downstream gene regulation. These insights offer a new perspective on riboswitch structure and enhance our understanding of the role that a conformational ensemble plays in riboswitch activity and regulation.
Dynamic programming algorithms within the NUPACK software suite enable analysis of equilibrium base-pairing properties for complex and test tube ensembles containing arbitrary numbers of interacting nucleic acid strands. Currently, calculations are limited to single-material systems that are either all-RNA or all-DNA. Here, to enable analysis of mixed-material systems that are critical for modern applications in vitro, in situ, and in vivo, we develop physical models and dynamic programming algorithms that allow the material of the system to be specified at nucleotide resolution. Free energy parameter sets are constructed for both RNA/DNA and RNA/2'OMe-RNA mixed-material systems by combining available empirical mixed-material parameters with single-material parameter sets to enable treatment of the full complex and test tube ensembles. New dynamic programming recursions account for the material of each nucleotide throughout the recursive process. For a complex with N nucleotides, the mixed-material dynamic programming algorithms maintain the O(N 3 ) time complexity of the single-material algorithms, enabling efficient calculation of diverse physical quantities over complex and test tube ensembles (e.g., complex partition function, equilibrium complex concentrations, equilibrium base-pairing probabilities, minimum free energy secondary structure(s), and Boltzmann-sampled secondary structures) at a cost increase of roughly 2.0-3.5×. The results of existing single-material algorithms are exactly reproduced when applying the new mixed-material algorithms to single-material systems. Accuracy is significantly enhanced using mixed-material models and algorithms to predict RNA/DNA and RNA/2'OMe-RNA duplex melting temperatures from the experimental literature as well as RNA/DNA melt profiles from new experiments. In conclusion, mixed-material analyses can be performed online using the NUPACK web app (www.nupack.org) or locally using the NUPACK Python module.
Although advances in nanotechnology have enabled the construction of complex and functional synthetic nucleic acid–based nanoarchitectures, high-resolution discrete structures are lacking because of the difficulty in obtaining good diffracting crystals. Here, we report the design and construction of RNA nanostructures based on homooligomerizable one-stranded tiles for x-ray crystallographic determination. We solved three structures to near-atomic resolution: a 2D parallelogram, a 3D nanobracelet unexpectedly formed from an RNA designed for a nanocage, and, eventually, a bona fide 3D nanocage designed with the guidance of the two previous structures. Structural details of their constituent motifs, such as kissing loops, branched kissing loops, and T-junctions, that resemble natural RNA motifs and resisted x-ray determination are revealed, providing insights into those natural motifs. This work unveils the largely unexplored potential of crystallography in gaining high-resolution feedback for nanoarchitectural design and suggests a route to investigate RNA motif structures by configuring them into nanoarchitectures.
hnRNPA2 is a human ribonucleoprotein (RNP) involved in RNA metabolism. It forms fibrils both under cellular stress and in mutated form in neurodegenerative conditions. Previous work established that the C-terminal low-complexity domain (LCD) of hnRNPA2 fibrillizes under stress, and missense mutations in this domain are found in the disease multisystem proteinopathy (MSP). However, little is known at the atomic level about the hnRNPA2 LCD structure that is involved in those processes and how disease mutations cause structural change. Here we present the cryo-electron microscopy (cryoEM) structure of the hnRNPA2 LCD fibril core and demonstrate its capability to form a reversible hydrogel in vitro containing amyloid-like fibrils. Whereas these fibrils, like pathogenic amyloid, are formed from protein chains stacked into β-sheets by backbone hydrogen bonds, they display distinct structural differences: the chains are kinked, enabling non-covalent cross-linking of fibrils and disfavoring formation of pathogenic steric zippers. Both reversibility and energetic calculations suggest these fibrils are less stable than pathogenic amyloid. Moreover, the crystal structure of the disease-mutation-containing segment (D290V) of hnRNPA2 suggests that the replacement fundamentally alters the fibril structure to a more stable energetic state. These findings illuminate how molecular interactions promote protein fibril networks and how mutation can transform fibril structure from functional to a pathogenic form.
Beetroot is a homodimeric in vitro selected RNA that binds and activates DFAME, a conditional fluorophore derived from GFP. It is 70% sequence-identical to the previously characterized homodimeric aptamer Corn, which binds one molecule of its cognate fluorophore DFHO at its interprotomer interface. We have now determined the Beetroot-DFAME co-crystal structure at 1.95 Å resolution, discovering that this RNA homodimer binds two molecules of the fluorophore, at sites separated by ~30 Å. In addition to this overall architectural difference, the local structures of the non-canonical, complex quadruplex cores of Beetroot and Corn are distinctly different, underscoring how subtle RNA sequence differences can give rise to unexpected structural divergence. Through structure-guided engineering, we generated a variant that has a 12-fold fluorescence activation selectivity switch toward DFHO. Beetroot and this variant form heterodimers and constitute the starting point for engineered tags whose through-space inter-fluorophore interaction could be used to monitor RNA dimerization.
Abstract CRISPR-Cas12a systems are becoming an attractive genome editing tool for cell engineering due to their broader editing capabilities compared to CRISPR-Cas9 counterparts. As opposed to Cas9, the Cas12a endonucleases are characterized by a lack of trans-activating crRNA (tracrRNA), which reduces the complexity of the editing system and simultaneously makes CRISPR RNA (crRNA) engineering a promising approach toward further improving and modulating editing activity of the CRISPR-Cas12a systems. Here, we design and validate sixteen types of structurally engineered Cas12a crRNAs targeting various immunologically relevant loci in-vitro and in-cellulo . We show that all our structural modifications in the loop region, ranging from engineered breaks (STAR-crRNAs) to large gaps (Gap-crRNAs), as well as nucleotide substitutions, enable gene-cutting in the presence of various Cas12a nucleases. Moreover, we observe similar insertion rates of short HDR templates using the engineered crRNAs compared to the wild-type crRNAs, further demonstrating that the introduced modifications in the loop region led to comparable genome editing efficiencies. In conclusion, we show that Cas12a nucleases can broadly utilize structurally engineered crRNAs with breaks or gaps in the otherwise highly-conserved loop region, which could further facilitate a wide range of genome editing applications.
CRISPR-Cas systems generally provide adaptive immunity in prokaryotes through RNA-guided degradation of foreign genetic elements like bacteriophages and plasmids. Recently, however, transposon-encoded and nuclease-deficient CRISPR-Cas systems were characterized and shown to be co-opted by Tn7-like transposons for CRISPR RNA-guided DNA transposition. As a genome engineering tool, these CRISPR-Cas systems and their associated transposon proteins can be deployed for programmable, site-specific integration of sizable cargo DNA, circumventing the need for DNA cleavage and homology-directed repair involving endogenous repair machinery. Here, we selected a diverse set of type I-F3 CRISPR-associated transposon systems derived from Gammaproteobacteria, predicted all components essential for transposition activity, and deployed them for functionality testing within Escherichia coli. Our results demonstrate that these systems possess a significant range of integration efficiencies with regards to temperature, transposon size, and flexible PAM requirements. Additionally, our findings support the categorization of these systems into functional compatibility groups for efficient and orthogonal RNA-guided DNA integration. This work expands the CRISPR-based toolbox with new CRISPR RNA-guided DNA integrases that can be applied to complex and extensive genome engineering efforts.
Coronavirus nucleocapsid protein (NP) of SARS-CoV-2 plays a central role in many functions important for virus proliferation including packaging and protecting genomic RNA. The protein shares sequence, structure, and architecture with nucleocapsid proteins from betacoronaviruses. The N-terminal domain (NP RBD ) binds RNA and the C-terminal domain is responsible for dimerization. After infection, NP is highly expressed and triggers robust host immune response. The anti-NP antibodies are not protective and not neutralizing but can effectively detect viral proliferation soon after infection. Two structures of SARS-CoV-2 NP RBD were determined providing a continuous model from residue 48 to 173, including RNA binding region and key epitopes. Five structures of NP RBD complexes with human mAbs were isolated using an antigen-bait sorting. Complexes revealed a distinct complement-determining regions and unique sets of epitope recognition. This may assist in the early detection of pathogens and designing peptide-based vaccines. Mutations that significantly increase viral load were mapped on developed, full length NP model, likely impacting interactions with host proteins and viral RNA.
Bacterial mRNA sequencing is inefficient due to the abundance of ribosomal RNA that is challenging to deplete. While commercial kits target rRNA from common bacterial species, they are frequently inefficient when applied to divergent species, including those from environmental isolates. Similarly, other methods typically employ large probe sets that tile the entire length of rRNAs; however, such approaches are infeasible when applied to many species. Therefore, we present EMBR-seq+, which requires fewer than 10 oligonucleotides per rRNA by combining rRNA blocking primers with RNase H-mediated depletion to achieve rRNA removal efficiencies of up to 99% in diverse bacterial species. Furthermore, in more complex microbial co-cultures between Fibrobacter succinogenes strain UWB7 and anaerobic fungi, EMBR-seq+ depleted both bacterial and fungal rRNA, with a fourfold improvement in bacterial rRNA depletion compared with a commercial kit, thereby demonstrating that the method can be applied to non-model microbial mixtures. Notably, for microbes with unknown rRNA sequences, EMBR-seq+ enables rapid iterations in probe design without requiring to start experiments from total RNA. Finally, efficient depletion of rRNA enabled systematic quantification of the reprogramming of the bacterial transcriptome when cultured in the presence of the anaerobic fungi Anaeromyces robustus or Caecomyces churrovis. We observed that F. succinogenes strain UWB7 downregulated several lignocellulose-degrading carbohydrate-active enzymes in the presence of anaerobic gut fungi, suggesting close interactions between two cellulolytic species that specialize in different aspects of biomass breakdown. Thus, EMBR-seq+ enables efficient, cost-effective, and rapid quantification of the transcriptome to gain insights into non-model microbial systems.
The application of RNA-sequencing has led to numerous breakthroughs related to investigating gene expression levels in complex biological systems. Among these are knowledge of how organisms, such as the vertebrate model organism zebrafish (Danio rerio), respond to toxicant exposure. Recently, the development of 3' RNA-seq has allowed for the determination of gene expression levels with a fraction of the required reads compared to standard RNA-seq. While 3' RNA-seq has many advantages, a comparison to standard RNA-seq has not been performed in the context of whole organism toxicity and sparse data. Here, we examined samples from zebrafish exposed to perfluorobutane sulfonamide (FBSA) with either 3' or standard RNA-seq to determine the advantages of each with regards to the identification of functionally enriched pathways. We found that 3' and standard RNA-seq showed specific advantages when focusing on annotated or unannotated regions of the genome. We also found that standard RNA-seq identified more differentially expressed genes (DEGs), but that this advantage disappeared under conditions of sparse data. We also found that standard RNA-seq had a significant advantage in identifying functionally enriched pathways via analysis of DEG lists but that this advantage was minimal when identifying pathways via gene set enrichment analysis of all genes. These results show that each approach has experimental conditions where they may be advantageous. Our observations can help guide others in the choice of 3' RNA-seq vs standard RNA sequencing to query gene expression levels in a range of biological systems.
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.