Difference between revisions of "Team:Peshawar/Model"

(Prototype team page)
 
 
(11 intermediate revisions by the same user not shown)
Line 1: Line 1:
{{Peshawar}}
 
 
<html>
 
<html>
  
 +
<head>
 +
    <link rel="icon" type="image/png" href="https://static.igem.org/mediawiki/2016/f/f5/T-Peshawar--favicon-v1-compressor.png">
 +
    <title>Modeling - iGEM Peshawar 2016</title>
 +
    <!--Mind the order of loader and CSS-->
 +
    <script type="text/javascript" src="https://www.gstatic.com/charts/loader.js"></script>
 +
<link rel="stylesheet" type="text/css" href="https://2016.igem.org/Template:Peshawar/CSS?action=raw&ctype=text/css" />
 +
    <script type="text/x-mathjax-config">
 +
        MathJax.Hub.Config({
 +
        tex2jax: {
 +
          inlineMath: [ ['$','$'], ["\\(","\\)"] ],
 +
          processEscapes: true
 +
        }
 +
        });
 +
    </script>
 +
    <script type="text/javascript" src="http://cdn.mathjax.org/mathjax/latest/MathJax.js?config=TeX-AMS-MML_HTMLorMML"></script>
  
<div class="column full_size judges-will-not-evaluate">
+
</head>
<h3>★  ALERT! </h3>
+
<p>This page is used by the judges to evaluate your team for the <a href="https://2016.igem.org/Judging/Awards#SpecialPrizes">Best Model award</a>. </p>
+
  
 +
<body>
 +
    <div class="outermost-box">
 +
        <div id="menu-button" class="floating-button-wrapper">
 +
    <div id="menu-background" class="js-toggle-left-slidebar">
 +
    </div>
 +
</div>
 +
<!--End of floating menu button wrapper -->        <div off-canvas="slidebar-1 left reveal">
 +
    <nav>
 +
        <ul class="navigation-list">
 +
            <li><a href="https://2016.igem.org/Team:Peshawar">Home</a></li>
 +
            <hr>
 +
            <li><a href="https://2016.igem.org/Team:Peshawar/Description">Project</a></li>
 +
            <li><a href="https://2016.igem.org/Team:Peshawar/Model">Modeling</a></li>
 +
            <li><a href="https://2016.igem.org/Team:Peshawar/Notebook">Notebook</a></li>
 +
            <li><a href="https://2016.igem.org/Team:Peshawar/Safety">Safety</a></li>
 +
            <li><a href="https://2016.igem.org/Team:Peshawar/Parts">Parts</a></li>
 +
            <li><a href="https://2016.igem.org/Team:Peshawar/Basic_Part">Basic Part</a></li>
 +
            <li><a href="https://2016.igem.org/Team:Peshawar/Composite_Part">Composite Part</a></li>
 +
            <hr>
 +
            <li><a href="https://2016.igem.org/Team:Peshawar/Team">Team</a></li>
 +
            <li><a href="https://2016.igem.org/Team:Peshawar/Attributions">Attributions</a></li>
 +
            <li><a href="https://2016.igem.org/Team:Peshawar/Collaborations">Collaborations</a></li>
 +
            <hr>
 +
            <li><a href="https://2016.igem.org/Team:Peshawar/Human_Practices">Human Practices</a></li>
 +
            <li>
 +
                <a href="https://2016.igem.org/Team:Peshawar/HP/Silver">Silver</a>
 +
                <span>  |  </span>
 +
                <a href="https://2016.igem.org/Team:Peshawar/HP/Gold">Gold</a>
 +
            </li>
 +
            <li><a href="https://2016.igem.org/Team:Peshawar/Integrated_Practices">Integrated Practices</a></li>
 +
            <hr>
 +
            <li><a href="https://2016.igem.org/Team:Peshawar/Engagement">Public Engagement</a></li>
 +
            <li><a href="https://2016.igem.org/Team:Peshawar/Awards">Awards</a></li>
 +
            <hr>
 +
            <li><a href="https://2016.igem.org/Team:Peshawar/Wiki">How to design an iGEM Wiki</a></li>
  
<p> Delete this box in order to be evaluated for this medal. See more information at <a href="https://2016.igem.org/Judging/Pages_for_Awards/Instructions"> Instructions for Pages for awards</a>.</p>
+
        </ul>
 +
    </nav>
 +
    <article id="social-main-footer">
 +
    <!--FACEBOOK-->
 +
    <a href="https://www.facebook.com/iGEMPeshawar2016/" target="_blank">
 +
        <svg class="social-main-footer-svg facebook" width="61" height="60" viewBox="0 0 61 60" xmlns="http://www.w3.org/2000/svg">
 +
            <title>
 +
                T--Peshawar--facebook-round-v2
 +
            </title>
 +
            <g id="Page-1" fill="none" fill-rule="evenodd">
 +
                <g id="T--Peshawar--facebook-round-v2" fill="#808184">
 +
                    <path d="M30.594 1.465c-15.735 0-28.536 12.8-28.536 28.537 0 15.734 12.8 28.536 28.536 28.536S59.13 45.736 59.13 30.002c0-15.736-12.8-28.537-28.536-28.537m0 58.536c-16.542 0-30-13.455-30-29.997 0-16.543 13.458-30 30-30s30 13.457 30 30c0 16.542-13.458 30-30 30" id="facebook-circle" />
 +
                    <path d="M22.91 25.905h3.72V22.29c0-1.597.042-4.057 1.2-5.58 1.22-1.614 2.895-2.71 5.777-2.71 4.696 0 6.672.67 6.672.67l-.933 5.514s-1.55-.45-2.997-.45c-1.447 0-2.743.52-2.743 1.966v4.205h5.936l-.414 5.386h-5.524V50H26.63V31.29h-3.72v-5.385z" id="facebook" />
 +
                </g>
 +
            </g>
 +
        </svg>
 +
    </a>
 +
    <!--TWITTER-->
 +
    <a href="https://twitter.com/igem_peshawar" target="_blank">
 +
        <svg class="social-main-footer-svg twitter" width="60" height="60" viewBox="0 0 60 60" xmlns="http://www.w3.org/2000/svg">
 +
            <title>
 +
                T--Peshawar--twitter-round-v2
 +
            </title>
 +
            <g id="Page-1" fill="none" fill-rule="evenodd">
 +
                <g id="T--Peshawar--twitter-round-v2" fill="#808184">
 +
                    <path d="M30 1.463C14.266 1.463 1.464 14.265 1.464 30S14.264 58.537 30 58.537 58.538 45.735 58.538 30 45.738 1.463 30 1.463M30 60C13.46 60 0 46.542 0 30S13.46 0 30 0c16.542 0 30 13.458 30 30S46.542 60 30 60" id="twitter-circle" />
 +
                    <path d="M49.642 21.783c-1.325.587-2.75.984-4.242 1.163 1.527-.914 2.698-2.36 3.247-4.085-1.428.848-3.006 1.46-4.692 1.79-1.343-1.433-3.265-2.33-5.39-2.33-4.077 0-7.384 3.308-7.384 7.385 0 .58.065 1.142.19 1.683-6.134-.31-11.576-3.247-15.222-7.718-.636 1.093-.998 2.362-.998 3.715 0 2.56 1.303 4.82 3.283 6.148-1.208-.038-2.35-.373-3.347-.923v.09c0 3.58 2.547 6.566 5.927 7.242-.62.174-1.27.26-1.946.26-.476 0-.94-.044-1.39-.13.94 2.933 3.666 5.07 6.9 5.126-2.53 1.982-5.715 3.163-9.173 3.163-.596 0-1.185-.032-1.762-.1 3.27 2.093 7.15 3.316 11.32 3.316 13.59 0 21.015-11.25 21.015-21.013 0-.32-.006-.64-.02-.955 1.444-1.043 2.7-2.344 3.685-3.824" id="twitter" />
 +
                </g>
 +
            </g>
 +
        </svg>
 +
    </a>
 +
    <!--INSTAGRAM-->
 +
    <a href="https://www.instagram.com/igempeshawar/" target="_blank">
 +
        <svg class="social-main-footer-svg instagram" width="61" height="61" viewBox="0 0 61 61" xmlns="http://www.w3.org/2000/svg">
 +
            <title>
 +
                T--Peshawar--instagram-round-v2
 +
            </title>
 +
            <g id="Page-1" fill="none" fill-rule="evenodd">
 +
                <g id="T--Peshawar--instagram-round-v2" fill="#808184">
 +
                    <path d="M30.802 1.64c-15.736 0-28.537 12.802-28.537 28.537s12.8 28.537 28.537 28.537c15.735 0 28.536-12.802 28.536-28.537S46.538 1.64 30.802 1.64m0 58.537c-16.542 0-30-13.458-30-30s13.458-30 30-30 30 13.458 30 30-13.458 30-30 30" id="instagram-circle" />
 +
                    <path d="M45.3 42.212c0 1.902-1.548 3.446-3.45 3.446H19.763c-1.903 0-3.45-1.544-3.45-3.446v-14.73h5.375c-.464 1.144-.727 2.39-.727 3.695 0 5.426 4.417 9.837 9.845 9.837 5.428 0 9.843-4.41 9.843-9.837 0-1.306-.264-2.55-.725-3.694H45.3v14.73zm-19.63-14.73c1.15-1.593 3.022-2.64 5.137-2.64 2.115 0 3.985 1.047 5.136 2.64.75 1.042 1.203 2.314 1.203 3.695 0 3.492-2.845 6.33-6.34 6.33-3.494 0-6.336-2.838-6.336-6.33 0-1.38.45-2.653 1.2-3.694zm18.174-10.145l.795-.004v6.097l-6.083.022-.022-6.098 5.308-.017zM41.85 13.19H19.763c-3.835 0-6.957 3.118-6.957 6.952v22.07c0 3.834 3.122 6.952 6.957 6.952H41.85c3.834 0 6.957-3.118 6.957-6.952v-22.07c0-3.834-3.123-6.952-6.958-6.952z" id="instagram" />
 +
                </g>
 +
            </g>
 +
        </svg>
 +
    </a>
 +
    <!--BLOG-->
 +
    <a href="http://www.synbiokp.com/blog/" target="_blank">
 +
        <svg class="social-main-footer-svg blog" width="60" height="60" viewBox="0 0 60 60" xmlns="http://www.w3.org/2000/svg">
 +
            <title>
 +
                T--Peshawar--blog-round-v2
 +
            </title>
 +
            <g id="Page-1" fill="none" fill-rule="evenodd">
 +
                <g id="T--Peshawar--blog-round-v2" fill="#808184">
 +
                    <path d="M30 1.463C14.265 1.463 1.463 14.263 1.463 30c0 15.735 12.802 28.536 28.537 28.536S58.537 45.736 58.537 30 45.735 1.463 30 1.463M0 30C0 13.458 13.458 0 30 0s30 13.458 30 30-13.458 30-30 30S0 46.542 0 30z" id="blog-circle" />
 +
                    <path d="M34.95 19.916c7.654 0 14 6.357 14 14.025 0 7.64-6.375 14.03-14 14.03h-7.942C19.393 47.968 13 41.57 13 33.94V15.22c0-1.78 1.47-3.255 3.25-3.255 1.768 0 3.247 1.482 3.247 3.256 0 6.24.01 12.48.01 18.72 0 4.11 3.4 7.52 7.5 7.52h7.943c4.084 0 7.504-3.424 7.504-7.514 0-4.108-3.403-7.516-7.504-7.516h-8.568c-1.777 0-3.248-1.475-3.248-3.255 0-1.773 1.477-3.254 3.248-3.254h8.568z" id="blog" />
 +
                </g>
 +
            </g>
 +
        </svg>
 +
    </a>
 +
</article>
 
</div>
 
</div>
 +
<!--End of off-canvas menu-->        <div canvas="container">
 +
            <section class="inner-container">
 +
                <div id="header-oc" class="header-outer-container">
 +
    <div class="sticky-header">
 +
        <div class="logo-wrapper">
 +
            <a href="https://2016.igem.org/Team:Peshawar">
 +
                <img id="id-logo" class="logo" src="https://static.igem.org/mediawiki/2016/7/7b/T--Peshawar--iGEM-Peshawar-Logo---New-Pakol-v2-2.png">
 +
            </a>                           
 +
        </div>
 +
    </div>
 +
</div>
 +
<!--End of full width header containing responsive logo  -->                <script type="text/javascript">
 +
 +
    var button = document.getElementById('menu-button');
 +
    button.style.position = "fixed";
 +
    button.style.zIndex = "-1";
 +
    var buttonImage = document.getElementById("menu-background");
 +
 +
    //Wait for the Logo to load and fill he div
 +
    document.getElementById('id-logo').onload = function() {
 +
        var header = document.getElementById('header-oc');
 +
        headerHeight = header.scrollHeight;
 +
        screenWidth = header.scrollWidth;
 +
        console.log('headerHeight',headerHeight, 'headerwidth', screenWidth);
 +
 +
        //Scale the button according to header height
 +
        var buttonheight = headerHeight / 3.5;
 +
        buttonImage.style.backgroundSize = buttonheight + 'px';
 +
        buttonImage.style.height = buttonheight +'px';
 +
        buttonImage.style.width = buttonheight*1.25 +'px';
 +
        button.style.zIndex = "2";
 +
 +
        // Calculate button coordinates
 +
        menuButtonHeight = 80;                       
 +
        var buttonTopPosition = headerHeight / 2 - buttonheight / 2;
 +
        var buttonLeftPosition = screenWidth / 36 ;
 +
 +
        // Poistion the button at optimal location
 +
        button.style.left = buttonLeftPosition + 'px';
 +
        button.style.top = buttonTopPosition + 'px';
 +
    };
 +
 +
    // Hack for menu button not appear when the back button is pressed in history
 +
    window.onpageshow = function(event) {
 +
        var button = document.getElementById('menu-button');
 +
        button.style.position = "fixed";
 +
        button.style.zIndex = "-1";
 +
        var buttonImage = document.getElementById("menu-background");
 +
       
 +
        var header = document.getElementById('header-oc');
 +
        headerHeight = header.scrollHeight;
 +
        screenWidth = header.scrollWidth;
 +
        console.log('headerHeight',headerHeight, 'headerwidth', screenWidth);
 +
 +
        //Scale the button according to header height
 +
        var buttonheight = headerHeight / 3.5;
 +
        buttonImage.style.backgroundSize = buttonheight + 'px';
 +
        buttonImage.style.height = buttonheight +'px';
 +
        buttonImage.style.width = buttonheight*1.25 +'px';
 +
        button.style.zIndex = "2";
 +
 +
        // Calculate button coordinates
 +
        menuButtonHeight = 80;                       
 +
        var buttonTopPosition = headerHeight / 2 - buttonheight / 2;
 +
        var buttonLeftPosition = screenWidth / 36 ;
 +
 +
        // Poistion the button at optimal location
 +
        button.style.left = buttonLeftPosition + 'px';
 +
        button.style.top = buttonTopPosition + 'px';
 +
    };
 +
 +
