about dynamical system modelling
play

About dynamical system modelling Illustration in micro-biology - PowerPoint PPT Presentation

About dynamical system modelling Illustration in micro-biology Alain Rapaport alain.rapaport@montpellier.inra.fr Research School AgreenSkills Toulouse, 29-31 october 2014 Contents Modeling bacterial growth in batch and continuous cultures


  1. About dynamical system modelling Illustration in micro-biology Alain Rapaport alain.rapaport@montpellier.inra.fr Research School AgreenSkills Toulouse, 29-31 october 2014

  2. Contents ◮ Modeling bacterial growth in batch and continuous cultures ◮ Simple representations of space ◮ Questions related to biodiversity ◮ Modeling and analysis of growth inhibition ◮ Facing models with experimental data

  3. Bacterial growth in batch number of cells per volume time

  4. Bacterial growth in batch number of cells per volume time

  5. The logistic growth (Verhulst, 1838) x max → x ( t ) = � � x max 1 + e − rt x 0 − 1 Equivalent formulation: x ( · ) is solution of the differential equation � � dx x dt = rx 1 − x max

  6. Jacques Monod’s experiments (1930)

  7. Jacques Monod’s experiments (1930) does not fit the logistic curves...

  8. Mathematical modeling Hypothesis H0. There exists y s.t. x + ys = m = cste � � Hypothesis H1a. dx dt = µ sx ⇒ dx dt = µ m 1 − x x y m ���� r Hypothesis H1b. dx dt = µ ( s ) x where µ ( · ) is not necessarily linear: µ S � How to identify with accuracy this specific growth curve?

  9. The chemostat device pump pump feed bootle culture vessel collection vessel Monod 1950 – Novick & Szilard 1950

  10. Dynamical modelling ◮ Mechanics: 2nd motion law: � mass × acceleration = forces ◮ Bio-reaction kinetics: mass balance (for biotic and abiotic): mass variation = � inputs − � outputs + growth − death � �� � for biotic only

  11. Dynamical modelling ◮ chemical kinetics: → P = A α B β ⇒ v = dP dt = k [ A ] α [ B ] β Ex . : α A + β B − ◮ microbial kinetics: → v = dB [ B ][ S ] k Ex . : B + kS − → B + B − dt = ✘✘✘ ❳❳❳ ✘ ❳ growth is not simply a matter of matching bacteria with molecules of substrate...

  12. The mathematical model (in concentrations)  dx  − Q µ ( s ) x V x dt   = +   − µ ( s ) Q   ds x V ( s in − s ) y dt growth dilution Remark. Q = 0 ⇒ x + ys = constant ˙= d D = Q Simplification and notations. y = 1 dt V s ˙ = − µ ( s ) x + D ( s in − s ) � x ˙ = µ ( s ) x − Dx

  13. Determination of equilibria x ⋆ = 0 µ ( s ⋆ ) = D s = − µ ( s ) x + D ( s in − s ) ˙ ⇒ or s ⋆ = s in x ⋆ = s in − s ⋆ x = µ ( s ) x − Dx ˙ wash-out positive equilibrium µ D S S* S in

  14. Null-clines s = − µ ( s ) x + D ( s in − s ) ˙ x = µ ( s ) x − Dx ˙ x = 0 ˙ s = 0 ˙ x s

  15. Vector field s = − µ ( s ) x + D ( s in − s ) ˙ x = µ ( s ) x − Dx ˙ x = 0 ˙ s = 0 ˙ x s

  16. Phase portrait s = − µ ( s ) x + D ( s in − s ) ˙ x = µ ( s ) x − Dx ˙ x = 0 ˙ s = 0 ˙ x s

  17. Back to experiments For each dilution rate D , one obtains a steady state s ⋆ with µ ( s ⋆ ) = D : µ D 6 D 5 D 4 D 3 D 2 D 1 s s*s* s* s s* s* s* 2 in 1 3 4 5 6

  18. The Monod’s law µ max µ µ ( s ) = µ max s µ K + s : max 2 K S

  19. The chemostat model s = − µ ( s ) x + D ( s in − s ) ˙ x = µ ( s ) x − Dx ˙ Ecology of mountain lakes Industrial bioreactors

  20. For various increasing dilution rates 0 0 0 0 0 0 0 0 0 0 0 0

  21. For various increasing input concentrations 0 0 0 0 0 0 0 0 0 0 0 0

  22. About conversion yield at equilibrium ◮ The mathematical model of the chemostat predicts that the substrate concentration at equilibrium is independant of the input concentration s in (provided that µ ( s in ) > D ). ◮ Micro-biologists report that this property is not verified when the tank is not homogeneous or in natural ecosystems such as soil ecosystems. Question: What is the influence of a spatial repartition on output substrate concentration at steady state?

Recommend


More recommend