Team:Imperial College London/Modelling/Protein Display/Detailed Description

Objectives; Detailed Description; Parameters & Constants; Results & Conclusion;Download MatLab Files;

1. Identification of all active and relevant elements of the isolated part of the system.
2. Identification of interactions between the identified elements.
3. Identification of threshold concentration of auto inducing peptide (AIP) needed for activation of the receptor.
4. Defining a control volume around the bacterial cell (after cleavage, the surface protein will float around in the extracellular environment).
5. Determining the surface protein production to estimate the maximum abundance on the bacterial surface.

1. Elements of the system

1. The surface protein consists of a cell wall binding domain, linker, AIP (Auto Inducing Peptide). It is expressed by constitutive gene expression. It was assumed that the bacteria would be fully grown before carrying out the detection so the cell wall would be covered with as many surface proteins as the cell can maintain.
2. Schistosoma elastase (enzyme released by the parasite) cleaves the AIP from the cell wall binding domain at the linker site. In the laboratory we used TEV protease as we could not get hold of Schistosoma elastase, so the model was adjusted appropriately (TEV enzyme kinetic parameters were used).
3. The ComD receptor is activated by a sufficiently high AIP concentration. Once activated, ComD signals into the cell to activate the colour response.

2. Interactions between elements

• E+S ↔ ES → E+P
• Substrate (S) = Protein
• Enzyme (E) = TEV (Protease)
• Product (P) = Peptide

3. Threshold concentration of AIP

4. Control volume selection

Initial Choice of Control Volume

Control volume initial choice
This control volume is considered to be wrong by us, but the details were kept for the reference.
The control volume: The inner boundary is determined by the bacterial cell (proteins after being displayed and cleaved cannot diffuse back into bacterium). The outer boundary is more time scale dependent. We have assumed that after mass cleavage of the display-proteins by TEV, many of these AIPs will bind to the receptors quite quickly (eg. 8 seconds). Our volume is determined by the distance that AIPs could travel outwards by diffusion within that short time. In this way, we are sure that the concentration of AIPs outside our control volume after a given time is approximately 0. This approach is not very accurate and can lead us to false negative conclusions (as in reality there will be a concentration gradient, with highest conentration on the cell wall).

Control volume (Volume of B.sub to be excluded. x indicates the distance travelled by AIPs from the surface in time t).

Difussion distance was calculated using Fick's 1st Law: x=√2Dt, where: x - diffusion distance, D - diffusion constant, t - time of diffusion
D_average = 10.7x10^-11 m²s^-1 for a protein in agarose gel for pH=5.6 [3]
t = 8s (arbitrarily chosen)
This gives: x = 4.14x10^-5m
The control volume can be calculated by adding 2x to the length and the diamter of the cell.
This gives a control volume (CV) = 4.81x10^-7m³

We realized that our initial choice of control volume was not accurate because we had assumed that the bacteria was the medium. However, in reality bacteria live in colonies very close to each other. Since our bacteria was meant to be used in suspension we had to reconsider this issue.

CFU/ml vs. Optical Density at 600nm (OD600).

The graph shows values of CFU/ml for different optical densities. The range of CFU/ml is therefore between 0.5x10^⁸ - 5x10⁸.
In this calculation, we assumed that only one cell will grow and become one colony (i.e. no more than one cell will form no more than one colony). Therefore, the maximum number of cells in 1ml of solution is 5x10⁸. Taking the volume of 1 ml = 10^-3 dm³ and dividing by the (maximum) number of cells in 1ml gives the average control volume (CV) around each cell: 2x10^-12 dm³/cell. For simplicity, we choose the control volume to be cubic. Taking the third root of this value gives the length of each side of the control volume.

Side length of CV = y = 1.26x10^-4 dm = 1.26x10^-5 m.

