Ensemble Modeling for Robustness Analysis in engineering non-native metabolic pathways

Engineering Escherichia coli for D-allulose biosynthesis from glycerol

D-allulose, a naturally occurring monosaccharide, is present in small quantities in nature. It is considered a valuable low-calorie sweetener due to its low absorption in the digestive tract and zero energy for growth. Most of the recent efforts to produce D-allulose have focused on in vitro enzyme catalysis. However, microbial fermentation is emerging as a promising alternative that offers the advantage of combining enzyme manufacturing and product synthesis within a single bioreactor. Here, a novel approach was proposed for the efficient biosynthesis of D-allulose from glycerol using metabolically engineered Escherichia coli. FbaA, Fbp, AlsE, and A6PP were used to construct the D-allulose synthesis pathway. Subsequently, PfkA, PfkB, and Pgi were disrupted to block the entry of the intermediate fructose-6-phosphate (F6P) into the Embden-Meyerhof-Parnas (EMP) and pentose phosphate (PP) pathways. Additionally, GalE and FryA were inactivated to reduce D-allulose consumption by the cells. Finally, a fed-batch fermentation process was implemented to optimize the performance of the cell factory. As a result, the titer of D-allulose reached 7.02 g/L with a maximum yield of 0.287 g/g.

Model-guided metabolic rewiring to bypass pyruvate oxidation for pyruvate derivative synthesis by minimizing carbon loss

Engineering microbial hosts to synthesize pyruvate derivatives depends on blocking pyruvate oxidation, thereby causing severe growth defects in aerobic glucose-based bioprocesses. To decouple pyruvate metabolism from cell growth to improve pyruvate availability, a genome-scale metabolic model combined with constraint-based flux balance analysis, geometric flux balance analysis, and flux variable analysis was used to identify genetic targets for strain design. Using translation elements from a ~3,000 cistronic library to modulate fxpK expression in a bicistronic cassette, a bifido shunt pathway was introduced to generate three molecules of non-pyruvate-derived acetyl-CoA from one molecule of glucose, bypassing pyruvate oxidation and carbon dioxide generation. The dynamic control of flux distribution by T7 RNAP-mediated synthetic small RNA decoupled pyruvate catabolism from cell growth. Adaptive laboratory evolution and multi-omics analysis revealed that a mutated isocitrate dehydrogenase functioned as a metabolic switch to activate the glyoxylate shunt as the only C4 anaplerotic pathway to generate malate from two molecules of acetyl-CoA input and bypass two decarboxylation reactions in the tricarboxylic acid cycle. A chassis strain for pyruvate derivative synthesis was constructed to reduce carbon loss by using the glyoxylate shunt as the only C4 anaplerotic pathway and the bifido shunt as a non-pyruvate-derived acetyl-CoA synthetic pathway and produced 22.46, 27.62, and 6.28 g/L of l-leucine, l-alanine, and l-valine by a controlled small RNA switch, respectively. Our study establishes a novel metabolic pattern of glucose-grown bacteria to minimize carbon loss under aerobic conditions and provides valuable insights into cell design for manufacturing pyruvate-derived products.IMPORTANCEBio-manufacturing from biomass-derived carbon sources using microbes as a cell factory provides an eco-friendly alternative to petrochemical-based processes. Pyruvate serves as a crucial building block for the biosynthesis of industrial chemicals; however, it is different to improve pyruvate availability in vivo due to the coupling of pyruvate-derived acetyl-CoA with microbial growth and energy metabolism via the oxidative tricarboxylic acid cycle. A genome-scale metabolic model combined with three algorithm analyses was used for strain design. Carbon metabolism was reprogrammed using two genetic control tools to fine-tune gene expression. Adaptive laboratory evolution and multi-omics analysis screened the growth-related regulatory targets beyond rational design. A novel metabolic pattern of glucose-grown bacteria is established to maintain growth fitness and minimize carbon loss under aerobic conditions for the synthesis of pyruvate-derived products. This study provides valuable insights into the design of a microbial cell factory for synthetic biology to produce industrial bio-products of interest.

The Calvin Benson cycle in bacteria: New insights from systems biology

The Calvin Benson cycle in phototrophic and chemolithoautotrophic bacteria has ecological and biotechnological importance, which has motivated study of its regulation. I review recent advances in our understanding of how the Calvin Benson cycle is regulated in bacteria and the technologies used to elucidate regulation and modify it, and highlight differences between and photoautotrophic and chemolithoautotrophic models. Systems biology studies have shown that in oxygenic phototrophic bacteria, Calvin Benson cycle enzymes are extensively regulated at post-transcriptional and post-translational levels, with multiple enzyme activities connected to cellular redox status through thioredoxin. In chemolithoautotrophic bacteria, regulation is primarily at the transcriptional level, with effector metabolites transducing cell status, though new methods should now allow facile, proteome-wide exploration of biochemical regulation in these models. A biotechnological objective is to enhance CO2 fixation in the cycle and partition that carbon to a product of interest. Flux control of CO2 fixation is distributed over multiple enzymes, and attempts to modulate gene Calvin cycle gene expression show a robust homeostatic regulation of growth rate, though the synthesis rates of products can be significantly increased. Therefore, de-regulation of cycle enzymes through protein engineering may be necessary to increase fluxes. Non-canonical Calvin Benson cycles, if implemented with synthetic biology, could have reduced energy demand and enzyme loading, thus increasing the attractiveness of these bacteria for industrial applications.

Creating new-to-nature carbon fixation: A guide

Synthetic biology aims at designing new biological functions from first principles. These new designs allow to expand the natural solution space and overcome the limitations of naturally evolved systems. One example is synthetic CO2-fixation pathways that promise to provide more efficient ways for the capture and conversion of CO2 than natural pathways, such as the Calvin Benson Bassham (CBB) cycle of photosynthesis. In this review, we provide a practical guideline for the design and realization of such new-to-nature CO2-fixation pathways. We introduce the concept of "synthetic CO2-fixation", and give a general overview over the enzymology and topology of synthetic pathways, before we derive general principles for their design from their eight naturally evolved analogs. We provide a comprehensive summary of synthetic carbon-assimilation pathways and derive a step-by-step, practical guide from the theoretical design to their practical implementation, before ending with an outlook on new developments in the field.

High-level and -yield production of L-leucine in engineered Escherichia coli by multistep metabolic engineering