</script>
 +
                <div class="full-content-wrapper-container">
 +
                    <div class="full-width-content-wrapper">
 +
                    </div>
 +
                    <!--End of full-width-content-wrapper-->
 +
                </div>
 +
 +
                <div class="full-content-wrapper-container">
 +
                    <div class="full-width-content-wrapper modeling-full">
 +
                        <div id="js-parallax-window" class="parallax-window">
 +
                          <div class="parallax-static-content">
 +
                            <h1>Modeling and Simulation</h1>
 +
                            <div class="video-overlay" style="background-image: url(&quot;https://d3k5xyayaartr5.cloudfront.net/_assets/pattern-overlays/patterns/black-medium-checks.png&quot;);"></div>
 +
                          </div>
 +
                          <div id="js-parallax-background" class="parallax-background"></div>
 +
                        </div>
 +
                    </div>
 +
                    <!--End of full-width-content-wrapper-->
 +
                </div>
 +
                <!--End of full-width-content-wrapper-container-->
 +
 +
 +
 +
 +
 +
                <!--End of full-width-content-wrapper-container-->
 +
                <div class="centered-content-wrapper-container model">
 +
                    <div class="centered-content-wrapper">
 +
                        <h2>Why model?</h2>
 +
                        <p>Before building a part or a device in lab, it is probably a good idea to simulate it in a simulation package &mdash; such as SimBiology, Cello, or GenoCAD &mdash; to check if the part or device one intends to make would work the way one wants it to. Although simulation can't solve all the problems that one might encounter during the actual physical design, it can certainly reduce them, and help in getting a better understanding of how a design behaves under different parameters and conditions.</p>
 +
 +
                        <p>Our project has two biological ystems: one for detection of corbon monoxide, and another for detection for oxides of nitrogen. Below we show you how we modelled and simulated both these systems.</p>
 +
                       
 +
                        <h2>Modeling of the CO-sensing system</h2>
 +
 +
                        <p>In our CO-sensing system, the CO binds to CooA (an activator). It activates the activator protein, and so the activated activator protein (CooA)2CO binds to promoter and initiates transcription.</p>
 +
                       
 +
                        <p>The reactions which take place can be given as:</p>
 +
 +
                        $$\textrm{[CO]}+\textrm{[CooA]} \longleftrightarrow \textrm{[(CooA)2CO]} \label{R1} \tag{R1}$$
 +
 +
                        <p>The equation above represents the first reaction that takes place in our system in which the CO present in the environment is taken up into the cell where it binds to the CooA protein. The CooA protein is constantly expressed in our system using a constitutive promoter. The coefficient 2 represents the number of CO which bind to CooA; this number was taken from the literature.</p>
 +
                       
 +
                        <p>To develop a dissociation constant of this equation, we can write it as:</p>
 +
 +
                        $$KI=\frac{\textrm{[CO].[CooA]}}{\textrm{[(CooA)2CO]}} \label{1} \tag{1}$$
 +
 +
                        <p>Note that the concentrations of $\textrm{[CO]}$, $\textrm{[CooA]}$, and $\textrm{[(CooA)2CO]}$ are at equilibrium before the start of the reaction R1.</p>
 +
                        <p>The second reaction that takes place in our system is the binding of the activated activator protein complex of $\textrm{[(CooA)2CO]}$ to the promoter Pr.</p>
 +
 +
                        <p>It can be described as:</p>
 +
 +
                        $$\textrm{[(CooA)2CO]} \longleftrightarrow \textrm{[(CooA)2CO.Pr]} \label{R2} \tag{R2}$$
 +
                       
 +
                        <p>Here, $\textrm{[(CooA)2CO.Pr]}$ represents the activated protein $\textrm{[(CooA)2CO]}$ bound to the promoter.</p>
 +
 +
                        <p>To develop a dissociation constant of this equation, we can write as:</p>
 +
                       
 +
                        $$KP=\frac{\textrm{[(CooA)2CO]}}{\textrm{[(CooA)2CO.Pr]}} \label{2} \tag{2}$$
 +
 +
                        <p>Now, in order to get an equation which describes the change in the concentration of mRNA with time we develop the following differential equation:</p>
 +
                        $$\frac{d(mRNA)}{dt}=k_{on1}.P_{exp}-\gamma_1.[mRNA] \label{3} \tag{3}$$
 +
 +
                        <p>where, $k_{on1}$ is represents the formation rate constant of the mRNA, and $\gamma_1$ represents the degradation rate constant of mRNA.</p>
 +
 +
                        <p>In Eq. (\ref{3}) it can be seen that the change in the $\textrm{[mRNA]}$ with time is equal to the difference of the formation rate of mRNA and  its degradation rate.</p>
 +
 +
                        <p>The initial conditions for the Eq. (\ref{3}) are assumed to be:</p>
 +
 +
                        $$\textrm{[mRNA(t=0)] = 0 molar}$$
 +
 +
                        <p>This means that at the time $t=0$ sec (initial time), there are no molecules of mRNA encoding for our desired chromoprotein in the cell.</p>
 +
 +
                        <p>In Eq. (\ref{3}), $P_{exp}$ represents the probability that the promoter will be activated by the activated activator protein complex $\textrm{[(CooA)2CO]}$. It can be calculated by the following equation:</p>
 +
                        $$P_{exp} = \frac{\textrm{number of situations in which the promoter will be activated}}{\textrm{total number of situations}} = \frac{w}{z} \label{4} \tag{4}$$
 +
 +
                        <p>In our modeling of the CO-sensing system, we have assumed that at any moment only one molecule of the activated activator protein complex $\textrm{[(CooA)2CO]}$ can bind to the promoter, and that the \ref{R1} reaches equilibrium before the transcription of the system begins.</p>
 +
 +
                        <p>To develop our equations further, we say that the number of situations in which a promoter will be activated is equal to:</p>
 +
                        $$ w=\textrm{[(CooA)2CO.Pr]}=\frac{\textrm{[(CooA)2CO]}}{K_p} \label{5} \tag{5}$$
 +
 +
                        <p>Eq. (\ref{5}) demonstrates the relation between the binding of the activated activator protein complex $\textrm{[(CooA)2CO]}$, showing that it is equal to the amount of $\textrm{[(CooA)2CO]}$ upon the dissociation constant of the binding of $\textrm{[(CooA)2CO]}$ to the promoter.</p>
 +
 +
                        <p>To simplify this equation, we assume that the number of situations in which the promoter will not be activated by $\textrm{[(CooA)2CO]}$ is equal to $1$.</p>
 +
 +
                        <p>Note that this assumption is made because the only situation in which transcription will not take place is when there are no molecules bound to the promoter.</p>
 +
 +
                        <p>Substituting the assumptions in Eq. (\ref{5}), we get:</p>
 +
                        $$ z=w+1=\frac{\textrm{[(CooA)2CO]}}{K_p}+1  \label{6} \tag{6}$$
 +
                        <p>where, $z$ represents the total number of situations, $w$ represents the number of situations in which the promoter will be activated, and $1$ represents the number of situations in which the promoter will not be activated.</p>
 +
                        <p>Sustituting Eq. (\ref{5}) and (\ref{6}) into Eq. (\ref{4}), we get:</p>
 +
 +
                        $$P_{exp} = \frac{w}{z} = \frac{\frac{\textrm{[(CooA)2CO]}}{K_p}}{\frac{\textrm{[(CooA)2CO]}}{K_p}+1} \label{7} \tag{7}$$
 +
 +
                        <p>As we want our $P_{exp}$ to be expressed in the terms of $\textrm{[CO]}$ and $\textrm{[CooA]}$ separately instead of $\textrm{[(CooA)2CO]}$, we repeat these above equations by taking into consideration \ref{R1} and its dissociation constant $\textrm{KI}$.</p>
 +
 +
                        <p>By doing so, we get:</p>
 +
                        $$P_{exp} = \frac{\frac{\textrm{[CooA]}.\frac{\frac{\textrm{CO}}{K_I}}{1+\frac{\textrm{CO}}{K_I}}}{K_p}} {1+\frac{\textrm{[CooA]}.\frac{\frac{\textrm{CO}}{K_I}}{1+\frac{\textrm{CO}}{K_I}}}{K_p}} \label{8} \tag{8}$$
 +
 +
                        <p>As the Eq. (\ref{8}) demonstrates, $P_{exp}$ depends only upon the concentration of $\textrm{[CO]}$ and $\textrm{[CooA]}$.</p>
 +
                        <p>Recall Eq. (\ref{3})</p>
 +
                        $$\frac{d(\textrm{mRNA})}{dt}=k_{on1}.P_{exp}-\gamma_1.[\textrm{mRNA}]$$
 +
 +
                        <p>The initial condition for this equation was $\textrm{[mRNA(t=0)] = 0 molar}$.</p>
 +
 +
                        <p>The analytical solution of this differential equation is equal to:</p>
 +
                        $$[mRNA]=\frac{k_{on1}.P_{exp}}{\gamma_1}(1-e^{-\gamma_1t})$$
 +
 +
                        <p>Now that we have got the equation describing the change in concentration of mRNA, we can use the same method to derive an equation which describes the change in the concentration of the chromoprotein in time as following:</p>
 +
                        $$\frac{d\textrm{[CP]}}{dt} = k_{on2}.\textrm{[mRNA]}-\gamma_2\textrm{[CP]} \label{9} \tag{9}$$
 +
 +
                        <p>Here, $k_{on2}$  is the formation rate constant of the chromoprotein and $\gamma_2$  is the degradation rate constant of the chromoprotein.</p>
 +
 +
                        <br>
 +
                        <figure id="co-block">
 +
                            <img src="https://static.igem.org/mediawiki/2016/f/f2/T--Peshawar--co-block-diagram.png">
 +
                            <figcaption>
 +
                            <b>Figure 1: </b> The central theme considered for mathematical modeling of CO system, where the chromoprotein DNA leads to the production of mRNA of chromoprotein and then subsequently of chromoprotein. $K_{on1}$ is the rate of transcription of chromoprotein DNA, while $K_{on2}$ is the translation rate of mRNA of chromoprotein. Also, $\gamma_1$ is the degradation rate of mRNA of cheomoprotein and $\gamma_2$ is the degradation rate of chromoprotein.           
 +
                            </figcaption>
 +
                        </figure>
 +
                        <br>
 +
 +
 +
                        <p>Again, the initial conditions for this equation are assumed to be:</p>
 +
                        $$\textrm{[CP(t=0)] = 0 molar}$$
 +
 +
                        <p>which states that at the initial time, the amount of chromoprotein in the cell is equal to $0$.</p>
 +
 +
                        <p>Here on forth, we will repeat the steps from Eq. (\ref{4}) to Eq. (\ref{8}) and arrive to the analytical solution of the differential equation of Eq. (\ref{9}), which is:</p>
 +
                        $$\textrm{[CP]} = k_{on2}.\frac{k_{on1}.P_{exp}}{\gamma_1.\gamma_2}(1-e^{-\gamma_2t}) + k_{on2}.\frac{k_{on1}.P_{exp}}{\gamma_1(\gamma_1-\gamma_2)}(e^{-\gamma_2t}-e^{-\gamma_1t}) \label{10} \tag{10}$$
 +
 +
                        <p>Eq. (\ref{10}) is the final equation which describes the amount of $\textrm{[CP]}$ along with the amount of $\textrm{[mRNA]}$.</p>
 +
 +
                        <h3>Assumptions</h3>
 +
 +
                        <ol>
 +
                            <li>1. Reactions \ref{R1} and \ref{R2} reach equilibrium before the transcription begins.</li>
 +
                            <li>2. Only one molecule of the complex $\textrm{[(CooA)2CO]}$ can be bound to the promoter at a time.</li>
 +
                            <li>3. $\textrm{[mRNA(t=0)] = 0 molar}$</li>
 +
                            <li>4. $\textrm{[CO]}$ and $\textrm{[CooA]}$ concentrations inside the cell are constant.</li>
 +
                        </ol>
 +
 +
                        <h2>Simulation of the CO-sensing system</h2>
 +
                        <h3>Simulation parameters</h3>                           
 +
                        <br>
 +
                        <table class="table-minimal">
 +
                            <caption><b>Table 1:</b> Parameters, their description, and their values used in MATLAB simulation of CO-sensing system.</caption>
 +
                            <thead>
 +
                                <tr>
 +
                                    <th>Parameter</th>
 +
                                    <th>Description</th>
 +
                                    <th>Value</th>
 +
                                </tr>
 +
                            </thead>
 +
                            <tbody>
 +
                                <tr>
 +
                                    <td>$\textrm{[CO]}$</td>
 +
                                    <td>Carbon monoxide concentration at equilibrium</td>
 +
                                    <td>$\textrm{10}$</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>$\textrm{[CooA]}$</td>
 +
                                    <td>CooA protein concentration at equilibrium</td>
 +
                                    <td>$11$</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>$\textrm{KI}$</td>
 +
                                    <td>Dissociation constant of the binding of CO to CooA</td>
 +
                                    <td>$\textrm{0.011 mM}$</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>$\textrm{KP}$</td>
 +
                                    <td>Dissociation Constant of the binding of $\textrm{(CO2CooA)}$ to the promoter</td>
 +
                                    <td>$\textrm{0.000009}$ $\textrm{min}^{-1}$</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>Kon1</td>
 +
                                    <td>Formation rate constant of mRNA</td>
 +
                                    <td>$\textrm{0.001}$ $\textrm{min}^{-1}$</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>Kon2</td>
 +
                                    <td>Formation rate constant of the chromoprotein</td>
 +
                                    <td>0.006 $\textrm{min}^{-1}$ / 0.0051 $\textrm{min}^{-1}$</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>Gama1 ($\gamma$)</td>
 +
                                    <td>Degradation rate constant of the mRNA</td>
 +
                                    <td>$\textrm{0.00039}$ $\textrm{min}^{-1}$</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>Gama2 ($\gamma$)</td>
 +
                                    <td>Degradation rate constant of the chromoprotein</td>
 +
                                    <td>$\textrm{0.001}$ $\textrm{min}^{-1}$</td>
 +
                                </tr>
 +
                            </tbody>
 +
                        </table>
 +
                        <br>
 +
 +
                        <h3>Simulation results</h3>
 +
                        <br><br>
 +
 +
                        <figure>
 +
                        <!--Demogrpahics-->
 +
