Team:Imperial/CBD Kinetics Model
From 2014.igem.org
(11 intermediate revisions not shown) | |||
Line 11: | Line 11: | ||
<div id="subNav"> | <div id="subNav"> | ||
<ul> | <ul> | ||
- | <li><a data-scroll href="#introduction" >Introduction</a> | + | <li><a data-scroll href="#introduction" ">Introduction</a> |
</li> | </li> | ||
- | + | ||
<li><a data-scroll href="#maths" >Mathematical Background</a> | <li><a data-scroll href="#maths" >Mathematical Background</a> | ||
</li> | </li> | ||
Line 20: | Line 20: | ||
<li><a data-scroll href="#results" >Results & Conclusions</a> | <li><a data-scroll href="#results" >Results & Conclusions</a> | ||
</li> | </li> | ||
- | + | <li><a data-scroll href="#references" >References</a> | |
</li> | </li> | ||
- | + | <li><a data-scroll href="#code" >Appendix: Code</a> | |
</li> | </li> | ||
+ | |||
</ul> | </ul> | ||
<br class="clearLeft" /> | <br class="clearLeft" /> | ||
Line 31: | Line 32: | ||
<section id="introduction"> | <section id="introduction"> | ||
<h2>Introduction</h2> | <h2>Introduction</h2> | ||
- | <p>Based on a set of ordinary differential equations describing the kinetics of | + | <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> | ||
<section id="maths"> | <section id="maths"> | ||
<h2>Mathematical Background</h2> | <h2>Mathematical Background</h2> | ||
- | |||
<ul> | <ul> | ||
- | <li>Variables: | + | <li>Dependent Variables: |
<ul> | <ul> | ||
- | <li>CBS( | + | |
- | + | <li>CBS(t): concentration of unoccupied cellulose binding sites (CBSs) at time t</li> | |
- | + | ||
- | <li> | + | <li>CBD(t): concentration of unbound CBD at time t</li> |
- | <li> | + | <li>CBS_CBD(t): concentration of binding complex (CBS bounded with CBD) at time t</li> |
- | + | ||
</ul> | </ul> | ||
</li> | </li> | ||
<li>Parameters: | <li>Parameters: | ||
- | + | <ul> | |
- | + | <li>CBS(0): initial concentration of CBS</li> | |
- | <li>k<sub>off</sub>: dissociation rate constant</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> | </ul> | ||
</li> | </li> | ||
<li>Binding Kinetics: | <li>Binding Kinetics: | ||
<ul> | <ul> | ||
- | <li><img class="content-image " src="https://static.igem.org/mediawiki/2014/8/89/IC14-modelling-kinetics-binding1.png" height= | + | <li> |
- | <li>Assume | + | <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>where k<sub>on</sub> | + | <li>CBS_CBD(∞) = CBS(0), assume that BC binding sites can be fully saturated given enough time, as k<sub>on</sub>CBD(0) >> 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) >> k<sub>off</sub></li> | ||
</ul></li> | </ul></li> | ||
</ul> | </ul> | ||
- | + | ||
</ul> | </ul> | ||
Line 86: | Line 99: | ||
</section> | </section> | ||
<section id="lab"> | <section id="lab"> | ||
- | <h2> | + | <h2>Interaction with wetlab experiments</h2> |
<figure class="content-image image-right image-small"> | <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"> | <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> | <figcaption>How our model and the wetlab interact with each other</figcaption> | ||
</figure> | </figure> | ||
- | <p>Parameters | + | <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> | |
- | </p> | + | |
- | + | ||
- | + | ||
- | + | ||
</section> | </section> | ||
<section id="results"> | <section id="results"> | ||
<h2>Results and Conclusions</h2> | <h2>Results and Conclusions</h2> | ||
- | + | ||
<figure class="content-image image-half"> | <figure class="content-image image-half"> | ||
<img class="image-full" src="https://static.igem.org/mediawiki/2014/7/7c/IC14-modelling-kinetics-results1.png"> | <img class="image-full" src="https://static.igem.org/mediawiki/2014/7/7c/IC14-modelling-kinetics-results1.png"> | ||
- | <figcaption> | + | <figcaption>Figure 1: CBD-BC complex concentration against time |
- | </figcaption> | + | </figcaption> |
</figure> | </figure> | ||
<figure class="content-image image-half"> | <figure class="content-image image-half"> | ||
<img class="image-full" src="https://static.igem.org/mediawiki/2014/8/81/IC14-modelling-kinetics-results2.png"> | <img class="image-full" src="https://static.igem.org/mediawiki/2014/8/81/IC14-modelling-kinetics-results2.png"> | ||
- | <figcaption> | + | <figcaption>Figure 2: Time for 90% saturation of BC binding sites against initial concentration of CBD</figcaption> |
</figure> | </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> | ||
- | + | ||
<section id="references"> | <section id="references"> | ||
- | < | + | <h2>References</h2> |
<ol> | <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>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> | + | </li> |
- | + | ||
Line 131: | Line 145: | ||
</ol> | </ol> | ||
</section> | </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> | ||
Latest revision as of 00:13, 18 October 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
![](https://static.igem.org/mediawiki/2014/8/8c/IC14-modelling-kinetics-synergy.png)
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
![](https://static.igem.org/mediawiki/2014/7/7c/IC14-modelling-kinetics-results1.png)
![](https://static.igem.org/mediawiki/2014/8/81/IC14-modelling-kinetics-results2.png)
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
- 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