Team:Imperial/CBD Kinetics Model

From 2014.igem.org

(Difference between revisions)
(Created page with "{{:Team:Imperial/Templates:header}} <html> <body> </body> </html> {{:Team:Imperial/Templates:footer}}")
 
(13 intermediate revisions not shown)
Line 2: Line 2:
<html>
<html>
<body>
<body>
 +
<div class="content-wrapper">
 +
 +
            <div class="pure-g">
 +
 +
                <div class="pure-u-1-1 main">
 +
                    <h1>CBD Kinetics Model</h1>
 +
                    <div id="navWrap">
 +
                        <div id="subNav">
 +
                            <ul>
 +
                                <li><a data-scroll href="#introduction" ">Introduction</a>
 +
                                </li>
 +
 +
                                <li><a data-scroll href="#maths" >Mathematical Background</a>
 +
                                </li>
 +
                                <li><a data-scroll href="#lab" >Wet-Lab Interaction</a>
 +
                                </li>
 +
                                <li><a data-scroll href="#results" >Results & Conclusions</a>
 +
                                </li>
 +
                                <li><a data-scroll href="#references" >References</a>
 +
                                </li>
 +
<li><a data-scroll href="#code" >Appendix: Code</a>
 +
                                </li>
 +
                               
 +
                            </ul>
 +
                            <br class="clearLeft" />
 +
                        </div>
 +
                    </div>
 +
 +
                    <section id="introduction">
 +
                        <h2>Introduction</h2>
 +
                        <p>Based on a set of ordinary differential equations (ODEs) describing the kinetics of carbohydrate-binding domains (CBDs) attaching to bacterial cellulose (BC), the model is designed to predict the efficiency and affinity of the binding process. Given an initial concentration of CBD injected into the system, the model calculates the time necessary for a certain percentage of protein binding sites on BC to be saturated.</p>
 +
                        <p>With a better understanding of the interaction between CBD and BC, <a href="https://2014.igem.org/Team:Imperial/Functionalisation">wetlab experiments</a> can be informed about the functionalisation capability of BC and therefore would be able to approximate efficiency of the contaminants-capturing process.</p>
 +
 +
                    </section>
 +
                    <section id="maths">
 +
                        <h2>Mathematical Background</h2>
 +
                        <ul>
 +
                            <li>Dependent Variables:
 +
                                <ul>
 +
 +
                                    <li>CBS(t): concentration of unoccupied cellulose binding sites (CBSs) at time t</li>
 +
 +
                                    <li>CBD(t): concentration of unbound CBD at time t</li>
 +
                                    <li>CBS_CBD(t): concentration of binding complex (CBS bounded with CBD) at time t</li>
 +
 +
 +
                                </ul>
 +
                            </li>
 +
                            <li>Parameters:
 +
                                <ul>
 +
                                    <li>CBS(0): initial concentration of CBS</li>
 +
                                    <li>CBD(0): initial concentration of CBD injected into the system</li>
 +
                                    <li>k<sub>on</sub>: association rate constant</li>
 +
                                    <li>k<sub>off</sub>: dissociation rate constant</li>
 +
                                    <li>K<sub>D</sub> = k<sub>off</sub>/k<sub>on</sub>: equilibrium dissociation constant</li>
 +
 +
 +
 +
                                </ul>
 +
                            </li>
 +
                            <li>Binding Kinetics:
 +
                                <ul>
 +
                                    <li>
 +
                                        <img class="content-image " src="https://static.igem.org/mediawiki/2014/8/89/IC14-modelling-kinetics-binding1.png" height=150>
 +
                                    </li>
 +
                                    <li>Assume that CBD(t) = CBD(0) since the initial concentration of CBD injected into the system is typically very large compared to the initial concentration of cellulose binding sites, i.e. CBS(0).</li>
 +
                            </ul>
 +
                            </li>
 +
                           
 +
                            <li>Initial and Boundary Conditions:
 +
                                <ul>
 +
                                    <li>CBD(0) = M<sub>0</sub>, where M<sub>0</sub> is a value determined as per experiment.</li>
 +
                                    <li>CBS_CBD(0) = 0</li>
 +
                                    <li>CBS_CBD(∞) = CBS(0), assume that BC binding sites can be fully saturated given enough time, as k<sub>on</sub>CBD(0) &gt;&gt; k<sub>off</sub>.</li>
 +
 +
 +
                                </ul>
 +
                            </li>
 +
                            <li>Solution:
 +
                                <ul>
 +
                                    <li>
 +
                                        <img class="content-image " src="https://static.igem.org/mediawiki/2014/a/ac/IC14-modelling-kinetics-binding2.png" height=60>
 +
                                    </li>
 +
                                        <li>where k<sub>on</sub> CBD(0) &gt;&gt; k<sub>off</sub></li>
 +
                                        </ul></li>
 +
                               
 +
                               
 +
                                </ul>
 +
                           
 +
 +
                        </ul>
 +
                     
 +
 +
 +
 +
                    </section>
 +
                    <section id="lab">
 +
                        <h2>Interaction with wetlab experiments</h2>
 +
                        <figure class="content-image image-right image-small">
 +
                            <img class="image-full" src="https://static.igem.org/mediawiki/2014/8/8c/IC14-modelling-kinetics-synergy.png">
 +
                            <figcaption>How our model and the wetlab interact with each other</figcaption>
 +
                        </figure>
 +
                        <p>Parameters k<sub>off</sub> and k<sub>d</sub> are determined by wetlab assays (wash experiments), from which k<sub>on</sub> is calculated: k<sub>on</sub>=k<sub>off</sub>/k<sub>d</sub>.</p>
 +
 +