<script type="text/javascript">
 +
google.charts.load('current', {
 +
    packages: ['corechart', 'line']
 +
});
 +
google.charts.setOnLoadCallback(drawCurveTypes);
 +
 +
function drawCurveTypes() {
 +
    var data = new google.visualization.DataTable();
 +
    data.addColumn('number', 'Time');
 +
    data.addColumn('number', "Conc. of chromoprotein with pCooF promoter");
 +
    data.addColumn('number', "Conc. of chromoprotein with pCooM promoter");
 +
    data.addColumn('number', "Conc. of mRNA");
 +
 +
    data.addRows([
 +
        [0, 0, 0, 0],
 +
        [0.00084492, 7.64009e-09, 6.49407e-09, 5.02372e-05],
 +
        [0.0016898, 3.05591e-08, 2.59752e-08, 0.000100473],
 +
        [0.0025347, 6.87552e-08, 5.84419e-08, 0.000150709],
 +
        [0.0033797, 1.22227e-07, 1.03893e-07, 0.000200943],
 +
        [0.0076042, 6.18647e-07, 5.2585e-07, 0.000452099],
 +
        [0.011829, 1.49667e-06, 1.27217e-06, 0.000703231],
 +
        [0.016053, 2.75607e-06, 2.34266e-06, 0.000954338],
 +
        [0.020278, 4.3966e-06, 3.73711e-06, 0.00120542],
 +
        [0.041401, 1.83083e-05, 1.55621e-05, 0.00246046],
 +
        [0.062524, 4.1714e-05, 3.54569e-05, 0.00371487],
 +
        [0.083647, 7.45849e-05, 6.33972e-05, 0.00496867],
 +
        [0.10477, 0.000116892, 9.93584e-05, 0.00622185],
 +
        [0.21038, 0.000468979, 0.000398632, 0.0124785],
 +
        [0.316, 0.00105272, 0.000894814, 0.0187196],
 +
        [0.42161, 0.00186462, 0.00158493, 0.0249454],
 +
        [0.52723, 0.00290122, 0.00246604, 0.0311558],
 +
        [1.0048, 0.0103016, 0.00875638, 0.0590466],
 +
        [1.4823, 0.0219241, 0.0186355, 0.0866276],
 +
        [1.9599, 0.0374825, 0.0318601, 0.113902],
 +
        [2.4375, 0.0567059, 0.0482, 0.140873],
 +
        [3.2995, 0.0998687, 0.0848884, 0.188803],
 +
        [4.1616, 0.15278, 0.129863, 0.235775],
 +
        [5.0236, 0.21421, 0.182079, 0.281809],
 +
        [5.8856, 0.283037, 0.240581, 0.326924],
 +
        [7.1129, 0.391805, 0.333035, 0.389605],
 +
        [8.3402, 0.511068, 0.434408, 0.450511],
 +
        [9.5675, 0.638707, 0.542901, 0.509693],
 +
        [10.7948, 0.772827, 0.656903, 0.5672],
 +
        [12.4036, 0.95574, 0.812379, 0.640122],
 +
        [14.0124, 1.14424, 0.972602, 0.71035],
 +
        [15.6212, 1.33605, 1.13564, 0.777984],
 +
        [17.2299, 1.52913, 1.29976, 0.843118],
 +
        [19.2775, 1.77411, 1.50799, 0.922551],
 +
        [21.3251, 2.0164, 1.71394, 0.998267],
 +
        [23.3727, 2.2544, 1.91624, 1.07044],
 +
        [25.4203, 2.48659, 2.1136, 1.13924],
 +
        [28.0161, 2.77081, 2.35519, 1.22185],
 +
        [30.6119, 3.04327, 2.58678, 1.29959],
 +
        [33.2078, 3.30361, 2.80807, 1.37275],
 +
        [35.8036, 3.55121, 3.01853, 1.4416],
 +
        [39.1395, 3.85025, 3.27271, 1.52415],
 +
        [42.4754, 4.12908, 3.50972, 1.60051],
 +
        [45.8113, 4.38884, 3.73051, 1.67113],
 +
        [49.1472, 4.63002, 3.93552, 1.73645],
 +
        [53.5666, 4.92212, 4.1838, 1.81549],
 +
        [57.986, 5.18634, 4.40839, 1.88677],
 +
        [62.4054, 5.42549, 4.61167, 1.95104],
 +
        [66.8248, 5.64133, 4.79513, 2.009],
 +
        [72.9821, 5.90656, 5.02057, 2.0804],
 +
        [79.1394, 6.13655, 5.21607, 2.14221],
 +
        [85.2967, 6.33635, 5.3859, 2.1957],
 +
        [91.4539, 6.50928, 5.53289, 2.24202],
 +
        [98.9539, 6.68778, 5.68461, 2.29017],
 +
        [106.4539, 6.83798, 5.81228, 2.33056],
 +
        [113.9539, 6.96511, 5.92034, 2.36441],
 +
        [121.4539, 7.07154, 6.01081, 2.3928],
 +
        [128.9539, 7.15888, 6.08504, 2.41667],
 +
        [136.4539, 7.23294, 6.148, 2.43669],
 +
        [143.9539, 7.29714, 6.20257, 2.45346],
 +
        [151.4539, 7.35063, 6.24803, 2.46753],
 +
        [158.9539, 7.39199, 6.28319, 2.47936],
 +
        [166.4539, 7.42804, 6.31384, 2.48929],
 +
        [173.9539, 7.46188, 6.3426, 2.4976],
 +
        [181.4539, 7.48963, 6.36619, 2.50458],
 +
        [188.9539, 7.50688, 6.38085, 2.51044],
 +
        [196.4539, 7.52364, 6.3951, 2.51536],
 +
        [203.9539, 7.54382, 6.41225, 2.51948],
 +
        [211.4539, 7.55968, 6.42573, 2.52294],
 +
        [218.4274, 7.5641, 6.42948, 2.52565],
 +
        [225.4009, 7.5706, 6.43501, 2.52796],
 +
        [232.3744, 7.58281, 6.44539, 2.52992],
 +
        [239.3479, 7.59222, 6.45339, 2.53158],
 +
        [245.7054, 7.59267, 6.45377, 2.53288],
 +
        [252.0629, 7.59504, 6.45578, 2.534],
 +
        [258.4203, 7.60131, 6.46111, 2.53497],
 +
        [264.7778, 7.6063, 6.46535, 2.5358],
 +
        [271.1393, 7.60623, 6.46529, 2.53651],
 +
        [277.5008, 7.60741, 6.4663, 2.53713],
 +
        [283.8624, 7.61107, 6.46941, 2.53766],
 +
        [290.2239, 7.61396, 6.47186, 2.53812],
 +
        [292.6679, 7.61412, 6.472, 2.53828],
 +
        [295.1119, 7.61436, 6.47221, 2.53843],
 +
        [297.556, 7.61466, 6.47246, 2.53857],
 +
        [300, 7.61498, 6.47273, 2.5387]
 +
 +
    ]);
 +
 +
    var options = {
 +
        title: 'Chromoprotein concentration with different promoters (pCooF, pCooM)',
 +
        colors: ['#185FB1', '#702EA6', '#ED2B72'],
 +
        lineWidth: 5,
 +
 +
        hAxis: {
 +
            gridlines: {color: '#F1F1F1'},
 +
            title: 'Time (minutes)',
 +
            textStyle: {
 +
                fontSize: 16,
 +
                //fontName: 'Arial',
 +
                bold: false,
 +
                italic: false
 +
            },
 +
            titleTextStyle: {
 +
                fontSize: 17,
 +
                bold: false,
 +
                italic: false
 +
 +
            }
 +
        },
 +
        vAxis: {
 +
            gridlines: {color: '#F1F1F1'},
 +
            title: 'Concentration [nM]',
 +
            textStyle: {
 +
                fontSize: 16,
 +
                //fontName: 'Arial',
 +
                bold: false,
 +
                italic: false
 +
            },
 +
            titleTextStyle: {
 +
                fontSize: 17,
 +
                bold: false,
 +
                italic: false
 +
 +
            }
 +
        },
 +
        titleTextStyle: {
 +
            fontSize: 19
 +
        },
 +
 +
 +
 +
    };
 +
 +
    var chart = new google.visualization.LineChart(document.getElementById('chart16_div'));
 +
    chart.draw(data, options);
 +
}
 +
</script>
 +
<div id="chart16_div" style="width: 900px; height: 600px;"></div>
 +
                            <figcaption>
 +
                                <b>Figure 2:</b> Activity of the different promoters after stimulation by carbon monoxide. The strong promoter pCooF, shown in blue, is yielding higher level of expression than the weak promoter pCooM, shown in violet. The difference between the Strong promoter (pCooF) and the weak promoter (pCooM) is meaningful, and both of these promoters will express enough chromoprotein to be visual. Hence, our goal is achieved.
 +
                            </figcaption>
 +
                        </figure>
 +
                        <h3>Conclusion</h3>
 +
                        <p>It is clear from the results that our model will achieve the steady state and give us the amount of chromoprotein that can be visualized.</p>
 +
 +
                        <h3>MATLAB Code</h3>
 +
<pre>
 +
<code class="language-matlab">
 +
function dx = ODEfun_CooM(t,x)
 +
    dx      = zeros(1,2);
 +
    %x is a vector (x(1),x(2))
 +
    %x(1) = [mRNA]
 +
    %x(2) = [Chromoprotein]
 +
    KI      = 0.011;    %Dissociation constant of the binding of CO to CooA
 +
    KP      = 0.1;      %Dissociation constant of the binding of
 +
                        %(CO2CooA) to the promoter.
 +
    Kon1    = 0.001;    %Formation rate constant of mRNA
 +
    Kon2    = 0.0051;  %Formation rate constant of the CP
 +
    gama1  = 0.00039;  %Degradation rate constant of mRNA
 +
    gama2  = 0.002;    %Degradation rate constant of the CP
 +
    co      = 10;      %Concentration of CO
 +
    cooa    = 11;      %Concentration of CooA
 +
    w      = (cooa*co/(KI*(1+co/KI)))/KP;
 +
    z      = (cooa*co/(KI*(1+co/KI)))/KP+1;
 +
    P_exp  = w/z;
 +
    dx1    = (Kon1*P_exp)-(gama1*x(1));    %Rate equation of mRNA
 +
    dx2    = (Kon2*x(1))-(gama2*x(2));    %Rate equation of CP
 +
    dx      = [dx1;dx2];
 +
end
 +
 +
%Save this code as ODEfun_CooM.m file
 +
</code>
 +
</pre>
 +
 +
<pre>
 +
<code class="language-matlab">
 +
function dx = ODEfun_CooF(t,x)
 +
    dx      = zeros(1,2);
 +
    %x is a vector [x(1),x(2)]
 +
    %x(1) = [mRNA]
 +
    %x(2) = [chromoprotein]
 +
    KI      = 0.011;    %Dissociation constant of the binding of CO to CooA
 +
    KP      = 0.1;      %Dissociation constant of the binding of
 +
                        %(CO2CooA) to the promoter.
 +
    Kon1    = 0.001;    %Formation rate constant of mRNA
 +
    Kon2    = 0.006;    %Formation rate constant of the CP
 +
    gama1  = 0.00039;  %Degradation rate constant of mRNA
 +
    gama2  = 0.002;    %Degradation rate constant of the CP
 +
    co      = 10;      %Concentration of CO
 +
    cooa    = 11;      %Concentration of CooA
 +
    w      = (cooa*co/(KI*(1+co/KI)))/KP;
 +
    z      = (cooa*co/(KI*(1+co/KI)))/KP+1;
 +
    P_exp  = w/z;
 +
    dx1    = (Kon1*P_exp)-(gama1*x(1));    %Rate equation of mRNA
 +
    dx2    = (Kon2*x(1))-(gama2*x(2));    %Rate equation of CP
 +
    dx      = [dx1;dx2];
 +
end
 +
 +
%Save this code as ODEfun_CooF.m file
 +
</code>
 +
</pre>
 +
                       
 +
<pre>
 +
<code class="language-matlab">
 +
% Clear worskpace
 +
clc
 +
clear all
 +
% Code for calling the previous function:
 +
time_range      = 18000;
 +
[T1,X1]        = ode45(@ODEfun_CooF,[0 time_range],[0 0]);
 +
[T2,X2]        = ode45(@ODEfun_CooM,[0 time_range],[0 0]);
 +
figure
 +
plot(T1,X1(:,2),'r', T2,X2(:,2),'b', T1,X1(:,1),'g')
 +
legend('Conc. of chromoprotein with pCooF promoter ',...
 +
    'Conc. of chromoprotein with pCooM promoter',...
 +
    'mRNA concentration')
 +
xlabel('Time [sec]');
 +
ylabel('Concentration [molar]');
 +
title('Chromoprotein concentration with different promoters (pCooF, pCooM)')
 +
 +
%Save this code as co_sensing_system_simulation.m file
 +
</code>
 +
