Team:Carnegie Mellon/Mod-Matlab
From 2012.igem.org
Matlab Documentation
The inputs to the model are the measurement tables of concentration of dye vs. time. Optional inputs to the model include an in vitro measurement of saturation of the dye, and measurements of the fluorescence of the dye with mRNA and protein synthesis turned off. The first optional measurement can be used to compare the in vitro fluorescence saturation levels with the in vivo fluorescence saturation levels in order to give a scaling factor for the all the measurements in the input. Estimations can be used in place of these to simplify the number of inputs. The second optional measurement can be used to determine the degradation rates of mRNA and protein.
Fluoro2.m
This function is the function that is called to run the entire program. In addition, it takes the mRNA titration tables (modeldata) and converts it into fluorescent mRNA concentrations. It then passes the degradation data to the degradation functions, Degradation.m and DegradationP.m to return alpha2 and beta2, the degradation coefficients.
1 function [ PoPSans ] = Fluoro2( matrix, modeldata ) 2 %Fluoro 2: This function takes in a matrix of titrations and determines 3 %both the possible percentage for bound mRNA as well as the actual 4 %fluorescent mRNA concentration from fluorescent input values. 5 6 DNA = 1 * 10^(-9); 7 8 controlc = matrix(end,:); %concentration and fluorescence of the control 9 controlconc = controlc(1); %concentration of dye for the control 10 controldat = controlc(1:end); %fluorescence of the control 11 maxc = max(controldat); 12 concs = matrix(:,1); %concentrations of the dye in the wells 13 concs1 = concs(1:(end - 2)); 14 controlmax = maxc; 15 16 Rf = []; 17 18 for i = 2:size(matrix, 2); 19 fluordat = matrix(:,i); %fluoroscence data at some time point 20 fluordat1 = fluordat(1:(end - 2)); 21 s = fitoptions('Method', 'NonlinearLeastSquares', 'Startpoint',... 22 [fluordat1(end), 1/(controlconc)]); 23 g = fittype('a * (1 - exp(b * (-x)))', 'coefficients', {'a', 'b'},... 24 'options', s); 25 h = fit(concs1, fluordat1, g); 26 coeffvals = coeffvalues(h); 27 %figure(i); 28 %plot(h, concs1, fluordat1) 29 factorScale = coeffvals(1) / controlmax; 30 %scaling factor from in vitro to in vivo 31 for j = 1:size(concs1); 32 if abs(h(concs1(j)) - h(1)) <= (.2 * h(1)); 33 Rf(i - 1) = concs1(j) * factorScale; 34 break 35 end 36 end 37 38 end 39 time = matrix(1,:); 40 fluortime = time(2:size(matrix, 2)); 41 42 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 43 %Section for Transcriptional Efficiency % 44 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 45 46 Rt = FluoroLR(Rf, fluortime); 47 48 49 alpha2 = Degradation(modeldata); 50 51 ETF = mRNAexpress(fluortime, DNA, Rt, alpha2); 52 %ETF = mean(ETF); 53 54 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 55 %Section for Translational Efficiency % 56 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 57 58 Pt = ProteinFunctions(matrix); 59 60 beta2 = DegradationP(modeldata); 61 62 Tl = proteinexpress(DNA, Pt, ETF, alpha2, beta2); 63 64 %Tl = proteinexpress(DNA, 1.6, ETF, alpha2, beta2) 65 66 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 67 %Section for Polymerase per Second % 68 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 69 70 PoPSans = PoPS(alpha2, beta2, Pt, Tl); 71 72 end
Degradation.m
This function takes in mRNA fluorescence data with mRNA synthesis turned off. This makes determining degradation rates easier, as all the change (as we will define degradation) in the mRNA concentration will be due to degradation. The function takes the data and fits a curve to the data, in the process calculating the degradation rate.
1 function [ alpha2 ] = Degradation ( modeldata ) 2 dRi = modeldata(:,3); 3 time = modeldata(:,1); 4 %dRi(length(dRi),:) = []; 5 C = max(dRi); 6 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 7 %Differential Equation: dRi./dt = alpha * Rt% 8 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 9 10 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 11 %dRi is the fluorescence measurements of RNA during degradation only% 12 %alpha is the desired result, so we solve the diff eq % 13 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 14 15 for i = 1:(length(dRi)); 16 if dRi(i) ~= 0 && time(i) ~= 0 17 alpha(i) = log(dRi(i) ./ C) ./ time(i); 18 else 19 alpha(i) = log(dRi(2) ./ C) ./ time(2); 20 end 21 end 22 alpha; 23 alpha2 = -mean(alpha); 24 %Rdegra = Rt * alpha2; 25 26 end
DegradationP.m
This function has a similar role to Degradation.m. This function takes in protein fluorescence data with protein synthesis turned off. This function, similarly to Degradation.m, takes the data and fits a curve to the data.
1 function [ beta2 ] = DegradationP ( modeldata ) 2 %DegradationP 3 4 dPi = modeldata(:,6); 5 timeP = modeldata(:,1); 6 %dRi ./ dt = alpha * Ri, solution: R = Ce^(alpha t) ./ alpha 7 C1 = max(dPi); 8 %ln(Ri ./ C) ./ t = alpha 9 for i = 1:length(dPi) 10 if dPi(i) ~= 0 && timeP(i) ~= 0 11 beta(i) = log(dPi(i) ./ C1) ./ timeP(i); 12 else 13 beta(i) = log(dPi(2) ./ C1) ./ timeP(2); 14 end 15 end 16 17 beta2 = -mean(beta); 18 %Pdegra = beta2 * Pt; 19 end
The degradation functions return alpha2 and beta2 to Fluoro2.m. Fluoro2.m then calls ProteinFunctions.m to convert the protein titration data to fluorescent protein concentrations.
ProteinFunctions.m
This function does the same thing as Fluoro2.m with the mRNA titration data. It returns fluorescent protein concentrations over time.
Fluoro2.m takes the fluorescent protein and fluorescent mRNA concentrations and passes them to FluoroLR.m and FluoroLP.m to convert to total protein and total mRNA concentrations.
1 function [ Pt ] = ProteinFunctions( matrix) 2 %UNTITLED Summary of this function goes here 3 % Detailed explanation goes here 4 5 controlc = matrix(end,:); %concentration and fluorescence of the control 6 controlconc = controlc(1); %concentration of dye for the control 7 controldat = controlc(1:end); %fluorescence of the control 8 maxc = max(controldat); 9 concs = matrix(:,1); %concentrations of the dye in the wells 10 concs1 = concs(1:(end - 2)); 11 controlmax = maxc; 12 13 Pf = []; 14 15 for i = 2:size(matrix, 2); 16 fluordat = matrix(:,i); %fluoroscence data at some time point 17 fluordat1 = fluordat(1:(end - 2)); 18 19 s = fitoptions('Method', 'NonlinearLeastSquares', 'Startpoint',... 20 [fluordat1(end), 1/(controlconc)]); 21 22 g = fittype('a * (1 - exp(b * (-x)))', 'coefficients', {'a', 'b'},... 23 'options', s); 24 25 h = fit(concs1, fluordat1, g); 26 coeffvals = coeffvalues(h); 27 %figure(i); 28 %plot(h, concs1, fluordat1) 29 factorScale = coeffvals(1) / controlmax;... 30 %scaling factor from in vitro to in vivo 31 for j = 1:size(concs1); 32 if abs(h(concs1(j)) - h(1)) <= (.2 * h(1)); 33 Pf(i - 1) = concs1(j) * factorScale; 34 break 35 end 36 end 37 38 end 39 time = matrix(1,:); 40 fluortime = time(2:size(matrix, 2)); 41 42 Pt = FluoroLP (Pf, fluortime); 43 44 end 45
FluoroLR.m
This function takes in fluorescent mRNA concentrations and converts it to mRNA concentrations using first-order chemical reactions. One dye molecule will bond to one mRNA molecule, creating an mRNA-dye complex. This leads to a rather simple conversion using the known dye concentration.
\begin{equation}K_{D_R} = \frac{[R_f]}{([R]_0 - [R_f])([D_R]_0 - [R_f])}\end{equation}Where KDR is the dissociation constant, $[Rf]$ is the fluorescent mRNA concentration, $[R]$ is the mRNA concentration, and $[D_R]$ is the dye concentration.
1 function [Rt] = FluoroLR (Rf, concentrations) 2 %DFHBI = 1*10^(-9); %concentration of dye 3 %Rfmax = 1*10^(-10);%max concentration of the mRNA 4 KD = 464*10^(-9); 5 6 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 7 %Relates the mRNA fluorescence levels with the total mRNA levels% 8 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 9 10 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 11 %fluorescence is measured, Rf is the fluorescent mRNA concentration% 12 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 13 14 %a = max(fluorescence) ./ Rfmax; 15 %Rf = zeros(length(fluorescence) - 1, 1); 16 %for i = 1:(length(fluorescence) - 1) 17 % Rf(i) = fluorescence(i) ./ a; 18 %end 19 20 %%%%%%%%%%%%%%%%%%%%%% 21 %Rt is the total mRNA% 22 %%%%%%%%%%%%%%%%%%%%%% 23 24 Rt = []; 25 for j = 1:(length(Rf)) 26 Rt(j) = Rf(j) * (1 + KD ./ (concentrations(j) - Rf(j))); 27 end 28 29 end
FluoroLP.m
This function has a similar role to FluoroLR.m, except using fluorescent protein concentrations and converting it to protein concentration.
1 function [Pt] = FluoroLP (Pf, concentrationsP) 2 MAG = 1*10^(-9); %concentration of dye 3 Pfmax = 1*10^(-10); %max concentration of the protein 4 KD2 = 464*10^(-9); 5 6 7 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 8 %Relates the protein fluorescence levels with the total protein levels% 9 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 10 11 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 12 %fluorescence is measured, Pf is the fluorescent protein concentration% 13 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 14 15 %b = max(fluorescenceP) ./ Pfmax; 16 %Pf = zeros(length(fluorescenceP), 1); 17 %for i = 1:(length(fluorescenceP)) 18 % Pf(i) = fluorescenceP(i) ./ b; 19 %end 20 21 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 22 %Pt is the total protein concentration% 23 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 24 25 Pt = []; 26 for j = 1:(length(concentrationsP)) 27 Pt(j) = Pf(j) * (1 + KD2 ./ (concentrationsP(j) - Pf(j))); 28 end 29 end
mRNAexpress.m
Fluoro2.m then takes the total mRNA concentrations passed by FluoroLR.m and passes it to mRNAexpress.m, which calculates the transcriptional efficiency. This is done via the differential equation
\begin{equation}\frac{d[R]}{dt} = Ts \cdot [D] - \alpha \cdot [R]\end{equation}to which the solution is
\begin{equation} Ts = \frac{[R] \cdot \alpha}{[D] \cdot (1 - e^{-\alpha \cdot t})}\end{equation}
where $[R]$ is mRNA concentration, $Ts$ is transcriptional efficiency, $[D]$ is DNA concentration, and $\alpha$ is the mRNA degradation coefficient.
proteinexpress.m
Fluoro2.m takes the transcriptional efficiency from mRNAexpress.m, total protein concentrations from FluoroLP.m, and alpha2 and beta2 from Degradation.m and DegradationP.m, and passes them to proteinexpress.m. proteinexpress.m computes the translational efficiency using the differential equation
\begin{equation}\frac{d[P]}{dt} = [R] \cdot Tl - \beta \cdot [P]\end{equation}
to which the solution is
\begin{equation}Tl = \frac{[P]}{\frac{Ts \cdot [D]}{(\alpha \cdot \beta)} \cdot (1 - e^{-\beta \cdot t}) - \frac{Ts \cdot [D]}{\alpha \cdot (-\alpha + \beta)} \cdot (e^{-\alpha \cdot t} - e^{-\beta \cdot t})} \label{eq:Tl}\end{equation}
where Tl is translational efficiency, beta is the protein degradation coefficient, and [P] is the protein concentration.
1 function [ Tl ] = proteinexpress (DNA, Pt, ETF, alpha2, beta2) 2 %proteinexpress 3 4 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 5 %differential equation comes from: dP./dt = Trans * RNA - beta * P% 6 %want to solve for Trans, the translational efficiency % 7 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 8 9 for k = 1:length(Pt) 10 Trans(k) = Pt(k) ./ (ETF * DNA ./ (alpha2 * beta2) *... 11 (1 - exp(-beta2 * (k - 1))) - ETF * DNA ./ (alpha2 *... 12 (-alpha2 + beta2)) * (exp(-alpha2 * (k - 1)) - exp(-beta2 * (k - 1)))); 13 end 14 15 %for k = 1:length(ETF) 16 % Pt(k) = Trans * (ETF(k) * DNA ./ (alpha2 * beta2) * (1 -... 17 % exp(-beta2 * (k - 1))) - ETF(k) * DNA ./ (alpha2 *... 18 % (-alpha2 + beta2)) * (exp(-alpha2 * (k - 1)) - ... 19 % exp(-beta2 * (k - 1)))); 20 % 21 %end 22 23 Tl = mean(Trans); 24 end
PoPS.m
Fluoro2.m passes to the final function alpha2 and beta2 from Degradation.m and DegradationP.m, total protein concentration from FluoroLP.m, and translational efficiency from proteinexpress.m. PoPS.m calculates the approximate polymerase per second using the equation
\begin{equation}PoPS = \frac{\alpha \cdot \beta \cdot [P]}{n \cdot Tl} \label{eq:PoPS}\end{equation}
where n is the approximate number of the promoters of interest in a cell (i.e. plasmid copy).
1 function [ PoPS ] = PoPS( alpha2, beta2, Pt, Tl ) 2 %This function calculates Polymerase per second given a few parameters 3 %This equation is valid at steady state. 4 5 n = 1; 6 PoPS = alpha2 * beta2 * Pt ./ (n * Tl); 7 8 end
The model outputs polymerase per second, although transcriptional efficiency and translational efficiency are also important factors in the model. Derivations of these equations can be found on the derivations page.