L-leucine is an essential amino acid widely used in food and pharmaceutical industries. However, the relatively low production efficiency limits its large-scale application. In this study, we rationally developed an efficient L-leucine-producing Escherichia coli strain. Initially, the L-leucine synthesis pathway was enhanced by overexpressing feedback-resistant 2-isopropylmalate synthase and acetohydroxy acid synthase both derived from Corynebacterium glutamicum, along with two other native enzymes. Next, the pyruvate and acetyl-CoA pools were enriched by deleting competitive pathways, employing the nonoxidative glycolysis pathway, and dynamically modulating the citrate synthase activity, which significantly promoted the L-leucine production and yield to 40.69 g/L and 0.30 g/g glucose, respectively. Then, the redox flux was improved by substituting the native NADPH-dependent acetohydroxy acid isomeroreductase, branched chain amino acid transaminase, and glutamate dehydrogenase with their NADH-dependent equivalents. Finally, L-leucine efflux was accelerated by precise overexpression of the exporter and deletion of the transporter. Under fed-batch conditions, the final strain LXH-21 produced 63.29 g/L of L-leucine, with a yield and productivity of 0.37 g/g glucose and 2.64 g/(L h), respectively. To our knowledge, this study achieved the highest production efficiency of L-leucine to date. The strategies presented here will be useful for engineering E. coli strains for producing L-leucine and related products on an industrial scale.

Thermodynamic and Kinetic Modeling Directs Pathway Optimization for Isopropanol Production in a Gas-Fermenting Bacterium

Highly efficient bioproduction from gaseous substrates (e.g., hydrogen and carbon oxides) will require systematic optimization of the host microbes. To date, the rational redesign of gas-fermenting bacteria is still in its infancy, due in part to the lack of quantitative and precise metabolic knowledge that can direct strain engineering.

Computer-Based Design of a Cell Factory for High-Yield Cytidine Production

Pyrimidine ribonucleotide de novo biosynthesis pathway (PRdnBP) is an important pathway to produce pyrimidine nucleosides. We attempted to systematically investigate PRdnBP in Escherichia coli with genome-scale metabolic models and utilized the models to guide strain design. The balance of central carbon metabolism and PRdnBP affected the production of cytidine from glucose. Using Bayesian metabolic flux analysis, the effect of modified PRdnBP on the metabolic network was analyzed. The acetate overflow became coupled with PRdnBP flux, while they were originally independent under oxygen-sufficient conditions. The coupling between cytidine production and acetate secretion in the modified strain was weakened by arcA deletion, which resulted in further improving the efficient accumulation of cytidine. In total, 1.28 g/L of cytidine with a yield of 0.26 g/g glucose was produced. The yield of cytidine produced by E. coli is higher than previous reports. Our strategy provides an effective attempt to find metabolic bottlenecks in genetically engineered bacteria by using flux coupling analysis.

Reconstructing Kinetic Models for Dynamical Studies of Metabolism using Generative Adversarial Networks

Kinetic models of metabolism relate metabolic fluxes, metabolite concentrations and enzyme levels through mechanistic relations, rendering them essential for understanding, predicting and optimizing the behaviour of living organisms. However, due to the lack of kinetic data, traditional kinetic modelling often yields only a few or no kinetic models with desirable dynamical properties, making the analysis unreliable and computationally inefficient. We present REKINDLE (Reconstruction of Kinetic Models using Deep Learning), a deep-learning-based framework for efficiently generating kinetic models with dynamic properties matching the ones observed in cells. We showcase REKINDLE's capabilities to navigate through the physiological states of metabolism using small numbers of data with significantly lower computational requirements. The results show that data-driven neural networks assimilate implicit kinetic knowledge and structure of metabolic networks and generate kinetic models with tailored properties and statistical diversity. We anticipate that our framework will advance our understanding of metabolism and accelerate future research in biotechnology and health.

Toward Methanol-Based Biomanufacturing: Emerging Strategies for Engineering Synthetic Methylotrophy in Saccharomyces cerevisiae

The global expansion of biomanufacturing is currently limited by the availability of sugar-based microbial feedstocks, which require farmland for cultivation and therefore cannot support large increases in production without impacting the human food supply. One-carbon feedstocks, such as methanol, present an enticing alternative to sugar because they can be produced independently of arable farmland from organic waste, atmospheric carbon dioxide, and hydrocarbons such as biomethane, natural gas, and coal. The development of efficient industrial microorganisms that can convert one-carbon feedstocks into valuable products is an ongoing challenge. This review discusses progress in the field of synthetic methylotrophy with a focus on how it pertains to the important industrial yeast, Saccharomyces cerevisiae. Recent insights generated from engineering synthetic methylotrophic xylulose- and ribulose-monophosphate cycles, reductive glycine pathways, and adaptive laboratory evolution studies are critically assessed to generate novel strategies for the future engineering of methylotrophy in S. cerevisiae.

Construction of Multiscale Genome-Scale Metabolic Models: Frameworks and Challenges

Genome-scale metabolic models (GEMs) are effective tools for metabolic engineering and have been widely used to guide cell metabolic regulation. However, the single gene–protein-reaction data type in GEMs limits the understanding of biological complexity. As a result, multiscale models that add constraints or integrate omics data based on GEMs have been developed to more accurately predict phenotype from genotype. This review summarized the recent advances in the development of multiscale GEMs, including multiconstraint, multiomic, and whole-cell models, and outlined machine learning applications in GEM construction. This review focused on the frameworks, toolkits, and algorithms for constructing multiscale GEMs. The challenges and perspectives of multiscale GEM development are also discussed.

Metabolic Engineering Interventions for Sustainable 2,3-Butanediol Production in Gas-Fermenting Clostridium autoethanogenum

Gas fermentation provides a promising platform to turn low-cost and readily available single-carbon waste gases into commodity chemicals, such as 2,3-butanediol. Clostridium autoethanogenum is usually used as a robust and flexible chassis for gas fermentation. Here, we leveraged constraint-based stoichiometric modeling and kinetic ensemble modeling of the C. autoethanogenum metabolic network to provide a systematic in silico analysis of metabolic engineering interventions for 2,3-butanediol overproduction and low carbon substrate loss in dissipated CO₂. Our analysis allowed us to identify and to assess comparatively the expected performances for a wide range of single, double, and triple interventions. Our analysis managed to individuate bottleneck reactions in relevant metabolic pathways when suggesting intervening strategies. Besides recapitulating intuitive and/or previously attempted genetic modifications, our analysis neatly outlined that interventions-at least partially-impinging on by-products branching from acetyl coenzyme A (acetyl-CoA) and pyruvate (acetate, ethanol, amino acids) offer valuable alternatives to the interventions focusing directly on the specific branch from pyruvate to 2,3-butanediol. IMPORTANCE Envisioning value chains inspired by environmental sustainability and circularity in economic models is essential to counteract the alterations in the global natural carbon cycle induced by humans. Recycling carbon-based waste gas streams into chemicals by devising gas fermentation bioprocesses mediated by acetogens of the genus Clostridium is one component of the solution. Carbon monoxide originates from multiple biogenic and abiogenic sources and bears a significant environmental impact. This study aims at identifying metabolic engineering interventions for increasing 2,3-butanediol production and avoiding carbon loss in CO₂ dissipation via C. autoethanogenum fermenting a substrate comprising CO and H₂. 2,3-Butanediol is a valuable biochemical by-product since, due to its versatility, can be transformed quite easily into chemical compounds such as butadiene, diacetyl, acetoin, and methyl ethyl ketone. These compounds are usable as building blocks to manufacture a vast range of industrially produced chemicals.

