Catherine L. Winder

Global metabolic profiling of Escherichia coli cultures: An evaluation of methods for quenching and extraction of intracellular metabolites

Metabolomics and systems biology require the acquisition of reproducible, robust, reliable, and homogeneous biological data sets. Therefore, we developed and validated standard operating procedures (SOPs) for quenching and efficient extraction of metabolites from Escherichia coli to determine the best methods to approach global analysis of the metabolome. E. coli was grown in chemostat culture so that cellular metabolism could be held in reproducible, steady-state conditions under a range of precisely defined growth conditions, thus enabling sufficient replication of samples. The metabolome profiles were generated using gas chromatography/time-of-flight mass spectrometry (GC/TOF-MS). We employed univariate and multivariate statistical analyses to determine the most suitable method. This investigation indicates that 60% cold (-48 degrees C) methanol solution is the most appropriate method to quench metabolism, and we recommend 100% methanol, also at -48 degrees C, with multiple freeze-thaw cycles for the extraction of metabolites. However, complementary extractions would be necessary for coverage of the entire complement of metabolites as detected by GC/TOF-MS. Finally, the observation that metabolite leakage was significant and measurable whichever quenching method is used indicates that methods should be incorporated into the experiment to facilitate the accurate quantification of intracellular metabolites.

Mass spectrometry tools and metabolite-specific databases for molecular identification in metabolomics

The chemical identification of mass spectrometric signals in metabolomic applications is important to provide conversion of analytical data to biological knowledge about metabolic pathways. The complexity of electrospray mass spectrometric data acquired from a range of samples (serum, urine, yeast intracellular extracts, yeast metabolic footprints, placental tissue metabolic footprints) has been investigated and has defined the frequency of different ion types routinely detected. Although some ion types were expected (protonated and deprotonated peaks, isotope peaks, multiply charged peaks) others were not expected (sodium formate adduct ions). In parallel, the Manchester Metabolomics Database (MMD) has been constructed with data from genome scale metabolic reconstructions, HMDB, KEGG, Lipid Maps, BioCyc and DrugBank to provide knowledge on 42,687 endogenous and exogenous metabolite species. The combination of accurate mass data for a large collection of metabolites, theoretical isotope abundance data and knowledge of the different ion types detected provided a greater number of electrospray mass spectrometric signals which were putatively identified and with greater confidence in the samples studied. To provide definitive identification metabolite-specific mass spectral libraries for UPLC-MS and GC-MS have been constructed for 1,065 commercially available authentic standards. The MMD data are available at http://dbkgroup.org/MMD/.

Improving metabolic flux predictions using absolute gene expression data.

BackgroundConstraint-based analysis of genome-scale metabolic models typically relies upon maximisation of a cellular objective function such as the rate or efficiency of biomass production. Whilst this assumption may be valid in the case of microorganisms growing under certain conditions, it is likely invalid in general, and especially for multicellular organisms, where cellular objectives differ greatly both between and within cell types. Moreover, for the purposes of biotechnological applications, it is normally the flux to a specific metabolite or product that is of interest rather than the rate of production of biomass per se.ResultsAn alternative objective function is presented, that is based upon maximising the correlation between experimentally measured absolute gene expression data and predicted internal reaction fluxes. Using quantitative transcriptomics data acquired from Saccharomyces cerevisiae cultures under two growth conditions, the method outperforms traditional approaches for predicting experimentally measured exometabolic flux that are reliant upon maximisation of the rate of biomass production.ConclusionDue to its improved prediction of experimentally measured metabolic fluxes, and of its lack of a requirement for knowledge of the biomass composition of the organism under the conditions of interest, the approach is likely to be of rather general utility. The method has been shown to predict fluxes reliably in single cellular systems. Subsequent work will investigate the method’s ability to generate condition- and tissue-specific flux predictions in multicellular organisms.

Systems Biology: The elements and principles of Life

Systems Biology has a mission that puts it at odds with traditional paradigms of physics and molecular biology, such as the simplicity requested by Occams razor and minimum energy/maximal efficiency. By referring to biochemical experiments on control and regulation, and on flux balancing in yeast, we show that these paradigms are inapt. Systems Biology does not quite converge with biology either: Although it certainly requires accurate ‘stamp collecting’, it discovers quantitative laws. Systems Biology is a science of its own, discovering own fundamental principles, some of which we identify here.

A model of yeast glycolysis based on a consistent kinetic characterisation of all its enzymes

We present an experimental and computational pipeline for the generation of kinetic models of metabolism, and demonstrate its application to glycolysis in Saccharomyces cerevisiae. Starting from an approximate mathematical model, we employ a “cycle of knowledge” strategy, identifying the steps with most control over flux. Kinetic parameters of the individual isoenzymes within these steps are measured experimentally under a standardised set of conditions. Experimental strategies are applied to establish a set of in vivo concentrations for isoenzymes and metabolites. The data are integrated into a mathematical model that is used to predict a new set of metabolite concentrations and reevaluate the control properties of the system. This bottom‐up modelling study reveals that control over the metabolic network most directly involved in yeast glycolysis is more widely distributed than previously thought.

Absolute Quantification of the Glycolytic Pathway in Yeast: DEPLOYMENT OF A COMPLETE QconCAT APPROACH

