A kinetic model of Escherichia coli core metabolism satisfying multiple sets of mutant flux data

Use of genomics & proteomics in studying lipase producing microorganisms & its application

In biotechnological applications, lipases are recognized as the most widely utilized and versatile enzymes, pivotal in biocatalytic processes, predominantly produced by various microbial species. Utilizing omics technology, natural sources can be meticulously screened to find microbial flora which are responsible for oil production. Lipases are versatile biocatalysts. They are used in a variety of bioconversion reactions and are receiving a lot of attention because of the quick development of enzyme technology and its usefulness in industrial operations. This article offers recent insights into microbial lipase sources, including fungi, bacteria, and yeast, alongside traditional and modern methods of purification such as precipitation, immunopurification and chromatographic separation. Additionally, it explores innovative methods like the reversed micellar system, aqueous two-phase system (ATPS), and aqueous two-phase flotation (ATPF). The article deals with the use of microbial lipases in a variety of sectors, including the food, textile, leather, cosmetics, paper, detergent, while also critically analyzing lipase-producing microbes. Moreover, it highlights the role of lipases in biosensors, biodiesel production, tea processing, bioremediation, and racemization. This review provides the concept of the use of omics technique in the mechanism of screening of microbial species those are capable of producing lipase and also find the potential applications.

Purification and biochemical characterization of pullulanase produced from Bacillus sp. modified by ethyl-methyl sulfonate for improved applications

Strain improvement via chemical mutagen could impart traits with better enzyme production or improved characteristics. The present study sought to investigate the physicochemical properties of pullulanase produced from the wild Bacillus sp and the mutant. The pullulanases produced from the wild and the mutant Bacillus sp. (obtained via induction with ethyl methyl sulfonate) were purified in a-three step purification procedure and were also characterized. The wild and mutant pullulanases, which have molecular masses of 40 and 43.23 kDa, showed yields of 2.3% with 6.0-fold purification and 2.0% with 5.0-fold purification, respectively, and were most active at 50 and 40 °C and pH 7 and 8, respectively. The highest stability of the wild and mutant was between 40 and 50 °C after 1 h, although the mutant retained greater enzymatic activity between pH 6 and 9 than the wild. The mutant had a decreased Km of 0.03 mM as opposed to the wild type of 1.6 mM. In comparison to the wild, the mutant demonstrated a better capacity for tolerating metal ions and chelating agents. These exceptional characteristics of the mutant pullulanase may have been caused by a single mutation, which could improve its utility in industrial and commercial applications.

Iterative design of training data to control intricate enzymatic reaction networks

Kinetic modeling of in vitro enzymatic reaction networks is vital to understand and control the complex behaviors emerging from the nonlinear interactions inside. However, modeling is severely hampered by the lack of training data. Here, we introduce a methodology that combines an active learning-like approach and flow chemistry to efficiently create optimized datasets for a highly interconnected enzymatic reactions network with multiple sub-pathways. The optimal experimental design (OED) algorithm designs a sequence of out-of-equilibrium perturbations to maximize the information about the reaction kinetics, yielding a descriptive model that allows control of the output of the network towards any cost function. We experimentally validate the model by forcing the network to produce different product ratios while maintaining a minimum level of overall conversion efficiency. Our workflow scales with the complexity of the system and enables the optimization of previously unobtainable network outputs.

Local flux coordination and global gene expression regulation in metabolic modeling

Genome-scale metabolic networks (GSMs) are fundamental systems biology representations of a cell's entire set of stoichiometrically balanced reactions. However, such static GSMs do not incorporate the functional organization of metabolic genes and their dynamic regulation (e.g., operons and regulons). Specifically, there are numerous topologically coupled local reactions through which fluxes are coordinated; the global growth state often dynamically regulates many gene expression of metabolic reactions via global transcription factor regulators. Here, we develop a GSM reconstruction method, Decrem, by integrating locally coupled reactions and global transcriptional regulation of metabolism by cell state. Decrem produces predictions of flux and growth rates, which are highly correlated with those experimentally measured in both wild-type and mutants of three model microorganisms Escherichia coli, Saccharomyces cerevisiae, and Bacillus subtilis under various conditions. More importantly, Decrem can also explain the observed growth rates by capturing the experimentally measured flux changes between wild-types and mutants. Overall, by identifying and incorporating locally organized and regulated functional modules into GSMs, Decrem achieves accurate predictions of phenotypes and has broad applications in bioengineering, synthetic biology, and microbial pathology.

Monitoring and modelling the glutamine metabolic pathway: a review and future perspectives

BACKGROUND

Analysis of the glutamine metabolic pathway has taken a special place in metabolomics research in recent years, given its important role in cell biosynthesis and bioenergetics across several disorders, especially in cancer cell survival. The science of metabolomics addresses the intricate intracellular metabolic network by exploring and understanding how cells function and respond to external or internal perturbations to identify potential therapeutic targets. However, despite recent advances in metabolomics, monitoring the kinetics of a metabolic pathway in a living cell in situ, real-time and holistically remains a significant challenge.

AIM

This review paper explores the range of analytical approaches for monitoring metabolic pathways, as well as physicochemical modeling techniques, with a focus on glutamine metabolism. We discuss the advantages and disadvantages of each method and explore the potential of label-free Raman microspectroscopy, in conjunction with kinetic modeling, to enable real-time and in situ monitoring of the cellular kinetics of the glutamine metabolic pathway.

KEY SCIENTIFIC CONCEPTS

Given its important role in cell metabolism, the ability to monitor and model the glutamine metabolic pathways are highlighted. Novel, label free approaches have the potential to revolutionise metabolic biosensing, laying the foundation for a new paradigm in metabolomics research and addressing the challenges in monitoring metabolic pathways in living cells.