</pre>
 +
 +
                        <!--NOx circuit-->
 +
                        <h2>Modeling of the NOx-sensing system</h2>
 +
                        <p>To gain an insight into the stability of NO-sensing circuit, we consider the cooperativity of NO binding with NsrR and the law of mass action. The goal of the modeling is to determine the stable expression of chromoprotein by considering a specified concentration of different parameters of the dynamic system.</p>
 +
                        <h3>Rate of change of chromoprotein expression</h3>
 +
                        <p>The following ordinary differential equations are used for determining the expression levels of chromoproteins (CP) and its mRNA (mCP):</p>
 +
 +
                        $$\frac{d\textrm{[mCP]}}{dt} = k_{transcription}\textrm{[PyeaR]} - \gamma_m \textrm{[mCP]} \label{11} \tag{11}$$
 +
                        $$\frac{d\textrm{[CP]}}{dt}=k_{translation}\textrm{[mCP]}-\gamma_p\textrm{[CP]} \label{12} \tag{12}$$
 +
 +
                        <br>
 +
                        <figure id="nox-block">
 +
                            <img src="https://static.igem.org/mediawiki/2016/f/f7/T--Peshawar--nox-block-diagram.png">
 +
                            <figcaption>
 +
                            <b>Figure 3: </b> The core idea used in mathematical modeling of NOx system, where the activation of PyeaR promotor leads to the production of mRNA of chromoprotein and then subsequently of chromoprotein. $K_{transcription}$ is the transcription rate of PyeaR promotor, while $K_{translation}$ is the translation rate of mRNA of chromoprotein. Additionally, $\gamma_m$ is the degradation rate of mCP (mRNA of cheomoprotein) and $\gamma_p$ is the degradation rate of CP (chromoprotein). 
 +
                            </figcaption>
 +
                        </figure>
 +
                        <br>
 +
 +
                        <h3>Binding cooperativity betweeen NO AND NsrR</h3>
 +
                        <p>The cooperativity of $\textrm{NO}$ binding with NsrR can be found using the Hill Equation, where free $\textrm{NO}^*$ combines with the $\textrm{n}$ amount of NsrR, leading to the formation of $\textrm{n[NsrR][NO]}$ complex. The total amount of inducer $NO_T$ concentration can be described by the following equation:</p>
 +
                        $$\textrm{NO}_T = \textrm{n[NsrR][NO]} + [\textrm{NO}^*]$$
 +
 +
                        <p>Additionally, $\textrm{NsrR}$ dimerization with respect to unbound $\textrm{NsrR}$ can be represented by the following equation by considering the degree of cooperativity of $\textrm{NO}$ binding:</p>   
 +
                        $$\frac{[\textrm{NO}^*]}{[\textrm{NO}_T]} = \frac{1}{1+\left(\frac{NsrR}{K_{NO}} \right)^n} \label{13} \tag{13}$$
 +
 +
                        <h3>Equation for PyeaR activity</h3>
 +
                        <p>The correlation between PyeaR activity and rate of transcription of mCP (\beta) can represented as: </p>
 +
                        $$\textrm{[PyeaR]}activity = \frac{\beta}{1+ \left( \frac{\textrm{NO}^*}{K_{d(\textrm{NsrR})}} \right)} \label{14} \tag{14}$$
 +
 +
                        <h3>Solving equations for chromoprotein expressions</h3>
 +
                        <p>Combining the Eqs. (\ref{13}) and (\ref{14}) into equation (\ref{11}), yields:</p>
 +
                        $$\frac{d[\textrm{mCP}]}{dt} = \frac{k_{transcription}[\textrm{PyeaR}]}{1 + \left(\frac{[\textrm{NsrR}]}{\left(1+\left[\frac{[\textrm{NO}]}{K_{NO}}\right]^n \right)}\right) k_{d(\textrm{NsrR})} } - \gamma_m [\textrm{mCP}]$$
 +
 +
                        <p>In order to further simplify the system, let</p>
 +
                        $$K=\frac{k_{transcription}[\textrm{PyeaR}]}{1 + \left(\frac{[\textrm{NsrR}]}{\left(1+\left[\frac{[\textrm{NO}]}{K_{NO}}\right]^n \right)}\right) k_{d(\textrm{NsrR})} }$$
 +
                        <p>Additionally, the value of $\textrm{K}$ is inversely proportional to the amount of $\textrm{NsrR}$ and directly proportional to the $\textrm{NO}$ concentration in the system. Substituting $\textrm{K}$ into the differential equation of $\textrm{mCP}$ generates the following equation:</p>
 +
                        $$\frac{d[\textrm{mCP}]}{dt} = K - \gamma_m [\textrm{mCP}]$$
 +
                        <p>After solving the equation through homogenous and particular solution, the following equation is obtained for steady&mdash;state concentration of mCP generated after PyeaR is triggered:</p>
 +
                        $$[\textrm{mCP}]=\frac{K}{\gamma_m} - \frac{K}{\gamma_m}e^{-\gamma_mt}$$
 +
                        <p>The steady&mdash;state concentration of mCP is then substituted into Eq. (\ref{12}) and is also solved through homogenous and particular solution. The resulting equation gives the steady&mdash;state concentration of CP in the desired system according to the specified formation and degradation rates of CP in response to promotor stimulation:</p>
 +
                        $$[\textrm{CP}]=\frac{K k_{translation}}{\gamma_p \gamma_m} + \frac{K k_{translation}}{{\gamma_m}^2-\gamma_p \gamma_m} e^{-\gamma_mt} - \left(\frac{K k_{translation}}{\gamma_p \gamma_m} + \frac{K k_{translation}}{{\gamma_m}^2-\gamma_p \gamma_m} \right)e^{-\gamma_pt} $$
 +
                        <br><br>
 +
 +
                        <h2>Simulation of the NOx-sensing system</h2>
 +
                        <h3>Simulation parameters</h3>
 +
                        <table class="table-minimal">
 +
                            <caption><b>Table 2:</b> Parameters, their description, and their values used in MATLAB simulation of the NOx-sensing system.</caption>
 +
                            <thead>
 +
                                <tr>
 +
                                    <th>Parameter </th>
 +
                                    <th>Description</th>
 +
                                    <th>Value</th>
 +
                                </tr>
 +
                            </thead>
 +
                            <tbody>
 +
                                <tr>
 +
                                    <td>$[\textrm{PyeaR}]$</td>
 +
                                    <td>Concentration of PyeaR</td>
 +
                                    <td>1 nM</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>$k_{transcription}$</td>
 +
                                    <td>Rate of chromoprotein mRNA synthesis</td>
 +
                                    <td>0.167$\textrm{min}^{-1}$</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>$k_{translation}$</td>
 +
                                    <td>Rate of chromoprotein synthesis</td>
 +
                                    <td>$0.0011 \textrm{min}^{-1}$</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>$\gamma_m$</td>
 +
                                    <td>mRNA degradation rate</td>
 +
                                    <td>$0.19 \textrm{min}^{-1}$</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>$\gamma_p$</td>
 +
                                    <td>Chromoprotein degradation rate</td>
 +
                                    <td>$0.18 \textrm{min}^{-1}$</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>$K_{\textrm{NO}}$</td>
 +
                                    <td>Dissociation constant of NO</td>
 +
                                    <td>$0.12 \textrm{min}^{-1}$</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>$K_{\textrm{d(NsrR)}}$</td>
 +
                                    <td>Dissociation constant of NsrR</td>
 +
                                    <td>$\textrm{35}$ nM</td>
 +
                                </tr>
 +
                                <tr>
 +
                                    <td>n</td>
 +
                                    <td>Cooperativity of NO binding to NsrR</td>
 +
                                    <td>2</td>
 +
                                </tr>
 +
                            </tbody>
 +
                        </table>
 +
                       
 +
                        <h3>Simulation results</h3>
 +
                        <br><br>
 +
 +
                        <figure>
 +
                            <!--Demogrpahics-->
 +
<script type="text/javascript">
 +
google.charts.load('current', {
 +
    packages: ['corechart', 'line']
 +
});
 +
google.charts.setOnLoadCallback(drawBasic);
 +
 +
function drawBasic() {
 +
 +
 +
    var data = new google.visualization.DataTable();
 +
    data.addColumn('number', 'Time (minutes)');
 +
    data.addColumn('number', 'mCP per cell (nM)');
 +
 +
    data.addRows([
 +
        [0, 0],
 +
        [1, 0.15209],
 +
        [2, 0.27787],
 +
        [3, 0.38188],
 +
        [4, 0.46789],
 +
        [5, 0.53902],
 +
        [6, 0.59784],
 +
        [7, 0.64649],
 +
        [8, 0.68671],
 +
        [9, 0.71998],
 +
        [10, 0.74748],
 +
        [11, 0.77023],
 +
        [12, 0.78904],
 +
        [13, 0.8046],
 +
        [14, 0.81747],
 +
        [15, 0.82811],
 +
        [16, 0.8369],
 +
        [17, 0.84418],
 +
        [18, 0.85019],
 +
        [19, 0.85517],
 +
        [20, 0.85928],
 +
        [21, 0.86269],
 +
        [22, 0.8655],
 +
        [23, 0.86783],
 +
        [24, 0.86975],
 +
        [25, 0.87134],
 +
        [26, 0.87266],
 +
        [27, 0.87375],
 +
        [28, 0.87465],
 +
        [29, 0.87539],
 +
        [30, 0.87601],
 +
        [31, 0.87652],
 +
        [32, 0.87694],
 +
        [33, 0.87728],
 +
        [34, 0.87757],
 +
        [35, 0.87781],
 +
        [36, 0.87801],
 +
        [37, 0.87817],
 +
        [38, 0.8783],
 +
        [39, 0.87842],
 +
        [40, 0.87851],
 +
        [41, 0.87858],
 +
        [42, 0.87865],
 +
        [43, 0.8787],
 +
        [44, 0.87874],
 +
        [45, 0.87878],
 +
        [46, 0.87881],
 +
        [47, 0.87883],
 +
        [48, 0.87885],
 +
        [49, 0.87887],
 +
        [50, 0.87888],
 +
        [51, 0.87889],
 +
        [52, 0.8789],
 +
        [53, 0.87891],
 +
        [54, 0.87892],
 +
        [55, 0.87892],
 +
        [56, 0.87893],
 +
        [57, 0.87893],
 +
        [58, 0.87893],
 +
        [59, 0.87894],
 +
        [60, 0.87894],
 +
        [61, 0.87894],
 +
        [62, 0.87894],
 +
        [63, 0.87894],
 +
        [64, 0.87894],
 +
        [65, 0.87894],
 +
        [66, 0.87894],
 +
        [67, 0.87894],
 +
        [68, 0.87895],
 +
        [69, 0.87895],
 +
        [70, 0.87895],
 +
        [71, 0.87895],
 +
        [72, 0.87895],
 +
        [73, 0.87895],
 +
        [74, 0.87895],
 +
        [75, 0.87895],
 +
        [76, 0.87895],
 +
        [77, 0.87895],
 +
        [78, 0.87895],
 +
        [79, 0.87895],
 +
        [80, 0.87895],
 +
        [81, 0.87895],
 +
        [82, 0.87895],
 +
        [83, 0.87895],
 +
        [84, 0.87895],
 +
        [85, 0.87895],
 +
        [86, 0.87895],
 +
        [87, 0.87895],
 +
        [88, 0.87895],
 +
        [89, 0.87895],
 +
        [90, 0.87895],
 +
        [91, 0.87895],
 +
        [92, 0.87895],
 +
        [93, 0.87895],
 +
        [94, 0.87895],
 +
        [95, 0.87895],
 +
        [96, 0.87895],
 +
        [97, 0.87895],
 +
        [98, 0.87895],
 +
        [99, 0.87895],
 +
        [100, 0.87895],
 +
        [101, 0.87895],
 +
        [102, 0.87895],
 +
        [103, 0.87895],
 +
        [104, 0.87895],
 +
        [105, 0.87895],
 +
        [106, 0.87895],
 +
        [107, 0.87895],
 +
        [108, 0.87895],
 +
        [109, 0.87895],
 +
        [110, 0.87895],
 +
        [111, 0.87895],
 +
        [112, 0.87895],
 +
        [113, 0.87895],
 +
        [114, 0.87895],
 +
        [115, 0.87895],
 +
        [116, 0.87895],
 +
        [117, 0.87895],
 +
        [118, 0.87895],
 +
        [119, 0.87895],
 +
        [120, 0.87895],
 +
        [121, 0.87895],
 +
        [122, 0.87895],
 +
        [123, 0.87895],
 +
        [124, 0.87895],
 +
        [125, 0.87895],
 +
        [126, 0.87895],
 +
        [127, 0.87895],
 +
        [128, 0.87895],
 +
        [129, 0.87895],
 +
        [130, 0.87895],
 +
        [131, 0.87895],
 +
        [132, 0.87895],
 +
        [133, 0.87895],
 +
        [134, 0.87895],
 +
        [135, 0.87895],
 +
        [136, 0.87895],
 +
        [137, 0.87895],
 +
        [138, 0.87895],
 +
        [139, 0.87895],
 +
        [140, 0.87895],
 +
        [141, 0.87895],
 +
        [142, 0.87895],
 +
        [143, 0.87895],
 +
        [144, 0.87895],
 +
        [145, 0.87895],
 +
        [146, 0.87895],
 +
        [147, 0.87895],
 +
        [148, 0.87895],
 +
        [149, 0.87895],
 +
        [150, 0.87895],
 +
        [151, 0.87895],
 +
        [152, 0.87895],
 +
        [153, 0.87895],
 +
        [154, 0.87895],
 +
        [155, 0.87895],
 +
        [156, 0.87895],
 +
        [157, 0.87895],
 +
        [158, 0.87895],
 +
        [159, 0.87895],
 +
        [160, 0.87895],
 +
        [161, 0.87895],
 +
        [162, 0.87895],
 +
        [163, 0.87895],
 +
        [164, 0.87895],
 +
        [165, 0.87895],
 +
        [166, 0.87895],
 +
        [167, 0.87895],
 +
        [168, 0.87895],
 +
        [169, 0.87895],
 +
        [170, 0.87895],
 +
        [171, 0.87895],
 +
        [172, 0.87895],
 +
        [173, 0.87895],
 +
        [174, 0.87895],
 +
        [175, 0.87895],
 +
        [176, 0.87895],
 +
        [177, 0.87895],
 +
        [178, 0.87895],
 +
        [179, 0.87895],
 +
        [180, 0.87895],
 +
        [181, 0.87895],
 +
        [182, 0.87895],
 +
        [183, 0.87895],
 +
        [184, 0.87895],
 +
        [185, 0.87895],
 +
        [186, 0.87895],
 +
        [187, 0.87895],
 +
        [188, 0.87895],
 +
        [189, 0.87895],
 +
        [190, 0.87895],
 +
        [191, 0.87895],
 +
        [192, 0.87895],
 +
        [193, 0.87895],
 +
        [194, 0.87895],
 +
        [195, 0.87895],
 +
        [196, 0.87895],
 +
        [197, 0.87895],
 +
        [198, 0.87895],
 +
        [199, 0.87895],
 +
        [200, 0.87895],
 +
        [201, 0.87895],
 +
        [202, 0.87895],
 +
        [203, 0.87895],
 +
        [204, 0.87895],
 +
        [205, 0.87895],
 +
        [206, 0.87895],
 +
        [207, 0.87895],
 +
        [208, 0.87895],
 +
        [209, 0.87895],
 +
        [210, 0.87895],
 +
        [211, 0.87895],
 +
        [212, 0.87895],
 +
        [213, 0.87895],
 +
        [214, 0.87895],
 +
        [215, 0.87895],
 +
        [216, 0.87895],
 +
        [217, 0.87895],
 +
        [218, 0.87895],
 +
        [219, 0.87895],
 +
        [220, 0.87895],
 +
        [221, 0.87895],
 +
        [222, 0.87895],
 +
        [223, 0.87895],
 +
        [224, 0.87895],
 +
        [225, 0.87895],
 +
        [226, 0.87895],
 +
        [227, 0.87895],
 +
        [228, 0.87895],
 +
        [229, 0.87895],
 +
        [230, 0.87895],
 +
        [231, 0.87895],
 +
        [232, 0.87895],
 +
        [233, 0.87895],
 +
        [234, 0.87895],
 +
        [235, 0.87895],
 +
        [236, 0.87895],
 +
        [237, 0.87895],
 +
        [238, 0.87895],
 +
        [239, 0.87895],
 +
        [240, 0.87895]
 +
    ]);
 +
 +
    var options = {
 +
        title: 'Expression of mRNA of chromoprotein (mCP)',
 +
        colors: ['#833789'],
 +
        lineWidth: 5,
 +
 +
        hAxis: {
 +
            gridlines: {color: '#F1F1F1'},
 +
            title: 'Time (minutes)',
 +
            textStyle: {
 +
                fontSize: 16,
 +
                //fontName: 'Arial',
 +
                bold: false,
 +
                italic: false
 +
            },
 +
            titleTextStyle: {
 +
                fontSize: 17,
 +
                bold: false,
 +
                italic: false
 +
 +
            }
 +
        },
 +
        vAxis: {
 +
            gridlines: {color: '#F1F1F1'},
 +
            title: 'mCP per cell (nM)',
 +
            textStyle: {
 +
                fontSize: 16,
 +
                //fontName: 'Arial',
 +
                bold: false,
 +
                italic: false
 +
            },
 +
            titleTextStyle: {
 +
                fontSize: 17,
 +
                bold: false,
 +
                italic: false
 +
 +
            }
 +
        },
 +
        legend: {
 +
            position: "none"
 +
        },
 +
        titleTextStyle: {
 +
            fontSize: 19
 +
        },
 +
 +
 +
 +
 +
    };
 +
 +
    var chart = new google.visualization.LineChart(document.getElementById('chart_div'));
 +
 +
    chart.draw(data, options);
 +
 +
}
 +