Unravelling Formaldehyde Metabolism in Bacteria: Road towards Synthetic Methylotrophy

Formaldehyde metabolism is prevalent in all organisms, where the accumulation of formaldehyde can be prevented through the activity of dissimilation pathways. Furthermore, formaldehyde assimilatory pathways play a fundamental role in many methylotrophs, which are microorganisms able to build biomass and obtain energy from single- and multicarbon compounds with no carbon–carbon bonds. Here, we describe how formaldehyde is formed in the environment, the mechanisms of its toxicity to the cells, and the cell’s strategies to circumvent it. While their importance is unquestionable for cell survival in formaldehyde rich environments, we present examples of how the modification of native formaldehyde dissimilation pathways in nonmethylotrophic bacteria can be applied to redirect carbon flux toward heterologous, synthetic formaldehyde assimilation pathways introduced into their metabolism. Attempts to engineer methylotrophy into nonmethylotrophic hosts have gained interest in the past decade, with only limited successes leading to the creation of autonomous synthetic methylotrophy. Here, we discuss how native formaldehyde assimilation pathways can additionally be employed as a premise to achieving synthetic methylotrophy. Lastly, we discuss how emerging knowledge on regulation of formaldehyde metabolism can contribute to creating synthetic regulatory circuits applied in metabolic engineering strategies.

Assessing the impact of substrate-level enzyme regulations limiting ethanol titer in Clostridium thermocellum using a core kinetic model

Clostridium thermocellum is a promising candidate for consolidated bioprocessing because it can directly ferment cellulose to ethanol. Despite significant efforts, achieved yields and titers fall below industrially relevant targets. This implies that there still exist unknown enzymatic, regulatory, and/or possibly thermodynamic bottlenecks that can throttle back metabolic flow. By (i) elucidating internal metabolic fluxes in wild-type C. thermocellum grown on cellobiose via ¹³C-metabolic flux analysis (¹³C-MFA), (ii) parameterizing a core kinetic model, and (iii) subsequently deploying an ensemble-docking workflow for discovering substrate-level regulations, this paper aims to reveal some of these factors and expand our knowledgebase governing C. thermocellum metabolism. Generated ¹³C labeling data were used with ¹³C-MFA to generate a wild-type flux distribution for the metabolic network. Notably, flux elucidation through MFA alluded to serine generation via the mercaptopyruvate pathway. Using the elucidated flux distributions in conjunction with batch fermentation process yield data for various mutant strains, we constructed a kinetic model of C. thermocellum core metabolism (i.e. k-ctherm138). Subsequently, we used the parameterized kinetic model to explore the effect of removing substrate-level regulations on ethanol yield and titer. Upon exploring all possible simultaneous (up to four) regulation removals we identified combinations that lead to many-fold model predicted improvement in ethanol titer. In addition, by coupling a systematic method for identifying putative competitive inhibitory mechanisms using K-FIT kinetic parameterization with the ensemble-docking workflow, we flagged 67 putative substrate-level inhibition mechanisms across central carbon metabolism supported by both kinetic formalism and docking analysis.

In-Depth Computational Analysis of Natural and Artificial Carbon Fixation Pathways

In the recent years, engineering new-to-nature CO2- and C1-fixing metabolic pathways made a leap forward. New, artificial pathways promise higher yields and activity than natural ones like the Calvin-Benson-Bassham (CBB) cycle. The question remains how to best predict their in vivo performance and what actually makes one pathway "better" than another. In this context, we explore aerobic carbon fixation pathways by a computational approach and compare them based on their specific activity and yield on methanol, formate, and CO2/H2 considering the kinetics and thermodynamics of the reactions. Besides pathways found in nature or implemented in the laboratory, this included two completely new cycles with favorable features: the reductive citramalyl-CoA cycle and the 2-hydroxyglutarate-reverse tricarboxylic acid cycle. A comprehensive kinetic data set was collected for all enzymes of all pathways, and missing kinetic data were sampled with the Parameter Balancing algorithm. Kinetic and thermodynamic data were fed to the Enzyme Cost Minimization algorithm to check for respective inconsistencies and calculate pathway-specific activities. The specific activities of the reductive glycine pathway, the CETCH cycle, and the new reductive citramalyl-CoA cycle were predicted to match the best natural cycles with superior product-substrate yield. However, the CBB cycle performed better in terms of activity compared to the alternative pathways than previously thought. We make an argument that stoichiometric yield is likely not the most important design criterion of the CBB cycle. Still, alternative carbon fixation pathways were paretooptimal for specific activity and product-substrate yield in simulations with C1 substrates and CO2/H2 and therefore hold great potential for future applications in Industrial Biotechnology and Synthetic Biology.

Metabolome and proteome analyses reveal transcriptional misregulation in glycolysis of engineered E. coli

Synthetic metabolic pathways are a burden for engineered bacteria, but the underlying mechanisms often remain elusive. Here we show that the misregulated activity of the transcription factor Cra is responsible for the growth burden of glycerol overproducing E. coli. Glycerol production decreases the concentration of fructose-1,6-bisphoshate (FBP), which then activates Cra resulting in the downregulation of glycolytic enzymes and upregulation of gluconeogenesis enzymes. Because cells grow on glucose, the improper activation of gluconeogenesis and the concomitant inhibition of glycolysis likely impairs growth at higher induction of the glycerol pathway. We solve this misregulation by engineering a Cra-binding site in the promoter controlling the expression of the rate limiting enzyme of the glycerol pathway to maintain FBP levels sufficiently high. We show the broad applicability of this approach by engineering Cra-dependent regulation into a set of constitutive and inducible promoters, and use one of them to overproduce carotenoids in E. coli.

Synthetic pathways represent a metabolic burden on host cells. Here the authors engineer Cra-binding sites to prevent misregulation in glycerol and carotenoid overproducing E. coli strains.

Enhancing in silico strain design predictions through next generation metabolic modeling approaches

