stationary states of genome scale metabolic networks in
play

Stationary states of genome scale metabolic networks in continuous - PowerPoint PPT Presentation

Stationary states of genome scale metabolic networks in continuous cell cultures Roberto Mulet Group of Complex Systems and Statistical Physics, Physics Faculty. University of Havana, CUBA With the support of the RISE-H2020 project, 2017-2021


  1. Stationary states of genome scale metabolic networks in continuous cell cultures Roberto Mulet Group of Complex Systems and Statistical Physics, Physics Faculty. University of Havana, CUBA With the support of the RISE-H2020 project, 2017-2021

  2. reservoir e ffl uent culture vessel

  3. Requirements ◮ It must include the internal metabolism ◮ It must include the chemostat ◮ Be computationally scalable to Genome scale metabolic networks ◮ Be flexible ◮ Toxicity ◮ Heterogeneity

  4. Outline Homogeneous chemostat Mathematical framework Stationary States From a Toy model to Genome Scale Heterogeneus chemostat Maximum Entropy Principle The Toy model again Genome Scale Metabolic Network Conclusions

  5. Homogeneous chemostat Mathematical framework Stationary States From a Toy model to Genome Scale Heterogeneus chemostat Maximum Entropy Principle The Toy model again Genome Scale Metabolic Network Conclusions

  6. Chemostat dX dt = ( µ − σ − D ) X (1) µ = µ ( ν ) σ = σ ( s ) (2) ds i dt = − u i X − ( s i − c i ) D (3)

  7. The cell lb k ≤ r k ≤ ub k

  8. The cell lb k ≤ r k ≤ ub k X } = min { V i , c i D − L i ≤ u i ≤ min { V i , c i ξ }

  9. The cell lb k ≤ r k ≤ ub k X } = min { V i , c i D − L i ≤ u i ≤ min { V i , c i ξ } � r k < K k

  10. The cell lb k ≤ r k ≤ ub k X } = min { V i , c i D − L i ≤ u i ≤ min { V i , c i ξ } � r k < K k � S ik r k − e i − y i µ + u i = 0 (4) k

  11. The cell lb k ≤ r k ≤ ub k X } = min { V i , c i D − L i ≤ u i ≤ min { V i , c i ξ } � r k < K k � S ik r k − e i − y i µ + u i = 0 (4) k This is a polytope in very high dimensions

  12. The cell lb k ≤ r k ≤ ub k X } = min { V i , c i D − L i ≤ u i ≤ min { V i , c i ξ } � r k < K k � S ik r k − e i − y i µ + u i = 0 (4) k The cell maximizes biomass production µ Linear Programming LP

  13. Mathematical framework lb k ≤ r q ≤ ub k dX dt = ( µ − σ − D ) X D − L i ≤ u i ≤ min { V i , c i X } � r k < K µ = µ ( ν ) σ = σ ( s ) k � S ik r k − e i − y i µ + u i = 0 k ds i dt = − u i X − ( s i − c i ) D The cell maximizes biomass production LP

  14. Flux Balance � S ik r k − e i − y i µ + u i = 0 k lb k ≤ r k ≤ ub k X } = min { V i , c i D − L i ≤ u i ≤ min { V i , c i ξ } � r i < K i

  15. Flux Balance � S ik r k − e i − y i µ + u i = 0 k lb k ≤ r k ≤ ub k X } = min { V i , c i D − L i ≤ u i ≤ min { V i , c i ξ } � r i < K i The cell maximizes biomass production µ LP u ∗ i ( ξ ) . . . µ ( ξ )

  16. Equilibrium in metabolite’s concentration ds i dt = − u ∗ i X − ( s i − c i ) D

  17. Equilibrium in metabolite’s concentration ds i dt = − u ∗ i X − ( s i − c i ) D s ∗ i ( ξ ) = c i − u ∗ i ( ξ ) ξ

  18. Stationarity in cell’s concentration dX dt = ( µ − σ − D ) X

  19. Stationarity in cell’s concentration dX dt = ( µ − σ − D ) X 0 = ( µ ∗ ( ξ ) − σ ∗ ( ξ ) − D ) X ∗

  20. Stationarity in cell’s concentration dX dt = ( µ − σ − D ) X D = µ ∗ ( ξ ) − σ ∗ ( ξ )

  21. Stationarity in cell’s concentration dX dt = ( µ − σ − D ) X X ∗ = µ ∗ ( ξ ) − σ ∗ ( ξ ) ξ

  22. Stationarity equations r ∗ k . . . u ∗ i ( ξ ) . . . µ ∗ ( ξ ) s ∗ i ( ξ ) = c i − u ∗ i ( ξ ) ξ X ∗ ( ξ ) = µ ∗ ( ξ ) − σ ∗ ( ξ ) ξ

  23. Small Network S E P W E Vazquez et al.. Macromolecular crowding explains overflow metabolism in cells. Scientific Reports 6, 31007 (2016)

  24. Toxicity is the key point bistable regime (a) (b)

  25. General Picture S S E E P W P W E E respiration over fl ow (a) Overflow. At high enough nutrient uptake the respiratory flux hit s the upper bound r max and the remaining nutrients are exported as W . (b) Respiration. The nutrient is completely oxidized with a large energy yield. (c) Threshold values of ξ . ξ 0 delimits the nutrient excess regime ( ξ < ξ 0 ) from the competition regime ( ξ > ξ 0 ). ξ sec delimits the transition between overflow metabolism ( ξ < ξ sec and respiration ( ξ > ξ sec ). Finally, maintenance demand cannot be met beyond ξ > ξ m .

  26. Genome Scale: CHO-K1 line ◮ 6663 reactions ◮ V glc = 0 . 5 mmol /gDW/h ◮ V i = . 1 V glc

  27. Metabolite uptakes and concentration

  28. General picture of the transitions limiting: limiting: gly, tyr, limiting: max. yield nutrient excess ser, asp, pro trp, his, arg, lys, phe glc, gln, asn aas aas glc aas aas aas asp asp glc glc asp glc asp glc asp succ for succ for acald acald for acald lac ala pyr lac for ala ala lac ala

  29. Steady state and bifurcation bistable regime (a) (b) J. Fernandez-de-Cossio Diaz, K. Le´ on and R. M. , Characterizing stationary states of genome scale metabolic networks in continuous culture, PLOS Computational Biology. 13 (11): e1005835 (2017)

  30. Homogeneous chemostat Mathematical framework Stationary States From a Toy model to Genome Scale Heterogeneus chemostat Maximum Entropy Principle The Toy model again Genome Scale Metabolic Network Conclusions

  31. Constraints � S ik r k − e i − y i µ + u i = 0 k lb k ≤ r q ≤ ub k D − L i ≤ u i ≤ min { V i , c i X } � r i < K i We must explore this polytope

  32. Stationarity: Dealing with the heterogeneity d � X dt = ( µ − σ − D ) � X

  33. Stationarity: Dealing with the heterogeneity d � X dt = ( µ − σ − D ) � X 0 = ( µ ( ν ) − σ ( s ) − D )

  34. Stationarity: Dealing with the heterogeneity d � X dt = ( µ − σ − D ) � X 0 = ( µ ( ν ) − σ ( s ) − D ) Effetive Growth rate = µ ( ν ) − σ ( s ) = D

  35. Maximum Entropy Principle If s is fixed, µ ( ν ) − σ ( s ∗ ) = D

  36. Maximum Entropy Principle P s ∗ ( ν ) ∼ e β ( µ ( ν ) − σ ( s ∗ ))

  37. Maximum Entropy Principle P s ∗ ( ν ) ∼ e β ( µ ( ν ) − σ ( s ∗ )) s i = c i − 1 � u a i D a

  38. Maximum Entropy Principle P s ∗ ( ν ) ∼ e β ( µ ( ν ) − σ ( s ∗ )) i = c i − X � s ∗ u i ( ν ) P s ∗ ( ν ) d ν D Π

  39. In short � � � e β ( µ ( ν ) − σ ( s ∗ )) µ ( ν ) − σ ( s ∗ ) D = X Π d ν ξ = � Π d ν e β ( µ ( ν ) − σ ( s ∗ ))

  40. In short � � � e β ( µ ( ν ) − σ ( s ∗ )) µ ( ν ) − σ ( s ∗ ) D = X Π d ν ξ = � Π d ν e β ( µ ( ν ) − σ ( s ∗ )) � s ∗ i = c i − ξ u i ( ν ) P s ∗ ( ν ) d ν Π

  41. Homogeneous vs Heterogeneous Chemostat D = X ξ = � µ ( ν ) − σ ( s ∗ ) � P s ∗ X ∗ ( ξ ) = µ ∗ ( ξ ) − σ ∗ ( ξ ) ξ � s ∗ i = c i − ξ u i ( ν ) P s ∗ ( ν ) d ν s ∗ i ( ξ ) = c i − ξ u ∗ i ( ξ ) Π

  42. Summarizing � S ik r k − e i − y i µ + u i = 0 k lb k ≤ r q ≤ ub k D − L i ≤ u i ≤ min { V i , c i X } D = � µ ( ν ) − σ ( s ∗ ) � P s ∗ � r i < K i = c i − X � s ∗ u i ( ν ) P s ∗ ( ν ) d ν i D Π

  43. Small Network again S E P W E

  44. Effect of the heterogeneity ∞

  45. Effect of the heterogeneity X (10 6 cells/mL) unfeasible 2.0 1.5 stable 1.0 0.5 unstable 0

  46. Genome Scale: CHO-K1 line ◮ 6663 reactions ◮ V glc = 0 . 5 mmol /gDW/h ◮ V i = . 1 V glc

  47. Exploring the space Π � S ik r k − e i − y i µ + u i = 0 k lb k ≤ r q ≤ ub k D − L i ≤ u i ≤ min { V i , c i X } � r i < K i

  48. Exploring the space Π � S ik r k − e i − y i µ + u i = 0 k lb k ≤ r q ≤ ub k D − L i ≤ u i ≤ min { V i , c i X } � r i < K i For β = ∞ : Expectation Propagation Alfredo Braunstein, Anna Paola Muntoni, Andrea Pagnani, An analytic approximation of the feasible space of metabolic networks, Nat. Comm. 8 , 14915 (2017) Here generalized for finite β

  49. Genome Scale Metabolic Networks λ m β = ∞ λ m β =0 λ m β =0 a) b) c) λ m β >970 s glc (mM) s nh4 (mM) λ m β =970 s lac (mM) λ m β = ∞ β =0 λ m β =970 0.01 0.1 1 10 100 1000 0.01 0.1 1 10 100 1000 0.01 0.1 1 10 100 1000 ξ (10 6 cells day/mL) ξ (10 6 cells day/mL) ξ (10 6 cells day/mL)

  50. Genome Scale Metabolic Networks 1.5 λ m β =970 X (10 6 cells/mL) 1.0 λ m β = ∞ 0.5 λ m β =97 0.0 0.5 1.0 1.5 2.0 D (1/day) J. Fernandez-de-Cossio Diaz, and R. M. , Maximum Entropy and Population Heterogeneity in continuos cell cultures , arXiv:1807.04218

  51. Homogeneous chemostat Mathematical framework Stationary States From a Toy model to Genome Scale Heterogeneus chemostat Maximum Entropy Principle The Toy model again Genome Scale Metabolic Network Conclusions

  52. Conclusions ◮ We developed a mathematical framework to determine the stationary states in a chemostat ◮ The presence of toxic waste: ◮ drives the appareance of many stationary states ◮ makes relevant the history of the system ◮ We provided a scheme to estimate the metabolic flux distribution of an heterogeneous culture in a chemostat ◮ The presence of heterogeneity in the culture ◮ changes the concentration of metabolites ◮ allows stationary states with a larger number of cells ◮ Everything is computationally tractable in Genome Scale metabolic networks

Recommend


More recommend