Bacteria within the gut microbiota possess the ability to metabolize a wide array of human drugs, foods, and toxins, but the responsible enzymes for these chemical events remain largely uncharacterized due to the time-consuming nature of current experimental approaches. Attempts have been made in the past to computationally predict which bacterial species and enzymes are responsible for chemical transformations in the gut environment, but with low accuracy due to minimal chemical representation and sequence similarity search schemes. Here, we present an in silico approach that employs chemical and protein Similarity algorithms that Identify MicrobioMe Enzymatic Reactions (SIMMER). We show that SIMMER accurately predicts the responsible species and enzymes for a queried reaction, unlike previous methods. We demonstrate SIMMER use cases in the context of drug metabolism by predicting previously uncharacterized enzymes for 88 drug transformations known to occur in the human gut. We validate these predictions on external datasets and provide an in vitro validation of SIMMER's predictions for metabolism of methotrexate, an anti-arthritic drug. After demonstrating its utility and accuracy, we made SIMMER available as both a command-line and web tool, with flexible input and output options for determining chemical transformations within the human gut. We present SIMMER as a computational addition to the microbiome researcher's toolbox, enabling them to make informed hypotheses before embarking on the lengthy laboratory experiments required to characterize novel bacterial enzymes that can alter human ingested compounds.
Data generated and analyzed during this study are provided in Figures 2-10 source data files, Table 1 source data file, supplemental files, and at https://github.com/aebustion/SIMMER. Accession numbers of previously published datasets are provided in the Materials and Methods section. SIMMER code can either be run at the SIMMER website (https://simmer.pollard.gladstone.org/) or downloaded directly from the above-linked GitHub.
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
© 2023, Bustion et al.
This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Copper is an essential enzyme cofactor in bacteria, but excess copper is highly toxic. Bacteria can cope with copper stress by increasing copper resistance and initiating chemorepellent response. However, it remains unclear how bacteria coordinate chemotaxis and resistance to copper. By screening proteins that interacted with the chemotaxis kinase CheA, we identified a copper-binding repressor CsoR that interacted with CheA in Pseudomonas putida. CsoR interacted with the HPT (P1), Dimer (P3), and HATPase_c (P4) domains of CheA and inhibited CheA autophosphorylation, resulting in decreased chemotaxis. The copper-binding of CsoR weakened its interaction with CheA, which relieved the inhibition of chemotaxis by CsoR. In addition, CsoR bound to the promoter of copper-resistance genes to inhibit gene expression, and copper-binding released CsoR from the promoter, leading to increased gene expression and copper resistance. P. putida cells exhibited a chemorepellent response to copper in a CheA-dependent manner, and CsoR inhibited the chemorepellent response to copper. Besides, the CheA-CsoR interaction also existed in proteins from several other bacterial species. Our results revealed a mechanism by which bacteria coordinately regulated chemotaxis and resistance to copper by CsoR.
5-Methylcytosine (m5C) is one of the posttranscriptional modifications in mRNA and is involved in the pathogenesis of various diseases. However, the capacity of existing assays for accurately and comprehensively transcriptome-wide m5C mapping still needs improvement. Here, we develop a detection method named DRAM (deaminase and reader protein assisted RNA methylation analysis), in which deaminases (APOBEC1 and TadA-8e) are fused with m5C reader proteins (ALYREF and YBX1) to identify the m5C sites through deamination events neighboring the methylation sites. This antibody-free and bisulfite-free approach provides transcriptome-wide editing regions which are highly overlapped with the publicly available bisulfite-sequencing (BS-seq) datasets and allows for a more stable and comprehensive identification of the m5C loci. In addition, DRAM system even supports ultralow input RNA (10 ng). We anticipate that the DRAM system could pave the way for uncovering further biological functions of m5C modifications.