A generalized computational framework to streamline thermodynamics and kinetics analysis of metabolic pathways

Metabolic engineering is a critical biotechnological approach in addressing global energy and environment challenges. Most engineering efforts, however, consist of laborious and inefficient trial-and-error of target pathways, due in part to the lack of methodologies that can comprehensively assess pathway properties in thermodynamics and kinetics. Metabolic engineering can benefit from computational tools that evaluate feasibility, expense and stability of non-natural metabolic pathways. Such tools can also help us understand natural pathways and their regulation at systems level. Here we introduce a computational toolbox, PathParser, which, for the first time, integrates multiple important functions for pathway analysis including thermodynamics analysis, kinetics-based protein cost optimization and robustness analysis. Specifically, PathParser enables optimization of the driving force of a pathway by minimizing the Gibbs free energy of least thermodynamically favorable reaction. In addition, based on reaction thermodynamics and enzyme kinetics, it can compute the minimal enzyme protein cost that supports metabolic flux, and evaluate pathway stability and flux in response to enzyme concentration perturbations. In a demo analysis of the Calvin–Benson–Bassham cycle and photorespiration pathway in the model cyanobacterium Synechocystis PCC 6803, the computation results are corroborated by experimental proteomics data as well as metabolic engineering outcomes. This toolbox may have broad application in metabolic engineering and systems biology in other microbial systems.     

Thermodynamic analysis of trinitrotoluene biodegradation and mineralization pathways

Biodegradation of 2,4,6-trinitrotoluene (TNT) proceeds through several different metabolic pathways. However, the reaction steps which are considered rate-controlling have not been fully determined. Glycolysis and other biological pathways contain biochemical reactions which are acutely rate-limiting due to enzyme control. These rate-limiting steps also have large negative Gibbs free energy changes. Because xenobiotic compounds such as TNT can be used by biological systems as nitrogen, carbon, and energy sources, it is likely that their degradation pathways also contain acutely rate-limiting steps. Identification of these rate-controlling reactions will enhance and better direct genetic engineering techniques to increase specific enzyme levels.This article identifies likely rate-controlling steps (or sets of steps) in reported TNT biodegradation pathways by estimating the Gibbs free energy change for each step and for the overall pathways. The biological standard Gibbs free energy change of reaction was calculated for each pathway step using a group contribution method specifically tailored for biomolecules. The method was also applied to hypothetical "pathways" constructed to mineralize TNT using several different microorganisms. Pathways steps that have large negative Gibbs free energy changes are postulated to be potentially rate-controlling. The microorganisms which utilize degradation pathways with the largest overall (from TNT to citrate) negatiave Gibbs free energy changes were also determined. Such microorganisms can extract more energy from the starting substrate and are thus assumed to have a competitive advantage over other microorganisms. Results from this modeling-based research are consistent with much experimental work available in the literature. (c) 1996 John Wiley & Sons, Inc.     

Thermodynamic constraints shape the structure of carbon fixation pathways

Thermodynamics impose a major constraint on the structure of metabolic pathways. Here, we use carbon fixation pathways to demonstrate how thermodynamics shape the structure of pathways and determine the cellular resources they consume. We analyze the energetic profile of prototypical reactions and show that each reaction type displays a characteristic change in Gibbs energy. Specifically, although carbon fixation pathways display a considerable structural variability, they are all energetically constrained by two types of reactions: carboxylation and carboxyl reduction. In fact, all adenosine triphosphate (ATP) molecules consumed by carbon fixation pathways - with a single exception - are used, directly or indirectly, to power one of these unfavorable reactions. When an indirect coupling is employed, the energy released by ATP hydrolysis is used to establish another chemical bond with high energy of hydrolysis, e.g. a thioester. This bond is cleaved by a downstream enzyme to energize an unfavorable reaction. Notably, many pathways exhibit reduced ATP requirement as they couple unfavorable carboxylation or carboxyl reduction reactions to exergonic reactions other than ATP hydrolysis. In the most extreme example, the reductive acetyl coenzyme A (acetyl-CoA) pathway bypasses almost all ATP-consuming reactions. On the other hand, the reductive pentose phosphate pathway appears to be the least ATP-efficient because it is the only carbon fixation pathway that invests ATP in metabolic aims other than carboxylation and carboxyl reduction. Altogether, our analysis indicates that basic thermodynamic considerations accurately predict the resource investment required to support a metabolic pathway and further identifies biochemical mechanisms that can decrease this requirement.     

