Dynamic Modeling of Cell-Free Biochemical Networks Using Effective Kinetic Models

Cell-free systems offer many advantages for the study, manipulation and modeling of metabolism compared to in vivo processes. Many of the challenges confronting genome-scale kinetic modeling can potentially be overcome in a cell-free system. For example, there is no complex transcriptional regulation to consider, transient metabolic measurements are easier to obtain, and we no longer have to consider cell growth. Thus, cell-free operation holds several significant advantages for model development, identification and validation. Theoretically, genome-scale cell-free kinetic models may be possible for industrially important organisms, such as E. coli, if a simple, tractable framework for integrating allosteric regulation with enzyme kinetics can be formulated. Toward this unmet need, we present an effective biochemical network modeling framework for building dynamic cell-free metabolic models. The key innovation of our approach is the integration of simple effective rules encoding complex allosteric regulation with traditional kinetic pathway modeling. We tested our approach by modeling the time evolution of several hypothetical cell-free metabolic networks. We found that simple effective rules, when integrated with traditional enzyme kinetic expressions, captured complex allosteric patterns such as ultrasensitivity or non-competitive inhibition in the absence of mechanistic information. Second, when integrated into network models, these rules captured classic regulatory patterns such as product-induced feedback inhibition. Lastly, we showed, at least for the network architectures considered here, that we could simultaneously estimate kinetic parameters and allosteric connectivity from synthetic data starting from an unbiased collection of possible allosteric structures using particle swarm optimization. However, when starting with an initial population that was heavily enriched with incorrect structures, our particle swarm approach could converge to an incorrect structure. While only an initial proof-of-concept, the framework presented here could be an important first step toward genome-scale cell-free kinetic modeling of the biosynthetic capacity of industrially important organisms.

[1]  Diran Basmadjian,et al.  A mathematical model of thrombin production in blood coagulation, Part I: The sparsely covered membrane case , 1994, Annals of Biomedical Engineering.

[2]  Michael R Pinsky,et al.  Current Evidence Based Guidelines for Factor VIIa Use in Trauma: The Good, the Bad, and the Ugly , 2008, The American surgeon.

[3]  J. Delâge [Coagulation]. , 2018, Laval medical.

[4]  M A Savageau,et al.  Biochemical systems theory: operational differences among variant representations and their significance. , 1991, Journal of theoretical biology.

[5]  A. Fogelson,et al.  Surface-mediated control of blood coagulation: the role of binding site densities and platelet deposition. , 2001, Biophysical journal.

[6]  Rustem F Ismagilov,et al.  Effects of shear rate on propagation of blood clotting determined using microfluidics and numerical simulations. , 2008, Journal of the American Chemical Society.

[7]  M M Domach,et al.  Computer model for glucose‐limited growth of a single cell of Escherichia coli B/r‐A , 1984, Biotechnology and bioengineering.

[8]  B. Kholodenko,et al.  Computational Approaches for Analyzing Information Flow in Biological Networks , 2012, Science Signaling.

[9]  Manolis Gavaises,et al.  A simplified mathematical model for thrombin generation. , 2014, Medical engineering & physics.

[10]  M. Shuler,et al.  A Single‐Cell Model for CHO Cells a , 1992, Annals of the New York Academy of Sciences.

[11]  K. Mann,et al.  Biochemistry and Physiology of Blood Coagulation , 1999, Thrombosis and Haemostasis.

[12]  James H Morrissey,et al.  Polyphosphate modulates blood coagulation and fibrinolysis. , 2006, Proceedings of the National Academy of Sciences of the United States of America.

[13]  M A Khanin,et al.  A mathematical model of the kinetics of blood coagulation. , 1989, Journal of theoretical biology.

[14]  Hisashi Kobayashi,et al.  Modeling and analysis , 1978 .

[15]  H C Hemker,et al.  Simulation model for thrombin generation in plasma. , 1991, Haemostasis.

[16]  Riccardo Poli,et al.  Particle swarm optimization , 1995, Swarm Intelligence.

