Search

Abstract The conversion of CO2 into methanol depicts one of the most promising emerging renewable routes for the chemical and biotech industry. Under this regard, native methylotrophs have a large potential for converting methanol into value-added products but require targeted engineering approaches to enhance their performances and to widen their product spectrum. Here we use a systems-based approach to analyze and engineer M. extorquens TK 0001 for production of glycolic acid. Application of constraint-based metabolic modeling reveals the great potential of M. extorquens for that purpose, which is not yet described in literature. In particular, a superior theoretical product yield of 1.0 C-molGlycolic acid C-molMethanol −1 is predicted by our model, surpassing theoretical yields of sugar fermentation. Following this approach, we show here that strain engineering is viable and present 1st generation strains producing glycolic acid via a heterologous NADPH-dependent glyoxylate reductase. It was found that lactic acid is a surprising by-product of glycolic acid formation in M. extorquens, most likely due to a surplus of available NADH upon glycolic acid synthesis. Finally, the best performing strain was tested in a fed-batch fermentation producing a mixture of up to total 1.2 g L−1 glycolic acid and lactic acid. Several key performance indicators of our glycolic acid producer strain are superior to state-of-the-art synthetic methylotrophs. The presented results open the door for further strain engineering of the native methylotroph M. extorquens and pave the way to produce two promising biopolymer building blocks from green methanol, i.e., glycolic acid and lactic acid.

Abstract Background Due to increasing ecological concerns, microbial production of biochemicals from sustainable carbon sources like acetate is rapidly gaining importance. However, to successfully establish large-scale production scenarios, a solid understanding of metabolic driving forces is required to inform bioprocess design. To generate such knowledge, we constructed isopropanol-producing Escherichia coli W strains. Results Based on strain screening and metabolic considerations, a 2-stage process was designed, incorporating a growth phase followed by a nitrogen-starvation phase. This process design yielded the highest isopropanol titers on acetate to date (13.3 g L−1). Additionally, we performed shotgun and acetylated proteomics, and identified several stress conditions in the bioreactor scenarios, such as acid stress and impaired sulfur uptake. Metabolic modeling allowed for an in-depth characterization of intracellular flux distributions, uncovering cellular demand for ATP and acetyl-CoA as limiting factors for routing carbon toward the isopropanol pathway. Moreover, we asserted the importance of a balance between fluxes of the NADPH-providing isocitrate dehydrogenase (ICDH) and the product pathway. Conclusions Using the newly gained system-level understanding for isopropanol production from acetate, we assessed possible engineering approaches and propose process designs to maximize production. Collectively, our work contributes to the establishment and optimization of acetate-based bioproduction systems. Graphical Abstract

Abstract Anaerobic microbial fermentations provide high product yields and are a cornerstone of industrial bio-based processes. However, the need for redox balancing limits the array of fermentable substrate-product combinations. To overcome this limitation, here we design an aerobic fermentative metabolism that allows the introduction of selected respiratory modules. These can use oxygen to re-balance otherwise unbalanced fermentations, hence achieving controlled respiro-fermentative growth. Following this design, we engineer and characterize an obligate fermentative Escherichia coli strain that aerobically ferments glucose to stoichiometric amounts of lactate. We then re-integrate the quinone-dependent glycerol 3-phosphate dehydrogenase and demonstrate glycerol fermentation to lactate while selectively transferring the surplus of electrons to the respiratory chain. To showcase the potential of this fermentation mode, we direct fermentative flux from glycerol towards isobutanol production. In summary, our design permits using oxygen to selectively re-balance fermentations. This concept is an advance freeing highly efficient microbial fermentation from the limitations imposed by traditional redox balancing.

Abstract Background Zymomonas mobilis is well known for its outstanding ability to produce ethanol with both high specific productivity and with high yield close to the theoretical maximum. The key enzyme in the ethanol production pathway is the pyruvate decarboxylase (PDC) which is converting pyruvate to acetaldehyde. Since it is widely considered that its gene pdc is essential, metabolic engineering strategies aiming to produce other compounds derived from pyruvate need to find ways to reduce PDC activity. Results Here, we present a new platform strain (sGB027) of Z. mobilis in which the native promoter of pdc was replaced with the IPTG-inducible PT7A1, allowing for a controllable expression of pdc. Expression of lactate dehydrogenase from E. coli in sGB027 allowed the production of D-lactate with, to the best of our knowledge, the highest reported specific productivity of any microbial lactate producer as well as with the highest reported lactate yield for Z. mobilis so far. Additionally, by expressing the L-alanine dehydrogenase of Geobacillus stearothermophilus in sGB027 we produced L-alanine, further demonstrating the potential of sGB027 as a base for the production of compounds other than ethanol. Conclusion We demonstrated that our new platform strain can be an excellent starting point for the efficient production of various compounds derived from pyruvate with Z. mobilis and can thus enhance the establishment of this organism as a workhorse for biotechnological production processes.