Optimal enzyme utilization suggests that concentrations and thermodynamics determine binding mechanisms and enzyme saturations

Deciphering the metabolic functions of organisms requires understanding the dynamic responses of living cells upon genetic and environmental perturbations, which in turn can be inferred from enzymatic activity. In this work, we investigate the optimal modes of operation for enzymes in terms of the evolutionary pressure driving them toward increased catalytic efficiency. We develop a framework using a mixed-integer formulation to assess the distribution of thermodynamic forces and enzyme states, providing detailed insights into the enzymatic mode of operation. We use this framework to explore Michaelis-Menten and random-ordered multi-substrate mechanisms. We show that optimal enzyme utilization is achieved by unique or alternative operating modes dependent on reactant concentrations. We find that in a bimolecular enzyme reaction, the random mechanism is optimal over any other ordered mechanism under physiological conditions. Our framework can investigate the optimal catalytic properties of complex enzyme mechanisms. It can further guide the directed evolution of enzymes and fill in the knowledge gaps in enzyme kinetics.

Comparative study of two Saccharomyces cerevisiae strains with kinetic models at genome-scale

The parameterization of kinetic models requires measurement of fluxes and/or metabolite levels for a base strain and a few genetic perturbations thereof. Unlike stoichiometric models that are mostly invariant to the specific strain, it remains unclear whether kinetic models constructed for different strains of the same species have similar or significantly different kinetic parameters. This important question underpins the applicability range and prediction limits of kinetic reconstructions. To this end, herein we parameterize two separate large-scale kinetic models using K-FIT with genome-wide coverage corresponding to two distinct strains of Saccharomyces cerevisiae: CEN.PK 113-7D strain (model k-sacce306-CENPK), and growth-deficient BY4741 (isogenic to S288c; model k-sacce306-BY4741). The metabolic network for each model contains 306 reactions, 230 metabolites, and 119 substrate-level regulatory interactions. The two models (for CEN.PK and BY4741) recapitulate, within one standard deviation, 77% and 75% of the fitted dataset fluxes, respectively, determined by ¹³C metabolic flux analysis for wild-type and eight single-gene knockout mutants of each strain. Strain-specific kinetic parameterization results indicate that key enzymes in the TCA cycle, glycolysis, and arginine and proline metabolism drive the metabolic differences between these two strains of S. cerevisiae. Our results suggest that although kinetic models cannot be readily used across strains as stoichiometric models, they can capture species-specific information through the kinetic parameterization process.

A dynamic kinetic model captures cell-free metabolism for improved butanol production

Cell-free systems are useful tools for prototyping metabolic pathways and optimizing the production of various bioproducts. Mechanistically-based kinetic models are uniquely suited to analyze dynamic experimental data collected from cell-free systems and provide vital qualitative insight. However, to date, dynamic kinetic models have not been applied with rigorous biological constraints or trained on adequate experimental data to the degree that they would give high confidence in predictions and broadly demonstrate the potential for widespread use of such kinetic models. In this work, we construct a large-scale dynamic model of cell-free metabolism with the goal of understanding and optimizing butanol production in a cell-free system. Using a combination of parameterization methods, the resultant model captures experimental metabolite measurements across two experimental conditions for nine metabolites at timepoints between 0 and 24 h. We present analysis of the model predictions, provide recommendations for butanol optimization, and identify the aldehyde/alcohol dehydrogenase as the primary bottleneck in butanol production. Sensitivity analysis further reveals the extent to which various parameters are constrained, and our approach for probing valid parameter ranges can be applied to other modeling efforts.

Synthetic Biology Meets Machine Learning

This chapter outlines the myriad applications of machine learning (ML) in synthetic biology, specifically in engineering cell and protein activity, and metabolic pathways. Though by no means comprehensive, the chapter highlights several prominent computational tools applied in the field and their potential use cases. The examples detailed reinforce how ML algorithms can enhance synthetic biology research by providing data-driven insights into the behavior of living systems, even without detailed knowledge of their underlying mechanisms. By doing so, ML promises to increase the efficiency of research projects by modeling hypotheses in silico that can then be tested through experiments. While challenges related to training dataset generation and computational costs remain, ongoing improvements in ML tools are paving the way for smarter and more streamlined synthetic biology workflows that can be readily employed to address grand challenges across manufacturing, medicine, engineering, agriculture, and beyond.

Computational design and engineering of an Escherichia coli strain producing the nonstandard amino acid para-aminophenylalanine

Introducing heterologous pathways into host cells constitutes a promising strategy for synthesizing nonstandard amino acids (nsAAs) to enable the production of proteins with expanded chemistries. However, this strategy has proven challenging, as the expression of heterologous pathways can disrupt cellular homeostasis of the host cell. Here, we sought to optimize the heterologous production of the nsAA para-aminophenylalanine (pAF) in Escherichia coli. First, we incorporated a heterologous pAF biosynthesis pathway into a genome-scale model of E. coli metabolism and computationally identified metabolic interventions in the host’s native metabolism to improve pAF production. Next, we explored different approaches of imposing these flux interventions experimentally and found that the upregulation of flux in the chorismate biosynthesis pathway through the elimination of feedback inhibition mechanisms could significantly raise pAF titers (∼20-fold) while maintaining a reasonable pAF production-growth rate trade-off. Overall, this study provides a promising strategy for the biosynthesis of nsAAs in engineered cells.

•

Sought to optimize para-aminophenylalanine (pAF) production and growth in E. coli

•

Identified interventions in the host native metabolism using genome-scale models

•

Constructed multiple mutant strains involving gene knockouts and/or overexpressions

•

Flux modification in chorismate biosynthesis pathway significantly raised pAF titer

Impact of Influenza A Virus Infection on Growth and Metabolism of Suspension MDCK Cells Using a Dynamic Model

