Binding of R2 to C

SCBs (C₂) bind to ScbR homo-dimer (R₂) and inactivate its repressing activity.

Chemical equation

The exact mechanism is still unclear, however in our model we assumed that two SCBs bind to the ScbR homo-dimer.

$C_{2} + R_{2} \rightleftharpoons C_{2}-R_{2}$

Rate equation

$r= \frac{k^{-}_{4}}{K_{d4}}\cdot [C_{2}]\cdot [R_{2}] - k^{-}_{4}\cdot [C_{2}-R_{2}]$

Parameters

The parameters of this reaction are the dissociation constant for binding of SCB to ScbR ($K_{d4}$) and the dissociation rate for binding of SCB to ScbR ($k^{-}_{4}$). ScbR is a member of the TetR family of repressors, named after the member of this group which is the most completely characterized, the TetR repressor protein. TetR binds to the operator tetO, repressing its own expression and that of the efflux determinant tetA. However, [MgTc]⁺ binds to TetR and thus the affinity of the later for the operator tetO is 9-fold reduced. This procedure is similar to ScbR binding to O_R and O_A and repressing its own expression and the expression of ScbA, while binding to SCBs reduces its affinity for the two operators. Therefore parameter values were derived from published data on the TetR-[MgTc]⁺ and TetR-Tc (without Mg⁺) interactions.

Name	Value	Units	Value in previous GBL models	Remarks-Reference
$K_{d4}$	$10^{-2} - 10^{4}$	$nM$	$0.083 nM^{-1} s^{-1}$ (Range tested: $10^{-7}-0.1 nM^{-1} s^{-1}$) (Bistability range: $0.077-0.17 nM^{-1} s^{-1}$ and $0.0042-0.253 nM^{-1} s^{-1}$)	According to Hillen et al. and Orth et al. the association constant for TetR-Tc binding is in the ~$10^{9} M^{-1}$ range in presence of Mg2+, as determined by in vitro measurements, therefore a Kd= ~$10^{-9} M = ~ 1 nM$ Orth et al. 2011 Hillen et al. 2010 On the other hand, Kedracka-Krok et al. conducted stopped-flow measurements using TetR overproduced in Escherichia coli strain RB 791. They consequently reported an association constant Ka=$0.96 \cdot 10^{5} M^{-1}$ in absence of Mg2+ (therefore a Kd=$1.04 \cdot 10^{-5} M = 1.04 \cdot 10^{4} nM$) and an association constant Ka=$6.3 \cdot 10^{6} M^{-1}$ for binding of TetR to [Tc-Mg]+ (therefore a Kd=$0.16 \cdot 10^{-6} M = 160 nM$). Kedracka-Krok et al. 2005 Kedracka-Krok et al. 2005 Finally, Schubert et al. reported Mg2+ independent KAs in the ~$10^{7} M^{-1} (10^{2} nM)$ range and Mg2+ dependent KAs in the ~$10^{11} M^{-1} (10^{-2} nM)$ range. The data was derived from in vitro and in vivo measurements in E. coli K12, strains DH5a and WH207. Schubert et al. 2004
$k^{-}_{4}$	$3 \cdot 10^{-4}-126$	$min^{-1}$	$630 s^{-1}$ (Range tested: $0-10^{3} s^{-1}$) (Bistability range: $460-630 s^{-1}$ and $8.5-195 s^{-1}$)	According to Kedracka-Krok et al. the unbinding rate for the TetR-Tc complex is $2.1 s^{-1} = 126 min^{-1}$ in absence of Mg2+ and $2.2 \cdot 10^{-2} s^{-1} = 1.32 min^{-1}$ in presence of Mg2+ (see figure above). However, Schubert et al. reported a dissociation rate of $5 \cdot 10^{-6} s^{-1} (3 \cdot 10^{-4} min^{-1})$ (see table above). In order to quantify the dependency between the parameters, a distribution for $k_{on4}$ needs to be defined. From the information reported in the tables above by Kedracka-Krok et al. and by Schubert et al., a range of values can be derived for this parameter as well. According to the experimental data, the value of $k_{on4}$ was found to be between $1.4 \cdot 10^{5}- 8 \cdot 10^{6} M^{-1} s^{-1} (0.0084-0.48 nM^{-1} min^{-1})$. These values will be used to define the corresponding probability distribution.

Parameters with uncertainty

When deciding how to describe the uncertainty for each parameter there are a few points to be taken into consideration. Firstly, the values reported in literature are spread in a relatively large range and correspond to TetR and TetR mutant proteins. Additionally, most of the values were acquired by in vitro testing. This means that there might be a notable difference between actual parameter values and the ones reported in literature. These facts influence the quantification of the parameter uncertainty and therefore the shape of the corresponding distributions.

With regards to the $K_{d4}$ the value that is mostly reported in different publications is $1 nM$, therefore we put the weight of the distribution in the range $1-10 nM$ and we consider as least likely the larger values. Therefore, the mode of the log-normal distribution is set to $5 nM$ and the confidence interval factor is $1000$. Thus the range where 95.45% of the values are found is between $5 \cdot 10^{-3}$ and $5 \cdot 10^{3} nM$.

Since the mode chosen for $K_{d4}$ is set in the $10^{-9} M$ range and the two parameters are interconnected, the smaller values reported by Schubert et al. ($10^{-4} min^{-1}$) are considered least likely for $k^{-}_{4}$. Therefore, the mode of the log-normal distribution for $k^{-}_{4}$ is set to $1.3 min^{-1}$ and the confidence interval factor is $100$. In this way the range where 95.45% of the values are found is between 0.013 and 130 $min^{-1}$.

Finally, the probability distribution for $k_{on4}$ is defined accordingly, in order to allow the exploration of the full range of the values retrieved from literature. Therefore, the mode is set to $0.06 nM^{-1} min^{-1}$ and the confidence interval factor is $7$. In this way the range where 95.45% of the values are found is between 0.0086 and 0.42 $nM^{-1} min^{-1}$.

The probability distributions for the three parameters, adjusted accordingly in order to reflect the above values, are the following:

500px 500px 500px

The location and scale parameters of the distributions are:

Parameter	μ	σ
$K_{d4}$	Failed to parse (Cannot store math image on filesystem.): 5.27937	Failed to parse (Cannot store math image on filesystem.): 1.9157
$k^{-}_{4}$	$2.3977$	$1.46129$