Team:UNAM Genomics Mexico/Modeling/Heavy Metal AND

From 2012.igem.org

Revision as of 07:56, 25 October 2012 by Abieltega (Talk | contribs)


UNAM-Genomics_Mexico


AND's Results


Metal AND


Our system consists of a hybrid promoter which contains binding sites for two transcription factors that act as repressors, CzrA and ArsR. In the presence of Met1, CzrA losses affinity for its binding site.

Similarly, if there is Met2 in the cell, ArsR also dissociates itself from its binding site. Taking into account these facts, unless Met1 and Met2 are present in the culture, we can say that the genes under this hybrid promoter will be repressed, just like a AND gate should function!

UNAM GENOMICS MEXICO AND bf hm.png

Repressors

UGM tv mtf.png

Inputs

UGM tv mi.png


Kappa



There is no metal catabolism, so metals just reach an equilibrium and non-toxic concentration within the cell through exporting it out of the cytosol, maintaining the concentration constant. This makes the behavior of a single burst or a continuous dose similar, because the metal levels in just one dose will be constant since metals won’t be degraded.


6000metales.png Single an conti.png


Transfer Function



We were interested in modelling the output of our logic gates based on the input. This is called the "Transfer function". Imagine this process like a black box that will give you the dynamic concentration through time of whatever is downstream of the "AND" when you feed it with input data. It doesn't matter if it’s a single input burst or a continuous input.

Our logic AND gate is fully dependent on the intracellular concentrations of the inducers (meaning metal ions).

Then, to accomplish our duty in model, we use the "Transfer Function", for which it is necessary to take into account the regulation and dynamics of the endogenous B. subtilis intake-efflux system that controls the intracellular concentrations of Heavy Metals and Sugars. All of this is going to be the "Black Box".

So, our first task was to reconstruct the regulatory network of B. subtilis for the intake-efflux systems. The regulation data was retrieved from several papers and databases like http://bsubcyc.org/


UNAM GENOMICS MEXICO HM AND.png

Heavy metal AND: Some metals are introduced to the cell by MntH. These metals, in the event of Met 1 being in the cell, are expelled from the cell by a transporter of the CadA operon and those of Met2 set by ArsB. The master regulators of each set are CzrA and ArsR. They control the production of the transporters by repressing them in the absence of metals. Both, CzrA and ArsR repress the expression of LasR or P4 and RFP under the hybrid promoter of CadA+ArsR.


Metal homeostasis in the cell involves an effective response in order to withstand possible stress conditions in which high concentrations of heavy metals are involved. Both approaches, ODEs with Matlab and the stochastic one with kappa, give us similar results. First, we were able to see that the AND response is not as sensitive to higher concentrations. For the AND to have a good performance, low concentrations of Met 1 and 2 are required. As we said, we infer that metal response is a fast process. The accumulation of the output protein shows a sudden increase in production resulting from a very small change in metal concentrations. Another remarkable fact is about the role each TF (CzrA and ArsR) has. Since ArsR represses itself, its production has a succession of spikes. Nevertheless, CzrA increases its production if the metal concentration within the cell also increases, but it doesn’t have a drop in production because their decay depends only on the dilution and degradation of the protein. That way, the effects of CzrA are more notable than ArsR. This is reflected in a greater contribution of Met1 in the regulation of the AND promoter than that of Met2.


Metal burst 2.png
UNAM GENOMICS MEXICO Metal burst 3D 2.png


Deterministic model



Our team made both, a deterministic and a stochastic approach. Here, we are going to describe the deterministic model.
General assumptions:
●This model's scope is single cell. While the regulation steps were modeled through Hill equations, other reactions were treated under Michaelis-Menten assumptions, notably isomerizations and transport processes (i.e. those that involved energy consumption).

●The maximal transcription rate for the system species was obtained by dividing the maximum number of nucleotides processed by the RNA polymerase (i.e. maximum polymerase activity), by the length of the nucleotide sequence.

●Translational rates were obtained with the maximum number of amino acids added by the ribosome divided by the length of the protein in amino acids.

●The protein degradation rate is considered as equal for all the protein species in the model.

●The degradation rate for mRNAs is equal for all the different species of mRNAs.

●The different binding sites of the same transcription factor have the same affinity to the TF.

●For each mRNA species, its concentration will be the sum of the production as affected by its respective promoter and the transcription factors that regulate it, minus the degradation of the mRNA at that time.

●For each protein, the change in its concentration depends on the amount of protein produced by translation minus the degradation rate of the protein.

●The Hill coefficients tend to be 2.

●The parameters of concentration were expressed in terms of molecules per cell(Vol of a B.subtilis cell ~10^-14) and the time units were converted sec.

In Bacillus subtilis, many metal ions, both essential, and others like arsenic, are introduced in the cell by the ion pump MntH. We make two assumptions, first, that the used metals were imported byMntH and second, that all the metals introduced by MntH (especially those that we use for our AND gate) have the same intake rate. Met1 is pumped out of the cell by a transporter of the CadA operon (most of them require ATP hydrolysis), and Met2 by ArsB (Moore C, 2005[21]). In that way, each metal’s concentration in the cell was defined by Michaelis-Menten equations as the amount of introduced metal by each MntH molecule in the cell, minus the amount of exported metal by its respective transporter, assuming that all the metals imported or exported by the same transporter had the same affinity to this transporter and had the same chance to pass into the cell.

UNAM GENOMICS MEXICO O1.png

The transcription factors which are involved in the AND processing are CzrA and ArsR. ArsR is repressed by itself when there are no Met2 ions. We use a Hill function to express the mRNA production. The function shows the effect of the binding of ArsR in its promoter. This way, the repression effect will be affected by the affinity of ArsR to the binding site and the concentration of ArsR. About CzrA, there are studies that suggest that the total amount of the protein in the cell is correlated with the levels of the metals that it tends to expel (Moore C, 2005 [9]). To describe this phenomenon, we defined the production of the mRNA of CzrA as the production due to intracellular Met1 mediated induction into the cell through a Michaelis-Menten relationship.

UGM O2.png


Finally, the proportion of active transcription factors (unbound to their respective metals) are defined with a hill function as the fraction of the total protein plus the relationship between the affinity constant of the binding and the amount of metal. The higher amount of metal or the affinity of it to the binding with TF corresponds with the lower amounts of active TF. The transporters for Met1 and Met2 are also under the regulation of ArsR and CzrA; the CadA operon is repressed by CzrA, and ArsB is repressed by ArsR.

UGM O3.png

The output of the AND operation is the production of one of two reporters, P4 or LasR respectively. The production of the reporter mRNA depends on the status of its promoter, which has a state space of being bound to ArsR, to CzrA, to both, or bound to none. This approach was made supposing that the binding-mediated repression of each TF is independent, and therefore the cumulative intensity of repression by ArsR and CzrA is the product of their independent repressions.

UGM O4.png

Conclusions

After making the model just for our AND, me coupled the endogenous Bacillus subtilis system to our predictions. After this, we found that the predictions given by both assumptions were not the same, proving that the endogenous system was, in fact, an important matter to consider, so we considered it.

Obviously, as the model that did not take into account the endogenous system was farther away from the actual system, we opted for using the model that did take it into account, making our transfer function a more realistic representation of the system, in spite of the fact that we could not find all the parameters.

Some of these parameters, unfortunately, were the diffusion within the nanotubes, and the amount, size, form, frequency of creation, or any other useful data of nanotubes you could think of (luckily for us, we could not think of that many). Even though we had to face all these problems, with the little light shone by people before us onto the spectral essence of the nanotubes, it was possible for us to say that our system can, in theory, simulate Boolean operations.