RNA Base Pairs Research Articles

Nonequilibrium phases of a biomolecular condensate facilitated by enzyme activity.

Biomolecular condensates represent a frontier in cellular organization, existing as dynamic materials driven out of equilibrium by active cellular processes. Here we explore active mechanisms of condensate regulation by examining the interplay between DEAD-box helicase activity and RNA base-pairing interactions within ribonucleoprotein condensates. We demonstrate how the ATP-dependent activity of DEAD-box helicases-a key class of enzymes in condensate regulation-acts as a nonequilibrium driver of condensate properties through the continuous remodeling of RNA interactions. By combining the LAF-1 DEAD-box helicase with a designer RNA hairpin concatemer, we unveil a complex landscape of dynamic behaviors, including time-dependent alterations in RNA partitioning, evolving condensate morphologies, and shifting condensate dynamics. Importantly, we reveal an antagonistic relationship between RNA secondary structure and helicase activity which promotes condensate homogeneity via a nonequilibrium steady state. By elucidating these nonequilibrium mechanisms, we gain a deeper understanding of cellular organization and expand the potential for active synthetic condensate systems.

Discovery and Quantification of Long-Range RNA Base Pairs in Coronavirus Genomes with SEARCH-MaP and SEISMIC-RNA.

RNA molecules perform a diversity of essential functions for which their linear sequences must fold into higher-order structures. Techniques including crystallography and cryogenic electron microscopy have revealed 3D structures of ribosomal, transfer, and other well-structured RNAs; while chemical probing with sequencing facilitates secondary structure modeling of any RNAs of interest, even within cells. Ongoing efforts continue increasing the accuracy, resolution, and ability to distinguish coexisting alternative structures. However, no method can discover and quantify alternative structures with base pairs spanning arbitrarily long distances - an obstacle for studying viral, messenger, and long noncoding RNAs, which may form long-range base pairs. Here, we introduce the method of Structure Ensemble Ablation by Reverse Complement Hybridization with Mutational Profiling (SEARCH-MaP) and software for Structure Ensemble Inference by Sequencing, Mutation Identification, and Clustering of RNA (SEISMIC-RNA). We use SEARCH-MaP and SEISMIC-RNA to discover that the frameshift stimulating element of SARS coronavirus 2 base-pairs with another element 1 kilobase downstream in nearly half of RNA molecules, and that this structure competes with a pseudoknot that stimulates ribosomal frameshifting. Moreover, we identify long-range base pairs involving the frameshift stimulating element in other coronaviruses including SARS coronavirus 1 and transmissible gastroenteritis virus, and model the full genomic secondary structure of the latter. These findings suggest that long-range base pairs are common in coronaviruses and may regulate ribosomal frameshifting, which is essential for viral RNA synthesis. We anticipate that SEARCH-MaP will enable solving many RNA structure ensembles that have eluded characterization, thereby enhancing our general understanding of RNA structures and their functions. SEISMIC-RNA, software for analyzing mutational profiling data at any scale, could power future studies on RNA structure and is available on GitHub and the Python Package Index.

Discovery and Quantification of Long-Range RNA Base Pairs in Coronavirus Genomes with SEARCH-MaP and SEISMIC-RNA.

