Team:Freiburg/Project/modeling

Introduction

Effectors

Effector Control

Targeting

uniBAss

Manual

Modeling

Introduction

The Networks

1. dCAS-VP16

Our network includes four different components dCas-VP16, a RNA complex (tracr/cr RNA), a RNA-dCas-VP16 complex and the Secreted alkaline phosphatase (SEAP). DCas-VP16 binds the RNA-complex and the whole complex binds the DNA, which leads to the production of SEAP.

Setting up the ODE

According to the graphical reaction network the ODE can be set up.

Cas9 is constitutively expressed by the CBh promoter and degraded proportional to the current concentration. It is used to build the DNA recognition complex and produced during complex decay.

The RNA-complex is build linearly. The production constant can be seen as production constant of the lower expressed RNA, because this expression limits the complex building. It is assumed that the RNA is degraded after DNA recognition complex decay and therefore the complex decay does not lead to more RNA.

The DNA recognition complex is built, when Cas9 and RNA meets and degraded proportional to the current DNA recognition complex concentration.

There is a leaky SEAP production and one that depends on the current concentration of the Cas9/RNA Complex. This dependency is assumed to follow the Monod-kinetic. Because of the long half time (T₂ > 500 h) of SEAP we can neglect the SEAP decay [2, 3].

The parameters are:
k1: linear production rate of Cas9 k2: Cas9 degradation rate k3: tracr/crRNA production rate k4: tracr/crRNA degradation rate k5: gene recognition complex building rate	k6: cr/trRNA /Cas9 degradation rate k7: SEAPs leaky production rate k8: Complex dependent SEAP production rate k9:

Finding the parameters

By setting up the ODE a n-dimensional hypothesis space (n is the number of parameters) is generated and finding the right parameter combination means finding a point in the space which fits the data best.
To find these parameters we used the maximum likelihood approach. The maximum likelihood hypothesis is the hypothesis which has the highest probability to generate the measured data. It is shown [3], that using the maximum likelihood approach and assuming gaussion noise in the data (an assumption that holds in our case) leads to a least-square error minimization problem.

A minimization problem is an optimization problem. You search for parameters (p0) for which holds, that the value of the function (f) at the point of the parameters is smaller than all other values. (f(p0)<=f(p)). In three dimensions the function can be thought as a landscape and minimization is finding the deepest valley. Depending on the method you use different problems arise. The most common problem is finding only a local minimum and not the global one (Fig. 3).

To avoid this and to be sure to have found a global minimum we started our minimization procedure using different start values for our parameters. To sample these parameters we used the latin hypercube sampling on a logarithmic scale (Fig. 4).

The number of different initial parameter settings is set to N and thus the parameter space is divided in N*N subspaces. For the initial parameter the values are chosen so that there is only one parameter in each row and column.

The resulting errors we plotted in an increasing order to be sure to have found a global minimum.

Data generation

Cas9 is quantified by using Western blot and we used SEAP as target protein that can be quantified by a SEAP assay. For more detailed information refer our modeling notebook.

Fitting Procedure and Results

1. dCas-VP16

Assuming the given ODE and using the fminsearch-function implemented in matlab with various initial parameter vectors the fitting process results in one optimal parameter composition. Our measurement time starts with the transfection, however the first data we got at time point 6, therefore all 14 parameters including the initial concentrations were set variable. This results in unlikely high initial concentration. Unlikely, because although the transfection end point is not clearly to define after transfection there should only be a small amount of the components. A second process in which the initial concentrations were assumed as zero followed the first one. The change between fixed and variable parameters was easily to perform, because of an additional vector (qfit). This vector contains boolean values depending on wether the parameter is fixed or flexible during the fitting process. Moreover another parameter is required for adjusting the absolute values of the duplicates. The results of this second process are shown. (Table 1).

The model reflects the general construction of the network (Fig. 5). As assumed Cas9 converges asymptotically to a stable state and there is an exponential increase in the SEAP concentration, both processes are reflected in the data. The model also shows a potential behaviour of the not measured components, the free tracr/crRNA-complex and the gene recognition complex. The tracr/crRNA-complex reachs it steady state quickly. There is however no possibility to distinguish between the two different RNAs. There might be some differences in their expressions, especially because of the different promoters (crRNA expressed under U6-promoter; tracrRNA expressed under h1-promoter), however the model won't show them.
The free gene recognition complex has, because of the small degradation rate, not reached its steady state yet.

Because of the fact that the fminsearch algorithm is not proved to converge to a minimum [1], we started at different points in the parameter space and therefore the probability of having found the global minimum is high Fig. 6 .

Table 1: Overview of the parameter results of the dCas9-VP16 fitting process

parameter	value
linear production rate of Cas9	8.78E+02 [M/h]
Cas9 degradation rate	0.54 [1/h]
Tracr/crRNA production rate	0.73
Tracr/crRNA degradation rate	1.21E+08
Gene recognition complex building rate	1.61 [1/M]
cr/trRNA /Cas9 degradation rate	8.59E-26
SEAPs leaky production rate	4.05E-05 [U/l*h]
k8	1.86E+08 [U/l*h]
k9	1.06
	2.06E-06
Cas9 degradation rate	0.00 [M]
Tracr/crRNA -Komplex	0.00
gene recognition complex	0.00
SEAP	2.73E-09 [U/l]
least square error	159.3423

The Code Files

References