[17]  Julia Mitchell,et al.  Haemophilia and inhibitors. 1: Diagnosis and treatment. , 2008, Nursing times.

[18]  Hisashi Adachi,et al.  Inflammation, atherosclerosis, and coronary artery disease. , 2005, The New England journal of medicine.

[19]  M. Schenone,et al.  The blood coagulation cascade , 2004, Current opinion in hematology.

[20]  F. Doyle,et al.  A benchmark for methods in reverse engineering and model discrimination: problem formulation and solutions. , 2004, Genome research.

[21]  C. Danforth,et al.  Defining the Boundaries of Normal Thrombin Generation: Investigations into Hemostasis , 2012, PloS one.

[22]  Rustem F Ismagilov,et al.  Minimal functional model of hemostasis in a biomimetic microfluidic system. , 2004, Angewandte Chemie.

[23]  Ken Lo,et al.  Stochastic Modeling of Blood Coagulation Initiation , 2006, Pathophysiology of Haemostasis and Thrombosis.

[24]  D. Ramkrishna,et al.  Systematic development of hybrid cybernetic models: Application to recombinant yeast co‐consuming glucose and xylose , 2009, Biotechnology and bioengineering.

[25]  K. Mann,et al.  Thrombin functions during tissue factor-induced blood coagulation. , 2002, Blood.

[26]  J. Heijnen,et al.  Dynamic simulation and metabolic re-design of a branched pathway using linlog kinetics. , 2003, Metabolic engineering.

[27]  Sang Ok Song,et al.  Ensembles of signal transduction models using Pareto Optimal Ensemble Techniques (POETs). , 2010, Biotechnology journal.

[28]  Marshall W. Nirenberg,et al.  The dependence of cell-free protein synthesis in E. coli upon naturally occurring or synthetic polyribonucleotides , 1961, Proceedings of the National Academy of Sciences.

[29]  Deyan Luan,et al.  Computationally Derived Points of Fragility of a Human Cascade Are Consistent with Current Therapeutic Strategies , 2007, PLoS Comput. Biol..

[30]  Jalal K. Siddiqui,et al.  Modelling and analysis of an ensemble of eukaryotic translation initiation models. , 2011, IET systems biology.

[31]  J. Sethna,et al.  Parameter Space Compression Underlies Emergent Theories and Predictive Models , 2013, Science.

[32]  G. Broze,et al.  Factor XI activation in a revised model of blood coagulation , 1991, Science.

[33]  B. Palsson,et al.  The Escherichia coli MG1655 in silico metabolic genotype: its definition, characteristics, and capabilities. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[34]  A G Fredrickson,et al.  Formulation of structured growth models. , 2000, Biotechnology and bioengineering.

[35]  M. L. Shuler,et al.  Structured model for Saccharomyces cerevisiae , 1989 .

[36]  T. Orfeo,et al.  The Prothrombotic Phenotypes in Familial Protein C Deficiency Are Differentiated by Computational Modeling of Thrombin Generation , 2012, PloS one.

[37]  A. Shapiro,et al.  Single-dose recombinant activated factor VII for the treatment of joint bleeds in hemophilia patients with inhibitors. , 2008, Clinical advances in hematology & oncology : H&O.

[38]  Peter K. Sorger,et al.  Logic-Based Models for the Analysis of Cell Signaling Networks† , 2010, Biochemistry.

[39]  Wolfgang Wiechert,et al.  Translating biochemical network models between different kinetic formats. , 2009, Metabolic engineering.

[40]  Thomas Orfeo,et al.  The Significance of Circulating Factor IXa in Blood* , 2004, Journal of Biological Chemistry.

[41]  J. Reifman,et al.  Kinetic modeling sheds light on the mode of action of recombinant factor VIIa on thrombin generation. , 2011, Thrombosis research.

[42]  Nathan E Lewis,et al.  Analysis of omics data with genome-scale models of metabolism. , 2013, Molecular bioSystems.

[43]  S. Diamond,et al.  Hierarchical organization in the hemostatic response and its relationship to the platelet-signaling network. , 2013, Blood.

