Team:Edinburgh/Phage Replication

Phage Replication

A basic activity in biorefinery consists of the degradation of cellulose, due to the presence of enzymes. We are not only concerned with the activities and the amount of enzymes, but also with metabolism and activities of bacteriophage.

M13 Replication

M13 is a filamentous bacteriophage: a worm-like virus approximately 1 um long with a 10 nm diameter that infects only E.coli. The viral particle consists of a single-stranded, closed circular DNA core surround by a protein coat. Prior to virus assembly, the coat proteins are fixed in the bacterial membrane by transmembrane domains. During assembly, viral DNA is extruded through the membrane and concomitantly enveloped by coat proteins. The ends of the assembled virus are capped by four minor coat proteins, and the length of the filament is covered by several thousand copies of the major coat protein(P8).

From Slonczewski and Foster (2010).

• The M13 phage attacks E. coli (host), multiplies in the host cell cytoplasm, and is released without causing the bacteria’s death (non-lytic).

Phage dynamic model

• dx/dt=ax-bvx

(Rate of change of quantity of uninfected E. coli equals to the uninfected E. coli replicate itself minus the E. coli infected by M13 phage.）

• dy/dt=ay+bvx

（Rate of change of quantity of infected E. coli equals to the quantity of infected E. coli replicate itself plus the E. coli infected by M13 phage.）

• dv/dt=cy-bvx-mv

（Rate of change of quantity of free phage equals to the phage released by infected E. coli minus the phage which is to infect an E. coli and the decayed phage.）

X(t) — uninfected E. coli

Y(t) — infected E. coli

V(t) — free phage

a — replication coefficient of E. coli

b — transmission coefficient of phage

c — replication coefficient of phage

m — decay rate of phage

Simulations

The MATLAB code uses a Runge-Kutta method of order four to solve the system.

simulation value: x0=y0=v0=2.00E5

The above figure shows a simulation going over 15 hours.

The simulation shows the infected E. coli population dominates. And the phage population decreases at first then increases.

Synergy on each phage

The synergy function, which means having enzymes closer together, is supposed to increase the efficiency of cellulose degradation. We attempted to construct a model of cellulase with regard of synergy.

For the model we construct, we assume that

1. Endoglucanase cuts cellulose chains in the middle,exoglucanase chews away at the end of a cellulase chain, they work together to produce cellobiose.

2. Because of the proximity of the β-glucosidase and the other two enzymes, as well as the sufficient amount of β-glucosidase on each phage, we assume the mdiate production, cellobiose, is all converted to glucose.

3. Production, which refers to glucose here, inhibits the action of the above enzymes.

Cellulase model

Substrate reactivity

• Rs = S/S0

Rate of reaction

• dG/dt = 1.056*2*(k1r*E1B*S*Rs) /(1+(G/K1IG))
• dS/dt = -(k1r*E1B*S*Rs) /(1+(G/K1IG))

Simulations

in biorefinery

The idea here is to bring the two models above together, which gives us a panaramic view of cellulose degradation using displayed phage, as well as make it possible to calculate the cost and profit in a biorefinery.

References

• Slonczewski JL, Foster JW (2010) Microbiology: An Evolving Science, 2nd edition. W. W. Norton & Company
• Robert J.H Payne, Vincent A. A. Jansen(2011)Understanding Bacteriaphage therapy as a density-dependent kinetic process
• Cattoen C (2003) Bacteriaphage mathematical model applied to the cheese industry
• A.O. Converse, J.D. Optekar(1993) A synergistic Kinetics Model for Enzymatic Cellulose Hydrolysis Compared to degree-of-synergism Experimental Results
• Kadam KL, Rydholm EC, McMillan JD (2004) Development and Validation of a Kinetic Model for Enzymatic Saccharification of Lignocellulosic Biomass