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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Lexis and Grammar of Mitochondrial RNA Processing in Trypanosomes.

Trypanosoma brucei spp. cause African human and animal trypanosomiasis, a burden on health and economy in Africa. These hemoflagellates are distinguished by a kinetoplast nucleoid containing mitochondrial DNAs of two kinds: maxicircles encoding ribosomal RNAs (rRNAs) and proteins and minicircles bearing guide RNAs (gRNAs) for mRNA editing. All RNAs are produced by a phage-type RNA polymerase as 3' extended precursors, which undergo exonucleolytic trimming. Most pre-mRNAs proceed through 3' adenylation, uridine insertion/deletion editing, and 3' A/U-tailing. The rRNAs and gRNAs are 3' uridylated. Historically, RNA editing has attracted major research effort, and recently essential pre- and postediting processing events have been discovered. Here, we classify the key players that transform primary transcripts into mature molecules and regulate their function and turnover.

The RNA binding activity of the first identified trypanosome protein with Z-DNA-binding domains.

RNA-binding proteins play a particularly important role in regulating gene expression in trypanosomes. A map of the network of protein complexes in Trypanosoma brucei uncovered an essential protein (Tb927.10.7910) that is postulated to be an RNA-binding protein implicated in the regulation of the mitochondrial post-transcriptional gene regulatory network by its association with proteins that participate in a multi-protein RNA editing complex. However, the mechanism by which this protein interacts with its multiple target transcripts remained unknown. Using sensitive database searches and experimental data, we identify Z-DNA-binding domains in T. brucei in the N- and C-terminal regions of Tb927.10.7910. RNA-binding studies of the wild-type protein, now referred to as RBP7910 (RNA binding protein 7910), and site-directed mutagenesis of residues important for the Z-DNA binding domains show that it preferentially interacts with RNA molecules containing poly(U) and poly(AU)-rich sequences. The interaction of RBP7910 with these regions may be involved in regulation of RNA editing of mitochondrial transcripts.

Gene Function Discovery for Kinetoplastid Pathogens.

We propose to integrate the existing and new experimental data with computational tools to model interaction networks for the most prominent kinetoplastid pathogens. These interaction networks will vastly expand the functional annotation of the kinetoplastid genomes, which in turn are critical for identifying new routes of disease intervention.

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.

Tail characteristics of Trypanosoma brucei mitochondrial transcripts are developmentally altered in a transcript-specific manner.

The intricate life cycle of Trypanosoma brucei requires extensive regulation of gene expression levels of the mtRNAs for adaptation. Post-transcriptional gene regulatory programs, including unencoded mtRNA 3' tail additions, potentially play major roles in this adaptation process. Intriguingly, T. brucei mitochondrial transcripts possess two distinct unencoded 3' tails, each with a differing functional role; i.e., while one type is implicated in RNA stability (in-tails), the other type appears associated with translation (ex-tails). We examined the degree to which tail characteristics differ among cytochrome c oxidase subunits I and III (CO1 and CO3), and NADH dehydrogenase subunit 1 (ND1) transcripts, and to what extent these characteristics differ developmentally. We found that CO1, CO3 and ND1 transcripts possess longer in-tails in the mammalian life stage. By mathematically modelling states of in-tail and ex-tail addition, we determined that the typical length at which an in-tail is extended to become an ex-tail differs by transcript and, in the case of ND1, by life stage. To the best of our knowledge, we provide the first evidence that developmental differences exist in tail length distributions of mtRNAs, underscoring the potential involvement of in-tail and ex-tail populations in mitochondrial post-transcriptional regulation mechanisms.

TrypsNetDB: An integrated framework for the functional characterization of trypanosomatid proteins.

Trypanosomatid parasites cause serious infections in humans and production losses in livestock. Due to the high divergence from other eukaryotes, such as humans and model organisms, the functional roles of many trypanosomatid proteins cannot be predicted by homology-based methods, rendering a significant portion of their proteins as uncharacterized. Recent technological advances have led to the availability of multiple systematic and genome-wide datasets on trypanosomatid parasites that are informative regarding the biological role(s) of their proteins. Here, we report TrypsNetDB (http://trypsNetDB.org), a web-based resource for the functional annotation of 16 different species/strains of trypanosomatid parasites. The database not only visualizes the network context of the queried protein(s) in an intuitive way but also examines the response of the represented network in more than 50 different biological contexts and its enrichment for various biological terms and pathways, protein sequence signatures, and potential RNA regulatory elements. The interactome core of the database, as of Jan 23, 2017, contains 101,187 interactions among 13,395 trypanosomatid proteins inferred from 97 genome-wide and focused studies on the interactome of these organisms.

The interaction of a Trypanosoma brucei KH-domain protein with a ribonuclease is implicated in ribosome processing.

Ribosomal RNA maturation is best understood in yeast. While substantial efforts have been made to explore parts of these essential pathways in animals, the similarities and uniquenesses of rRNA maturation factors in non-Opisthokonts remain largely unexplored. Eukaryotic ribosome synthesis requires the coordinated activities of hundreds of Assembly Factors (AFs) that transiently associate with pre-ribosomes, many of which are essential. Pno1 and Nob1 are two of six AFs that are required for the cytoplasmic maturation of the 20S pre-rRNA to 18S rRNA in yeast where it has been almost exclusively analyzed. Specifically, Nob1 ribonucleolytic activity generates the mature 3'-end of 18S rRNA. We identified putative Pno1 and Nob1 homologues in the protist Trypanosoma brucei, named TbPNO1 and TbNOB1, and set out to explore their rRNA maturation role further as they are both essential for normal growth. TbPNO1 is a nuclear protein with limited cytosolic localization relative to its yeast homologue. Like in yeast, it interacts directly with TbNOB1, with indications of associations with a larger AF-containing complex. Interestingly, in the absence of TbPNO1, TbNOB1 exhibits non-specific degradation activity on RNA substrates, and its cleavage activity becomes specific only in the presence of TbPNO1, suggesting that TbPNO1-TbNOB1 interaction is essential for regulation and site-specificity of TbNOB1 activity. These results highlight a conserved role of the TbPNO1-TbNOB1 complex in 18S rRNA maturation across eukaryotes; yet reveal a novel role of their interaction in regulation of TbNOB1 enzymatic activity.