Abstract Background The microbial production of isobutanol holds promise to become a sustainable alternative to fossil-based synthesis routes for this important chemical. Escherichia coli has been considered as one production host, however, due to redox imbalance, growth-coupled anaerobic production of isobutanol from glucose in E. coli is only possible if complex media additives or small amounts of oxygen are provided. These strategies have a negative impact on product yield, productivity, reproducibility, and production costs. Results In this study, we propose a strategy based on acetate as co-substrate for resolving the redox imbalance. We constructed the E. coli background strain SB001 (ΔldhA ΔfrdA ΔpflB) with blocked pathways from glucose to alternative fermentation products but with an enabled pathway for acetate uptake and subsequent conversion to ethanol via acetyl-CoA. This strain, if equipped with the isobutanol production plasmid pIBA4, showed robust exponential growth (µ = 0.05 h−1) under anaerobic conditions in minimal glucose medium supplemented with small amounts of acetate. In small-scale batch cultivations, the strain reached a glucose uptake rate of 4.8 mmol gDW−1 h−1, a titer of 74 mM and 89% of the theoretical maximal isobutanol/glucose yield, while secreting only small amounts of ethanol synthesized from acetate. Furthermore, we show that the strain keeps a high metabolic activity also in a pulsed fed-batch bioreactor cultivation, even if cell growth is impaired by the accumulation of isobutanol in the medium. Conclusions This study showcases the beneficial utilization of acetate as a co-substrate and redox sink to facilitate growth-coupled production of isobutanol under anaerobic conditions. This approach holds potential for other applications with different production hosts and/or substrate–product combinations.

Abstract The ubiquitous coexistence of the redox cofactors NADH and NADPH is widely considered to facilitate an efficient operation of cellular redox metabolism. However, it remains unclear what shapes the NAD(P)H specificity of specific redox reactions. Here, we present a computational framework to analyze the effect of redox cofactor swaps on the maximal thermodynamic potential of a metabolic network and use it to investigate key aspects of redox cofactor redundancy in Escherichia coli. As one major result, our analysis suggests that evolved NAD(P)H specificities are largely shaped by metabolic network structure and associated thermodynamic constraints enabling thermodynamic driving forces that are close or even identical to the theoretical optimum and significantly higher compared to random specificities. Furthermore, while redundancy of NAD(P)H is clearly beneficial for thermodynamic driving forces, a third redox cofactor would require a low standard redox potential to be advantageous. Our approach also predicts trends of redox-cofactor concentration ratios and could facilitate the design of optimal redox cofactor specificities.

Abstract One long‐standing question in microbiology is how microbes buffer perturbations in energy metabolism. In this study, we systematically analyzed the impact of different levels of ATP demand in Escherichia coli under various conditions (aerobic and anaerobic, with and without cell growth). One key finding is that, under all conditions tested, the glucose uptake increases with rising ATP demand, but only to a critical level beyond which it drops markedly, even below wild‐type levels. Focusing on anaerobic growth and using metabolomics and proteomics data in combination with a kinetic model, we show that this biphasic behavior is induced by the dual dependency of the phosphofructokinase on ATP (substrate) and ADP (allosteric activator). This mechanism buffers increased ATP demands by a higher glycolytic flux but, as shown herein, it collapses under very low ATP concentrations. Model analysis also revealed two major rate‐controlling steps in the glycolysis under high ATP demand, which could be confirmed experimentally. Our results provide new insights on fundamental mechanisms of bacterial energy metabolism and guide the rational engineering of highly productive cell factories.

Abstract Background The alcohol 2,3-butanediol (2,3-BDO) is an important chemical and an Escherichia coli producer strain was recently engineered for bio-based production of 2,3-BDO. However, further improvements are required for realistic applications. Results Here we report that enforced ATP wasting, implemented by overexpressing the genes of the ATP-hydrolyzing F1-part of the ATPase, leads to significant increases of yield and especially of productivity of 2,3-BDO synthesis in an E. coli producer strain under various cultivation conditions. We studied aerobic and microaerobic conditions as well as growth-coupled and growth-decoupled production scenarios. In all these cases, the specific substrate uptake and 2,3-BDO synthesis rate (up to sixfold and tenfold higher, respectively) were markedly improved in the ATPase strain compared to a control strain. However, aerobic conditions generally enable higher productivities only with reduced 2,3-BDO yields while high product yields under microaerobic conditions are accompanied with low productivities. Based on these findings we finally designed and validated a three-stage process for optimal conversion of glucose to 2,3-BDO, which enables a high productivity in combination with relatively high yield. The ATPase strain showed again superior performance and finished the process twice as fast as the control strain and with higher 2,3-BDO yield. Conclusions Our results demonstrate the high potential of enforced ATP wasting as a generic metabolic engineering strategy and we expect more applications to come in the future.

