Team:Groningen/Modelling/Arsenic

From 2009.igem.org

(Difference between revisions)
m (Lijstje in Tabel ergens onderaan)
m
Line 279: Line 279:
The following variables play an important role in our system (these can be concentrations of substances, the density of the cell, etc.):
The following variables play an important role in our system (these can be concentrations of substances, the density of the cell, etc.):
-
 
+
{|class="ourtable" style="clear:right;"
 +
|+ Core Substances
 +
!colspan="2"|Name
 +
!Description
 +
!Derivative to time
 +
|-
 +
|colspan="4"|''Extracellular''
 +
|-
 +
| ||As(III)<sub>ex</sub>T || As(III) in the solution. || (Vc/Vs) k8 ArsB<sub>As</sub> - (Vc/Vs) v5<sup>&dagger;</sup> As(III)<sub>ex</sub>T / (K5+As(III)<sub>ex</sub>T)
 +
|-
 +
|colspan="4"|''Membrane (all naturally occurring, but we plan to bring GlpF to overexpression)''
 +
|-
 +
| ||GlpFT || Importer of As(III) (concentration w.r.t. the exterior of the cell). || (not used directly in model, assumed to be constant)
 +
|-
 +
| ||ArsBT || Exporter of As(III) (concentration w.r.t. the interior of the cell). || βB ars1 - (ln(2)/τB) ArsB
 +
|-
 +
|colspan="4"|''Intracellular (ars2, pro and GV are introduced)''
 +
|-
 +
| ||As(III)<sub>in</sub>T || As(III) (bound and unbound) in the cell. || v5 As(III)<sub>ex</sub>T / (K5+As(III)<sub>ex</sub>T) - k8 ArsB<sub>As</sub>
 +
|-class="estimate"
 +
| ||ars1T || ArsR repressed promoters (bound and unbound) naturally occurring in E. coli. || (concentration is constant = 1.6605nM, one per cell)
 +
|-class="estimate"
 +
| ||ars2T || ArsR repressed promoters in front of gas vesicle genes. || (concentration is constant = 0-166.05nM)
 +
|-class="estimate"
 +
| ||proR || Constitutive promoters in front of arsR. || (concentration is constant = 0-166.05nM)
 +
|-class="estimate"
 +
| ||proM || Constitutive promoters in front of mbp-arsR. || (concentration is constant = 0-166.05nM)
 +
|-class="estimate"
 +
| ||proF || Constitutive promoters in front of fMT. || (concentration is constant = 0-166.05nM)
 +
|-
 +
| ||ArsRT || ArsR in the cell. || βRN ars1 + βR proR - (ln(2)/τR) ArsR
 +
|-
 +
| ||MBPArsRT || MBPArsR in the cell. || βM proM - (ln(2)/τM) MBPArsR
 +
|-
 +
| ||fMTT || fMT in the cell. || βF proF - (ln(2)/τF) fMT
 +
|-
 +
| ||GV || Concentration of gas vesicles. || βG ars2 - (ln(2)/τG) GV
 +
|-style="border:none;"
 +
|colspan="4"|
 +
{|class="ourtable" style="width:100%"
 +
!colspan="5"|
 +
|- style="text-align:center;"
 +
|class="fromPaper" style="padding:0;"|Directly from paper.
 +
|class="selfDerived" style="padding:0;"|Based on data from paper.
 +
|class="experimental" style="padding:0;"|Based on experiment.
 +
|class="estimate" style="padding:0;"|Rough estimate.
 +
|class="unknown" style="padding:0;"|Totally unknown.
 +
|}
 +
|}
* Extracellular:
* Extracellular:
** As(III)<sub>ex</sub>
** As(III)<sub>ex</sub>

Revision as of 16:59, 14 October 2009