The reconstruction and analysis of metabolic models has garnered increasing attention due to the multitude of applications in which these have proven to be practical. The growing number of generated metabolic models has been accompanied by an exponentially expanding arsenal of tools used to analyze them. In this work, we discussed the biological relevance of a number of promising modeling frameworks, focusing on the questions and hypotheses each method is equipped to address. To this end, we critically analyzed the steady-state modeling approaches focusing on resource allocation and incorporation of thermodynamic considerations which produce promising results and aid in the generation and experimental validation of numerous predictions. For smaller networks involving more complex regulation, we addressed kinetic modeling techniques which show encouraging results in addressing questions outside the scope of steady-state modeling. Finally, we discussed the potential application of the discussed frameworks within the field of strain design. Adoption of such methodologies is believed to significantly enhance the accuracy of in silico predictions and hence decrease the number of design-build-test cycles required.

Bottom-Up Synthesis of Artificial Cells: Recent Highlights and Future Challenges

The bottom-up approach in synthetic biology aims to create molecular ensembles that reproduce the organization and functions of living organisms and strives to integrate them in a modular and hierarchical fashion toward the basic unit of life-the cell-and beyond. This young field stands on the shoulders of fundamental research in molecular biology and biochemistry, next to synthetic chemistry, and, augmented by an engineering framework, has seen tremendous progress in recent years thanks to multiple technological and scientific advancements. In this timely review of the research over the past decade, we focus on three essential features of living cells: the ability to self-reproduce via recursive cycles of growth and division, the harnessing of energy to drive cellular processes, and the assembly of metabolic pathways. In addition, we cover the increasing efforts to establish multicellular systems via different communication strategies and critically evaluate the potential applications.

Biotechnology for secure biocontainment designs in an emerging bioeconomy

Genetically modified organisms (GMOs) have emerged as an integral component of a sustainable bioeconomy, with an array of applications in agriculture, bioenergy, and biomedicine. However, the rapid development of GMOs and associated synthetic biology approaches raises a number of biosecurity concerns related to environmental escape of GMOs, detection thereof, and impact upon native ecosystems. A myriad of genetic safeguards have been deployed in diverse microbial hosts, ranging from classical auxotrophies to global genome recoding. However, to realize the full potential of microbes as biocatalytic platforms in the bioeconomy, a deeper understanding of the fundamental principles governing microbial responsiveness to biocontainment constraints, and interactivity of GMOs with the environment, is required. Herein, we review recent analytical biotechnological advances and strategies to assess biocontainment and microbial bioproductivity, as well as opportunities for predictive systems biodesigns towards securing a viable bioeconomy.

Role of cyanobacterial phosphoketolase in energy regulation and glucose secretion under dark anaerobic and osmotic stress conditions

Primary metabolism in cyanobacteria is built on the Calvin-Benson-Bassham (CBB) cycle, oxidative pentose phosphate (OPP) pathway, Embden-Meyerhof-Parnas (EMP) pathway, and the tricarboxylic acid (TCA) cycle. Phosphoketolase (Xpk), commonly found in cyanobacteria, is an enzyme that is linked to all these pathways. However, little is known about its physiological role. Here, we show that most of the cyanobacterial Xpk surveyed are inhibited by ATP, and both copies of Xpk in nitrogen-fixing Cyanothece ATCC51142 are further activated by ADP, suggesting their role in energy regulation. Moreover, Xpk in Synechococcus elongatus PCC7942 and Cyanothece ATCC51142 show that their expressions are dusk-peaked, suggesting their roles in dark conditions. Finally, we find that Xpk in S. elongatus PCC7942 is responsible for survival using ATP produced from the glycogen-to-acetate pathway under dark, anaerobic condition. Interestingly, under this condition, xpk deletion causes glucose secretion in response to osmotic shock such as NaHCO₃, KHCO₃ and NaCl as part of incomplete glycogen degradation. These findings unveiled the role of this widespread enzyme and open the possibility for enhanced glucose secretion from cyanobacteria.

Systems biology approaches integrated with artificial intelligence for optimized metabolic engineering

Metabolic engineering aims to maximize the production of bio-economically important substances (compounds, enzymes, or other proteins) through the optimization of the genetics, cellular processes and growth conditions of microorganisms. This requires detailed understanding of underlying metabolic pathways involved in the production of the targeted substances, and how the cellular processes or growth conditions are regulated by the engineering. To achieve this goal, a large system of experimental techniques, compound libraries, computational methods and data resources, including multi-omics data, are used. The recent advent of multi-omics systems biology approaches significantly impacted the field by opening new avenues to perform dynamic and large-scale analyses that deepen our knowledge on the manipulations. However, with the enormous transcriptomics, proteomics and metabolomics available, it is a daunting task to integrate the data for a more holistic understanding. Novel data mining and analytics approaches, including Artificial Intelligence (AI), can provide breakthroughs where traditional low-throughput experiment-alone methods cannot easily achieve. Here, we review the latest attempts of combining systems biology and AI in metabolic engineering research, and highlight how this alliance can help overcome the current challenges facing industrial biotechnology, especially for food-related substances and compounds using microorganisms.

Converting Escherichia coli to a Synthetic Methylotroph Growing Solely on Methanol

Methanol, being electron rich and derivable from methane or CO₂, is a potentially renewable one-carbon (C1) feedstock for microorganisms. Although the ribulose monophosphate (RuMP) cycle used by methylotrophs to assimilate methanol differs from the typical sugar metabolism by only three enzymes, turning a non-methylotrophic organism to a synthetic methylotroph that grows to a high cell density has been challenging. Here we reprogrammed E. coli using metabolic robustness criteria followed by laboratory evolution to establish a strain that can efficiently utilize methanol as the sole carbon source. This synthetic methylotroph alleviated a so far uncharacterized hurdle, DNA-protein crosslinking (DPC), by insertion sequence (IS)-mediated copy number variations (CNVs) and balanced the metabolic flux by mutations. Being capable of growing at a rate comparable with natural methylotrophs in a wide range of methanol concentrations, this synthetic methylotrophic strain illustrates genome editing and evolution for microbial tropism changes and expands the scope of biological C1 conversion.

Mechanistic Modeling of Biochemical Systems without A Priori Parameter Values Using the Design Space Toolbox v.3.0

Mechanistic models of biochemical systems provide a rigorous description of biological phenomena. They are indispensable for making predictions and elucidating biological design principles. To date, mathematical analysis and characterization of these models encounter a bottleneck consisting of large numbers of unknown parameter values. Here, we introduce the Design Space Toolbox v.3.0 (DST3), a software implementation of the Design Space formalism enabling mechanistic modeling without requiring previous knowledge of parameter values. This is achieved by using a phenotype-centric modeling approach, in which the system is first decomposed into a series of biochemical phenotypes. Parameter values realizing phenotypes of interest are subsequently predicted. DST3 represents the most generally applicable implementation of the Design Space formalism and offers unique advantages over earlier versions. By expanding the Design Space formalism and streamlining its distribution, DST3 represents a valuable tool for elucidating biological design principles and designing novel synthetic circuits.

