Open-Source Chromatographic Data Analysis for Reaction Optimization and Screening.

Automation and digitalization solutions in the field of small molecule synthesis face new challenges for chemical reaction analysis, especially in the field of high-performance liquid chromatography (HPLC). Chromatographic data remains locked in vendors' hardware and software components, limiting their potential in automated workflows and data science applications. In this work, we present an open-source Python project called MOCCA for the analysis of HPLC-DAD (photodiode array detector) raw data. MOCCA provides a comprehensive set of data analysis features, including an automated peak deconvolution routine of known signals, even if overlapped with signals of unexpected impurities or side products. We highlight the broad applicability of MOCCA in four studies: (i) a simulation study to validate MOCCA's data analysis features; (ii) a reaction kinetics study on a Knoevenagel condensation reaction demonstrating MOCCA's peak deconvolution feature; (iii) a closed-loop optimization study for the alkylation of 2-pyridone without human control during data analysis; (iv) a well plate screening of categorical reaction parameters for a novel palladium-catalyzed cyanation of aryl halides employing O-protected cyanohydrins. By publishing MOCCA as a Python package with this work, we envision an open-source community project for chromatographic data analysis with the potential of further advancing its scope and capabilities.

A Pareto approach to resolve the conflict between information gain and experimental costs: Multiple-criteria design of carbon labeling experiments.

Science revolves around the best way of conducting an experiment to obtain insightful results. Experiments with maximal information content can be found by computational experimental design (ED) strategies that identify optimal conditions under which to perform the experiment. Several criteria have been proposed to measure the information content, each emphasizing different aspects of the design goal, i.e., reduction of uncertainty. Where experiments are complex or expensive, second sight is at the budget governing the achievable amount of information. In this context, the design objectives cost and information gain are often incommensurable, though dependent. By casting the ED task into a multiple-criteria optimization problem, a set of trade-off designs is derived that approximates the Pareto-frontier which is instrumental for exploring preferable designs. In this work, we present a computational methodology for multiple-criteria ED of information-rich experiments that accounts for virtually any set of design criteria. The methodology is implemented for the case of 13C metabolic flux analysis (MFA), which is arguably the most expensive type among the 'omics' technologies, featuring dozens of design parameters (tracer composition, analytical platform, measurement selection etc.). Supported by an innovative visualization scheme, we demonstrate with two realistic showcases that the use of multiple criteria reveals deep insights into the conflicting interplay between information carriers and cost factors that are not amendable to single-objective ED. For instance, tandem mass spectrometry turns out as best-in-class with respect to information gain, while it delivers this information quality cheaper than the other, routinely applied analytical technologies. Therewith, our Pareto approach to ED offers the investigator great flexibilities in the conception phase of a study to balance costs and benefits.

Stage-specific metabolic features of differentiating neurons: Implications for toxicant sensitivity.

Developmental neurotoxicity (DNT) may be induced when chemicals disturb a key neurodevelopmental process, and many tests focus on this type of toxicity. Alternatively, DNT may occur when chemicals are cytotoxic only during a specific neurodevelopmental stage. The toxicant sensitivity is affected by the expression of toxicant targets and by resilience factors. Although cellular metabolism plays an important role, little is known how it changes during human neurogenesis, and how potential alterations affect toxicant sensitivity of mature vs. immature neurons. We used immature (d0) and mature (d6) LUHMES cells (dopaminergic human neurons) to provide initial answers to these questions. Transcriptome profiling and characterization of energy metabolism suggested a switch from predominantly glycolytic energy generation to a more pronounced contribution of the tricarboxylic acid cycle (TCA) during neuronal maturation. Therefore, we used pulsed stable isotope-resolved metabolomics (pSIRM) to determine intracellular metabolite pool sizes (concentrations), and isotopically non-stationary 13C-metabolic flux analysis (INST 13C-MFA) to calculate metabolic fluxes. We found that d0 cells mainly use glutamine to fuel the TCA. Furthermore, they rely on extracellular pyruvate to allow continuous growth. This metabolic situation does not allow for mitochondrial or glycolytic spare capacity, i.e. the ability to adapt energy generation to altered needs. Accordingly, neuronal precursor cells displayed a higher sensitivity to several mitochondrial toxicants than mature neurons differentiated from them. In summary, this study shows that precursor cells lose their glutamine dependency during differentiation while they gain flexibility of energy generation and thereby increase their resistance to low concentrations of mitochondrial toxicants.

