Evolutionary Trends in RNA Base Selectivity Within the RNase A Superfamily

Dynamics of non-covalent interactions during the P-O bond cleavage reaction by ribonuclease A

In this work, an atomistic-scale investigation of the phosphodiester P-O bond cleavage reaction by the enzyme ribonuclease A was carried out using computer simulation techniques. It is shown that during the reaction the network of non-covalent interactions in the active center of the ribonuclease changes significantly, while the role of these non-covalent interactions is different: coordination of the corresponding groups, electron density transfer, and ligand holding in the active center. It is shown that the process of proton transfer from Asp121 to His119 is the first stage of this reaction; at the same time, the hydrogen bond between the phosphate ligand and the imino group of Arg39 is broken, which, although keeping the ligand in the active center, does not allow the ligand to orient itself more conveniently for subsequent proton transfers. Furthermore, the key step of this reaction occurs: proton transfer with the participation of imidazole rings His12 and His119, in which the guiding role is played by several hydrogen bonds with the participation of Phe120, and the role of an electron density carrier is played by the pnictogen bond between the oxygen of the phosphate ligand and the pyridine-like nitrogen of the imidazole ring His119, which was detected for the first time.

Analysis of RNA and Its Modifications

Ribonucleic acids (RNAs) are key biomolecules responsible for the transmission of genetic information, the synthesis of proteins, and modulation of many biochemical processes. They are also often the key components of viruses. Synthetic RNAs or oligoribonucleotides are becoming more widely used as therapeutics. In many cases, RNAs will be chemically modified, either naturally via enzymatic systems within a cell or intentionally during their synthesis. Analytical methods to detect, sequence, identify, and quantify RNA and its modifications have demands that far exceed requirements found in the DNA realm. Two complementary platforms have demonstrated their value and utility for the characterization of RNA and its modifications: mass spectrometry and next-generation sequencing. This review highlights recent advances in both platforms, examines their relative strengths and weaknesses, and explores some alternative approaches that lie at the horizon.

Selection on synonymous sites: the unwanted transcript hypothesis

Although translational selection to favour codons that match the most abundant tRNAs is not readily observed in humans, there is nonetheless selection in humans on synonymous mutations. We hypothesize that much of this synonymous site selection can be explained in terms of protection against unwanted RNAs - spurious transcripts, mis-spliced forms or RNAs derived from transposable elements or viruses. We propose not only that selection on synonymous sites functions to reduce the rate of creation of unwanted transcripts (for example, through selection on exonic splice enhancers and cryptic splice sites) but also that high-GC content (but low-CpG content), together with intron presence and position, is both particular to functional native mRNAs and used to recognize transcripts as native. In support of this hypothesis, transcription, nuclear export, liquid phase condensation and RNA degradation have all recently been shown to promote GC-rich transcripts and suppress AU/CpG-rich ones. With such 'traps' being set against AU/CpG-rich transcripts, the codon usage of native genes has, in turn, evolved to avoid such suppression. That parallel filters against AU/CpG-rich transcripts also affect the endosomal import of RNAs further supports the unwanted transcript hypothesis of synonymous site selection and explains the similar design rules that have enabled the successful use of transgenes and RNA vaccines.

Structural determinants for tRNA selective cleavage by RNase 2/EDN

tRNA-derived fragments (tRFs) have emerged as key players of immunoregulation. Some RNase A superfamily members participate in the shaping of the tRFs population. By comparing wild-type and knockout macrophage cell lines, our previous work revealed that RNase 2 can selectively cleave tRNAs. Here, we confirm the in vitro protein cleavage pattern by screening of synthetic tRNAs, single-mutant variants, and anticodon-loop DNA/RNA hairpins. By sequencing of tRF products, we identified the cleavage selectivity of recombinant RNase 2 with base specificity at B1 (U/C) and B2 (A) sites, consistent with a previous cellular study. Lastly, protein-hairpin complexes were predicted by MD simulations. Results reveal the contribution of the α1, loop 3 and loop 4, and β6 RNase 2 regions, where residues Arg36/Asn39/Gln40/Asn65/Arg68/Arg132 provide interactions, spanning from P-1 to P2 sites that are essential for anticodon loop recognition. Knowledge of RNase 2-specific tRFs generation might guide new therapeutic approaches for infectious and immune-related diseases.

A Common Polymorphism in RNASE6 Impacts Its Antimicrobial Activity toward Uropathogenic Escherichia coli