• Firstly, assume that the cells will be placed in the centre of the CV. Hence, after cleavage the protein will have an average distance of y/2 to travel in order to cross the boundary of the CV. This is calculated to happen within 0.18s. Even if the bacterium was not placed in the centre of the CV, the protein will travel from one end of the cube to the other in less than one second (~0.74s). Hence, it will take between 0.18 and 0.74s for the concentration of AIPs around the cell to be uniform. Noticing that these time values are very small, we can approximate our model to have a uniform concentration across the volume. Since we are underestimating the value of AIP concentration right next to the cell's surface, we are overestimating the time required for the AIP concentration to reach the threshold level.
• We can neglect the diffusive fluxes across the CV border (see figure below). Assuming that adjacent cells are producing the peptide at the same rate and that the concentration of TEV is the same around the cell, then the fluxes should be of the same value giving a net flux of zero. Hence, we can neglect diffusion and have our model limited to one bacterium.

Figure showing two cells with their control volumes.

Limitations resulting from our assumptions
Most of our assumptions concerning the control volume were plausible due to the careful choice of cell density = 5x10⁸ CFU/ml.
If the density is varied significantly, then our simplifications might not hold any longer. However, this does not mean that our system cannot function for lower cell densities. Our model might not be very accurate for predicting situations with cell densities that are much higher or lower than 5x10⁸ CFU/ml.
It was decided that the model should not be used for cell densities lower 10⁷ CFU/ml. Below this value AIPs take more than 1 second to diffuse accross half of a side-length of the control volume (assuming that the cell is inside control volume). We agreed that below 10^⁷ CFU/ml the approximation about uniform concentration throughout the control volume could be wrong and that the concentration gradients could become more significant. If our model was applied to this particular situation, it would possibly overestimate the time taken to activate the receptor.
It is not possible to increase the cell density by more than 10⁹ CFU/ml because of the size of a cell.

5. Protein production

• The paper mentions that each cell displays 2.4x10⁵ peptides. [2]
• 2.4x10⁵ molecules = 2.4x10⁵/6.02x10²³ mol = 0.398671x10^-18 mol
• Volume of B.sub: 2.79x10^-15 dm³
• Concentration = [mol/L]
• c = 1.4289x10^-4 mol/dm³. This is the concentration of protein that will be displayed on a single cell of B.sub.

Production of protein by transcription and translation.

The degradation rate was kept constant, and the production rate was changed according to the final concentration.

References

1. Havarstein, L., Coomaraswamy, G. & Morrison, D. (1995) An unmodified heptadecapeptide pheromone induces competence for genetic transformation in Streptococcus pneumoniae. Proc. Natl. [Online] 92, 11140-11144. Available from: http://ukpmc.ac.uk/backend/ptpmcrender.cgi?accid=PMC40587&blobtype=pdf [Accessed 27th August 2010]
2. Kobayashi, G. et al (2000) Accumulation of an artificial cell wall-binding lipase by Bacillus subtilis wprA and/or sigD mutants. FEMS Microbiology Letters. [Online] 188(2000), 165-169. Available from: http://onlinelibrary.wiley.com/doi/10.1111/j.1574-6968.2000.tb09188.x/pdf [Accessed 27th August 2010]
3. Gutenwik, J., Nilsson, B. & Axelsson, A. (2003) Determination of protein diffusion coefficients in agarose gel with a diffusion cell. Biochemical Engineering Journal. [Online] 19(2004), 1-7. Available from: http://www.sciencedirect.com/science?_ob=MImg&_imagekey=B6V5N-4B3MXDC-2-K&_cdi=5791&_user=217827&_pii=S1369703X03002377&_origin=search&_coverDate=07%2F01%2F2004&_sk=999809998&view=c&wchp=dGLzVtb-zSkzS&md5=c17d0e7320f03931006f9b1a10a438b9&ie=/sdarticle.pdf [Accessed August 20th 2010]
4. Imperial College London (2008) Biofabricator Subtilis - Designer Genes. [Online] Available from: https://2008.igem.org/Imperial_College/18_September_2008 [Accessed 1st September 2010]