Microbial transformation of dietary xenobiotics shapes gut microbiome composition.

Diet is a major determinant of gut microbiome composition, and variation in diet-microbiome interactions may contribute to variation in their health consequences. To mechanistically understand these relationships, here we map interactions between ∼150 small-molecule dietary xenobiotics and the gut microbiome, including the impacts of these compounds on community composition, the metabolic activities of human gut microbes on dietary xenobiotics, and interindividual variation in these traits. Microbial metabolism can toxify and detoxify these compounds, producing emergent interactions that explain community-specific remodeling by dietary xenobiotics. We identify the gene and enzyme responsible for detoxification of one such dietary xenobiotic, resveratrol, and demonstrate that this enzyme contributes to interindividual variation in community remodeling by resveratrol. Together, these results systematically map interactions between dietary xenobiotics and the gut microbiome and connect toxification and detoxification to interpersonal differences in microbiome response to diet.

Structural and Mechanistic Principles of ABC Transporters.

ATP-binding cassette (ABC) transporters constitute one of the largest and most ancient protein superfamilies found in all living organisms. They function as molecular machines by coupling ATP binding, hydrolysis, and phosphate release to translocation of diverse substrates across membranes. The substrates range from vitamins, steroids, lipids, and ions to peptides, proteins, polysaccharides, and xenobiotics. ABC transporters undergo substantial conformational changes during substrate translocation. A comprehensive understanding of their inner workings thus requires linking these structural rearrangements to the different functional state transitions. Recent advances in single-particle cryogenic electron microscopy have not only delivered crucial information on the architecture of several medically relevant ABC transporters and their supramolecular assemblies, including the ATP-sensitive potassium channel and the peptide-loading complex, but also made it possible to explore the entire conformational space of these nanomachines under turnover conditions and thereby gain detailed mechanistic insights into their mode of action.

A Circadian Clock in the Blood-Brain Barrier Regulates Xenobiotic Efflux.

Endogenous circadian rhythms are thought to modulate responses to external factors, but mechanisms that confer time-of-day differences in organismal responses to environmental insults/therapeutic treatments are poorly understood. Using a xenobiotic, we find that permeability of the Drosophila "blood"-brain barrier (BBB) is higher at night. The permeability rhythm is driven by circadian regulation of efflux and depends on a molecular clock in the perineurial glia of the BBB, although efflux transporters are restricted to subperineurial glia (SPG). We show that transmission of circadian signals across the layers requires cyclically expressed gap junctions. Specifically, during nighttime, gap junctions reduce intracellular magnesium ([Mg2+]i), a positive regulator of efflux, in SPG. Consistent with lower nighttime efflux, nighttime administration of the anti-epileptic phenytoin is more effective at treating a Drosophila seizure model. These findings identify a novel mechanism of circadian regulation and have therapeutic implications for drugs targeted to the central nervous system.

Boron catalysis in a designer enzyme.

Enzymes play an increasingly important role in improving the benignity and efficiency of chemical production, yet the diversity of their applications lags heavily behind chemical catalysts as a result of the relatively narrow range of reaction mechanisms of enzymes. The creation of enzymes containing non-biological functionalities facilitates reaction mechanisms outside nature's canon and paves the way towards fully programmable biocatalysis1-3. Here we present a completely genetically encoded boronic-acid-containing designer enzyme with organocatalytic reactivity not achievable with natural or engineered biocatalysts4,5. This boron enzyme catalyses the kinetic resolution of hydroxyketones by oxime formation, in which crucial interactions with the protein scaffold assist in the catalysis. A directed evolution campaign led to a variant with natural-enzyme-like enantioselectivities for several different substrates. The unique activation mode of the boron enzyme was confirmed using X-ray crystallography, high-resolution mass spectrometry (HRMS) and 11B NMR spectroscopy. Our study demonstrates that genetic-code expansion can be used to create evolvable enantioselective enzymes that rely on xenobiotic catalytic moieties such as boronic acids and access reaction mechanisms not reachable through catalytic promiscuity of natural or engineered enzymes.

Inactivation of conserved C. elegans genes engages pathogen- and xenobiotic-associated defenses.