Abstract Background The concept of minimal cut sets (MCS) has become an important mathematical framework for analyzing and (re)designing metabolic networks. However, the calculation of MCS in genome-scale metabolic models is a complex computational problem. The development of duality-based algorithms in the last years allowed the enumeration of thousands of MCS in genome-scale networks by solving mixed-integer linear problems (MILP). A recent advancement in this field was the introduction of the MCS2 approach. In contrast to the Farkas-lemma-based dual system used in earlier studies, the MCS2 approach employs a more condensed representation of the dual system based on the nullspace of the stoichiometric matrix, which, due to its reduced dimension, holds promise to further enhance MCS computations. Results In this work, we introduce several new variants and modifications of duality-based MCS algorithms and benchmark their effects on the overall performance. As one major result, we generalize the original MCS2 approach (which was limited to blocking the operation of certain target reactions) to the most general case of MCS computations with arbitrary target and desired regions. Building upon these developments, we introduce a new MILP variant which allows maximal flexibility in the formulation of MCS problems and fully leverages the reduced size of the nullspace-based dual system. With a comprehensive set of benchmarks, we show that the MILP with the nullspace-based dual system outperforms the MILP with the Farkas-lemma-based dual system speeding up MCS computation with an averaged factor of approximately 2.5. We furthermore present several simplifications in the formulation of constraints, mainly related to binary variables, which further enhance the performance of MCS-related MILP. However, the benchmarks also reveal that some highly condensed formulations of constraints, especially on reversible reactions, may lead to worse behavior when compared to variants with a larger number of (more explicit) constraints and involved variables. Conclusions Our results further enhance the algorithmic toolbox for MCS calculations and are of general importance for theoretical developments as well as for practical applications of the MCS framework.

Abstract Background Enforced ATP wasting has been recognized as a promising metabolic engineering strategy to enhance the microbial production of metabolites that are coupled to ATP generation. It also appears to be a suitable approach to improve production of ethanol by Saccharomyces cerevisiae. In the present study, we constructed different S. cerevisiae strains with heterologous expression of genes of the ATP-hydrolyzing F1-part of the ATPase enzyme to induce enforced ATP wasting and quantify the resulting effect on biomass and ethanol formation. Results In contrast to genomic integration, we found that episomal expression of the αβγ subunits of the F1-ATPase genes of Escherichia coli in S. cerevisiae resulted in significantly increased ATPase activity, while neither genomic integration nor episomal expression of the β subunit from Trichoderma reesei could enhance ATPase activity. When grown in minimal medium under anaerobic growth-coupled conditions, the strains expressing E. coli’s F1-ATPase genes showed significantly improved ethanol yield (increase of 10% compared to the control strain). However, elevated product formation reduces biomass formation and, therefore, volumetric productivity. We demonstrate that this negative effect can be overcome under growth-decoupled (nitrogen-starved) operation with high and constant biomass concentration. Under these conditions, which mimic the second (production) phase of a two-stage fermentation process, the ATPase-expressing strains showed significant improvement in volumetric productivity (up to 111%) compared to the control strain. Conclusions Our study shows that expression of genes of the F1-portion of E. coli’s ATPase induces ATPase activity in S. cerevisiae and can be a promising way to improve ethanol production. This ATP-wasting strategy can be easily applied to other metabolites of interest, whose formation is coupled to ATP generation.

Abstract Systems biology has experienced dramatic growth in the number, size, and complexity of computational models. To reproduce simulation results and reuse models, researchers must exchange unambiguous model descriptions. We review the latest edition of the Systems Biology Markup Language (SBML), a format designed for this purpose. A community of modelers and software authors developed SBML Level 3 over the past decade. Its modular form consists of a core suited to representing reaction‐based models and packages that extend the core with features suited to other model types including constraint‐based models, reaction‐diffusion models, logical network models, and rule‐based models. The format leverages two decades of SBML and a rich software ecosystem that transformed how systems biologists build and interact with models. More recently, the rise of multiscale models of whole cells and organs, and new data sources such as single‐cell measurements and live imaging, has precipitated new ways of integrating data with models. We provide our perspectives on the challenges presented by these developments and how SBML Level 3 provides the foundation needed to support this evolution.