The availability of label-free data derived from yeast cells (based on the summed intensity of the three strongest, isoform-specific peptides) permitted a preliminary assessment of protein abundances for glycolytic proteins. Following this analysis, we demonstrate successful application of the QconCAT technology, which uses recombinant DNA techniques to generate artificial concatamers of large numbers of internal standard peptides, to the quantification of enzymes of the glycolysis pathway in the yeast Saccharomyces cerevisiae. A QconCAT of 88 kDa (59 tryptic peptides) corresponding to 27 isoenzymes was designed and built to encode two or three analyte peptides per protein, and after stable isotope labeling of the standard in vivo, protein levels were determined by LC-MS, using ultra high performance liquid chromatography-coupled mass spectrometry. We were able to determine absolute protein concentrations between 14,000 and 10 million molecules/cell. Issues such as efficiency of extraction and completeness of proteolysis are addressed, as well as generic factors such as optimal quantotypic peptide selection and expression. In addition, the same proteins were quantified by intensity-based label-free analysis, and both sets of data were compared with other quantification methods.

TARDIS-based microbial metabolomics: time and relative differences in systems.

Metabolomics can play a particularly important role in elucidating novel anabolic and catabolic pathways in bacteria and fungi, and in understanding the dynamics of metabolism. In these approaches, an isotopically labelled substrate, with an artificially high abundance of isotopic label, is fed to the microorganism under study. The products become isotopically labelled, and can be measured using a combination of mass spectrometry and nuclear magnetic resonance spectroscopy. This mass isotopomer analysis is referred to as time and relative differences in systems (TARDIS)-based analysis, as it measures and quantifies the temporal sequential emergence of these labelled products. In this review, we cover this topic from an experimental point of view in relation to the study of metabolism, and summarise how the application of radioactive and stable isotopes is being used in pathway elucidation and metabolic flux determination (fluxomics).

Raman spectroscopy: lighting up the future of microbial identification.

Over the last decade Raman spectroscopy has become established as a physicochemical technique for the rapid identification of microbes. This powerful analytical method generates a spectroscopic fingerprint from the microbial sample, which provides quantitative and qualitative information that can be used to characterize, discriminate and identify microorganisms, in both bacteria slurry and at the single-cell level. Recent developments in Raman spectroscopy have dramatically increased in recent years due to the enhancement of the signal by techniques including tip-enhanced Raman spectroscopy and coherent anti-Stokes Raman spectroscopy and due to the availability of user-friendly instrumentation and software. The result of this has been reduced cost and rapid collection time, and it has allowed the nonspecialist access to this physical sciences approach for biological applications. In this article, we will briefly explain the technique of Raman spectroscopy and discuss enhancement techniques, including the recent application of tip-enhanced Raman spectroscopy to microbiology, as well as the move towards rapid microbial identification with Raman spectroscopy. Furthermore, recent studies have combined Raman spectroscopy with microfluidic devices, giving greater control of sample conditions, which will no doubt have an important impact in the future development of Raman spectroscopy for microbial identification.

VOC-based metabolic profiling for food spoilage detection with the application to detecting Salmonella typhimurium-contaminated pork

In this study, we investigated the feasibility of using a novel volatile organic compound (VOC)-based metabolic profiling approach with a newly devised chemometrics methodology which combined rapid multivariate analysis on total ion currents with in-depth peak deconvolution on selected regions to characterise the spoilage progress of pork. We also tested if such approach possessed enough discriminatory information to differentiate natural spoiled pork from pork contaminated with Salmonella typhimurium, a food poisoning pathogen commonly recovered from pork products. Spoilage was monitored in this study over a 72-h period at 0-, 24-, 48- and 72-h time points after the artificial contamination with the salmonellae. At each time point, the VOCs from six individual pork chops were collected for spoiled vs. contaminated meat. Analysis of the VOCs was performed by gas chromatography/mass spectrometry (GC/MS). The data generated by GC/MS analysis were initially subjected to multivariate analysis using principal component analysis (PCA) and multi-block PCA. The loading plots were then used to identify regions in the chromatograms which appeared important to the separation shown in the PCA/multi-block PCA scores plot. Peak deconvolution was then performed only on those regions using a modified hierarchical multivariate curve resolution procedure for curve resolution to generate a concentration profiles matrix C and the corresponding pure spectra matrix S. Following this, the pure mass spectra (S) of the peaks in those region were exported to NIST 02 mass library for chemical identification. A clear separation between the two types of samples was observed from the PCA models, and after deconvolution and univariate analysis using N-way ANOVA, a total of 16 significant metabolites were identified which showed difference between natural spoiled pork and those contaminated with S. typhimurium.

Systematic integration of experimental data and models in systems biology

BackgroundThe behaviour of biological systems can be deduced from their mathematical models. However, multiple sources of data in diverse forms are required in the construction of a model in order to define its components and their biochemical reactions, and corresponding parameters. Automating the assembly and use of systems biology models is dependent upon data integration processes involving the interoperation of data and analytical resources.ResultsTaverna workflows have been developed for the automated assembly of quantitative parameterised metabolic networks in the Systems Biology Markup Language (SBML). A SBML model is built in a systematic fashion by the workflows which starts with the construction of a qualitative network using data from a MIRIAM-compliant genome-scale model of yeast metabolism. This is followed by parameterisation of the SBML model with experimental data from two repositories, the SABIO-RK enzyme kinetics database and a database of quantitative experimental results. The models are then calibrated and simulated in workflows that call out to COPASIWS, the web service interface to the COPASI software application for analysing biochemical networks. These systems biology workflows were evaluated for their ability to construct a parameterised model of yeast glycolysis.ConclusionsDistributed information about metabolic reactions that have been described to MIRIAM standards enables the automated assembly of quantitative systems biology models of metabolic networks based on user-defined criteria. Such data integration processes can be implemented as Taverna workflows to provide a rapid overview of the components and their relationships within a biochemical system.