<p>With the parameters determined from the experiments, the model is then instantiated and used to predict the time period necessary for a certain percentage of the binding sites to be saturated. The validity of the model is <em>a posteriori</em> checked through further wetlab experiments (please refer to <a href="https://2014.igem.org/Team:Imperial/Functionalisation">CBD expression experiments</a> for further information).</p>
 +
 +
<p>The chart below illustrates how our model and the wetlab experiments interact with each other.</p>
 +
 +
 +
 +
 +
                    </section>
 +
 +
                    <section id="results">
 +
                        <h2>Results and Conclusions</h2>
 +
 +
                        <figure class="content-image image-half">
 +
                            <img class="image-full" src="https://static.igem.org/mediawiki/2014/7/7c/IC14-modelling-kinetics-results1.png">
 +
                            <figcaption>Figure 1: CBD-BC complex concentration against time
 +
                            </figcaption>
 +
                        </figure>
 +
                        <figure class="content-image image-half">
 +
                            <img class="image-full" src="https://static.igem.org/mediawiki/2014/8/81/IC14-modelling-kinetics-results2.png">
 +
                            <figcaption>Figure 2: Time for 90% saturation of BC binding sites against initial concentration of CBD</figcaption>
 +
                        </figure>
 +
                        <p>It can be observed from the second graph that the time necessary for 90% CBS saturation decreases significantly when initial CBD concentration increases beyond 0.5 µM; however, a further increase in initial CBD concentration does not contribute much to effectively reduce the time necessary for 90% CBS saturation. Therefore, an initial CBD concentration between 0.5 µM and 1 µM will be considered as optimal to efficiently functionalise BC.</p>
 +
 +
                    </section>
 +
 +
 +
                    <section id="references">
 +
                        <h2>References</h2>
 +
                        <ol>
 +
                            <li>Zhang, Mengmeng; Wang, Bin; Xu, Bingqian Mapping Single Molecular Binding Kinetics of Carbohydrate-Binding Module with Crystalline Cellulose by Atomic Force Microscopy Recognition Imaging J. Phys. Chem. B 2014, 118, 6714−6720s
 +
                            </li>
 +
 +
 +
 +
 +
 +
 +
                        </ol>
 +
                    </section>
 +