[http://2009.igem.org/Team:Groningen http://2009.igem.org/wiki/images/f/f1/Igemhomelogo.png]

Based on the quasi-steady-state derivation below we have made the simplified version of our model shown below. The simplification is based on two key assumptions (which are also illustrated below, next to the table "Breakdown of core substances"):

  • Binding and unbinding of arsenic to/from the transporters occurs on a much smaller time scale than changes in the concentration of arsenic inside and outside the cell. And similarly, we assume that (un)binding of ArsR to/from the ars promoter is much faster than the production of ArsR (for example).
  • The concentration of transporters is insignificant compared to the concentration of arsenic inside and outside the cell, and similarly for the concentration of ars promoters compared to the concentration of ArsR. The concentration of promoters might not be insignificant.

This leads to the Michaelis-Menten equation for import, but also some more general equations for export using ArsB and accumulation with ArsR (for example, the Hill equation can be recognized in the activity of the ars promoter).

A schematic representation of the processes involved in arsenic filtering (keep in mind that ArsR represses the expression of the genes behind ars). Note that MBPArsR and fMT are not shown for clarity.
Reactions
Reaction Description Rate
Transport
As(III)exT → As(III)inTImport of arsenic.(Vc/Vs) v5 As(III)exT / (K5+As(III)exT)
As(III)inT → As(III)exTExport of arsenic. k8 ArsBAs
ars1T → ars1T + ArsBTProduction of ArsB. βB ars1
ArsBT → nullDegradation of ArsB (ln(2)/τB) ArsB
Accumulation
ars1T → ars1T + ArsRTFrom chromosomal operon. βRN ars1
proR → proR + ArsRTProduction of ArsR. βR pro
proM → proM + MBPArsRTProduction of MBPArsR. βM pro
proF → proF + fMTTProduction of fMT. βF pro
ArsRT → nullDegradation of ArsR. (ln(2)/τR) ArsR
MBPArsRT → nullDegradation of MBPArsR. (ln(2)/τM) MBPArsR
fMTT → nullDegradation of fMT. (ln(2)/τF) fMT
Gas vesicles
ars2T → ars2T + GVTranscription + translation. βG ars2
GV → nullDegradation of gas vesicles. (ln(2)/τG) GV
Core Substances
Name Description Derivative to time
Extracellular
As(III)exT As(III) in the solution. (Vc/Vs) k8 ArsBAs - (Vc/Vs) v5 As(III)exT / (K5+As(III)exT)
Membrane (all naturally occurring, but we plan to bring GlpF to overexpression)
GlpFT Importer of As(III) (concentration w.r.t. the exterior of the cell). (not used directly in model, assumed to be constant)
ArsBT Exporter of As(III) (concentration w.r.t. the interior of the cell). βB ars1 - (ln(2)/τB) ArsB
Intracellular (ars2, pro and GV are introduced)
As(III)inT As(III) (bound and unbound) in the cell. v5 As(III)exT / (K5+As(III)exT) - k8 ArsBAs
ars1T ArsR repressed promoters (bound and unbound) naturally occurring in E. coli. (concentration is constant = 1.6605nM, one per cell)
ars2T ArsR repressed promoters in front of gas vesicle genes. (concentration is constant = 0-166.05nM)
proR Constitutive promoters in front of arsR. (concentration is constant = 0-166.05nM)
proM Constitutive promoters in front of mbp-arsR. (concentration is constant = 0-166.05nM)
proF Constitutive promoters in front of fMT. (concentration is constant = 0-166.05nM)
ArsRT ArsR in the cell. βRN ars1 + βR proR - (ln(2)/τR) ArsR
MBPArsRT MBPArsR in the cell. βM proM - (ln(2)/τM) MBPArsR
fMTT fMT in the cell. βF proF - (ln(2)/τF) fMT
GV Concentration of gas vesicles. βG ars2 - (ln(2)/τG) GV
Directly from paper. Based on data from paper. Based on experiment. Rough estimate. Totally unknown.
Note that the "constant" v5 depends on the concentration of GlpF transporters in the cell, and this can depend on whether we bring GlpF to overexpression or not. For simplicity the production/degradation of GlpF is not included explicitly in the model, instead we can vary the constant v5 relative to the value found for wild-type E. coli.
Breakdown of core substances
Core substance Component Relative abundance
ArsBT ArsB K7
ArsBAs As(III)in
As(III)inT As(III)in 1
ArsRAs ArsRT / (KRd + As(III)in)
MBPArsRAs MBPArsRT / (KMd + As(III)in)
fMTAs fMTT As(III)innf-1 / (KFdnf + As(III)innf)
arsT ars K3d²
ArsRars ArsR²
ars ars1 ars1T
ars2 ars2T
ArsRT ArsR KRd
ArsRAs As(III)in
MBPArsRT MBPArsR KMd
MBPArsRAs As(III)in
fMTT fMT KFdnf
fMTAs As(III)innf
Circles correspond to core substances. We consider the reactions between the overlapping substances so fast that we model them by determining the ratios between the substances when the reactions between them are in equilibrium. Also, the complexes formed with ars, GlpF and ArsB (the small circles) are considered to have such a low concentration that they are of no importance to the concentrations of As(III)in/-ex and ArsR (the large circles).
Constants
Name Units Value Description
k8 1/s Reaction rate constant representing how fast ArsB can export arsenic.
KRd M 6µM Dissociation constant for ArsR and As(III). Assumed to be about an order of magnitude smaller than KDd = 60µM, the corresponding constant for the similar protein ArsD from Chen1997.
KMd M 6µM Dissociation constant for MBPArsR and As(III). We assume this to be roughly equal to KRd.
KFd M Dissociation constant for fMT and As(III).
nf Hill coefficient for the formation of the complex fMTAs. This is related to the number of arsenic ions that bind to fMT.
K3d M 0.33µM Dissociation constants for ArsR and ars.
  • K3d² = k3off/k3on = (0.33µM)²? (Chen1997, suspect as the relevant reference doesn't actually seem to give any value for this)
v5 mol/(s·L) 3.1863µmol/(s·L) Maximum import rate per liter of cells (see Michaelis-Menten equation). Note that we have purposefully chosen to write the units as mol/(s·L) instead of M/s, to emphasize the fact that the rate is per liter of cells.
  • v5 = k6 GlpFT (Vs/Vc)
K5 M 27.718µM Concentration at which import reaches half its maximum import rate (see Michaelis-Menten equation).
  • K5 = (k5off+k6) / k5on
K7 M Concentration at which export reaches half its maximum export rate (see Michaelis-Menten equation).
  • K7 = (k7off+k8) / k7on
τB, τR, τG, etc. s Half-lifes (of ArsB, ArsR and GV, respectively). Degradation rate = ln(2)/τ
If you take just the degradation into account you will have the equation dC/dt = -k*C, which leads to C(t) = C(0) e-k t. So if k = ln(2)/τ we get C(t) = C(0) e-ln(2)/τ t = C(0) 2-t/τ. In other words τ is the time it takes for the concentration to half.
 i 
βB, βR, etc. 1/s Production rates.
  • βRN = the production rate for ArsR behind the ars1 promoter
  • βB = the production rate for ArsB behind the ars1 promoter
  • βG = the production rate for GV behind the ars2 promoter
  • βR = the production rate for ArsR behind a constitutive promoter
  • βM = the production rate for MBPArsR behind a constitutive promoter
  • βF = the production rate for fMT behind a constitutive promoter
Vs L Volume of solution (excluding cells).
Vc L Total volume of cells (in solution) (so Vs+Vc is the total volume).
Directly from paper. Based on data from paper. Based on experiment. Rough estimate. Totally unknown.

The raw model

The following variables play an important role in our system (these can be concentrations of substances, the density of the cell, etc.):

Core Substances
Name Description Derivative to time
Extracellular
As(III)exT As(III) in the solution. (Vc/Vs) k8 ArsBAs - (Vc/Vs) v5 As(III)exT / (K5+As(III)exT)
Membrane (all naturally occurring, but we plan to bring GlpF to overexpression)
GlpFT Importer of As(III) (concentration w.r.t. the exterior of the cell). (not used directly in model, assumed to be constant)
ArsBT Exporter of As(III) (concentration w.r.t. the interior of the cell). βB ars1 - (ln(2)/τB) ArsB
Intracellular (ars2, pro and GV are introduced)
As(III)inT As(III) (bound and unbound) in the cell. v5 As(III)exT / (K5+As(III)exT) - k8 ArsBAs
ars1T ArsR repressed promoters (bound and unbound) naturally occurring in E. coli. (concentration is constant = 1.6605nM, one per cell)
ars2T ArsR repressed promoters in front of gas vesicle genes. (concentration is constant = 0-166.05nM)
proR Constitutive promoters in front of arsR. (concentration is constant = 0-166.05nM)
proM Constitutive promoters in front of mbp-arsR. (concentration is constant = 0-166.05nM)
proF Constitutive promoters in front of fMT. (concentration is constant = 0-166.05nM)
ArsRT ArsR in the cell. βRN ars1 + βR proR - (ln(2)/τR) ArsR
MBPArsRT MBPArsR in the cell. βM proM - (ln(2)/τM) MBPArsR
fMTT fMT in the cell. βF proF - (ln(2)/τF) fMT
GV Concentration of gas vesicles. βG ars2 - (ln(2)/τG) GV
Directly from paper. Based on data from paper. Based on experiment. Rough estimate. Totally unknown.
  • Extracellular:
    • As(III)ex
  • Membrane (all naturally occurring, but we plan to bring GlpF to overexpression):
    • GlpF (concentration w.r.t. the exterior of the cell)
    • GlpFAs (concentration w.r.t. the exterior of the cell)
    • ArsB (concentration w.r.t. the interior of the cell)
    • ArsBAs (concentration w.r.t. the interior of the cell)
  • Intracellular (ars2, pro and GV are introduced):
    • As(III)ex
    • ars1 (concentration of unbound promoters naturally occurring in E. coli)
    • ars2 (concentration of unbound promoters in front of gas vesicle genes)
    • proR (concentration of constitutive promoters in front of arsR)
    • proM (concentration of constitutive promoters in front of mbp-arsR)
    • proF (concentration of constitutive promoters in front of fMT)
    • ArsR
      ArsR binds to ars to repress production of the genes they regulate, and binds to As(III) to make it less of a problem for the cell.
       i 
    • ArsRAs (bound to As(III))
    • ArsRars1 (bound to ars1)
    • ArsRars2 (bound to ars2)
    • MBPArsR
      A fusion of maltose binding protein and ArsR.
       i 
    • MBPArsRAs (bound to As(III))
    • fMT
    • fMTAs (bound to multiple As(III))
    • ArsRAs (bound to As(III))
    • GV (concentration of gas vesicles)

The variables above can be related to each other through the following "reactions" (color coding is continued below to show which parts of the differential equations refer to which groups of reactions):

  • Transport (based on Rosen1996, Meng2004 and Rosen2009)
    • As(III)ex + GlpF ↔ GlpFAs
    • GlpFAs → GlpF + As(III)
    • As(III)in + ArsB ↔ ArsBAs
    • ArsBAs → ArsB + As(III)ex
    • ArsB → null (degradation)
  • Accumulation (mostly based on Chen1997)
    • As(III)in + ArsR ↔ ArsRAs
    • As(III)in + MBPArsR ↔ MBPArsRAs
    • nf As(III)in + fMT ↔ fMTAs
    • ars1 + 2 ArsR ↔ ArsRars1
    • ars2 + 2 ArsR ↔ ArsRars2
    • ars1 → ars1 + ArsR + ArsB (transcription + translation)
    • ars2 → ars2 + GV (transcription + translation)
    • proR → proR + ArsR (transcription + translation)
    • proM → proM + MBPArsR (transcription + translation)
    • proF → proF + fMT (transcription + translation)
    • ArsR → null (degradation)
    • GV → null (degradation)

Resulting in the following differential equations (please note that some can be formed by linear combinations of the others), using color coding to show the correspondence to the reactions above:

  • (d/dt) As(III)ex = - (d/dt) GlpFAs - k6 GlpFAs + (Vc/Vs) k8 ArsBAs
  • (d/dt) GlpF = - (d/dt) GlpFAs
  • (d/dt) GlpFAs = k5on As(III)ex GlpF - (k5off+k6) GlpFAs
  • (d/dt) ArsB = - (d/dt) ArsBAs + β4 ars1 - ln(2)/τB ArsB
  • (d/dt) ArsBAs = k7on As(III)in ArsB - (k7off+k8) ArsBAs
  • (d/dt) As(III)in = - (d/dt) ArsRAs - (d/dt) MBPArsRAs - nf (d/dt) fMTAs - (d/dt) ArsBAs - k8 ArsBAs + (Vs/Vc) k6 GlpFAs
  • (d/dt) ars1 = - (d/dt) ArsRars1
  • (d/dt) ars2 = - (d/dt) ArsRars2
  • (d/dt) ArsR = βRN ars1 + βR proR - (ln(2)/τR) ArsR - (d/dt) ArsRAs - 2 (d/dt) ArsRars1 - 2 (d/dt) ArsRars2
  • (d/dt) ArsRAs = kRon ArsR As(III)in - kRoff ArsRAs
  • (d/dt) ArsRars1 = k3on ArsR² ars1 - k3off ArsRars1
  • (d/dt) ArsRars2 = k3on ArsR² ars2 - k3off ArsRars2
  • (d/dt) MBPArsR = βM proM - (ln(2)/τM) MBPArsR - (d/dt) MBPArsRAs
  • (d/dt) MBPArsRAs = kMon MBPArsR As(III)in - kMoff MBPArsRAs
  • (d/dt) fMT = βF proF - (ln(2)/τF) fMT - (d/dt) fMTAs
  • (d/dt) fMTAs = kFon fMT As(III)innf - kFoff fMTAs
  • (d/dt) GV = βG ars2 - ln(2)/τG GV

Using the following constants/definitions:

Name Units Description
kRon, kMon, k5on, etc. 1/(M·s) Reaction rate constants for reactions to a complex.
k3on 1/(M²·s) Reaction rate constants for reactions to a complex.
kFon 1/(Mnf·s) Reaction rate constants for reactions to a complex.
kRoff, kMoff, kFoff, k3off, k5off, etc. 1/s Reaction rate constants for reactions from a complex.
k6, k8 1/s Reaction rate constants representing how fast transporters transport their cargo to "the other side".
τB, τR, τM, τF, τG s Half-lifes (of ArsB, ArsR, MBPArsR, fMT and GV, respectively). Degradation rate = ln(2)/τ
If you take just the degradation into account you will have the equation dC/dt = -k*C, which leads to C(t) = C(0) e-k t. So if k = ln(2)/τ we get C(t) = C(0) e-ln(2)/τ t = C(0) 2-t/τ. In other words τ is the time it takes for the concentration to half.
 i 
βRN, βR, etc. 1/s Production rates.
  • βRN = the production rate for ArsR behind the ars1 promoter
  • βB = the production rate for ArsB behind the ars1 promoter
  • βG = the production rate for GV behind the ars2 promoter
  • βR = the production rate for ArsR behind a constitutive promoter
  • βM = the production rate for MBPArsR behind a constitutive promoter
  • βF = the production rate for fMT behind a constitutive promoter
Vs L Volume of solution (excluding cells).
Vc L Total volume of cells (in solution) (so Vs+Vc is the total volume).

See Chen1997 for the interplay between ArsR and ArsD (the latter has a role similar to ArsR, but we do not treat it, as it is not present in our system).

Quasi steady state

First of all, we assume the concentration of transporters is quite low compared to the concentration of the transported substances. After all, if this were not the case the transporters would act more like "storage" proteins than transporters (note that this can be even more rigorously justified if, for example, GlpFT<<K5). Similarly, there will generally only be a few ars promoters, compared to very many ArsR molecules. This leads to:

As(III)exT ≈ As(III)ex
As(III)inT ≈ As(III)in + ArsRAs + MBPArsRAs + fMTAs
ArsRT ≈ ArsR + ArsRAs

Also, we assume the binding and unbinding of molecules to the transporters occurs on a much finer time-scale than any actual changes to the concentrations inside and outside the cell. Similarly, within the cell we assume diffusion processes are very fast and binding/unbinding of substances is quite fast compared to the production of proteins. This leads us to assume that the following ratios between substances are constantly in equilibrium:

We use the following when grouping the ars promoters:

arsT = ars + ArsRars
ars1 / ars1T = ars2 / ars2T

 ars = ars1 + ars2
 ars = ars1 (1 + ars2T / ars1T)
ars1 = ars / (1 + ars2T / ars1T)
ars1 = ars ars1T / arsT

ars2 = ars ars2T / arsT
As(III)ex : GlpFAs ≈ As(III)ex : 0
GlpF : GlpFAs
ArsB : ArsBAs
As(III)in : ArsRAs : MBPArsRAs : fMTAs : ArsBAs ≈ As(III)in : ArsRAs : MBPArsRAs : fMTAs : 0
ArsR : ArsRAs : 2 ArsRars ≈ ArsR : ArsRAs : 0
ars : ArsRars

To determine what the unknown ratios are we can set the following derivatives to zero (these are the derivatives of the complexes corresponding to the four overlapping regions in the diagram):

0 = (d/dt) GlpFAs = k5on As(III)ex GlpF - (k5off+k6) GlpFAs
0 = (d/dt) ArsBAs = k7on As(III)in ArsB - (k7off+k8) ArsBAs
0 = (d/dt) ArsRars = k3on ArsR² ars - k3off ArsRars
0 = (d/dt) ArsRAs = kRon ArsR As(III)in - kRoff ArsRAs
0 = (d/dt) MBPArsRAs = kMon MBPArsR As(III)in - kMoff MBPArsRAs
0 = (d/dt) fMTAs = kFon fMT As(III)in^nf - kFoff fMTAs

The first two derivates let us determine the ratios between bound and unbound transporters:

0 = (d/dt) GlpFAs = k5on As(III)ex GlpF - (k5off+k6) GlpFAs

k5on As(III)ex GlpF = (k5off+k6) GlpFAs
               GlpF = (k5off+k6)/k5on GlpFAs / As(III)ex
               GlpF = K5 GlpFAs / As(III)ex

                 GlpF : GlpFAs
K5 GlpFAs / As(III)ex : GlpFAs
                   K5 : As(III)ex

0 = (d/dt) ArsBAs = k7on As(III)in ArsB - (k7off+k8) ArsBAs

k7on As(III)in ArsB = (k7off+k8) ArsBAs
               ArsB = (k7off+k8)/k7on ArsBAs / As(III)in
               ArsB = K7 ArsBAs / As(III)in

                 ArsB : ArsBAs
K7 ArsBAs / As(III)in : ArsBAs
                   K7 : As(III)in

The next two differential equations can be used to determine the relative abundances of ArsR and ArsRAs, and ars and ArsRars:

0 = (d/dt) ArsRAs = kRon ArsR As(III)in - kRoff ArsRAs

kRon ArsR As(III)in = kRoff ArsRAs
           ArsRAs = kRon/kRoff ArsR As(III)in
           ArsRAs = ArsR As(III)in / KRd

ArsR : ArsRAs
ArsR : ArsR As(III)in / KRd
KRd  : As(III)in

0 = (d/dt) ArsRars = k3on ArsR² ars - k3off ArsRars

k3on ArsR² ars = k3off ArsRars
       ArsRars = k3on/k3off ArsR² ars
       ArsRars = ArsR² ars / K3d²

ars : ArsRars
ars : ArsR² ars / K3d²
K3d² : ArsR²

For MBPArsR and fMT we find:

0 = (d/dt) MBPArsRAs = kMon MBPArsR As(III)in - kMoff MBPArsRAs

MBPArsR : MBPArsRAs = KMd : As(III)in

0 = (d/dt) fMTAs = kFon fMT As(III)in^nf - kFoff fMTAs

fMT : fMTAs = KFd^nf : As(III)in^nf

And finally the relative abundances of arsenic:

ArsRAs = ArsR As(III)in / KRd

As(III)in : ArsRAs                            : MBPArsRAs                            : fMTAs
As(III)in : ArsRT As(III)in / (KRd+As(III)in) : MBPArsRT As(III)in / (KMd+As(III)in) : fMTT As(III)in^nf / (KFd^nf+As(III)in^nf)
1         : ArsRT / (KRd+As(III)in)           : MBPArsRT / (KMd+As(III)in)           : fMTT As(III)in^(nf-1) / (KFd^nf+As(III)in^nf)

Summarizing:

GlpF : GlpFAs = K5 : As(III)ex
ArsB : ArsBAs = K7 : As(III)in
As(III)in : ArsRAs : MBPArsRAs : fMTAs ≈ 1 : ArsRT / (KRd+As(III)in) : MBPArsRT / (KMd+As(III)in) : fMTT As(III)in^(nf-1) / (KFd^nf+As(III)in^nf)
ars : ArsRars = K3d² : ArsR²
ArsR : ArsRAs ≈ KRd : As(III)in
MBPArsR : MBPArsRAs = KMd : As(III)in
fMT : fMTAs = KFd^nf : As(III)in^nf

Now we can look at the differential equations for the totals of ArsB (so ArsBT=ArsB+ArsBAs), ArsR, As(III)in and As(III)ex (GlpFT and arsT are assumed to be constant):

(d/dt) As(III)exT = (d/dt) As(III)ex + (d/dt) GlpFAs
                  = - (d/dt) GlpFAs - k6 GlpFAs + (Vc/Vs) k8 ArsBAs + (d/dt) GlpFAs
                  = (Vc/Vs) k8 ArsBAs - k6 GlpFAs
                  = (Vc/Vs) k8 ArsBAs - (Vc/Vs) v5 GlpFAs / GlpFT
                  = (Vc/Vs) k8 ArsBAs - (Vc/Vs) v5 As(III)ex / (K5+As(III)ex)
                  = (Vc/Vs) k8 ArsBAs - (Vc/Vs) v5 As(III)exT / (K5+As(III)exT)
(d/dt) ArsBT = (d/dt) ArsB + (d/dt) ArsBAs
             = - (d/dt) ArsBAs + βB ars1 - ln(2)/τB ArsB + (d/dt) ArsBAs
             = βB ars1 - ln(2)/τB ArsB
(d/dt) As(III)inT = -(Vs/Vc) (d/dt) As(III)exT
                  = v5 As(III)exT / (K5+As(III)exT) - k8 ArsBT As(III)in / (K7+As(III)in)
(d/dt) ArsRT = (d/dt) ArsR + (d/dt) ArsRAs + 2 (d/dt) ArsRars
             = βRN ars1 + βR proR - (ln(2)/τR) ArsR - (d/dt) ArsRAs - 2 (d/dt) ArsRars + (d/dt) ArsRAs + 2 (d/dt) ArsRars
             = βRN ars1 + βR proR - (ln(2)/τR) ArsR
(d/dt) MBPArsRT = (d/dt) MBPArsR + (d/dt) MBPArsRAs
                = βM proM - (ln(2)/τM) MBPArsR
(d/dt) fMTT = (d/dt) fMT + (d/dt) fMTAs
            = βF proF - (ln(2)/τF) fMT

Steady state

By looking at the steady state of the system we can say something about its long-term behaviour. This also makes it easier to analyze relations between variables. To derive the steady state solution we take the quasi steady state solution and simplify it further by setting additional derivatives to zero:

0 = (d/dt) ArsBT = βB ars1 - ln(2)/τB ArsB
0 = (d/dt) As(III)inT = v5 As(III)exT / (K5+As(III)exT) - k8 ArsBAs
0 = (d/dt) ArsRT = βRN ars1 + βR pro - (ln(2)/τR) ArsR
0 = (d/dt) MBPArsRT = βM proM - (ln(2)/τM) MBPArsR
0 = (d/dt) fMTT = βF proF - (ln(2)/τF) fMT
0 = (d/dt) GV = βG ars2 - ln(2)/τG GV

This directly leads to:

   0 = βB ars1 - ln(2)/τB ArsB
ArsB = βB (τB/ln(2)) ars1
ArsB = βB (τB/ln(2)) ars1T K3d²/(K3d²+ArsR²)

      0 = βM proM - (ln(2)/τM) MBPArsR
MBPArsR = βM (τM/ln(2)) proM

  0 = βF proF - (ln(2)/τF) fMT
fMT = βF (τF/ln(2)) proF

 0 = βG ars2 - ln(2)/τG GV
GV = βG (τB/ln(2)) ars2
GV = βG (τB/ln(2)) ars2T K3d²/(K3d²+ArsR²)
For the intra- and extracellular concentrations we can find the following equation, giving a maximum for As(III)in of K7 v5/(k8 ArsB) (as As(III)exT cannot be negative)
Conveniently the function x/(c-x) is non-negative and non-decreasing for x∈[0,c⟩.
 i 
:
         0 = v5 As(III)exT / (K5+As(III)exT) - k8 ArsBAs
         0 = v5 As(III)exT / (K5+As(III)exT) - k8 ArsB As(III)in / K7
         0 = v5 As(III)exT - k8 ArsB As(III)in / K7 (K5+As(III)exT)
         0 = v5 As(III)exT - k8 ArsB As(III)in As(III)exT / K7 - k8 ArsB As(III)in K5 / K7
         0 = As(III)exT (v5 - k8 ArsB As(III)in / K7) - k8 ArsB As(III)in K5 / K7
As(III)exT = k8 ArsB As(III)in K5 / (v5 K7 - k8 ArsB As(III)in)
As(III)exT = K5 As(III)in / (K7 v5/(k8 ArsB) - As(III)in)

As we can safely assume arsenic neither disappears into nothingness nor appears from nothingness, we can use this to derive (As(III)T is the total amount of arsenic):

As(III)inT = As(III)in (1 + ArsR/KRd + MBPArsR/KMd + fMT As(III)in^(nf-1)/KFd^nf)

As(III)T = Vs As(III)exT + Vc As(III)inT
0 = Vs As(III)exT + Vc As(III)inT - As(III)T
0 = Vs K5 As(III)in / (K7 v5/(k8 ArsB) - As(III)in) + Vc As(III)in (1 + ArsR/KRd + MBPArsR/KMd + fMT As(III)in^(nf-1)/KFd^nf) - As(III)T

As the function on the right-hand side is non-decreasing for As(III)in∈[0,K7 v5/(k8 ArsB)⟩ it at most has one zero on this interval (and it has one, as it starts at a negative value and gets arbitrarily large as As(III)in approaches the end of its range). So this zero can safely be found using any number of numerical methods.

Finally, for ArsR we can find the following third-order equation:

0 = βRN ars1 + βR pro - (ln(2)/τR) ArsR
0 = βRN ars1T K3d²/(K3d²+ArsR²) + βR pro - (ln(2)/τR) ArsR
0 = βRN ars1T K3d² + βR pro (K3d²+ArsR²) - (ln(2)/τR) ArsR (K3d²+ArsR²)
0 = βRN ars1T K3d² + βR pro K3d² + βR pro ArsR² - (ln(2)/τR) ArsR K3d² - (ln(2)/τR) ArsR³
0 = (βRN ars1T + βR pro) K3d² - (ln(2)/τR) K3d² ArsR + βR pro ArsR² - (ln(2)/τR) ArsR³
0 = (βRN ars1T + βR pro) (τR/ln(2)) K3d² - K3d² ArsR + βR (τR/ln(2)) pro ArsR² - ArsR³

According to Mathematica's solution of Reduce[eq && K3d > 0 && arsT >= 0 && pro >= 0 && β1 > 0 && β3 > 0 && τR > 0, ArsR, Reals] (where eq is the equation shown above) there is only one real solution (examining the discriminant of eq confirms this), so we can solve the equation safely using Newton's (or Halley's) method.


Variables
Variable Description
Extracellular
As(III)ex
Membrane (all naturally occurring, but we plan to bring GlpF to overexpression)
GlpFconcentration w.r.t. the exterior of the cell
GlpFAsconcentration w.r.t. the exterior of the cell
ArsBconcentration w.r.t. the interior of the cell
ArsBAs concentration w.r.t. the interior of the cell
Intracellular (ars2, pro and GV are introduced)
As(III)ex
ars1 concentration of unbound promoters naturally occurring in E. coli
ars2 concentration of unbound promoters in front of gas vesicle genes
proR concentration of constitutive promoters in front of arsR
proM concentration of constitutive promoters in front of mbp-arsR
proF concentration of constitutive promoters in front of fMT
ArsR
ArsR binds to ars to repress production of the genes they regulate, and binds to As(III) to make it less of a problem for the cell.
 i 
ArsRAs bound to As(III)
ArsRars1 bound to ars1
ArsRars2 bound to ars2
MBPArsR
A fusion of maltose binding protein and ArsR.
 i 
MBPArsRAs bound to As(III)
fMT
fMTAs bound to multiple As(III)
ArsRAs bound to As(III)
GV concentration of gas vesicles