Team:TU Delft-Leiden/Modeling

Modeling Overview

We modeled all three different modules of our project, namely the conductive curli module, the extracellular electron transport (EET) module and the landmine detection module. In order to achieve this, we had to use all kinds of different modeling methods. One of these methods is based on the percolation theory and is used in the conductive curli modeling module. To the best of our knowledge, this theory has never before been used to model a biological system. Because we use this theorem and for other reasons, the model for the conductive curli module is a stochastic model. The model for the EET and landmine modules are deterministic models.
We used Matlab for most of the calculations; the scripts we made can be found in the code repository. Also, we had great interactions with the Life Science and Microfluidics departments, which for the conductive curli module can be read here, for the EET module can be read here and for the landmine detection module can be read here.

Curli Module

The goal of our project for the conductive curli module is to produce a biosensor that consists of E. coli that are able to build a conductive biofilm, induced by any promoter, in our case a promoter that gets activated in the presence of DNT/TNT, see the gadget section of our wiki and the extracellular electron transport (EET) module. The biofilm consists of curli containing His-tags that can connect to gold nanoparticles, see the conductive curli module. When the curli density is sufficiently high, a dense network of connected curli fibrils is present around the cells. Further increasing the amount of curli results in a conductive pathway connecting the cells, thereby forming conductive clusters. Increasing the amount of curli even further, sufficiently curli fibrils are present to have a cluster that connects the two electrodes and thus have a conducting system.
The goal of the modeling of the curli module is to prove that our biosensor system works as expected and to capture the dynamics of our system. So, we want to answer the question: "Does a conductive path between the two electrodes arise at a certain point in time and at which time does this happen?" However, we not only want to answer the question if our system works as expected qualitatively, but we also want to make quantitative predictions about the resistance between the two electrodes of our system in time.

To capture the dynamics of our system, we have implemented a three-layered model, consisting of the gene level layer, the cell level layer and the colony level layer:

• At the gene level, we calculate the curli subunits production rates and curli subunit growth that will be used in the cell level.
• At the cell level, we use these production and growth rates to calculate the curli growth in time, which we will use at the colony level.
• At the colony level layer, we determine if our system works as expected, ie. determine if a conductive path between the two electrodes arises at a certain point in time and at which time this happens. We also determine the change of the resistance between the two electrodes of our system in time.

A figure of our three-layered model is displayed below.

conclusions has to be written

EET Module

Flux Balance Analysis

In the wet lab, we integrated the Extracellular Electron Transport (EET) module of S. oneidensis into E. coli, see the extracellular electron transport (EET) module. For the modeling of the EET module, we wanted at first to gain insight in the consequences of the integration of the EET module into E. coli. To achieve this, we simulated the cell metabolism of E. coli including the EET module using the Flux Balance Analysis (FBA) method. Our goal was to investigate the carbon metabolism providing the electrons for the EET module [reference]. Also, we want the EET pathway used by the cells in order to have a measurable electrical signal for our biosensor, see the gadget section of our wiki.

From the FBA method, we conclude that in aerobic conditions the cell does not use the EET pathway, but oxygen gets reduced instead, as it is a stronger oxidizing agent. However, in anaerobic growth the cell does use the EET pathway to export electrons out of the cell. When the cell is grown on glucose, the growth rate will be higher than when the cell is grown on D-lactate.
We also conclude that in an experimental setting the EET pathway has a higher chance of being used when the cells are grown on D-lactate as the EET pathway is necessary in order for the cells to grow, while when grown on glucose and the EET pathway is turned off (represented by \(0 \ mmol \ (gDW)^{-1} \ hr^{-1}\) (per gram dry weight per hour) maximum EET flux) growth is still possible, see figure 2. From Flux Variability Analysis (FVA) we conclude that for maximum growth for each specific combination of carbon source uptake flux and maximum EET flux, only one possible EET flux is possible for both growth on glucose and growth on D-lactate, namely the EET flux shown in figure 2.

As performing the FBA method while maximizing the EET flux yielded no growth, we wondered if there are pathways possible that would yield growth. So, we performed FVA, the results can be found in figure 3. From this figure, we conclude that when maximizing the EET flux, there are pathways possible that yield growth, as the figure displays the difference between the optimized maximum value and the optimized minimum value of growth. We see that this value is not equal to zero everywhere in the figure, thus there are multiple pathways possible when maximizing EET flux that all yield different values for growth. Note that the EET flux is 2.8 times higher and 2 times higher in comparison to maximizing for growth rate, for glucose and D-lactate as a carbon source, respectively.

From figures 2 and 3 we conclude that there are different regions in which the cell can operate. In an experimental setting, it can be investigated in which region the cell actually operates and if it maximizes its growth rate or its EET flux. To be able to do this, the experimental observed pathway has to be compared to the possible pathways when maximizing the EET flux and to the pathway when maximizing the growth rate. From these regions, it can be deduced if the experimentally observed EET flux and growth rate are carbon source limited or limited by the maximum possible EET flux.
Finally, we investigated an extended model of E. coli metabolism. This model contains, in contrast to the previously used core model, L-lactate as a metabolite. Using the extended model, we found that for glucose and D-lactate as carbon sources, the maximized growth rate agreed quite well to the previous analysis's using the core model. Using L-lactate as a carbon source, we conclude that a steady state solution in which E. coli can grow on L-lactate and use the EET pathway is not possible. A possible way to obtain information about the EET flux when the cells are not in steady state as observed by Goldbeck et al. [2], would be by the use of dynamic flux balance analysis (dFBA), which can also model the dynamics of a system before it reaches steady state [3].

