Team:SydneyUni Australia/Modelling Principles

From 2013.igem.org

Revision as of 05:18, 28 October 2013 by Robbie.oppenheimer (Talk | contribs)

SydneyUniversity Top Banner.jpg SydneyUniversity Bottom Banner.jpg



A Schematic of the Engineered Metabolic Pathway:

Pathway RRR.png


The above figure is a simplified schematic of the two metabolic pathways we have modelled. The symbols in squares signify the intracellular concentrations of the associated metabolite (Greek letters) or enzyme (capital English letters). These symbols are used throughout the analysis.


General Information regarding the Enzymes Involved in the Metabolic Pathway

Enzyme Gene Symbol Constants Substrate Product Ref
1,2-Dichloroethane Dechlorinase dhlA A KM A = 0.53 mM, kcat A = 3.3 s-1 DCA 2-chloroethanol [1]
Alcohol Dehydrogenase Adhlb 1/2* B KM B = 0.94 mM, kcat B = 0.0871 s-1 2-chloroethanol chloroacetaldehyde [2]
p450 p450 C KM C = 0.120 mM, kcat C = 0.0113 s-1 DCA chloroacetaldehyde [3]
Chloroacetaldehyde Dehydrogenase aldA D KM D = 0.06mM, kcat D = 0.60 s-1 chloroacetaldehyde chloroacetate [4]
Haloacetate Dehydrogenase dhlB E KM E = 20 mM, kcat E = 25.4 s-1 chloroacetate glycolate [5]

A table that summarises the enzymes that we used in both of the engineered metabolic pathways. The constants presented were (painstakingly) obtained from the literature (referenced).


ODE Model

Shown below are the two systems of ODEs we used to model each pathway.


The non-monooxygenase pathway (with dhlA (A) and adh1b1*2 (B)).

Igem ode 111.png


The monooxygenase pathway (with p450 (C)).

Igem ode 22.png

Generally, each line represents how the metabolic intermediate changes over time; this is dictated by the relative rate at which the intermediate is formed vs. the rate at which it is used/removed.

The two above systems of ODEs were used to model two things:

1. The rate at which extracellular DCA is removed from solution (given the number of cells in solution) . This is considering the rate at which DCA crosses the plasma membrane (from solution to cytoplasm) of all cells in solution.

2. The rate at which the intracellular concentrations of the metabolic intermediates change over time within a single cell. This considers the relative activities of the enzymes of the metabolic pathway.


The symbols A, B, C, D & E represent the intracellular enzyme concentrations of 1,2-Dichloroethane Dechlorinase, Alcohol Dehydrogenase, Cytochrome p450, Chloroacetaldehyde Dehydrogenase and Haloacetate Dehydrogenase respectively. They have units of mM.


The symbols αin, β, γ, δ and ε represent the intracellular concentrations of the metabolic intermediates DCA, 2-chloroethanol, chloroacetaldehyde, chloroacetate and glycolate respectively. They have units mM. αout represents the extracellular concentration of DCA, i.e. the concentration of DCA in solution. It has units mM.


The symbol Ψ represents the concentration of cells in solution and has units cells/mL.


The function J(αoutin) represents the rate which extracellular DCA moves across the plasma membrane and into the single cell (the flux rate). It is a function of the extra and intracellular concentrations of DCA and has a value of surface area per time m2s-1. The value S represents the average surface area of the cellular membrane of E. coli. By multiplying S by J one achieves the total amount of DCA flowing into a single cell per unit time. The total amount of DCA flowing from solution is achieved by multiplying Ψ by S and J.


Constructing a model for the metabolic pathway

The rate of change for the intracellular concentration of metabolites (β, γ, δ and ε) is purely determined through the relative activity (described by MM equations) of the enzymes creating or removing the metabolite.

Luckily the enzymes of our metabolic pathway are accurately described by MM kinetics and all constants were available in the literature. The system is based on the forward and reverse kinetic rate (k1 and k-1 respectively) of substrate-enzyme binding followed by the irreversible catalytic step of product formation (k2) as depicted in the figure below:

Igem pathway 1.png

The rate of catalysis for an enzyme is a function of the enzyme and substrate concentration and incorporates the catalytic constant kcat and the Michaelis constant KM, which are given in most enzyme kinetics studies.

Igem Re.jpg

The symbol RE denotes the reaction rate: i.e. the rate of product formation (which is equal to the rate of substrate removal).

Igem Re= dproduct.png

So, by observing:

i) that the rate at which a substrate is removed is directly equal to the rate at which the associated product is created (e.g. the rate at which γ is removed is the rate at which δ is formed) ii) and how each intermediate simultaneously acts as a substrate and product (e.g. how D creates δ, and E uses δ to create ε)

it is easy to see how this process gives rise to a connected system of ODEs.


Estimating Protein Concentration

The accuracy of this model is strongly dependent on the intracellular enzyme concentrations; these cannot be determined theoretically as there is a weak correlation with the level of gene expression and protein abundance.

Ishihama et al. ([11]) profiled the protein concentration in E. coli as:

Why was ET brown?.png