Cell cultured-based influenza virus production is a viable option for vaccine manufacturing. In order to achieve a high concentration of viable cells, is requirement to have not only optimal process conditions, but also an active metabolism capable of intracellular synthesis of viral components. Experimental metabolic data collected in such processes are complex and difficult to interpret, for which mathematical models are an appropriate way to simulate and analyze the complex and dynamic interaction between the virus and its host cell. A dynamic model with 35 states was developed in this study to describe growth, metabolism, and influenza A virus production in shake flask cultivations of suspension Madin-Darby Canine Kidney (MDCK) cells. It considers cell growth (concentration of viable cells, mean cell diameters, volume of viable cells), concentrations of key metabolites both at the intracellular and extracellular level and virus titers. Using one set of parameters, the model accurately simulates the dynamics of mock-infected cells and correctly predicts the overall dynamics of virus-infected cells for up to 60 h post infection (hpi). The model clearly suggests that most changes observed after infection are related to cessation of cell growth and the subsequent transition to apoptosis and cell death. However, predictions do not cover late phases of infection, particularly for the extracellular concentrations of glutamate and ammonium after about 12 hpi. Results obtained from additional in silico studies performed indicated that amino acid degradation by extracellular enzymes resulting from cell lysis during late infection stages may contribute to this observed discrepancy.

Combining Kinetic and Constraint-Based Modelling to Better Understand Metabolism Dynamics

To understand the phenotypic capabilities of organisms, it is useful to characterise cellular metabolism through the analysis of its pathways. Dynamic mathematical modelling of metabolic networks is of high interest as it provides the time evolution of the metabolic components. However, it also has limitations, such as the necessary mechanistic details and kinetic parameters are not always available. On the other hand, large metabolic networks exhibit a complex topological structure which can be studied rather efficiently in their stationary regime by constraint-based methods. These methods produce useful predictions on pathway operations. In this review, we present both modelling techniques and we show how they bring complementary views of metabolism. In particular, we show on a simple example how both approaches can be used in conjunction to shed some light on the dynamics of metabolic networks.

Genome Scale Modeling to Study the Metabolic Competition between Cells in the Tumor Microenvironment

The tumor's physiology emerges from the dynamic interplay of numerous cell types, such as cancer cells, immune cells and stromal cells, within the tumor microenvironment. Immune and cancer cells compete for nutrients within the tumor microenvironment, leading to a metabolic battle between these cell populations. Tumor cells can reprogram their metabolism to meet the high demand of building blocks and ATP for proliferation, and to gain an advantage over the action of immune cells. The study of the metabolic reprogramming mechanisms underlying cancer requires the quantification of metabolic fluxes which can be estimated at the genome-scale with constraint-based or kinetic modeling. Constraint-based models use a set of linear constraints to simulate steady-state metabolic fluxes, whereas kinetic models can simulate both the transient behavior and steady-state values of cellular fluxes and concentrations. The integration of cell- or tissue-specific data enables the construction of context-specific models that reflect cell-type- or tissue-specific metabolic properties. While the available modeling frameworks enable limited modeling of the metabolic crosstalk between tumor and immune cells in the tumor stroma, future developments will likely involve new hybrid kinetic/stoichiometric formulations.

Recent advances in metabolic engineering-integration of in silico design and experimental analysis of metabolic pathways

Microorganisms are widely used to produce valuable compounds. Because thousands of metabolic reactions occur simultaneously and many metabolic reactions are related to target production and cell growth, the development of a rational design method for metabolic pathway modification to optimize target production is needed. In this paper, recent advances in metabolic engineering are reviewed, specifically considering computational pathway modification design and experimental evaluation of metabolic fluxes by ¹³C-metabolic flux analysis. Computational tools for seeking effective gene deletion targets and recruiting heterologous genes are described in flux balance analysis approaches. A kinetic model and adaptive laboratory evolution were applied to identify and eliminate the rate-limiting step in metabolic pathways. Data science-based approaches for process monitoring and control are described to maximize the performance of engineered cells in bioreactors.

A structural property for reduction of biochemical networks

Large-scale biochemical models are of increasing sizes due to the consideration of interacting organisms and tissues. Model reduction approaches that preserve the flux phenotypes can simplify the analysis and predictions of steady-state metabolic phenotypes. However, existing approaches either restrict functionality of reduced models or do not lead to significant decreases in the number of modelled metabolites. Here, we introduce an approach for model reduction based on the structural property of balancing of complexes that preserves the steady-state fluxes supported by the network and can be efficiently determined at genome scale. Using two large-scale mass-action kinetic models of Escherichia coli, we show that our approach results in a substantial reduction of 99% of metabolites. Applications to genome-scale metabolic models across kingdoms of life result in up to 55% and 85% reduction in the number of metabolites when arbitrary and mass-action kinetics is assumed, respectively. We also show that predictions of the specific growth rate from the reduced models match those based on the original models. Since steady-state flux phenotypes from the original model are preserved in the reduced, the approach paves the way for analysing other metabolic phenotypes in large-scale biochemical networks.

The effects of model complexity and size on metabolic flux distribution and control: case study in Escherichia coli

Background

Significant efforts have been made in building large-scale kinetic models of cellular metabolism in the past two decades. However, most kinetic models published to date, remain focused around central carbon pathways or are built around ad hoc reduced models without clear justification on their derivation and usage. Systematic algorithms exist for reducing genome-scale metabolic reconstructions to build thermodynamically feasible and consistently reduced stoichiometric models. However, it is important to study how network complexity affects conclusions derived from large-scale kinetic models built around consistently reduced models before we can apply them to study biological systems.

Results