</script>
 +
<div id="chart_div" style="width: 900px; height: 400px;"></div>
 +
                            <figcaption>
 +
                                <b>Figure 4.</b> The graph shows the steady&ndash;state concentration/expression of mCP over time with respect to the mRNA degradation rates of chromoprotein. It is clear from the graph that presence of 0.87 nM of mCP is important in the cell in order to achieve the steady&ndash;state concentration of mCP.
 +
                            </figcaption>
 +
                        </figure>
 +
 +
                        <br><br>
 +
                        <figure>
 +
                        <!--Demogrpahics-->
 +
<script type="text/javascript">
 +
google.charts.load('current', {
 +
    packages: ['corechart', 'line']
 +
});
 +
google.charts.setOnLoadCallback(drawBasic);
 +
 +
function drawBasic() {
 +
 +
 +
    var data = new google.visualization.DataTable();
 +
    data.addColumn('number', 'Time (minutes)');
 +
    data.addColumn('number', 'CP per cell (nM)');
 +
 +
    data.addRows([
 +
        [0, 0],
 +
        [1, 8.1271e-05],
 +
        [2, 0.0002882],
 +
        [3, 0.00057603],
 +
        [4, 0.00091153],
 +
        [5, 0.0012704],
 +
        [6, 0.0016352],
 +
        [7, 0.0019937],
 +
        [8, 0.0023375],
 +
        [9, 0.0026615],
 +
        [10, 0.0029626],
 +
        [11, 0.0032392],
 +
        [12, 0.003491],
 +
        [13, 0.0037186],
 +
        [14, 0.0039229],
 +
        [15, 0.0041053],
 +
        [16, 0.0042674],
 +
        [17, 0.0044108],
 +
        [18, 0.0045372],
 +
        [19, 0.0046483],
 +
        [20, 0.0047457],
 +
        [21, 0.0048308],
 +
        [22, 0.004905],
 +
        [23, 0.0049695],
 +
        [24, 0.0050255],
 +
        [25, 0.0050741],
 +
        [26, 0.0051161],
 +
        [27, 0.0051524],
 +
        [28, 0.0051837],
 +
        [29, 0.0052107],
 +
        [30, 0.0052339],
 +
        [31, 0.0052539],
 +
        [32, 0.005271],
 +
        [33, 0.0052857],
 +
        [34, 0.0052983],
 +
        [35, 0.0053091],
 +
        [36, 0.0053183],
 +
        [37, 0.0053262],
 +
        [38, 0.0053329],
 +
        [39, 0.0053386],
 +
        [40, 0.0053435],
 +
        [41, 0.0053477],
 +
        [42, 0.0053513],
 +
        [43, 0.0053543],
 +
        [44, 0.0053569],
 +
        [45, 0.0053591],
 +
        [46, 0.0053609],
 +
        [47, 0.0053625],
 +
        [48, 0.0053639],
 +
        [49, 0.005365],
 +
        [50, 0.005366],
 +
        [51, 0.0053668],
 +
        [52, 0.0053675],
 +
        [53, 0.0053681],
 +
        [54, 0.0053686],
 +
        [55, 0.005369],
 +
        [56, 0.0053694],
 +
        [57, 0.0053697],
 +
        [58, 0.0053699],
 +
        [59, 0.0053702],
 +
        [60, 0.0053703],
 +
        [61, 0.0053705],
 +
        [62, 0.0053706],
 +
        [63, 0.0053707],
 +
        [64, 0.0053708],
 +
        [65, 0.0053709],
 +
        [66, 0.005371],
 +
        [67, 0.005371],
 +
        [68, 0.0053711],
 +
        [69, 0.0053711],
 +
        [70, 0.0053712],
 +
        [71, 0.0053712],
 +
        [72, 0.0053712],
 +
        [73, 0.0053712],
 +
        [74, 0.0053713],
 +
        [75, 0.0053713],
 +
        [76, 0.0053713],
 +
        [77, 0.0053713],
 +
        [78, 0.0053713],
 +
        [79, 0.0053713],
 +
        [80, 0.0053713],
 +
        [81, 0.0053713],
 +
        [82, 0.0053713],
 +
        [83, 0.0053713],
 +
        [84, 0.0053713],
 +
        [85, 0.0053713],
 +
        [86, 0.0053713],
 +
        [87, 0.0053713],
 +
        [88, 0.0053713],
 +
        [89, 0.0053713],
 +
        [90, 0.0053713],
 +
        [91, 0.0053713],
 +
        [92, 0.0053713],
 +
        [93, 0.0053713],
 +
        [94, 0.0053713],
 +
        [95, 0.0053713],
 +
        [96, 0.0053713],
 +
        [97, 0.0053713],
 +
        [98, 0.0053713],
 +
        [99, 0.0053713],
 +
        [100, 0.0053713],
 +
        [101, 0.0053713],
 +
        [102, 0.0053713],
 +
        [103, 0.0053713],
 +
        [104, 0.0053713],
 +
        [105, 0.0053713],
 +
        [106, 0.0053713],
 +
        [107, 0.0053713],
 +
        [108, 0.0053713],
 +
        [109, 0.0053713],
 +
        [110, 0.0053713],
 +
        [111, 0.0053713],
 +
        [112, 0.0053713],
 +
        [113, 0.0053713],
 +
        [114, 0.0053713],
 +
        [115, 0.0053713],
 +
        [116, 0.0053713],
 +
        [117, 0.0053713],
 +
        [118, 0.0053713],
 +
        [119, 0.0053713],
 +
        [120, 0.0053713],
 +
        [121, 0.0053713],
 +
        [122, 0.0053713],
 +
        [123, 0.0053713],
 +
        [124, 0.0053713],
 +
        [125, 0.0053713],
 +
        [126, 0.0053713],
 +
        [127, 0.0053713],
 +
        [128, 0.0053713],
 +
        [129, 0.0053713],
 +
        [130, 0.0053713],
 +
        [131, 0.0053713],
 +
        [132, 0.0053713],
 +
        [133, 0.0053713],
 +
        [134, 0.0053713],
 +
        [135, 0.0053713],
 +
        [136, 0.0053713],
 +
        [137, 0.0053713],
 +
        [138, 0.0053713],
 +
        [139, 0.0053713],
 +
        [140, 0.0053713],
 +
        [141, 0.0053713],
 +
        [142, 0.0053713],
 +
        [143, 0.0053713],
 +
        [144, 0.0053713],
 +
        [145, 0.0053713],
 +
        [146, 0.0053713],
 +
        [147, 0.0053713],
 +
        [148, 0.0053713],
 +
        [149, 0.0053713],
 +
        [150, 0.0053713],
 +
        [151, 0.0053713],
 +
        [152, 0.0053713],
 +
        [153, 0.0053713],
 +
        [154, 0.0053713],
 +
        [155, 0.0053713],
 +
        [156, 0.0053713],
 +
        [157, 0.0053713],
 +
        [158, 0.0053713],
 +
        [159, 0.0053713],
 +
        [160, 0.0053713],
 +
        [161, 0.0053713],
 +
        [162, 0.0053713],
 +
        [163, 0.0053713],
 +
        [164, 0.0053713],
 +
        [165, 0.0053713],
 +
        [166, 0.0053713],
 +
        [167, 0.0053713],
 +
        [168, 0.0053713],
 +
        [169, 0.0053713],
 +
        [170, 0.0053713],
 +
        [171, 0.0053713],
 +
        [172, 0.0053713],
 +
        [173, 0.0053713],
 +
        [174, 0.0053713],
 +
        [175, 0.0053713],
 +
        [176, 0.0053713],
 +
        [177, 0.0053713],
 +
        [178, 0.0053713],
 +
        [179, 0.0053713],
 +
        [180, 0.0053713],
 +
        [181, 0.0053713],
 +
        [182, 0.0053713],
 +
        [183, 0.0053713],
 +
        [184, 0.0053713],
 +
        [185, 0.0053713],
 +
        [186, 0.0053713],
 +
        [187, 0.0053713],
 +
        [188, 0.0053713],
 +
        [189, 0.0053713],
 +
        [190, 0.0053713],
 +
        [191, 0.0053713],
 +
        [192, 0.0053713],
 +
        [193, 0.0053713],
 +
        [194, 0.0053713],
 +
        [195, 0.0053713],
 +
        [196, 0.0053713],
 +
        [197, 0.0053713],
 +
        [198, 0.0053713],
 +
        [199, 0.0053713],
 +
        [200, 0.0053713],
 +
        [201, 0.0053713],
 +
        [202, 0.0053713],
 +
        [203, 0.0053713],
 +
        [204, 0.0053713],
 +
        [205, 0.0053713],
 +
        [206, 0.0053713],
 +
        [207, 0.0053713],
 +
        [208, 0.0053713],
 +
        [209, 0.0053713],
 +
        [210, 0.0053713],
 +
        [211, 0.0053713],
 +
        [212, 0.0053713],
 +
        [213, 0.0053713],
 +
        [214, 0.0053713],
 +
        [215, 0.0053713],
 +
        [216, 0.0053713],
 +
        [217, 0.0053713],
 +
        [218, 0.0053713],
 +
        [219, 0.0053713],
 +
        [220, 0.0053713],
 +
        [221, 0.0053713],
 +
        [222, 0.0053713],
 +
        [223, 0.0053713],
 +
        [224, 0.0053713],
 +
        [225, 0.0053713],
 +
        [226, 0.0053713],
 +
        [227, 0.0053713],
 +
        [228, 0.0053713],
 +
        [229, 0.0053713],
 +
        [230, 0.0053713],
 +
        [231, 0.0053713],
 +
        [232, 0.0053713],
 +
        [233, 0.0053713],
 +
        [234, 0.0053713],
 +
        [235, 0.0053713],
 +
        [236, 0.0053713],
 +
        [237, 0.0053713],
 +
        [238, 0.0053713],
 +
        [239, 0.0053713],
 +
        [240, 0.0053713]
 +
    ]);
 +
 +
    var options = {
 +
        title: 'Expression of chromoprotein (CP)',
 +
        colors: ['#248689'],
 +
        lineWidth: 5,
 +
 +
        hAxis: {
 +
            gridlines: {color: '#F1F1F1'},
 +
            title: 'Time (minutes)',
 +
            textStyle: {
 +
                fontSize: 16,
 +
                //fontName: 'Arial',
 +
                bold: false,
 +
                italic: false
 +
            },
 +
            titleTextStyle: {
 +
                fontSize: 17,
 +
                bold: false,
 +
                italic: false
 +
 +
            }
 +
        },
 +
        vAxis: {
 +
            gridlines: {color: '#F1F1F1'},
 +
            title: 'CP per cell (nM)',
 +
            textStyle: {
 +
                fontSize: 16,
 +
                //fontName: 'Arial',
 +
                bold: false,
 +
                italic: false
 +
            },
 +
            titleTextStyle: {
 +
                fontSize: 17,
 +
                bold: false,
 +
                italic: false
 +
 +
            }
 +
        },
 +
        legend: {
 +
            position: "none"
 +
        },
 +
        titleTextStyle: {
 +
            fontSize: 19
 +
        },
 +
 +
 +
 +
 +
    };
  
 +
    var chart = new google.visualization.LineChart(document.getElementById('chart2_div'));
  
 +
    chart.draw(data, options);
  
 +
}
 +