Human Ribonuclease (RNase) 6 is a monocyte and macrophage-derived protein with potent antimicrobial activity toward uropathogenic bacteria. The RNASE6 gene is heterogeneous in humans due to the presence of single nucleotide polymorphisms (SNPs). RNASE6 rs1045922 is the most common non-synonymous SNP, resulting in a G to A substitution that determines an arginine (R) to glutamine (Q) transversion at position 66 in the protein sequence. By structural analysis we observed that R66Q substitution significantly reduces the positive electrostatic charge at the protein surface. Here, we generated both recombinant RNase 6-R66 and -Q66 protein variants and determined their antimicrobial activity toward uropathogenic Escherichia coli (UPEC), the most common cause of UTI. We found that the R66 variant, encoded by the major SNP rs1045922 allele, exhibited superior bactericidal activity in comparison to the Q66 variant. The higher bactericidal activity of R66 variant correlated with an increase in the protein lipopolysaccharide binding and bacterial agglutination abilities, while retaining the same enzymatic efficiency. These findings encourage further work to evaluate RNASE6 SNP distribution and its impact in UTI susceptibility.

Paternal methotrexate exposure affects sperm small RNA content and causes craniofacial defects in the offspring

Folate is an essential vitamin for vertebrate embryo development. Methotrexate (MTX) is a folate antagonist that is widely prescribed for autoimmune diseases, blood and solid organ malignancies, and dermatologic diseases. Although it is highly contraindicated for pregnant women, because it is associated with an increased risk of multiple birth defects, the effect of paternal MTX exposure on their offspring has been largely unexplored. Here, we found MTX treatment of adult medaka male fish (Oryzias latipes) causes cranial cartilage defects in their offspring. Small non-coding RNA (sncRNAs) sequencing in the sperm of MTX treated males identify differential expression of a subset of tRNAs, with higher abundance for specific 5' tRNA halves. Sperm RNA methylation analysis on MTX treated males shows that m5C is the most abundant and differential modification found in RNAs ranging in size from 50 to 90 nucleotides, predominantly tRNAs, and that it correlates with greater testicular Dnmt2 methyltransferase expression. Injection of sperm small RNA fractions from MTX-treated males into normal fertilized eggs generated cranial cartilage defects in the offspring. Overall, our data suggest that paternal MTX exposure alters sperm sncRNAs expression and modifications that may contribute to developmental defects in their offspring.

Impact of food restriction on the medio-basal hypothalamus of intact ewes as revealed by a large-scale transcriptomics study

In mammals, the medio-basal hypothalamus (MBH) integrates photoperiodic and food-related cues to ensure timely phasing of physiological functions, including seasonal reproduction. The current human epidemics of obesity and associated reproductive disorders exemplifies the tight link between metabolism and reproduction. Yet, how food-related cues impact breeding at the level of the MBH remains unclear. In this respect, the sheep, which is a large diurnal mammal with a marked dual photoperiodic/metabolic control of seasonal breeding, is a relevant model. Here, we present a large-scale study in ewes (n = 120), which investigated the impact of food restriction (FRes) on the MBH transcriptome using unbiased RNAseq, followed by RT-qPCR. Few genes (~100) were impacted by FRes and the transcriptional impact was very modest (<2-fold increase or < 50% decrease for most genes). As anticipated, FRes increased expression of Npy/AgRP/LepR and decreased expression of Pomc/Cartpt, while Kiss1 expression was not impacted. Of particular interest, Eya3, Nmu and Dio2, genes involved in photoperiodic decoding within the MBH, were also affected by FRes. Finally, we also identified a handful of genes not known to be regulated by food-related cues (e.g., RNase6, HspA6, Arrdc2). In conclusion, our transcriptomics study provides insights into the impact of metabolism on the MBH in sheep, which may be relevant to human, and identifies possible molecular links between metabolism and (seasonal) reproduction.

Exploring the RNase A scaffold to combine catalytic and antimicrobial activities. Structural characterization of RNase 3/1 chimeras

