Team:McGill/Modeling

From 2009.igem.org

Revision as of 14:58, 21 October 2009 by Fnaqib (Talk | contribs)


Home Home Team Team Project Notebook Results Sponsors

Contents

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.

One Oscillator

We first explored the different potential dynamics when the separation distance between an activating and inhibitory site was increased.

Figure 2 – One Oscillator – The red bar represents the inhibitory site, which remains fixed in position while the activating site, blue bar, is sequentially moved around the ring.

The following is an example of the dynamics observed when the two sites are at a distance where oscillations occur.

Figure 3 – One Oscillator – One activation and one inhibitory site separated by 11 intervals. The concentration of the activating molecule is measured at the activating site and the inhibitory molecule at the inhibitory site.

By analyzing curves similar to those in figure 3 the frequency of oscillations was calculated. Figure 4 illustrates the sudden appearance of high frequency oscillations when the sites are near each other and the decreasing frequency as they are moved apart until oscillations disappear.

Figure 4 - Frequency of Activation Molecule at Activation Site vs Separation Distance – The symmetrical appearance is a result of the activation site at first moving away from the inhibitory site, making its way around the ring, and then getting closer to the inhibitory site along the opposite side. The graph of frequency of the inhibitory molecule at the inhibitory site vs separation distance was equivalent.

The observation that the steady state becomes unstable for a range of separations has been reported previously by Shymko and Glass (1974).

Two Oscillators

We next looked at a system consisting of two oscillators, where each consists of an activation and inhibitory site.

Figure 5 – Two Oscillators – Each oscillator consists of one activation and one inhibitory site. There are two arrangements of the system: BA AB, termed the symmetrical system, and BA BA, the unsymmetrical system. In this document we only discuss the symmetrical system. Oscillator 2 travels around the ring.

The distance between the two oscillators was varied while the distance between two sites within an oscillator was held fixed at 5 intervals. This value was chosen for demonstration purposes, however the dynamics to be described have been observed at various separation distances.

We first look at the change in frequency as the oscillators are moved apart. [[Image: McGill09Frequency_vs_Distance_for_Osc2_Geometry1.png|frame|center|Figure 6 - Frequency of Activation Molecule at Activation Site of Oscillator 2 vs Separation Distance – There are four separation distances that resulted in a frequency of -1. This value is used as a marker to denote a curve who's frequency did not converge to a single value. As will be seen later, these distances resulted in complex dynamics.

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.