</script>
 +
<div id="chart2_div" style="width: 900px; height: 400px;"></div>
 +
                            <figcaption>
 +
                                <b>Figure 5.</b> The graph shows the relationship between the steady&ndash;state concentration over time of CP with respect to its degradation rate. It can be seen from the graph that 0.005 nM of CP is required in the cell in order to acquire the steady&ndash;state concentration of CP.
 +
                            </figcaption>
 +
                        </figure>
 +
                        <h3>Matlab Code</h3>
 +
<pre>
 +
    <code class="language-matlab">
 +
%Defining the parameters
 +
ktr    = 0.167;    %mRNA synthesis rate of
 +
                    %Chromoprotein per minute (rate of transcription)
 +
ktrans  = 0.0011;    %Chromoprotein synthesis rate
 +
                    %per minute (rate of translation)
 +
ycp    = 0.18;      %Chromoproteins degradation rate per minute
 +
ym      = 0.19;      %mRNA degradation rate of chromoproteins per minute
 +
Kno    = 0.12;      %Dissociation constant for inducer (NO)
 +
Kdn    = 35;        %Dissociation constant for repressor (NsrR)
 +
pyear  = 1;        %Number of plasmids per cell
 +
n      = 2;        %Cooperativity of inducer (NO)
  
 +
%Assume that the repressor is removed
 +
NsrR    = 0;
 +
NO      = 2;
 +
t      = 0:0.1:240;
  
 +
%Calculate mRNA and Chromoprotein values
 +
K      = ktr*pyear/(1+(NsrR/((1+(NO/Kno)^n)*Kdn)));
 +
mCP    = K/ym-K/ym*exp(-ym*t);
 +
CP      = K*ktrans/(ycp*ym)+K*ktrans/(ym^2-ycp*ym)*exp(-ym*t)...
 +
    -(K*ktrans/(ycp*ym)+K*ktrans/(ym^2-ycp*ym))*exp(-ycp*t);
  
 +
disp('Steady State mCP Concentration:')
 +
disp(mCP(end)) %Steady state mRNA value
 +
disp('Steady State CP Concentration:')
 +
disp(CP(end)) %Steady state chromoprotein value
 +
 +
%Plotting mRNA concentration over time
 +
figure(1)
 +
plot(t, mCP,'r','LineWidth',2)
 +
hold on
 +
grid on
 +
xlabel('Time (minutes)')
 +
ylabel('mCP per cell (nM)')
 +
title('Expression of mRNA of chromoprotein (mCP)')
  
<div class="column full_size">
+
%Plotting the chromoprotein concentrations over time
<h2> Modeling</h2>
+
figure(2)
<p>Mathematical models and computer simulations provide a great way to describe the function and operation of BioBrick Parts and Devices. Synthetic Biology is an engineering discipline, and part of engineering is simulation and modeling to determine the behavior of your design before you build it. Designing and simulating can be iterated many times in a computer before moving to the lab. This award is for teams who build a model of their system and use it to inform system design or simulate expected behavior in conjunction with experiments in the wetlab.</p>
+
plot(t,CP,'r', 'LineWidth',2)
 +
grid on
 +
hold on
 +
xlabel('Time (minutes)')
 +
ylabel('CP per cell (nM)')
 +
title('Expression of chromoprotein (CP)')
  
 +
%Save this as nox_sensing_system_simulation.m file
 +
    </code>
 +
</pre>
 +
                        <h3>Conclusion</h3>
 +
                        <p>The system is considered as stable, because when NsrR detaches itself from pYeaR, the promoter then initiate the production of mCP and CP, until a steady&ndash;state concentration is reached. This steady&ndash;state concentration is the maximum production rate of mCP and CP with respect to their degradation rates. However, it is clear from the graphs that the degradation rates are only related to the maximum production rate/ number of molecules in the cell; it is not responsible for altering the time span required to reach the steady&ndash;state concentration. Furthermore, this simplified solution for NOx&ndash;sensing system demonstrates the migration of system towards steady&ndash;state under assumed parameters and prescribed conditions. Therefore, this mathematical model can be considered as an appropriate model for determining the behavior of NOx&ndash;sensing system. </p>
 +
                       
 +
                        <h2>References</h2>
 +
                        <ol>
 +
                        <li><p>Li, J., Zhao, Z., Kazakov, A., Chaos, M., Dryer, F.L. and Scire, J.J. (2007) 'A comprehensive kinetic mechanism for CO, CH2O, and CH3OH combustion', International Journal of Chemical Kinetics, 39(3), pp. 109&ndash;136. doi: 10.1002/kin.20218.</p></li>
 +
                        <li><p>Chen, H., Shiroguchi, K., Ge, H. and Xie, X.S. (2015) 'Genome-wide study of mRNA degradation and transcript elongation in Escherichia coli', Molecular Systems Biology, 11(1), pp. 781&ndash;781. doi: 10.15252/msb.20145794.</p></li>
 +
                        <li><p>Leduc, J., Thorsteinsson, M.V., Gaal, T. and Roberts, G.P. (2001) 'Mapping CooA.RNA polymerase interactions: Identification of activating regions 2 and 3 in CooA, the co&ndash;sensing transcriptional activator', Journal of Biological Chemistry, 276(43), pp. 39968&ndash;39973. doi: 10.1074/jbc.m105758200.</p></li>
 +
                        <li><p>Benabbas, A., Karunakaran, V., Youn, H., Poulos, T.L. and Champion, P.M. (2012) 'Effect of DNA binding on Geminate CO recombination Kinetics in CO-sensing transcription factor CooA', Journal of Biological Chemistry, 287(26), pp. 21729&ndash;21740. doi: 10.1074/jbc.m112.345090.</p></li>
 +
                        <li><p>Kuchinskas, M., Li, H., Conrad, M., Roberts, G. and Poulos, T.L. (2006) 'The role of the DNA-Binding domains in CooA activation', Biochemistry, 45(23), pp. 7148&ndash;27153. doi: 10.1021/bi052609o.</p></li>
 +
                        <li><p>Bergmann, John E, and Harvey F Lodish. 'A kinetic model of protein synthesis. Application to hemoglobin synthesis and translational control.' The Journal of biological chemistry 254.23 (1979) : 11927&ndash;37.</p></li>
 +
                        <li><p>Nick Csicsery, and Ricky O'Laughlin. 'A Mathematical Model of a Synthetically Constructed Genetic Toggle Switch' (2015) </p></li>
 +
                        </ol>
  
<h5> Inspiration </h5>
 
<p>
 
Here are a few examples from previous teams:
 
</p>
 
<ul>
 
<li><a href="https://2014.igem.org/Team:ETH_Zurich/modeling/overview">ETH Zurich 2014</a></li>
 
<li><a href="https://2014.igem.org/Team:Waterloo/Math_Book">Waterloo 2014</a></li>
 
</ul>
 
  
 +
                    </div>
 +
                    <!--End of centered-content-wrapper-->
 +
                </div>
 +
                <!--End of centered-content-wrapper-container-->
 +
               
 +
                <div class="all-in-one-footer-container">
 +
    <div class="all-in-one-footer-address">
 +
        <h3>Reach Us</h3>
 +
        <h4><img id="location" src="https://static.igem.org/mediawiki/2016/7/7c/T--Peshawar--location-compressor.png">  Institute of Integrative Biosciences,
 +
        <br>CECOS University of IT &amp; Emerging Sciences,
 +
        <br>Phase 6, F-5, Hayatabad, Peshawar, Pakistan
 +
        </h4>
 +
        <h4 id="h4-email"><img id="email" src="https://static.igem.org/mediawiki/2016/0/04/T--Peshawar--email.png">
 +
        <a href="mailto:iib@cecos.edu.pk">iib@cecos.edu.pk</a>
 +
        </h4>
 +
        <h4 id="h4-phone"><img id="phone" src="https://static.igem.org/mediawiki/2016/c/c6/T--Peshawar--phone-compressor.png">  +92 91 5860291</h4>
 +
    </div>
  
 +
    <div class="all-in-one-sponsors">
 +
        <h3>Our Partners</h3>
 +
        <br>
 +
        <div class="sponsors-icons govt-kpk">
 +
            <a href="http://kp.gov.pk/" target="_blank">
 +
                <img src="https://static.igem.org/mediawiki/2016/a/a1/T--Peshawar--govt-of-kpk-logo-small-gs-compressed.png">
 +
            </a>
 +
        </div>
 +
        <div class="sponsors-icons dost">
 +
            <a href="http://www.dost.gov.pk/" target="_blank">
 +
                <img src="https://static.igem.org/mediawiki/2016/e/e2/T--peshawar--dost-logo-small-gs-compresed-compressor.png">
 +
            </a>
 +
        </div>
 +
        <div class="sponsors-icons cecos">
 +
            <a href="http://cecos.edu.pk/" target="_blank">
 +
                <img src="https://static.igem.org/mediawiki/2016/0/03/T--Peshawar--cecos-logo-small-gs-compressed-compressor.png">
 +
            </a>
 +
        </div>
 +
        <div class="sponsors-icons iib">
 +
            <a href="http://iib.cecos.edu.pk/" target="_blank">
 +
                <img src="https://static.igem.org/mediawiki/2016/6/60/T--Peshawar--iib-logo-small-gs-compressed-compressor.png">
 +
            </a>
 +
        </div>
 +
    </div>
 +
    <div class="all-in-one-footer-social">
 +
        <h3>Connect with us</h3>
 +
        <br>
 +
        <article id="social-page-footer">
 +
    <!--FACEBOOK-->
 +
    <a href="https://www.facebook.com/iGEMPeshawar2016/" target="_blank">
 +
        <svg class="social-page-footer-svg facebook" width="61" height="60" viewBox="0 0 61 60" xmlns="http://www.w3.org/2000/svg">
 +
            <title>
 +
                T--Peshawar--facebook-round-v2
 +
            </title>
 +
            <g id="Page-1" fill="none" fill-rule="evenodd">
 +
                <g id="T--Peshawar--facebook-round-v2" fill="#808184">
 +
                    <path d="M30.594 1.465c-15.735 0-28.536 12.8-28.536 28.537 0 15.734 12.8 28.536 28.536 28.536S59.13 45.736 59.13 30.002c0-15.736-12.8-28.537-28.536-28.537m0 58.536c-16.542 0-30-13.455-30-29.997 0-16.543 13.458-30 30-30s30 13.457 30 30c0 16.542-13.458 30-30 30" id="facebook-circle" />
 +
                    <path d="M22.91 25.905h3.72V22.29c0-1.597.042-4.057 1.2-5.58 1.22-1.614 2.895-2.71 5.777-2.71 4.696 0 6.672.67 6.672.67l-.933 5.514s-1.55-.45-2.997-.45c-1.447 0-2.743.52-2.743 1.966v4.205h5.936l-.414 5.386h-5.524V50H26.63V31.29h-3.72v-5.385z" id="facebook" />
 +
                </g>
 +
            </g>
 +
        </svg>
 +
    </a>
 +
    <!--TWITTER-->
 +
    <a href="https://twitter.com/igem_peshawar" target="_blank">
 +
        <svg class="social-page-footer-svg twitter" width="60" height="60" viewBox="0 0 60 60" xmlns="http://www.w3.org/2000/svg">
 +
            <title>
 +
                T--Peshawar--twitter-round-v2
 +
            </title>
 +
            <g id="Page-1" fill="none" fill-rule="evenodd">
 +
                <g id="T--Peshawar--twitter-round-v2" fill="#808184">
 +
                    <path d="M30 1.463C14.266 1.463 1.464 14.265 1.464 30S14.264 58.537 30 58.537 58.538 45.735 58.538 30 45.738 1.463 30 1.463M30 60C13.46 60 0 46.542 0 30S13.46 0 30 0c16.542 0 30 13.458 30 30S46.542 60 30 60" id="twitter-circle" />
 +
                    <path d="M49.642 21.783c-1.325.587-2.75.984-4.242 1.163 1.527-.914 2.698-2.36 3.247-4.085-1.428.848-3.006 1.46-4.692 1.79-1.343-1.433-3.265-2.33-5.39-2.33-4.077 0-7.384 3.308-7.384 7.385 0 .58.065 1.142.19 1.683-6.134-.31-11.576-3.247-15.222-7.718-.636 1.093-.998 2.362-.998 3.715 0 2.56 1.303 4.82 3.283 6.148-1.208-.038-2.35-.373-3.347-.923v.09c0 3.58 2.547 6.566 5.927 7.242-.62.174-1.27.26-1.946.26-.476 0-.94-.044-1.39-.13.94 2.933 3.666 5.07 6.9 5.126-2.53 1.982-5.715 3.163-9.173 3.163-.596 0-1.185-.032-1.762-.1 3.27 2.093 7.15 3.316 11.32 3.316 13.59 0 21.015-11.25 21.015-21.013 0-.32-.006-.64-.02-.955 1.444-1.043 2.7-2.344 3.685-3.824" id="twitter" />
 +
                </g>
 +
            </g>
 +
        </svg>
 +
    </a>
 +
    <!--INSTAGRAM-->
 +
    <a href="https://www.instagram.com/igempeshawar/" target="_blank">
 +
        <svg class="social-page-footer-svg instagram" width="61" height="61" viewBox="0 0 61 61" xmlns="http://www.w3.org/2000/svg">
 +
            <title>
 +
                T--Peshawar--instagram-round-v2
 +
            </title>
 +
            <g id="Page-1" fill="none" fill-rule="evenodd">
 +
                <g id="T--Peshawar--instagram-round-v2" fill="#808184">
 +
                    <path d="M30.802 1.64c-15.736 0-28.537 12.802-28.537 28.537s12.8 28.537 28.537 28.537c15.735 0 28.536-12.802 28.536-28.537S46.538 1.64 30.802 1.64m0 58.537c-16.542 0-30-13.458-30-30s13.458-30 30-30 30 13.458 30 30-13.458 30-30 30" id="instagram-circle" />
 +
                    <path d="M45.3 42.212c0 1.902-1.548 3.446-3.45 3.446H19.763c-1.903 0-3.45-1.544-3.45-3.446v-14.73h5.375c-.464 1.144-.727 2.39-.727 3.695 0 5.426 4.417 9.837 9.845 9.837 5.428 0 9.843-4.41 9.843-9.837 0-1.306-.264-2.55-.725-3.694H45.3v14.73zm-19.63-14.73c1.15-1.593 3.022-2.64 5.137-2.64 2.115 0 3.985 1.047 5.136 2.64.75 1.042 1.203 2.314 1.203 3.695 0 3.492-2.845 6.33-6.34 6.33-3.494 0-6.336-2.838-6.336-6.33 0-1.38.45-2.653 1.2-3.694zm18.174-10.145l.795-.004v6.097l-6.083.022-.022-6.098 5.308-.017zM41.85 13.19H19.763c-3.835 0-6.957 3.118-6.957 6.952v22.07c0 3.834 3.122 6.952 6.957 6.952H41.85c3.834 0 6.957-3.118 6.957-6.952v-22.07c0-3.834-3.123-6.952-6.958-6.952z" id="instagram" />
 +
                </g>
 +
            </g>
 +
        </svg>
 +
    </a>
 +
    <!--BLOG-->
 +
    <a href="http://www.synbiokp.com/blog/" target="_blank">
 +
        <svg class="social-page-footer-svg blog" width="60" height="60" viewBox="0 0 60 60" xmlns="http://www.w3.org/2000/svg">
 +
            <title>
 +
                T--Peshawar--blog-round-v2
 +
            </title>
 +
            <g id="Page-1" fill="none" fill-rule="evenodd">
 +
                <g id="T--Peshawar--blog-round-v2" fill="#808184">
 +
                    <path d="M30 1.463C14.265 1.463 1.463 14.263 1.463 30c0 15.735 12.802 28.536 28.537 28.536S58.537 45.736 58.537 30 45.735 1.463 30 1.463M0 30C0 13.458 13.458 0 30 0s30 13.458 30 30-13.458 30-30 30S0 46.542 0 30z" id="blog-circle" />
 +
                    <path d="M34.95 19.916c7.654 0 14 6.357 14 14.025 0 7.64-6.375 14.03-14 14.03h-7.942C19.393 47.968 13 41.57 13 33.94V15.22c0-1.78 1.47-3.255 3.25-3.255 1.768 0 3.247 1.482 3.247 3.256 0 6.24.01 12.48.01 18.72 0 4.11 3.4 7.52 7.5 7.52h7.943c4.084 0 7.504-3.424 7.504-7.514 0-4.108-3.403-7.516-7.504-7.516h-8.568c-1.777 0-3.248-1.475-3.248-3.255 0-1.773 1.477-3.254 3.248-3.254h8.568z" id="blog" />
 +
                </g>
 +
            </g>
 +
        </svg>
 +
    </a>
 +