Design of novel antibiotics to fight antimicrobial resistance is one of the first global health priorities. Novel protein-based strategies come out as alternative therapies. Based on the structure-function knowledge of the RNase A superfamily we have engineered a chimera that combines RNase 1 highest catalytic activity with RNase 3 unique antipathogen properties. A first construct (RNase 3/1-v1) was successfully designed with a catalytic activity 40-fold higher than RNase 3, but alas in detriment of its anti-pathogenic activity. Next, two new versions of the original chimeric protein were created showing improvement in the antimicrobial activity. Both second generation versions (RNases 3/1-v2 and -v3) incorporated a loop characteristic of RNase 3 (L7), associated to antimicrobial activity. Last, removal of an RNase 1 flexible loop (L1) in the third version enhanced its antimicrobial properties and catalytic efficiency. Here we solved the 3D structures of the three chimeras at atomic resolution by X-ray crystallography. Structural analysis outlined the key functional regions. Prediction by molecular docking of the protein chimera in complex with dinucleotides highlighted the contribution of the C-terminal region to shape the substrate binding cavity and determine the base selectivity and catalytic efficiency. Nonetheless, the structures that incorporated the key features related to RNase 3 antimicrobial activity retained the overall RNase 1 active site conformation together with the essential structural elements for binding to the human ribonuclease inhibitor (RNHI), ensuring non-cytotoxicity. Results will guide us in the design of the best RNase pharmacophore for anti-infective therapies.

RNA Cleavage Properties of Nucleobase-Specific RNase MC1 and Cusativin Are Determined by the Dinucleotide-Binding Interactions in the Enzyme-Active Site

Knowledge of the cleavage specificity of ribonucleases is critical for their application in RNA modification mapping or RNA-protein binding studies. Here, we detail the cleavage specificity and efficiency of ribonuclease MC1 and cusativin using a customized RNA sequence that contained all dinucleotide combinations and homopolymer sequences. The sequencing of the oligonucleotide digestion products by a semi-quantitative liquid chromatography coupled with mass spectrometry (LC-MS) analysis documented as little as 0.5-1% cleavage levels for a given dinucleotide sequence combination. While RNase MC1 efficiently cleaved the [A/U/C]pU dinucleotide bond, no cleavage was observed for the GpU bond. Similarly, cusativin efficiently cleaved Cp[U/A/G] dinucleotide combinations along with UpA and [A/U]pU, suggesting a broader specificity of dinucleotide preferences. The molecular interactions between the substrate and active site as determined by the dinucleotide docking studies of protein models offered additional evidence and support for the observed substrate specificity. Targeted alteration of the key amino acid residues in the nucleotide-binding site confirms the utility of this in silico approach for the identification of key interactions. Taken together, the use of bioanalytical and computational approaches, involving LC-MS and ligand docking of tertiary structural models, can form a powerful combination to help explain the RNA cleavage behavior of RNases.

Selective cleavage of ncRNA and antiviral activity by RNase2/EDN in THP1-induced macrophages

RNase2 is the member of the RNaseA family most abundant in macrophages. Here, we knocked out RNase2 in THP-1 cells and analysed the response to Respiratory Syncytial Virus (RSV). RSV induced RNase2 expression, which significantly enhanced cell survival. Next, by cP-RNAseq sequencing, which amplifies the cyclic-phosphate endonuclease products, we analysed the ncRNA population. Among the ncRNAs accumulated in WT vs KO cells, we found mostly tRNA-derived fragments (tRFs) and second miRNAs. Differential sequence coverage identified tRFs from only few parental tRNAs, revealing a predominant cleavage at anticodon and d-loops at U/C (B1) and A (B2) sites. Selective tRNA cleavage was confirmed in vitro using the recombinant protein. Likewise, only few miRNAs were significantly more abundant in WT vs RNase2-KO cells. Complementarily, by screening of a tRF & tiRNA array, we identified an enriched population associated to RNase2 expression and RSV exposure. The results confirm the protein antiviral action and provide the first evidence of its cleavage selectivity on ncRNAs.

Graphical abstract

Characterization of SARS2 Nsp15 nuclease activity reveals it's mad about U

Nsp15 is a uridine specific endoribonuclease that coronaviruses employ to cleave viral RNA and evade host immune defense systems. Previous structures of Nsp15 from across Coronaviridae revealed that Nsp15 assembles into a homo-hexamer and has a conserved active site similar to RNase A. Beyond a preference for cleaving RNA 3′ of uridines, it is unknown if Nsp15 has any additional substrate preferences. Here, we used cryo-EM to capture structures of Nsp15 bound to RNA in pre- and post-cleavage states. The structures along with molecular dynamics and biochemical assays revealed critical residues involved in substrate specificity, nuclease activity, and oligomerization. Moreover, we determined how the sequence of the RNA substrate dictates cleavage and found that outside of polyU tracts, Nsp15 has a strong preference for purines 3′ of the cleaved uridine. This work advances our understanding of how Nsp15 recognizes and processes viral RNA, and will aid in the development of new anti-viral therapeutics.

