The discovery and characterization of two novel structural motifs on the carboxy-terminal domain of kinetoplastid RNA editing ligases.

Parasitic protozoans of the Trypanosoma and Leishmania species have a uniquely organized mitochondrial genome, the kinetoplast. Most kinetoplast-transcribed mRNAs are cryptic and encode multiple subunits for the electron transport chain following maturation through a uridine insertion/deletion process called RNA editing. This process is achieved through an enzyme cascade by an RNA editing catalytic complex (RECC), where the final ligation step is catalyzed by the kinetoplastid RNA editing ligases, KREL1 and KREL2. While the amino-terminal domain (NTD) of these proteins is highly conserved with other DNA ligases and mRNA capping enzymes, with five recognizable motifs, the functional role of their diverged carboxy-terminal domain (CTD) has remained elusive. In this manuscript, we assayed recombinant KREL1 in vitro to unveil critical residues from its CTD to be involved in protein-protein interaction and dsRNA ligation activity. Our data show that the α-helix (H)3 of KREL1 CTD interacts with the αH1 of its editosome protein partner KREPA2. Intriguingly, the OB-fold domain and the zinc fingers on KREPA2 do not appear to influence the RNA ligation activity of KREL1. Moreover, a specific KWKE motif on the αH4 of KREL1 CTD is found to be implicated in ligase auto-adenylylation analogous to motif VI in DNA ligases. In summary, we present in the KREL1 CTD a motif VI for auto-adenylylation and a KREPA2 binding motif for RECC integration.

Hammerhead ribozyme-based U-insertion and deletion RNA editing assays for multiplexing in HTS applications.

Untranslatable mitochondrial transcripts in kinetoplastids are decrypted post-transcriptionally through an RNA editing process that entails uridine insertion/deletion. This unique stepwise process is mediated by the editosome, a multiprotein complex that is a validated drug target of considerable interest in addressing the unmet medical needs for kinetoplastid diseases. With that objective, several in vitro RNA editing assays have been developed, albeit with limited success in discovering potent inhibitors. This manuscript describes the development of three hammerhead ribozyme (HHR) FRET reporter-based RNA editing assays for precleaved deletion, insertion, and ligation assays that bypass the rate-limiting endonucleolytic cleavage step, providing information on U-deletion, U-insertion, and ligation activities. These assays exhibit higher editing efficiencies in shorter incubation times while requiring significantly less purified editosome and 10,000-fold less ATP than the previously published full round of in vitro RNA editing assay. Moreover, modifications in the reporter ribozyme sequence enable the feasibility of multiplexing a ribozyme-based insertion/deletion editing (RIDE) assay that simultaneously surveils U-insertion and deletion editing suitable for HTS. These assays can be used to find novel chemical compounds with chemotherapeutic applications or as probes for studying the editosome machinery.

A Selective and ''Off-On'' Fluorescent Chemosensor Based on Fluorescein for Al3+: Synthesis, Characterization, Spectroscopy Analyses, and DFT Calculation.

An efficient fluorescent cation chemosensor based on fluorescein L4 was well prepared and identified with spectroscopy analyses. UV-vis and fluorescence measurements examined the analyte complexation of the L4 with various cations, demonstrating a clear tendency to Al3+ ion. In the Job plot study, a stoichiometry ratio of a complex between L4 and Al3+ ion was determined to be 1: 2 (L4: Al3+). A stoichiometry ratio of complex between L4 and Al3+ ion was determined to be 1: 2 (L4: Al3+) using the Job plot. The association constant (Ka) of the L4-Al3+ complex was found 2.8 × 107 M-2. The obtained limit of detection (LOD) value (1.37 × 10-6 M for Al3+) exhibited the considerable sensitivity of the chemosensor L4 to Al3+ ion. DFT/TD-DFT calculations have also been employed to support the binding mode and photophysical properties of the complexation of chemosensor L4 to Al3+ ion and also to investigate the enhancement of L4 fluorescence by Al3+ ion.

Structure-Function Analysis of RBP7910: An Editosome Z-Binding Protein in Trypanosomatids.