The nematode C. elegans is attracted to nutritious bacteria and is repelled by pathogens and toxins. Here we show that RNAi and toxin-mediated disruption of core cellular activities, including translation, respiration, and protein turnover, stimulate behavioral avoidance of normally attractive bacteria. RNAi of these and other essential processes induces expression of detoxification and innate immune effectors, even in the absence of toxins or pathogens. Disruption of core processes in non-neuronal tissues was sufficient to stimulate aversion behavior, revealing a neuroendocrine axis of control that additionally required serotonergic and Jnk kinase signaling pathways. We propose that surveillance pathways overseeing core cellular activities allow animals to detect invading pathogens that deploy toxins and virulence factors to undermine vital host functions. Variation in cellular surveillance and endocrine pathways controlling behavior, detoxification, and immunity selected by past toxin or microbial interactions could underlie aberrant responses to foods, medicines, and microbes.

Flavin-containing monooxygenase (FMO): Beyond xenobiotics.

Flavin-containing monooxygenases (FMOs), traditionally known for detoxifying xenobiotics, are now recognized for their involvement in endogenous metabolism. We recently discovered that an isoform of FMO, fmo-2 in Caenorhabditis elegans, alters endogenous metabolism to impact longevity and stress tolerance. Increased expression of fmo-2 in C. elegans modifies the flux through the key pathway known as One Carbon Metabolism (OCM). This modified flux results in a decrease in the ratio of S-adenosyl-methionine (SAM) to S-adenosyl-homocysteine (SAH), consequently diminishing methylation capacity. Here we discuss how FMO-2-mediated formate production during tryptophan metabolism may serve as a trigger for changing the flux in OCM. We suggest formate bridges tryptophan and OCM, altering metabolic flux away from methylation during fmo-2 overexpression. Additionally, we highlight how these metabolic results intersect with the mTOR and AMPK pathways, in addition to mitochondrial metabolism. In conclusion, the goal of this essay is to bring attention to the central role of FMO enzymes but lack of understanding of their mechanisms. We justify a call for a deeper understanding of FMO enzyme's role in metabolic rewiring through tryptophan/formate or other yet unidentified substrates. Additionally, we emphasize the identification of novel drugs and microbes to induce FMO activity and extend lifespan.

XRE transcription factors conserved in Caulobacter and φCbK modulate adhesin development and phage production.

The xenobiotic response element (XRE) family of transcription factors (TFs), which are commonly encoded by bacteria and bacteriophage, regulate diverse features of bacterial cell physiology and impact phage infection dynamics. Through a pangenome analysis of Caulobacter species isolated from soil and aquatic ecosystems, we uncovered an apparent radiation of a paralogous XRE TF gene cluster, several of which have established functions in the regulation of holdfast adhesin development and biofilm formation in C. crescentus. We further discovered related XRE TFs throughout the class Alphaproteobacteria and its phages, including the φCbK Caulophage, suggesting that members of this cluster impact host-phage interactions. Here we show that a closely related group of XRE transcription factors encoded by both C. crescentus and φCbK can physically interact and function to control the transcription of a common gene set, influencing processes including holdfast development and the production of φCbK virions. The φCbK-encoded XRE paralog, tgrL, is highly expressed at the earliest stages of infection and can directly inhibit transcription of host genes including hfiA, a potent holdfast inhibitor, and gafYZ, an activator of prophage-like gene transfer agents (GTAs). XRE proteins encoded from the C. crescentus chromosome also directly repress gafYZ transcription, revealing a functionally redundant set of host regulators that may protect against spurious production of GTA particles and inadvertent cell lysis. Deleting the C. crescentus XRE transcription factors reduced φCbK burst size, while overexpressing these host genes or φCbK tgrL rescued this burst defect. We conclude that this XRE TF gene cluster, shared by C. crescentus and φCbK, plays an important role in adhesion regulation under phage-free conditions, and influences host-phage dynamics during infection.

Aryl Hydrocarbon Receptor and Its Diverse Ligands and Functions: An Exposome Receptor.

The aryl hydrocarbon receptor (AhR) is a transcriptional factor that regulates multiple functions following its activation by a variety of ligands, including xenobiotics, natural products, microbiome metabolites, and endogenous molecules. Because of this diversity, the AhR constitutes an exposome receptor. One of its main functions is to regulate several lines of defense against chemical insults and bacterial infections. Indeed, in addition to its well-established detoxication function, it has several functions at physiological barriers, and it plays a critical role in immunomodulation. The AhR is also involved in the development of several organs and their homeostatic maintenance. Its activity depends on the type of ligand and on the time frame of the receptor activation, which can be either sustained or transient, leading in some cases to opposite modes of regulations as illustrated in the regulation of different cancer pathways. The development of selective modulators and their pharmacological characterization are important areas of research.