[44]  D. Ramkrishna,et al.  Metabolic engineering from a cybernetic perspective: aspartate family of amino acids. , 1999, Metabolic engineering.

[45]  R. Ismagilov,et al.  Modular chemical mechanism predicts spatiotemporal dynamics of initiation in the complex network of hemostasis , 2006, Proceedings of the National Academy of Sciences.

[46]  Doraiswami Ramkrishna,et al.  Prediction of dynamic behavior of mutant strains from limited wild-type data. , 2012, Metabolic engineering.

[47]  Edward G. D. Tuddenham,et al.  The Molecular Genetics of Haemostasis and Its Inherited Disorders , 1994 .

[48]  Kenneth G. Mann,et al.  Inhibitory Mechanism of the Protein C Pathway on Tissue Factor-induced Thrombin Generation , 1997, The Journal of Biological Chemistry.

[49]  Adam M. Feist,et al.  Basic and applied uses of genome-scale metabolic network reconstructions of Escherichia coli , 2013, Molecular systems biology.

[50]  Jordan C. Atlas,et al.  Incorporating genome-wide DNA sequence information into a dynamic whole-cell model of Escherichia coli: application to DNA replication. , 2008, IET systems biology.

[51]  R. Bird,et al.  Review article: Coagulation cascade and therapeutics update: Relevance to nephrology. Part 1: Overview of coagulation, thrombophilias and history of anticoagulants , 2009, Nephrology.

[52]  B. Palsson,et al.  Escherichia coli K-12 undergoes adaptive evolution to achieve in silico predicted optimal growth , 2002, Nature.

[53]  H. Tien,et al.  The use of recombinant factor VIIa in trauma patients. , 2009, The Journal of the American Academy of Orthopaedic Surgeons.

[54]  Aaron Fogelson,et al.  An overview of mathematical modeling of thrombus formation under flow. , 2014, Thrombosis research.

[55]  Scott L. Diamond,et al.  Systems Biology of Coagulation Initiation: Kinetics of Thrombin Generation in Resting and Activated Human Blood , 2010, PLoS Comput. Biol..

[56]  Adam M. Feist,et al.  A genome-scale metabolic reconstruction for Escherichia coli K-12 MG1655 that accounts for 1260 ORFs and thermodynamic information , 2007, Molecular systems biology.

[57]  J. Bailey Complex biology with no parameters , 2001, Nature Biotechnology.

[58]  Kapil G. Gadkar,et al.  Cybernetic Model Predictive Control of a Continuous Bioreactor with Cell Recycle , 2003, Biotechnology progress.

[59]  G. T. Tsao,et al.  Investigation of bacterial growth on mixed substrates: Experimental evaluation of cybernetic models , 1986, Biotechnology and bioengineering.

[60]  B. Palsson,et al.  Regulation of gene expression in flux balance models of metabolism. , 2001, Journal of theoretical biology.

[61]  J. Peyman,et al.  Regulation of gene Expression , 2012 .

[62]  B. Palsson,et al.  Constraining the metabolic genotype–phenotype relationship using a phylogeny of in silico methods , 2012, Nature Reviews Microbiology.

[63]  Kenichi A Tanaka,et al.  Blood Coagulation: Hemostasis and Thrombin Regulation , 2009, Anesthesia and analgesia.

[64]  Igor Goryanin,et al.  Kinetic Model of phosphofructokinase-1 from Escherichia coli , 2008, J. Bioinform. Comput. Biol..

[65]  Jeffrey D Varner,et al.  Modeling and analysis of retinoic acid induced differentiation of uncommitted precursor cells. , 2011, Integrative biology : quantitative biosciences from nano to macro.

[66]  Brittany E. Bannish,et al.  Modelling fibrinolysis: a 3D stochastic multiscale model. , 2014, Mathematical medicine and biology : a journal of the IMA.

[67]  U. Hedner,et al.  Factor VIIa and its potential therapeutic use in bleeding-associated pathologies , 2008, Thrombosis and Haemostasis.