RNA editing, a unique post-transcriptional modification, is observed in trypanosomatid parasites as a crucial procedure for the maturation of mitochondrial mRNAs. The editosome protein complex, involving multiple protein components, plays a key role in this process. In Trypanosoma brucei, a putative Z-DNA binding protein known as RBP7910 is associated with the editosome. However, the specific Z-DNA/Z-RNA binding activity and the interacting interface of RBP7910 have yet to be determined. In this study, we conducted a comparative analysis of the binding behavior of RBP7910 with different potential ligands using microscale thermophoresis (MST). Additionally, we generated a 3D model of the protein, revealing potential Z-α and Z-ß nucleic acid-binding domains of RBP7910. RBP7910 belongs to the winged-helix-turn-helix (HTH) superfamily of proteins with an α1α2α3ß1ß2 topology. Finally, using docking techniques, potential interacting surface regions of RBP7910 with notable oligonucleotide ligands were identified. Our findings indicate that RBP7910 exhibits a notable affinity for (CG)n Z-DNA, both in single-stranded and double-stranded forms. Moreover, we observed a broader interacting interface across its Z-α domain when bound to Z-DNA/Z-RNA compared to when bound to non-Z-form nucleic acid ligands.

Sulfonated inhibitors of the RNA editing ligases validate the essential role of the MRP1/2 proteins in kinetoplastid RNA editing.

The RNA editing core complex (RECC) catalyzes mitochondrial U-insertion/deletion mRNA editing in trypanosomatid flagellates. Some naphthalene-based sulfonated compounds, such as C35 and MrB, competitively inhibit the auto-adenylylation activity of an essential RECC enzyme, kinetoplastid RNA editing ligase 1 (KREL1), required for the final step in editing. Previous studies revealed the ability of these compounds to interfere with the interaction between the editosome and its RNA substrates, consequently affecting all catalytic activities that comprise RNA editing. This observation implicates a critical function for the affected RNA binding proteins in RNA editing. In this study, using the inhibitory compounds, we analyzed the composition and editing activities of functional editosomes and identified the mitochondrial RNA binding proteins 1 and 2 (MRP1/2) as their preferred targets. While the MRP1/2 heterotetramer complex is known to bind guide RNA and promote annealing to its cognate pre-edited mRNA, its role in RNA editing remained enigmatic. We show that the compounds affect the association between the RECC and MRP1/2 heterotetramer. Furthermore, RECC purified post-treatment with these compounds exhibit compromised in vitro RNA editing activity that, remarkably, recovers upon the addition of recombinant MRP1/2 proteins. This work provides experimental evidence that the MRP1/2 heterotetramer is required for in vitro RNA editing activity and substantiates the hypothesized role of these proteins in presenting the RNA duplex to the catalytic complex in the initial steps of RNA editing.

Genetic and morphometric categorization of Taenia ovis from Sheep in Iran.

Little is known about the genetic and morphological characters of Taenia ovis. The purpose of the present study was to characterize sheep isolates of T. ovis using rostellar hook morphometry as well as mitochondrial genes sequence analysis. Ninety sheep specimens of Cysticercus ovis were collected from 18 slaughterhouses in Iran. The mean ± s.d. for total length of large and small hooks were 174.1 ± 6.4 and 116.7 ± 5.4 µm, respectively. CO1 and 12S rRNA sequence analysis showed 11 and nine haplotypes, respectively. The level of pairwise nucleotide variations between individual haplotypes of CO1 and 12S rRNA genes were 0.3-1.1 and 0.2-1.0%, respectively. Level of nucleotide variation in CO1 and 12S rRNA between T. ovis haplotypes from present study and eight other Taenia species was found to be 11.3-17.8 and 5.3-16.3%, respectively. Phylogenetic analysis clustered all T. ovis isolates into a single clade comprised of the all CO1 and 12S rRNA haplotypes. CO1 nucleotide difference between T. ovis ovis and T. asiatica was 13.6% that is lesser than the corresponding difference between T. ovis ovis and T. ovis krabbei, warranting the designation of two separate species as T. ovis and T. krabbei. Interclass correlation coefficients showed that there was no significant association between rostellar hook length variation and the variability of the mitochondrial genes.

