Team:Stanford/Research/Modeling

From 2010.igem.org

Contents

Goals

Our intuition for what makes a good ratio sensor could only take us so far. From the very first stages of design, we wanted to back up and test our ideas with mathematical tools. Luckily, we found that solving the equations of mass action kinetics at steady-state was enough to give us clear design criteria. We present the mathematical basis for sensors that are capable of sensing a single ratio digitally, or many ratios in an analog fashion.

Kinase/Phosphatase System

Throughout this derivation, we refer to Goldbeter and Koshland's seminal work on the kinetics of emzyme pairs [1]. For an opposing pair of enzymes X and Y modifying and unmodifying a substrate Z with modification p,

MainGBK.JPG

From Michaelis-Menten kinetics we know that the rate at which Zp is dephosphorylated is

R1.JPG [2]

and the rate at which Z is phosphorylated is

R2.JPG [2]

We also know that Z can be modified an unmodified, but its total is conserved: [Z]o = [Zp] + [Z]

At steady state, modification and unmodification rates are equal, leading us to

1stblock.JPG [2]

where

2nd block.JPG [2]


Rearranging gives

3rdblock.JPG

Now, assume that the two opposing enzymes are nowhere near being saturated with substrate:

4thblock.JPG

which is the equation for a Michaelis-Menten-shaped curve as shown below:

Langmuir.JPG

So, the fraction of modified substrate is a saturating function of the ratio of inputs, the half-max value is the ratio of the strengths of the two enzymes, and the fraction of unmodified substrate is a linear function of the ratio of inputs,

6thblock.JPG

We now know that the requirements for this analog ratio sensor are

  1. Consistent relationships between inputs and enzyme activities (non-trivial)
  2. A non-saturating amount of substrate

But what happens if we disobey requirement 1? Let's increase the amount of substrate to saturation so that

7thblock.JPG
8thblock.JPG

When the enzymes are saturated, the only steady-state solution is v1 = v2. Any ratio of v1/v2 below this, and all the substrate becomes modified. Any ratio above and all the substrate becomes unmodified. Simply by changing the concentration of substrate, we have converted our analog ratio sensor to a digital ratio sensor.

The requirements for this digital ratio sensor are:

  1. Consistent relationships between inputs and enzyme activities (non-trivial)
  2. An amount of substrate sufficient to saturate both opposing enzymes

sRNA System

This sRNA system is bases its analysis on the work of Nicolas Buchler. [3]

Consider an sRNA that can bind to an mRNA and inhibit its translation. The pair binds and dissociates with dissociation constant Kd.

For a moment, imagine that this interaction is infinitely strong, such that Kd -> 0. In this case, all the mRNA available will bind to the sRNA, until every sRNA is occupied by one mRNA. Once this threshold, defined by the total sRNA concentration, is reached, additional mRNA added is free to be translated. Here we show the accumulation of translatable free mRNA, and hence Gene Expression vs. total mRNA, in the presence of varying concentrations of infinitely strong sRNA

SRNA1.JPG

At the point that the sRNA becomes saturated, two interesting things occur.

  1. There is a sharp transition between no gene expression and a linearly-increasing amount of gene expression.
  2. This transition occurs only when the concentration of mRNA is equal to the concentration of sRNA.

We have designed a sensor that responds to a 1:1 mRNA:sRNA ratio with a sharp change in output.

But will this system have the same desirable properties when the sRNA is only realistically strong? Here we assume a 1 nM interaction between sRNA and mRNA

SRNA2.JPG

We find that the transition is still quite sharp when the sRNA has a realistic dissociation constant.

References

[1] [http://www.ncbi.nlm.nih.gov/pubmed/6947258 Goldbeter A, Koshland DE Jr.: An amplified sensitivity arising from covalent modification in biological systems. PNAS. 1981 Nov; 78(11) 6840-4.]

[2] [http://en.wikipedia.org/wiki/Goldbeter-Koshland_kinetics Goldbeter-Koshland kinetics In Wikipedia. Retrieved July 13, 2010, from http://en.wikipedia.org/wiki/Goldbeter-Koshland_kinetics]

[3] [http://www.ncbi.nlm.nih.gov/pubmed/18938177 Buchler NE, Luis M. Molecular titration and ultrasensitivity in regulatory networks. J. Mol. Biol. 2008 Dec 31; 384(5):1106-19]