Quantitative analysis of energy metabolic pathways in MCF-7 breast cancer cells by selected reaction monitoring assay

To investigate the quantitative response of energy metabolic pathways in human MCF-7 breast cancer cells to hypoxia, glucose deprivation, and estradiol stimulation, we developed a targeted proteomics assay for accurate quantification of protein expression in glycolysis/gluconeogenesis, TCA cycle, and pentose phosphate pathways. Cell growth conditions were selected to roughly mimic the exposure of cells in the cancer tissue to the intermittent hypoxia, glucose deprivation, and hormonal stimulation. Targeted proteomics assay allowed for reproducible quantification of 76 proteins in four different growth conditions after 24 and 48 h of perturbation. Differential expression of a number of control and metabolic pathway proteins in response to the change of growth conditions was found. Elevated expression of the majority of glycolytic enzymes was observed in hypoxia. Cancer cells, as opposed to near-normal MCF-10A cells, exhibited significantly increased expression of key energy metabolic pathway enzymes (FBP1, IDH2, and G6PD) that are known to redirect cellular metabolism and increase carbon flux through the pentose phosphate pathway. Our quantitative proteomic protocol is based on a mass spectrometry-compatible acid-labile detergent and is described in detail. Optimized parameters of a multiplex selected reaction monitoring (SRM) assay for 76 proteins, 134 proteotypic peptides, and 401 transitions are included and can be downloaded and used with any SRM-compatible mass spectrometer. The presented workflow is an integrated tool for hypothesis-driven studies of mammalian cells as well as functional studies of proteins, and can greatly complement experimental methods in systems biology, metabolic engineering, and metabolic transformation of cancer cells.     

Analysis of metabolic networks using a pathway distance metric through linear programming

The solution of the shortest path problem in biochemical systems constitutes an important step for studies of their evolution. In this paper, a linear programming (LP) algorithm for calculating minimal pathway distances in metabolic networks is studied. Minimal pathway distances are identified as the smallest number of metabolic steps separating two enzymes in metabolic pathways. The algorithm deals effectively with circularity and reaction directionality. The applicability of the algorithm is illustrated by calculating the minimal pathway distances for Escherichia coli small molecule metabolism enzymes, and then considering their correlations with genome distance (distance separating two genes on a chromosome) and enzyme function (as characterised by enzyme commission number). The results illustrate the effectiveness of the LP model. In addition, the data confirm that propinquity of genes on the genome implies similarity in function (as determined by co-involvement in the same region of the metabolic network), but suggest that no correlation exists between pathway distance and enzyme function. These findings offer insight into the probable mechanism of pathway evolution.     

On the role of the axial ligand in heme-based catalysis of the peroxidase and P450 type

 The present commentary focusses on the role of the axial ligand in peroxidase- and P450-type catalysis. Based on molecular orbital calculations and the experimental evidence available, it is argued that the ligand of a heme-containing enzyme may be a factor in setting the relative chance, although not the intrinsic capability, of the enzyme to catalyse a specific type of heme-based reaction chemistry. The ligand can do so by influencing the electrophilicity, i.e. the redox potential of the high-valency iron-oxo complex, and also by influencing the energy barrier for a reaction pathway through delocalization of valence electrons along the axial ligands, thereby, in the case of a cysteinate but not a histidine axial ligand, stabilizing oxygen transfer pathways. Received and accepted: 7 May 1996     

New insights on trehalose: a multifunctional molecule