circTAIL-seq, a targeted method for deep analysis of RNA 3' tails, reveals transcript-specific differences by multiple metrics.

Post-transcriptionally added RNA 3' nucleotide extensions, or tails, impose numerous regulatory effects on RNAs, including effects on RNA turnover and translation. However, efficient methods for in-depth tail profiling of a transcript of interest are still lacking, hindering available knowledge particularly of tail populations that are highly heterogeneous. Here, we developed a targeted approach, termed circTAIL-seq, to quantify both major and subtle differences of heterogeneous tail populations. As proof-of-principle, we show that circTAIL-seq quantifies the differences in tail qualities between two selected Trypanosoma brucei mitochondrial transcripts. The results demonstrate the power of the developed method in identification, discrimination, and quantification of different tail states that the population of one transcript can possess. We further show that circTAIL-seq can detect the tail characteristics for variants of transcripts that are not easily detectable by conventional approaches, such as degradation intermediates. Our findings are not only well supported by previous knowledge, but they also expand this knowledge and provide experimental evidence for previous hypotheses. In the future, this approach can be used to determine changes in tail qualities in response to environmental or internal stimuli, or upon silencing of genes of interest in mRNA-processing pathways. In summary, circTAIL-seq is an effective tool for comparing nonencoded RNA tails, especially when the tails are extremely variable or transcript of interest is low abundance.

Systematic discovery of structural elements governing stability of mammalian messenger RNAs.

Decoding post-transcriptional regulatory programs in RNA is a critical step towards the larger goal of developing predictive dynamical models of cellular behaviour. Despite recent efforts, the vast landscape of RNA regulatory elements remains largely uncharacterized. A long-standing obstacle is the contribution of local RNA secondary structure to the definition of interaction partners in a variety of regulatory contexts, including--but not limited to--transcript stability, alternative splicing and localization. There are many documented instances where the presence of a structural regulatory element dictates alternative splicing patterns (for example, human cardiac troponin T) or affects other aspects of RNA biology. Thus, a full characterization of post-transcriptional regulatory programs requires capturing information provided by both local secondary structures and the underlying sequence. Here we present a computational framework based on context-free grammars and mutual information that systematically explores the immense space of small structural elements and reveals motifs that are significantly informative of genome-wide measurements of RNA behaviour. By applying this framework to genome-wide human mRNA stability data, we reveal eight highly significant elements with substantial structural information, for the strongest of which we show a major role in global mRNA regulation. Through biochemistry, mass spectrometry and in vivo binding studies, we identified human HNRPA2B1 (heterogeneous nuclear ribonucleoprotein A2/B1, also known as HNRNPA2B1) as the key regulator that binds this element and stabilizes a large number of its target genes. We created a global post-transcriptional regulatory map based on the identity of the discovered linear and structural cis-regulatory elements, their regulatory interactions and their target pathways. This approach could also be used to reveal the structural elements that modulate other aspects of RNA behaviour.

Structural Studies of Trypanosoma brucei RNA Editing Ligases and Their Binding Partner Proteins.

To study the mechanism of ligating nicked RNA strands, we conducted molecular dynamics simulations of Trypanosoma brucei RNA editing ligases L1 and L2 complexed with double-stranded RNA (dsRNA) fragments. In each resulting model, a Mg(2+) ion coordinates the 5'-PO4 of the nicked nucleotide and the 3'-OH of the terminal nucleotide for a nucleophilic reaction consistent with the postulated step 3 chemistry of the ligation mechanism. Moreover, coordination of the 3'-OH to the Mg(2+) ion may lower its pKa, thereby rendering it a more effective nucleophile as an oxyanion. Thus, Mg(2+) may play a twofold role: bringing the reactants into the proximity of each other and activating the nucleophile. We also conducted solvated interaction energy calculations to explore whether ligation specificities can be correlated to ligase-dsRNA binding affinity changes. The calculated dsRNA binding affinities are stronger for both L1 and L2 when the terminal nucleotide is changed from cytosine to guanine, in line with their experimentally measured ligation specificities. Because the ligation mechanism is also influenced by interactions of the ligase with partner proteins from the editosome subcomplex, we also modeled the structure of the RNA-bound L2 in complex with the oligonucleotide binding (OB) domain of largest editosome interacting protein A1. The resulting L2-dsRNA-A1 model, which is consistent with mutagenesis and binding data recorded to date, provides the first atomic-level glimpse of plausible interactions around the RNA ligation site in the presence of an OB domain presented in-trans to a nucleic acid ligase.