RNA molecules perform a diversity of essential functions for which their linear sequences must fold into higher-order structures. Techniques including crystallography and cryogenic electron microscopy have revealed 3D structures of ribosomal, transfer, and other well-structured RNAs; while chemical probing with sequencing facilitates secondary structure modeling of any RNAs of interest, even within cells. Ongoing efforts continue increasing the accuracy, resolution, and ability to distinguish coexisting alternative structures. However, no method can discover and quantify alternative structures with base pairs spanning arbitrarily long distances - an obstacle for studying viral, messenger, and long noncoding RNAs, which may form long-range base pairs. Here, we introduce the method of Structure Ensemble Ablation by Reverse Complement Hybridization with Mutational Profiling (SEARCH-MaP) and software for Structure Ensemble Inference by Sequencing, Mutation Identification, and Clustering of RNA (SEISMIC-RNA). We use SEARCH-MaP and SEISMIC-RNA to discover that the frameshift stimulating element of SARS coronavirus 2 base-pairs with another element 1 kilobase downstream in nearly half of RNA molecules, and that this structure competes with a pseudoknot that stimulates ribosomal frameshifting. Moreover, we identify long-range base pairs involving the frameshift stimulating element in other coronaviruses including SARS coronavirus 1 and transmissible gastroenteritis virus, and model the full genomic secondary structure of the latter. These findings suggest that long-range base pairs are common in coronaviruses and may regulate ribosomal frameshifting, which is essential for viral RNA synthesis. We anticipate that SEARCH-MaP will enable solving many RNA structure ensembles that have eluded characterization, thereby enhancing our general understanding of RNA structures and their functions. SEISMIC-RNA, software for analyzing mutational profiling data at any scale, could power future studies on RNA structure and is available on GitHub and the Python Package Index.

ChimericFragments: computation, analysisand visualization of global RNA networks.

RNA-RNA interactions are a key feature of post-transcriptional gene regulation in all domains of life. While ever more experimental protocols are being developed to study RNA duplex formation on a genome-wide scale, computational methods for the analysis and interpretation of the underlying data are lagging behind. Here, we present ChimericFragments, an analysis framework for RNA-seq experiments that produce chimeric RNA molecules. ChimericFragments implements a novel statistical method based on the complementarity of the base-pairing RNAs around their ligation site and provides an interactive graph-based visualization for data exploration and interpretation. ChimericFragments detects true RNA-RNA interactions with high precision and is compatible with several widely used experimental procedures such as RIL-seq, LIGR-seqor CLASH. We further demonstrate that ChimericFragments enables the systematic detection of novel RNA regulators and RNA-target pairs with crucial roles in microbial physiology and virulence. ChimericFragments is written in Julia and available at: https://github.com/maltesie/ChimericFragments.

The Tautomeric State of N4-Hydroxycytidine within Base-Paired RNA.

Antiviral nucleoside analogues (e.g., Molnupiravir, Remdesivir) played key roles in the treatment of COVID-19 by targeting SARS-CoV-2 RNA-dependent RNA polymerase (RdRp). The nucleoside of Molnupiravir, N4-hydroxycytidine (NHC), exists in two tautomeric forms that pair either with G or A within the RdRp active site, causing an accumulation of viral RNA mutations during replication. Detailed insights into the tautomeric states within base pairs and the structural influence of NHC in RNA are still missing. In this study, we investigate the properties of NHC:G and NHC:A base pairs in a self-complementary RNA duplex by UV thermal melting and NMR spectroscopy using atom-specifically 15N-labeled versions of NHC that were incorporated into oligonucleotides by solid-phase synthesis. NMR analysis revealed that NHC forms a Watson-Crick base pair with G via its amino form, whereas two equally populated conformations were detected for the NHC:A base pair: a weakly hydrogen-bonded Watson-Crick base pair with NHC in the imino form and another conformation with A shifted toward the minor groove. Moreover, we found a variable influence of NHC:G and NHC:A base pairs on the neighboring duplex environment. This study provides conclusive experimental evidence for the existence of two tautomeric forms of NHC within RNA base pairs.

An Integrated Affinity Chromatography-Based Approach to Unravel the sRNA Interactome in Nitrogen-Fixing Rhizobia.

The activity mechanism and function of bacterial base-pairing small non-coding RNA regulators (sRNAs) are largely shaped by their main interacting cellular partners, i.e., proteins and mRNAs. We describe here an MS2 affinity chromatography-based procedure adapted to unravel the sRNA interactome in nitrogen-fixing legume endosymbiotic bacteria. The method consists of tagging of the bait sRNA at its 5'-end with the MS2 aptamer followed by pulse overexpression and immobilization of the chimeric transcript from cell lysates by an MS2-MBP fusion protein conjugated to an amylose resin. The sRNA-binding proteins and target mRNAs are further profiled by mass spectrometry and RNAseq, respectively.