We reduced the iJO1366 Escherichia Coli genome-scale metabolic reconstruction systematically to build three stoichiometric models of different size. Since the reduced models are expansions around the core subsystems for which the reduction was performed, the models are nested. We present a method for scaling up the flux profile and the concentration vector reference steady-states from the smallest model to the larger ones, whilst preserving maximum equivalency. Populations of kinetic models, preserving similarity in kinetic parameters, were built around the reference steady-states and their metabolic sensitivity coefficients (MSCs) were computed. The MSCs were sensitive to the model complexity. We proposed a metric for measuring the sensitivity of MSCs to these structural changes.

Conclusions

We proposed for the first time a workflow for scaling up the size of kinetic models while preserving equivalency between the kinetic models. Using this workflow, we demonstrate that model complexity in terms of networks size has significant impact on sensitivity characteristics of kinetic models. Therefore, it is essential to account for the effects of network complexity when constructing kinetic models. The presented metric for measuring MSC sensitivity to structural changes can guide modelers and experimentalists in improving model quality and guide synthetic biology and metabolic engineering. Our proposed workflow enables the testing of the suitability of a kinetic model for answering certain study-specific questions. We argue that the model-based metabolic design targets that are common across models of different size are of higher confidence, while those that are different could be the objective of investigations for model improvement.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04066-y.

Altered Metabolic Strategies: Elaborate Mechanisms Adopted by Oenococcus oeni in Response to Acid Stress

Oenococcus oeni plays a key role in inducing malolactic fermentation in wine. Acid stress is often encountered under wine conditions. However, the lack of systematic studies of acid resistance mechanisms limits the downstream fermentation applications. In this study, the acid responses of O. oeni were investigated by combining transcriptome, metabolome, and genome-scale metabolic modeling approaches. Metabolite profiling highlighted the decreased abundance of nucleotides under acid stress. The gene-metabolite bipartite network showed negative correlations between nucleotides and genes involved in ribosome assembly, translation, and post-translational processes, suggesting that stringent response could be activated under acid stress. Genome-scale metabolic modeling revealed marked flux rerouting, including reallocation of pyruvate, attenuation of glycolysis, utilization of carbon sources other than glucose, and enhancement of nucleotide salvage and the arginine deiminase pathway. This study provided novel insights into the acid responses of O. oeni, which will be useful for designing strategies to address acid stress in wine malolactic fermentation.

MASSpy: Building, simulating, and visualizing dynamic biological models in Python using mass action kinetics

Mathematical models of metabolic networks utilize simulation to study system-level mechanisms and functions. Various approaches have been used to model the steady state behavior of metabolic networks using genome-scale reconstructions, but formulating dynamic models from such reconstructions continues to be a key challenge. Here, we present the Mass Action Stoichiometric Simulation Python (MASSpy) package, an open-source computational framework for dynamic modeling of metabolism. MASSpy utilizes mass action kinetics and detailed chemical mechanisms to build dynamic models of complex biological processes. MASSpy adds dynamic modeling tools to the COnstraint-Based Reconstruction and Analysis Python (COBRApy) package to provide an unified framework for constraint-based and kinetic modeling of metabolic networks. MASSpy supports high-performance dynamic simulation through its implementation of libRoadRunner: the Systems Biology Markup Language (SBML) simulation engine. Three examples are provided to demonstrate how to use MASSpy: (1) a validation of the MASSpy modeling tool through dynamic simulation of detailed mechanisms of enzyme regulation; (2) a feature demonstration using a workflow for generating ensemble of kinetic models using Monte Carlo sampling to approximate missing numerical values of parameters and to quantify biological uncertainty, and (3) a case study in which MASSpy is utilized to overcome issues that arise when integrating experimental data with the computation of functional states of detailed biological mechanisms. MASSpy represents a powerful tool to address challenges that arise in dynamic modeling of metabolic networks, both at small and large scales.

Building kinetic models for metabolic engineering

Kinetic formalisms of metabolism link metabolic fluxes to enzyme levels, metabolite concentrations and their allosteric regulatory interactions. Though they require the identification of physiologically relevant values for numerous parameters, kinetic formalisms uniquely establish a mechanistic link across heterogeneous omics datasets and provide an overarching vantage point to effectively inform metabolic engineering strategies. Advances in computational power, gene annotation coverage, and formalism standardization have led to significant progress over the past few years. However, careful interpretation of model predictions, limited metabolic flux datasets, and assessment of parameter sensitivity remain as challenges. In this review we highlight fundamental considerations which influence model quality and prediction, advances in methodologies, and success stories of deploying kinetic models to guide metabolic engineering.

The Expanding Computational Toolbox for Engineering Microbial Phenotypes at the Genome Scale

Microbial strains are being engineered for an increasingly diverse array of applications, from chemical production to human health. While traditional engineering disciplines are driven by predictive design tools, these tools have been difficult to build for biological design due to the complexity of biological systems and many unknowns of their quantitative behavior. However, due to many recent advances, the gap between design in biology and other engineering fields is closing. In this work, we discuss promising areas of development of computational tools for engineering microbial strains. We define five frontiers of active research: (1) Constraint-based modeling and metabolic network reconstruction, (2) Kinetics and thermodynamic modeling, (3) Protein structure analysis, (4) Genome sequence analysis, and (5) Regulatory network analysis. Experimental and machine learning drivers have enabled these methods to improve by leaps and bounds in both scope and accuracy. Modern strain design projects will require these tools to be comprehensively applied to the entire cell and efficiently integrated within a single workflow. We expect that these frontiers, enabled by the ongoing revolution of big data science, will drive forward more advanced and powerful strain engineering strategies.

Recent advances in constraint and machine learning-based metabolic modeling by leveraging stoichiometric balances, thermodynamic feasibility and kinetic law formalisms