</article>
 +
        <div class="wiki-design">
 +
            <!-- <h4>Wiki designed by <a href="https://2016.igem.org/Team:Peshawar/Team">Adnan Niazi</a></h4> -->
 +
            <br>
 +
            <button><a href="https://2016.igem.org/Team:Peshawar/Wiki">Learn iGEM Wiki design</a></button>
 +
        </div>
 +
    </div>
 
</div>
 
</div>
 +
<!--End of all-in-one-footer-container-->
 +
            </section>
 +
            <!--End of inner-container-->
 +
        </div>
 +
        <!--End of on-canvas container for all the content-->
 +
    </div>
 +
    <!--End of outermost box-->
 +
<script type="text/javascript" src="https://2016.igem.org/Template:Peshawar/Javascript?action=raw&ctype=text/javascript"></script>
 +
    <!-- Start of StatCounter Code for Default Guide -->
 +
<script type="text/javascript">
 +
var sc_project=11100754;
 +
var sc_invisible=1;
 +
var sc_security="135c4f8d";
 +
var scJsHost = (("https:" == document.location.protocol) ?
 +
"https://secure." : "http://www.");
 +
document.write("<sc"+"ript type='text/javascript' src='" +
 +
scJsHost+
 +
"statcounter.com/counter/counter.js'></"+"script>");
 +
</script>
 +
<noscript><div class="statcounter"><a title="shopify
 +
analytics ecommerce" href="http://statcounter.com/shopify/"
 +
target="_blank"><img class="statcounter"
 +
src="//c.statcounter.com/11100754/0/135c4f8d/1/"
 +
alt="shopify analytics ecommerce"></a></div></noscript>
 +
<!-- End of StatCounter Code for Default Guide -->
 +
</body>
  
 
</html>
 
</html>

Latest revision as of 22:24, 19 October 2016

Modeling - iGEM Peshawar 2016

Modeling and Simulation

Why model?

Before building a part or a device in lab, it is probably a good idea to simulate it in a simulation package — such as SimBiology, Cello, or GenoCAD — to check if the part or device one intends to make would work the way one wants it to. Although simulation can't solve all the problems that one might encounter during the actual physical design, it can certainly reduce them, and help in getting a better understanding of how a design behaves under different parameters and conditions.

Our project has two biological ystems: one for detection of corbon monoxide, and another for detection for oxides of nitrogen. Below we show you how we modelled and simulated both these systems.

Modeling of the CO-sensing system

In our CO-sensing system, the CO binds to CooA (an activator). It activates the activator protein, and so the activated activator protein (CooA)2CO binds to promoter and initiates transcription.

The reactions which take place can be given as:

$$\textrm{[CO]}+\textrm{[CooA]} \longleftrightarrow \textrm{[(CooA)2CO]} \label{R1} \tag{R1}$$

The equation above represents the first reaction that takes place in our system in which the CO present in the environment is taken up into the cell where it binds to the CooA protein. The CooA protein is constantly expressed in our system using a constitutive promoter. The coefficient 2 represents the number of CO which bind to CooA; this number was taken from the literature.

To develop a dissociation constant of this equation, we can write it as:

$$KI=\frac{\textrm{[CO].[CooA]}}{\textrm{[(CooA)2CO]}} \label{1} \tag{1}$$

Note that the concentrations of $\textrm{[CO]}$, $\textrm{[CooA]}$, and $\textrm{[(CooA)2CO]}$ are at equilibrium before the start of the reaction R1.

The second reaction that takes place in our system is the binding of the activated activator protein complex of $\textrm{[(CooA)2CO]}$ to the promoter Pr.

It can be described as:

$$\textrm{[(CooA)2CO]} \longleftrightarrow \textrm{[(CooA)2CO.Pr]} \label{R2} \tag{R2}$$

Here, $\textrm{[(CooA)2CO.Pr]}$ represents the activated protein $\textrm{[(CooA)2CO]}$ bound to the promoter.

To develop a dissociation constant of this equation, we can write as:

$$KP=\frac{\textrm{[(CooA)2CO]}}{\textrm{[(CooA)2CO.Pr]}} \label{2} \tag{2}$$

Now, in order to get an equation which describes the change in the concentration of mRNA with time we develop the following differential equation:

$$\frac{d(mRNA)}{dt}=k_{on1}.P_{exp}-\gamma_1.[mRNA] \label{3} \tag{3}$$

where, $k_{on1}$ is represents the formation rate constant of the mRNA, and $\gamma_1$ represents the degradation rate constant of mRNA.

In Eq. (\ref{3}) it can be seen that the change in the $\textrm{[mRNA]}$ with time is equal to the difference of the formation rate of mRNA and its degradation rate.

The initial conditions for the Eq. (\ref{3}) are assumed to be:

$$\textrm{[mRNA(t=0)] = 0 molar}$$

This means that at the time $t=0$ sec (initial time), there are no molecules of mRNA encoding for our desired chromoprotein in the cell.

In Eq. (\ref{3}), $P_{exp}$ represents the probability that the promoter will be activated by the activated activator protein complex $\textrm{[(CooA)2CO]}$. It can be calculated by the following equation:

$$P_{exp} = \frac{\textrm{number of situations in which the promoter will be activated}}{\textrm{total number of situations}} = \frac{w}{z} \label{4} \tag{4}$$

In our modeling of the CO-sensing system, we have assumed that at any moment only one molecule of the activated activator protein complex $\textrm{[(CooA)2CO]}$ can bind to the promoter, and that the \ref{R1} reaches equilibrium before the transcription of the system begins.

To develop our equations further, we say that the number of situations in which a promoter will be activated is equal to:

$$ w=\textrm{[(CooA)2CO.Pr]}=\frac{\textrm{[(CooA)2CO]}}{K_p} \label{5} \tag{5}$$

Eq. (\ref{5}) demonstrates the relation between the binding of the activated activator protein complex $\textrm{[(CooA)2CO]}$, showing that it is equal to the amount of $\textrm{[(CooA)2CO]}$ upon the dissociation constant of the binding of $\textrm{[(CooA)2CO]}$ to the promoter.

To simplify this equation, we assume that the number of situations in which the promoter will not be activated by $\textrm{[(CooA)2CO]}$ is equal to $1$.

Note that this assumption is made because the only situation in which transcription will not take place is when there are no molecules bound to the promoter.

Substituting the assumptions in Eq. (\ref{5}), we get:

$$ z=w+1=\frac{\textrm{[(CooA)2CO]}}{K_p}+1 \label{6} \tag{6}$$

where, $z$ represents the total number of situations, $w$ represents the number of situations in which the promoter will be activated, and $1$ represents the number of situations in which the promoter will not be activated.

Sustituting Eq. (\ref{5}) and (\ref{6}) into Eq. (\ref{4}), we get:

$$P_{exp} = \frac{w}{z} = \frac{\frac{\textrm{[(CooA)2CO]}}{K_p}}{\frac{\textrm{[(CooA)2CO]}}{K_p}+1} \label{7} \tag{7}$$

As we want our $P_{exp}$ to be expressed in the terms of $\textrm{[CO]}$ and $\textrm{[CooA]}$ separately instead of $\textrm{[(CooA)2CO]}$, we repeat these above equations by taking into consideration \ref{R1} and its dissociation constant $\textrm{KI}$.

By doing so, we get:

$$P_{exp} = \frac{\frac{\textrm{[CooA]}.\frac{\frac{\textrm{CO}}{K_I}}{1+\frac{\textrm{CO}}{K_I}}}{K_p}} {1+\frac{\textrm{[CooA]}.\frac{\frac{\textrm{CO}}{K_I}}{1+\frac{\textrm{CO}}{K_I}}}{K_p}} \label{8} \tag{8}$$

As the Eq. (\ref{8}) demonstrates, $P_{exp}$ depends only upon the concentration of $\textrm{[CO]}$ and $\textrm{[CooA]}$.

Recall Eq. (\ref{3})

$$\frac{d(\textrm{mRNA})}{dt}=k_{on1}.P_{exp}-\gamma_1.[\textrm{mRNA}]$$

The initial condition for this equation was $\textrm{[mRNA(t=0)] = 0 molar}$.

The analytical solution of this differential equation is equal to:

$$[mRNA]=\frac{k_{on1}.P_{exp}}{\gamma_1}(1-e^{-\gamma_1t})$$

Now that we have got the equation describing the change in concentration of mRNA, we can use the same method to derive an equation which describes the change in the concentration of the chromoprotein in time as following:

$$\frac{d\textrm{[CP]}}{dt} = k_{on2}.\textrm{[mRNA]}-\gamma_2\textrm{[CP]} \label{9} \tag{9}$$

Here, $k_{on2}$ is the formation rate constant of the chromoprotein and $\gamma_2$ is the degradation rate constant of the chromoprotein.


Figure 1: The central theme considered for mathematical modeling of CO system, where the chromoprotein DNA leads to the production of mRNA of chromoprotein and then subsequently of chromoprotein. $K_{on1}$ is the rate of transcription of chromoprotein DNA, while $K_{on2}$ is the translation rate of mRNA of chromoprotein. Also, $\gamma_1$ is the degradation rate of mRNA of cheomoprotein and $\gamma_2$ is the degradation rate of chromoprotein.

Again, the initial conditions for this equation are assumed to be:

$$\textrm{[CP(t=0)] = 0 molar}$$

which states that at the initial time, the amount of chromoprotein in the cell is equal to $0$.

Here on forth, we will repeat the steps from Eq. (\ref{4}) to Eq. (\ref{8}) and arrive to the analytical solution of the differential equation of Eq. (\ref{9}), which is:

$$\textrm{[CP]} = k_{on2}.\frac{k_{on1}.P_{exp}}{\gamma_1.\gamma_2}(1-e^{-\gamma_2t}) + k_{on2}.\frac{k_{on1}.P_{exp}}{\gamma_1(\gamma_1-\gamma_2)}(e^{-\gamma_2t}-e^{-\gamma_1t}) \label{10} \tag{10}$$

Eq. (\ref{10}) is the final equation which describes the amount of $\textrm{[CP]}$ along with the amount of $\textrm{[mRNA]}$.

Assumptions

  1. 1. Reactions \ref{R1} and \ref{R2} reach equilibrium before the transcription begins.
  2. 2. Only one molecule of the complex $\textrm{[(CooA)2CO]}$ can be bound to the promoter at a time.
  3. 3. $\textrm{[mRNA(t=0)] = 0 molar}$
  4. 4. $\textrm{[CO]}$ and $\textrm{[CooA]}$ concentrations inside the cell are constant.

Simulation of the CO-sensing system

Simulation parameters


Table 1: Parameters, their description, and their values used in MATLAB simulation of CO-sensing system.
Parameter Description Value
$\textrm{[CO]}$ Carbon monoxide concentration at equilibrium $\textrm{10}$
$\textrm{[CooA]}$ CooA protein concentration at equilibrium $11$
$\textrm{KI}$ Dissociation constant of the binding of CO to CooA $\textrm{0.011 mM}$
$\textrm{KP}$ Dissociation Constant of the binding of $\textrm{(CO2CooA)}$ to the promoter $\textrm{0.000009}$ $\textrm{min}^{-1}$
Kon1 Formation rate constant of mRNA $\textrm{0.001}$ $\textrm{min}^{-1}$
Kon2 Formation rate constant of the chromoprotein 0.006 $\textrm{min}^{-1}$ / 0.0051 $\textrm{min}^{-1}$
Gama1 ($\gamma$) Degradation rate constant of the mRNA $\textrm{0.00039}$ $\textrm{min}^{-1}$
Gama2 ($\gamma$) Degradation rate constant of the chromoprotein $\textrm{0.001}$ $\textrm{min}^{-1}$