Trehalose is a nonreducing disaccharide in which the two glucose units are linked in an alpha,alpha-1,1-glycosidic linkage. This sugar is present in a wide variety of organisms, including bacteria, yeast, fungi, insects, invertebrates, and lower and higher plants, where it may serve as a source of energy and carbon. In yeast and plants, it may also serve as a signaling molecule to direct or control certain metabolic pathways or even to affect growth. In addition, it has been shown that trehalose can protect proteins and cellular membranes from inactivation or denaturation caused by a variety of stress conditions, including desiccation, dehydration, heat, cold, and oxidation. Finally, in mycobacteria and corynebacteria, trehalose is an integral component of various glycolipids that are important cell wall structures. There are now at least three different pathways described for the biosynthesis of trehalose. The best known and most widely distributed pathway involves the transfer of glucose from UDP-glucose (or GDP-glucose in some cases) to glucose 6-phosphate to form trehalose-6-phosphate and UDP. This reaction is catalyzed by the trehalose-P synthase (TPS here, or OtsA in Escherichia coli ). Organisms that use this pathway usually also have a trehalose-P phosphatase (TPP here, or OtsB in E. coli) that converts the trehalose-P to free trehalose. A second pathway that has been reported in a few unusual bacteria involves the intramolecular rearrangement of maltose (glucosyl-alpha1,4-glucopyranoside) to convert the 1,4-linkage to the 1,1-bond of trehalose. This reaction is catalyzed by the enzyme called trehalose synthase and gives rise to free trehalose as the initial product. A third pathway involves several different enzymes, the first of which rearranges the glucose at the reducing end of a glycogen chain to convert the alpha1,4-linkage to an alpha,alpha1,1-bond. A second enzyme then releases the trehalose disaccharide from the reducing end of the glycogen molecule. Finally, in mushrooms there is a trehalose phosphorylase that catalyzes the phosphorolysis of trehalose to produce glucose-1-phosphate and glucose. This reaction is reversible in vitro and could theoretically give rise to trehalose from glucose-1-P and glucose. Another important enzyme in trehalose metabolism is trehalase (T), which may be involved in energy metabolism and also have a regulatory role in controlling the levels of trehalose in cells. This enzyme may be important in lowering trehalose concentrations once the stress is alleviated. Recent studies in yeast indicate that the enzymes involved in trehalose synthesis (TPS, TPP) exist together in a complex that is highly regulated at the activity level as well as at the genetic level.     

Sequential Activation of Metabolic Pathways: a Dynamic Optimization Approach

The regulation of cellular metabolism facilitates robust cellular operation in the face of changing external conditions. The cellular response to this varying environment may include the activation or inactivation of appropriate metabolic pathways. Experimental and numerical observations of sequential timing in pathway activation have been reported in the literature. It has been argued that such patterns can be rationalized by means of an underlying optimal metabolic design. In this paper we pose a dynamic optimization problem that accounts for time-resource minimization in pathway activation under constrained total enzyme abundance. The optimized variables are time-dependent enzyme concentrations that drive the pathway to a steady state characterized by a prescribed metabolic flux. The problem formulation addresses unbranched pathways with irreversible kinetics. Neither specific reaction kinetics nor fixed pathway length are assumed. In the optimal solution, each enzyme follows a switching profile between zero and maximum concentration, following a temporal sequence that matches the pathway topology. This result provides an analytic justification of the sequential activation previously described in the literature. In contrast with the existent numerical approaches, the activation sequence is proven to be optimal for a generic class of monomolecular kinetics. This class includes, but is not limited to, Mass Action, Michaelis–Menten, Hill, and some Power-law models. This suggests that sequential enzyme expression may be a common feature of metabolic regulation, as it is a robust property of optimal pathway activation.     

Catalytic pathways of Euphorbia characias peroxidase reacting with hydrogen peroxide

The reaction of Euphorbia characias latex peroxidase (ELP) with hydrogen peroxide as the sole substrate was studied by conventional and stopped-flow spectrophotometry. The reaction mechanism occurs via three distinct pathways. In the first (pathway I), ELP shows catalase-like activity: H2O2 oxidizes the native enzyme to compound I and subsequently acts as a reducing substrate, again converting compound I to the resting ferric enzyme. In the presence of an excess of hydrogen peroxide, compound I is still formed and further reacts in two other pathways. In pathway II, compound I initiates a series of cyclic reactions leading to the formation of compound II and compound III, and then returns to the native resting state. In pathway III, the enzyme is inactivated and compound I is converted into a bleached inactive species; this reaction proceeds faster in samples illuminated with bright white light, demonstrating that at least one of the intermediates is photosensitive. Calcium ions decrease the rate of pathway I and accelerate the rate of pathways II and III. Moreover, in the presence of calcium the inactive stable verdohemochrome P670 species accumulates. Thus, Ca2+ ions seem to be the key for all catalytic pathways of Euphorbia peroxidase.     

QM/MM investigation on the catalytic mechanism of Bacteroides thetaiotaomicron α-glucosidase BtGH97a