[68]  U. Sauer,et al.  Systematic evaluation of objective functions for predicting intracellular fluxes in Escherichia coli , 2007, Molecular systems biology.

[69]  Doraiswami Ramkrishna,et al.  Exacting predictions by cybernetic model confirmed experimentally: Steady state multiplicity in the chemostat , 2012, Biotechnology progress.

[70]  U. Sauer,et al.  Systematic identification of allosteric protein-metabolite interactions that control enzyme activity in vivo , 2013, Nature Biotechnology.

[71]  Markus J. Herrgård,et al.  Integrating high-throughput and computational data elucidates bacterial networks , 2004, Nature.

[72]  Paola Annoni,et al.  Variance based sensitivity analysis of model output. Design and estimator for the total sensitivity index , 2010, Comput. Phys. Commun..

[73]  Douglas A Lauffenburger,et al.  Querying quantitative logic models (Q2LM) to study intracellular signaling networks and cell-cytokine interactions , 2011, Biotechnology journal.

[74]  Mingsheng Zhang,et al.  Comparing signaling networks between normal and transformed hepatocytes using discrete logical models. , 2011, Cancer research.

[75]  Rudiyanto Gunawan,et al.  Iterative approach to model identification of biological networks , 2005, BMC Bioinformatics.

[76]  M. Jewett,et al.  Cell-free synthetic biology: thinking outside the cell. , 2012, Metabolic engineering.

[77]  Peter D. Karp,et al.  EcoCyc: fusing model organism databases with systems biology , 2012, Nucleic Acids Res..

[78]  J. Varner,et al.  Large-scale prediction of phenotype: concept. , 2000, Biotechnology and bioengineering.

[79]  Morten H. Christiansen,et al.  A computational model , 2014 .

[80]  K. Mann,et al.  Surface-dependent reactions of the vitamin K-dependent enzyme complexes. , 1990, Blood.

[81]  A. Imbriano [Blood coagulation]. , 1955, La Semana medica.

[82]  E G Tuddenham,et al.  The hemophilias--from royal genes to gene therapy. , 2001, The New England journal of medicine.

[83]  Michael C Jewett,et al.  An integrated cell-free metabolic platform for protein production and synthetic biology , 2008, Molecular systems biology.

[84]  J. Reifman,et al.  Kinetic model facilitates analysis of fibrin generation and its modulation by clotting factors: implications for hemostasis-enhancing therapies. , 2014, Molecular bioSystems.

[85]  F J Doyle,et al.  Model identification of signal transduction networks from data using a state regulator problem. , 2005, Systems biology.

[86]  J. Stelling,et al.  Ensemble modeling for analysis of cell signaling dynamics , 2007, Nature Biotechnology.

[87]  J. Heijnen Approximative kinetic formats used in metabolic network modeling , 2005, Biotechnology and bioengineering.

[88]  Rustem F Ismagilov,et al.  Propagation of blood clotting in the complex biochemical network of hemostasis is described by a simple mechanism. , 2007, Journal of the American Chemical Society.

[89]  Dougald M Monroe,et al.  Manipulation of prothrombin concentration improves response to high‐dose factor VIIa in a cell‐based model of haemophilia , 2006, British journal of haematology.

[90]  K. Mann,et al.  A model for the tissue factor pathway to thrombin. I. An empirical study. , 1994, The Journal of biological chemistry.

[91]  Adam M. Feist,et al.  Reconstruction of biochemical networks in microorganisms , 2009, Nature Reviews Microbiology.

[92]  James R. Swartz,et al.  Production and stabilization of the trimeric influenza hemagglutinin stem domain for potentially broadly protective influenza vaccines , 2013, Proceedings of the National Academy of Sciences.

[93]  K. C. Jones,et al.  A Model for the Stoichiometric Regulation of Blood Coagulation* , 2002, The Journal of Biological Chemistry.

[94]  D. Ramkrishna,et al.  Cybernetic models based on lumped elementary modes accurately predict strain‐specific metabolic function , 2011, Biotechnology and bioengineering.