<section id="code">
 +
<h2>Appendix: Code</h2>
 +
 +
<p>The code for this model can be found here: <a href="https://static.igem.org/mediawiki/2014/1/16/CBD_kinetics_model_imperial_iGEM_2014.m">CBD_kinetics_model_imperial_iGEM_2014.m</a></p>
 +
</section>
 +
 +
                </div>
 +
 +
 +
 +
            </div>
 +
 +
 +
 +
 +
 +
 +
 +
 +
 +
        </div>
</body>
</body>
</html>
</html>
{{:Team:Imperial/Templates:footer}}
{{:Team:Imperial/Templates:footer}}

Latest revision as of 00:13, 18 October 2014

Imperial iGEM 2014

CBD Kinetics Model

Introduction

Based on a set of ordinary differential equations (ODEs) describing the kinetics of carbohydrate-binding domains (CBDs) attaching to bacterial cellulose (BC), the model is designed to predict the efficiency and affinity of the binding process. Given an initial concentration of CBD injected into the system, the model calculates the time necessary for a certain percentage of protein binding sites on BC to be saturated.

With a better understanding of the interaction between CBD and BC, wetlab experiments can be informed about the functionalisation capability of BC and therefore would be able to approximate efficiency of the contaminants-capturing process.

Mathematical Background

  • Dependent Variables:
    • CBS(t): concentration of unoccupied cellulose binding sites (CBSs) at time t
    • CBD(t): concentration of unbound CBD at time t
    • CBS_CBD(t): concentration of binding complex (CBS bounded with CBD) at time t
  • Parameters:
    • CBS(0): initial concentration of CBS
    • CBD(0): initial concentration of CBD injected into the system
    • kon: association rate constant
    • koff: dissociation rate constant
    • KD = koff/kon: equilibrium dissociation constant
  • Binding Kinetics:
    • Assume that CBD(t) = CBD(0) since the initial concentration of CBD injected into the system is typically very large compared to the initial concentration of cellulose binding sites, i.e. CBS(0).
  • Initial and Boundary Conditions:
    • CBD(0) = M0, where M0 is a value determined as per experiment.
    • CBS_CBD(0) = 0
    • CBS_CBD(∞) = CBS(0), assume that BC binding sites can be fully saturated given enough time, as konCBD(0) >> koff.
  • Solution:
    • where kon CBD(0) >> koff

Interaction with wetlab experiments

How our model and the wetlab interact with each other

Parameters koff and kd are determined by wetlab assays (wash experiments), from which kon is calculated: kon=koff/kd.

With the parameters determined from the experiments, the model is then instantiated and used to predict the time period necessary for a certain percentage of the binding sites to be saturated. The validity of the model is a posteriori checked through further wetlab experiments (please refer to CBD expression experiments for further information).

The chart below illustrates how our model and the wetlab experiments interact with each other.

Results and Conclusions

Figure 1: CBD-BC complex concentration against time
Figure 2: Time for 90% saturation of BC binding sites against initial concentration of CBD

It can be observed from the second graph that the time necessary for 90% CBS saturation decreases significantly when initial CBD concentration increases beyond 0.5 µM; however, a further increase in initial CBD concentration does not contribute much to effectively reduce the time necessary for 90% CBS saturation. Therefore, an initial CBD concentration between 0.5 µM and 1 µM will be considered as optimal to efficiently functionalise BC.

References

  1. Zhang, Mengmeng; Wang, Bin; Xu, Bingqian Mapping Single Molecular Binding Kinetics of Carbohydrate-Binding Module with Crystalline Cellulose by Atomic Force Microscopy Recognition Imaging J. Phys. Chem. B 2014, 118, 6714−6720s

Appendix: Code

The code for this model can be found here: CBD_kinetics_model_imperial_iGEM_2014.m