Abstract Background In order to improve the accuracy of constraint-based metabolic models, several approaches have been developed which intend to integrate additional biological information. Two of these methods, MOMENT and GECKO, incorporate enzymatic (kcat) parameters and enzyme mass constraints to further constrain the space of feasible metabolic flux distributions. While both methods have been proven to deliver useful extensions of metabolic models, they may considerably increase size and complexity of the models and there is currently no tool available to fully automate generation and calibration of such enzyme-constrained models from given stoichiometric models. Results In this work we present three major developments. We first conceived short MOMENT (sMOMENT), a simplified version of the MOMENT approach, which yields the same predictions as MOMENT but requires significantly fewer variables and enables direct inclusion of the relevant enzyme constraints in the standard representation of a constraint-based model. When measurements of enzyme concentrations are available, these can be included as well leading in the extreme case, where all enzyme concentrations are known, to a model representation that is analogous to the GECKO approach. Second, we developed the AutoPACMEN toolbox which allows an almost fully automated creation of sMOMENT-enhanced stoichiometric metabolic models. In particular, this includes the automatic read-out and processing of relevant enzymatic data from different databases and the reconfiguration of the stoichiometric model with embedded enzymatic constraints. Additionally, tools have been developed to adjust (kcat and enzyme pool) parameters of sMOMENT models based on given flux data. We finally applied the new sMOMENT approach and the AutoPACMEN toolbox to generate an enzyme-constrained version of the E. coli genome-scale model iJO1366 and analyze its key properties and differences with the standard model. In particular, we show that the enzyme constraints improve flux predictions (e.g., explaining overflow metabolism and other metabolic switches) and demonstrate, for the first time, that these constraints can markedly change the spectrum of metabolic engineering strategies for different target products. Conclusions The methodological and tool developments presented herein pave the way for a simplified and routine construction and analysis of enzyme-constrained metabolic models.

Microbial communities have become a major research focus due to their importance for biogeochemical cycles, biomedicine and biotechnological applications. While some biotechnological applications, such as anaerobic digestion, make use of naturally arising microbial communities, the rational design of microbial consortia for bio-based production processes has recently gained much interest. One class of synthetic microbial consortia is based on specifically designed strains of one species. A common design principle for these consortia is based on division of labor, where the entire production pathway is divided between the different strains to reduce the metabolic burden caused by product synthesis. We first show that classical division of labor does not automatically reduce the metabolic burden when metabolic flux per biomass is analyzed. We then present ASTHERISC (Algorithmic Search of THERmodynamic advantages in Single-species Communities), a new computational approach for designing multi-strain communities of a single-species with the aim to divide a production pathway between different strains such that the thermodynamic driving force for product synthesis is maximized. ASTHERISC exploits the fact that compartmentalization of segments of a product pathway in different strains can circumvent thermodynamic bottlenecks arising when operation of one reaction requires a metabolite with high and operation of another reaction the same metabolite with low concentration. We implemented the ASTHERISC algorithm in a dedicated program package and applied it on E. coli core and genome-scale models with different settings, for example, regarding number of strains or demanded product yield. These calculations showed that, for each scenario, many target metabolites (products) exist where a multi-strain community can provide a thermodynamic advantage compared to a single strain solution. In some cases, a production with sufficiently high yield is thermodynamically only feasible with a community. In summary, the developed ASTHERISC approach provides a promising new principle for designing microbial communities for the bio-based production of chemicals.

Flux balance analysis (FBA) is a key method for the constraint-based analysis of metabolic networks. A technical problem may occur in FBA when known (e.g., measured) fluxes of certain reactions are integrated into an FBA scenario rendering the underlying linear program (LP) infeasible, for example, due to inconsistencies between some of the measured fluxes causing a violation of the steady-state or other constraints. Here, we present and compare two methods, one based on an LP and one on a quadratic program (QP), to find minimal corrections for the given flux values so that the FBA problem becomes feasible. We provide a general guide on how to treat infeasible FBA systems in practice and discuss relevant examples of potentially infeasible scenarios in core and genome-scale metabolic models. Finally, we also highlight and clarify the relationships to classical metabolic flux analysis, where solely algebraic approaches are used to compute unknown metabolic rates from measured fluxes and to balance infeasible flux scenarios.

