Team:Chalmers-Gothenburg/Model

Chalmers Gothenburg iGEM 2017

Introduction

The modeling part of our project was to simulate our biosensor in silico to gain insight in how the output might look, and what may be changed or modified in it to produce a better result.

One of the main parts of our biosensor uses the yeast pheromone pathway with custom GPCRs that will make the cells sense volatile organic compounds (VOCs) that are lung cancer biomarkers. These receptors, when activated by the VOCs in the breath of patients will trigger a proteolytic cascade, the MAPK cascade, in the pheromone pathway which ultimately leads to increased levels of a transcription factor called Ste12 [1]. Ste12 acts as a “first checkpoint” in our biosensor as its values will correlate with the downstream processes. This transcription factor will then promote expression of the protein Fus1, whose function involves coordinating cell polarization and fusion in order to make the yeast mate with the opposing mating type [3]. Therefore, the promoter region (PFUS1) was used as a way to signal when the pheromone pathway is activated.

In parallel to the normal pheromone pathway, a construct containing Cre recombinase with PFUS1 in front of it was added, which causes it to be co-expressed by pheromone pathway activation. Cre in turn flips the constitutive promoter PTEF by recombining two loxP-sites surrounding the promoter. This leads to expression of Cas9 in one mating type and a gRNA in the other. When the cells of opposing mating types then mate to form a diploid, the Cas9 and gRNA will bind to each other. This will in turn cause the Cas9-gRNA complex to cut the ADE2 gene in yeast and turn the cells red.

Modeling and results

To model this system, we decided to use a set of coupled ordinary differential equations (ODEs) that correspond to the different species involved in the pheromone pathway. The model relies upon the work of Kofahl and Klipp [1], whose model of this pathway has been used extensively by many iGEM teams throughout the years and has produced good results [4,5,6]. The model is very detailed and involves 36 different species with almost 50 reactions between them [1]. The modeling part of the project uses all of the ODEs and parameter from their original paper, along with other ODEs stated in Table 3.

Figure 1.The reaction network used in the model by Kofahl and Klipp [1].

As seen in the bottom right corner of Figure 1, the cascade ends with an activation of Ste12 which makes it possible for it to bind into DNA and begin transcriptional activation. This was therefore used as a first step into validating the function of the biosensor. The concentration of Ste12 should be higher in lung cancer patients compared to healthy people, due to the fact that they produce a higher amount of VOCs that can activate the pheromone pathway [7]. To check this, data on different VOC levels was used, both in in lung cancer patients and compared to a healthy control group and a group consisting of smokers [7]. The data of the median level of each group is presented in Table 1 and consists only of n-octanal levels. The choice of only studying the n-octanal levels is due to the fact that no data (kinetic rates v1 to v5 in Kofahl and Klipp’s model, see Figure 1) is available for the two receptors used in this project, thus the 2-butanone levels were omitted. The original yeast parameters of these reactions were used together with the n-octanal concentrations from the three patient groups.

Table 1. Median n-octanal level values for different patient groups [7].
Patient group (n=12) Controls Smokers Lung cancer patients
Median n-octanal values (nmol/L) 0,011 0,006 0,052

Given these assumptions, the simulation for the Ste12 concentration can be run on the three different patient groups, which gives the result presented in Figure 2a. Even though the n-octanal level is about five times higher in lung cancer patients compared to the controls, this only leads to a 50% increase in peak Ste12 concentration. To study whether this affected expression of Cre recombinase another ODE was added into the system, see Table 2, that models binding and gene expression of Ste12 onto PFUS1.

Figure 2. Time course plots for the different species in our model.
Table 2. ODEs describing the temporal behavior of different species in the model. The parameter values are either randomly chosen or from literature.
ODEs added to Kofahl & Klipp's model Parameter values
Cre/Fus1 production: [2,8]
Cas9/gRNA formation: Arbitrary value
Diploid formation: Arbitrary value
Haploid mating: Arbitrary value
Cas9:gRNA complex formation: Arbitrary value

