Team:Imperial/CBD Kinetics Model

From 2014.igem.org

(Difference between revisions)
 
(6 intermediate revisions not shown)
Line 32: 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 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>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, wetlab experiments can be informed about the functionalisation capability of BC and therefore would be able to approximate efficiency of the contaminants-capturing process.</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>
Line 42: Line 42:
                                 <ul>
                                 <ul>
-
                                     <li>CBS(t): concentration of unoccupied cellulose binding sites(CBS) at time t</li>
+
                                     <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>CBD(t): concentration of unbound CBD at time t</li>
Line 52: Line 52:
                             <li>Parameters:
                             <li>Parameters:
                                 <ul>
                                 <ul>
-
                                     CBS(0): initial concentration of CBS
+
                                     <li>CBS(0): initial concentration of CBS</li>
-
                                     <li>CBD(0): initial concentration of CBM injected into the system</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>on</sub>: association rate constant</li>
                                     <li>k<sub>off</sub>: dissociation rate constant</li>
                                     <li>k<sub>off</sub>: dissociation rate constant</li>
Line 67: Line 67:
                                         <img class="content-image " src="https://static.igem.org/mediawiki/2014/8/89/IC14-modelling-kinetics-binding1.png" height=150>
                                         <img class="content-image " src="https://static.igem.org/mediawiki/2014/8/89/IC14-modelling-kinetics-binding1.png" height=150>
                                     </li>
                                     </li>
-
                                     <li>Assume CBM(t) = CBM(0) due to large concentration of CBM injected into the system compared to initial concentration of cellulose binding sites.</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>
 +
                           
                             <li>Initial and Boundary Conditions:
                             <li>Initial and Boundary Conditions:
                                 <ul>
                                 <ul>
-
                                     <li>CBD(0)=M<sub>0</sub>, where M0 is a value determined as per experiment.</li>
+
                                     <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(0) = 0</li>
-
                                     <li>CBS_CBD(∞) = CBS(0), assume that BC binding sites can be fully saturated given enough time, as konCBD(0) &gt;&gt; koff.</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>
Line 83: Line 85:
                                         <img class="content-image " src="https://static.igem.org/mediawiki/2014/a/ac/IC14-modelling-kinetics-binding2.png" height=60>
                                         <img class="content-image " src="https://static.igem.org/mediawiki/2014/a/ac/IC14-modelling-kinetics-binding2.png" height=60>
                                     </li>
                                     </li>
-
                                    <li>where k<sub>on</sub> CBD(0) &gt;&gt; k<sub>>off</sub</li>
+
                                        <li>where k<sub>on</sub> CBD(0) &gt;&gt; k<sub>off</sub></li>
                                         </ul></li>
                                         </ul></li>
                                  
                                  
Line 102: Line 104:
                             <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 koff and kd 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>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 a posteriori checked through further wetlab experiments (please refer to CBD expression experiments for further information).</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>The chart below illustrates how our model and the wetlab experiments interact with each other.</p>
Line 123: Line 125:
                         <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>Figure 2: Time for 90% saturation of BC binding sites against initial solution concentration CBD</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>
                         <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>
Line 146: Line 148:
<h2>Appendix: Code</h2>
<h2>Appendix: Code</h2>
-
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>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>
</section>

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