Team:Potsdam Bioware/Project/Details

Details

Lorem ipsum dolor sit amet, consetetur sadipscing elitr, sed diam nonumy eirmod tempor invidunt ut labore et dolore magna aliquyam erat, sed diam voluptua. At vero eos et accusam et justo duo dolores et ea rebum. Stet clita kasd gubergren, no sea takimata sanctus est Lorem ipsum dolor sit amet. Lorem ipsum dolor sit amet, consetetur sadipscing elitr, sed diam nonumy eirmod tempor invidunt ut labore et dolore magna aliquyam erat, sed diam voluptua. At vero eos et accusam et justo duo dolores et ea rebum. Stet clita kasd gubergren, no sea takimata sanctus est Lorem ipsum dolor sit amet.

Microviridin

Lorem ipsum dolor sit amet, consetetur sadipscing elitr, sed diam nonumy eirmod tempor invidunt ut labore et dolore magna aliquyam erat, sed diam voluptua. At vero eos et accusam et justo duo dolores et ea rebum. Stet clita kasd gubergren, no sea takimata sanctus est Lorem ipsum dolor sit amet. Lorem ipsum dolor sit amet, consetetur sadipscing elitr, sed diam nonumy eirmod tempor invidunt ut labore et dolore magna aliquyam erat, sed diam voluptua. At vero eos et accusam et justo duo dolores et ea rebum. Stet clita kasd gubergren, no sea takimata sanctus est Lorem ipsum dolor sit amet.

Phage Display

In Vivo Selection

Lorem ipsum dolor sit amet, consetetur sadipscing elitr, sed diam nonumy eirmod tempor invidunt ut labore et dolore magna aliquyam erat, sed diam voluptua. At vero eos et accusam et justo duo dolores et ea rebum. Stet clita kasd gubergren, no sea takimata sanctus est Lorem ipsum dolor sit amet. Lorem ipsum dolor sit amet, consetetur sadipscing elitr, sed diam nonumy eirmod tempor invidunt ut labore et dolore magna aliquyam erat, sed diam voluptua. At vero eos et accusam et justo duo dolores et ea rebum. Stet clita kasd gubergren, no sea takimata sanctus est Lorem ipsum dolor sit amet.

Modelling

Motivation

There is no synthetic biology without modeling, of course. In principle there is structure modeling and system modeling. In structure-modeling the conformation and structure of proteins is examined and steric consequences for reactions or the whole system can be estimated. We focused on the second sort of modeling: The systems modeling in which the reaction kinetics of the whole system is analysed and outcomes are predicted. Thus a synthetic biology approach can be chosen because a better understanding of the system is achieved and further changes can be planed - just like in engineering.

Model

The following schema shows the major reactions taking place in our cell system. The Romanic numbers indicate the place of reactions that were written down as equations and then numerically propagated through time. We were able to see that our system works very well in theory. We learned about correct time-scales for our triggering and we were able to identify expected cell-division rates as reference for lab work.

Figure xx: Simplified schema of our cell system including labels and markers for the chemical reactions in it. Marked reactions are considered in our system modeling.

The schema shows the triggered expression of mircoviridin (our inhibitor), the protease (that needs to be inhibited for medical reasons) and the lactamase (that is destroyed by the protease and works as an indicator for our system because it is a protection against ampicillin).
There are three important trigger activation times:
(t0) - beginn (microviridin added already)
t1 - start expression of protease
t2 - start expression of lactamase
t3 - ampicillin added into medium
(t4) - end of the experiment: cellcultures survive or die.
Keeping these times in mind, the reactions can be written down in seven chemical reaction equations of different sort and order.

Figure xxx: Reaction equations between relevant molecules in the microviridin-inhibitor-concept including indications about the triggered time period (t1,t2,t3). The Romanic numbers correspond to the reactions marked in the above schema.

Concentration equations

From the above reaction equations differential equations can be derived that describe the change of substance concentrations. Three concentrations are fixed, however:

Figure xxx: Given substance concentrations for differential equations. Those concentrations are added to the system and will not increase over time.

All other concentrations can be represented in form of differential equations. Between time t1 and t2 four factors are introduced:
k+1 in (1/s*#molecules)- factor for association of mdna (microviridin) and Prot (protease)
k-1 in (1/s) - factor for dissociation of the inhibited protease
kexpr.prot in (#molecules/s) - factor for the expression of protease
kdeg1 in (1/s) -factor for degradation of protease

Figure xxx: Differential equations for substance concentrations from t1 until t2. At t1 the expression of protease begins.

Between time t2 and t3 six additional factors are introduced:
k+2 in (1/s*#molecules)- factor for association of Prot (protease) and TorABla (lactamase)
k-2 in (1/s) - factor for dissociation of protease and substrate
kcat in (1/s) – factor for the catalytic enzyme reaction that cuts lactamase in inactive parts
kexpr.Tor in (#molecules/s) - factor for the expression of lactamase
kdeg2 in (1/s) – factor for degradation of lactamase
ktransTor (1/s) – factor by which lactamase in the cytoplasm is able to pass the membrane and get into the periplasm

Figure xxx: Differential equations for substance concentrations from t1 until t2. At t2 the expression of lactamase begins.

After t3 there are only two more factors that we should introduce:
ktransAmp (1/s) – factor by which the added ampicillin in the medium is able to pass the outer membrane and get into the periplasm
kcat2 in (1/s*#molecules) – factor for the catalytic enzyme reaction that cuts lactamase in inactive parts

Figure xxx: Differential equations for substance concentrations from t1 until t2. At t1 ampicillin is added.

Results

The equations above were solved using MATLAB. Using the Avogadro constant and the volume of a E.coli cell (10^-15L) and a periplam-volume of about 5%, following konstants were calculated and compared with similar literature values:
k+1 = 1.e-6 (1/molecules*s)
k-1 = 2.e-4 (1/s)
kexpr.Prot = 1 (molecules/s)
kdegr1 = 3.e-3 (1/s)
k+2 = 5.e-7 (1/molecules*s)
k-2 = 4.e-4 (1/s)
kcat = 8 (1/s)
kexpr.Tor = 1.4 (1/s)
kdegr2 = 3.e-3 (1/s)
ktransTor = 3.e-3 (1/s)
ktransAmp = 3.e-3 (1/s)
Kcat2 = 0.3 (1/molecules*s)
MIC(0) = 8000 (molecules)
AMP(0) = 5000 (molecules)

Figure xxx: MATLAB numerical solution to our system of differential equations for substance concentrations (above) over all time segments.

It can be seen that there is only a negligible amount of ampicillin in the periplam if the expression works as indicated. There is also a wide tolerance left for suboptimal conditions. The triggertime t1, t2 and t3 shal be about one hour after one another.
In the first section the effect of the inhibition can clearly be seen: Even though the continuous expression of protease, only a very small number remains active in the cytoplasm. In the second section lactamase is expressed. Most of the lactamase in the cytoplasm is destroyed by the protease, however the number is by far large enough to steadily release molecules into the periplam where it cannot be affected by the protease any more. Because the volume of the periplasm is very small, the concentration there is even higher than the absolute amount of molecules in the graph suggests compared to the amount of molecules in the cytoplam. The third section shows that the concentrations remain very steady and the cells are ampicillin resistent.
Cells passing this procedure are very ampicillin resistent and grow fast: They double their cell volume about every 20 minutes. If an error would appear and the ampicillin concentration inside the periplasma would increase over 2µg/ml, the growth rate would slow down drastically and the cells might die. This way defect cell cultures can easyly be seperated from the good ones.
The ratio of ampicillin concentration in n µg/ml to the growth rate (cells double during this timespan) is approximatelly this: T(groth) = 20min+10*2^(n-1).

MATLAB code

Media:conzt1.m
Media:conzt2.m
Media:conzt3.m
Media:plottfuntionendreiaufeinmal2.m