Natural isotope correction of MS/MS measurements for metabolomics and (13)C fluxomics.

Fluxomics and metabolomics are crucial tools for metabolic engineering and biomedical analysis to determine the in vivo cellular state. Especially, the application of (13)C isotopes allows comprehensive insights into the functional operation of cellular metabolism. Compared to single MS, tandem mass spectrometry (MS/MS) provides more detailed and accurate measurements of the metabolite enrichment patterns (tandem mass isotopomers), increasing the accuracy of metabolite concentration measurements and metabolic flux estimation. MS-type data from isotope labeling experiments is biased by naturally occurring stable isotopes (C, H, N, O, etc.). In particular, GC-MS(/MS) requires derivatization for the usually non-volatile intracellular metabolites introducing additional natural isotopes leading to measurements that do not directly represent the carbon labeling distribution. To make full use of LC- and GC-MS/MS mass isotopomer measurements, the influence of natural isotopes has to be eliminated (corrected). Our correction approach is analyzed for the two most common applications; (13)C fluxomics and isotope dilution mass spectrometry (IDMS) based metabolomics. Natural isotopes can have an impact on the calculated flux distribution which strongly depends on the substrate labeling and the actual flux distribution. Second, we show that in IDMS based metabolomics natural isotopes lead to underestimated concentrations that can and should be corrected with a nonlinear calibration. Our simulations indicate that the correction for natural abundance in isotope based fluxomics and quantitative metabolomics is essential for correct data interpretation.

Fluxome study of Pseudomonas fluorescens reveals major reorganisation of carbon flux through central metabolic pathways in response to inactivation of the anti-sigma factor MucA.

BACKGROUND: The bacterium Pseudomonas fluorescens switches to an alginate-producing phenotype when the pleiotropic anti-sigma factor MucA is inactivated. The inactivation is accompanied by an increased biomass yield on carbon sources when grown under nitrogen-limited chemostat conditions. A previous metabolome study showed significant changes in the intracellular metabolite concentrations, especially of the nucleotides, in mucA deletion mutants compared to the wild-type. In this study, the P. fluorescens SBW25 wild-type and an alginate non-producing mucA- ΔalgC double-knockout mutant are investigated through model-based (13)C-metabolic flux analysis ((13)C-MFA) to explore the physiological consequences of MucA inactivation at the metabolic flux level. Intracellular metabolite extracts from three carbon labelling experiments using fructose as the sole carbon source are analysed for (13)C-label incorporation in primary metabolites by gas and liquid chromatography tandem mass spectrometry. RESULTS: From mass isotopomer distribution datasets, absolute intracellular metabolic reaction rates for the wild type and the mutant are determined, revealing extensive reorganisation of carbon flux through central metabolic pathways in response to MucA inactivation. The carbon flux through the Entner-Doudoroff pathway was reduced in the mucA- ΔalgC mutant, while flux through the pentose phosphate pathway was increased. Our findings also indicated flexibility of the anaplerotic reactions through down-regulation of the pyruvate shunt in the mucA- ΔalgC mutant and up-regulation of the glyoxylate shunt. CONCLUSIONS: Absolute metabolic fluxes and metabolite levels give detailed, integrated insight into the physiology of this industrially, medically and agriculturally important bacterial species and suggest that the most efficient way of using a mucA- mutant as a cell factory for alginate production would be to use non-growing conditions and nitrogen deprivation.

Metabolic engineering of Corynebacterium glutamicum for methanol metabolism.

Methanol is already an important carbon feedstock in the chemical industry, but it has found only limited application in biotechnological production processes. This can be mostly attributed to the inability of most microbial platform organisms to utilize methanol as a carbon and energy source. With the aim to turn methanol into a suitable feedstock for microbial production processes, we engineered the industrially important but nonmethylotrophic bacterium Corynebacterium glutamicum toward the utilization of methanol as an auxiliary carbon source in a sugar-based medium. Initial oxidation of methanol to formaldehyde was achieved by heterologous expression of a methanol dehydrogenase from Bacillus methanolicus, whereas assimilation of formaldehyde was realized by implementing the two key enzymes of the ribulose monophosphate pathway of Bacillus subtilis: 3-hexulose-6-phosphate synthase and 6-phospho-3-hexuloisomerase. The recombinant C. glutamicum strain showed an average methanol consumption rate of 1.7 ± 0.3 mM/h (mean ± standard deviation) in a glucose-methanol medium, and the culture grew to a higher cell density than in medium without methanol. In addition, [(13)C]methanol-labeling experiments revealed labeling fractions of 3 to 10% in the m + 1 mass isotopomers of various intracellular metabolites. In the background of a C. glutamicum Δald ΔadhE mutant being strongly impaired in its ability to oxidize formaldehyde to CO2, the m + 1 labeling of these intermediates was increased (8 to 25%), pointing toward higher formaldehyde assimilation capabilities of this strain. The engineered C. glutamicum strains represent a promising starting point for the development of sugar-based biotechnological production processes using methanol as an auxiliary substrate.