Human RNase3 immune modulation by catalytic-dependent and independent modes in a macrophage-cell line infection model

The human RNase3 is a member of the RNaseA superfamily involved in host immunity. RNase3 is expressed by leukocytes and shows broad-spectrum antimicrobial activity. Together with a direct antimicrobial action, RNase3 exhibits immunomodulatory properties. Here, we have analysed the transcriptome of macrophages exposed to the wild-type protein and a catalytic-defective mutant (RNase3-H15A). The analysis of differently expressed genes (DEGs) in treated THP1-derived macrophages highlighted a common pro-inflammatory "core-response" independent of the protein ribonucleolytic activity. Network analysis identified the epidermal growth factor receptor (EGFR) as the main central regulatory protein. Expression of selected DEGs and MAPK phosphorylation were inhibited by an anti-EGFR antibody. Structural analysis suggested that RNase3 activates the EGFR pathway by direct interaction with the receptor. Besides, we identified a subset of DEGs related to the protein ribonucleolytic activity, characteristic of virus infection response. Transcriptome analysis revealed an early pro-inflammatory response, not associated to the protein catalytic activity, followed by a late activation in a ribonucleolytic-dependent manner. Next, we demonstrated that overexpression of macrophage endogenous RNase3 protects the cells against infection by Mycobacterium aurum and the human respiratory syncytial virus. Comparison of cell infection profiles in the presence of Erlotinib, an EGFR inhibitor, revealed that the receptor activation is required for the antibacterial but not for the antiviral protein action. Moreover, the DEGs related and unrelated to the protein catalytic activity are associated to the immune response to bacterial and viral infection, respectively. We conclude that RNase3 modulates the macrophage defence against infection in both catalytic-dependent and independent manners.

Immune Sensing of Synthetic, Bacterial, and Protozoan RNA by Toll-like Receptor 8 Requires Coordinated Processing by RNase T2 and RNase 2

Human toll-like receptor 8 (TLR8) activation induces a potent T helper-1 (Th1) cell response critical for defense against intracellular pathogens, including protozoa. The receptor harbors two distinct binding sites, uridine and di- and/or trinucleotides, but the RNases upstream of TLR8 remain poorly characterized. We identified two endolysosomal endoribonucleases, RNase T2 and RNase 2, that act synergistically to release uridine from oligoribonucleotides. RNase T2 cleaves preferentially before, and RNase 2 after, uridines. Live bacteria, P. falciparum-infected red blood cells, purified pathogen RNA, and synthetic oligoribonucleotides all required RNase 2 and T2 processing to activate TLR8. Uridine supplementation restored RNA recognition in RNASE2^-/- or RNASET2^-/- but not RNASE2^-/-RNASET2^-/- cells. Primary immune cells from RNase T2-hypomorphic patients lacked a response to bacterial RNA but responded robustly to small-molecule TLR8 ligands. Our data identify an essential function of RNase T2 and RNase 2 upstream of TLR8 and provide insight into TLR8 activation.

An ancient evolutionary connection between Ribonuclease A and EndoU families

The ribonuclease A family of proteins is well studied from the biochemical and biophysical points of view, but its evolutionary origins are obscure, as no sequences homologous to this family have been reported outside of vertebrates. Recently, the spatial structure of the ribonuclease domain from a bacterial polymorphic toxin was shown to be closely similar to the structure of vertebrate ribonuclease A. The absence of sequence similarity between the two structures prompted a speculation of convergent evolution of bacterial and vertebrate ribonuclease A-like enzymes. We show that bacterial and homologous archaeal polymorphic toxin ribonucleases with a known or predicted ribonuclease A-like fold are distant homologs of the ribonucleases from the EndoU family, found in all domains of cellular life and in viruses. We also detected a homolog of vertebrate ribonucleases A in the transcriptome assembly of the sea urchin Mesocentrotus franciscanus These observations argue for the common ancestry of prokaryotic ribonuclease A-like and ubiquitous EndoU-like ribonucleases, and suggest a better-grounded scenario for the origin of animal ribonucleases A, which could have emerged in the deuterostome lineage, either by an extensive modification of a copy of an EndoU gene, or, more likely, by a horizontal acquisition of a prokaryotic immunity-mediating ribonuclease gene.