Understanding the governing principles behind organisms' metabolism and growth underpins their effective deployment as bioproduction chassis. A central objective of metabolic modeling is predicting how metabolism and growth are affected by both external environmental factors and internal genotypic perturbations. The fundamental concepts of reaction stoichiometry, thermodynamics, and mass action kinetics have emerged as the foundational principles of many modeling frameworks designed to describe how and why organisms allocate resources towards both growth and bioproduction. This review focuses on the latest algorithmic advancements that have integrated these foundational principles into increasingly sophisticated quantitative frameworks.

KiMoSys 2.0: an upgraded database for submitting, storing and accessing experimental data for kinetic modeling

The KiMoSys (https://kimosys.org), launched in 2014, is a public repository of published experimental data, which contains concentration data of metabolites, protein abundances and flux data. It offers a web-based interface and upload facility to share data, making it accessible in structured formats, while also integrating associated kinetic models related to the data. In addition, it also supplies tools to simplify the construction process of ODE (Ordinary Differential Equations)-based models of metabolic networks. In this release, we present an update of KiMoSys with new data and several new features, including (i) an improved web interface, (ii) a new multi-filter mechanism, (iii) introduction of data visualization tools, (iv) the addition of downloadable data in machine-readable formats, (v) an improved data submission tool, (vi) the integration of a kinetic model simulation environment and (vii) the introduction of a unique persistent identifier system. We believe that this new version will improve its role as a valuable resource for the systems biology community. Database URL:  www.kimosys.org.

Novel allosteric inhibition of phosphoribulokinase identified by ensemble kinetic modeling of Synechocystis sp. PCC 6803 metabolism

The present study attempted a computer simulation of the metabolism of a model cyanobacteria, Synechocystis sp. PCC 6803 (PCC 6803) to predict allosteric inhibitions that are likely to occur in photoautotrophic and mixotrophic conditions as well as in a metabolically engineered strain. PCC 6803 is a promising host for direct biochemical production from CO₂; however, further investigation of allosteric regulation is required for rational metabolic engineering to produce target compounds. Herein, ensemble modeling of microbial metabolism was applied to build accurate predictive models by synthesizing the results of multiple models with different parameter sets into a single score to identify plausible allosteric inhibitions. The data driven-computer simulation using metabolic flux, enzyme abundance, and metabolite concentration data successfully identified candidates for allosteric inhibition. The enzyme assay experiment using the recombinant protein confirmed isocitrate was a non-competitive inhibitor of phosphoribulokinase as a novel allosteric regulation of cyanobacteria metabolism.

Characterization of maximal enzyme catalytic rates in central metabolism of Arabidopsis thaliana

Availability of plant-specific enzyme kinetic data is scarce, limiting the predictive power of metabolic models and precluding identification of genetic factors of enzyme properties. Enzyme kinetic data are measured in vitro, often under non-physiological conditions, and conclusions elicited from modeling warrant caution. Here we estimate maximal in vivo catalytic rates for 168 plant enzymes, including photosystems I and II, cytochrome-b6f complex, ATP-citrate synthase, sucrose-phosphate synthase as well as enzymes from amino acid synthesis with previously undocumented enzyme kinetic data in BRENDA. The estimations are obtained by integrating condition-specific quantitative proteomics data, maximal rates of selected enzymes, growth measurements from Arabidopsis thaliana rosette with and fluxes through canonical pathways in a constraint-based model of leaf metabolism. In comparison to findings in Escherichia coli, we demonstrate weaker concordance between the plant-specific in vitro and in vivo enzyme catalytic rates due to a low degree of enzyme saturation. This is supported by the finding that concentrations of nicotinamide adenine dinucleotide (phosphate), adenosine triphosphate and uridine triphosphate, calculated based on our maximal in vivo catalytic rates, and available quantitative metabolomics data are below reported KM values and, therefore, indicate undersaturation of respective enzymes. Our findings show that genome-wide profiling of enzyme kinetic properties is feasible in plants, paving the way for understanding resource allocation.

Microbial lipases and their industrial applications: a comprehensive review

Lipases are very versatile enzymes, and produced the attention of the several industrial processes. Lipase can be achieved from several sources, animal, vegetable, and microbiological. The uses of microbial lipase market is estimated to be USD 425.0 Million in 2018 and it is projected to reach USD 590.2 Million by 2023, growing at a CAGR of 6.8% from 2018. Microbial lipases (EC 3.1.1.3) catalyze the hydrolysis of long chain triglycerides. The microbial origins of lipase enzymes are logically dynamic and proficient also have an extensive range of industrial uses with the manufacturing of altered molecules. The unique lipase (triacylglycerol acyl hydrolase) enzymes catalyzed the hydrolysis, esterification and alcoholysis reactions. Immobilization has made the use of microbial lipases accomplish its best performance and hence suitable for several reactions and need to enhance aroma to the immobilization processes. Immobilized enzymes depend on the immobilization technique and the carrier type. The choice of the carrier concerns usually the biocompatibility, chemical and thermal stability, and insolubility under reaction conditions, capability of easy rejuvenation and reusability, as well as cost proficiency. Bacillus spp., Achromobacter spp., Alcaligenes spp., Arthrobacter spp., Pseudomonos spp., of bacteria and Penicillium spp., Fusarium spp., Aspergillus spp., of fungi are screened large scale for lipase production. Lipases as multipurpose biological catalyst has given a favorable vision in meeting the needs for several industries such as biodiesel, foods and drinks, leather, textile, detergents, pharmaceuticals and medicals. This review represents a discussion on microbial sources of lipases, immobilization methods increased productivity at market profitability and reduce logistical liability on the environment and user.

Autocatalytic networks in biology: structural theory and algorithms

Self-sustaining autocatalytic networks play a central role in living systems, from metabolism at the origin of life, simple RNA networks and the modern cell, to ecology and cognition. A collectively autocatalytic network that can be sustained from an ambient food set is also referred to more formally as a 'reflexively autocatalytic food-generated' (RAF) set. In this paper, we first investigate a simplified setting for studying RAFs, which is nevertheless relevant to real biochemistry and which allows an exact mathematical analysis based on graph-theoretic concepts. This, in turn, allows for the development of efficient (polynomial-time) algorithms for questions that are computationally intractable (NP-hard) in the general RAF setting. We then show how this simplified setting for RAF systems leads naturally to a more general notion of RAFs that are 'generative' (they can be built up from simpler RAFs) and for which efficient algorithms carry over to this more general setting. Finally, we show how classical RAF theory can be extended to deal with ensembles of catalysts as well as the assignment of rates to reactions according to which catalysts (or combinations of catalysts) are available.

A dynamic model linking cell growth to intracellular metabolism and extracellular by-product accumulation

Mathematical modeling of animal cell growth and metabolism is essential for the understanding and improvement of the production of biopharmaceuticals. Models can explain the dynamic behavior of cell growth and product formation, support the identification of the most relevant parameters for process design, and significantly reduce the number of experiments to be performed for process optimization. Few dynamic models have been established that describe both extracellular and intracellular dynamics of growth and metabolism of animal cells. In this study, a model was developed, which comprises a set of 33 ordinary differential equations to describe batch cultivations of suspension AGE1.HN.AAT cells considered for the production of α1-antitrypsin. This model combines a segregated cell growth model with a structured model of intracellular metabolism. Overall, it considers the viable cell concentration, mean cell diameter, viable cell volume, concentration of extracellular substrates, and intracellular concentrations of key metabolites from the central carbon metabolism. Furthermore, the release of metabolic by-products such as lactate and ammonium was estimated directly from the intracellular reactions. Based on the same set of parameters, this model simulates well the dynamics of four independent batch cultivations. Analysis of the simulated intracellular rates revealed at least two distinct cellular physiological states. The first physiological state was characterized by a high glycolytic rate and high lactate production. Whereas the second state was characterized by efficient adenosine triphosphate production, a low glycolytic rate, and reactions of the TCA cycle running in the reverse direction from α-ketoglutarate to citrate. Finally, we show possible applications of the model for cell line engineering and media optimization with two case studies.

Bayesian inference of metabolic kinetics from genome-scale multiomics data

Modern biological tools generate a wealth of data on metabolite and protein concentrations that can be used to help inform new strain designs. However, learning from these data to predict how a cell will respond to genetic changes, a key need for engineering, remains challenging. A promising technique for leveraging omics measurements in metabolic modeling involves the construction of kinetic descriptions of the enzymatic reactions that occur within a cell. Parameterizing these models from biological data can be computationally difficult, since methods must also quantify the uncertainty in model parameters resulting from the observed data. While the field of Bayesian inference offers a wide range of methods for efficiently estimating distributions in parameter uncertainty, such techniques are poorly suited to traditional kinetic models due to their complex rate laws and resulting nonlinear dynamics. In this paper, we employ linear-logarithmic kinetics to simplify the calculation of steady-state flux distributions and enable efficient sampling and inference methods. We demonstrate that detailed information on the posterior distribution of parameters can be obtained efficiently at a variety of problem scales, including nearly genome-scale kinetic models trained on multiomics datasets. These results allow modern Bayesian machine learning tools to be leveraged in understanding biological data and in developing new, efficient strain designs.

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.

Particle-Based Simulation Reveals Macromolecular Crowding Effects on the Michaelis-Menten Mechanism

Many computational models for analyzing and predicting cell physiology rely on in vitro data collected in dilute and controlled buffer solutions. However, this can mislead models because up to 40% of the intracellular volume—depending on the organism, the physiology, and the cellular compartment—is occupied by a dense mixture of proteins, lipids, polysaccharides, RNA, and DNA. These intracellular macromolecules interfere with the interactions of enzymes and their reactants and thus affect the kinetics of biochemical reactions, making in vivo reactions considerably more complex than the in vitro data indicates. In this work, we present a new, to our knowledge, type of kinetics that captures and quantifies the effect of volume exclusion and other spatial phenomena on the kinetics of elementary reactions. We further developed a framework that allows for the efficient parameterization of these kinetics using particle simulations. Our formulation, entitled generalized elementary kinetics, can be used to analyze and predict the effect of intracellular crowding on enzymatic reactions and was herein applied to investigate the influence of crowding on phosphoglycerate mutase in Escherichia coli, which exhibits prototypical reversible Michaelis-Menten kinetics. Current research indicates that many enzymes are reaction limited and not diffusion limited, and our results suggest that the influence of fractal diffusion is minimal for these reaction-limited enzymes. Instead, increased association rates and decreased dissociation rates lead to a strong decrease in the effective maximal velocities V_max and the effective Michaelis-Menten constants K_M under physiologically relevant volume occupancies. Finally, the effects of crowding were explored in the context of a linear pathway, with the finding that crowding can have a redistributing effect on the effective flux responses in the case of twofold enzyme overexpression. We suggest that this framework, in combination with detailed kinetics models, will improve our understanding of enzyme reaction networks under nonideal conditions.

Human Systems Biology and Metabolic Modelling: A Review-From Disease Metabolism to Precision Medicine

In cell and molecular biology, metabolism is the only system that can be fully simulated at genome scale. Metabolic systems biology offers powerful abstraction tools to simulate all known metabolic reactions in a cell, therefore providing a snapshot that is close to its observable phenotype. In this review, we cover the 15 years of human metabolic modelling. We show that, although the past five years have not experienced large improvements in the size of the gene and metabolite sets in human metabolic models, their accuracy is rapidly increasing. We also describe how condition-, tissue-, and patient-specific metabolic models shed light on cell-specific changes occurring in the metabolic network, therefore predicting biomarkers of disease metabolism. We finally discuss current challenges and future promising directions for this research field, including machine/deep learning and precision medicine. In the omics era, profiling patients and biological processes from a multiomic point of view is becoming more common and less expensive. Starting from multiomic data collected from patients and N-of-1 trials where individual patients constitute different case studies, methods for model-building and data integration are being used to generate patient-specific models. Coupled with state-of-the-art machine learning methods, this will allow characterizing each patient's disease phenotype and delivering precision medicine solutions, therefore leading to preventative medicine, reduced treatment, and in silico clinical trials.

Approaches to Computational Strain Design in the Multiomics Era

Modern omics analyses are able to effectively characterize the genetic, regulatory, and metabolic phenotypes of engineered microbes, yet designing genetic interventions to achieve a desired phenotype remains challenging. With recent developments in genetic engineering techniques, timelines associated with building and testing strain designs have been greatly reduced, allowing for the first time an efficient closed loop iteration between experiment and analysis. However, the scale and complexity associated with multi-omics datasets complicates manual biological reasoning about the mechanisms driving phenotypic changes. Computational techniques therefore form a critical part of the Design-Build-Test-Learn (DBTL) cycle in metabolic engineering. Traditional statistical approaches can reduce the dimensionality of these datasets and identify common motifs among high-performing strains. While successful in many studies, these methods do not take full advantage of known connections between genes, proteins, and metabolic networks. There is therefore a growing interest in model-aided design, in which modeling frameworks from systems biology are used to integrate experimental data and generate effective and non-intuitive design predictions. In this mini-review, we discuss recent progress and challenges in this field. In particular, we compare methods augmenting flux balance analysis with additional constraints from fluxomic, genomic, and metabolomic datasets and methods employing kinetic representations of individual metabolic reactions, and machine learning. We conclude with a discussion of potential future directions for improving strain design predictions in the omics era and remaining experimental and computational hurdles.

Biochemical engineering in China

Chinese biochemical engineering is committed to supporting the chemical and food industries, to advance science and technology frontiers, and to meet major demands of Chinese society and national economic development. This paper reviews the development of biochemical engineering, strategic deployment of these technologies by the government, industrial demand, research progress, and breakthroughs in key technologies in China. Furthermore, the outlook for future developments in biochemical engineering in China is also discussed.

Cellular determinants of metabolite concentration ranges

Cellular functions are shaped by reaction networks whose dynamics are determined by the concentrations of underlying components. However, cellular mechanisms ensuring that a component’s concentration resides in a given range remain elusive. We present network properties which suffice to identify components whose concentration ranges can be efficiently computed in mass-action metabolic networks. We show that the derived ranges are in excellent agreement with simulations from a detailed kinetic metabolic model of Escherichia coli. We demonstrate that the approach can be used with genome-scale metabolic models to arrive at predictions concordant with measurements from Escherichia coli under different growth scenarios. By application to 14 genome-scale metabolic models from diverse species, our approach specifies the cellular determinants of concentration ranges that can be effectively employed to make predictions for a variety of biotechnological and medical applications.

We present a computational approach for inferring concentration ranges from genome-scale metabolic models. The approach specifies a determinant and molecular mechanism underling facile control of concentration ranges for components in large-scale cellular networks. Most importantly, the predictions about concentration ranges do not require knowledge of kinetic parameters (which are difficult to specify at a genome scale), provided measurements of concentrations in a reference state. The approach assumes that reaction rates follow the mass action law used in the derivations of other types of kinetics. We apply the approach with large-scale kinetic and stoichiometric metabolic models of organisms from different kingdoms of life to show that we can identify a proportion of metabolites to which our approach is applicable. By challenging the predictions of concentration ranges in the genome-scale metabolic network of E. coli with real-world data sets, we further demonstrate the prediction power and limitations of the approach.

Assessing Escherichia coli metabolism models and simulation approaches in phenotype predictions: Validation against experimental data

Over the last years, several genome-scale metabolic models (GEMs) and kinetic models of Escherichia coli were published. Their predictive performance varies according to the evaluation metric considered, the computational simulation methods used, and the type/quality of experimental data available. However, the GEM approach is often not compared with the kinetic modeling framework. Also, the different genome-scale reconstruction versions and simulation methods of mutant phenotypes are usually not validated to predict intracellular fluxes using large experimental datasets. Here, we intended to (i) systematically evaluate the prediction performance of three E. coli GEMs (iJR904, iAF1260, and iJO1366) available in the literature according to predictive growth metrics (intracellular flux distribution); (ii) assess the reliability of a E. coli GEM in the prediction of gene knockout phenotypes when different simulation methods (parsimonious flux balance analysis, Minimization of Metabolic Adjustment, linear version of MoMA, Regulatory on/off minimization, and Minimization of Metabolites Balance) are used; and finally (iii) investigate the flux distribution predictive power of the constrained-based modeling approach (selected stoichiometric GEM) and compare it with the kinetic modeling approach (two published kinetic models) for E. coli central metabolism, in order to assess their accuracy. Results show that the phenotype predictions were not significantly sensitive to the metabolic models, although the GEM iAF1260 was more accurate in the prediction of central carbon fluxes at low dilution rates. Furthermore, we observed that the choice of the appropriate simulation method of mutant phenotypes depends on the biological question to be addressed. In terms of the two modeling approaches, none outperformed the other for all the tested scenarios. © 2018 American Institute of Chemical Engineers Biotechnol. Prog., 34:1344-1354, 2018.

Finite-time parametric identification for the model representing the metabolic and genetic regulatory effects of sequential aerobic respiration and anaerobic fermentation processes in Escherichia coli

Mathematical modelling applied to biological systems allows for the inferring of changes in the dynamic behaviour of organisms associated with variations in the environment. Models based on ordinary differential equations are most commonly used because of their ability to describe the mechanisms of biological systems such as transcription. The disadvantage of using this approach is that there is a large number of parameters involved and that it is difficult to obtain them experimentally. This study presents an algorithm to obtain a finite-time parameter characterization of the model used to describe changes in the metabolic behaviour of Escherichia coli associated with environmental changes. In this scheme, super-twisting algorithm was proposed to recover the derivative of all the proteins and mRNA of E. coli associated to changes in the concentration of oxygen available in the growth media. The 75 identified parameters in this study maintain the biological coherence of the system and they were estimated with no more than 20% error with respect to the real ones included in the proposed model.

Metabolic engineering to increase crop yield: From concept to execution

Although the return on investment over the last 20 years for mass screening of individual plant genes to improve crop performance has been low, the investment in these activities was essential to establish the infrastructure and tools of modern plant genomics. Complex traits such as crop yield are likely multigenic, and the exhaustive screening of random gene combinations to achieve yield gains is not realistic. Clearly, smart approaches must be developed. In silico analyses of plant metabolism and gene networks can move a trait discovery program beyond trial-and-error approaches and towards rational design strategies. Metabolic models employing flux-balance analysis are useful to determine the contribution of individual genes to a trait, or to compare, optimize, or even design metabolic pathways. Regulatory association networks provide a transcriptome-based view of the plant and can lead to the identification of transcription factors that control expression of multiple genes affecting a trait. In this review, the use of these models from the perspective of an Ag innovation company's trait discovery and development program will be discussed. Important decisions that can have significant impacts on the cost and timeline to develop a commercial trait will also be presented.

A Review of Dynamic Modeling Approaches and Their Application in Computational Strain Optimization for Metabolic Engineering

Mathematical modeling is a key process to describe the behavior of biological networks. One of the most difficult challenges is to build models that allow quantitative predictions of the cells' states along time. Recently, this issue started to be tackled through novel in silico approaches, such as the reconstruction of dynamic models, the use of phenotype prediction methods, and pathway design via efficient strain optimization algorithms. The use of dynamic models, which include detailed kinetic information of the biological systems, potentially increases the scope of the applications and the accuracy of the phenotype predictions. New efforts in metabolic engineering aim at bridging the gap between this approach and other different paradigms of mathematical modeling, as constraint-based approaches. These strategies take advantage of the best features of each method, and deal with the most remarkable limitation—the lack of available experimental information—which affects the accuracy and feasibility of solutions. Parameter estimation helps to solve this problem, but adding more computational cost to the overall process. Moreover, the existing approaches include limitations such as their scalability, flexibility, convergence time of the simulations, among others. The aim is to establish a trade-off between the size of the model and the level of accuracy of the solutions. In this work, we review the state of the art of dynamic modeling and related methods used for metabolic engineering applications, including approaches based on hybrid modeling. We describe approaches developed to undertake issues regarding the mathematical formulation and the underlying optimization algorithms, and that address the phenotype prediction by including available kinetic rate laws of metabolic processes. Then, we discuss how these have been used and combined as the basis to build computational strain optimization methods for metabolic engineering purposes, how they lead to bi-level schemes that can be used in the industry, including a consideration of their limitations.

DMPy: a Python package for automated mathematical model construction of large-scale metabolic systems

Background

Models of metabolism are often used in biotechnology and pharmaceutical research to identify drug targets or increase the direct production of valuable compounds. Due to the complexity of large metabolic systems, a number of conclusions have been drawn using mathematical methods with simplifying assumptions. For example, constraint-based models describe changes of internal concentrations that occur much quicker than alterations in cell physiology. Thus, metabolite concentrations and reaction fluxes are fixed to constant values. This greatly reduces the mathematical complexity, while providing a reasonably good description of the system in steady state. However, without a large number of constraints, many different flux sets can describe the optimal model and we obtain no information on how metabolite levels dynamically change. Thus, to accurately determine what is taking place within the cell, finer quality data and more detailed models need to be constructed.

Results

In this paper we present a computational framework, DMPy, that uses a network scheme as input to automatically search for kinetic rates and produce a mathematical model that describes temporal changes of metabolite fluxes. The parameter search utilises several online databases to find measured reaction parameters. From this, we take advantage of previous modelling efforts, such as Parameter Balancing, to produce an initial mathematical model of a metabolic pathway. We analyse the effect of parameter uncertainty on model dynamics and test how recent flux-based model reduction techniques alter system properties. To our knowledge this is the first time such analysis has been performed on large models of metabolism. Our results highlight that good estimates of at least 80% of the reaction rates are required to accurately model metabolic systems. Furthermore, reducing the size of the model by grouping reactions together based on fluxes alters the resulting system dynamics.

Conclusion

The presented pipeline automates the modelling process for large metabolic networks. From this, users can simulate their pathway of interest and obtain a better understanding of how altering conditions influences cellular dynamics. By testing the effects of different parameterisations we are also able to provide suggestions to help construct more accurate models of complete metabolic systems in the future.

Electronic supplementary material

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

A machine learning approach to predict metabolic pathway dynamics from time-series multiomics data

New synthetic biology capabilities hold the promise of dramatically improving our ability to engineer biological systems. However, a fundamental hurdle in realizing this potential is our inability to accurately predict biological behavior after modifying the corresponding genotype. Kinetic models have traditionally been used to predict pathway dynamics in bioengineered systems, but they take significant time to develop, and rely heavily on domain expertise. Here, we show that the combination of machine learning and abundant multiomics data (proteomics and metabolomics) can be used to effectively predict pathway dynamics in an automated fashion. The new method outperforms a classical kinetic model, and produces qualitative and quantitative predictions that can be used to productively guide bioengineering efforts. This method systematically leverages arbitrary amounts of new data to improve predictions, and does not assume any particular interactions, but rather implicitly chooses the most predictive ones.