The enzyme-encoding-genes are under the regulation of an artificial promoter termed Psyn. The promoter is designed to be constitutive and drive high expression of the enzymes. We estimated the protein concentration as 10,000 protein units per cell - the upper quartile limit in the ‘Highly abundant group’ box in the figure above.

By taking the cell volume as 0.65µ3 = 0.63 x 10-18 m3 [7]

Pikachu is simply a mouse that sneezes electricity.png

So the intracellular concentration for all proteins is predicted to be 25.55 mM


DCA Diffusion Across the Plasma Membrane

We were presented with a scenario where it was necessary to model how DCA would move from the solution into the cell where it’s metabolised: i.e. linking αout and αin.

This was tackled by modelling the rate of DCA diffusion across a cell membrane through Fick’s first law of diffusion:

J=p.png

Fick’s first law of diffusion is justified since 1,2-DCA is non-polar and the cellular membrane is thin. The law states that the flux, J, of DCA across the membrane is equal to the permeability coefficient, P, times the concentration difference of DCA across the cell membrane. The flux has units m2 s-1 The permeability coefficient of DCA is not reported in the literature, so it had to be estimated from the definition of the permeability constant.


P=kow.png

Where D is the diffusion constant for DCA diffusion across the plasma membrane and d is the length of the membrane. The partition coefficient for DCA across a plasma membrane is not documented, but can be estimated to be similar to the octanol-water partition coefficient (Kow). This constant determines the equilibrium ratio of DCA in octanol and water. Like many properties of DCA, the diffusion constant, D, is not documented in the literature. It was determined through the famous Stokes-Einstein equation:


D=kb.png

Where kB, T, η & r represent the Boltzmann constant (1.3806488 × 10-23 m2 kg s-2 K-1), temperature, viscosity of the membrane and radius of DCA respectively. Here we must assume DCA is spherical. The ‘radius’ of DCA isn’t described in the literature. But the van der Waal constant, b, can be used to calculate the van der Waal volume, VW , and hence the van der Waal radius, rW3. Note DCA is not spherical, and this method is used to calculate atoms.

Igem b=.png

By combining all of the above, the flux can be described as:

Igem summary J=.png


Summary of Constants
Symbol Name Value (units) Ref
η Cellular Membrane Viscosity 1.9 kg/m/s [6]
S E. coli Membrane Surface Area 6x10-12 m2 [7]
r DCA radius 0.3498 nm [8]
Kow Octanol-water Partition Coefficient for DCA 28.2 [9]
d Length of Cellular Membrane 2 nm [10]

Extending the System to Cell Cultures

The model describes the process of DCA diffusion and metabolism of a single cell. So given a concentration of our engineered cells in solution one can calculate the rate at which the DCA concentration in solution (ααout) decreases by simply multiply the flux rate of DCA across the membrane by the total concentration of cells in the solution. So, by letting Ψ represent the concentration of cells in solution. The overall rate of decrease of DCA in solution is:

Igem datot=daout1213.png

Since the latter portion of the system of ODEs describes the intracellular concentrations within a single cell, the rate at which the DCA comes into the cell and is metabolised (dαin/dt) is simply:

Simply is a simply simple word.png


The motivation for doing this as it seems more natural to consider the intracellular metabolite concentrations of a single cell, rather than averaged across many. It also allows one to gauge whether chloroacetaldehyde reaches cytotoxic levels. The assumption made here is that all the cells in the solution are exactly the same.

Factoring in Cell Growth

The cells are expected to grow due to the production of glycolate (used as a carbon source for growth). So one can model the rate at which the cell concentration, Ψ, increases. Ideally one would want to experimentally determine the cell growth as function of time, but unfortunately we did not have enough time to do this. We turned to other, more approximate measures: The growth of E. coli over time due to the presence of glycolate was initially described by Lord [12]:


Ned flanders.png

Firstly it must be noted that this analysis is highly approximate, and is used only to obtain a very rough estimate of the growth rate as a function of glycolate produced. By using the approximation that an OD550 of 0.1 = 1 x 108 cells/mL one can use the above graph to generate rates of E. coli growth due to glycolate production. From graph 1, it can be seen that a 2.2 x 108cells us 3.33 mM of glycolate in 95 seconds. Or in other words, 1 cell uses 1.5789 x 10-10 mM of glycolate per second:

Jakethedog.png

This assumes that only the original cells existing in solution from the graphs above (the 2.2E8 cells) use glycolate for growth. From the graph b, the OD550 increases at a rate of 0.25 per 50 seconds (gradient of graph 2) which shows that the cells increase at a maximum rate of 5x106 cells / mL / s when growing in saturated glycolate. From the graphs, it can be seen that the point where maximal growth no longer occurs (where the straight line becomes curved) is when glycolate is at a very low concentration. We will take it that the saturation point is so small that it is negligible – i.e. the growth rate is always maximal when in the presence of glycolate. So the cellular growth rate can be described as:

Finthehuman.png

Using this approach, the model for growth rate is most appropriate for cell concentrations in the range of 2 x 108 to 6 x 108 cells/mL.


With thanks to: