Absolute quantification of intracellular metabolite pools is a prerequisite for modeling and in-depth biological interpretation of metabolomics data. It is the final step of an elaborate metabolomics workflow, with challenges associated with all steps—from sampling to quantifying the physicochemically diverse metabolite pool. Chromatographic separation combined with mass spectrometric (MS) detection is the superior platform for high coverage, selective, and sensitive detection of metabolites. Herein, we apply our quantitative MS-metabolomics workflow to measure and present the central carbon metabolome of a panel of commonly applied biological model systems. The workflow includes three chromatographic methods combined with isotope dilution tandem mass spectrometry to allow for absolute quantification of 68 metabolites of glycolysis, the pentose phosphate pathway, the tricarboxylic acid cycle, and the amino acid and (deoxy) nucleoside pools. The biological model systems; Bacillus subtilis, Saccharomyces cerevisiae, two microalgal species, and four human cell lines were all cultured in commonly applied culture media and sampled in exponential growth phase. Both literature and databases are scarce with comprehensive metabolite datasets, and existing entries range over several orders of magnitude. The workflow and metabolite panel presented herein can be employed to expand the list of reference metabolomes, as encouraged by the metabolomics community, in a continued effort to develop and refine high-quality quantitative metabolomics workflows.

Precise and accurate quantification is a prerequisite for interpretation of targeted metabolomics data, but this task is challenged by the inherent instability of the analytes. The sampling, quenching, extraction, and sample purification conditions required to recover and stabilize metabolites in representative extracts have also been proven highly dependent on species-specific properties. For Escherichia coli, unspecific leakage has been demonstrated for conventional microbial metabolomics sampling protocols. We herein present a fast filtration-based sampling protocol for this widely applied model organism, focusing on pitfalls such as inefficient filtration, selective loss of biomass, matrix contamination, and membrane permeabilization and leakage. We evaluate the effect of and need for removal of extracellular components and demonstrate how residual salts can challenge analytical accuracy of hyphenated mass spectrometric analyses, even when sophisticated correction strategies are applied. Laborious extraction procedures are bypassed by direct extraction in cold acetonitrile:water:methanol (3:5:2, v/v%), ensuring compatibility with sample concentration and thus, any downstream analysis. By applying this protocol, we achieve and demonstrate high precision and low metabolite turnover, and, followingly, minimal perturbation of the inherent metabolic state. This allows us to herein report absolute intracellular concentrations in E. coli and explore its central carbon metabolome at several commonly applied cultivation conditions.

Genome-scale metabolic models (GEMs) are mathematical representations of metabolism that allow for in silico simulation of metabolic phenotypes and capabilities. A prerequisite for these predictions is an accurate representation of the biomolecular composition of the cell necessary for replication and growth, implemented in GEMs as the so-called biomass objective function (BOF). The BOF contains the metabolic precursors required for synthesis of the cellular macro- and micromolecular constituents (e.g. protein, RNA, DNA), and its composition is highly dependent on the particular organism, strain, and growth condition. Despite its critical role, the BOF is rarely constructed using specific measurements of the modeled organism, drawing the validity of this approach into question. Thus, there is a need to establish robust and reliable protocols for experimental condition-specific biomass determination. Here, we address this challenge by presenting a general pipeline for biomass quantification, evaluating its performance on Escherichia coli K-12 MG1655 sampled during balanced exponential growth under controlled conditions in a batch-fermentor set-up. We significantly improve both the coverage and molecular resolution compared to previously published workflows, quantifying 91.6% of the biomass. Our measurements display great correspondence with previously reported measurements, and we were also able to detect subtle characteristics specific to the particular E. coli strain. Using the modified E. coli GEM i ML1515a, we compare the feasible flux ranges of our experimentally determined BOF with the original BOF, finding that the changes in BOF coefficients considerably affect the attainable fluxes at the genome-scale.