Difference between revisions of "Binding of R2 to C"

From ISMOC
Jump to: navigation, search
(Parameters with uncertainty)
(Parameters with uncertainty)
Line 79: Line 79:
 
The probability distributions for the three parameters, adjusted accordingly in order to reflect the above values, are the following:
 
The probability distributions for the three parameters, adjusted accordingly in order to reflect the above values, are the following:
  
[[Image:Kd4.png|500px]] [[Image:K4.png|500px]] [[Image:Kon4.png|500px]]
+
[[Image:KD4u.png|500px]] [[Image:K4u.png|500px]] [[Image:Kon4u.png|500px]]
  
 
The correlation matrix which is necessary to define the relationship between the two marginal distributions (<math>K_{d4}</math>,<math>k^{-}_{4}</math>) of the bivariate system is derived by employing random values generated by the two distributions.  
 
The correlation matrix which is necessary to define the relationship between the two marginal distributions (<math>K_{d4}</math>,<math>k^{-}_{4}</math>) of the bivariate system is derived by employing random values generated by the two distributions.  

Revision as of 18:33, 18 January 2017

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

Go back to overview
About this image

Chemical equation

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

2C + 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 OR and OA 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 [1] [2] Remarks-Reference
K_{d4} 10^{-2} - 10^{4} [3] [4] [5] [6] [7] 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}[1]

and 0.0042-0.253 nM^{-1} s^{-1}[2])

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[5]
  • Hillen et al. 2010[6]

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[3]
  • Kedracka-Krok et al. 2005[3]

Finally, Schubert et al. reported Mg2+ independent KAs of 1.95 \pm 0.11 \cdot 10^{7} M^{-1} and 10 \pm 0.31 \cdot 10^{7} M^{-1} (K_d=51.28 \pm 2.89 nM and K_d=10 \pm 0.31 nM) and Mg2+ dependent KAs of 2.3 \pm 0.21 \cdot 10^{11} M^{-1} and 9.11 \pm 0.3 \cdot 10^{11} M^{-1} (K_d=0.435 \pm 0.04 nM and K_d=0.11 \pm 0.0036 nM). The data was derived from in vitro and in vivo measurements in E. coli K12, strains DH5a and WH207.

Schubert et al. 2004[7]
k^{-}_{4} 3 \cdot 10^{-4}-126 [3] [7] min^{-1} 630 s^{-1}

(Range tested: 0-10^{3} s^{-1})

(Bistability range: 460-630 s^{-1}[1]

and 8.5-195 s^{-1}[2])

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 dissociation rates of 5 \pm 0.5 \cdot 10^{-6} s^{-1} (3 \pm 0.3 \cdot 10^{-4} min^{-1}) and 7 \pm 0.4 \cdot 10^{-6} s^{-1} (4.2 \pm 0.24 \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. Therefore, by assigning the appropriate weights to the parameter values and using the method described here, the appropriate probability distributions were designed.

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  9.6 nM and the confidence interval factor is  37232 . Thus the range where 95.45% of the values are found is between  2.6 \cdot 10^{-4} and  3.6 \cdot 10^{5} nM. (Note: This distribution represents our initial beliefs about the system but is redefined as explained below, in order to account for the thermodynamic consistency of the system.)

The smaller values reported by Schubert et al. ( 10^{-4} min^{-1}) are also considered least likely for k^{-}_{4}. The mode of the log-normal distribution for k^{-}_{4} is calculated as 1.23 min^{-1} and the confidence interval factor is 24779. In this way the range where 95.45% of the values are found is between  4.9 \cdot 10^{-5} and  3.05 \cdot 10^{4} 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.071 nM^{-1} min^{-1} and the confidence interval factor is 19.4. In this way the range where 95.45% of the values are found is between 0.0037 and 1.3802 nM^{-1} min^{-1}.

Since the three parameters are interdependent, thermodynamic consistency also needs to be taken into account. This is achieved by creating a bivariate system as described here. Since K_{d4} is the parameter with the largest geometric coefficient of variation, this is set as the dependent parameter as per: K_{d4}=\frac{k^{-}_{4}}{k_{on4}}, and an updated probability distribution is defined. The location and scale parameters of K_{d4} (μ=7.684 and σ=2.6711) were calculated from those of k_{on4} and k^{-}_{4}. The probability distributions for the three parameters, adjusted accordingly in order to reflect the above values, are the following:

KD4u.png K4u.png Kon4u.png

The correlation matrix which is necessary to define the relationship between the two marginal distributions (K_{d4},k^{-}_{4}) of the bivariate system is derived by employing random values generated by the two distributions.

The final parameters of the distributions of the multivariate system are:

Parameter μ σ Correlation matrix
k_{on4} -1.4918 1.0729 N/A
K_{d4} 7.684 2.6711 \begin{pmatrix} 1 & 0.5955 \\
0.5955  & 1 \end{pmatrix}
k^{-}_{4} 6.1922 2.4461

The multivariate system of the normal distributions (ln(k_{d4}) and ln(k^{-}_{4})) and the resulting samples of values are presented in the following figure:

Error creating thumbnail: Unable to save thumbnail to destination

In this way, a system of distributions is created where each distribution is described and constrained by the other two. Therefore, the parameters will be sampled by the two marginal distributions in a way consistent with our beliefs and with the relevant thermodynamic constraints.

References