Team:DTU-Denmark/Modeling

From 2011.igem.org

(Difference between revisions)
(General kinetic model)
Line 20: Line 20:
\begin{eqnarray}
\begin{eqnarray}
-
\leftrightharpoons
+
\leftrightharpoons_k^l
\label{eq:r1}
\label{eq:r1}
&\ce{\textit{\color{blue}m} + \textit{\color{red}s} <=>[k_{1,s}][k_{-1,s}] $c_{ms}$ ->[k_{2,s}]} (1 - p_s) {\color{red}s}
&\ce{\textit{\color{blue}m} + \textit{\color{red}s} <=>[k_{1,s}][k_{-1,s}] $c_{ms}$ ->[k_{2,s}]} (1 - p_s) {\color{red}s}

Revision as of 11:44, 12 September 2011

Project

On the design of life

The way humans do design is fundamentally different than the evolutionary principles from which life origins. Life has evolved due to an incredible long series of trial-and-error driven and governed by selectional pressure. As this approach have led to an enormous variety of solutions to ways of being alive on earth, from the most sophisticated machine known to mankind, the human brain, to the most resilient microorganisms, there is no doubt to the success of evolution. But what principles may humans apply if we where to do biological design with similar success? One possibility is to guide and facilitate selection, a principle adhered to by breeding programs of the last millennia. This way we might change wolves into dogs or improve crop yields, but there might be and probably are impossibilities of this approach. What if something really useful is only achievable by completely redesigning whole organisms or specific pathways as opposed to simply gradually changing them? A wonderful example of redesign of organisms is the production of human insulin from microorganisms, an application made possible by novel understanding of the components of life and recombinant DNA technology. But from an engineering perspective introducing one recombinant gene is a lot simpler than introducing whole pathways not to speak of designing pathways a priori. The reason microbial insulin production was possible is that the problem reduces to one well characterized component namely the gene encoding insulin. This way the topic of complexity is avoided by utilizing the abilities of the cells to perform the complex task of protein synthesis. Thus one of the most successful cases of applied biotechnology and the redesign of life highlights our current limited ability to understand and control life itself.

To broaden the scope of synthetic biology and biology in general it is necessary to address complexity and not simply reduce scientific inquires to individual components. The emergent properties of life which we ultimately want to understand and control arises from interactions between components. These interactions are often so complex that the unaided human brain might be unable to perceive what is going on. For example the large subunit of the ribosome in simple prokaryotes consists of 2 rRNAs and 34 different proteins all exerting some function. Even if one had the precise structure and function of each component it is unlikely that we would have the ability to fully grasp and understand this subunit of the ribosome. Thus there are clearly limits to human understanding of which one must be aware to successfully conduct research. It is clear that some knowledge of components are required but one might continue to gather knowledge of components for all eternity without gaining understanding of the phenomena of protein synthesis. This indicates that not all knowledge is equal and some knowledge might even be useless to practical applications in the life sciences. Today biology is generating discoveries and insights at an unprecedented rate but despite of the massive gain in knowledge promises of biotechnology are often repeatedly postponed. This indicates that life is often more complicated than anticipated and that new paradigms might be needed to generate useful knowledge of biological systems.

The increasing computing power are emerging as a decisive factor for modern biology. The effect of faster computation have already transformed biology enabling the vast field of bioinformatics and as computing continues to improve its impact on biology are likely to be even greater. More advanced modeling and simulations also known as in silico experiments are enabled, which may guide researchers and designers of biological systems. To construct accurate models or simulations it necessary to have a framework representing the actual system. In terms of mathematics this framework could be a network of interactions and a set of parameters governing the magnitudes of interaction. In developing tools for synthetic biology, such as the trap-RNA system, is important to provide means of incorporating the tools into both living organisms as well as advanced modeling and simulation.

Purpose of modeling

In this project we developed a model of a trap-RNA system with the dual purpose of both developing hypotheses about the system and providing synthetic biologists the framework to incorporate the trap-RNA systems into advanced modeling of their designs. Furthermore we provide the ability to predict the caused gene repression based on key parameters some of which are changeable by altering the primary sequences involved...

biologist who wish to employ a trap-RNA system with the means of how to assess the tunability of gene repression by altering...

General kinetic model

The modeling of the trap-RNA system is based on reaction kinetics and multiple models are proposed based on two general reaction schemes


\begin{eqnarray} \leftrightharpoons_k^l \label{eq:r1} &\ce{\textit{\color{blue}m} + \textit{\color{red}s} <=>[k_{1,s}][k_{-1,s}] $c_{ms}$ ->[k_{2,s}]} (1 - p_s) {\color{red}s} \\ \label{eq:r2} &\ce{\textit{\color{red}s} + \textit{\color{green}r} <=>[k_{1,r}][k_{-1,r}] $c_{sr}$ ->[k_{2,r}]} (1-p_r) {\color{green}r} \end{eqnarray}

In the first reaction ${\color{blue} m}$RNA binds to a ${\color{red}s}$RNA forming a RNA:RNA complex called $c_{ms}$. The RNAs of the duplex is then irreversibly degraded in an RNAse-E dependent reaction with stoichiometries defined by $p_s$, which denotes the probability that $s$ is codegraded in the reaction. With this definition $(1 - p_s) {\color{red}s}$ represents sRNA that is not codegraded or released following the degradation of mRNA. The majority of investigated small regulatory RNA (srRNA) acts stoichiometrically, they are degraded 1:1 with their target mRNA corresponding to $p_s = 1$. Interestingly studies indicate that the trap-RNA system acts catalytically with $p_s$ close to zero[{\color{red}citation}], leading to different regulatory properties some of which we try to explore mathematically. For the second reaction the main experimental observation is that the t{\color{green}r}ap RNA inhibits the activity of the ${\color{red}s}$RNA but the mechanism is largely unknown. The general scheme allows different hypotheses; either trap-RNA works by competing for the sRNA binding site or it could mediate RNAse-E dependent degradation of the sRNA either catalytically or stoichiometrically.

\emph{In vivo} genes end their RNAs are constantly being expressed and degraded giving rise to finite lifetimes of RNA molecules. To model the trap-RNA system in context of living cells the expression of each RNA is described by a production term $\alpha_i$ and the degradation is described by a first order degradation term $\beta_i[RNA]_i$. \begin{equation} \label{proddeg} \ce{->[production] [RNA]_i ->[degradation]} \end{equation} Eventually the amount of RNA will settle into a steady state where production equals degradation providing stability for mathematical analysis.

\begin{eqnarray} \label{test} y & = & ax^2 + bx + c \\ f(x) & = & x^2 + 2xy + y^2 \end{eqnarray} \eqref{test}