The concept of minimal cut sets (MCS) provides a flexible framework for analyzing properties of metabolic networks and for computing metabolic intervention strategies. In particular, it has been used to support the targeted design of microbial strains for bio-based production processes. Herein we present a number of major extensions that generalize the existing MCS approach and broaden its scope for applications in metabolic engineering. We first introduce a modified approach to integrate gene-protein-reaction associations (GPR) in the metabolic network structure for the computation of gene-based intervention strategies. In particular, we present a set of novel compression rules for GPR associations, which effectively speedup the computation of gene-based MCS by a factor of up to one order of magnitude. These rules are not specific for MCS and as well applicable to other computational strain design methods. Second, we enhance the MCS framework by allowing the definition of multiple target (undesired) and multiple protected (desired) regions. This enables precise tailoring of the metabolic solution space of the designed strain with unlimited flexibility. Together with further generalizations such as individual cost factors for each intervention, direct combinations of reaction/gene deletions and additions as well as the possibility to search for substrate co-feeding strategies, the scope of the MCS framework could be broadly extended. We demonstrate the applicability and performance benefits of the described developments by computing (gene-based) Escherichia coli strain designs for the bio-based production of 2,3-butanediol, a chemical, that has recently received much attention in the field of metabolic engineering. With our extended framework, we could identify promising strain designs that were formerly unpredictable, including those based on substrate co-feeding.

Interventions in metabolic networks can improve microbes for chemical production. Here, the authors develop a tool to identify metabolic valves that can decouple growth and production to systematically improve the rate and yield of biochemical production processes.

Coupling of growth and product synthesis is an important principle in metabolic engineering, but its range of applicability is unclear. Here, the authors use a dedicated computational framework to study the feasibility of coupling the production of metabolites to growth in the genome-scale metabolic models of five production organisms, and show that coupling can be achieved for most metabolites.

The process of anaerobic digestion in which waste biomass is transformed to methane by complex microbial communities has been modeled for more than 16 years by parametric gray box approaches that simplify process biology and do not resolve intracellular microbial activity. Information on such activity, however, has become available in unprecedented detail by recent experimental advances in metatranscriptomics and metaproteomics. The inclusion of such data could lead to more powerful process models of anaerobic digestion that more faithfully represent the activity of microbial communities. We augmented the Anaerobic Digestion Model No. 1 (ADM1) as the standard kinetic model of anaerobic digestion by coupling it to Flux-Balance-Analysis (FBA) models of methanogenic species. Steady-state results of coupled models are comparable to standard ADM1 simulations if the energy demand for non-growth associated maintenance (NGAM) is chosen adequately. When changing a constant feed of maize silage from continuous to pulsed feeding, the final average methane production remains very similar for both standard and coupled models, while both the initial response of the methanogenic population at the onset of pulsed feeding as well as its dynamics between pulses deviates considerably. In contrast to ADM1, the coupled models deliver predictions of up to 1,000s of intracellular metabolic fluxes per species, describing intracellular metabolic pathway activity in much higher detail. Furthermore, yield coefficients which need to be specified in ADM1 are no longer required as they are implicitly encoded in the topology of the species’ metabolic network. We show the feasibility of augmenting ADM1, an ordinary differential equation-based model for simulating biogas production, by FBA models implementing individual steps of anaerobic digestion. While cellular maintenance is introduced as a new parameter, the total number of parameters is reduced as yield coefficients no longer need to be specified. The coupled models provide detailed predictions on intracellular activity of microbial species which are compatible with experimental data on enzyme synthesis activity or abundance as obtained by metatranscriptomics or metaproteomics. By providing predictions of intracellular fluxes of individual community members, the presented approach advances the simulation of microbial community driven processes and provides a direct link to validation by state-of-the-art experimental techniques.

Constraint-based modeling (CBM) is increasingly used to analyze the metabolism of complex microbial communities involved in ecology, biomedicine, and various biotechnological processes. While CBM is an established framework for studying the metabolism of single species with linear stoichiometric models, CBM of communities with balanced growth is more complicated, not only due to the larger size of the multi-species metabolic network but also because of the bilinear nature of the resulting community models. Moreover, the solution space of these community models often contains biologically unrealistic solutions, which, even with model linearization and under application of certain objective functions, cannot easily be excluded. Here we present RedCom, a new approach to build reduced community models in which the metabolisms of the participating organisms are represented by net conversions computed from the respective single-species networks. By discarding (single-species) net conversions that violate a minimality criterion in the exchange fluxes, it is ensured that unrealistic solutions in the community model are excluded where a species altruistically synthesizes large amounts of byproducts (instead of biomass) to fulfill the requirements of other species. We employed the RedCom approach for modeling communities of up to nine organisms involved in typical degradation steps of anaerobic digestion in biogas plants. Compared to full (bilinear and linearized) community models, we found that the reduced community models obtained with RedCom are not only much smaller but allow, also in the largest model with nine species, extensive calculations required to fully characterize the solution space and to reveal key properties of communities with maximum methane yield and production rates. Furthermore, the predictive power of the reduced community models is significantly larger because they predict much smaller ranges of feasible community compositions and exchange fluxes still being consistent with measurements obtained from enrichment cultures. For an enrichment culture for growth on ethanol, we also used metaproteomic data to further constrain the solution space of the community models. Both model and proteomic data indicated a dominance of acetoclastic methanogens (Methanosarcinales) and Desulfovibrionales being the least abundant group in this microbial community.