Massively parallel characterization of CYP2C9 variant enzyme activity and abundance.

CYP2C9 encodes a cytochrome P450 enzyme responsible for metabolizing up to 15% of small molecule drugs, and CYP2C9 variants can alter the safety and efficacy of these therapeutics. In particular, the anti-coagulant warfarin is prescribed to over 15 million people annually and polymorphisms in CYP2C9 can affect individual drug response and lead to an increased risk of hemorrhage. We developed click-seq, a pooled yeast-based activity assay, to test thousands of variants. Using click-seq, we measured the activity of 6,142 missense variants in yeast. We also measured the steady-state cellular abundance of 6,370 missense variants in a human cell line by using variant abundance by massively parallel sequencing (VAMP-seq). These data revealed that almost two-thirds of CYP2C9 variants showed decreased activity and that protein abundance accounted for half of the variation in CYP2C9 function. We also measured activity scores for 319 previously unannotated human variants, many of which may have clinical relevance.

Detection and Imaging of Exposure-Related Metabolites and Xenobiotics in Hard Tissues by Laser Sampling and Mass Spectrometry.

The utility of two novel laser-based methods, laser ablation electrospray ionization (LAESI) and laser desorption ionization (LDI) from silicon nanopost array (NAPA), is explored via local analysis and mass spectrometry imaging (MSI) of hard tissues (tooth and hair) for the detection and mapping of organic components. Complex mass spectra are recorded in local analysis mode from tooth dentin and scalp hair samples. Nicotine and its metabolites (cotinine, hydroxycotinine, norcotinine, and nicotine) are detected by LAESI-MS in the teeth of rats exposed to tobacco smoke. The intensities of the detected metabolite peaks are proportional to the degree of exposure. Incorporating ion mobility separation in the LAESI-MS analysis of scalp hair enables the detection of cotinine in smoker hair along with other common molecular species, including endogenous steroid hormones and some lipids. Single hair strands are imaged by MALDI-MSI and NAPA-LDI-MSI to explore longitudinal variations in the level of small molecules. Comparing spectra integrated from NAPA-LDI-MSI and MALDI-MSI images reveals that the two techniques provide complementary information. There were 105 and 82 sample-related peaks for MALDI and NAPA, respectively, with an overlap of only 16 peaks, indicating a high degree of complementarity. Enhanced molecular coverage and spatial resolution offered by LAESI-MS and NAPA-LDI-MSI can reveal the distributions of known and potential biomarkers in hard tissues, facilitating exposome research.

MetSim: Integrated Programmatic Access and Pathway Management for Xenobiotic Metabolism Simulators.

Xenobiotic metabolism is a key consideration in evaluating the hazards and risks posed by environmental chemicals. A number of software tools exist that are capable of simulating metabolites, but each reports its predictions in a different format and with varying levels of detail. This makes comparing the performance and coverage of the tools a practical challenge. To address this shortcoming, we developed a metabolic simulation framework called MetSim, which comprises three main components. A graph-based schema was developed to allow metabolism information to be harmonized. The schema was implemented in MongoDB to store and retrieve metabolic graphs for subsequent analysis. MetSim currently includes an application programming interface for four metabolic simulators: BioTransformer, the OECD Toolbox, EPA's chemical transformation simulator (CTS), and tissue metabolism simulator (TIMES). Lastly, MetSim provides functions to help evaluate simulator performance for specific data sets. In this study, a set of 112 drugs with 432 reported metabolites were compiled, and predictions were made using the 4 simulators. Fifty-nine of the 112 drugs were taken from the Small Molecule Pathway Database, with the remainder sourced from the literature. The human models within BioTransformer and CTS (Phase I only) and the rat models within TIMES and the OECD Toolbox (Phase I only) were used to make predictions for the chemicals in the data set. The recall and precision (recall, precision) ranked in order of highest recall for each individual tool were CTS (0.54, 0.017), BioTransformer (0.50, 0.008), Toolbox in vitro (0.40, 0.144), TIMES in vivo (0.40, 0.133), Toolbox in vivo (0.40, 0.118), and TIMES in vitro (0.39, 0.128). Combining all of the model predictions together increased the overall recall (0.73, 0.008). MetSim enabled insights into the performance and coverage of in silico metabolic simulators to be more efficiently derived, which in turn should aid future efforts to evaluate other data sets.

