Current machine learning techniques enable robust association of biological signals with measured phenotypes, but these approaches are incapable of identifying causal relationships. Here, we develop an integrated ''white-box'' biochemical screening, network modeling, and machine learning approach for revealing causal mechanisms and apply this approach to understanding antibiotic efficacy. We counter-screen diverse metabolites against bactericidal antibiotics in Escherichia coli and simulate their corresponding metabolic states using a genomescale metabolic network model. Regression of the measured screening data on model simulations reveals that purine biosynthesis participates in antibiotic lethality, which we validate experimentally. We show that antibiotic-induced adenine limitation increases ATP demand, which elevates central carbon metabolism activity and oxygen consumption, enhancing the killing effects of antibiotics. This work demonstrates how prospective network modeling can couple with machine learning to identify complex causal mechanisms underlying drug efficacy.

Monoterpene indole alkaloids (MIAs) are a diverse family of complex plant secondary metabolites with many medicinal properties, including the essential anti-cancer therapeutics vinblastine and vincristine1. As MIAs are difficult to chemically synthesize, the world’s supply chain for vinblastine relies on low-yielding extraction and purification of the precursors vindoline and catharanthine from the plant Catharanthus roseus, which is then followed by simple in vitro chemical coupling and reduction to form vinblastine at an industrial scale2,3. Here, we demonstrate the de novo microbial biosynthesis of vindoline and catharanthine using a highly engineered yeast, and in vitro chemical coupling to vinblastine. The study showcases a very long biosynthetic pathway refactored into a microbial cell factory, including 30 enzymatic steps beyond the yeast native metabolites geranyl pyrophosphate and tryptophan to catharanthine and vindoline. In total, 56 genetic edits were performed, including expression of 34 heterologous genes from plants, as well as deletions, knock-downs and overexpression of ten yeast genes to improve precursor supplies towards de novo production of catharanthine and vindoline, from which semisynthesis to vinblastine occurs. As the vinblastine pathway is one of the longest MIA biosynthetic pathways, this study positions yeast as a scalable platform to produce more than 3,000 natural MIAs and a virtually infinite number of new-to-nature analogues.

A B S T R A C TFast metabolite quantification methods are required for high throughput screening of microbial strains obtained by combinatorial or evolutionary engineering approaches. In this study, a rapid RIP-LC-MS/MS (RapidRIP) method for high-throughput quantitative metabolomics was developed and validated that was capable of quantifying 102 metabolites from central, amino acid, energy, nucleotide, and cofactor metabolism in less than 5 minutes. The method was shown to have comparable sensitivity and resolving capability as compared to a full length RIP-LC-MS/MS method (FullRIP). The RapidRIP method was used to quantify the metabolome of seven industrial strains of E. coli revealing significant differences in glycolytic, pentose phosphate, TCA cycle, amino acid, and energy and cofactor metabolites were found. These differences translated to statistically and biologically significant differences in thermodynamics of biochemical reactions between strains that could have implications when choosing a host for bioprocessing.