Team:McGill/Modeling

From 2009.igem.org

(Difference between revisions)
(Introduction)
(Introduction)
Line 65: Line 65:
where Ψ1 and Ψ2 represent the concentrations of the activating and inhibiting molecules, respectively, γi the degradation constant, Di the diffusion constant, λi the maximal synthesis rate of molecule i, and δ the Dirac function. fi represents the Hill function describing the dependence on the opposing molecule:
where Ψ1 and Ψ2 represent the concentrations of the activating and inhibiting molecules, respectively, γi the degradation constant, Di the diffusion constant, λi the maximal synthesis rate of molecule i, and δ the Dirac function. fi represents the Hill function describing the dependence on the opposing molecule:
[[Image: Mcgill09HillFunctions.png|frame|center]]
[[Image: Mcgill09HillFunctions.png|frame|center]]
 +
where n, b, and θ are positive.
 +
To simplify the analysis, we chose equal parameters between the activating and inhibiting sites (Appendix A).
 +
 +
=='''Numerical Simulation'''==
 +
The above system was solved numerically using a forward Euler scheme in time and a centered difference scheme in space. Cyclical boundary conditions were assumed; meaning the spatial dimension formed a ring. This was chosen since simulating an approximate infinite line is computationally costly. However, this also allowed us to investigate two site geometries during one simulation (explained later). The ring was given a physical length of 50 and discretized into 500 intervals. For simplicity, separation distances will be reported in terms of numerical intervals rather than physical distance.
 +
 +
 +
=='''Appendix A - Parameters'''==
 +
The standard set of parameters used to observe oscillations were taken directly from Shymko and Glass (1974).
 +
γ = 2
 +
D = 2
 +
λ = 54
 +
θ = 1
 +
b = 0
 +
N = 8
 +
Both strains were assumed to have identical parameters in order to simplify the model as well as explicitly observe the dependence of dynamics on separation distance.

Revision as of 05:57, 21 October 2009


Home Home Team Team Project Notebook Results Sponsors

Introduction

Many models examining intercellular signaling do not take into account the separation distances of the signaling bodies. We use a partial differential equation (PDE) based model to gain insight into spatially heterogeneous activation-inhibition intercellular signaling.

Two types of signaling molecules exist: activating and inhibiting. Each molecule is synthesized by a unique strain of cells and affects the synthesis rate of the other strain.

Figure 1 – Activation-inhibition intercellular signaling – Activating molecule (A) synthesized and diffuses to increase synthesis of inhibiting molecule (B) in secondary strain. Inhibiting molecule also diffuses back to initial cell and decreases synthesis of activating molecule.

This is modeled using the following system of PDEs:

Mcgill09PDEs.png

where Ψ1 and Ψ2 represent the concentrations of the activating and inhibiting molecules, respectively, γi the degradation constant, Di the diffusion constant, λi the maximal synthesis rate of molecule i, and δ the Dirac function. fi represents the Hill function describing the dependence on the opposing molecule:

Mcgill09HillFunctions.png

where n, b, and θ are positive. To simplify the analysis, we chose equal parameters between the activating and inhibiting sites (Appendix A).

Numerical Simulation

The above system was solved numerically using a forward Euler scheme in time and a centered difference scheme in space. Cyclical boundary conditions were assumed; meaning the spatial dimension formed a ring. This was chosen since simulating an approximate infinite line is computationally costly. However, this also allowed us to investigate two site geometries during one simulation (explained later). The ring was given a physical length of 50 and discretized into 500 intervals. For simplicity, separation distances will be reported in terms of numerical intervals rather than physical distance.


Appendix A - Parameters

The standard set of parameters used to observe oscillations were taken directly from Shymko and Glass (1974). γ = 2 D = 2 λ = 54 θ = 1 b = 0 N = 8 Both strains were assumed to have identical parameters in order to simplify the model as well as explicitly observe the dependence of dynamics on separation distance.