Difference between revisions of "Team:Rice/Modeling"

Violacein is a fluorescent reporter with anticancer activity (Carvalho et. al, 2006) that has been used in several other iGEM projects (Cambridge 2009, Slovenia 2010, Johns Hopkins 2011, UCSF 2012). Although it would be a good pigment candidate for our project, it has a complex synthetic pathway requiring five specialized enzymes and oxygen (Fig. 2). It also presents multiple off-path reactions that can reduce the efficiency of the pathway. Before building constructs to use for violacein production, we needed to find a way to determine which promoters to use for the five genes involved in the pathway. Although there are studies focused on the optimization of the production of violacein (Lee et. al, 2013), none of the studies give a biochemical model of the rates of the reactions that take place in the bacteria.

Create a biochemical model of the violacein production based on the synthetic pathway and violacein production data from bacteria with different promoters for each of the five genes involved in the pathway.

1. The rate of dilution of the enzymes and the intermediaries is much greater than its degradation (for example by ubiquitination for the proteins or by conversion to products not included on the pathway)
2. There is no saturation of the enzymes and all the reactions will follow the law of mass action
3. Independence of external factors such as oxygen and NADH in the reactions
4. None of the reactions are reversible

We use the mass action kinetics because this type of equation only requires one parameter for reaction and is less susceptible to overdosing

Because a major goal of the model is to predict the effects of the selection of promoters on the final production of violacein, we decided to find a way to characterize promoters first. To simplify the computation, we used the promoter strength as a single standard to characterize the promoters. Moreover, we assumed the degradation rate of proteins only depends on the growth rate of E.coli. Then, every enzyme has the same degradation rate. The bacteriophage T7 promoter has been widely used for protein expression and purification (Jones et al., 2013), so we used data of five mutant T7 promoters to create a proof-of-concept model. If this model was functional, we could implement the same modeling technique to the promoters we were working with.The five mutant T7 promoters have distinct promoter strength over time after induction. The experimental data from the literature are shown in the figure below (Jones et al., 2013).



The first step of our model is to describe the rate of change of enzymes based on promoter strength. Here we assumed that the enzyme production rate is directly proportional to strength of the promoter. Therefore, we were able to use a mass-action kinetics equation of promoters to describe the enzyme concentration. The equation is shown below:



In this equation, Ai is the concentration of enzyme i, ki­ is the production rate of each enzyme i, kd is the degradation rate of all enzymes, and t is time. By solving this equation, we derived the equation of enzyme concentration against time.



Since we assumed that the promoter strength is proportional to the promoter concentration, we would use the equation to fit our data using least squares method (Fig. 1).



Figure 1. Linear regressions fitted to normalized fluorescence vs time. The circles represent data from Jones et al., 2013. The solid lines are our regression lines. The colors indicate with which promoters the circles and lines correspond.
In general, the regression lines are able to capture the change of strength of each enzyme over time. In this way, the parameters are determined. The table below lists the parameter values.



In the table, ki­ (i = 1,2,3,4,5) are the production rate coefficients of promoter I (i = 1,2,3,4,5), and kd is the degradation rate coefficient of all promoters.
After we finished the regression model of each promoter, we created a second model to describe the violacein biosynthetic pathway. The pathway (Fig. 2) involves five enzyme-catalyzed reactions and one non-enzymatic reaction (Lee et al, 2013).



Figure 2. Violacein synthetic pathway. The purple arrows highlight the five enzymatic and one non-enzymatic steps of violacein production from two molecules of tryptophan. The five enzymes are indicated by bolding (VioA, VioB, etc.).
The model was developed as three major parts. A pseudocode of this model is provided here.
Define ODE System

1. Calculate the production and degradation rate of each molecule in the pathway from the concentration of reagents and parameters.
2. Obtain the rate of change of each molecule based on the production and degradation rates.

Solve the System of Nonlinear Equations at Steady State

1. Solve the system of nonlinear equations at steady state starting at an initial guess X0.
2. Use the result as a new initial guess; repeat the numerical method to solve the system of equations again.
3. Calculate the relative error of each chemical in the new result.
4. If the maximum error is smaller than 0.0001%, output violacein concentration at steady state as the final result.

Optimize Parameters to Fit Experimental Data

1. Set the initial guess of the parameters.
2. Load the data from literature, which include the choice of promoter for each gene and the corresponding violacein yield determined experimentally.
3. For each promoter selection scenario, pass each promoter numbers and the temporary parameters to the steady-state model.
4. Obtain the violacein yield predicted by the steady-state model for each promoter selection scenario.
5. Compute the residual sum of squares (RSS) of between the predicted violacein yields and the violacein yields given by experiment.
6. Determine the optimal parameters by minimizing the RSS (least square method).