•

DST3 extends the Design Space formalism by identifying additional phenotypes

•

Additional phenotypes arise from cycles, conservations, and metabolic imbalances

•

DST3 enables mechanistic modeling without previous knowledge of parameter values

•

It fully unlocks the potential of the novel phenotype-centric modeling strategy

Understanding biochemical design principles with ensembles of canonical non-linear models

Systems biology applies concepts from engineering in order to understand biological networks. If such an understanding was complete, biologists would be able to design ad hoc biochemical components tailored for different purposes, which is the goal of synthetic biology. Needless to say that we are far away from creating biological subsystems as intricate and precise as those found in nature, but mathematical models and high throughput techniques have brought us a long way in this direction. One of the difficulties that still needs to be overcome is finding the right values for model parameters and dealing with uncertainty, which is proving to be an extremely difficult task. In this work, we take advantage of ensemble modeling techniques, where a large number of models with different parameter values are formulated and then tested according to some performance criteria. By finding features shared by successful models, the role of different components and the synergies between them can be better understood. We will address some of the difficulties often faced by ensemble modeling approaches, such as the need to sample a space whose size grows exponentially with the number of parameters, and establishing useful selection criteria. Some methods will be shown to reduce the predictions from many models into a set of understandable “design principles” that can guide us to improve or manufacture a biochemical network. Our proposed framework formulates models within standard formalisms in order to integrate information from different sources and minimize the dimension of the parameter space. Additionally, the mathematical properties of the formalism enable a partition of the parameter space into independent subspaces. Each of these subspaces can be paired with a set of criteria that depend exclusively on it, thus allowing a separate sampling/screening in spaces of lower dimension. By applying tests in a strict order where computationally cheaper tests are applied first to each subspace and applying computationally expensive tests to the remaining subset thereafter, the use of resources is optimized and a larger number of models can be examined. This can be compared to a complex database query where the order of the requests can make a huge difference in the processing time. The method will be illustrated by analyzing a classical model of a metabolic pathway with end-product inhibition. Even for such a simple model, the method provides novel insight.

From Escherichia coli mutant 13C labeling data to a core kinetic model: A kinetic model parameterization pipeline

Kinetic models of metabolic networks offer the promise of quantitative phenotype prediction. The mechanistic characterization of enzyme catalyzed reactions allows for tracing the effect of perturbations in metabolite concentrations and reaction fluxes in response to genetic and environmental perturbation that are beyond the scope of stoichiometric models. In this study, we develop a two-step computational pipeline for the rapid parameterization of kinetic models of metabolic networks using a curated metabolic model and available 13C-labeling distributions under multiple genetic and environmental perturbations. The first step involves the elucidation of all intracellular fluxes in a core model of E. coli containing 74 reactions and 61 metabolites using 13C-Metabolic Flux Analysis (13C-MFA). Here, fluxes corresponding to the mid-exponential growth phase are elucidated for seven single gene deletion mutants from upper glycolysis, pentose phosphate pathway and the Entner-Doudoroff pathway. The computed flux ranges are then used to parameterize the same (i.e., k-ecoli74) core kinetic model for E. coli with 55 substrate-level regulations using the newly developed K-FIT parameterization algorithm. The K-FIT algorithm employs a combination of equation decomposition and iterative solution techniques to evaluate steady-state fluxes in response to genetic perturbations. k-ecoli74 predicted 86% of flux values for strains used during fitting within a single standard deviation of 13C-MFA estimated values. By performing both tasks using the same network, errors associated with lack of congruity between the two networks are avoided, allowing for seamless integration of data with model building. Product yield predictions and comparison with previously developed kinetic models indicate shifts in flux ranges and the presence or absence of mutant strains delivering flux towards pathways of interest from training data significantly impact predictive capabilities. Using this workflow, the impact of completeness of fluxomic datasets and the importance of specific genetic perturbations on uncertainties in kinetic parameter estimation are evaluated.

Uncertainty reduction in biochemical kinetic models: Enforcing desired model properties

A persistent obstacle for constructing kinetic models of metabolism is uncertainty in the kinetic properties of enzymes. Currently, available methods for building kinetic models can cope indirectly with uncertainties by integrating data from different biological levels and origins into models. In this study, we use the recently proposed computational approach iSCHRUNK (in Silico Approach to Characterization and Reduction of Uncertainty in the Kinetic Models), which combines Monte Carlo parameter sampling methods and machine learning techniques, in the context of Bayesian inference. Monte Carlo parameter sampling methods allow us to exploit synergies between different data sources and generate a population of kinetic models that are consistent with the available data and physicochemical laws. The machine learning allows us to data-mine the a priori generated kinetic parameters together with the integrated datasets and derive posterior distributions of kinetic parameters consistent with the observed physiology. In this work, we used iSCHRUNK to address a design question: can we identify which are the kinetic parameters and what are their values that give rise to a desired metabolic behavior? Such information is important for a wide variety of studies ranging from biotechnology to medicine. To illustrate the proposed methodology, we performed Metabolic Control Analysis, computed the flux control coefficients of the xylose uptake (XTR), and identified parameters that ensure a rate improvement of XTR in a glucose-xylose co-utilizing S. cerevisiae strain. Our results indicate that only three kinetic parameters need to be accurately characterized to describe the studied physiology, and ultimately to design and control the desired responses of the metabolism. This framework paves the way for a new generation of methods that will systematically integrate the wealth of available omics data and efficiently extract the information necessary for metabolic engineering and synthetic biology decisions.

Kinetic models are the most promising tool for understanding the complex dynamic behavior of living cells. The primary goal of kinetic models is to capture the properties of the metabolic networks as a whole, and thus we need large-scale models for dependable in silico analyses of metabolism. However, uncertainty in kinetic parameters impedes the development of kinetic models, and uncertainty levels increase with the model size. Tools that will address the issues with parameter uncertainty and that will be able to reduce the uncertainty propagation through the system are therefore needed. In this work, we applied a method called iSCHRUNK that combines parameter sampling and machine learning techniques to characterize the uncertainties and uncover intricate relationships between the parameters of kinetic models and the responses of the metabolic network. The proposed method allowed us to identify a small number of parameters that determine the responses in the network regardless of the values of other parameters. As a consequence, in future studies of metabolism, it will be sufficient to explore a reduced kinetic space, and more comprehensive analyses of large-scale and genome-scale metabolic networks will be computationally tractable.