The binding kinetics are assumed to be of Hill type with values found in the literature [2,8]. The initial concentrations of Fus1 and Cre recombinase are assumed to be zero. The time courses for Cre/Fus1 are shown in Figure 2b. Notice that we now have about double the Fus1/Cre recombinase amounts for the lung cancer patients compared to the controls. The Cre recombinase will flip a constitutive promoter PTEF1 flanked between two LoxP sites. This will in turn activate transcription of the gRNA for ADE2 disruption in one mating type and transcription of the protein Cas9 in the other. Since the production rate of PTEF1 is unknown, we set it to an arbitrary value of 1mmol/min/L. Setting an arbitrary production rate does not allow us to make any quantitative predictions about the gRNA/Cas9 levels. However, due to the fact that the original yeast parameters are used for the receptors, quantitative information is lost anyway. Although we can still estimate relative concentrations of these species between our patient groups.

The mating rate in this model is assumed to be proportional to the amount of haploid cells in the culture and the Fus1 concentration. Given an initial 1000 haploids, we produce the corresponding amount of diploids for each patient group, shown in Figure 2d. For the concentration of activated Cas9-gRNA complex per diploid cell we assumed it to be proportional to the amount of haploids, the gRNA concentration in each cell and the Cas9 concentration in each cell. The result from this simulation is shown in Figure 2e. After 24 hours there is about a 50% increase in Cas9:gRNA dimers compared with the controls. This increase will in turn promote deletion of the ADE2 gene, causing the cell to turn red. Although this can be modeled, since there were not any parameters at all for the ADE2 pathway in the literature, we decided to only model up until the Cas9:gRNA dimer formation.

From this we can conclude that it might be possible to differentiate between lung cancer patients and controls since we get a 50% higher level of Cas9:dimers after 24 hours.

MATLAB code

The MATLAB code used in this project can be found here:

pheromonepathway.m

This is the main .m-file that should be run. Available here.

f.m

This function is called by the ODE solver and contains all of the ODEs. Available here.

References

  • [1]   Kofahl B,Klipp E Modelling the dynamics of the yeast pheromone pathway. Yeast 2004; 21: 831–850.
  • [2]   Dyjack N, Azeredo-Tseng C, Yildirim N Mathematical modeling reveals differential regulation of MAPK activity by phosphatase proteins in the yeast pheromone response pathway. Mol. BioSyst. Mol. BioSyst 2017; 13: 1323-1335.
  • [3]   Yeastgenome.org. (2017). FUS1 | SGD. [online] Available at: https://www.yeastgenome.org/locus/S000000532 [Accessed 17 Sep. 2017].
  • [4]   2012.igem.org. (2017). Team:TU-Delft/Modeling/SingleCellModel - 2012.igem.org. [online] Available at: https://2012.igem.org/Team:TU-Delft/Modeling/SingleCellModel [Accessed 17 Sep. 2017].
  • [5]   2015.igem.org. (2017). Team:CGU Taiwan/Modeling - 2015.igem.org. [online] Available at: https://2015.igem.org/Team:CGU_Taiwan/Modeling [Accessed 17 Sep. 2017].
  • [6]  2015.igem.org. (2017). Team:Chalmers-Gothenburg/Detection - 2015.igem.org. [online] Available at: https://2015.igem.org/Team:Chalmers-Gothenburg/Detection [Accessed 17 Sep. 2017].
  • [7]   Fuchs P, Loeseken C, Schubert J. K. , Miekisch W Breath gas aldehydes as biomarkers of lung cancer. International Journal of Cancer 2010.
  • [8]   Groves B, Khakhar A, Nadel C. M. , Gardner R. G. ,Seelig G, Rewiring MAP kinases in Saccharomyces cerevisiae to regulate novel targets through ubiquitination. eLife 2016.