Simulation results



Figure 2: Activity of the different promoters after stimulation by carbon monoxide. The strong promoter pCooF, shown in blue, is yielding higher level of expression than the weak promoter pCooM, shown in violet. The difference between the Strong promoter (pCooF) and the weak promoter (pCooM) is meaningful, and both of these promoters will express enough chromoprotein to be visual. Hence, our goal is achieved.

Conclusion

It is clear from the results that our model will achieve the steady state and give us the amount of chromoprotein that can be visualized.

MATLAB Code


function dx = ODEfun_CooM(t,x)
    dx      = zeros(1,2);
    %x is a vector (x(1),x(2))
    %x(1) = [mRNA]
    %x(2) = [Chromoprotein]
    KI      = 0.011;    %Dissociation constant of the binding of CO to CooA
    KP      = 0.1;      %Dissociation constant of the binding of 
                        %(CO2CooA) to the promoter.
    Kon1    = 0.001;    %Formation rate constant of mRNA
    Kon2    = 0.0051;   %Formation rate constant of the CP
    gama1   = 0.00039;  %Degradation rate constant of mRNA
    gama2   = 0.002;    %Degradation rate constant of the CP
    co      = 10;       %Concentration of CO
    cooa    = 11;       %Concentration of CooA
    w       = (cooa*co/(KI*(1+co/KI)))/KP;
    z       = (cooa*co/(KI*(1+co/KI)))/KP+1;
    P_exp   = w/z;
    dx1     = (Kon1*P_exp)-(gama1*x(1));    %Rate equation of mRNA
    dx2     = (Kon2*x(1))-(gama2*x(2));     %Rate equation of CP
    dx      = [dx1;dx2];
end

%Save this code as ODEfun_CooM.m file


function dx = ODEfun_CooF(t,x)
    dx      = zeros(1,2);
    %x is a vector [x(1),x(2)]
    %x(1) = [mRNA]
    %x(2) = [chromoprotein]
    KI      = 0.011;    %Dissociation constant of the binding of CO to CooA
    KP      = 0.1;      %Dissociation constant of the binding of 
                        %(CO2CooA) to the promoter.
    Kon1    = 0.001;    %Formation rate constant of mRNA
    Kon2    = 0.006;    %Formation rate constant of the CP
    gama1   = 0.00039;  %Degradation rate constant of mRNA
    gama2   = 0.002;    %Degradation rate constant of the CP
    co      = 10;       %Concentration of CO
    cooa    = 11;       %Concentration of CooA
    w       = (cooa*co/(KI*(1+co/KI)))/KP;
    z       = (cooa*co/(KI*(1+co/KI)))/KP+1;
    P_exp   = w/z;
    dx1     = (Kon1*P_exp)-(gama1*x(1));    %Rate equation of mRNA
    dx2     = (Kon2*x(1))-(gama2*x(2));     %Rate equation of CP
    dx      = [dx1;dx2];
end

%Save this code as ODEfun_CooF.m file


% Clear worskpace
clc 
clear all
% Code for calling the previous function:
time_range      = 18000; 
[T1,X1]         = ode45(@ODEfun_CooF,[0 time_range],[0 0]); 
[T2,X2]         = ode45(@ODEfun_CooM,[0 time_range],[0 0]); 
figure 
plot(T1,X1(:,2),'r', T2,X2(:,2),'b', T1,X1(:,1),'g')
legend('Conc. of chromoprotein with pCooF promoter ',...
    'Conc. of chromoprotein with pCooM promoter',...
    'mRNA concentration') 
xlabel('Time [sec]'); 
ylabel('Concentration [molar]'); 
title('Chromoprotein concentration with different promoters (pCooF, pCooM)')

%Save this code as co_sensing_system_simulation.m file

Modeling of the NOx-sensing system

To gain an insight into the stability of NO-sensing circuit, we consider the cooperativity of NO binding with NsrR and the law of mass action. The goal of the modeling is to determine the stable expression of chromoprotein by considering a specified concentration of different parameters of the dynamic system.

Rate of change of chromoprotein expression

The following ordinary differential equations are used for determining the expression levels of chromoproteins (CP) and its mRNA (mCP):

$$\frac{d\textrm{[mCP]}}{dt} = k_{transcription}\textrm{[PyeaR]} - \gamma_m \textrm{[mCP]} \label{11} \tag{11}$$ $$\frac{d\textrm{[CP]}}{dt}=k_{translation}\textrm{[mCP]}-\gamma_p\textrm{[CP]} \label{12} \tag{12}$$
Figure 3: The core idea used in mathematical modeling of NOx system, where the activation of PyeaR promotor leads to the production of mRNA of chromoprotein and then subsequently of chromoprotein. $K_{transcription}$ is the transcription rate of PyeaR promotor, while $K_{translation}$ is the translation rate of mRNA of chromoprotein. Additionally, $\gamma_m$ is the degradation rate of mCP (mRNA of cheomoprotein) and $\gamma_p$ is the degradation rate of CP (chromoprotein).

Binding cooperativity betweeen NO AND NsrR

The cooperativity of $\textrm{NO}$ binding with NsrR can be found using the Hill Equation, where free $\textrm{NO}^*$ combines with the $\textrm{n}$ amount of NsrR, leading to the formation of $\textrm{n[NsrR][NO]}$ complex. The total amount of inducer $NO_T$ concentration can be described by the following equation:

$$\textrm{NO}_T = \textrm{n[NsrR][NO]} + [\textrm{NO}^*]$$

Additionally, $\textrm{NsrR}$ dimerization with respect to unbound $\textrm{NsrR}$ can be represented by the following equation by considering the degree of cooperativity of $\textrm{NO}$ binding:

$$\frac{[\textrm{NO}^*]}{[\textrm{NO}_T]} = \frac{1}{1+\left(\frac{NsrR}{K_{NO}} \right)^n} \label{13} \tag{13}$$

Equation for PyeaR activity

The correlation between PyeaR activity and rate of transcription of mCP (\beta) can represented as:

$$\textrm{[PyeaR]}activity = \frac{\beta}{1+ \left( \frac{\textrm{NO}^*}{K_{d(\textrm{NsrR})}} \right)} \label{14} \tag{14}$$

Solving equations for chromoprotein expressions

Combining the Eqs. (\ref{13}) and (\ref{14}) into equation (\ref{11}), yields:

$$\frac{d[\textrm{mCP}]}{dt} = \frac{k_{transcription}[\textrm{PyeaR}]}{1 + \left(\frac{[\textrm{NsrR}]}{\left(1+\left[\frac{[\textrm{NO}]}{K_{NO}}\right]^n \right)}\right) k_{d(\textrm{NsrR})} } - \gamma_m [\textrm{mCP}]$$

In order to further simplify the system, let

$$K=\frac{k_{transcription}[\textrm{PyeaR}]}{1 + \left(\frac{[\textrm{NsrR}]}{\left(1+\left[\frac{[\textrm{NO}]}{K_{NO}}\right]^n \right)}\right) k_{d(\textrm{NsrR})} }$$

Additionally, the value of $\textrm{K}$ is inversely proportional to the amount of $\textrm{NsrR}$ and directly proportional to the $\textrm{NO}$ concentration in the system. Substituting $\textrm{K}$ into the differential equation of $\textrm{mCP}$ generates the following equation:

$$\frac{d[\textrm{mCP}]}{dt} = K - \gamma_m [\textrm{mCP}]$$

After solving the equation through homogenous and particular solution, the following equation is obtained for steady—state concentration of mCP generated after PyeaR is triggered:

$$[\textrm{mCP}]=\frac{K}{\gamma_m} - \frac{K}{\gamma_m}e^{-\gamma_mt}$$

The steady—state concentration of mCP is then substituted into Eq. (\ref{12}) and is also solved through homogenous and particular solution. The resulting equation gives the steady—state concentration of CP in the desired system according to the specified formation and degradation rates of CP in response to promotor stimulation:

$$[\textrm{CP}]=\frac{K k_{translation}}{\gamma_p \gamma_m} + \frac{K k_{translation}}{{\gamma_m}^2-\gamma_p \gamma_m} e^{-\gamma_mt} - \left(\frac{K k_{translation}}{\gamma_p \gamma_m} + \frac{K k_{translation}}{{\gamma_m}^2-\gamma_p \gamma_m} \right)e^{-\gamma_pt} $$

Simulation of the NOx-sensing system

Simulation parameters

Table 2: Parameters, their description, and their values used in MATLAB simulation of the NOx-sensing system.
Parameter Description Value
$[\textrm{PyeaR}]$ Concentration of PyeaR 1 nM
$k_{transcription}$ Rate of chromoprotein mRNA synthesis 0.167$\textrm{min}^{-1}$
$k_{translation}$ Rate of chromoprotein synthesis $0.0011 \textrm{min}^{-1}$
$\gamma_m$ mRNA degradation rate $0.19 \textrm{min}^{-1}$
$\gamma_p$ Chromoprotein degradation rate $0.18 \textrm{min}^{-1}$
$K_{\textrm{NO}}$ Dissociation constant of NO $0.12 \textrm{min}^{-1}$
$K_{\textrm{d(NsrR)}}$ Dissociation constant of NsrR $\textrm{35}$ nM
n Cooperativity of NO binding to NsrR 2

Simulation results



Figure 4. The graph shows the steady–state concentration/expression of mCP over time with respect to the mRNA degradation rates of chromoprotein. It is clear from the graph that presence of 0.87 nM of mCP is important in the cell in order to achieve the steady–state concentration of mCP.


Figure 5. The graph shows the relationship between the steady–state concentration over time of CP with respect to its degradation rate. It can be seen from the graph that 0.005 nM of CP is required in the cell in order to acquire the steady–state concentration of CP.

Matlab Code

    
%Defining the parameters
ktr     = 0.167;     %mRNA synthesis rate of 
                     %Chromoprotein per minute (rate of transcription)
ktrans  = 0.0011;    %Chromoprotein synthesis rate 
                     %per minute (rate of translation)
ycp     = 0.18;      %Chromoproteins degradation rate per minute
ym      = 0.19;      %mRNA degradation rate of chromoproteins per minute 
Kno     = 0.12;      %Dissociation constant for inducer (NO)
Kdn     = 35;        %Dissociation constant for repressor (NsrR)
pyear   = 1;         %Number of plasmids per cell
n       = 2;         %Cooperativity of inducer (NO)

%Assume that the repressor is removed 
NsrR    = 0;
NO      = 2;
t       = 0:0.1:240; 

%Calculate mRNA and Chromoprotein values
K       = ktr*pyear/(1+(NsrR/((1+(NO/Kno)^n)*Kdn)));
mCP     = K/ym-K/ym*exp(-ym*t);
CP      = K*ktrans/(ycp*ym)+K*ktrans/(ym^2-ycp*ym)*exp(-ym*t)...
    -(K*ktrans/(ycp*ym)+K*ktrans/(ym^2-ycp*ym))*exp(-ycp*t);

disp('Steady State mCP Concentration:')
disp(mCP(end)) %Steady state mRNA value
disp('Steady State CP Concentration:')
disp(CP(end)) %Steady state chromoprotein value
 
%Plotting mRNA concentration over time
figure(1)
plot(t, mCP,'r','LineWidth',2)
hold on
grid on
xlabel('Time (minutes)')
ylabel('mCP per cell (nM)')
title('Expression of mRNA of chromoprotein (mCP)')

%Plotting the chromoprotein concentrations over time
figure(2)
plot(t,CP,'r', 'LineWidth',2)
grid on
hold on
xlabel('Time (minutes)')
ylabel('CP per cell (nM)')
title('Expression of chromoprotein (CP)')

%Save this as nox_sensing_system_simulation.m file
    

Conclusion

The system is considered as stable, because when NsrR detaches itself from pYeaR, the promoter then initiate the production of mCP and CP, until a steady–state concentration is reached. This steady–state concentration is the maximum production rate of mCP and CP with respect to their degradation rates. However, it is clear from the graphs that the degradation rates are only related to the maximum production rate/ number of molecules in the cell; it is not responsible for altering the time span required to reach the steady–state concentration. Furthermore, this simplified solution for NOx–sensing system demonstrates the migration of system towards steady–state under assumed parameters and prescribed conditions. Therefore, this mathematical model can be considered as an appropriate model for determining the behavior of NOx–sensing system.

References

  1. Li, J., Zhao, Z., Kazakov, A., Chaos, M., Dryer, F.L. and Scire, J.J. (2007) 'A comprehensive kinetic mechanism for CO, CH2O, and CH3OH combustion', International Journal of Chemical Kinetics, 39(3), pp. 109–136. doi: 10.1002/kin.20218.

  2. Chen, H., Shiroguchi, K., Ge, H. and Xie, X.S. (2015) 'Genome-wide study of mRNA degradation and transcript elongation in Escherichia coli', Molecular Systems Biology, 11(1), pp. 781–781. doi: 10.15252/msb.20145794.

  3. Leduc, J., Thorsteinsson, M.V., Gaal, T. and Roberts, G.P. (2001) 'Mapping CooA.RNA polymerase interactions: Identification of activating regions 2 and 3 in CooA, the co–sensing transcriptional activator', Journal of Biological Chemistry, 276(43), pp. 39968–39973. doi: 10.1074/jbc.m105758200.

  4. Benabbas, A., Karunakaran, V., Youn, H., Poulos, T.L. and Champion, P.M. (2012) 'Effect of DNA binding on Geminate CO recombination Kinetics in CO-sensing transcription factor CooA', Journal of Biological Chemistry, 287(26), pp. 21729–21740. doi: 10.1074/jbc.m112.345090.

  5. Kuchinskas, M., Li, H., Conrad, M., Roberts, G. and Poulos, T.L. (2006) 'The role of the DNA-Binding domains in CooA activation', Biochemistry, 45(23), pp. 7148–27153. doi: 10.1021/bi052609o.

  6. Bergmann, John E, and Harvey F Lodish. 'A kinetic model of protein synthesis. Application to hemoglobin synthesis and translational control.' The Journal of biological chemistry 254.23 (1979) : 11927–37.

  7. Nick Csicsery, and Ricky O'Laughlin. 'A Mathematical Model of a Synthetically Constructed Genetic Toggle Switch' (2015)