Seeing the wood for the trees: a forest of methods for optimization and omic-network integration in metabolic modelling

Metabolic modelling has entered a mature phase with dozens of methods and software implementations available to the practitioner and the theoretician. It is not easy for a modeller to be able to see the wood (or the forest) for the trees. Driven by this analogy, we here present a 'forest' of principal methods used for constraint-based modelling in systems biology. This provides a tree-based view of methods available to prospective modellers, also available in interactive version at http://modellingmetabolism.net, where it will be kept updated with new methods after the publication of the present manuscript. Our updated classification of existing methods and tools highlights the most promising in the different branches, with the aim to develop a vision of how existing methods could hybridize and become more complex. We then provide the first hands-on tutorial for multi-objective optimization of metabolic models in R. We finally discuss the implementation of multi-view machine learning approaches in poly-omic integration. Throughout this work, we demonstrate the optimization of trade-offs between multiple metabolic objectives, with a focus on omic data integration through machine learning. We anticipate that the combination of a survey, a perspective on multi-view machine learning and a step-by-step R tutorial should be of interest for both the beginner and the advanced user.

Defining informative priors for ensemble modeling in systems biology

Ensemble modeling in molecular systems biology requires the reproducible translation of kinetic parameter data into informative probability distributions (priors), as well as approaches that sample parameters from these distributions without violating the thermodynamic consistency of the overall model. Although a number of pioneering frameworks for ensemble modeling have been published, the issue of generating informative priors has not yet been addressed. Here, we present a protocol that aims to fill this gap. This protocol discusses the collection of parameter values from a diverse range of sources (literature, databases and experiments), assessment of their plausibility, and creation of log-normal probability distributions that can be used as informative priors in ensemble modeling. Furthermore, the protocol enables sampling from the generated distributions while maintaining thermodynamic consistency. Once all parameter values have been retrieved from literature and databases, the protocol can be implemented within ~5-10 min per parameter. The aim of this protocol is to facilitate the design and use of informative distributions for ensemble modeling, especially in fields such as synthetic biology and systems medicine.

Metabolic kinetic modeling provides insight into complex biological questions, but hurdles remain

Metabolic models containing kinetic information can answer unique questions about cellular metabolism that are useful to metabolic engineering. Several kinetic modeling frameworks have recently been developed or improved. In addition, techniques for systematic identification of model structure, including regulatory interactions, have been reported. Each framework has advantages and limitations, which can make it difficult to choose the most appropriate framework. Common limitations are data availability and computational time, especially in large-scale modeling efforts. However, recently developed experimental techniques, parameter identification algorithms, as well as model reduction techniques help alleviate these computational bottlenecks. Opportunities for additional improvements may come from the rich literature in catalysis and chemical networks. In all, kinetic models are positioned to make significant impact in cellular engineering.

Allosteric Feedback Inhibition Enables Robust Amino Acid Biosynthesis in E. coli by Enforcing Enzyme Overabundance

Microbes must ensure robust amino acid metabolism in the face of external and internal perturbations. This robustness is thought to emerge from regulatory interactions in metabolic and genetic networks. Here, we explored the consequences of removing allosteric feedback inhibition in seven amino acid biosynthesis pathways in Escherichia coli (arginine, histidine, tryptophan, leucine, isoleucine, threonine, and proline). Proteome data revealed that enzyme levels decreased in five of the seven dysregulated pathways. Despite that, flux through the dysregulated pathways was not limited, indicating that enzyme levels are higher than absolutely needed in wild-type cells. We showed that such enzyme overabundance renders the arginine, histidine, and tryptophan pathways robust against perturbations of gene expression, using a metabolic model and CRISPR interference experiments. The results suggested a sensitive interaction between allosteric feedback inhibition and enzyme-level regulation that ensures robust yet efficient biosynthesis of histidine, arginine, and tryptophan in E. coli.

•

Amino acid biosynthesis enzymes do not normally operate at maximum capacity

•

Allosteric feedback inhibition ensures that enzymes are overabundant

•

Enzyme overabundance provides robustness against decreases in gene expression

An ensemble of mathematical models showing diauxic growth behaviour

Background

Carbon catabolite repression (CCR) controls the order in which different carbon sources are metabolised. Although this system is one of the paradigms of regulation in bacteria, the underlying mechanisms remain controversial. CCR involves the coordination of different subsystems of the cell - responsible for the uptake of carbon sources, their breakdown for the production of energy and precursors, and the conversion of the latter to biomass. The complexity of this integrated system, with regulatory mechanisms cutting across metabolism, gene expression, and signalling, has motivated important modelling efforts over the past four decades, especially in the enterobacterium Escherichia coli.

Results

Starting from a simple core model with only four intracellular metabolites, we develop an ensemble of model variants, all showing diauxic growth behaviour during a batch process. The model variants fall into one of the four categories: flux balance models, kinetic models with growth dilution, kinetic models with regulation, and resource allocation models. The model variants differ from one another in only a single aspect, each breaking the symmetry between the two substrate assimilation pathways in a different manner, and can be quantitatively compared using a so-called diauxic growth index. For each of the model variants, we predict the behaviour in two new experimental conditions, namely a glucose pulse for a culture growing in minimal medium with lactose and a batch culture with different initial concentrations of the components of the transport systems. When qualitatively comparing these predictions with experimental data for these two conditions, a number of models can be excluded while other model variants are still not discriminable. The best-performing model variants are based on inducer inclusion and activation of enzymatic genes by a global transcription factor, but the other proposed factors may complement these well-known regulatory mechanisms.

Conclusions

The model ensemble presented here offers a better understanding of the variety of mechanisms that have been proposed to play a role in CCR. In addition, it provides an educational resource for systems biology, as it gives an introduction to a broad range of modeling approaches in the context of a simple but biologically relevant example.

Electronic supplementary material

The online version of this article (10.1186/s12918-018-0604-8) contains supplementary material, which is available to authorized users.

Discovery and Evaluation of Biosynthetic Pathways for the Production of Five Methyl Ethyl Ketone Precursors