Global identification of conserved post-transcriptional regulatory programs in trypanosomatids.

While regulatory programs are extensively studied at the level of transcription, elements that are involved in regulation of post-transcriptional processes are largely unknown, and methods for systematic identification of these elements are in early stages. Here, using a novel computational framework, we have integrated sequence information with several functional genomics data sets to characterize conserved regulatory programs of trypanosomatids, a group of eukaryotes that almost entirely rely on post-transcriptional processes for regulation of mRNA abundance. This analysis revealed a complex network of linear and structural RNA elements that potentially govern mRNA abundance across different life stages and environmental conditions. Furthermore, we show that the conserved regulatory network that we have identified is responsive to chemical perturbation of several biological functions in trypanosomatids. We have further characterized one of the most abundant regulatory RNA elements that we discovered, an AU-rich element (ARE) that can be found in 3' untranslated region of many trypanosomatid genes. Using bioinformatics approaches as well as in vitro and in vivo experiments, we have identified three ELAV-like homologs, including the developmentally critical protein TbRBP6, which regulate abundance of a large number of trypanosomatid ARE-containing transcripts. Together, these studies lay out a roadmap for characterization of mechanisms that modulate development and metabolic pathways in trypanosomatids.

Genetic variability of Taenia saginata inferred from mitochondrial DNA sequences.

Taenia saginata is an important tapeworm, infecting humans in many parts of the world. The present study was undertaken to identify inter- and intraspecific variation of T. saginata isolated from cattle in different parts of Iran using two mitochondrial CO1 and 12S rRNA genes. Up to 105 bovine specimens of T. saginata were collected from 20 slaughterhouses in three provinces of Iran. DNA were extracted from the metacestode Cysticercus bovis. After PCR amplification, sequencing of CO1 and 12S rRNA genes were carried out and two phylogenetic analyses of the sequence data were generated by Bayesian inference on CO1 and 12S rRNA sequences. Sequence analyses of CO1 and 12S rRNA genes showed 11 and 29 representative profiles respectively. The level of pairwise nucleotide variation between individual haplotypes of CO1 gene was 0.3-2.4% while the overall nucleotide variation among all 11 haplotypes was 4.6%. For 12S rRNA sequence data, level of pairwise nucleotide variation was 0.2-2.5% and the overall nucleotide variation was determined as 5.8% among 29 haplotypes of 12S rRNA gene. Considerable genetic diversity was found in both mitochondrial genes particularly in 12S rRNA gene.

Prioritization of Trypanosoma brucei editosome protein interactions interfaces at residue resolution through proteome-scale network analysis.

BACKGROUND: Trypanosoma brucei is the causative agent for trypanosomiasis in humans and livestock, which presents a growing challenge due to drug resistance. While identifying novel drug targets is vital, the process is delayed due to a lack of functional information on many of the pathogen's proteins. Accordingly, this paper presents a computational framework for prioritizing drug targets within the editosome, a vital molecular machinery responsible for mitochondrial RNA processing in T. brucei. Importantly, this framework may eliminate the need for prior gene or protein characterization, potentially accelerating drug discovery efforts. RESULTS: By integrating protein-protein interaction (PPI) network analysis, PPI structural modeling, and residue interaction network (RIN) analysis, we quantitatively ranked and identified top hub editosome proteins, their key interaction interfaces, and hotspot residues. Our findings were cross-validated and further prioritized by incorporating them into gene set analysis and differential expression analysis of existing quantitative proteomics data across various life stages of T. brucei. In doing so, we highlighted PPIs such as KREL2-KREPA1, RESC2-RESC1, RESC12A-RESC13, and RESC10-RESC6 as top candidates for further investigation. This includes examining their interfaces and hotspot residues, which could guide drug candidate selection and functional studies. CONCLUSION: RNA editing offers promise for target-based drug discovery, particularly with proteins and interfaces that play central roles in the pathogen's life cycle. This study introduces an integrative drug target identification workflow combining information from the PPI network, PPI 3D structure, and reside-level information of their interface which can be applicable to diverse pathogens. In the case of T. brucei, via this pipeline, the present study suggested potential drug targets with residue-resolution from RNA editing machinery. However, experimental validation is needed to fully realize its potential in advancing urgently needed antiparasitic drug development.

