Team:ULB-Brussels/Modelling

From 2014.igem.org

(Difference between revisions)
m
m
Line 229: Line 229:
</section></section>
</section></section>
 +
<section style="margin: -5px"></section>
</tr>
</tr>

Revision as of 19:32, 14 October 2014

$~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ \newcommand{\MyColi}{{\small Mighty\hspace{0.12cm}Coli}} \newcommand{\Stabi}{\small Stabi}$ $\newcommand{\EColi}{\small E.coli} \newcommand{\SCere}{\small S.cerevisae}\\[0cm] ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ \newcommand{\PI}{\small PI}$ $\newcommand{\Igo}{\Large\mathcal{I}} \newcommand{\Tgo}{\Large\mathcal{T}} \newcommand{\Ogo}{\Large\mathcal{O}} ~$ Example of a hierarchical menu in CSS

carousel slider




- Université Libre de Bruxelles -


Pop Dyn > TA Sys > 2A Pep >

$Overview$


$\Igo$n this section, we will introduce the modelling part with an overview of our biological system and its expected effects with numerical simulations of bacterial populations. Structural components of the chosen plasmids will also be covered. In next sections, the mathematical procedure will be detailed.

Afterwards, these simulations will be compared with experimental results. In parallel, an estimation of the production of sub-populations in bioreactors will be made by coupling the recombinant protein with an essential protein and a numerical estimation will be runned. Now, we're continuing to perform it.

The following diagram was drawn to identify the principal components of our biological system (i.e. the effects related to the sequences induced by a promotor, to the growth of bacteria, to the transcription and translation of a specific protein) :

Figure m1 : This diagram illustrates the production of GFP fluorescent proteins controlled by PSK (post-segregational killing), in the case of TA CcdBA system.


Population Dynamics Model

~Pop~

The growth of bacteria involves ...

a Population Dynamics Model can be fitted in our system. Theorically, two approaches have been planned:

$1.1)$ $By$ $Probabilities$

When some new plasmids are introduced into the cytoplasm of E.Coli bacteria, it does not garantee that the daughter cells will contain it. Indeed, these plasmids can be lost after cell division or replication. In this regard, it is interesting to study a model based on the different possibilities of plasmid combinations in bacteria, as in the study of mutations in animals. A typical example of a similar case is found if we study the mutations of the eyes color in a family, by vertical genes transfer. Of course, in the case of plasmids, one must also take horizontal genes transfer into account.

A Probabilistic Model is useful because easily undertstandable, but requires some assumptions. Actually, models based on genetical alphabets are continuing to be suggested in the literature [19]. A Monte Carlo rejection sampling method for DNA using conditional probabilities was recently designed for a 4-letter alphabet [20]. The irreversibility and spatial asymmetry implic that the DNA molecule can be viewed as an out-of-equilibrium structure.

$1.2)$ $By$ $Logistic$ $\small\&\normalsize$ $Lotka$ $Equations$

The Logistic Equation was initially introduced during the beginning of the XIXth Century, by the belgian mathematician P.F. Verhulst. Now, this equation is mainly used in Population Dynamics Models, especially in Biological Sciences. Mathematicians currently finish Ph.D thesis using it, and the analytical Lotka-Volterra model is directly associated with the Verhulst theory. Another interesting model is obtained from Euler-Lotka equation to the Leslie matrix coefficients.

Other models exist, f.e. by Monod equation, but these ideas are less consistent with our global and partial systems.

-

Toxin-Antitoxin Systems

~TA~

Two type II TA systems are investigating in our $\MyColi$ project.

The $\hspace{0.12cm}\small\mathtt{1}\normalsize^{st}$ consists of ccdB (the toxin, T) and ccdA (the antitoxin, A) and for the $\hspace{0.12cm}\small\mathtt{2}\normalsize^{nd}$ these are Kid (T) and Kis (A):

$\newcommand{\AA}{\mathbb{A}} \newcommand{\CC}{\mathbb{C}} \newcommand{\TT}{\mathbb{T}} \newcommand{\GG}{\mathbb{G}} \newcommand{\KK}{\mathcal{K}}$

$2.1)$ $CcdBA$

One of the most studied and characterized TA systems, CcdBA involves two principal components : ccdB and its antidote ccdA. As we explained in the introduction page of our project, ccdB is an inhibitor of the DNA gyrase, so it binds the subunit A of the DNA gyrase complex when it's bound to DNA. When DNA double strand is broken, there is activation of emergency signals (SOS system blocks cellular division in bacteria). If the DNA gyrase is not protected by a mutation (such events are possible, but excessively rare) or if the antidote is degraded (very frequent because ccdA is unstable in comparison with ccdA), the death of a bacterium in unavoidable [13-18].

This TA system will be implemented for $\EColi$ bacteria.

$2.2)$ $Kis/Kid$

The same is relevant about the second TA system studied. In this case, the two principal components are Kid and its antidote Kis and the parameters are chosen a little bit different than in the first system.

This TA system will be implemented for $\SCere$ yeasts. By modelling and by comparison with experiments, we would obtain finally a model that correctly describes our TA system.

-

Bioreactor and 2A Peptide

~2A~

The p2A peptide is an important building block in our system.

We think that its use can enhance the production of bioreactors .

In this perspective, we will compare the classical $\Stabi$ system with our $\MyColi$ system and identify when Mighty Coli should be applied in a bioreactor. We will also investigate the effective gain in production of a protein with the Mighty Coli system.

We would like to mention that different 2A peptids exists (f.e. F2A and P2A), but these peptids are not completely characterized. A future continuation of our project could be a study of this aspect to refine our current model and to build new Mighty Coli Biobricks that use different 2A peptids [14,16].

$\hspace{0.12cm}$Bibliography

  • [13] A.V. Hill, (1910). The possible effects of the Aggregation of the molecules of hemoglobin on its Dissociation curves, J. Physiol, No.40, iv-vii.
  • [14] T. Ogura, S. Hiraga, (1983). Mini-F plasmids genes that couple Host cell division to Plasmid proliferation, Proc. Natl. Acad. Sci. USA, 80, 4784-4788.
  • [15] M. Santillán, (2008). On the use of the Hill functions in Mathematical models od Gene regulatory networks, Math. Model. Nat. Phenom., Vol.3, No.2, 85-97.
  • [16] P. Wang, R.E. Dalbey, (2010). In Vitro and in Vivo approaches to studying the Bacterial signal Peptide processing, Springer Protocols, A. Economou ed. Humana Press, Protein Secretion, Methods in Molecular Biology 619, 21-37.
  • [17] L. Gelens, L. Hill, A. Vandervelde, J. Danckaert, R. Loris, (2013). A general model for Toxin-antitoxin module dynamics can explain Persister cell formation in E.coli, PLOS Computational Biology, Vol.9, Iss.8, e1003190.
  • [18] N. Goeders, L. Van Melderen, (2014). Toxin-antitoxin systems as Multilevel interaction systems, Toxins, 6, 304-324, ISSN 2072-6651.
  • [19] D.A. Malyshev, K. Dhami, T. Lavergne, T. Chen, N. Dai, J.M. Foster, I.R. Corrêa Jr & F.E. Romesberg, (2014). A semi-synthetic organism with an expanded genetic alphabet, Research Letter, Nature 13314, 1-17.
  • [20] A. Provata, C. Nicolis & G. Nicolis, (2014). DNA viewed as an out-of-equilibrium structure, Phys. Rev. E 89, 052105.

Population Dynamics >