Interaction with B-type lamin reveals the function of Drosophila Keap1 xenobiotic response factor in nuclear architecture.

BACKGROUND: The Keap1-Nrf2 pathway serves as a central regulator that mediates transcriptional responses to xenobiotic and oxidative stimuli. Recent studies have shown that Keap1 and Nrf2 can regulate transcripts beyond antioxidant and detoxifying genes, yet the underlying mechanisms remain unclear. Our research has uncovered that Drosophila Keap1 (dKeap1) and Nrf2 (CncC) proteins can control high-order chromatin structure, including heterochromatin. METHODS AND RESULTS: In this study, we identified the molecular interaction between dKeap1 and lamin Dm0, the Drosophila B-type lamin responsible for the architecture of nuclear lamina and chromatin. Ectopic expression of dKeap1 led to an ectopic localization of lamin to the intra-nuclear area, corelated with the spreading of the heterochromatin marker H3K9me2 into euchromatin regions. Additionally, mis-regulated dKeap1 disrupted the morphology of the nuclear lamina. Knocking down of dKeap1 partially rescued the lethality induced by lamin overexpression, suggesting their genetic interaction during development. CONCLUSIONS: The discovered dKeap1-lamin interaction suggests a novel role for the Keap1 oxidative/xenobiotic response factor in regulating chromatin architecture.

A serum-induced gene signature in hepatocytes is associated with pediatric nonalcoholic fatty liver disease.

OBJECTIVE: Pediatric nonalcoholic fatty liver disease (NAFLD) is a growing problem, but its underlying mechanisms are poorly understood. We used transcriptomic reporter cell assays to investigate differences in transcriptional signatures induced in hepatocyte reporter cells by the sera of children with and without NAFLD. METHODS: We studied serum samples from 45 children with NAFLD and 28 children without NAFLD. The sera were used to induce gene expression in cultured HepaRG cells and RNA-sequencing was used to determine gene expression. Computational techniques were used to compare gene expression patterns. RESULTS: Sera from children with NAFLD induced the expression of 195 genes that were significantly differentially expressed in hepatocytes compared to controls with obesity. NAFLD was associated with increased expression of genes promoting inflammation, collagen synthesis, and extracellular matrix remodeling. Additionally, there was lower expression of genes involved in endobiotic and xenobiotic metabolism, and downregulation of peroxisome function, oxidative phosphorylation, and xenobiotic, bile acid, and fatty acid metabolism. A 13-gene signature, including upregulation of TREM1 and MMP1 and downregulation of CYP2C9, was consistently associated with all diagnostic categories of pediatric NAFLD. CONCLUSION: The extracellular milieu of sera from children with NAFLD induced specific gene profiles distinguishable by a hepatocyte reporter system. Circulating factors may contribute to inflammation and extracellular matrix remodeling and impair xenobiotic and endobiotic metabolism in pediatric NAFLD.

Meeting report of the 4th European biotransformation workshop.

1. Challenges, strategies and new technologies in the field of biotransformation were presented and discussed at the fourth European Biotransformation Workshop which was held in collaboration with the joint ISSX/DMDG meeting on June 15, 2023 at the University of Hertfordshire in Hatfield, UK.2. In this meeting report we summarise the presentations and discussions from this workshop.3. The topics covered are listed below: Unusual biotransformation reactionsBiotransformation Workflows in Discovery utilising various softwares for structure elucidationBiotransformation software for the identification of peptide metabolitesAccelerator Mass Spectrometry (AMS) for endogenous and xenobiotic metabolite profilingMetabolite profiling using quantitative Nuclear magnetic resonance (NMR) and liquid chromatography coupled to inductively coupled plasma-mass spectrometry (LC-ICP-MS).

Exploring precision-cut liver slices for comparative xenobiotic metabolism profiling in swine and cattle.