Bacteroides thetaiotaomicron α-glucosidase BtGH97a is an inverting enzyme. In this paper, the hydrolysis mechanism of p-nitro-phenyl α-d-glucopyranoside (pNP-Glc) catalyzed by BtGH97a was firstly studied by using quantum mechanical/molecular mechanical (QM/MM) approach. Two possible reaction pathways were considered. In the first pathway, a water molecule deprotonated by a nucleophilic base (here E439 or E508) attacks firstly on the anomeric carbon of pNP-Glc, then a proton from an acid residue (E532) attacks on the glycosidic oxygen to finish the hydrolysis reaction (named as nucleophilic attack-first pathway). In the second pathway, the proton from E532 attacks firstly on the glycosidic oxygen, then the water deprotonated by the nucleophilic base attacks on the anomeric carbon of pNP-Glc (named as proton attack-first pathway). Our calculation results indicate that the nucleophilic attack-first pathway is favorable in energy, in which the nucleophilic attack process is the rate-determining step with an energy barrier of 15.4 kcal/mol in the case of residue E508 as nucleophilic base. In this rate-determining step, the deprotonation of water and the attack on the anomeric carbon are concerted. In the proton attack-first pathway, the proton attack on the glycosidic oxygen is the rate-determining step, and the energy barrier is 24.1 kcal/mol. We conclude that the hydrolysis mechanism would follow nucleophilic attack-first pathway.     

Control in channelled pathways. A matrix method calculating the enzyme control coefficients

The usual equations expressing the enzyme control coefficients (quantitative indicators of 'global' control properties of a pathway) via the elasticity coefficients (reflecting local kinetic properties of an enzyme reaction), cannot be applied to a variety of 'non-ideal' pathways, in particular to pathways with metabolic channelling. Here we show that the relationship between the control and elasticity coefficients can be obtained by considering such a metabolic pathway as a network of elemental chemical conversions (steps). To calculate the control coefficients of enzymes one should first determine the elasticity coefficients of such elemental steps and then take their appropriate combinations. Although the method is illustrated for a channelled pathway it can be used for any non-ideal pathway including those with high enzyme concentrations where the sequestration of metabolites by enzymes cannot be neglected.     

Pathway Choice in Glutamate Synthesis in Escherichia coli

Escherichia coli has two primary pathways for glutamate synthesis. The glutamine synthetase-glutamate synthase (GOGAT) pathway is essential for synthesis at low ammonium concentration and for regulation of the glutamine pool. The glutamate dehydrogenase (GDH) pathway is important during glucose-limited growth. It has been hypothesized that GDH is favored when the organism is stressed for energy, because the enzyme does not use ATP as does the GOGAT pathway. The results of competition experiments between the wild-type and a GDH-deficient mutant during glucose-limited growth in the presence of the nonmetabolizable glucose analog α-methylglucoside were consistent with the hypothesis. Enzyme measurements showed that levels of the enzymes of the glutamate pathways dropped as the organism passed from unrestricted to glucose-restricted growth. However, other conditions influencing pathway choice had no substantial effect on enzyme levels. Therefore, substrate availability and/or modulation of enzyme activity are likely to be major determinants of pathway choice in glutamate synthesis.     

QM/MM investigation on the catalytic mechanism of Bacteroides thetaiotaomicron α-glucosidase BtGH97a

Bacteroides thetaiotaomicron α-glucosidase BtGH97a is an inverting enzyme. In this paper, the hydrolysis mechanism of p-nitro-phenyl α-d-glucopyranoside (pNP-Glc) catalyzed by BtGH97a was firstly studied by using quantum mechanical/molecular mechanical (QM/MM) approach. Two possible reaction pathways were considered. In the first pathway, a water molecule deprotonated by a nucleophilic base (here E439 or E508) attacks firstly on the anomeric carbon of pNP-Glc, then a proton from an acid residue (E532) attacks on the glycosidic oxygen to finish the hydrolysis reaction (named as nucleophilic attack-first pathway). In the second pathway, the proton from E532 attacks firstly on the glycosidic oxygen, then the water deprotonated by the nucleophilic base attacks on the anomeric carbon of pNP-Glc (named as proton attack-first pathway). Our calculation results indicate that the nucleophilic attack-first pathway is favorable in energy, in which the nucleophilic attack process is the rate-determining step with an energy barrier of 15.4kcal/mol in the case of residue E508 as nucleophilic base. In this rate-determining step, the deprotonation of water and the attack on the anomeric carbon are concerted. In the proton attack-first pathway, the proton attack on the glycosidic oxygen is the rate-determining step, and the energy barrier is 24.1kcal/mol. We conclude that the hydrolysis mechanism would follow nucleophilic attack-first pathway.     

