A universal framework for 13C metabolic flux analysis.

A general methodology is presented for the modeling, simulation, design, evaluation, and statistical analysis of (13)C-labeling experiments for metabolic flux analysis. The universal software framework 13C-FLUX was implemented to support all steps of this process. Guided by the example of anaplerotic flux determination in Corynebacterium glutamicum, the technical details of the model setup, experimental design, and data evaluation are discussed. It is shown how the network structure, the input substrate composition, the assumptions about fluxes, and the measurement configuration are specified within 13C-FLUX. Based on the network model, different experimental designs are computed depending on the goal of the investigations. Finally, a specific experiment is evaluated and the various statistical methods used to analyze the results are briefly explained. The appendix gives some details about the software implementation and availability.

In vivo quantification of parallel and bidirectional fluxes in the anaplerosis of Corynebacterium glutamicum.

The C(3)-C(4) metabolite interconversion at the anaplerotic node in many microorganisms involves a complex set of reactions. C(3) carboxylation to oxaloacetate can originate from phosphoenolpyruvate and pyruvate, and at the same time multiple C(4)-decarboxylating enzymes may be present. The functions of such parallel reactions are not yet fully understood. Using a (13)C NMR-based strategy, we here quantify the individual fluxes at the anaplerotic node of Corynebacterium glutamicum, which is an example of a bacterium possessing multiple carboxylation and decarboxylation reactions. C. glutamicum was grown with a (13)C-labeled glucose isotopomer mixture as the main carbon source and (13)C-labeled lactate as a cosubstrate. 58 isotopomers as well as 15 positional labels of biomass compounds were quantified. Applying a generally applicable mathematical model to include metabolite mass and carbon labeling balances, it is shown that pyruvate carboxylase contributed 91 +/- 7% to C(3) carboxylation. The total in vivo carboxylation rate of 1.28 +/- 0.14 mmol/g dry weight/h exceeds the demand of carboxylated metabolites for biosyntheses 3-fold. Excess oxaloacetate was recycled to phosphoenolpyruvate by phosphoenolpyruvate carboxykinase. This shows that the reactions at the anaplerotic node might serve additional purposes other than only providing C(4) metabolites for biosynthesis.

Determination of full 13C isotopomer distributions for metabolic flux analysis using heteronuclear spin echo difference NMR spectroscopy.

13C-isotopomer labeling experiments play an increasingly important role in the analysis of intracellular metabolic fluxes for genetic engineering purposes. 13C NMR spectroscopy is a key technique in the experimental determination of isotopomer distributions. However, only subsets of isotopomers can be quantitated using this technique due to redundancies in the scalar coupling patterns and due to invisibility of the 12C isotope in NMR. Therefore, we developed and describe in this paper a 1H NMR spectroscopy method that allows to determine the complete isotopomer distribution in metabolites having a backbone consisting of up to at least four carbons. The proposed pulse sequences employ up to three alternately applied frequency-selective inversion pulses in the 13C channel. In a first application study, the complete isotopomer distribution of aspartate isolated from [1-13C]ethanol-grown Ashbya gossypii was determined. A tentative model of the central metabolism of this organism was constructed and used for metabolic flux analysis. The aspartate isotopomer NMR data played a key role in the successful determination of the flux through the peroxisomal glyoxylate pathway. The new NMR method can be highly instrumental in generating the data upon which isotopomer labeling experiments for flux analysis, that are becoming increasingly important, are based.

Bidirectional reaction steps in metabolic networks: III. Explicit solution and analysis of isotopomer labeling systems.

The last few years have brought tremendous progress in experimental methods for metabolic flux determination by carbon-labeling experiments. A significant enlargement of the available measurement data set has been achieved, especially when isotopomer fractions within intracellular metabolite pools are quantitated. This information can be used to improve the statistical quality of flux estimates. Furthermore, several assumptions on bidirectional intracellular reaction steps that were hitherto indispensable may now become obsolete. To make full use of the complete measurement information a general mathematical model for isotopomer systems is established in this contribution. Then, by introducing the important new concept of cumomers and cumomer fractions, it is shown that the arising nonlinear isotopomer balance equations can be solved analytically in all cases. In particular, the solution of the metabolite flux balances and the positional carbon-labeling balances presented in part I of this series turn out to be just the first two steps of the general solution procedure for isotopomer balances. A detailed analysis of the isotopomer network structure then opens up new insights into the intrinsic structure of isotopomer systems. In particular, it turns out that isotopomer systems are not as complex as they appear at first glance. This enables some far-reaching conclusions to be drawn on the information potential of isotopomer experiments with respect to flux identification. Finally, some illustrative examples are examined to show that an information increase is not guaranteed when isotopomer measurements are used in addition to positional enrichment data.

Bidirectional reaction steps in metabolic networks: IV. Optimal design of isotopomer labeling experiments.

This article generalizes the statistical tools for the evaluation of carbon-labeling experiments that have been developed for the case of positional enrichment systems in part II of this series to the general case of isotopomer systems. For this purpose, a new generalized measurement equation is introduced that can describe all kinds of measured data, like positional enrichments, relative (13)C nuclear magnetic resonance ((13)C NMR) multiplet intensities, or mass isotopomer fractions produced with mass spectroscopy (MS) instruments. Then, to facilitate the specification of the various measurement procedures available, a new flexible textual notation is introduced from which the complicated generalized measurement equations are generated automatically. Based on these measurement equations, a statistically optimal flux estimator is established and parameter covariance matrices for the flux estimation are computed. Having implemented these tools, different kinds of labeling experiments can be compared by using statistical quality measures. A general framework for the optimal design of carbon-labeling experiments is established on the basis of this method. As an example it is applied to the Corynebacterium network from part II extended by various NMR and MS measurements. In particular, the positional enrichment, multiplet, or mass isotopomer measurements with the greatest information content for flux estimation are computed (measurement design) and various differently labeled input substrates are compared with respect to flux estimation (input design). It is examined in detail how the measurement procedure influences the estimation quality of specific fluxes like the pentose phosphate pathway influx.

Object-oriented programming for the biosciences.

The development of software systems for the biosciences is always closely connected to experimental practice. Programs must be able to handle the inherent complexity and heterogeneous structure of biological systems in combination with the measuring equipment. Moreover, a high degree of flexibility is required to treat rapidly changing experimental conditions. Object-oriented methodology seems to be well suited for this purpose. It enables an evolutionary approach to software development that still maintains a high degree of modularity. This paper presents experience with object-oriented technology gathered during several years of programming in the fields of bioprocess development and metabolic engineering. It concentrates on the aspects of experimental support, data analysis, interaction and visualization. Several examples are presented and discussed in the general context of the experimental cycle of knowledge acquisition, thus pointing out the benefits and problems of object-oriented technology in the specific application field of the biosciences. Finally, some strategies for future development are described.