In vitro systems are useful tools for unravelling species differences in xenobiotic metabolism.The current work aimed to validate the technique of precision-cut liver slices (PCLS) for comparative studies on xenobiotic metabolism in swine and cattle.PCLS from swine (n = 3) and cattle (n = 3) were produced using a Brendel-VitronTM Tissue Slicer and cultured for 6 h. Tissue viability was preserved throughout the whole culture period.Metabolic viability was evaluated using the anthelmintics albendazole (ABZ) and fenbendazole (FBZ) as model drugs, as well as other substrates of hepatic monooxygenases: benzydamine (BZ) N-oxygenase (FMO-dependent), and the O-dealkylations of 7-ethoxyresorufin (EROD, CYP1A1-dependent) and 7-methoxyresorufin (MROD, CYP1A2-dependent).ABZ S-oxygenation resulted 6-fold (cattle) and 13.6-fold (swine) higher (p = 0.001) compared to FBZ S-oxygenation.Similar BZ N-oxygenation and EROD activities were observed in PCLS cultures from both species. MROD was 2.5-fold higher (p = 0.033) in swine than in cattle. Similarly, ABZ S-oxygenation was 1.7-fold higher (p = 0.0002) in swine than in cattle. Conversely, a 82% higher (p = 0.0003) rate of FBZ S-oxygenation was evidenced in PCLS cultures from cattle compared to those from swine.Overall, this work shows that PCLS cultures are useful to obtain relevant information on species differences in xenobiotic metabolism.

Mechanism-based inactivation of cytochromes P450: implications in drug interactions and pharmacotherapy.

Cytochrome P40 (CYP) enzymes dominate the metabolism of numerous endogenous and xenobiotic substances. While it is commonly believed that CYP-catalysed reactions result in the detoxication of foreign substances, these reactions can also yield reactive intermediates that can bind to cellular macromolecules to cause cytotoxicity or irreversibly inactivate CYPs that create them.Mechanism-based inactivation (MBI) produces either irreversible or quasi-irreversible inactivation and is commonly caused by CYP metabolic bioactivation to an electrophilic reactive intermediate. Many drugs that have been known to cause MBI in CYPs have been discovered as perpetrators in drug-drug interactions throughout the last 20-30 years.This review will highlight the key findings from the recent literature about the mechanisms of CYP enzyme inhibition, with a focus on the broad mechanistic elements of MBI for widely used drugs linked to the phenomenon. There will also be a brief discussion of the clinical or pharmacokinetic consequences of CYP inactivation with regard to drug interaction and toxicity risk.Gaining knowledge about the selective inactivation of CYPs by common therapeutic drugs helps with the assessment of factors that affect the systemic clearance of co-administered drugs and improves comprehension of anticipated interactions with other drugs or xenobiotics.

Harnessing machine learning to predict cytochrome P450 inhibition through molecular properties.

Cytochrome P450 enzymes are a superfamily of enzymes responsible for the metabolism of a variety of medicines and xenobiotics. Among the Cytochrome P450 family, five isozymes that include 1A2, 2C9, 2C19, 2D6, and 3A4 are most important for the metabolism of xenobiotics. Inhibition of any of these five CYP isozymes causes drug-drug interactions with high pharmacological and toxicological effects. So, the inhibition or non-inhibition prediction of these isozymes is of great importance. Many techniques based on machine learning and deep learning algorithms are currently being used to predict whether these isozymes will be inhibited or not. In this study, three different molecular or substructural properties that include Morgan, MACCS and Morgan (combined) and RDKit of the various molecules are used to train a distinct SVM model against each isozyme (1A2, 2C9, 2C19, 2D6, and 3A4). On the independent dataset, Morgan fingerprints provided the best results, while MACCS and Morgan (combined) achieved comparable results in terms of balanced accuracy (BA), sensitivity (Sn), and Mathews correlation coefficient (MCC). For the Morgan fingerprints, balanced accuracies (BA), Mathews correlation coefficients (MCC), and sensitivities (Sn) against each CYPs isozyme, 1A2, 2C9, 2C19, 2D6, and 3A4 on an independent dataset ranged between 0.81 and 0.85, 0.61 and 0.70, 0.72 and 0.83, respectively. Similarly, on the independent dataset, MACCS and Morgan (combined) fingerprints achieved competitive results in terms of balanced accuracies (BA), Mathews correlation coefficients (MCC), and sensitivities (Sn) against each CYPs isozyme, 1A2, 2C9, 2C19, 2D6, and 3A4, which ranged between 0.79 and 0.85, 0.59 and 0.69, 0.69 and 0.82, respectively.

Deciphering Toxic Pollutants Breakdown Potential in Microbial Community of Chumathang Hot Spring, Ladakh, India via Shotgun Metagenome Sequencing.

