Team:UNOTT/Modelling

Constitutive Gene Expression For Protein and mRNA Expression over Time

Biological insight had told us we need a model with constant gene expression. Investigating models from literature ¹ so see which model would satisfy these conditions, and it was found the constitutive gene expression model was suitable to guide the model.

The first step was to take the general model from literature and apply it in our scenario using the proteins (GFP, ECHP, RFP.)

^{Figure 1} $$ sfGFP \underset{Transcriptin}{\rightarrow} mRNA \underset{Translation}{\rightarrow} sfGFP $$

The equation above describes the process of which the gene undergoes transcription to produce mRNA. The mRNA carries the genetic information copied from the DNA which codes for protein. The expression of protein, can therefore, be measured by the fluorescence which is the desired output of the system.

^{Figure 2} $$ mRNA \underset{Degradation}{\rightarrow} \oslash $$ $$ sfGFP \underset{Degradation}{\rightarrow} \oslash $$

The two equations above state the same time, the concentration of protein and mRNA would undergo degradation which means the concentration would drop. However, since there is always protein and mRNA being created, over time, the creation and degradation keep the concentration constant. ²

We can apply Law of Mass Action combine both equations for the concentration of protein and mRNA over time. This model can be described as:

^{Figure 3} $$ mRNA = k_{1} -d _{1 } mRNA $$ $$ Protein = k_{2} \cdot mRNA - d_{2} \cdot Protein $$

Where...

mRNA is the concentration of mRNA
Protein is the concentration of Protein
k₁ is the constitutive transcription rate. This represents the number of mRNA molecules produced per gene, per unit of time.
d ₁ is the mRNA degradation rate
k₂ is the translation rate. This represents the number of protein molecules produced per mRNA molecule, per unit of time.
d ₂ is the protein degradation rate.

This is important because we can use this model to calculate the concentration of proteins we can expect over time. This is useful as we can use this information to calculate the total emitted light spectra during the time period which is what we are looking for in our system. However, the constants and variables are individual for each protein and which means parameters for each protein would need to be found. These constants were found using literature ³ (for GFP) and lab results (the rest.)

¹ GB Stan, 20137. Modeling in Biology. London, the United Kingdom: Imperial College London. p, pp.59-65.

² See Non-Inhibited conditions from Figure 5 Gene Transcription Regulation by Repressors (CRISPRi) - Concentration over Time

³ See Relationship between Max Fluorescence and Protein Concentration for more details

Gene Transcription Regulation by Repressors (CRISPRi) - Concentration over Time

The next step in developing our simulation was to calculate our protein concentration at any given time when using CRISPRi. Discussion with wet-lab revealed our method would be using CRISPRi as a repressor, which works by inhibiting the expression of one or more genes by binding to the promoter region ¹. The expanded mRNA and Protein concentration models from the Constitutive Gene Expression Model ² were modified to include the element of repression from the CRISPRi inhibition.

$$ Gene \overset{Repressor}{\rightarrow} mRNA \rightarrow Protein $$ $$ mRNA \underset{Degradation}{\rightarrow} \oslash $$ $$ sfGFP \underset{Degradation}{\rightarrow} \oslash $$

This change can be applied to the Law of Mass Action ³ :

$$ m = k_{1} \cdot \frac{k^{n}}{k^{n} + R^{n}}- d_{1}m $$ $$ p = k_{2} m - d_{2}p $$

Where...

m is mRNA concentration, p is Protein concentration, R is Repressor, k1 is Max Transcription Rate, k is the Repression Coefficient, n is number of repressors that need to cooperatively bind the promoter to trigger the inhibition of gene expression (Hill Coefficient), R is Repressor, d1 is mRNA degradation rate, d2 is Protein degradation rate

The value for these constants and variables were taken from literature and calculating them ⁴ but later, adjusted to the lab results.

^{Figure 6}

Figure 6 shows the structure which underwent CRISPRi inhibition are expected to produce lower concentration of the protein whose expression were are inhibiting. This is important as it means the team can calculate concentration of proteins which are inhibited and compare them to the control conditions as well as giving the correct concentration for the simulation.

Furthermore, by having a model which can calculate the protein concentration at any given time, we can deduce how much fluorescence is being emitted at that time period by the bacteria

⁴ See Relationship between Max Fluorescence and Protein Concentration

Relationship between Max Fluorescence and Protein Concentration

At any given time, it is expected that the proteins would be expressed so the bacteria would fluoresce. This can be confirmed by looking at the bacteria after being constructed and observing that they are giving off light.

This means that an equation must be developed to find out what the intensity of fluorescence would be at that certain time. This consisted of of calculating the protein concentration at the time and using real life lab data of the fluorescence at that time period, the team could map that intensity to the protein concentration at that time.

When the fluorescence data received from the wet lab were graphed, a model was constructed, refined and optimised to demonstrate the trends shown from the real data gained from the labs. Originally, the data from the lab was the Fluorescence against Time but by using the Gene Transcription Regulation by Repressors model developed earlier ¹, the team was able to estimate the protein concentration at a certain time periods.

^{Figure 7}

Due to time constraints, rather than implementing the relationship directly from lab data, the data was fitted using a Polynomial Fit of Order 3 using Excel and an equation was calculated from these. These equations were directly plugged into the simulation. However, this is inaccurate as the R squared value was ... , suggesting that it doesn't fully capture the data trend.

These relationships were implemented into the simulation to give the expected spectra produced by each protein. This highlights another use: by adding or subtracting values from our fit, we can create a threshold for our Keys. This was essential when developing the Raw Data Simulator. ²

¹ See Gene Transcription Regulation by Repressors (CRISPRi) - Concentration over Time ² See Software

MODELLING

Overview

Constitutive Gene Expression For Protein and mRNA Expression over Time

Gene Transcription Regulation by Repressors (CRISPRi) - Concentration over Time

This change can be applied to the Law of Mass Action ³ :

Relationship between Max Fluorescence and Protein Concentration

Are Our Constructions Random?

Team:UNOTT/Modelling

MODELLING

Overview

Constitutive Gene Expression For Protein and mRNA Expression over Time

Gene Transcription Regulation by Repressors (CRISPRi) - Concentration over Time

This change can be applied to the Law of Mass Action 3 :

Relationship between Max Fluorescence and Protein Concentration

Are Our Constructions Random?

This change can be applied to the Law of Mass Action ³ :