Constraint-based modeling techniques have become a standard tool for the in silico analysis of metabolic networks. To further improve their accuracy, recent methodological developments focused on integration of thermodynamic information in metabolic models to assess the feasibility of flux distributions by thermodynamic driving forces. Here we present OptMDFpathway, a method that extends the recently proposed framework of Max-min Driving Force (MDF) for thermodynamic pathway analysis. Given a metabolic network model, OptMDFpathway identifies both the optimal MDF for a desired phenotypic behavior as well as the respective pathway itself that supports the optimal driving force. OptMDFpathway is formulated as a mixed-integer linear program and is applicable to genome-scale metabolic networks. As an important theoretical result, we also show that there exists always at least one elementary mode in the network that reaches the maximal MDF. We employed our new approach to systematically identify all substrate-product combinations in Escherichia coli where product synthesis allows for concomitant net CO2 assimilation via thermodynamically feasible pathways. Although biomass synthesis cannot be coupled to net CO2 fixation in E. coli we found that as many as 145 of the 949 cytosolic carbon metabolites contained in the genome-scale model iJO1366 enable net CO2 incorporation along thermodynamically feasible pathways with glycerol as substrate and 34 with glucose. The most promising products in terms of carbon assimilation yield and thermodynamic driving forces are orotate, aspartate and the C4-metabolites of the tricarboxylic acid cycle. We also identified thermodynamic bottlenecks frequently limiting the maximal driving force of the CO2-fixing pathways. Our results indicate that heterotrophic organisms like E. coli hold a possibly underestimated potential for CO2 assimilation which may complement existing biotechnological approaches for capturing CO2. Furthermore, we envision that the developed OptMDFpathway approach can be used for many other applications within the framework of constrained-based modeling and for rational design of metabolic networks.

Abstract Liver regeneration is a tightly controlled process mainly achieved by proliferation of usually quiescent hepatocytes. The specific molecular mechanisms ensuring cell division only in response to proliferative signals such as hepatocyte growth factor (HGF) are not fully understood. Here, we combined quantitative time‐resolved analysis of primary mouse hepatocyte proliferation at the single cell and at the population level with mathematical modeling. We showed that numerous G1/S transition components are activated upon hepatocyte isolation whereas DNA replication only occurs upon additional HGF stimulation. In response to HGF, Cyclin:CDK complex formation was increased, p21 rather than p27 was regulated, and Rb expression was enhanced. Quantification of protein levels at the restriction point showed an excess of CDK2 over CDK4 and limiting amounts of the transcription factor E2F‐1. Analysis with our mathematical model revealed that T160 phosphorylation of CDK2 correlated best with growth factor‐dependent proliferation, which we validated experimentally on both the population and the single cell level. In conclusion, we identified CDK2 phosphorylation as a gate‐keeping mechanism to maintain hepatocyte quiescence in the absence of HGF.

Elementary flux modes (EFMs) emerged as a formal concept to describe metabolic pathways and have become an established tool for constraint-based modeling and metabolic network analysis. EFMs are characteristic (support-minimal) vectors of the flux cone that contains all feasible steady-state flux vectors of a given metabolic network. EFMs account for (homogeneous) linear constraints arising from reaction irreversibilities and the assumption of steady state; however, other (inhomogeneous) linear constraints, such as minimal and maximal reaction rates frequently used by other constraint-based techniques (such as flux balance analysis [FBA]), cannot be directly integrated. These additional constraints further restrict the space of feasible flux vectors and turn the flux cone into a general flux polyhedron in which the concept of EFMs is not directly applicable anymore. For this reason, there has been a conceptual gap between EFM-based (pathway) analysis methods and linear optimization (FBA) techniques, as they operate on different geometric objects. One approach to overcome these limitations was proposed ten years ago and is based on the concept of elementary flux vectors (EFVs). Only recently has the community started to recognize the potential of EFVs for metabolic network analysis. In fact, EFVs exactly represent the conceptual development required to generalize the idea of EFMs from flux cones to flux polyhedra. This work aims to present a concise theoretical and practical introduction to EFVs that is accessible to a broad audience. We highlight the close relationship between EFMs and EFVs and demonstrate that almost all applications of EFMs (in flux cones) are possible for EFVs (in flux polyhedra) as well. In fact, certain properties can only be studied with EFVs. Thus, we conclude that EFVs provide a powerful and unifying framework for constraint-based modeling of metabolic networks.

