Difference between revisions of "Team:UNITN-Trento/Modeling"

m
Line 83: Line 83:
 
<div class="content">
 
<div class="content">
 
<header>
 
<header>
<h3 class="wow fadeInDown">Functions used</h3>
+
<h3 class="wow fadeInDown">MATLAB code used</h3>
 
</header>
 
</header>
<p>ciao</p>
+
<div class="captionbox" style="width=80%;margin:auto;"><img src="https://static.igem.org/mediawiki/2015/1/10/UniTN_pics_2015_code2.png"/></div><br/>
 +
To produce the graph
 +
<div class="captionbox" style="width=80%;margin:auto;"><img src="https://static.igem.org/mediawiki/2015/0/0e/UniTN_pics_2015_code1.png"/></div>
 
</div>
 
</div>
 
</section>
 
</section>
 
<section class="wrapper style4 container" style="margin-top:1em;">
 
<div class="content">
 
<header>
 
<h3 class="wow fadeInDown">What we used to produce the graphs</h3>
 
</header>
 
<p>ciao</p>
 
</div>
 
</section>
 
 
 
<section class="wrapper style4 container">
 
<section class="wrapper style4 container">
 
<!-- Content -->  
 
<!-- Content -->  
Line 103: Line 95:
 
<h3 class="wow fadeInDown">References</h3>
 
<h3 class="wow fadeInDown">References</h3>
 
</header>
 
</header>
 
 
<ol type="1" class="sourcebox">  
 
<ol type="1" class="sourcebox">  
 
<a class="anchor-off" name="refs_1" id="refs_1"></a>  
 
<a class="anchor-off" name="refs_1" id="refs_1"></a>  

Revision as of 18:55, 17 September 2015

Modeling

Can we relate light incidence to ATP production?

Introduction

We really want to thank the team iGEM Kent 2015 for this modeling born from the collaboration between the two groups.

Proteorhodopsin is a light-powered proton pump. This protein has the property to use light energy to generate an outward proton flux that can subsequently power cellular processes, such as ATP synthesis.After having experimentally assessed that our engineered E. coli was able to produce more ATP thanks to the proteorhodopsin, we decided to quantify this effect by modeling the amount of ATP that can be produced in function of the proteorhodopsin’s activation.

How many ATP molecules?

Membrane as an electric circuitElectric circuit analogy for the membrane [2]

What is the number of ATP molecules that can be produced per second as a function of light irradiance that hits the bacterial membrane?
Once a photon is absorbed by proteorhodopsin (PR), PR must complete its photocycle before it can absorb another photon [1]. At high light irradiance, this leads to saturation. For this we choose to exploit the Michaelis-Menten kinetics, where V_max is the maximum rate of the system and the Michaelis-Menten constant, Km, is the substrate concentration at which the reaction rate is $\frac{1}{2}V\max$.


Walter et al. demonstrated that the system is analogous to a circuit (figure 1), in this circuit representation; the proteorhodopsin (PR) acts like a battery with internal resistance. [2][3]

The current through the system is inversely related to the PR resistor and is dependent on light irradiance.

$R_{PR=}\left( \frac{V_{\max }*I}{K_{m}+I} \right)^{-1}$

Walter et al. determined that $V_{\max}$ is fixed by the boundary condition that $R_{PR≈}\frac{R_{\sin k}}{10}$ at the highest light irradiance $I=\frac{160mW}{cm^{2}}$. $\; R_{\sin k}≈R_{\mbox{re}s}≈10^{15}\; \Omega$ and $K_{m=}\frac{60mW}{cm^{2}}$. Where light irradiance of $\frac{20mW}{cm^{2}}\;$ is roughly equivalent to PR absorption from solar illumination at sea level. [2]

At the boundary condition:

$Rpr=\frac{R_{\sin k}}{10}=10^{14}\Omega =\left( \frac{V_{\max }*I}{K_{m}+I} \right)^{-1}$

Hence:

$V\max \; =\; \frac{K_{m}+I}{R_{PR}*I}\; =\; 1.375*10^{-14}\; \Omega ^{-1}$

The rate of reaction, $v$, has units of $\Omega ^{\left( -1 \right)}$; through dimensional analysis we can see that $\Omega ^{-1}\; =\; \frac{Amps}{Volts}\; =\frac{coulombs}{\left( \sec ond*voltage \right)}$. The voltage across the PR, $V_{PR=}0.2\; Volts\;$ [2] and the charge of a proton is $q=1.6*10^{\left( -19 \right)}\; \mbox{C}$.

Therefore we can work out the number of protons pumped by the PR per second as

$N_{Proton}\; =\; \frac{V_{\max }\; I}{K_{m}+I}*\frac{V_{PR}}{q}$

If an electron pair is composed of 10 protons and there is a net gain of 2.5 ATP molecules per electron pair then the number of ATP molecules produced per second is simply:

$N_{ATP}\; =\; \frac{1}{4}N_{Proton}\; =\; \frac{1}{4}*\frac{V_{\max }\; I}{K_{m}+I}*\frac{V_{PR}}{q}$

The rate of ATP production per second per bacterium as a function of light irradiance has been plotted in figure 2. From the graph, most ATP production rates per second per bacterium are in the range 102-103, after 5 minutes of illumination each cell would have produced a net gain of about 105 ATP molecules, which agrees with experiment [4].

A Michaelis-Menten curveATP production per second per bacterium as a function of a irradiance of light.

MATLAB code used


To produce the graph

References

  1. Béja, O., Aravind, L., Koonin, E. V., Suzuki, M. T., Hadd, A., Nguyen, L. P., DeLong, E. F. (2000).
    Bacterial rhodopsin: evidence for a new type of phototrophy in the sea.
    Science, 289(5486), 1902-1906.
  2. Walter, J. M., Greenfield, D., Bustamante, C., & Liphardt, J. (2007).Light-powering Escherichia coli with proteorhodopsin.
    Proceedings of the National Academy of Sciences, 104(7), 2408-2412.
  3. Friedrich, T., Geibel, S., Kalmbach, R., Chizhov, I., Ataka, K., Heberle, J., & Bamberg, E. (2002).Proteorhodopsin is a light-driven proton pump with variable vectoriality.
    Journal of molecular biology, 321(5), 821-838.
  4. Martinez, A., Bradley, A. S., Waldbauer, J. R., Summons, R. E., & DeLong, E. F. (2007). Proteorhodopsin photosystem gene expression enables photophosphorylation in a heterologous host.
    Proceedings of the National Academy of Sciences, 104(13), 5590-5595.