Kinetic and thermodynamic principles determining the structural design of ATP-producing systems

It is theoretically analysed whether the structural design of ATP-producing pathways, in particular the design of glycolysis, may be explained by optimization principles. On the basis of kinetic and thermodynamic principles conclusions are derived concerning the stoichiometry of these pathways in states of high ATP production rates. One of the extensions to previous investigations is that the concentrations of the adenine nucleotides are taken into account as variable quantities. This necessitates the consideration of an interaction of the ATP-producing system I with an external ATP-consuming system II. A great variety of pathways is studied which differ in the number and location of ATP-consuming reactions, ATP-producing reactions and reactions involving inorganic phosphate. The corresponding number of possible pathways may be calculated in an explicit manner as a function of the number of those reactions which do not couple to ATP or inorganic phosphate. The kinetics of the individual reactions are described by linear or bilinear functions of reactant concentrations and all rate equations are expressed in terms of equilibrium constants and characteristic times. A thermodynamical analysis of the two coupled systems yields upper and lower limits for the concentration of ATP and an explicit expression for the maximal difference between the number of ATP-producing and ATP-consuming reactions of system I. The following results of the optimization are obtained. (i) The ATP production rate always increases if the ATP-producing reactions as well as those reactions characterized by an uptake of inorganic phosphate are shifted as far as possible towards the end of system I. (ii) Explicit conditions for the optimal location of the ATP-consuming reactions are presented. The results are discussed in the context of characteristic times as well as in terms of enzyme kinetic parameters. (iii) For two sets of characteristic times the resulting stoichiometries and their corresponding steady-state fluxes are investigated in detail. One of these stoichiometries shows a close correspondence to contemporary standard glycolysis. (iv) It is shown that most possible pathways result in a very low steady-state flux, that is, the optimal stoichiometry is characterized by a significant selective advantage. (v) The standard free energy profile of a pathway with an optimal stoichiometry is discussed. It differs significantly from the free energy profiles of nonoptimized pathways.     

Local kinetic measures of macromolecular structure reveal partitioning among multiple parallel pathways from the earliest steps in the folding of a large RNA molecule

At the heart of the RNA folding problem is the number, structures, and relationships among the intermediates that populate the folding pathways of most large RNA molecules. Unique insight into the structural dynamics of these intermediates can be gleaned from the time-dependent changes in local probes of macromolecular conformation (e.g. reports on individual nucleotide solvent accessibility offered by hydroxyl radical (()OH) footprinting). Local measures distributed around a macromolecule individually illuminate the ensemble of separate changes that constitute a folding reaction. Folding pathway reconstruction from a multitude of these individual measures is daunting due to the combinatorial explosion of possible kinetic models as the number of independent local measures increases. Fortunately, clustering of time progress curves sufficiently reduces the dimensionality of the data so as to make reconstruction computationally tractable. The most likely folding topology and intermediates can then be identified by exhaustively enumerating all possible kinetic models on a super-computer grid. The folding pathways and measures of the relative flux through them were determined for Mg(2+) and Na(+)-mediated folding of the Tetrahymena thermophila group I intron using this combined experimental and computational approach. The flux during Mg(2+)-mediated folding is divided among numerous parallel pathways. In contrast, the flux during the Na(+)-mediated reaction is predominantly restricted through three pathways, one of which is without detectable passage through intermediates. Under both conditions, the folding reaction is highly parallel with no single pathway accounting for more than 50% of the molecular flux. This suggests that RNA folding is non-sequential under a variety of different experimental conditions even at the earliest stages of folding. This study provides a template for the systematic analysis of the time-evolution of RNA structure from ensembles of local measures that will illuminate the chemical and physical characteristics of each step in the process. The applicability of this analysis approach to other macromolecules is discussed.     

Properties and Regulation of Microsomal PAF-Synthesizing Enzymes in Rat Brain Cortex