Abstract Large‐scale protein signalling networks are useful for exploring complex biochemical pathways but do not reveal how pathways respond to specific stimuli. Such specificity is critical for understanding disease and designing drugs. Here we describe a computational approach—implemented in the free CNO software—for turning signalling networks into logical models and calibrating the models against experimental data. When a literature‐derived network of 82 proteins covering the immediate‐early responses of human cells to seven cytokines was modelled, we found that training against experimental data dramatically increased predictive power, despite the crudeness of Boolean approximations, while significantly reducing the number of interactions. Thus, many interactions in literature‐derived networks do not appear to be functional in the liver cells from which we collected our data. At the same time, CNO identified several new interactions that improved the match of model to data. Although missing from the starting network, these interactions have literature support. Our approach, therefore, represents a means to generate predictive, cell‐type‐specific models of mammalian signalling from generic protein signalling networks.

Abstract Purple non‐sulfur bacteria (Rhodospirillaceae) have been extensively employed for studying principles of photosynthetic and respiratory electron transport phosphorylation and for investigating the regulation of gene expression in response to redox signals. Here, we use mathematical modeling to evaluate the steady‐state behavior of the electron transport chain (ETC) in these bacteria under different environmental conditions. Elementary‐modes analysis of a stoichiometric ETC model reveals nine operational modes. Most of them represent well‐known functional states, however, two modes constitute reverse electron flow under respiratory conditions, which has been barely considered so far. We further present and analyze a kinetic model of the ETC in which rate laws of electron transfer steps are based on redox potential differences. Our model reproduces well‐known phenomena of respiratory and photosynthetic operation of the ETC and also provides non‐intuitive predictions. As one key result, model simulations demonstrate a stronger reduction of ubiquinone when switching from high‐light to low‐light conditions. This result is parameter insensitive and supports the hypothesis that the redox state of ubiquinone is a suitable signal for controlling photosynthetic gene expression.

Signaling pathways are characterized by crosstalk, feedback and feedforward mechanisms giving rise to highly complex and cell-context specific signaling networks. Dissecting the underlying relations is crucial to predict the impact of targeted perturbations. However, a major challenge in identifying cell-context specific signaling networks is the enormous number of potentially possible interactions. Here, we report a novel hybrid mathematical modeling strategy to systematically unravel hepatocyte growth factor (HGF) stimulated phosphoinositide-3-kinase (PI3K) and mitogen activated protein kinase (MAPK) signaling, which critically contribute to liver regeneration. By combining time-resolved quantitative experimental data generated in primary mouse hepatocytes with interaction graph and ordinary differential equation modeling, we identify and experimentally validate a network structure that represents the experimental data best and indicates specific crosstalk mechanisms. Whereas the identified network is robust against single perturbations, combinatorial inhibition strategies are predicted that result in strong reduction of Akt and ERK activation. Thus, by capitalizing on the advantages of the two modeling approaches, we reduce the high combinatorial complexity and identify cell-context specific signaling networks.

One ultimate goal of metabolic network modeling is the rational redesign of biochemical networks to optimize the production of certain compounds by cellular systems. Although several constraint-based optimization techniques have been developed for this purpose, methods for systematic enumeration of intervention strategies in genome-scale metabolic networks are still lacking. In principle, Minimal Cut Sets (MCSs; inclusion-minimal combinations of reaction or gene deletions that lead to the fulfilment of a given intervention goal) provide an exhaustive enumeration approach. However, their disadvantage is the combinatorial explosion in larger networks and the requirement to compute first the elementary modes (EMs) which itself is impractical in genome-scale networks. We present MCSEnumerator, a new method for effective enumeration of the smallest MCSs (with fewest interventions) in genome-scale metabolic network models. For this we combine two approaches, namely (i) the mapping of MCSs to EMs in a dual network, and (ii) a modified algorithm by which shortest EMs can be effectively determined in large networks. In this way, we can identify the smallest MCSs by calculating the shortest EMs in the dual network. Realistic application examples demonstrate that our algorithm is able to list thousands of the most efficient intervention strategies in genome-scale networks for various intervention problems. For instance, for the first time we could enumerate all synthetic lethals in E.coli with combinations of up to 5 reactions. We also applied the new algorithm exemplarily to compute strain designs for growth-coupled synthesis of different products (ethanol, fumarate, serine) by E.coli. We found numerous new engineering strategies partially requiring less knockouts and guaranteeing higher product yields (even without the assumption of optimal growth) than reported previously. The strength of the presented approach is that smallest intervention strategies can be quickly calculated and screened with neither network size nor the number of required interventions posing major challenges.