Persistent Organic Pollutants (POPs) have been in focus of research due to their massive contamination of environment and bio-accumulation. Bioremediation and high-throughput research have gained momentum to curb the harmful effects of POPs. The present research has explored the microbial diversity of Chumathang Hot Spring, Ladakh, India, through Illumina metagenomic HiSeq 4000 sequencing platform and their potential to degrade persistent pollutants, especially xenobiotics. Taxonomic characterization based on raw metagenomic data illuminated the abundance of members of Pseudomonadota and Actinomyceota. The re-construction of the microbial genomes from assembled contigs and scaffolds using de novo assembler metaSPAdes and their further annotation through contig alignment with available reference genomes elucidated the landscape of the hot spring's microbes. The predominantly occupied key genera reported were Pannonibacter and Novosphingobium. Comparative genomic analysis established evolutionary relationships and functional diversities among hot spring microbial communities. The function annotation through MG-RAST has revealed their metabolic versatility of degrading a wide array of xenobiotic compounds, including caprolactam, dioxin, chlorobenzene, benzoate, and. Further, the hydroxylating dioxygenase (Saro_3901) was identified as a pivotal component in the aromatic degradation pathways, showcasing extensive metabolic interconnectivity. Interestingly, protein interaction network analysis identified hub genes like Saro_1233 (protocatechuate 4,5-dioxygenase alpha subunit), while Saro_3057 (amidase) was noted for its critical role in network communication and control. The resilience of thermal ecosystems, evidenced by robust enzymatic activity and degradation capability among organisms with < 95% genetic similarity, underscores their potential for industrial and bioremediation exploration, emphasizing the importance of preserving and studying biodiverse habitats.

The response of Anisakis simplex (s. s.) to anthelmintics - Specific changes in xenobiotic metabolic processes.

Anisakiasis is a parasitic disease transmitted through the consumption of raw or undercooked fish and cephalopods that are infected with larvae of Anisakis simplex (sensu stricto) or Anisakis pegreffii. The purpose of this study was to investigate how A. simplex (s. s.) responds to the influence of anthelmintics such as ivermectin (IVM) and pyrantel (PYR). In vitro experiments were conducted using larvae at two developmental stages of A. simplex (s. s.) (L3 and L4) obtained from Baltic herring (Clupea harengus membras). Larvae were cultured with different concentrations of IVM or PYR (1.56, 3.125, and 6.25 µg/mL) for various durations (3, 6, 9, and 12 h) under anaerobic conditions (37 °C, 5% CO2). The gene expression of actin, ABC transporter, antioxidant enzymes, γ-aminobutyric acid receptors, and nicotinic acetylcholine receptors, as well as the oxidative status were analyzed. The results showed that A. simplex (s. s.) L3 stage had lower mobility when cultured with PYR compared to IVM. The analysis of relative gene expression revealed significant differences in the mRNA level of ABC transporters after treatment with IVM and PYR, compared to the control group. Similar patterns were observed in the gene expression of antioxidant enzymes in response to both drugs. Furthermore, the total antioxidant capacity (TAC) and glutathione S-transferase (GST) activity were higher in the treatment groups than in the control group. These findings suggest a relationship between the expression of the studied genes, including those related to oxidative metabolism, and the effectiveness of the tested drugs.

BioTransformer 3.0-a web server for accurately predicting metabolic transformation products.

BioTransformer 3.0 (https://biotransformer.ca) is a freely available web server that supports accurate, rapid and comprehensive in silico metabolism prediction. It combines machine learning approaches with a rule-based system to predict small-molecule metabolism in human tissues, the human gut as well as the external environment (soil and water microbiota). Simply stated, BioTransformer takes a molecular structure as input (SMILES or SDF) and outputs an interactively sortable table of the predicted metabolites or transformation products (SMILES, PNG images) along with the enzymes that are predicted to be responsible for those reactions and richly annotated downloadable files (CSV and JSON). The entire process typically takes less than a minute. Previous versions of BioTransformer focused exclusively on predicting the metabolism of xenobiotics (such as plant natural products, drugs, cosmetics and other synthetic compounds) using a limited number of pre-defined steps and somewhat limited rule-based methods. BioTransformer 3.0 uses much more sophisticated methods and incorporates new databases, new constraints and new prediction modules to not only more accurately predict the metabolic transformation products of exogenous xenobiotics but also the transformation products of endogenous metabolites, such as amino acids, peptides, carbohydrates, organic acids, and lipids. BioTransformer 3.0 can also support customized sequential combinations of these transformations along with multiple iterations to simulate multi-step human biotransformation events. Performance tests indicate that BioTransformer 3.0 is 40-50% more accurate, far less prone to combinatorial 'explosions' and much more comprehensive in terms of metabolite coverage/capabilities than previous versions of BioTransformer.