Team:Hong Kong-CUHK/modelling-1.html

From 2014.igem.org

Contents

Model for protein expression system

1. Overview

Amino acids are the basic building blocks of protein. The rate of protein translation (V) increases with the concentration of amino acids (C) and bounds when the amount of amino acids is saturated.(i.e. there exists Vmax, such that V is smaller or equal to Vmax no matter how high C is.) Normally, the protein expression rates of different system can be compared by the transcription efficiency of the system only because we provide enough nitrogen sources to make V close to Vmax. However, this is not quite applicable to our protein expression system. The transcription process of our system is repressed by nitrogen sources in the cell, so we believe that it is not practical to make C reaches saturated. However, this system allows us to achieve a constant protein expression rate by controlling the duration for expressing T7 polymerase.

2. Analysis and Discussion

Let us review our protein expression system. The schematic diagram is shown below. <img src="">

2.1 Assumptions

First, we may assume that:
a. Concentration of σ-RNA polymerase is constant.
b. All proteins being expressed share the same nitrogen source.

Then, let:
[N] be the concentration of nitrogen source in cell;
[T7p] be the concentration of T7 polymerase;
[TP] be the concentration of target protein;
Knc be the rate of nitrogen consumption out of this protein expression system;
Vmax([T7p]) be the maximum rate for target protein expression;
A,B,C,D,E,F,G,H, be constant.

2.2 Discussion of the two stages

In the first stage:

In this stage, we want to optimize the production of T7 polymerase by setting a relatively low concentration of nitrogen source.
We have the following equation for the rate of change of [N]: <img src="">

Where in the right-hand side of the equation,
First term corresponds to the nitrogen consumption rate for T7 polymerase;
Second term corresponds to the nitrogen consumption rate for target protein;
Third term corresponds to the nitrogen consumption rate of the bacteria without the system;
A is the nitrogen needed to expression one T7 polymerase;
B is the nitrogen needed to expression one target protein;
C is the maximum speed of nitrogen consumption of unmodified bacteria;
D is the concentration needed for the rate of nitrogen consumption equals to half of maximum speed.

Now, we need to figure out the formulae of rate of change of [T7P] and [TP].
Since [TP] is not the major focus in this stage, let us ignore this temporarily.
Then, we need to predict the formula for rate of change of [T7P]. Below is the information we have:
i As [N] tends to 0, rate of change of [T7P] tends to 0.
ii When [N] is very large, rate of change of [T7P] tends to 0.
iii When [N] is saturated for the translation process (denoted by [Nsat]), the second derivative of [T7P] is negative.
vi The formula is continuous.
Base on the above information, we proposed the following two formulae:

<img src="">22:58, 17 October 2014 (CDT)~ (2)
<img src="">(3)

Reason for proposing (2): we assume that the repressive power of nitrogen sources grows exponentially.
<img src="">

Graph of (2) with different values of F. (red: 0.5; blue: 0.75; green: 1) Reason for proposing (3): we assume that the repressive power of nitrogen sources grows linearly.
<img src="">

Graph of (3) with different values of G. ([Nsat] = 1) (red: 0.5; blue: 1; green: 2)

Now, we want to find out the optimized value for [N] to produce polymerase as fast as possible. Therefore, we want to have the maximum of (2) and (3).
For (2), we have
<img src="">

Solving the equation, we have [N optimal] = 1/F,
For (3), we have [N optimal] = [Nsat]
If we want to set [N] = [N optimal] for a certain period, we need to find the rate of change of [N] at [N] = [N optimal]
When [N] is optimal, [N] is quite small. As we assume that the transcription efficiency of nifH promoter is very high with small [N] and [T7p] initially is very small, we have

 <img src="">

In the second stage:
In this stage, we want to optimize the production of target protein by setting an excess nitrogen source.
We have the equation for the rate of target protein production:

  <img src="">22:58, 17 October 2014 (CDT)~~(4)

The formula is similar to the Michaelis-Menten Equation with Vmax replaced by a function of [T7P].
When [N] is very large, we have rate of change of [N] closed to zero.
Also, rate of change of [T7P] closes to zero as the nifH promoter is highly repressed.
This implies that [T7P] is a constant
As a result, Vmax([T7P]) is a constant
From (4), when [N] is very large, rate of change of [TP] closes to Vmax([T7P]).
Hence, we have rate of target protein production is a constant.

2.3 Constant Search

The composition of T7 RNA polymerase has been investigated [1]. From the table above, the expected value of A can be found.
If the target protein chosen have been investigated, B should be able to be found in some literatures
C, D can be found experimentally:
i. We can grow wild type Azotobacter vinelandii and monitor the nitrogen concentration at different time points.
E, F, G can be found experimentally:
i. We can calculate [T7P] by measuring the RNA concentration. Here, we assume that RNA concentration is always directly proportional to the [T7P].
ii. We can modify T7 polymerase and find [T7P] with the aid of western blot.
It is not necessary to search H. We can predict without knowing H.

Reference:

[1] Moffatt, Barbara A., John J. Dunn, and F. William Studier. "Nucleotide sequence of the gene for bacteriophage T7 RNA polymerase." Journal of molecular biology 173.2 (1984): 265-269.