Team:SJTU-BioX-Shanghai/Project/project3

From 2012.igem.org

Revision as of 01:27, 27 September 2012 by Huanan1991 (Talk | contribs)

Acceleration Modeling

  • State of Art

In mathematics, a partial differential equation (PDE) is a differential equation that contains unknown multivariable functions and their partial derivatives. PDEs are used to formulate problems involving functions of several variables, and are either solved by hand, or used to create a relevant computer model.

Mathematical modeling is a useful tool for understanding how a bio-system works and how to improve it. We build a mathematical model mainly based on the mechanism of the involved biochemical reaction. Here we focus on the dynamic process of membrane protein cluster in order to reveal how the system works over time.

  • Aim

PDEs can be used to describe a wide variety of phenomena such as sound, heat, electrostatics, electrodynamics, fluid flow, or elasticity. Defined by the PDE, the variable describing the phenomena can be solved according to different boundary values.

  • Achievements

Based on diffusion equation, by solving it under different boundary values, we found that the cluster system can increase production significantly.

Basic Assumptions

1. Protein clusters locate on the membrane, which means the products of our device diffuse in an area of hemisphere. We assume that the membrane has no extra effect on the behavior of molecular diffusion. So the concentration will be doubled compared to the diffusion progress in a volume of sphere.
2. The primary substrate concentration remains the same in the cell cytoplasm.
3. The products of our cluster enzyme diffuse as the diffusion equation. The diffusion equation is a partial differential equation which describes density dynamics in a material undergoing diffusion. It is also used to describe processes exhibiting diffusive-like behavior, which means the diffusion rate is proportional to the gradient of density.
4. Each our membrane device occupies an area of membrane, adjacent to the others. Because each device is the same as others, the diffusion of the products out of the area can be regarded as reflected by the imaginary barrier, which means the normal derivative at the edge of the area is 0.
5. The enzyme can be regarded as a thermal source, producing the molecular at a rate of v.

Acceleration Model

Diffusion Model

The diffusion area is regarded as a hemisphere which radius is R with an enzyme located on the center which occupies a region of r0. So we have the equaction:

12SJTU formula1.png

While:

12SJTU formula11.png

By solving the equations we can show how the concentration floats with radius and time (Fig.1 ):
Fig.1 : Diffusion Model (left) and the comparison of concentration at different distance

Diffusion Model with Production Consumption

We place another downstream enzyme beside the current enzyme, so the production will be consumed. The consumption rate is considered to be Michaelis-Menten equation, and simplified as proportional to the product concentration.
  • ∂u/∂t-a22 u=f(r)
  • u(r,0)=0
  • ∂u/∂r(R,t)=0
While:
f(r)=(vmax [A])/(Km+[A]), r ≤ r0 or
f(r)=0, r > r0 or
f(r)=-βu, r∈Φ
Then we have three-dimensional graphics as following (Fig.2 ):
Fig.2 : Membrane Cluster Diffusion Model near/far
To compare the concentration at different distance, tracing the concentration aside the enzyme we have:
Fig.3 : Comparison of concentration at different distance in Membrane Cluster Diffusion Model

Enzymatic Reaction with Product Inhibition

  • d[P]/dt=k2 [ES]
  • [E][S]/[ES]=k-1/k1 =Ks
  • [E][P]/[EP]=KP
Solving these equations we have:
  • d[P]/dt=(vmax [S])/(Ks (1+[P]/KP)+[S])
  • [P]=(-(Ks+[S] )+((Ks+[S])2+2vmax[S]Ks t/KP )0.5)/(Ks/KP )
Fig.4 : Production in enzymatic reaction with product inhibition

Discussion

As we assumed that the area we chose is a closed region, the products cannot diffuse to the outer side of the area. However when we are using the system in the situation of industry, most of the products can diffuse through the membrane and enter the culture medium. Under this situation the equilibrium will be changed. Reaction beside the membrane barrier can enhance the process the trans-membrane transportation as the concentration of products is much higher than the cytoplasm. Consequently there will be a dramatic increase in product amount.