Functional domain annotation by structural similarity.

Traditional automated in silico functional annotation uses tools like Pfam that rely on sequence similarities for domain annotation. However, structural conservation often exceeds sequence conservation, suggesting an untapped potential for improved annotation through structural similarity. This approach was previously overlooked before the AlphaFold2 introduction due to the need for more high-quality protein structures. Leveraging structural information especially holds significant promise to enhance accurate annotation in diverse proteins across phylogenetic distances. In our study, we evaluated the feasibility of annotating Pfam domains based on structural similarity. To this end, we created a database from segmented full-length protein structures at their domain boundaries, representing the structure of Pfam seeds. We used Trypanosoma brucei, a phylogenetically distant protozoan parasite as our model organism. Its structome was aligned with our database using Foldseek, the ultra-fast structural alignment tool, and the top non-overlapping hits were annotated as domains. Our method identified over 400 new domains in the T. brucei proteome, surpassing the benchmark set by sequence-based tools, Pfam and Pfam-N, with some predictions validated manually. We have also addressed limitations and suggested avenues for further enhancing structure-based domain annotation.

From a bimodal to a multi-stage view on trypanosomes' differential RNA editing.

Significant variations in the abundance of mitochondrial RNA processing proteins and their target RNAs across trypanosome life stages present an opportunity to explore the regulatory mechanisms that drive these changes. Utilizing omics approaches can uncover unconventional targets, aiding our understanding of the parasites' adaptation and enabling targeted interventions for differentiation.

Pilot-Scale Screening of Clinically Approved Drugs to Identify Uridine Insertion/Deletion RNA Editing Inhibitors in Trypanosoma brucei.

RNA editing pathway is a validated target in kinetoplastid parasites (Trypanosoma brucei, Trypanosoma cruzi, and Leishmania spp.) that cause severe diseases in humans and livestock. An essential large protein complex, the editosome, mediates uridine insertion and deletion in RNA editing through a stepwise process. This study details the discovery of editosome inhibitors by screening a library of widely used human drugs using our previously developed in vitro biochemical Ribozyme Insertion Deletion Editing (RIDE) assay. Subsequent studies on the mode of action of the identified hits and hit expansion efforts unveiled compounds that interfere with RNA-editosome interactions and novel ligase inhibitors with IC50 values in the low micromolar range. Docking studies on the ligase demonstrated similar binding characteristics for ATP and our novel epigallocatechin gallate inhibitor. The inhibitors demonstrated potent trypanocidal activity and are promising candidates for drug repurposing due to their lack of cytotoxic effects. Further studies are necessary to validate these targets using more definitive gene-editing techniques and to enhance the safety profile.

High-throughput screening of compounds targeting RNA editing in Trypanosoma brucei: Novel molecular scaffolds with broad trypanocidal effects.

Mitochondrial uridine insertion/deletion RNA editing, catalyzed by a multiprotein complex (editosome), is essential for gene expression in trypanosomes and Leishmania parasites. As this process is absent in the human host, a drug targeting this mechanism promises high selectivity and reduced toxicity. Here, we successfully miniaturized our FRET-based full-round RNA editing assay, which replicates the complete RNA editing process, adapting it into a 1536-well format. Leveraging this assay, we screened over 100,000 compounds against purified editosomes derived from Trypanosoma brucei, identifying seven confirmed primary hits. We sourced and evaluated various analogs to enhance the inhibitory and parasiticidal effects of these primary hits. In combination with secondary assays, our compounds marked inhibition of essential catalytic activities, including the RNA editing ligase and interactions of editosome proteins. Although the primary hits did not exhibit any growth inhibitory effect on parasites, we describe eight analog compounds capable of effectively killing T. brucei and/or Leishmania donovani parasites within a low micromolar concentration. Whether parasite killing is - at least in part - due to inhibition of RNA editing in vivo remains to be assessed. Our findings introduce novel molecular scaffolds with the potential for broad antitrypanosomal effects.

MetaboAnalystR 4.0: a unified LC-MS workflow for global metabolomics.

The wide applications of liquid chromatography - mass spectrometry (LC-MS) in untargeted metabolomics demand an easy-to-use, comprehensive computational workflow to support efficient and reproducible data analysis. However, current tools were primarily developed to perform specific tasks in LC-MS based metabolomics data analysis. Here we introduce MetaboAnalystR 4.0 as a streamlined pipeline covering raw spectra processing, compound identification, statistical analysis, and functional interpretation. The key features of MetaboAnalystR 4.0 includes an auto-optimized feature detection and quantification algorithm for LC-MS1 spectra processing, efficient MS2 spectra deconvolution and compound identification for data-dependent or data-independent acquisition, and more accurate functional interpretation through integrated spectral annotation. Comprehensive validation studies using LC-MS1 and MS2 spectra obtained from standards mixtures, dilution series and clinical metabolomics samples have shown its excellent performance across a wide range of common tasks such as peak picking, spectral deconvolution, and compound identification with good computing efficiency. Together with its existing statistical analysis utilities, MetaboAnalystR 4.0 represents a significant step toward a unified, end-to-end workflow for LC-MS based global metabolomics in the open-source R environment.

Circular mitochondrial-encoded mRNAs are a distinct subpopulation of mitochondrial mRNA in Trypanosoma brucei.

Since the first identification of circular RNA (circRNA) in viral-like systems, reports of circRNAs and their functions in various organisms, cell types, and organelles have greatly expanded. Here, we report the first evidence of circular mRNA in the mitochondrion of the eukaryotic parasite, Trypanosoma brucei . While using a circular RT-PCR technique developed to sequence mRNA tails of mitochondrial transcripts, we found that some mRNAs are circularized without an in vitro circularization step normally required to produce PCR products. Starting from total in vitro circularized RNA and in vivo circRNA, we high-throughput sequenced three transcripts from the 3' end of the coding region, through the 3' tail, to the 5' start of the coding region. We found that fewer reads in the circRNA libraries contained tails than in the total RNA libraries. When tails were present on circRNAs, they were shorter and less adenine-rich than the total population of RNA tails of the same transcript. Additionally, using hidden Markov modelling we determined that enzymatic activity during tail addition is different for circRNAs than for total RNA. Lastly, circRNA UTRs tended to be shorter and more variable than those of the same transcript sequenced from total RNA. We propose a revised model of Trypanosome mitochondrial tail addition, in which a fraction of mRNAs is circularized prior to the addition of adenine-rich tails and may act as a new regulatory molecule or in a degradation pathway.

Circular mitochondrial-encoded mRNAs are a distinct subpopulation of mitochondrial mRNA in Trypanosoma brucei.

Since the first identification of circular RNA (circRNA) in viral-like systems, reports of circRNAs and their functions in various organisms, cell types, and organelles have greatly expanded. Here, we report the first evidence, to our knowledge, of circular mRNA in the mitochondrion of the eukaryotic parasite, Trypanosoma brucei. While using a circular RT-PCR technique developed to sequence mRNA tails of mitochondrial transcripts, we found that some mRNAs are circularized without an in vitro circularization step normally required to produce PCR products. Starting from total in vitro circularized RNA and in vivo circRNA, we high-throughput sequenced three transcripts from the 3' end of the coding region, through the 3' tail, to the 5' start of the coding region. We found that fewer reads in the circRNA libraries contained tails than in the total RNA libraries. When tails were present on circRNAs, they were shorter and less adenine-rich than the total population of RNA tails of the same transcript. Additionally, using hidden Markov modelling we determined that enzymatic activity during tail addition is different for circRNAs than for total RNA. Lastly, circRNA UTRs tended to be shorter and more variable than those of the same transcript sequenced from total RNA. We propose a revised model of Trypanosome mitochondrial tail addition, in which a fraction of mRNAs is circularized prior to the addition of adenine-rich tails and may act as a new regulatory molecule or in a degradation pathway.