Deterministic Model of EET Complex Assembly

The EET module consists of three proteins: MtrA, a cytochrome on the inside of the outer membrane, MtrB, a β-barrel protein located in the outer membrane, and MtrC, another cytochrome, located on the cell surface. This complex enables the cell to transport electrons from the cytoplasm of the cell to the extracellular environment, see the extracellular electron transport (EET) module.
The assembly of the trans-membrane EET complex depends on many factors other than transcriptional and translational control, as it requires a large amount of post-translational modifications. We set up a simplified model of this assembly process, largely based on section 1.3 of the thesis of Jensen [4]. With the use of deterministic modeling methods, our goal is to predict how many EET complexes are formed under different initial conditions.
In our modeling of the assembly of the EET complex, in addition to the assembly mechanism, we also focus on the apparent reduced cell viability. Jensen [4] proposes two possible explanations for this: the formation of cytosolic aggregates and reduced membrane integrity due to the high amount of trans-membrane protein complexes.

We constructed two different models, one based on an extensive model of EET complex assembly, which we based upon the work of Jensen [4], the other based on a very much simplified model of EET complex assembly, which includes only the most fundamental reactions of the assembly process.
From the extensive model, we concluded that the amount of δ-ALA (and therefore heme) is rate limiting and not the amount of available binding sites. We therefore predict that adding extra δ-ALA to the cells will increase the amount of EET complexes, see figure 4. This effect is also observed by Jensen [4].

Although the extensive model proved to be valuable in the investigation of the mechanism which assembles the EET pathway, it is not suitable for the quantitative prediction of the amount of EET complexes. The most important reason for this is the large number of unknown parameters. Therefore we decided not to aim at enhancing this model, and rather set up a more simplified model.

To sidestep the difficulties experienced with the extensive model of the assembly of the EET complex, we reduced the system to a bare minimum. For this simplified model, we only included the production of MtrCAB, the formation of cytosolic aggregates and the assembly of the EET complex.
The simplified model does not represent all the actual mechanisms of the EET complex assembly process that happens in nature, but it is able to match the experimental data of Goldbeck et al. [2], see figure 5. A maximum at low promoter strength is clearly visible. This corroborates the statement in [2] that maximum promoter strength does not result in maximum EET concentration due to reduced cell viability.

This is a confirmation that our assumptions concerning cell viability might be correct. Therefore, we conclude that reduced cell viability because of the implementation of the EET pathway is the consequence of three molecular processes, namely firstly, the amount of EET complexes reduces the transcription and translation of the MtrCAB proteins due to reduced membrane integrity, secondly, the forming of MtrCAB aggregates and thirdly, the clogging of the secretion system transporting the MtrCAB complexes.

Landmine Module

An important part of our iGEM project is a promoter sensitive to DNT/TNT. We will use two promoters that are sensitive to DNT/TNT, namely ybiJ and ybiFB2A1, in our project, see the landmine detection module. Of these promoters, not much is known other than the fact that they have a DNT/TNT-dependent response curve. Our goal was to find a model which would be able to reproduce the response curves of both promoters. To achieve this, we constructed two different models, both using deterministic modeling methods. One model is based on a simple binding model of DNT to the promoter, the other is based on cooperative binding of DNT to the promoter.
When based on the simple binding model, fits of promoter activation with respect to DNT concentration to the experimental data of [1] did not yield good results. However, when the fits were based on the cooperative binding model, we were able to match the experimental data in [1] really well, see figure 1.
From the fits shown in figure 6, it is clear that both the standard and cooperative activation model fail to describe the experimental data we obtained. This might have several reasons. Firstly, the response ratio of our measurements is very small compared to the measurements in [1]. Secondly, our experimental data suggest a detrimental effect of induction with high DNT concentration. This effect is not described in [1] and we don't know its reason. It is therefore not included in our model. Thirdly, both the data set from [1] and our experimental data set are small compared to the amount of parameters. This makes finding the right model difficult.
To improve the promoter activation model, a larger and more consistent data set has to be obtained. We need to conduct fluorescence measurements over a wide range of DNT concentrations. Besides that, possible detrimental effects of its solvent (acytonitrile) have to be investigated. Also, the induction of the N genes with rhamnose might influence the measurements in an unexpected way. Doing all this is unfortunately not possible in the time span of this iGEM project.

References

[1] S. Yagur-Kroll, S. Belkin et al., “Escherichia Coli bioreporters for the detection of 2,4-dinitrotoluene and 2,4,6-trinitrotoluene”, Appl. Microbiol. Biotechnol. 98, 885-895, 2014.

[2] C.P. Goldbeck, H.M. Jensen et al., “Tuning Promoter Strengths for Improved Synthesis and Function of Electron Conduits in Escherichia coli”, ACS Synth. Biol. 2, 150-159, 2013.

[3] R. Mahadevan, J.S. Edwards & F.J. Doyle, “Dynamic Flux Balanace Analysis of Diauxic Growth in Escherichia coli”, Biophys. J. 83, 1331-1340, 2002.

[4] H.M. Jensen, “Engineering Escherichia coli for molecularly defined electron transfer to metal oxides and electrodes”, PhD Thesis Chemistry UC Berkeley, 2013.