Platelet-activating factor (PAF) is a phospholipid mediator of long-term potentiation, synaptic plasticity and memory formation as well as of the development of brain damage. In brain, PAF is synthesized by two distinct pathways but their relative contribution to its productions, in various physiological and pathological conditions, is not established. We have further investigated on the properties of the two enzymes that catalyze the last step of the de novo or remodeling pathways in rat brain microsomes, PAF-synthesizing phosphocholinetransferase (PAF-PCT) and lysoPAF acetyltransferase (lysoPAF-AT), respectively. The latter enzyme is fully active at M Ca²⁺ concentration, inhibited by MgATP and activated by phosphorylation. Because the reversibility of the reaction catalyzed by PAF-PCT, its direction depends on the ratio [CDP-choline]/[CMP] which is related to the energy charge of the cell. These and other properties indicate that the de novo pathway should mainly contribute to PAF synthesis for maintaining its basal levels under physiological conditions. The remodeling pathway should be more involved in the production of PAF during ischemia. During reperfusion, the overproduction of PAF should be the result of the concomitant activation of both pathways.     

The alternative complement pathway revisited

Alternative pathway amplification plays a major role for the final effect of initial specific activation of the classical and lectin complement pathways, but the quantitative role of the amplification is insufficiently investigated. In experimental models of human diseases in which a direct activation of alternative pathway has been assumed, this interpretation needs revision placing a greater role on alternative amplification. We recently documented that the alternative amplification contributed to 80-90% of C5 activation when the initial activation was highly specific for the classical pathway. The recent identification of properdin as a recognition factor directly initiating alternative pathway activation, like C1q in the classical and mannose-binding lectin in the lectin pathway initiates a renewed interest in the reaction mechanisms of complement. Complement and Toll-like receptors, including the CD14 molecule, are two main upstream recognition systems of innate immunity, contributing to the inflammatory reaction in a number of conditions including ischemia-reperfusion injury and sepsis. These systems act as "double-edged swords", being protective against microbial invasion, but harmful to the host when activated improperly or uncontrolled. Combined inhibition of complement and Toll-like receptors/CD14 should be explored as a treatment regimen to reduce the overwhelming damaging inflammatory response during sepsis. The alternative pathway should be particularly considered in this regard, due to its uncontrolled amplification in sepsis. The alternative pathway should be regarded as a dual system, namely a recognition pathway principally similar to the classical and lectin pathways, and an amplification mechanism, well known, but quantitatively probably more important than generally recognized.     

Predicting microbial nitrogen pathways from basic principles

Nitrogen compounds are transformed by a complicated network of competing geochemical processes or microbial pathways, each performed by a different ecological guild of microorganisms. Complete experimental unravelling of this network requires a prohibitive experimental effort. Here we present a simple model that predicts relative rates of hypothetical nitrogen pathways, based only on the stoichiometry and energy yield of the performed redox reaction, assuming competition for resources between alternative pathways. Simulating competing pathways in hypothetical freshwater and marine sediment situations, we surprisingly found that much of the variation observed in nature can simply be predicted from these basic principles. Investigating discrepancies between observations and predictions led to two important biochemical factors that may create barriers for the viability of pathways: enzymatic costs for long pathways and high ammonium activation energy. We hypothesize that some discrepancies can be explained by non-equilibrium dynamics. The model predicted a pathway that has not been discovered in nature yet: the dismutation of nitrite to the level of nitrate and dinitrogen gas.     

Engineering cytosolic acetyl-coenzyme A supply in Saccharomyces cerevisiae: Pathway stoichiometry,free-energy conservation and redox-cofactor balancing

Saccharomyces cerevisiae is an important industrial cell factory and an attractive experimental model for evaluating novel metabolic engineering strategies. Many current and potential products of this yeast require acetyl coenzyme A (acetyl-CoA) as a precursor and pathways towards these products are generally expressed in its cytosol. The native S. cerevisiae pathway for production of cytosolic acetyl-CoA consumes 2 ATP equivalents in the acetyl-CoA synthetase reaction. Catabolism of additional sugar substrate, which may be required to generate this ATP, negatively affects product yields. Here, we review alternative pathways that can be engineered into yeast to optimize supply of cytosolic acetyl-CoA as a precursor for product formation. Particular attention is paid to reaction stoichiometry, free-energy conservation and redox-cofactor balancing of alternative pathways for acetyl-CoA synthesis from glucose. A theoretical analysis of maximally attainable yields on glucose of four compounds (n-butanol, citric acid, palmitic acid and farnesene) showed a strong product dependency of the optimal pathway configuration for acetyl-CoA synthesis. Moreover, this analysis showed that combination of different acetyl-CoA production pathways may be required to achieve optimal product yields. This review underlines that an integral analysis of energy coupling and redox-cofactor balancing in precursor-supply and product-formation pathways is crucial for the design of efficient cell factories.