Cross-referencing experimental data with our current knowledge of signaling network topologies is one central goal of mathematical modeling of cellular signal transduction networks. We present a new methodology for data-driven interrogation and training of signaling networks. While most published methods for signaling network inference operate on Bayesian, Boolean, or ODE models, our approach uses integer linear programming (ILP) on interaction graphs to encode constraints on the qualitative behavior of the nodes. These constraints are posed by the network topology and their formulation as ILP allows us to predict the possible qualitative changes (up, down, no effect) of the activation levels of the nodes for a given stimulus. We provide four basic operations to detect and remove inconsistencies between measurements and predicted behavior: (i) find a topology-consistent explanation for responses of signaling nodes measured in a stimulus-response experiment (if none exists, find the closest explanation); (ii) determine a minimal set of nodes that need to be corrected to make an inconsistent scenario consistent; (iii) determine the optimal subgraph of the given network topology which can best reflect measurements from a set of experimental scenarios; (iv) find possibly missing edges that would improve the consistency of the graph with respect to a set of experimental scenarios the most. We demonstrate the applicability of the proposed approach by interrogating a manually curated interaction graph model of EGFR/ErbB signaling against a library of high-throughput phosphoproteomic data measured in primary hepatocytes. Our methods detect interactions that are likely to be inactive in hepatocytes and provide suggestions for new interactions that, if included, would significantly improve the goodness of fit. Our framework is highly flexible and the underlying model requires only easily accessible biological knowledge. All related algorithms were implemented in a freely available toolbox SigNetTrainer making it an appealing approach for various applications.

The epidermal growth factor receptor (EGFR) signaling pathway is probably the best-studied receptor system in mammalian cells, and it also has become a popular example for employing mathematical modeling to cellular signaling networks. Dynamic models have the highest explanatory and predictive potential; however, the lack of kinetic information restricts current models of EGFR signaling to smaller sub-networks. This work aims to provide a large-scale qualitative model that comprises the main and also the side routes of EGFR/ErbB signaling and that still enables one to derive important functional properties and predictions. Using a recently introduced logical modeling framework, we first examined general topological properties and the qualitative stimulus-response behavior of the network. With species equivalence classes, we introduce a new technique for logical networks that reveals sets of nodes strongly coupled in their behavior. We also analyzed a model variant which explicitly accounts for uncertainties regarding the logical combination of signals in the model. The predictive power of this model is still high, indicating highly redundant sub-structures in the network. Finally, one key advance of this work is the introduction of new techniques for assessing high-throughput data with logical models (and their underlying interaction graph). By employing these techniques for phospho-proteomic data from primary hepatocytes and the HepG2 cell line, we demonstrate that our approach enables one to uncover inconsistencies between experimental results and our current qualitative knowledge and to generate new hypotheses and conclusions. Our results strongly suggest that the Rac/Cdc42 induced p38 and JNK cascades are independent of PI3K in both primary hepatocytes and HepG2. Furthermore, we detected that the activation of JNK in response to neuregulin follows a PI3K-dependent signaling pathway.

Cellular decisions are determined by complex molecular interaction networks. Large-scale signaling networks are currently being reconstructed, but the kinetic parameters and quantitative data that would allow for dynamic modeling are still scarce. Therefore, computational studies based upon the structure of these networks are of great interest. Here, a methodology relying on a logical formalism is applied to the functional analysis of the complex signaling network governing the activation of T cells via the T cell receptor, the CD4/CD8 co-receptors, and the accessory signaling receptor CD28. Our large-scale Boolean model, which comprises 94 nodes and 123 interactions and is based upon well-established qualitative knowledge from primary T cells, reveals important structural features (e.g., feedback loops and network-wide dependencies) and recapitulates the global behavior of this network for an array of published data on T cell activation in wild-type and knock-out conditions. More importantly, the model predicted unexpected signaling events after antibody-mediated perturbation of CD28 and after genetic knockout of the kinase Fyn that were subsequently experimentally validated. Finally, we show that the logical model reveals key elements and potential failure modes in network functioning and provides candidates for missing links. In summary, our large-scale logical model for T cell activation proved to be a promising in silico tool, and it inspires immunologists to ask new questions. We think that it holds valuable potential in foreseeing the effects of drugs and network modifications.