The limited supply of fossil fuels and the establishment of new environmental policies shifted research in industry and academia toward sustainable production of the second generation of biofuels, with methyl ethyl ketone (MEK) being one promising fuel candidate. MEK is a commercially valuable petrochemical with an extensive application as a solvent. However, as of today, a sustainable and economically viable production of MEK has not yet been achieved despite several attempts of introducing biosynthetic pathways in industrial microorganisms. We used BNICE.ch as a retrobiosynthesis tool to discover all novel pathways around MEK. Out of 1325 identified compounds connecting to MEK with one reaction step, we selected 3-oxopentanoate, but-3-en-2-one, but-1-en-2-olate, butylamine, and 2-hydroxy-2-methylbutanenitrile for further study. We reconstructed 3 679 610 novel biosynthetic pathways toward these 5 compounds. We then embedded these pathways into the genome-scale model of E. coli, and a set of 18 622 were found to be the most biologically feasible ones on the basis of thermodynamics and their yields. For each novel reaction in the viable pathways, we proposed the most similar KEGG reactions, with their gene and protein sequences, as candidates for either a direct experimental implementation or as a basis for enzyme engineering. Through pathway similarity analysis we classified the pathways and identified the enzymes and precursors that were indispensable for the production of the target molecules. These retrobiosynthesis studies demonstrate the potential of BNICE.ch for discovery, systematic evaluation, and analysis of novel pathways in synthetic biology and metabolic engineering studies.

Engineering an in vivo EP-bifido pathway in Escherichia coli for high-yield acetyl-CoA generation with low CO2 emission

The low carbon yield from native metabolic machinery produces unfavorable process economics during the biological conversion of substrates to desirable bioproducts. To obtain higher carbon yields, we constructed a carbon conservation pathway named EP-bifido pathway in Escherichia coli by combining Embden-Meyerhof-Parnas Pathway, Pentose Phosphate Pathway and "bifid shunt", to generate high yield acetyl-CoA from glucose. ¹³C-Metabolic flux analysis confirmed the successful and appropriate employment of the EP-bifido pathway. The CO₂ release during fermentation significantly reduced compared with the control strains. Then we demonstrated the in vivo effectiveness of the EP-bifido pathway using poly-β-hydroxybutyrate (PHB), mevalonate and fatty acids as example products. The engineered EP-bifido strains showed greatly improved PHB yield (from 26.0 mol% to 63.7 mol%), fatty acid yield (from 9.17% to 14.36%), and the highest mevalonate yield yet reported (64.3 mol% without considering the substrates used for cell mass formation). The synthetic pathway can be employed in the production of chemicals that use acetyl-CoA as a precursor and can be extended to other microorganisms.

Augmenting the Calvin-Benson-Bassham cycle by a synthetic malyl-CoA-glycerate carbon fixation pathway

The Calvin–Benson–Bassham (CBB) cycle is presumably evolved for optimal synthesis of C3 sugars, but not for the production of C2 metabolite acetyl-CoA. The carbon loss in producing acetyl-CoA from decarboxylation of C3 sugar limits the maximum carbon yield of photosynthesis. Here we design a synthetic malyl-CoA-glycerate (MCG) pathway to augment the CBB cycle for efficient acetyl-CoA synthesis. This pathway converts a C3 metabolite to two acetyl-CoA by fixation of one additional CO₂ equivalent, or assimilates glyoxylate, a photorespiration intermediate, to produce acetyl-CoA without net carbon loss. We first functionally demonstrate the design of the MCG pathway in vitro and in Escherichia coli. We then implement the pathway in a photosynthetic organism Synechococcus elongates PCC7942, and show that it increases the intracellular acetyl-CoA pool and enhances bicarbonate assimilation by roughly 2-fold. This work provides a strategy to improve carbon fixation efficiency in photosynthetic organisms.

Improving carbon fixation efficiency and reducing carbon loss have been long term goals for people working on photosynthetic organism improvement. Here, the authors design a synthetic malyl-CoA-glycerate pathway for efficient acetyl-CoA synthesis and verify its function in vitro, in E. coli and in cyanobacterium.

Construction and evolution of an Escherichia coli strain relying on nonoxidative glycolysis for sugar catabolism

We constructed an Escherichia coli strain that does not use glycolysis for sugar catabolism. Instead, it uses the synthetic nonoxidative glycolysis cycle to directly synthesize stoichiometric amounts of the two-carbon building block (acetyl-CoA), which is then converted to three-carbon metabolites to support growth. The resulting strain grows aerobically in glucose minimal medium and can achieve near-complete carbon conservation in the production of acetyl-CoA–derived products during anaerobic fermentation. This strain improves the theoretical carbon yield from 66.7% to 100% in acetyl-CoA–derived product formation.

The Embden–Meyerhoff–Parnas (EMP) pathway, commonly known as glycolysis, represents the fundamental biochemical infrastructure for sugar catabolism in almost all organisms, as it provides key components for biosynthesis, energy metabolism, and global regulation. EMP-based metabolism synthesizes three-carbon (C3) metabolites before two-carbon (C2) metabolites and must emit one CO₂ in the synthesis of the C2 building block, acetyl-CoA, a precursor for many industrially important products. Using rational design, genome editing, and evolution, here we replaced the native glycolytic pathways in Escherichia coli with the previously designed nonoxidative glycolysis (NOG), which bypasses initial C3 formation and directly generates stoichiometric amounts of C2 metabolites. The resulting strain, which contains 11 gene overexpressions, 10 gene deletions by design, and more than 50 genomic mutations (including 3 global regulators) through evolution, grows aerobically in glucose minimal medium but can ferment anaerobically to products with nearly complete carbon conservation. We confirmed that the strain metabolizes glucose through NOG by ¹³C tracer experiments. This redesigned E. coli strain represents a different approach for carbon catabolism and may serve as a useful platform for bioproduction.

Control of primary metabolism by a virulence regulatory network promotes robustness in a plant pathogen

Robustness is a key system-level property of living organisms to maintain their functions while tolerating perturbations. We investigate here how a regulatory network controlling multiple virulence factors impacts phenotypic robustness of a bacterial plant pathogen. We reconstruct a cell-scale model of Ralstonia solanacearum connecting a genome-scale metabolic network, a virulence macromolecule network, and a virulence regulatory network, which includes 63 regulatory components. We develop in silico methods to quantify phenotypic robustness under a broad set of conditions in high-throughput simulation analyses. This approach reveals that the virulence regulatory network exerts a control of the primary metabolism to promote robustness upon infection. The virulence regulatory network plugs into the primary metabolism mainly through the control of genes likely acquired via horizontal gene transfer, which results in a functional overlay with ancestral genes. These results support the view that robustness may be a selected trait that promotes pathogenic fitness upon infection.

Machine Learning Methods for Analysis of Metabolic Data and Metabolic Pathway Modeling

Machine learning uses experimental data to optimize clustering or classification of samples or features, or to develop, augment or verify models that can be used to predict behavior or properties of systems. It is expected that machine learning will help provide actionable knowledge from a variety of big data including metabolomics data, as well as results of metabolism models. A variety of machine learning methods has been applied in bioinformatics and metabolism analyses including self-organizing maps, support vector machines, the kernel machine, Bayesian networks or fuzzy logic. To a lesser extent, machine learning has also been utilized to take advantage of the increasing availability of genomics and metabolomics data for the optimization of metabolic network models and their analysis. In this context, machine learning has aided the development of metabolic networks, the calculation of parameters for stoichiometric and kinetic models, as well as the analysis of major features in the model for the optimal application of bioreactors. Examples of this very interesting, albeit highly complex, application of machine learning for metabolism modeling will be the primary focus of this review presenting several different types of applications for model optimization, parameter determination or system analysis using models, as well as the utilization of several different types of machine learning technologies.