Evaluation of 5'-End Phosphorylation for Small RNA Stability and Target Regulation In Vivo.

Bacterial small RNAs (sRNAs) can be equipped at the 5' end with triphosphate (5'PPP) or monophosphate (5'P) groups, depending on whether they are primary transcripts, undergo dephosphorylation or originate via processing. Often, 5' groups hallmark RNAs for rapid decay, but whether this also applies to sRNAs is little explored. Moreover, the sRNA 5'P group could activate endoribonuclease RNase E to cleave the base-paired target RNA, but a tool for investigation in vivo was lacking. Here, we describe a two-plasmid system suitable for the generation of 5' monophosphorylated RNAs on demand inside the cell. The sRNA gene of interest is fused to the 3' end of a fragment of sRNA GlmZ and transcribed from a plasmid in an IPTG-inducible manner. The fusion RNA gets cleaved upon arabinose-controlled expression of rapZ, provided on a compatible plasmid. Adaptor protein RapZ binds the GlmZ aptamer and directs RNase E to release the sRNA of choice with 5'P ends. An isogenic plasmid generating the same sRNA with a 5'PPP end allows for direct comparison. The fates of the sRNA variants and target RNA(s) are monitored by Northern blotting. This tool is applicable to E. coli and likely other enteric bacteria.

G-U base-paired hpRNA confers potent inhibition of small RNA function in plants.

MicroRNA (miRNA) target mimicry technologies, utilizing naturally occurring miRNA decoy molecules, represent a potent tool for analyzing miRNA function. In this study, we present a highly efficient small RNA (sRNA) target mimicry design based on G-U base-paired hairpin RNA (hpG:U), which allows for the simultaneous targeting of multiple sRNAs. The hpG:U constructs consistently generate high amounts of intact, polyadenylated stem-loop (SL) RNA outside the nuclei, in contrast to traditional hairpin RNA designs with canonical base pairing (hpWT), which were predominantly processed resulting in a loop. By incorporating a 460-bp G-U base-paired double-stranded stem and a 312-576 nt loop carrying multiple miRNA target mimicry sites (GUMIC), the hpG:U construct displayed effective repression of three Arabidopsis miRNAs, namely miR165/166, miR157, and miR160, both individually and in combination. Additionally, a GUMIC construct targeting a prominent cluster of siRNAs derived from cucumber mosaic virus (CMV) Y-satellite RNA (Y-Sat) effectively inhibited Y-Sat siRNA-directed silencing of the chlorophyll biosynthetic gene CHLI, thereby reducing the yellowing symptoms in infected Nicotiana plants. Therefore, the G-U base-paired hpRNA, characterized by differential processing compared to traditional hpRNA, acts as an efficient decoy for both miRNAs and siRNAs. This technology holds great potential for sRNA functional analysis and the management of sRNA-mediated diseases.

4.5SH RNA counteracts deleterious exonization of SINE B1 in mice

4.5SH RNA is a highly abundant, small rodent-specific noncoding RNA that localizes to nuclear speckles enriched in pre-mRNA-splicing regulators. To investigate the physiological functions of 4.5SH RNA, we have created mutant mice that lack the expression of 4.5SH RNA. The mutant mice exhibited embryonic lethality, suggesting that 4.5SH RNA is an essential species-specific noncoding RNA in mice. RNA-sequencing analyses revealed that 4.5SH RNA protects the transcriptome from abnormal exonizations of the antisense insertions of the retrotransposon SINE B1 (asB1), which would otherwise introduce deleterious premature stop codons or frameshift mutations. Mechanistically, 4.5SH RNA base pairs with complementary asB1-containing exons via the target recognition region and recruits effector proteins including Hnrnpm via its 5' stem loop region. The modular organization of 4.5SH RNA allows us to engineer a programmable splicing regulator to induce the skipping of target exons of interest. Our results also suggest the general existence of splicing regulatory noncoding RNAs.

In vivo-like nearest neighbor parameters improve prediction of fractional RNA base-pairing in cells.

We conducted a thermodynamic analysis of RNA stability in Eco80 artificial cytoplasm, which mimics in vivo conditions, and compared it to transcriptome-wide probing of mRNA. Eco80 contains 80% of Escherichia coli metabolites, with biological concentrations of metal ions, including 2 mM free Mg2+ and 29 mM metabolite-chelated Mg2+. Fluorescence-detected binding isotherms (FDBI) were used to conduct a thermodynamic analysis of 24 RNA helices and found that these helices, which have an average stability of -12.3 kcal/mol, are less stable by ΔΔGo37 ∼1 kcal/mol. The FDBI data was used to determine a set of Watson-Crick free energy nearestneighbor parameters (NNPs), which revealed that Eco80 reduces the stability of three NNPs. This information was used to adjust the NN model using the RNAstructure package. The in vivo-like adjustments have minimal effects on the prediction of RNA secondary structures determined in vitro and in silico, but markedly improve prediction of fractional RNA base pairing in E. coli, as benchmarked with our in vivo DMS and EDC RNA chemical probing data. In summary, our thermodynamic and chemical probing analyses of RNA helices indicate that RNA secondary structures are less stable in cells than in artificially stable in vitro buffer conditions.

A conserved isoleucine in the binding pocket of RIG-I controls immune tolerance to mitochondrial RNA.

RIG-I is a cytosolic receptor of viral RNA essential for the immune response to numerous RNA viruses. Accordingly, RIG-I must sensitively detect viral RNA yet tolerate abundant self-RNA species. The basic binding cleft and an aromatic amino acid of the RIG-I C-terminal domain(CTD) mediate high-affinity recognition of 5'triphosphorylated and 5'base-paired RNA(dsRNA). Here, we found that, while 5'unmodified hydroxyl(OH)-dsRNA demonstrated residual activation potential, 5'-monophosphate(5'p)-termini, present on most cellular RNAs, prevented RIG-I activation. Determination of CTD/dsRNA co-crystal structures and mutant activation studies revealed that the evolutionarily conserved I875 within the CTD sterically inhibits 5'p-dsRNA binding. RIG-I(I875A) was activated by both synthetic 5'p-dsRNA and endogenous long dsRNA within the polyA-rich fraction of total cellular RNA. RIG-I(I875A) specifically interacted with long, polyA-bearing, mitochondrial(mt) RNA, and depletion of mtRNA from total RNA abolished its activation. Altogether, our study demonstrates that avoidance of 5'p-RNA recognition is crucial to prevent mtRNA-triggered RIG-I-mediated autoinflammation.

Effect of Water Models on The Stability of RNA: Role of Counter-Ions

Various force fields and water model potentials influence significantly RNA conformations. The polyanionic nature of RNA attracts the water molecules and the counter ions which in turn affects their stability. The interfacial water's structural and dynamic aspects affect the RNA's base-pair opening and denaturation by breaking or making inter/intra-hydrogen bonds. Herein, we employed an MD simulations study using SPC/E and modified TIP3P water models in combination with different force fields CHARMM and AMBER to find their influence on the hydration shell of the SARS-CoV-2 RNA genome at different temperatures. AMBER-mTIP3P model was found to give more dynamic and transient conformations for RNA. The lower dielectric constant of the SPC/E model helps in the formation of the ion-contact pair near the negatively charged phosphate group (Na+-PO4−) leading to strong RNA-ion interaction and strong hydration shells having higher hydrogen bond lifetime. The Na+ ion survival probability at the interface was found to be more in the SPC/E model. At lower temperatures, the water molecules inside these hydration shells were found to be inhomogeneous, with lower void space, higher-coordinated, and non-tetrahedral. The higher dielectric constant of the mTIP3P model screened out the attraction between the ion pairs leading to a more homogenous solvation shell having a lesser hydrogen bond lifetime and more diffusive water. The distribution of the ions near the RNA structure is confirmed by metadynamics simulations. Both water models were found to disrupt the base pair orientation due to the formation of water bridges between the O2ʹ group of RNA and the water molecules.

Analyzes In Silico Indicate the lncRNAs MIR31HG and LINC00939 as Possible Epigenetic Inhibitors of the Osteogenic Differentiation in PDLCs.

Chromatin conformation, DNA methylation pattern, transcriptional profile, and non-coding RNAs (ncRNAs) interactions constitute an epigenetic pattern that influences the cellular phenotypic commitment and impacts the clinical outcomes in regenerative therapies. Here, we investigated the epigenetic landscape of the SP7 transcriptor factor (SP7) and Distal-Less Homeobox 4 (DLX4) osteoblastic transcription factors (TFs), in human periodontal ligament mesenchymal cells (PDLCs) with low (l-PDLCs) and high (h-PDLCs) osteogenic potential. Chromatin accessibility (ATAC-seq), genome DNA methylation (Methylome), and RNA sequencing (RNA-seq) assays were performed in l- and h-PDLCs, cultured at 10 days in non-induced (DMEM) and osteogenic (OM) medium in vitro. Data were processed in HOMER, Genome Studio, and edgeR programs, and metadata was analyzed by online bioinformatics tools and in R and Python environments. ATAC-seq analyses showed the TFs genomic regions are more accessible in l-PDLCs than in h-PDLCs. In Methylome analyses, the TFs presented similar average methylation intensities (AMIs), without differently methylated probes (DMPs) between l- and h-PDLCs; in addition, there were no differences in the expression profiles of TFs signaling pathways. Interestingly, we identified the long non-coding RNAs (lncRNAs), MIR31HG and LINC00939, as upregulated in l-PDLCs, in both DMEM and OM. In the following analysis, the web-based prediction tool LncRRIsearch predicted RNA:RNA base-pairing interactions between SP7, DLX4, MIR31HG, and LINC00939 transcripts. The machine learning program TriplexFPP predicted DNA:RNA triplex-forming potential for the SP7 DNA site and for one of the LINC00939 transcripts (ENST00000502479). PCR data confirmed the upregulation of MIR31HG and LINC00939 transcripts in l-PDLCs (× h-PDLCs) in both DMEM and OM (p < 0.05); conversely, SP7 and DLX4 were downregulated, confirming those results observed in the RNA-Seq analysis. Together, these results indicate the lncRNAs MIR31HG and LINC00939 as possible epigenetic inhibitors of the osteogenic differentiation in PDLCs by (post)transcriptional and translational repression of the SP7 and DLX4 TFs.

Crosstalk between Long Non-Coding RNA and Spliceosomal microRNA as a Novel Biomarker for Cancer.

Non-coding RNAs (ncRNAs) play diverse roles in regulating cellular processes and have been implicated in pathological conditions, including cancer, where interactions between ncRNAs play a role. Relevant here are (i) microRNAs (miRNAs), mainly known as negative regulators of gene expression in the cytoplasm. However, identification of miRNAs in the nucleus suggested novel nuclear functions, and (ii) long non-coding RNA (lncRNA) regulates gene expression at multiple levels. The recent findings of miRNA in supraspliceosomes of human breast and cervical cancer cells revealed new candidates of lncRNA targets. Here, we highlight potential cases of crosstalk between lncRNA and supraspliceosomal miRNA expressed from the same genomic region, having complementary sequences. Through RNA:RNA base pairing, changes in the level of one partner (either miRNA or lncRNA), as occur in cancer, could affect the level of the other, which might be involved in breast and cervical cancer. An example is spliceosomal mir-7704 as a negative regulator of the oncogenic lncRNA HAGLR. Because the expression of spliceosomal miRNA is cell-type-specific, the list of cis-interacting lncRNA:spliceosomal miRNA presented here is likely just the tip of the iceberg, and such interactions are likely relevant to additional cancers. We thus highlight the potential of lncRNA:spliceosomal miRNA interactions as novel targets for cancer diagnosis and therapies.

RNA covariation at helix-level resolution for the identification of evolutionarily conserved RNA structure.

Many biologically important RNAs fold into specific 3D structures conserved through evolution. Knowing when an RNA sequence includes a conserved RNA structure that could lead to new biology is not trivial and depends on clues left behind by conservation in the form of covariation and variation. For that purpose, the R-scape statistical test was created to identify from alignments of RNA sequences, the base pairs that significantly covary above phylogenetic expectation. R-scape treats base pairs as independent units. However, RNA base pairs do not occur in isolation. The Watson-Crick (WC) base pairs stack together forming helices that constitute the scaffold that facilitates the formation of the non-WC base pairs, and ultimately the complete 3D structure. The helix-forming WC base pairs carry most of the covariation signal in an RNA structure. Here, I introduce a new measure of statistically significant covariation at helix-level by aggregation of the covariation significance and covariation power calculated at base-pair-level resolution. Performance benchmarks show that helix-level aggregated covariation increases sensitivity in the detection of evolutionarily conserved RNA structure without sacrificing specificity. This additional helix-level sensitivity reveals an artifact that results from using covariation to build an alignment for a hypothetical structure and then testing the alignment for whether its covariation significantly supports the structure. Helix-level reanalysis of the evolutionary evidence for a selection of long non-coding RNAs (lncRNAs) reinforces the evidence against these lncRNAs having a conserved secondary structure.

Insight into the intercalating ligands affecting the interactions of ruthenium(II) polypyridyl complexes with RNA poly(A)*poly(U)

A series of ruthenium(II) complexes [Ru(phen)2dppp3]2+ (Ru1, phen = 1,10-phenanthrolline; dppp3 = pyrido[3′,4′:5,6]pyrazino[2,3-f][1,10]phenanthroline), [Ru(phen)2dppp2]2+ (Ru2, dppp2 = pyrido[2′,3′:5,6]pyrazino[2,3-f][1,10]phenanthroline) and [Ru(phen)2bppp]2+ (Ru3, bppp = 10-bromopyrido[2′,3′:5,6]pyrazino-[2,3-f][1,10]phenanthroline) have been synthesized and characterized that helps to further investigate RNA-complex interactions. The viscosity measurement, UV–vis absorption spectrum study as well as EB (ethidium bromide) fluorescence competition have been explored to prove that the binding modes of Ru1, Ru2 and Ru3 with double-stranded RNA were intercalation, and these binding affinities of Ru1, Ru2 and Ru3 toward poly(A)*poly(U) (where “*” denotes the Watson-Crick base) have followed the order of Ru1 > Ru2 > Ru3. The analysis of thermal denaturation results obtained the identical conclusion and revealed that three complexes have positive roles in the stabilization of RNA base pairs. To my knowledge, this study presents the first Ru(Ⅱ) polypyridine complexes with the most effective thermal stabilization of RNA poly(A)*poly(U). All results showed the structural effects of complex main ligands in RNA study, in which the RNA-affinities were impacted greatly by the positions of substituents, and the sizes of functional groups were also key factors controlling the interaction between the complexes and double-stranded RNA.

Thermodynamics and kinetics of an A-U RNA base pair under force studied by molecular dynamics simulations.

Mechanical force has been widely used to study RNA folding and unfolding. Understanding how the force affects the opening and closing of a single base pair, which is a basic step for RNA folding and unfolding and a fundamental behavior in some important biological activities, is crucial to understanding the mechanism of RNA folding and unfolding under mechanical force. In this work, we investigated the opening and closing process of an RNA base pair under mechanical force with constant-force stretching molecular dynamics simulations. It was found that high mechanical force results in overstretching, and the open state is a high-energy state. The enthalpy and entropy change of the base-pair opening-closing transition were obtained and the results at low forces were in good agreement with the nearest-neighbor model. The temperature and force dependence of the opening and closing rates were also obtained. The position of the transition state for the base-pair opening-closing transition under mechanical force was determined. The free energy barrier of opening a base pair without force is the enthalpy increase, and the work done by the force from the closed state to the transition state decreases the barrier and increases the opening rate. The free energy barrier of closing the base pair without force results from the entropy loss, and the work done by the force from the open state to the transition state increases the barrier and decreases the closing rate. The transition rates are strongly dependent on the temperature and force, while the transition path times are weakly dependent on force and temperature.

Molecular insights into the stereospecificity of arginine in RNA tetraloop folding.

One of the hypotheses for the homochirality of amino acids in the context of the origin of life is that only a particular stereoisomer provides preferential stability to RNA folding by acting as a chemical chaperon. However, the effect at the molecular level is not well understood. This study provides a molecular understanding of such preferential stability for a small GAAA RNA tetraloop in the presence of chiral arginine through a multidimensional free energy landscape constructed using a combination of umbrella sampling and parallel bias metadynamics (PBMetaD) simulations. We show that the origin of the chirality difference in RNA folding-unfolding dynamics is due to differences in the configurational diversity of RNA in adopting various non-natural conformations that accompany the diverse binding modes of D-arginine and L-arginine. We show that while D-arginine stabilizes the native folded state of RNA, L-arginine destabilizes it. Furthermore, free energy calculations on the binding of D- and L-arginine reveal a specific geometric constraint that helps D-arginine to stack with the terminal base pairs of RNA and pushes L-arginine for groove binding.

SNPD-CRISPR: Single Nucleotide Polymorphism-Distinguishable Repression or Enhancement of a Target Gene Expression by CRISPR System.

A wide range of diseases, including cancer, autoimmune diseases, or neurodegenerative diseases, have been associated with single nucleotide mutations in their causative genes. Clustered regularly interspaced short palindromic repeats/CRISPR-associated protein 9 (CRISPR/Cas9) system is a flexible and efficient genome engineering technology widely used for researches and therapeutic applications which offers immense opportunity to treat genetic diseases. The complex of Cas9 and the guide RNA acts as an RNA-guided endonuclease. Cas9 recognizes a sequence motif known as a protospacer adjacent motif (PAM), and then the guide RNA base pairs with its proximal target region of 20 nucleotides with sequence complementarity. Here we describe the procedure named single nucleotide polymorphism-distinguishable (SNPD)-CRISPR system which can suppress or enhance the expression of disease-causative gene with single nucleotide mutation distinguished from its wild-type. In this study, we used HRAS, one of most famous cancer-causative genes, as an example of a target gene.

Substrate Specificity and Kinetics of RNA Hydrolysis by SARS-CoV-2 NSP10/14 Exonuclease.

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), the virus that causes COVID-19, continues to evolve resistance to vaccines and existing antiviral therapies at an alarming rate, increasing the need for new direct-acting antiviral drugs. Despite significant advances in our fundamental understanding of the kinetics and mechanism of viral RNA replication, there are still open questions regarding how the proofreading exonuclease (NSP10/NSP14 complex) contributes to replication fidelity and resistance to nucleoside analogs. Through single turnover kinetic analysis, we show that the preferred substrate for the exonuclease is double-stranded RNA without any mismatches. Double-stranded RNA containing a 3'-terminal remdesivir was hydrolyzed at a rate similar to a correctly base-paired cognate nucleotide. Surprisingly, single-stranded RNA or duplex RNA containing a 3'-terminal mismatch was hydrolyzed at rates 125- and 45-fold slower, respectively, compared to the correctly base-paired double-stranded RNA. These results define the substrate specificity and rate of removal of remdesivir for the exonuclease and outline rigorous kinetic assays that could help in finding next-generation exonuclease inhibitors or nucleoside analogs that are able to evade excision. These results also raise important questions about the role of the polymerase/exonuclease complex in proofreading during viral replication. Addressing these questions through rigorous kinetic analysis will facilitate the search for desperately needed antiviral drugs to combat COVID-19.