How to measure metabolic fluxes: a taxonomic guide for (13)C fluxomics.

Metabolic reaction rates (fluxes) contribute fundamentally to our understanding of metabolic phenotypes and mechanisms of cellular regulation. Stable isotope-based fluxomics integrates experimental data with biochemical networks and mathematical modeling to 'measure' the in vivo fluxes within an organism that are not directly observable. In recent years, (13)C fluxomics has evolved into a technology with great experimental, analytical, and mathematical diversity. This review aims at establishing a unified taxonomy by means of which the various fluxomics methods can be compared to each other. By linking the developed modeling approaches to recent studies, their challenges and opportunities are put into perspective. The proposed classification serves as a guide for scientific 'travelers' who are striving to resolve research questions with the currently available (13)C fluxomics toolset.

13C-flux spectral analysis of host-pathogen metabolism reveals a mixed diet for intracellular Mycobacterium tuberculosis.

Whereas intracellular carbon metabolism has emerged as an attractive drug target, the carbon sources of intracellularly replicating pathogens, such as the tuberculosis bacillus Mycobacterium tuberculosis, which causes long-term infections in one-third of the world's population, remain mostly unknown. We used a systems-based approach--(13)C-flux spectral analysis (FSA) complemented with manual analysis-to measure the metabolic interaction between M. tuberculosis and its macrophage host cell. (13)C-FSA analysis of experimental data showed that M. tuberculosis obtains a mixture of amino acids, C1 and C2 substrates from its host cell. We experimentally confirmed that the C1 substrate was derived from CO2. (13)C labeling experiments performed on a phosphoenolpyruvate carboxykinase mutant revealed that intracellular M. tuberculosis has access to glycolytic C3 substrates. These findings provide constraints for developing novel chemotherapeutics.

13CFLUX2--high-performance software suite for (13)C-metabolic flux analysis.

SUMMARY: (13)C-based metabolic flux analysis ((13)C-MFA) is the state-of-the-art method to quantitatively determine in vivo metabolic reaction rates in microorganisms. 13CFLUX2 contains all tools for composing flexible computational (13)C-MFA workflows to design and evaluate carbon labeling experiments. A specially developed XML language, FluxML, highly efficient data structures and simulation algorithms achieve a maximum of performance and effectiveness. Support of multicore CPUs, as well as compute clusters, enables scalable investigations. 13CFLUX2 outperforms existing tools in terms of universality, flexibility and built-in features. Therewith, 13CFLUX2 paves the way for next-generation high-resolution (13)C-MFA applications on the large scale. AVAILABILITY AND IMPLEMENTATION: 13CFLUX2 is implemented in C++ (ISO/IEC 14882 standard) with Java and Python add-ons to run under Linux/Unix. A demo version and binaries are available at www.13cflux.net.

Visualizing multi-omics data in metabolic networks with the software Omix: a case study.

Systems Biology is a multi-disciplinary research field with the aim of understanding the function of complex processes in living organisms. These intracellular processes are described by biochemical networks. Experimental studies in alliance with computer simulation lead to a continually increasing amount of data in liaison with different layers of biochemical networks. Thus, visualization is very important for getting an overview of data in association with the network components. Omix is a software for the visualization of any data in biochemical networks. The unique feature of Omix is: the software is programmable by a scripting language called Omix Visualization Language (OVL). In Omix, the visualization of data coming from experiment or simulation is completely performed by the software user realized in concise OVL scripts. By this, visualization becomes most flexible and adaptable to the requirements of the user and can be adapted to new application fields. We present four case studies of visualizing data of diverse kind in biochemical networks on metabolic level by using Omix and the OVL scripting language. These worked examples demonstrate the power of OVL in conjunction with pleasing visualization, an important requirement for successful interdisciplinary communication in the interface between more experimental and more theoretical researchers.