Acceleration Strategies to Enhance Metabolic Ensemble Modeling Performance

Developing reliable, predictive kinetic models of metabolism is a difficult, yet necessary, priority toward understanding and deliberately altering cellular behavior. Constraint-based modeling has enabled the fields of metabolic engineering and systems biology to make great strides in interrogating cellular metabolism but does not provide sufficient insight into regulation or kinetic limitations of metabolic pathways. Moreover, the growth-optimized assumptions that constraint-based models often rely on do not hold when studying stationary or persistor cell populations. However, developing kinetic models provides many unique challenges, as many of the kinetic parameters and rate laws governing individual enzymes are unknown. Ensemble modeling (EM) was developed to circumnavigate this challenge and effectively sample the large kinetic parameter solution space using consistent experimental datasets. Unfortunately, EM, in its base form, requires long solve times to complete and often leads to unstable kinetic model predictions. Furthermore, these limitations scale prohibitively with increasing model size. As larger metabolic models are developed with increasing genetic information and experimental validation, the demand to incorporate kinetic information increases. Therefore, in this work, we have begun to tackle the challenges of EM by introducing additional steps to the existing method framework specifically through reducing computation time and optimizing parameter sampling. We first reduce the structural complexity of the network by removing dependent species, and second, we sample locally stable parameter sets to reflect realistic biological states of cells. Lastly, we presort the screening data to eliminate the most incorrect predictions in the earliest screening stages, saving further calculations in later stages. Our complementary improvements to this EM framework are easily incorporated into concurrent EM efforts and broaden the application opportunities and accessibility of kinetic modeling across the field.

A review of computational tools for design and reconstruction of metabolic pathways

Metabolic pathways reflect an organism's chemical repertoire and hence their elucidation and design have been a primary goal in metabolic engineering. Various computational methods have been developed to design novel metabolic pathways while taking into account several prerequisites such as pathway stoichiometry, thermodynamics, host compatibility, and enzyme availability. The choice of the method is often determined by the nature of the metabolites of interest and preferred host organism, along with computational complexity and availability of software tools. In this paper, we review different computational approaches used to design metabolic pathways based on the reaction network representation of the database (i.e., graph or stoichiometric matrix) and the search algorithm (i.e., graph search, flux balance analysis, or retrosynthetic search). We also put forth a systematic workflow that can be implemented in projects requiring pathway design and highlight current limitations and obstacles in computational pathway design.

Standardization in synthetic biology: an engineering discipline coming of age

BACKGROUND

Leaping DNA read-and-write technologies, and extensive automation and miniaturization are radically transforming the field of biological experimentation by providing the tools that enable the cost-effective high-throughput required to address the enormous complexity of biological systems. However, standardization of the synthetic biology workflow has not kept abreast with dwindling technical and resource constraints, leading, for example, to the collection of multi-level and multi-omics large data sets that end up disconnected or remain under- or even unexploited.

PURPOSE

In this contribution, we critically evaluate the various efforts, and the (limited) success thereof, in order to introduce standards for defining, designing, assembling, characterizing, and sharing synthetic biology parts. The causes for this success or the lack thereof, as well as possible solutions to overcome these, are discussed.

CONCLUSION

Akin to other engineering disciplines, extensive standardization will undoubtedly speed-up and reduce the cost of bioprocess development. In this respect, further implementation of synthetic biology standards will be crucial for the field in order to redeem its promise, i.e. to enable predictable forward engineering.

Formulation, construction and analysis of kinetic models of metabolism: A review of modelling frameworks

Kinetic models are critical to predict the dynamic behaviour of metabolic networks. Mechanistic kinetic models for large networks remain uncommon due to the difficulty of fitting their parameters. Recent modelling frameworks promise new ways to overcome this obstacle while retaining predictive capabilities. In this review, we present an overview of the relevant mathematical frameworks for kinetic formulation, construction and analysis. Starting with kinetic formalisms, we next review statistical methods for parameter inference, as well as recent computational frameworks applied to the construction and analysis of kinetic models. Finally, we discuss opportunities and limitations hindering the development of larger kinetic reconstructions.

Development of a core Clostridium thermocellum kinetic metabolic model consistent with multiple genetic perturbations

BACKGROUND

Clostridium thermocellum is a Gram-positive anaerobe with the ability to hydrolyze and metabolize cellulose into biofuels such as ethanol, making it an attractive candidate for consolidated bioprocessing (CBP). At present, metabolic engineering in C. thermocellum is hindered due to the incomplete description of its metabolic repertoire and regulation within a predictive metabolic model. Genome-scale metabolic (GSM) models augmented with kinetic models of metabolism have been shown to be effective at recapitulating perturbed metabolic phenotypes.

RESULTS

In this effort, we first update a second-generation genome-scale metabolic model (iCth446) for C. thermocellum by correcting cofactor dependencies, restoring elemental and charge balances, and updating GAM and NGAM values to improve phenotype predictions. The iCth446 model is next used as a scaffold to develop a core kinetic model (k-ctherm118) of the C. thermocellum central metabolism using the Ensemble Modeling (EM) paradigm. Model parameterization is carried out by simultaneously imposing fermentation yield data in lactate, malate, acetate, and hydrogen production pathways for 19 measured metabolites spanning a library of 19 distinct single and multiple gene knockout mutants along with 18 intracellular metabolite concentration data for a Δgldh mutant and ten experimentally measured Michaelis-Menten kinetic parameters.

CONCLUSIONS

The k-ctherm118 model captures significant metabolic changes caused by (1) nitrogen limitation leading to increased yields for lactate, pyruvate, and amino acids, and (2) ethanol stress causing an increase in intracellular sugar phosphate concentrations (~1.5-fold) due to upregulation of cofactor pools. Robustness analysis of k-ctherm118 alludes to the presence of a secondary activity of ketol-acid reductoisomerase and possible regulation by valine and/or leucine pool levels. In addition, cross-validation and robustness analysis allude to missing elements in k-ctherm118 and suggest additional experiments to improve kinetic model prediction fidelity. Overall, the study quantitatively assesses the advantages of EM-based kinetic modeling towards improved prediction of C. thermocellum metabolism and develops a predictive kinetic model which can be used to design biofuel-overproducing strains.