[95]  Eric Jones,et al.  SciPy: Open Source Scientific Tools for Python , 2001 .

[96]  Y. Nakatomi,et al.  A novel therapeutic approach combining human plasma‐derived Factors VIIa and X for haemophiliacs with inhibitors: evidence of a higher thrombin generation rate in vitro and more sustained haemostatic activity in vivo than obtained with Factor VIIa alone , 2003, Vox sanguinis.

[97]  K. Mann,et al.  A model for the tissue factor pathway to thrombin. II. A mathematical simulation. , 1994, Journal of Biological Chemistry.

[98]  Deyan Luan,et al.  Ensembles of uncertain mathematical models can identify network response to therapeutic interventions. , 2010, Molecular bioSystems.

[99]  I. Sobola,et al.  Global sensitivity indices for nonlinear mathematical models and their Monte Carlo estimates , 2001 .

[100]  J. Liao,et al.  Ensemble modeling of metabolic networks. , 2008, Biophysical journal.

[101]  M L Shuler,et al.  A modular minimal cell model: purine and pyrimidine transport and metabolism. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[102]  Jeffrey Varner,et al.  Modeling and Analysis of the Molecular Basis of Pain in Sensory Neurons , 2009, PloS one.

[103]  Jeffrey D. Varner,et al.  Analysis of the Molecular Networks in Androgen Dependent and Independent Prostate Cancer Revealed Fragile and Robust Subsystems , 2010, PloS one.

[104]  B. Palsson,et al.  Genome-scale Reconstruction of Metabolic Network in Bacillus subtilis Based on High-throughput Phenotyping and Gene Essentiality Data* , 2007, Journal of Biological Chemistry.

[105]  K. Mann,et al.  "Normal" thrombin generation. , 1999, Blood.

[106]  Shuai Li,et al.  ASD v2.0: updated content and novel features focusing on allosteric regulation , 2013, Nucleic Acids Res..

[107]  K. Mann,et al.  Mechanism of factor VIIa-dependent coagulation in hemophilia blood. , 2002, Blood.

[108]  Ali R. Zomorrodi,et al.  Mathematical optimization applications in metabolic networks. , 2012, Metabolic engineering.

[109]  J. H. Matthaei,et al.  Characteristics and stabilization of DNAase-sensitive protein synthesis in E. coli extracts. , 1961, Proceedings of the National Academy of Sciences of the United States of America.

[110]  R J Leipold,et al.  Mathematical Model of Serine Protease Inhibition in the Tissue Factor Pathway to Thrombin (*) , 1995, The Journal of Biological Chemistry.

[111]  John D. Hunter,et al.  Matplotlib: A 2D Graphics Environment , 2007, Computing in Science & Engineering.

[112]  C. Kessler,et al.  Newer concepts of blood coagulation , 1998, Haemophilia : the official journal of the World Federation of Hemophilia.

[113]  K. S. Brown,et al.  Statistical mechanical approaches to models with many poorly known parameters. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[114]  G. T. Tsao,et al.  A cybernetic view of microbial growth: Modeling of cells as optimal strategists , 1985, Biotechnology and bioengineering.

[115]  Julio Saez-Rodriguez,et al.  Training Signaling Pathway Maps to Biochemical Data with Constrained Fuzzy Logic: Quantitative Analysis of Liver Cell Responses to Inflammatory Stimuli , 2011, PLoS Comput. Biol..

[116]  K Fujikawa,et al.  Activation of human blood coagulation factor XI independent of factor XII. Factor XI is activated by thrombin and factor XIa in the presence of negatively charged surfaces. , 1991, The Journal of biological chemistry.

[117]  Jeffrey D. Varner,et al.  Computational Modeling and Analysis of Insulin Induced Eukaryotic Translation Initiation , 2011, PLoS Comput. Biol..

[118]  S. Diamond,et al.  Simulation of Intrathrombus Fluid and Solute Transport Using In Vivo Clot Structures with Single Platelet Resolution , 2013, Annals of Biomedical Engineering.