Growth models for cells in the chemostat
|
|
- Clinton Stanley
- 5 years ago
- Views:
Transcription
1 Growth models for cells in the chemostat V. Lemesle, J-L. Gouzé COMORE Project, INRIA Sophia Antipolis BP93, Sophia Antipolis, FRANCE Valerie.Lemesle, Abstract A chemostat is a continuous device where microorganisms consume a nutrient to grow. To describe microorganisms growth in the chemostat, Monod [12] introduced a very simple model. Droop [5] then presented another model considering the possible nutrient storage in a plankton cell. However, this work was mainly based on experimental observations and was not completely justified by biochemical arguments. In this paper, a new model of phytoplankton growth is proposed. Though our model is finally similar to the Droop one, we base our modelling on biochemical arguments. We study two versions of our model considering or not cell mortality. We give also some indications on the observation of such models, and related observers. Keywords biochemical mechanisms, structured model, ordinary differential equations, chemostat, plankton growth models, observers. 1 Introduction Unstructured mathematical models are often employed to describe cell growth in continuous cultures devices, so called chemostats. The most representative one has been proposed by Monod [12]. Its approach is based on very strong assumptions on the interaction between the two main components in the culture: the microorganisms, more precisely bacteria, and the environment. It has been assumed that the biomass concentration alone is able to characterize the physiological state of the whole cell population and that the biomass specific growth rate is controlled only by the concentration of a unique growth limiting nutrient. Although 1
2 this model could be successfully used to fit steady state data, its prediction, when phytoplankton is considered, is far from being satisfactory. Then, a new modelling approach is therefore required to adequately used this continuous culture model when phytoplankton is studied. In fact, structured models seems very efficient for this purpose. One of the most classical one is proposed by Droop [5]. Its approach is based on two steps: the storage and the metabolization of the limiting nutrient for phytoplankton cell. It is assumed that the intracellular limiting component quantity and the total biomass are able to characterize the cell population and that the biomass specific growth rate is given by a function depending on a minimum of internal nutrient quantity. All this assumptions are based on experimental observations. In this paper, we propose another structured model, more based on mechanistic biochemical explanations. Our approach is explained by biochemical mechanisms and not by biological data only. It is assumed that two components are limiting for the growth and that the metabolized and un-metabolized limiting intracellular components characterize the cell population (phytoplankton). This paper is organized as follows. First, we present the chemostat apparatus, recalling the two classical models and their main properties. Then, we explain the basis of our modelling approach and we compare this model to the Droop one. We give finally some indications on building observers for these models. 2 Growth models 2.1 The Monod model The most classical chemostat model has been introduced by Monod to give a synthetic picture of what happens when a microorganism culture (bacteria) is grown in a reactor. Then, only two variables are chosen to describe the various and complex reactions occurring in the reactor vessel: the biomass concentration and the limiting substrate concentration. The mathematical modelling is divided in two parts: the previously mentioned physical part due to the flow and the biological part, which describes the biological reactions in the vessel. This description leads to the well known mass balanced equations: (1) Biological part Physical part Based on biological data,, the specific growth rate, is often defined by a classical function, so-called Michaelis-Menten function.! 2
3 This model can be formulated as a biochemical reaction defined such that: Reactants and product are considered; the rate of reaction depends on biological mechanisms. Indeed, it is assumed some qualitative hypotheses about this specific rate considering some biological characteristics. Hypothesis H 2.1 and is, increasing and bounded. Hypothesis H 2.2 Moreover, the mass conservation principle is verified. Indeed, take the variable (the total mass in the chemostat) then the dynamical equation of follows equation: "! Then, the asymptotic behavior of the system (1) is given by the following proposition. Proposition 2.3 Under the hypotheses (H2.1) and (H2.2), the washout steady state is unstable and the non trivial steady state is globally asymptotically stable (wrt the positive orthant). This simple unstructured model is very intuitive and gives a good idea of the bacteria dynamics. The main problem of this model is that all the physiological characterization is contained in only one variable which is not enough to describe all the different microorganisms. Moreover, this model predicts, for example, that there is no growth if there is no more substrate ( % ). This contradicts experimental facts for phytoplankton. 2.2 The Droop model The Droop model has been built to describe the phytoplankton growth submitted to limiting vitamin B '& in the chemostat [5]. This model is more descriptive than the Monod one: indeed, it distinguishes the evolution of the biomass concentration and of the intracellular limiting nutrient concentration (. The variable ( is socalled the intracellular quota of limiting nutrient. Thus, considering the biological $# 3
4 phenomena and the physical part due to the through passing flow, the following variable yield model is obtained: ( ( ( Biological part Physical part The equation of ( is often replaced by the equation of ( (deduced of (2)): ( ( ( (3) Let us remark that the quantity ( does not makes a difference between the stored and the metabolized nutrient. This model can be illustrated by the following biochemical mechanisms: ( ( storage metabolization The phenomenon of metabolization cannot be described very precisely using the variables (,, : this is the main problem of this model. Indeed, the variable ( does not make the difference between the storage of nutrient and how this storage is used for metabolization. As in the Monod model, Droop assumed some qualitative hypotheses about the specific growth rate ( and storage rate functions based on biological characteristics. Hypotheses H 2.4 ' and is, increasing and bounded. ( is, non negative, increasing and bounded, such that there exists ( defined as follows (. Hypothesis H 2.5 The mass conservation principle is verified; indeed, take the variable ( (the total mass in the chemostat), the dynamical equation of verifies equation: "! The asymptotic behavior is given by the following proposition proved by Lange et al. in [9, 13]. Proposition 2.6 Under all the hypotheses above, the washout equilibrium is unstable and the non trivial steady state is globally asymptotically stable (wrt. the positive orthant). We propose now a structured model described by variables which have a more biological sense and where biochemical phenomena are more taken into account. 4 $# (2)
5 3 A structured phytoplankton growth model In this section, we consider phytoplankton growth limited by the two major components of photo-synthetic cells: nitrogen and carbon. 3.1 Model formulation It is well known that the phytoplanktonic cells grow via absorption of extracellular nitrogen and carbon. We simplify the complex absorption of extracellular nitrogen in two steps: storage and metabolization. The absorption of extracellular carbon is very simple: it is directly metabolized into intracellular organic carbon. The two input limiting components are denoted and. Moreover, to describe the two absorption phenomena, the following variables are chosen: the extracellular nitrogen substrate, the intracellular inorganic nitrogen, the intracellular organic nitrogen, the extracellular carbon substrate, the intracellular organic carbon. We consider the following biochemical mechanisms. storage mortality growth and respiration Let us focus on the mechanism of growth and respiration. The nitrogen metabolization and the carbon metabolization are both represented. More precisely, the nitrogen metabolization is explained by one unit of intracellular inorganic nitrogen metabolized into one unit of intracellular organic nitrogen with the rate of reaction ; the carbon metabolization is characterized by more than one unit of extracellular inorganic carbon metabolized into one unit of intracellular organic carbon with the rate of reaction, the rest being rejected by the respiration with the same rate of reaction. The cell mortality describes the possible cell death before its division; indeed, the dilution rate cannot justify alone the total biomass disappearance. It is modeled by a constant rate ; let us remark that the dynamics of this dead cells (, 5
6 , ) is not defined in this model. In the literature, it can be chosen to add a delay phenomenon in the extracellular inorganic components for the recycling of the dead cells ; indeed, the conversion of organic components into inorganic component is not immediate [2, 3]. Moreover, we take into account another phenomenon, which is not described by the biochemical mechanisms: the conversion of atmospheric carbon into dissolved carbon (extracellular inorganic carbon) with a rate (. Let us remark that this rate can be negative. Adding the physical part due to the flow passing through the chemostat, this description leads to the following model, denoted (S): (4) (5) (6) (7) ( (8) Biological part Physical part The variable represents the available stored nitrogen per organic biomass unit, since, in biology, the usual measure unit of organic biomass is the organic carbon. We make the same hypotheses as for the Monod model (H2.1) on the qualitative behavior of the growth rates. Hypotheses H 3.1 ' and is a increasing bounded function; "! ' and is a increasing bounded function.. Let us remark that the function is not defined for. This is not a problem because this function describe the evolution of the intracellular inorganic components, then has not a biological sense for this. 3.2 Model reduction In this section, we explain why the equations (4), (5), (7) are enough to describe the dynamics of the model (S). First, consider the variable carbon unit and nitrogen unit., means that a cell is composed by one 6
7 The dynamical equation associated to is The variable and have the same asymptotic behavior. More precisely, the equilibrium (i.e. ), is globally asymptotically stable. We deduce that, asymptotically, the organic nitrogen and the organic carbon are linked by a fixed ratio. This fact is often observed at equilibrium on biological data. Thus, the dynamical equation (6) does not give more informations than the one (7). The dynamical equation (9) only depends on the dynamics of (7), (5) and itself; moreover, all the other equations do not rely on the carbon dynamics. Then, in the sequel, we consider the reduced system: (9) 3.3 Interpretation of the limiting nutrient quota The intracellular quota of limiting nutrient is proposed with respect to experimental observations in the Droop model. It represents the intracellular inorganic limiting nutrient quantity per unit of biomass. With the variables of our model, the quota is the total nitrogen quantity in the cell per biomass unit, meaning that is given by (. The dynamical equation of ( is then as follows: ( (10) Compare this equation with the Droop equation (3). The model (2) has been built under equilibria conditions and without mortality phenomenon (i.e. ). If we consider the model (S) under equilibria conditions and with, we obtain the following conditions: and Thus, the quota equation (10) becomes: ( ( ( Though we have built model (9) in a different manner Droop had, we obtain the same formulation for the intracellular quota. Moreover, we are able to explain the 7
8 clear biological meaning of the Droop minimum intracellular quota ( defined by the minimum stored nutrient of a cell to stay alive and estimated from data: it is the meaning of the constant given in the previous subsection. 4 Model analysis The asymptotic behavior of the two versions of the model (9) can be studied: without any mortality rate, or with a fixed positive mortality rate. The second one is very similar to the Droop model. 4.1 Existence of steady states and invariant set First, we give some general properties of the system (9) with. For this system, is not acceptable since the function is not well defined at this point. Then, to make a complete study of this system, we consider the change of variables: which defines a diffeomorphism in the set. We obtain the new system: (11) This new system is well defined in the set Note that it can be seen as an extended system of (9) since it is defined for. Thus, we only consider initial conditions for the state variables belonging to the set which only has a biological sense for (11). Moreover, let us remark that this new system is a loop system [11]: for To ensure the existence of the non trivial equilibrium, we make the following hypotheses. Hypotheses H 4.1 "!. "!. 8
9 Proposition 4.2 Under the hypotheses (H3.1) and (H4.1), there exists two admissible stationary points for the system (11): the washout point and the non trivial equilibrium. Proposition 4.3 Under the hypotheses (H3.1), is invariant for the system (11). 4.2 Local analysis For the local behavior of the system (11), we use the classical method of linearization. Thus, we compute the associated Jacobian of system (11) at a generic point : Let us remark that as the number of non positive off diagonal terms is odd, the system is similar to a competitive system taking [14]. Proposition 4.4 Under the hypotheses (H3.1) and (H4.1), the washout point is unstable. Proposition 4.5 Under the hypotheses (H3.1) and (H4.1), the non trivial equilibrium is locally asymptotically stable. 4.3 Global analysis without a mortality rate We consider the system (11) in the case. We obtain the new system: (12) Proposition 4.6 Under hypotheses (H3.1) and (H4.1), the mass principle is verified and the state variables are bounded. Proof : Consider the total mass in the chemostat. The dynamical equation of verifies equation: $# 9
10 Meaning that: "! #. Thus, the mass principle is verified and the variables are bounded. To prove the global stability of this equilibrium, we use the same techniques as Lange et al. [9, 13] for the study of the classical Droop model. Proposition 4.7 Under hypotheses (H3.1) and (H4.1), the non trivial equilibrium is globally asymptotically stable (wrt. the positive orthant). 5 Observers It is often the case that the quota (either for the Droop model, or for our model) is not measured or is difficult to measure. It can be useful to build observers, that will estimate this quota from the available measurements. The problem here is that some parts of the model can also be not so well known, for example the growth rate. One possibility in this case is to built a kind of observer (sometimes called asymptotic observer, see [1]) independent on this unknown part. Let us suppose that and are the measurements. Then,! verifies the equation and the simple estimator: estimated by. converges towards. The quota can be High gain observers can also be built, if the model is well known (cf. [4]). We will present here an intermediate case (subject of ongoing research) where the model is not well known, but bounds are available. The idea is to build a hybrid observer, similar to a high-gain observer when far of the solution, and switching to an asymptotic observer when closer to the solution (cf. [10]). Let us suppose that the substrate is the only measurement and that the growth rate is such that a bound is known: With the variable, the system is: and an hybrid observer will be: (13) 10
11 & & If, the observer is an asymptotic observer; for, the gains can be adjusted to obtain a high gain observer (for known ). The switch is done when the output error is smaller than. A convergence towards the solution faster than the asymptotic one can be attained. 6 Conclusion In this paper, a model similar to the Droop model has been considered. The new structured model is characterized by biochemical mechanisms which describe the cell storage and how the cell uses this storage to grow. Thus, we obtain a descriptive model using variables which have a more biological sense and which do not come from experimental data only. However, if the same experimental conditions as in the Droop model are considered (equilibria conditions, no mortality rate), then the same formulation for the intracellular quota is obtained. This structured model contains, therefore, more informations than the Droop one. The building of observers for such models is the subject of ongoing research. References [1] G. Bastin and D.Dochain. On-line Estimation and Adaptative Control of Bioreactors, volume 1. Elseiver, Amsterdam, [2] E. Beretta and Y. Takeuchi. Qualitative properties of chemostat equations with time delays: boundedness, local and global asymptotic stability. Differential Equations and Dynamical Systems, 2(1):19 40, [3] E. Beretta and Y. Takeuchi. Qualitative properties of chemostat equations with time delays ii. Differential Equations and Dynamical Systems, 2(4): , [4] O. Bernard, G. Sallet, and A. Sciandra. Nonlinear observers for a class of biological systems. Application to validation of a phytoplanktonic growth model. IEEE Trans. Aut. Control, 43(8): ,
12 [5] M. Droop. Vitamin B12 and marine ecology. IV. The kinetics of uptake, growth and inhibition in Monochrysis Lutheri. Journal of Marine Biology, 48: , [6] J. Guckenheimer and P. Holmes. Nonlinear Oscillations, Dynamical Systems, and Bifurcations of Vector Fields. Springer-Verlag, New York, [7] M. Hirsch and S. Smale. Differential Equations, Dynamical Systems and Linear Algebra. Academic Press, Orlando, [8] J. Hofbauer and K. Sigmund. The theory of evolution and dynamical systems: mathematical aspects of selection. Cambridge University Press, London, [9] K. Lange and F. Oyarzun. The attractiveness of the Droop equations. Mathematical Biosciences, 111: , [10] V. Lemesle and J. L. Gouzé. A two dimensionnal bounded error observer for a class of bioreactor models. In Proceedings of ECC03, Cambridge, September , [11] J. Mallet-Paret and H. Smith. The Poincare-Bendixon Theorem for Monotone Cyclic Feedback system. Journal of Dynamics and Differential Equations, 2: , [12] J. Monod. Recherches sur la Croissance des Cultures Bactériennes. Herman, Paris, [13] F. Oyarzun and K. Lange. The attractiveness of the Droop equations. II. Generic uptake and growth functions. Mathematical Biosciences, 121: , [14] H. L. Smith. Monotone Dynamical Systems. An Introduction to the theory of competitive and cooperative systems. American Mathematical Society, United States of America, [15] M. Vidyasagar. Nonlinear Systems Analysis, second edition. Prentice Hall International, [16] F. Viel, E. Busvelle, and J. P. Gauthier. Stability of polymerization reactors using I/O linearization and a high-gain observer. Automatica, 31: ,
DYNAMICS OF A PREDATOR-PREY INTERACTION IN CHEMOSTAT WITH VARIABLE YIELD
Journal of Sustainability Science Management Volume 10 Number 2, December 2015: 16-23 ISSN: 1823-8556 Penerbit UMT DYNAMICS OF A PREDATOR-PREY INTERACTION IN CHEMOSTAT WITH VARIABLE YIELD SARKER MD SOHEL
More informationTwo simple growth models in the chemostat
Two simple growth models in the chemostat Jean-Luc Gouzé, Valérie Lemesle To cite this version: Jean-Luc Gouzé, Valérie Lemesle. Two simple growth models in the chemostat. Revue Africaine de la Recherche
More informationCompetitive Exclusion in a Discrete-time, Size-structured Chemostat Model
Competitive Exclusion in a Discrete-time, Size-structured Chemostat Model Hal L. Smith Department of Mathematics Arizona State University Tempe, AZ 85287 1804, USA E-mail: halsmith@asu.edu Xiao-Qiang Zhao
More informationI/O monotone dynamical systems. Germán A. Enciso University of California, Irvine Eduardo Sontag, Rutgers University May 25 rd, 2011
I/O monotone dynamical systems Germán A. Enciso University of California, Irvine Eduardo Sontag, Rutgers University May 25 rd, 2011 BEFORE: Santa Barbara, January 2003 Having handed to me a photocopied
More informationModeling Microbial Populations in the Chemostat
Modeling Microbial Populations in the Chemostat Hal Smith A R I Z O N A S T A T E U N I V E R S I T Y H.L. Smith (ASU) Modeling Microbial Populations in the Chemostat MBI, June 3, 204 / 34 Outline Why
More informationDESIGN OF PROBABILISTIC OBSERVERS FOR MASS-BALANCE BASED BIOPROCESS MODELS. Benoît Chachuat and Olivier Bernard
DESIGN OF PROBABILISTIC OBSERVERS FOR MASS-BALANCE BASED BIOPROCESS MODELS Benoît Chachuat and Olivier Bernard INRIA Comore, BP 93, 692 Sophia-Antipolis, France fax: +33 492 387 858 email: Olivier.Bernard@inria.fr
More informationFeedback-mediated oscillatory coexistence in the chemostat
Feedback-mediated oscillatory coexistence in the chemostat Patrick De Leenheer and Sergei S. Pilyugin Department of Mathematics, University of Florida deleenhe,pilyugin@math.ufl.edu 1 Introduction We study
More informationFeedback control for a chemostat with two organisms
Feedback control for a chemostat with two organisms Patrick De Leenheer and Hal Smith Arizona State University Department of Mathematics and Statistics Tempe, AZ 85287 email: leenheer@math.la.asu.edu,
More informationBifurcation Analysis of Continuous Biochemical Reactor Models
Biotechnol. Prog. 2001, 17, 647 660 647 Bifurcation Analysis of Continuous Biochemical Reactor Models Yongchun Zhang and Michael A. Henson* Department of Chemical Engineering, Louisiana State University,
More informationStationary phase. Time
An introduction to modeling of bioreactors Bengt Carlsson Dept of Systems and Control Information Technology Uppsala University August 19, 2002 Abstract This material is made for the course Wastewater
More informationGlobal Stability of Full Open Reversible Michaelis-Menten Reactions
Preprints of the 8th IFAC Symposium on Advanced Control of Chemical Processes The International Federation of Automatic Control Global Stability of Full Open Reversible Michaelis-Menten Reactions Ismail
More informationON THE MATHEMATICAL MODELLING OF MICROBIAL GROWTH: SOME COMPUTATIONAL ASPECTS. Svetoslav M. Markov
Serdica J. Computing 5 (2011), 153 168 Dedicated to the memories of my biology teachers Dr. Kalcho Markov and Dr. Roumen Tsanev ON THE MATHEMATICAL MODELLING OF MICROBIAL GROWTH: SOME COMPUTATIONAL ASPECTS
More informationPolytopic Lyapunov functions for persistence analysis of competing species
Polytopic Lyapunov functions for persistence analysis of competing species Alain Rapaport, Frédéric Grognard, Frédéric Mazenc To cite this version: Alain Rapaport, Frédéric Grognard, Frédéric Mazenc. Polytopic
More informationGlobal stability in a chemostat with multiple nutrients
J. Math. Biol. 2006) Mathematical Biology Digital Obect Identifier DOI): 10.1007/s00285-005-0344-4 Patrick De Leenheer Simon A. Levin Eduardo D. Sontag Christopher A. Klausmeier Global stability in a chemostat
More informationarxiv: v1 [math.ds] 20 Sep 2016
THE PITCHFORK BIFURCATION arxiv:1609.05996v1 [math.ds] 20 Sep 2016 Contents INDIKA RAJAPAKSE AND STEVE SMALE Abstract. We give development of a new theory of the Pitchfork bifurcation, which removes the
More informationDynamical Systems in Biology
Dynamical Systems in Biology Hal Smith A R I Z O N A S T A T E U N I V E R S I T Y H.L. Smith (ASU) Dynamical Systems in Biology ASU, July 5, 2012 1 / 31 Outline 1 What s special about dynamical systems
More informationNetwork Analysis of Biochemical Reactions in Complex Environments
1 Introduction 1 Network Analysis of Biochemical Reactions in Complex Environments Elias August 1 and Mauricio Barahona, Department of Bioengineering, Imperial College London, South Kensington Campus,
More informationDroop models of nutrient-plankton interaction with intratrophic predation
Droop models of nutrient-plankton interaction with intratrophic predation S. R.-J. Jang 1, J. Baglama 2 1. Department of Mathematics, University of Louisiana at Lafayette, Lafayette, LA 754-11 2. Department
More informationFeedback control for chemostat models
J. Math. Biol.46, 48 70 (2003) Mathematical Biology Digital Object Identifier (DOI): 10.1007/s00285-002-0170-x Patrick De Leenheer Hal Smith Feedback control for chemostat models Received: 1 November 2001
More informationExperiment selection for the discrimination of semi-quantitative models of dynamical systems
Artificial Intelligence 170 (2006) 472 506 www.elsevier.com/locate/artint Experiment selection for the discrimination of semi-quantitative models of dynamical systems Ivayla Vatcheva a, Hidde de Jong b,,
More informationDynamics of Modified Leslie-Gower Predator-Prey Model with Predator Harvesting
International Journal of Basic & Applied Sciences IJBAS-IJENS Vol:13 No:05 55 Dynamics of Modified Leslie-Gower Predator-Prey Model with Predator Harvesting K. Saleh Department of Mathematics, King Fahd
More informationGlobal Qualitative Analysis for a Ratio-Dependent Predator Prey Model with Delay 1
Journal of Mathematical Analysis and Applications 266, 401 419 (2002 doi:10.1006/jmaa.2001.7751, available online at http://www.idealibrary.com on Global Qualitative Analysis for a Ratio-Dependent Predator
More informationSimulated niche partitioning by bacteria
Simulated niche partitioning by bacteria Steven S. Andrews and Adam P. Arkin Physical Biosciences Division Lawrence Berkeley National Laboratory 1 Cyclotron Road Berkeley, CA 94720 ssandrews@lbl.gov 1.
More informationObservations on the Stability Properties of Cooperative Systems
1 Observations on the Stability Properties of Cooperative Systems Oliver Mason and Mark Verwoerd Abstract We extend two fundamental properties of positive linear time-invariant (LTI) systems to homogeneous
More informationProblem Set 2. 1 Competitive and uncompetitive inhibition (12 points) Systems Biology (7.32/7.81J/8.591J)
Problem Set 2 1 Competitive and uncompetitive inhibition (12 points) a. Reversible enzyme inhibitors can bind enzymes reversibly, and slowing down or halting enzymatic reactions. If an inhibitor occupies
More informationScience Online Instructional Materials Correlation to the 2010 Biology Standards of Learning and Curriculum Framework
and Curriculum Framework Provider York County School Divison Course Title Biology Last Updated 2010-11 Course Syllabus URL http://yorkcountyschools.org/virtuallearning/coursecatalog.aspx BIO.1 The student
More information2D-Volterra-Lotka Modeling For 2 Species
Majalat Al-Ulum Al-Insaniya wat - Tatbiqiya 2D-Volterra-Lotka Modeling For 2 Species Alhashmi Darah 1 University of Almergeb Department of Mathematics Faculty of Science Zliten Libya. Abstract The purpose
More informationControl of a Bioreactor with Quantized Measurements
Control of a Bioreactor with Quantized Measurements Francis Mairet, Jean-Luc Gouzé To cite this version: Francis Mairet, Jean-Luc Gouzé. Control of a Bioreactor with Quantized Measurements. Fages, François
More informationChapter One. Introduction
Chapter One Introduction With the ever-increasing influence of mathematical modeling and engineering on biological, social, and medical sciences, it is not surprising that dynamical system theory has played
More informationA game-theoretical approach to the analysis of microbial metabolic pathways
A game-theoretical approach to the analysis of microbial metabolic pathways Stefan Schuster Dept. of Bioinformatics Friedrich Schiller University Jena, Germany Introduction Widely used hypothesis: During
More informationObjectives INTRODUCTION TO METABOLISM. Metabolism. Catabolic Pathways. Anabolic Pathways 3/6/2011. How to Read a Chemical Equation
Objectives INTRODUCTION TO METABOLISM. Chapter 8 Metabolism, Energy, and Life Explain the role of catabolic and anabolic pathways in cell metabolism Distinguish between kinetic and potential energy Distinguish
More informationFurther Results on Stabilization of Periodic Trajectories for a Chemostat with Two Species
Further Results on Stabilization of Periodic Trajectories for a Chemostat with Two Species Frédéric Mazenc, Michael Malisoff, and Jérôme Harmand Abstract We discuss an important class of problems involving
More informationA canonical toolkit for modeling plant function
F S P M 0 4 A canonical toolkit for modeling plant function M. Renton, J. Hanan, K. Burrage ACMC, University of Queensland, Australia Introduction From seeds, forms emerge, growing and evolving and interacting
More informationDynamics on a General Stage Structured n Parallel Food Chains
Memorial University of Newfoundland Dynamics on a General Stage Structured n Parallel Food Chains Isam Al-Darabsah and Yuan Yuan November 4, 2013 Outline: Propose a general model with n parallel food chains
More informationScience Textbook and Instructional Materials Correlation to the 2010 Biology Standards of Learning and Curriculum Framework. Publisher Information
Publisher Information Copyright date 2013 Contact Carol Kornfeind Phone# 847-486-2065 E-mail carol.kornfeind@pearson.com Biology 1 of 12 Virginia Department of Education Text Miller Levine Biology, Virginia
More informationApplication demonstration. BifTools. Maple Package for Bifurcation Analysis in Dynamical Systems
Application demonstration BifTools Maple Package for Bifurcation Analysis in Dynamical Systems Introduction Milen Borisov, Neli Dimitrova Department of Biomathematics Institute of Mathematics and Informatics
More informationChemotherapy in heterogeneous cultures of cancer cells with interconversion
Chemotherapy in heterogeneous cultures of cancer cells with interconversion Rui DILãO Institut des Hautes Études Scientifiques 35, route de Chartres 91440 Bures-sur-Yvette (France) Juin 2013 IHES/M/13/20
More informationHOMOTOPY METHODS FOR COUNTING REACTION NETWORK GHEORGHE CRACIUN, J. WILLIAM HELTON, AND RUTH J. WILLIAMS. August 26, 2008
HOMOTOPY METHODS FOR COUNTING REACTION NETWORK EQUILIBRIA GHEORGHE CRACIUN, J. WILLIAM HELTON, AND RUTH J. WILLIAMS August 26, 2008 Abstract. Dynamical system models of complex biochemical reaction networks
More informationA First Course on Kinetics and Reaction Engineering. Class 9 on Unit 9
A First Course on Kinetics and Reaction Engineering Class 9 on Unit 9 Part I - Chemical Reactions Part II - Chemical Reaction Kinetics Where We re Going A. Rate Expressions - 4. Reaction Rates and Temperature
More informationQualitative analysis of a chemostat model with inhibitory exponential substrate uptake q
Chaos, Solitons and Fractals 23 (2005) 873 886 www.elsevier.com/locate/chaos Qualitative analysis of a chemostat model with inhibitory exponential substrate uptake q Guifang Fu a, *, Wanbiao Ma a, Shigui
More informationLecture Notes in Biomathematics
r 6 Lecture Notes in Biomathematics Managing Editor: S. Levin 68 The Dynamics of Physiologically Structured Populations Edited by J.A.J. Metz arid O. Diekmann Springer-Verlag Berlin Heidelberg New York
More informationMonotone Control System. Brad C. Yu SEACS, National ICT Australia And RSISE, The Australian National University June, 2005
Brad C. Yu SEACS, National ICT Australia And RSISE, The Australian National University June, 005 Foreword The aim of this presentation is to give a (primitive) overview of monotone systems and monotone
More informationNonlinear Control Lecture 1: Introduction
Nonlinear Control Lecture 1: Introduction Farzaneh Abdollahi Department of Electrical Engineering Amirkabir University of Technology Fall 2011 Farzaneh Abdollahi Nonlinear Control Lecture 1 1/15 Motivation
More informationRosenzweig-MacArthur Model. Considering the Function that Protects a Fixed. Amount of Prey for Population Dynamics
Contemporary Engineering Sciences, Vol. 11, 18, no. 4, 1195-15 HIKAI Ltd, www.m-hikari.com https://doi.org/1.1988/ces.18.8395 osenzweig-macarthur Model Considering the Function that Protects a Fixed Amount
More informationAn Evolutionary Game Perspective to ALOHA with power control
An Evolutionary Game Perspective to ALOHA with power control E. Altman a, N. Bonneau a, M. Debbah b and G. Caire b a MAESTRO, INRIA Sophia-Antipolis, 004 Route des Lucioles, B.P.93, 0690 Sophia-Antipolis,
More informationA Producer-Consumer Model With Stoichiometry
A Producer-Consumer Model With Stoichiometry Plan B project toward the completion of the Master of Science degree in Mathematics at University of Minnesota Duluth Respectfully submitted by Laura Joan Zimmermann
More informationLyapunov Stability Theory
Lyapunov Stability Theory Peter Al Hokayem and Eduardo Gallestey March 16, 2015 1 Introduction In this lecture we consider the stability of equilibrium points of autonomous nonlinear systems, both in continuous
More informationDIFFERENTIAL EQUATIONS, DYNAMICAL SYSTEMS, AND AN INTRODUCTION TO CHAOS
DIFFERENTIAL EQUATIONS, DYNAMICAL SYSTEMS, AND AN INTRODUCTION TO CHAOS Morris W. Hirsch University of California, Berkeley Stephen Smale University of California, Berkeley Robert L. Devaney Boston University
More informationPractical exercises INSA course: modelling integrated networks of metabolism and gene expression
Practical exercises INSA course: modelling integrated networks of metabolism and gene expression Hidde de Jong October 7, 2016 1 Carbon catabolite repression in bacteria All free-living bacteria have to
More informationMetabolic diversity is based on the Electron donors, acceptors, and carbon sources available - thermodynamics
To date you have covered microbial community sampling using molecular techniques to identify who is present in the environment. You have also looked at various genetic mechanisms to understand how organisms
More informationComputer Simulation and Data Analysis in Molecular Biology and Biophysics
Victor Bloomfield Computer Simulation and Data Analysis in Molecular Biology and Biophysics An Introduction Using R Springer Contents Part I The Basics of R 1 Calculating with R 3 1.1 Installing R 3 1.1.1
More informationProceedings of the FOSBE 2007 Stuttgart, Germany, September 9-12, 2007
Proceedings of the FOSBE 2007 Stuttgart, Germany, September 9-12, 2007 A FEEDBACK APPROACH TO BIFURCATION ANALYSIS IN BIOCHEMICAL NETWORKS WITH MANY PARAMETERS Steffen Waldherr 1 and Frank Allgöwer Institute
More informationSummer School on Mathematical Control Theory. Control of biological processes
united nations educational, scientific and cultural organization the i international centre for theoretical physics international atomic energy agency SMR1 327/25 Summer School on Mathematical Control
More informationOrdinary Differential Equations
Ordinary Differential Equations Michael H. F. Wilkinson Institute for Mathematics and Computing Science University of Groningen The Netherlands December 2005 Overview What are Ordinary Differential Equations
More informationAnalysis of Plankton Populations
Analysis of Plankton Populations Sean Li September 21, 25 Abstract I use numerical and analytical techniques to study the dynamics of various models of plankton food webs with and without resource fluctuations.
More informationSystems Biology in Photosynthesis INTRODUCTION
1 / 26 Systems Biology in Photosynthesis INTRODUCTION Rainer Machné Institute for Theoretical Chemistry, University of Vienna, Austria PSI - Photon System Instruments, Czech Republic Brno, April, 2011
More informationEQUILIBRIA AND STABILITY ANALYSIS OF A BRANCHED METABOLIC NETWORK WITH FEEDBACK INHIBITION. F. Grognard Y. Chitour G. Bastin
EQUILIBRIA AND STABILITY ANALYSIS OF A BRANCHED METABOLIC NETWORK WITH FEEDBACK INHIBITION F. Grognard Y. Chitour G. Bastin Projet COMORE. INRIA Sophia-Antipolis. BP 93 06902 Sophia-Antipolis Cedex, France
More informationBIBO STABILITY AND ASYMPTOTIC STABILITY
BIBO STABILITY AND ASYMPTOTIC STABILITY FRANCESCO NORI Abstract. In this report with discuss the concepts of bounded-input boundedoutput stability (BIBO) and of Lyapunov stability. Examples are given to
More informationStabilization of a 3D Rigid Pendulum
25 American Control Conference June 8-, 25. Portland, OR, USA ThC5.6 Stabilization of a 3D Rigid Pendulum Nalin A. Chaturvedi, Fabio Bacconi, Amit K. Sanyal, Dennis Bernstein, N. Harris McClamroch Department
More informationOn the status of the Michaelis-Menten equation and its implications for enzymology
1 On the status of the Michaelis-Menten equation and its implications for enzymology Sosale Chandrasekhar 1 Department of Organic Chemistry, Indian Institute of Science, Bangalore 560 012, India 1 E-mail:
More information11. S. Jang, Dynamics of a discrete host-parasitoid system with stocking, Discrete Dynamics
Sophia Jang Department of Mathematics and Statistics Texas Tech University Office: MA202 Phone: (806) 834-7006 Fax: (806) 742-1112 E-mail: sophia.jang@ttu.edu Publications 1. M. De Silva, S. Jang, Period-doubling
More informationBasic Synthetic Biology circuits
Basic Synthetic Biology circuits Note: these practices were obtained from the Computer Modelling Practicals lecture by Vincent Rouilly and Geoff Baldwin at Imperial College s course of Introduction to
More informationChapter 6 # METABOLISM PowerPoint Image Slideshow
COLLEGE BIOLOGY PHYSICS Chapter 6 # METABOLISM Chapter Title PowerPoint Image Slideshow Figure 8.1 Metabolism Figure 6.2 Energy from the sun. Plants photosynthesis Herbivores eat those plants Carnivores
More informationAnalysis and Simulation of Biological Systems
Analysis and Simulation of Biological Systems Dr. Carlo Cosentino School of Computer and Biomedical Engineering Department of Experimental and Clinical Medicine Università degli Studi Magna Graecia Catanzaro,
More informationDIFFERENTIAL EQUATIONS, DYNAMICAL SYSTEMS, AND AN INTRODUCTION TO CHAOS
DIFFERENTIAL EQUATIONS, DYNAMICAL SYSTEMS, AND AN INTRODUCTION TO CHAOS Morris W. Hirsch University of California, Berkeley Stephen Smale University of California, Berkeley Robert L. Devaney Boston University
More informationPredator-Prey Model with Ratio-dependent Food
University of Minnesota Duluth Department of Mathematics and Statistics Predator-Prey Model with Ratio-dependent Food Processing Response Advisor: Harlan Stech Jana Hurkova June 2013 Table of Contents
More informationB L U E V A L L E Y D I S T R I C T C U R R I C U L U M Science 7 th grade
B L U E V A L L E Y D I S T R I C T C U R R I C U L U M Science 7 th grade ORGANIZING THEME/TOPIC UNIT 1: CELLS Structure and Function of Cells MS-LS1-1. Conduct an investigation to provide evidence that
More informationAn Introduction to Metabolism
LECTURE PRESENTATIONS For CAMPBELL BIOLOGY, NINTH EDITION Jane B. Reece, Lisa A. Urry, Michael L. Cain, Steven A. Wasserman, Peter V. Minorsky, Robert B. Jackson Chapter 8 An Introduction to Metabolism
More informationCANONICAL LOSSLESS STATE-SPACE SYSTEMS: STAIRCASE FORMS AND THE SCHUR ALGORITHM
CANONICAL LOSSLESS STATE-SPACE SYSTEMS: STAIRCASE FORMS AND THE SCHUR ALGORITHM Ralf L.M. Peeters Bernard Hanzon Martine Olivi Dept. Mathematics, Universiteit Maastricht, P.O. Box 616, 6200 MD Maastricht,
More informationBiological Pathways Representation by Petri Nets and extension
Biological Pathways Representation by and extensions December 6, 2006 Biological Pathways Representation by and extension 1 The cell Pathways 2 Definitions 3 4 Biological Pathways Representation by and
More informationAn Introduction to Metabolism. Chapter 8
An Introduction to Metabolism Chapter 8 METABOLISM I. Introduction All of an organism s chemical reactions Thousands of reactions in a cell Example: digest starch use sugar for energy and to build new
More informationContinuous time population models
Continuous time population models Jaap van der Meer jaap.van.der.meer@nioz.nl Abstract Many simple theoretical population models in continuous time relate the rate of change of the size of two populations
More informationOn predator-prey systems and small-gain theorems
On predator-prey systems and small-gain theorems Patrick De Leenheer, David Angeli and Eduardo D. Sontag December 4, Abstract This paper deals with an almost global attractivity result for Lotka-Volterra
More informationEcological Modelling
Ecological Modelling 243 (212 81 88 Contents lists available at SciVerse ScienceDirect Ecological Modelling jo ur n al homep ag e: www.elsevier.com/locate/ecolmodel On the strict homeostasis assumption
More informationTHE DESIGN of stabilizing feedback control laws for
IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 42, NO. 4, APRIL 1997 473 Robust Feedback Stabilization of Chemical Reactors F. Viel, F. Jadot, and G. Bastin Abstract This paper deals the temperature stabilization
More informationChapter 8: An Introduction to Metabolism. 1. Energy & Chemical Reactions 2. ATP 3. Enzymes & Metabolic Pathways
Chapter 8: An Introduction to Metabolism 1. Energy & Chemical Reactions 2. ATP 3. Enzymes & Metabolic Pathways 1. Energy & Chemical Reactions 2 Basic Forms of Energy Kinetic Energy (KE) energy in motion
More informationPopulation Dynamics in a Microfounded Predator-Prey Model. Thomas Christiaans. University of Siegen. Discussion Paper No
Volkswirtschaftliche Diskussionsbeiträge V W L Population Dynamics in a Microfounded Predator-Prey Model Thomas Christiaans University of Siegen Discussion Paper No. 8-04 ISSN 433-058x UNIVERSITÄT SIEGEN
More informationLimit Cycles in High-Resolution Quantized Feedback Systems
Limit Cycles in High-Resolution Quantized Feedback Systems Li Hong Idris Lim School of Engineering University of Glasgow Glasgow, United Kingdom LiHonIdris.Lim@glasgow.ac.uk Ai Poh Loh Department of Electrical
More informationEnrichment in a General Class of Stoichiometric Producer-Consumer Population Growth Models
Enrichment in a General Class of Stoichiometric Producer-Consumer Population Growth Models Harlan Stech a, c, Bruce Peckham a and John Pastor b a Department of Mathematics and Statistics, University of
More informationOptimization of a photobioreactor biomass production using natural light
INSTITUT NATIONAL DE RECHERCHE EN INFORMATIQUE ET EN AUTOMATIQUE arxiv:1009.2137v1 [math.oc] 11 Sep 2010 Optimization of a photobioreactor biomass production using natural light F. Grognard A.R. Akhmetzhanov
More informationAn Introduction to Metabolism
Chapter 8 1 An Introduction to Metabolism PowerPoint Lecture Presentations for Biology Eighth Edition Neil Campbell and Jane Reece Lectures by Chris Romero, updated by Erin Barley with contributions from
More informationChapter 6- An Introduction to Metabolism*
Chapter 6- An Introduction to Metabolism* *Lecture notes are to be used as a study guide only and do not represent the comprehensive information you will need to know for the exams. The Energy of Life
More informationCellular Metabolic Models
Cellular Metabolic Models. Cellular metabolism. Modeling cellular metabolism. Flux balance model of yeast glycolysis 4. Kinetic model of yeast glycolysis Cellular Metabolic Models Cellular Metabolism Basic
More informationResearch Article Global Dynamics of a Competitive System of Rational Difference Equations in the Plane
Hindawi Publishing Corporation Advances in Difference Equations Volume 009 Article ID 1380 30 pages doi:101155/009/1380 Research Article Global Dynamics of a Competitive System of Rational Difference Equations
More informationNetworks in systems biology
Networks in systems biology Matthew Macauley Department of Mathematical Sciences Clemson University http://www.math.clemson.edu/~macaule/ Math 4500, Spring 2017 M. Macauley (Clemson) Networks in systems
More informationExtremum seeking via continuation techniques for optimizing biogas production in the chemostat
9th IFAC Symposium on Nonlinear Control Systems Toulouse, France, September 4-6, 2013 WeB2.1 Extremum seeking via continuation techniques for optimizing biogas production in the chemostat A. Rapaport J.
More informationCHEM-E3205 BIOPROCESS OPTIMIZATION AND SIMULATION
CHEM-E3205 BIOPROCESS OPTIMIZATION AND SIMULATION TERO EERIKÄINEN ROOM D416d tero.eerikainen@aalto.fi COURSE LECTURES AND EXERCISES Week Day Date Time Place Lectures/Execises 37 Mo 12.9.2016 10:15-11:45
More informationStability Analysis of Nonlinear Hybrid System in Microbial Fed-batch Culture
The Fourth International Conference on Computational Systems Biology (ISB010) Suzhou, China, September 9 11, 010 Copyright 010 ORSC & APORC, pp. 363 37 Stability Analysis of Nonlinear Hybrid System in
More informationFeedback stabilisation with positive control of dissipative compartmental systems
Feedback stabilisation with positive control of dissipative compartmental systems G. Bastin and A. Provost Centre for Systems Engineering and Applied Mechanics (CESAME Université Catholique de Louvain
More informationA Note on the Ramsey Growth Model with the von Bertalanffy Population Law
Applied Mathematical Sciences, Vol 4, 2010, no 65, 3233-3238 A Note on the Ramsey Growth Model with the von Bertalanffy Population aw uca Guerrini Department of Mathematics for Economic and Social Sciences
More informationMathematical Modeling of Competition for Light and Nutrients Between Phytoplankton Species in a Poorly Mixed Water Column
University of Wisconsin Milwaukee UWM Digital Commons Theses and Dissertations May 2014 Mathematical Modeling of Competition for Light and Nutrients Between Phytoplankton Species in a Poorly Mixed Water
More information1. (rate of) {production of / energy incorporated into / eq} {biomass / organic material / organic molecules / tissue} ;
1(a)(i) 1. (rate of) {production of / energy incorporated into / eq} {biomass / organic material / organic molecules / tissue} ; 2. reference to {losses in respiration / GPP- R } ; 3. in {producers / plants
More informationA LaSalle version of Matrosov theorem
5th IEEE Conference on Decision Control European Control Conference (CDC-ECC) Orlo, FL, USA, December -5, A LaSalle version of Matrosov theorem Alessro Astolfi Laurent Praly Abstract A weak version of
More informationVITA EDUCATION. Arizona State University, , Ph.D. Huazhong Univ. of Sci. and Tech., , M.S. Lanzhou University, , B.S.
VITA Bingtuan Li Department of Mathematics University of Louisville Louisville, KY 40292 Office phone: (502) 852-6149 Fax: (502) 852-7132 E-mail: bing.li@louisville.edu EDUCATION Arizona State University,
More informationA dynamic Model coupling Photoacclimation and Photoinhibition in Microalgae
A dynamic Model coupling Photoacclimation and Photoinhibition in Microalgae Philipp Hartmann, Andrea Nikolaou, Benoît Chachuat, Olivier Bernard To cite this version: Philipp Hartmann, Andrea Nikolaou,
More informationInflow Qin, Sin. S, X Outflow Qout, S, X. Volume V
UPPSALA UNIVERSITET AVDELNINGEN FÖR SYSTEMTEKNIK BC,PSA 9809, Last rev August 17, 2000 SIMULATION OF SIMPLE BIOREACTORS Computer laboratory work in Wastewater Treatment W4 1. Microbial growth in a "Monode"
More informationSemi-analytical solutions for cubic autocatalytic reaction-diffusion equations; the effect of a precursor chemical
ANZIAM J. 53 (EMAC211) pp.c511 C524, 212 C511 Semi-analytical solutions for cubic autocatalytic reaction-diffusion equations; the effect of a precursor chemical M. R. Alharthi 1 T. R. Marchant 2 M. I.
More informationBioinformatics 3. V18 Kinetic Motifs. Fri, Jan 8, 2016
Bioinformatics 3 V18 Kinetic Motifs Fri, Jan 8, 2016 Modelling of Signalling Pathways Curr. Op. Cell Biol. 15 (2003) 221 1) How do the magnitudes of signal output and signal duration depend on the kinetic
More informationSelf-Organization in Nonequilibrium Systems
Self-Organization in Nonequilibrium Systems From Dissipative Structures to Order through Fluctuations G. Nicolis Universite Libre de Bruxelles Belgium I. Prigogine Universite Libre de Bruxelles Belgium
More informationIntroduction on metabolism & refresher in enzymology
Introduction on metabolism & refresher in enzymology Daniel Kahn Laboratoire de Biométrie & Biologie Evolutive Lyon 1 University & INRA MIA Department Daniel.Kahn@univ-lyon1.fr General objectives of the
More information