Difference between revisions of "Team:Valencia UPV/Experiments"

Revision as of 14:59, 6 October 2018

Stack Multipurpose HTML Template

Index

Human Practices
Software
Hardware
Wet Lab
Modeling
Graphic Design
Wiki
Organization / Communication

Modeling

Do you think it is possible to mathematically describe a cell? Would you like to know the possibilities that Printeria modeling offers you?

One of the fundamental bases of the Printeria project has undoubtedly been mathematical modeling. Thanks to the development and application of new mathematical models, it is possible to quantify the expression of proteins in cells, and therefore characterize through different experiments the parts designed by Printeria. From the Printeria modeling team, we intend to reach different goals:

Design simple mathematical models based on differential equations that describe the biochemical processes of a cell. With them, we can simulate the different genetic circuits that Printeria allows us to build.
Develop a Simulation Tool that allows the user to visualize a prediction of the results of their experiment before running it in Printeria.
Optimize model parameters to match simulation results to experimental data obtained from Printeria constructions.
Characterize the parts of our Part Collection from the optimization results and provide the user with all the information about the Printeria kit.

Although in the development of the project we have dealt with all these aspects, all of them have a single purpose: to understand and describe in a mathematical way the biological processes that take place inside the cell.

Models & Experiments

Characterization procedures

Simulation Tool

Models & Experiments

The development of new and simple mathematical models has been one of the essential bases of the Printeria project. The models offer us a multitude of applications: to describe the basic biochemical reactions of the cell, to characterize the parts of our Printeria kit by means of different experiments or to elaborate new tools that facilitate to our user the learning of Synthetic Biology.

In Printeria Modeling team, we have grouped the models designed in two fields: constitutive expression models and inducible expression models. We have also established an experimental protocol and a multi-objective optimization protocol in order to contrast theoretical and experimental results. With them, we can describe practically all the cells modified with Printeria.

Want to find out more about how we did it?

Modeling process

At the beginning of the design of any mathematical model, we have started from a cell scheme in which all the cellular biochemical reactions we wish to model are represented. From the reactions, a set of equations is deduced. Each of these equations describes the temporal variation of the biochemical species in the cell (DNA, RNA, proteins...), and depends on a set of parameters with a physical meaning.

Because they are differential and generally non-linear equations, mathematical models are simulated using software tools such as MATLAB. The results obtained from the simulation reflect the evolution of the concentration of biochemical species over time. These theoretical results can be contrasted with the experimental data obtained in the laboratory and and thereby adjust the parameters so that the model adapts to the data in an optimal way. This process, called multiobjective parameter optimization, is the process that allows us to validate our models and see if they respond to what happens inside the cells.

In the development of the set of mathematical models that describe the genetic circuits printed by Printeria, we have prioritized simplicity over complexity of the model. A more complex model may be more precise, but it requires a large number of parameters whose value is often unknown. Therefore, our models are composed of few equations and parameters, and they are all deterministic models, that is, their parameters take a concrete, non-probabilistic value.

Experimental protocol

The Printeria Modeling and Lab team and have jointly designed an experimentation protocol for the laboratory experiments. Thanks to it, and from the colonies of the different UT the experimental data can be obtained, processed and ready to be optimized.

Materials:

Printeria transcriptional units (see our Printeria Part Collection)
Measuring equipment: Biotek Cytation3
96 well plate
MATLAB 2018a software

Protocol:

The TU colonies assembled by our Lab team are diluted in the medium (generally we have used LB or M9 medium with Kanamycin) so that they all have the same OD.

The experiment is designed on the measuring equipment. In our case, we use the Biotek Cytation3 equipment. We establish the equipment parameters.

Parameters	Description
Time	06:00:00 (HH:MM:SS) usually. Measurement interval: 05:00 (MM:SS)
Number of samples	We normally set 8 samples of reporter protein for each TU colony
Number of medium samples	We normally set 8 samples of medium. We normally use LB or M9 medium with Kanamycin
Temperature	37 ºC
Shake	Double Orbital. Continuously. We shake the plate before each measure
Absorbance. Optical Density (OD) measure	Wavelenght at 600 nm emission
Excitacion wavelength	We normally set 485 nm
Emission wavelength	We normally set 528 nm
Gain (G)	Normally the gain value is G = 60, although for proteins with lower fluorescence, it is recommended that G takes higher values.

Samples are introduced into the 96 well plate and the experiment begins.
After the experiment, fluorescence and absorbance data obtained are exported to an Excel file.
We run the MATLAB convert_data.m script. This script uses several additional MATLAB files with which:
1. We extract the fluorescence and absorbance matrices from the Excel file (see readExperiment3.m MATLAB function).
2. We apply the corrections to the data (see plotExperiment_yb.m MATLAB function).
3. We save the data of FOD, OD, molecules and time in a .mat format file.

Multi-objective parameter optimization

To what extent is our model valid? What values should the parameters take? To what extent do the theoretical results resemble the experimental ones? Is our model able to explain the behaviour of cells in reality?

One of the most important phases in the field of Printeria Modeling has been the validation of theoretical models with experimental data. In the Printeria team we wanted to confirm that the models designed are consistent with our experiments. However, we know that the parameters will not take a fixed value, but vary depending on the genetic construction with which we experiment. With this idea in mind, the need to apply an optimization process arises.

The process of multiobjective parameter optimization consists in finding those parameters of the model that best fit a particular objective, i.e., a particular set of experimental data. From a dataset, and applying a mathematical algorithm, we can find the optimal solutions, also known as Pareto-optimal or Pareto front solutions. These solutions are such that there are no more solutions that do not improve one objective without worsening the rest.

From the Printeria Modeling team we have established a protocol that allows us to obtain the optimal solutions of the model parameters for any experiment:

We define the objectives to be optimized. These will be the values of Fluorescence F/Unit of Absorbance OD or FOD, and the cellular growth (OD) of a Transcriptional Unit (TU). Therefore, if we have different n TU in an experiment, we will optimize 2n objectives.
We choose the parameters to optimize and establish the range of values they can take.
We select the identification and validation data. The identification data are those used in the optimization process. The validation data are used to compare the experiment with the simulated theoretical model and the optimized parameters.
We define the cost function, which specifies the error function to be minimized, simulates the model with different parameter values, and calculates the error between the simulation data and the experimental identification data.
The mathematical optimization algorithm is executed. In our case we use the Multi-objective Differential Evolution Algorithm with Spherical Pruning. This algorithm uses the cost function to test different parameter values and thus search for optimal parameter values.
The Pareto set and Pareto Front diagrams and the different solutions of our optimized parameters for the different objectives are obtained.
Decision Making is carried out by the designer. The optimal parameters are chosen for our 2n objectives.
The model is simulated with the optimized parameters and compared with the validation data.

Following this protocol, and particularizing it for each experiment, we have been able to determine the validity of our models and obtain the best parameters for our experiments.

Constituve expression models

Constitutive expression models are those in which the expression of the cell's messenger RNA (mRNA) is not regulated. The gene does not depend on any inducing biochemical species, so it will produce RNA molecules continuously, and therefore the translation of the protein will also be constant.

Model design

A simple scheme that summarizes the cellular transcription and translation processes in a constitutive expression model is shown in the following Figure. We only consider as biochemical species involved in reactions to mRNA and our protein of interest (PoI).

From the cell schema we can deduce the following biochemical reactions.

Parameter	Description	Units	Value
C_R	Constitutive transcription rate is calculated by multiplying the transcription rate per plasmid K_R by the mean number of plasmids in cell c_n . In our case, we are using the pMBI replication origin, so c_n ∼ 500, and C_R = K_R·c_n.	molecules.min^-1	---
p	Translation rate	min^-1	---
d_m	mRNA degradation rate	min^-1	---
d_p	PoI degradation rate	min^-1	---
μ	Dilution rate	min^-1	---
K_max	Maximum growth capacity	cells	---

With the reactions and approximations already established, we apply the Law of Mass Action kynetics (LMA) to deduce the differential equations. The LMA establishes that the variation of the species resulting from a reaction is proportional to the product of the reactants. In our case, to the set of differential equations deduced from the reactions we add a cell growth equation. Biochemical species and our model of equations obtained are shown below.

Variable	Biochemical species	Units
n₁	mRNA	Molecules
n₂	PoI	Molecules
n₃	Number of cells	Cells

By obtaining the equations, and seeking the simplicity of the model, a series of quasi-static approximations have been made. We have established that mRNA is generated with a constant C_R constitutive transcription rate. Another consideration has been to assume that the copy number of plasmids in cell is constant. We have also proposed that other species such as ribosomes, which also participate in transcription, are constant over time.

After developing the model and raising the reactions and equations, two files have been programmed in MATLAB: the mc_simple.m function, which describes the system of differential equations, and the script model_const_simple.m in which the parameters are defined and the model is solved by the ode45 MATLAB function. The result of the simulation of the model is shown in the following Figure.

In the simulation of mRNA and PoI we distinguish two phases in the temporal evolution of the species: the transient phase and the stationary phase. If we perform different simulations varying some parameters of the model, we can observe how the equilibrium value of the stationary phase varies. MATLAB script mc_sim_analysis.m repeatedly simulates the model by varying the C_R and p parameters, and graphs the results. The following Figure shows how an increase in both ratios implies a greater expression of the protein.

At this point anyone could ask us: Which are the advantages of a constitutive expression model like this? Why is a model like this really useful to us? Faced with these questions, the Printeria Modeling team has looked for the answers in the great characteristics that the model offers us.

A simple, easy-to-understand model that clearly explains the processes of cellular transcription and translation.
It is valid for any Printeria construction that has a constitutive expression promoter. The variation between different constructions occurs in the parameters values, not in the model equations!
It has few parameters, with a clear physical meaning, and easily optimizable.
As a compact model, simulations and optimizations are performed at high speed.

Experiments & optimization

The simple model of constitutive expression is a model that characterizes a large number of Printeria Transcriptional Units (TU), specifically, all those constructions whose promoter is constitutive. We know that all these TU follow the same model, but each one has different parameters, and therefore different experimental values of fluorescence and absorbance. Faced with this situation, the Printeria Lab and Modeling teams have designed some experiments changing RBS and promoters, in which we have applied the multi-objective parameter optimization, and thus check whether the designed constitutive expression model responds to experimental reality.

Experiment parameters	Description
Time	06:00:00 (HH:MM:SS). Measurement interval: 05:00 (MM:SS)
Number of samples	8 samples for each TU colony
Number of medium samples	8 samples of M9 medium
Temperature	37 ºC
Shake	Double Orbital. Continuously
Absorbance. Optical Density (OD) measure	Wavelenght at 600 nm emission
Excitacion wavelength	Wavelenght at 485 nm
Emission wavelength	Wavelenght at 528 nm
Gain (G)	60

Experiments changing RBS

We have designed two experiments following the same experimental protocol. In them we have assembled different Printeria TU with the same promoters, CDS (sfGFP reporter protein) and transcriptional terminator, but with different RBS. After obtaining the results, and following the optimization protocol, we have obtained the parameters of the model and have validated our model.

Printeria RBS:

Strong expression: BBa_K2656009.
Medium expression: BBa_K2656011.
Low expression: BBa_K2656010.
Very low expression: BBa_K2656008, BBa_K2656012.

Optimization specifications	Description
Parameters	Constitutive transcription rate C_R: fixed Translation rate p: to optimize PoI degradation rate d_p: to optimize mRNA degradation rate d_m: fixed Dilution rate μ: to optimize Maximum growth capacity K_max: experimental value
Objetives to optimize	For each TU, we set 2 objetives: FOD (Fluorescence) and OD (Absorbance). In each experiment we have measured 3 TU, so we must set 6 objectives par experiment
Parameter ranges	Translation rate p: [0.001 - 6] min^-1 PoI degradation rate d_p: [0.0058 - 0.0087] min^-1 Dilution rate μ: [0.0058 - 0.035] min^-1
MATLAB files	spMODEparam.m: it defines the parameters to be optimized, their value ranges, number of objectives, the cost function, the identification and validation experimental data and other spMODE algorithm parameters. When the script is executed, a spMODEDat structure variable is defined. This structure contains all the standardized optimization information that the spMODE algorithm will need for its execution. Experiment 1: spMODEparam_exp1sfGFP.m Experiment 2: spMODEparam_exp2sfGFP.m spMODE.m: it contains the Multi-objective Differential Evolution Algorithm with Spherical Pruning, which optimizes parameters for our experimental results (to execute spMODE.m file we need also SphPruning.m file). CostFunction_RBS_RMS_2n.m: it simulates the model with different vector parameters, and calculates the Root Mean Square Error with the identification dataset. Then, it returns the error for each objective and parameter vector. levelDiagram.m: plots Pareto front and Pareto set that gives us back the algorithm. execute_RBS_2n.m: this script launches the optimization by executing spMODEparam.m, spMODEm algorithm and levelDiagram.m. Then, allows the user to enter the best parameters, and plots validation experimental data and simulation results in the same graph.

The results of the graphs have been obtained after an optimization and decision-making process, in which the optimal parameters have been selected and the constituent model for each objective has been simulated and compared with the experimental data. The optimized parameters are summarized in the following table.

Optimized parameters

Values

Translation rate p

Experiment 1:

For strong RBS BBa_K2656009: p = 0.2059 min^-1
For medium RBS BBa_K2656011: p = 0.082 min^-1
For low RBS BBa_K2656010: p = 0.0072 min^-1

Experiment 2:

For strong RBS BBa_K2656009: p = - min^-1
For very low RBS BBa_K2656008: p = - min^-1
For very low RBS BBa_K2656012: p = - min^-1

PoI degradation rate d_p

Experiment 1: d_p = 0.0058 min^-1

Experiment 2: d_p = - min^-1

Dilution rate μ

Experiment 1:

For strong RBS BBa_K2656009: μ = 0.01492 min^-1
For medium RBS BBa_K2656011: μ = 0.01631 min^-1
For low RBS BBa_K2656010: μ = 0.0188 min^-1

Experiment 2:

For strong RBS BBa_K2656009: μ = - min^-1
For very low RBS BBa_K2656008: μ = - min^-1
For very low RBS BBa_K2656012: μ = - min^-1

In addition to simulating the optimized models, the Printeria Modeling team has also calculated the relative force between the different RBS, taking BBa_K2656009 strong RBS as a reference. The relative force has been defined as the quotient between the values of the protein in equilibrium of the results of the simulation of one RBS and another reference RBS. Likewise, a ratio between p parameters of the different RBS parts and p parameter of the reference RBS has been calculated. The characterization of the RBS parts by their relative strength is shown below.

Printeria RBS part	Relative strength	p parameter ratio (p_RBS/p_ref)
BBa_K2656009 (Reference)	1	1
BBa_K2656011	0.371	0.398
BBa_K2656010	0.029	0.035
BBa_K2656008	-	-
BBa_K2656012	-	-

Experiments changing promoters

We have designed two experiments following the same experimental protocol. In them we have assembled different Printeria TU with the same RBS, CDS (GFP reporter protein) and transcriptional terminator, and with different promoters. After obtaining the results, and following the optimization protocol, we have obtained the parameters of the model and have validated our model.

Printeria promoters:

Strong promoters: BBa_K2656005
Medium promoters: BBa_K2656007
Low promoters: BBa_K2656004

Optimization specifications	Description
Parameters	Constitutive transcription rate C_R: to optimize Translation rate p: fixed PoI degradation rate d_p: to optimize mRNA degradation rate d_m: fixed Dilution rate μ: to optimize Maximum growth capacity K_max: experimental value
Objetives to optimize	For each TU, we set 2 objetives: FOD (Fluorescence) and OD (Absorbance). In each experiment we have measured 3 TU, so we must set 6 objectives par experiment
Parameter ranges	Constitutive transcription rate C_R: [ - ] min^-1 PoI degradation rate d_p: [0.0058 - 0.0087] min^-1 Dilution rate μ: [0.0058 - 0.035] min^-1
MATLAB files	spMODEparam.m: it defines the parameters to be optimized, their value ranges, number of objectives, the cost function, the identification and validation experimental data and other spMODE algorithm parameters. When the script is executed, a spMODEDat structure variable is defined. This structure contains all the standardized optimization information that the spMODE algorithm will need for its execution. Experiment 1: -- spMODE.m: it contains the Multi-objective Differential Evolution Algorithm with Spherical Pruning, which optimizes parameters for our experimental results (to execute spMODE.m file we need also SphPruning.m file). CostFunction_Prom_RMS_2n.m: it simulates the model with different vector parameters, and calculates the Root Mean Square Error with the identification dataset. Then, it returns the error for each objective and parameter vector. levelDiagram.m: plots Pareto front and Pareto set that gives us back the algorithm. execute.m: this script launches the optimization by executing spMODEparam.m, spMODEm algorithm and levelDiagram.m. Then, allows the user to enter the best parameters, and plots validation experimental data and simulation results in the same graph.

GRÁFICAS RESULTADOS 2 EXPERIMENTOS PROMOTORES

The results of the graphs have been obtained after an optimization and decision-making process, in which the optimal parameters have been selected and the constituent model for each objective has been simulated and compared with the experimental data. The optimized parameters are summarized in the following table.

Optimized parameters	Values
Constitutive transcription rate C_R	For strong promoter BBa_K2656005: C_R = - min^-1 For medium promoter BBa_K2656007: C_R = - min^-1 For low promoter BBa_K2656004: C_R = - min^-1
PoI degradation rate d_p	Experiment 1: d_p = - min^-1
Dilution rate μ	For strong promoter BBa_K2656005: μ = - min^-1 For medium promoter BBa_K2656007: μ = - min^-1 For low promoter BBa_K2656004: μ = - min^-1

In addition to simulating the optimized models, the Printeria Modeling team has also calculated the relative force between the different promoters, taking BBa_K2656005 promoter as a reference. Likewise, a ratio between C_R parameters of the different promoter parts and C_R parameter of the reference promoter has been calculated. The characterization of the promoter parts by their relative strength is shown below.

Printeria RBS part	Relative strength	p parameter ratio (p_RBS/p_ref)
BBa_K2656009 (Reference)	1	1
BBa_K2656011	-	-
BBa_K2656010	-	-

Inducible expression models

Inducible expression models are those in which the expression of the cell's messenger RNA (mRNA) depends on one or more inducing molecules, which can act as activators or repressors. The concentration of the inducing species determines the level of expression of mRNA and protein in the cell.

Our Printeria device gives us the possibility to create some inducible genetic constructs. Of all of them, we have managed to design and experiment with some models...

P_BAD/araC inducible model
LuxR/LuxI inducible model

P_BAD/araC inducible model

The inducible promoter P_BAD/araC is a promoter whose expression depends on two inducing molecules: the L-arabinose (activator) and the araC (repressor). The dimeric protein araC binds to the DNA chain forming a DNA loop and thus prevents the binding between DNA and RNA polymerase. However, when two arabinose molecules bind to the araC dimer, the DNA loop is broken and the binding between the RNA polymerase to the promoting region is possible.

We have proposed a possible mathematical model describing the P_BAD/araC mechanism. This theoretical model has been included in the Simulation Tool, so that with Printeria it is possible to predict the behavior of P_BAD/araC inducible promoter circuits.

Model design

After analyzing the biochemical functioning of the Pbad/araC promoter, we propose the following cellular scheme that would describe its behavior. It is worth mentioning that we have assumed a very high concentration of araC in the cell, as we have considered that the gene encoding for the araC protein is found in the cell genome (e.g. in DH5α strains).

From the cell schema we can deduce the following biochemical reactions.

Parameter	Description	Units	Value
D	Diffusion coefficient	min^-1	---
V_c	Cell volume(μL)/External volume(μL)	adimensional	---
k_u	Dissociation rate of (P_BAD.araC.arabi) complex	min^-1	---
k_d	Dissociation constant of (P_BAD.araC.arabi) complex	molec²	---
k_b	Association rate of (P_BAD.araC.arabi) complex	molec^-2min^-1	---
K	Transcription rate is calculated by multiplying the transcription rate per plasmid K_t by the mean number of plasmids in cell c_n . In our case, we are using the pMBI replication origin, so c_n ∼ 500, and C_R = K_R·c_n.	min^-1	---
α	(P_BAD.araC) basal activity constant	Admimensional	---
p	Translation rate	min^-1	---
d_m	mRNA degradation rate	min^-1	---
d_p	PoI degradation rate	min^-1	---
μ	Dilution rate	min^-1	---
K_max	Maximum growth capacity	cells	---

With the reactions and approximations already established, we apply the Law of Mass Action kynetics (LMA) to deduce the differential equations. The LMA establishes that the variation of the species resulting from a reaction is proportional to the product of the reactants. In our case, to the set of differential equations deduced from the reactions we add a cell growth equation. Biochemical species and our model of equations obtained are shown below.

Variable	Biochemical species	Units
n₁	mRNA	Molecules
n₂	PoI	Molecules
n₃	Number of cells	Cells

LuxR LuxI inducible model

Characterization procedures

This year, Valencia UPV iGEM team has designed an extensive Part Collection in purpose of allowing the user to design multiple genetic constructions and experiments. One of our main objectives has been to show the user clear and structured information about the pieces that make up the Printeria kit. For this reason, we have considered the characterization of the parts as a priority when developing the project. In this way, we have elaborated some procedures which have allowed us to systematically obtain and structure information from the parts.

Procedure for obtaining protein spectra
Comparison between sfGFP and GFP relative fluorescence intensity

Procedure for obtaining protein spectra

Obtaining excitation and emission spectra is a fundamental aspect in the process of characterization of a fluorescent protein. Each protein has a characteristic spectrum, which indicates the energy in which the molecule is excited or emits at a certain wavelength.

The characterization of the reporter proteins by excitation and emission spectra are of great importance in experimentation. When contrasting experimental information with the theoretical results of mathematical models, we experiment with numerous reporter proteins, such as fluorescence proteins or chromoproteins. However, the fluorescence data obtained must be corrected applying diferent operations in order to obtain representative fluorescence data:

The subtraction of the medium fluorescence:

The quotient of the fluorescence with the gain:

The quotient of the fluorescence with excitation and emission efficiency:

Until now, the corrections applied to the experiments performed with these reporter proteins were the subtraction of the fluorescence of the medium, and the division by the gain factor of the measuring equipment. However, with the protein spectra we can also normalize the fluorescence data to values that would have been obtained with maximum excitation and emission.

Owing to this reason, a protocol has been established in the lab by Lab and Modeling team to obtain the spectrum of any reporter protein.

Materials:

Measuring equipment: Biotek Cytation3
96 well plate
MATLAB 2018a software

Procedure:

We look for the theoretical spectra of the protein to be measured or similar molecules in order to determine the wavelength at which the protein is excited or emitted at maximum energy, i.e. where the theoretical spectral peaks occurs.

We define the protocol of our equipment to get the absorbance and fluorescence dataset. In our protocol, the most important parameters to be established are summarized in the following Table.

Parameters	Description
Number of samples	From 3 to 6 samples of reporter protein and 3 samples of medium (LB or M9 medium)
Temperature	37 ºC
Shake	Double Orbital. 01:00 (MM:SS)
Absorbance. Optical Density (OD)	Wavelenght at 600 nm emission
Excitacion and emission scans	The scans occur between two wavelength limit values. The established range will depend on the theoretical spectrum of the protein.
Excitacion and emission wavelengths	These values will depend on the range, and therefore on the spectrum. Values far from the theoretical peak lead to more attenuated fluorescence curves, and values very close to the peak can lead to overlap and error in reading the data. Therefore, a compromise must be reached between curve resolution and reading overlap.
Gain (G)	Normally the gain value is G = 60, although for proteins with lower fluorescence, it is recommended that G takes higher values.

The experiment is introduced, and the experimental absorbance data and the fluorescence curves of the samples are obtained with the medium fluorescence correction applied.
The dataset is exported to an Excel file.
The MATLAB script spectrum.m for fluorescent proteins is executed:
1. Dataset is extracted from the Excel file. We discard readings that have suffered overlap, or that take negative values.
2. The fluorescence curves of all samples are averaged, and the result is normalised (from 0 to 100%).
3. Graphs of the normalized absorption and emission spectra are plotted. The X-axis represents the wavelength (nm), and the Y-axis represents the normalized fluorescence intensity (%).

In practical terms, the protocol has been applied to all of our reporter proteins: GFP, sfGFP, YFP and mRFP. All results can be found in our parts collection as well as in the iGEM catalog.

In the particular case of a reporter chromoprotein, such as amilCP, we do not measure fluorescence, but absorbance. In this case, in rder to obtain the corrected absorbance curve, we must subtract from the cell absorbance data with the reporter protein the absorbance of a medium with cells without chromoprotein. Once the data have been corrected, we normalize them between values of 0 and 100 and with this we elaborate the graph. The protocol used can be found in the MATLAB script amilCP_spectrum.m.

Finally, we have also obtained the spectra of the fluorescein molecule. These spectra have been used to correct the fluorescence data used in the Interlab Study to obtain the Relative Fluorescence Units (RFU) to Molecules of Equivalent Fluorochrome (MEFL) conversion factor. In addition, the comparison between GFP and sfGFP proteins RFU has been possible thanks to fluorescein spectra.

Parameters	Value
Number of samples	6 samples. 3 samples of medium (LB or M9 medium)
Temperature	37 ºC
Shake	Double Orbital. 01:00 (MM:SS)
Absorbance. Optical Density (OD)	Wavelenght at 600 nm emission
Emission range	[495 - 580] nm
Fixed excitation wavelength	480 nm
Excitation range	[430 - 520] nm
Fixed emission wavelength	545 nm
Gain (G)	60

Comparison between sfGFP and GFP relative fluorescence intensity

One of the main problems we encounter when processing the results of experiments in Synthetic Biology are the units of measurement of fluorescence data. Unlike absorbance, where there is a simple conversion between Optical Density (OD) and cell number, there is no direct relationship between the Relative Units of Fluorescence (RFU) and the number of protein molecules in the cell. Moreover, RFUs vary among reporter fluoroproteins: for example, an RFU of the GFP protein does not have to be equivalent to an RFU of the mRFP protein, sfGFP, etc.

Thanks to initiatives such as the Interlab Study, we have been able to go a step further and obtain a MEFL/cell_GFP factor of equivalence between the RFU of the GFP protein and the Molecules of Equivalent Fluorochrome (MEFL). This relationship is an important breakthrough, as it can give us a more accurate estimation of the amount of GFP molecules in the cell.

However, another reporter protein widely used in the experiments is the sfGFP. This protein has a much faster folding than GFP, which translates into a higher fluorescence intensity per molecule. In order to obtain the MEFL/cell_sfGFP factor from sfGFP, the Printeria Modeling and Lab teams have designed a comparative experiment between both proteins. The experiment consists of designing two identical transcriptional units (TU), changing only the CDS sequence so that each TU will produce GFP and sfGFP, respectively. It should be added that this experiment is based on two fundamental assumptions:

The number of GFP molecules produced in the cells is equivalent to the number of MEFL.
Given two TU with identical promoters, RBS and terminators, but with different CDS, under the same experimental conditions, the number of molecules produced by each TU is the same.

Taking these axioms into account, the materials and procedure followed to calculate the MEFL/cell_sfGFP factor were as follows.

Materials:

Measuring equipment: Biotek Cytation3
96 well plate
MATLAB 2018a software

Procedure:

From the fluorescein spectrum, and the fluorescence data obtained from the Interlab Study experiment, we apply medium, gain and efficiency corrections .

The data are introduced in the Excel file of the Interlab Study. From the Fluorescein standard curve sheet we can obtain the MEFL/RFU factor, and then calculate the MEFL/cell factor.

Parameters	Value
Cuture volume (96 well plate)	200 μL
Cells per μL per OD unit	200 cells/OD μ L
MEFL/RFU factor	4,38.10¹⁰ MEFL/RFU
MEFL/cell_GFP factor	273.7 MEFL OD/RFU cell

We establish the protocol of the experiment and the parameters of our equipment. The experiment consists in the measurement of the absorbance and fluorescence of two TU with identical promoter, RBS and transcriptional terminator, but whose CDS codifies for GFP and sfGFP proteins. The most relevant information of the experiment is detailed in the following Table:

Parameters	Description
Time	06:00:00 (HH:MM:SS). Measurement interval: 05:00 (MM:SS)
Number of samples	8 samples for each TU. Total samples: 16
Temperature	37 ºC
Shake	Double Orbital Continuous
Absorbance. Optical Density (OD)	Wavelenght at 600 nm emission
Excitacion wavelength	485 nm
Emission wavelength	528 nm
Gain	60

The experiment is introduced into the equipment and the absorbance and fluorescence curves are obtained, as well as the curve of Fluorescence F/Absorbance OD or FOD ratio with the applied medium fluorescence correction.
The absorbance and fluorescence data are exported to an Excel file.
Running the MATLAB GFP_sfGFP_comparison.m script:
1. We extract the FOD data from the GFP and sfGFP from the Excel file.
2. We apply the gain and efficiency corrections of the spectrum to the FOD.
3. We plot the FOD curves and look for a stationary equilibrium time interval in the expression of the GFP and sfGFP proteins. In our experiment, we have decided to chose the interval [145,290] min
4. We obtain the average value of both regions.
5. We calculate the number of GFP MEFLs by multiplying the FOD data by the MEFL/cell_GFP factor.
6. If we assume that the number of molecules expressed by both TU is the same, the number of MEFL for the calculated GFP is the same as the number of MEFL for the sfGFP.
7. We calculate the MEFL/cell_sfGFP factor by dividing the MEFL number by the average FOD value of the sfGFP in the stationary region. In our case, MEFL/cell_sfGFP = 121.4 MEFL·OD/RFU·cell.

Reporter protein	MEFL/cell factor (MEFL·OD/RFU·cell)
GFP	273.7
sfGFP	121.4

With MEFL/cell factors, FOD data obtained in any experiment in which GFP or sfGFP was used as reporter proteins can be transformed into equivalent fluorescein molecules by applying the following ratio: Molecules = (MEFL/cell factor) · FOD.

The new data give us a more accurate estimation of the number of molecules in the cell. Consequently, by relating the experimental results with the theoretical mathematical models in the optimization process, the parameters of the model acquire values more consistent with their physical significance, working in equivalent molecules and not in arbitrary units.

@@ Line 1: / Line 1: @@
-{{Valencia_UPV}}
+{{Valencia_UPV/bootstrap.css}}
+<!--{{Valencia_UPV/bootstrapMin.css}}-->
+{{Valencia_UPV/custom.css}}
+{{Valencia_UPV/flickity.css}}
+{{Valencia_UPV/iconsmind.css}}
+{{Valencia_UPV/jquery.steps.css}}
+{{Valencia_UPV/lightbox.min.css}}
+{{Valencia_UPV/stackinterface.css}}
+{{Valencia_UPV/theme.css}}
+{{Valencia_UPV/navBar}}
 <html>
+<style>
-<div class="column full_size">
+#home_logo, #sideMenu { display:none; }
+#sideMenu, #top_title, .patrollink  {display:none;}
+#content { margin-left:0px; margin-top:-30px; padding:0px; width:100%;}
-<h1>Experiments</h1>
+</style>
-<p>Describe the research, experiments, and protocols you used in your iGEM project. These should be detailed enough for another team to repeat your experiments.</p>
-<p>
+    <head>
-Please remember to put all characterization and measurement data for your parts on the corresponding Registry part pages.
+        <meta charset="utf-8">
-</p>
+        <title>Stack Multipurpose HTML Template</title>
+        <meta name="viewport" content="width=device-width, initial-scale=1.0">
+        <meta name="description" content="Site Description Here">
+        <link href="https://fonts.googleapis.com/css?family=Open+Sans:200,300,400,400i,500,600,700%7CMerriweather:300,300i" rel="stylesheet">
+        <link href="https://fonts.googleapis.com/icon?family=Material+Icons" rel="stylesheet">
+        <link href="https://fonts.googleapis.com/css?family=Lato:400,700,900" rel="stylesheet">
+        <link href="https://2018.igem.org/Template:Valencia_UPV/fontAwesome?action=raw&ctype=text/css" rel="stylesheet" type="text/css">
-</div>
+    </head>
+  <body class="" style="background-color:#353535;">
+        <a class="anchorOffset" id="start"></a>
-<div class="column two_thirds_size">
-<h3>What should this page contain?</h3>
-<ul>
-<li> Protocols </li>
-<li> Experiments </li>
-<li> Documentation of the development of your project </li>
-</ul>
-</div>
+<!--------- NAVBAR-------------------------------------------------------------------------------------------------------------------------------------------------------------->
-<div class="column third_size">
-<div class="highlight decoration_A_full">
-<h3>Inspiration</h3>
-<ul>
-<li><a href="https://2014.igem.org/Team:Colombia/Protocols">2014 Colombia </a></li>
-<li><a href="https://2014.igem.org/Team:Imperial/Protocols">2014 Imperial </a></li>
+<!------------CONTENEDOR PRINCIPAL ---------------------------------------------------------------------------------------------------------------------------------->
-<li><a href="https://2014.igem.org/Team:Caltech/Project/Experiments">2014 Caltech </a></li>
-</ul>
+<div class="main-container">
+            <section style="/* background-color: #353535; *//* background-image: url(&quot;https://previews.123rf.com/images/bogdanhoda/bogdanhoda1504/bogdanhoda150400014/38730978-test-tubes-in-clinic-pharmacy-and-medical-research-laboratory-with-male-scientist-using-pipette.jpg&quot;); *//* background-size: 100%; */height:44.2em;margin-top: 0em;padding-top: 1.1em;" class="height-90 parallax">
+               <div class="main-container" style="
+height: 100%;">
+               <center style="
+height: 100%;">
+               <div class="col-md-10" style="
+height: 100%;margin-top: 4em;">
+<div id="logoPrinteria" class="item" style="
+background-image: url(http://via.placeholder.com/1155x666);
+height: 100%;
+width: 100%;background-attachment: fixed; background-position: center; background-position-y: 100px; background-repeat: no-repeat">
+    <div style="/* position: absolute; *//* top: 81%; *//* padding-left: 31em; */">
+<a class="btn down inner-link active" href="#story" style="font-size: 82%;/* position: fixed; *//* bottom: 7%; */top: 541px;z-index: 99;background-color: white;position: relative;border-radius: 80%;width: 3.8em;height: 3.8em;padding: 0;padding-top: 14px;">
+            <i class="stack-interface stack-down-open-big"></i>
+        </a>
 </div>
+ </div>
+</div>
+               </center>
 </div>
-<div class="clear"></div>
+              </section></div>
+<a class="anchorOffset" id="story"></a>
+<section class="feature-large" style="padding-top: 6em; padding-bottom: 50px; outline: none;" data-scroll-id="story" tabindex="-1">
+   <div class="container">
+      <div class="row">
+         <div class="col-md-3" style="
+            margin-left: 0em;
+            padding-right: 0px;
+            padding-left: 0.2em;">
+            <div class="tabs-container tabs--vertical page-navigator" style="position: -webkit-sticky;position: sticky;padding-bottom: 0px;margin-top: 11.4em;">
+               <h4 style="
+                  padding-left: 1em;
+                  ">Index</h4>
+               <ul class="lateral">
+                  <li class="lateral" style="margin-bottom: 0px;padding-bottom: 1em;padding-top: initial;padding-left: 1em;">
+                     <div class="tab__title" style="
+                        line-height: 1.3em;
+                        ">
+                        <a href="#lmHuman" class="lateral inner-link" style="
+                           color: #353535;
+                           opacity: 1;
+                           list-style-type: circle;
+                           ">Human Practices</a>
+                     </div>
+                  </li>
+                  <li class="lateral">
+                     <div class="tab__title" style="
+                        line-height: 1.3em;
+                        ">
+                        <a href="#lmSoftware" class="lateral inner-link" style="
+                           color: #353535;
+                           opacity: 1;
+                           ">Software</a>
+                     </div>
+                  </li>
+                  <li class="lateral">
+                     <div class="tab__title" style="
+                        line-height: 1.3em;
+                        ">
+                        <a href="#lmHardware" class="lateral inner-link" style="
+                           color: #353535;
+                           opacity: 1;
+                           ">Hardware</a>
+                     </div>
+                  </li>
+                  <li class="lateral">
+                     <div class="tab__title" style="
+                        line-height: 1.3em;
+                        ">
+                        <a href="#lmWetLab" class="lateral inner-link" style="
+                           color: #353535;
+                           opacity: 1;
+                           ">Wet Lab</a>
+                     </div>
+                  </li>
+                  <li class="lateral">
+                     <div class="tab__title" style="
+                        line-height: 1.3em;
+                        ">
+                        <a href="#lmModeling" class="lateral inner-link" style="
+                           color: #353535;
+                           opacity: 1;
+                           ">Modeling</a>
+                     </div>
+                  </li>
+                  <li class="lateral">
+                     <div class="tab__title" style="
+                        line-height: 1.3em;
+                        ">
+                        <a href="#lmGraphicDesign" class="lateral inner-link" style="
+                           color: #353535;
+                           opacity: 1;
+                           ">Graphic Design</a>
+                     </div>
+                  </li>
+                  <li class="lateral">
+                     <div class="tab__title" style="
+                        line-height: 1.3em;
+                        ">
+                        <a href="#lmWiki" class="lateral inner-link" style="
+                           color: #353535;
+                           opacity: 1;
+                           ">Wiki</a>
+                     </div>
+                  </li>
+                  <li class="lateral">
+                     <div class="tab__title" style="
+                        line-height: 1.3em;
+                        ">
+                        <a href="#lmOrgCom" class="lateral inner-link" style="
+                           color: #353535;
+                           opacity: 1;
+                           ">Organization / Communication</a>
+                     </div>
+                  </li>
+               </ul>
+               <ul class="tabs-content">
+               </ul>
+            </div>
+         </div>
+         <div class="col-md-9 col9Attr" style="padding-left: 6em;padding-right: 2em;">
+            <div>
+               <a class="anchorOffset" id="lmModeling"></a>
+               <h2 class="h2Modeling">Modeling</h2>
+               <p>
+               	<b>Do you think it is possible to mathematically describe a cell? Would you like to know the possibilities that Printeria modeling offers you?</b>
+               	</p>
+               	<p>
+				One of the fundamental bases of the Printeria project has undoubtedly been <b>mathematical modeling</b>. Thanks to the development and application of new mathematical models, it is possible to <b>quantify the expression of proteins</b> in cells, and therefore <b>characterize</b> through different experiments the parts designed by Printeria. From the Printeria modeling team, we intend to reach different goals:
+				</p>
+				<ul>
+				<li><p>
+					<b>Design simple mathematical models</b> based on differential equations that describe the biochemical processes of a cell. With them, we can simulate the different genetic circuits that Printeria allows us to build.
+				</p></li>
+				<li><p>
+					Develop a <b>Simulation Tool</b> that allows the user to visualize a prediction of the results of their experiment before running it in Printeria.
+				</p></li>
+				<li><p>
+					<b>Optimize model parameters</b> to match simulation results to experimental data obtained from Printeria constructions.
+				</p></li>
+				<li><p>
+					<b>Characterize the parts of our <a href="https://2018.igem.org/Team:Valencia_UPV/Part_Collection" target="_blank">Part Collection</a></b> from the optimization results and provide the user with all the information about the Printeria kit.
+				</p></li>
+				</ul>
+				<p>
+				Although in the development of the project we have dealt with all these aspects, all of them have a single purpose: to understand and <b>describe in a mathematical way the biological processes</b> that take place inside the cell.
+                </p>
+               <!--BOTONES PARA MOVERSE POR LA PÁGINA -->
+               <!--Los enlaces están bien pero los ids a donde apuntan están descuadrados adrede para que al pulsar lleve un poco por encima del contenido -->
+               <div id="" data-scroll-id="" tabindex="-1" style="outline: none;margin-top: 1em;margin-bottom: 1em;justify-content: center;" class="row">
+                  <div class="col-md-3" style="text-align: center;margin-right: 2em;">
+                     <a class="inner-link" href="#models"><img class="fotosModeling" src="https://static.igem.org/mediawiki/2018/c/ca/T--Valencia_UPV--instagramUPV2018.png"></a>
+                     <p style="text-align: center !important; font-weight: bold;">Models &amp; Experiments</p>
+                  </div>
+                  <div class="col-md-3" style="text-align: center;margin-right: 2em;">
+                     <a class="inner-link" href="#charact"><img class="fotosModeling" src="https://static.igem.org/mediawiki/2018/c/ca/T--Valencia_UPV--instagramUPV2018.png"></a>
+                     <p style="text-align: center !important; font-weight: bold;">Characterization procedures</p>
+                  </div>
+                  <div class="col-md-3" style="text-align: center;margin-right: 2em;">
+                     <a class="inner-link" href="#simtool"><img class="fotosModeling" src="https://static.igem.org/mediawiki/2018/c/ca/T--Valencia_UPV--instagramUPV2018.png"></a>
+                     <p style="text-align: center !important; font-weight: bold;">Simulation Tool</p>
+                  </div>
+               </div>
+               <!--FIN BOTONES PARA MOVERSE POR LA PÁGINA -->
+               <div>
+               	<a class="anchorOffset" id="models"></a>
+                  <h3>Models & Experiments</h3>
+                    <p><a href=""></a>
+                    The development of <b>new and simple mathematical models</b> has been one of the essential bases of the Printeria project. The models offer us a <b>multitude of applications</b>: to describe the basic biochemical reactions of the cell, to <a href="#charact" class="inner-link">characterize</a> the parts of our Printeria kit by means of different experiments or to elaborate <a href="#simtool" class="inner-link">new tools</a> that facilitate to our user the learning of Synthetic Biology.
+                    </p>
+                    <p>
+                    In Printeria Modeling team, we have grouped the models designed in <b>two fields</b>: <a href="#const_models" class="inner-link"><b>constitutive expression models</b></a> and <a href="#ind_models" class="inner-link"><b>inducible expression models</b></a>. We have also established an <a href="#exp_protocol" class="inner-link"><b>experimental protocol</b></a> and a <a href="#optimization" class="inner-link"><b>multi-objective optimization protocol</b></a> in order to contrast theoretical and experimental results. With them, we can describe practically all the cells modified with Printeria.
+                    </p>
+                    <p align="center">
+                    	<b>Want to find out more about how we did it?</b>
+                    </p>
+                    <a class="anchorOffset" id="design_process"></a>
+                    <h4>Modeling process</h4>
+                    <p>
+                    At the beginning of the design of any mathematical model, we have started from a <b>cell scheme</b> in which all the <b>cellular biochemical reactions</b> we wish to model are represented. From the reactions, a set of equations is deduced. Each of these equations describes the <b>temporal variation of the biochemical species</b> in the cell (DNA, RNA, proteins...), and depends on a set of <b>parameters with a physical meaning</b>.
+                    </p>
+                    <img src="https://static.igem.org/mediawiki/2018/7/7e/T--Valencia_UPV--esq_modelingUPV2018.png">
+                    <p>
+                    Because they are <b>differential and generally non-linear equations</b>, mathematical models are simulated using software tools such as MATLAB. The results obtained from the <b>simulation</b> reflect the evolution of the <b>concentration of biochemical species over time</b>. These theoretical results can be contrasted with the experimental data obtained in the laboratory and and thereby adjust the parameters so that the model adapts to the data in an optimal way. This process, called <b>multiobjective parameter optimization</b>, is the process that allows us to <b>validate our models</b> and see if they respond to what happens inside the cells.
+                    </p>
+                    <p>
+                    In the development of the set of mathematical models that describe the genetic circuits printed by Printeria, we have prioritized <b>simplicity over complexity</b> of the model. A more complex model may be more precise, but it requires a large number of parameters whose value is often unknown. Therefore, our models are composed of <b>few equations and parameters</b>, and they are all <b>deterministic models</b>, that is, their parameters take a concrete, non-probabilistic value.
+                    </p>
+                    <a class="anchorOffset" id="exp_protocol"></a>
+                    <h4>Experimental protocol</h4>
+                    <p>
+                    The Printeria Modeling and Lab team and have jointly designed an experimentation protocol for the laboratory experiments. Thanks to it, and from the colonies of the different UT the experimental data can be obtained, processed and ready to be optimized.
+                    </p>
+                    <p><i>Materials:</i></p>
+					<ul>
+						<li><a href=""></a>
+							<p>Printeria transcriptional units</b> (see our Printeria <a href="https://2018.igem.org/Team:Valencia_UPV/Part_Collection" target="_blank">Part Collection</a>)</p>
+						</li>
+						<li>
+							<p>Measuring equipment: <b>Biotek Cytation3</b></p>
+						</li>
+						<li>
+							<p><b>96 well plate</b></p>
+						</li>
+						<li>
+							<p><b>MATLAB 2018a software</b></p>
+						</li>
+					</ul>
+					<p><i>Protocol:</i></p>
+					<ol>
+						<li><p>
+							The TU colonies assembled by our Lab team are diluted in the medium (generally we have used LB or M9 medium with Kanamycin) so that they all have the same OD.
+						</p></li>
+						<li><p>
+							The experiment is designed on the measuring equipment. In our case, we use the Biotek Cytation3 equipment. We establish the equipment parameters.
+							<table style="width:100%">
+							  <tr>
+							    <th><p>Parameters</p></th>
+							    <th><p>Description</p></th>
+							  </tr>
+							  <tr>
+							    <td><p>Time</p></td>
+							    <td><p><b>06:00:00 (HH:MM:SS)</b> usually. <b>Measurement interval: 05:00 (MM:SS)</b></p></td>
+							  </tr>
+							  <tr>
+							    <td><p>Number of samples</p></td>
+							    <td><p>We normally set <b>8 samples of reporter protein for each TU colony</b></p></td>
+							  </tr>
+							  <tr>
+							    <td><p>Number of medium samples</p></td>
+							    <td><p>We normally set <b>8 samples of medium</b>. We normally use <b>LB or M9 medium with Kanamycin</b></p></td>
+							  </tr>
+							  <tr>
+							    <td><p>Temperature</p></td>
+							    <td><p>37 ºC</p></td>
+							  </tr>
+							  <tr>
+							    <td><p>Shake</p></td>
+							    <td><p>Double Orbital. Continuously. We shake the plate before each measure</p></td>
+							  </tr>
+  							  <tr>
+							    <td><p>Absorbance. Optical Density (OD) measure</p></td>
+							    <td><p>Wavelenght at <b>600 nm</b> emission</p></td>
+							  </tr>
+							  <tr>
+							    <td><p>Excitacion wavelength</p></td>
+							    <td><p>We normally set <b>485 nm</b></p></td>
+							  </tr>
+							  <tr>
+							    <td><p>Emission wavelength</p></td>
+							    <td><p>We normally set <b>528 nm</b></p></td>
+							  </tr>
+							  <tr>
+							    <td><p>Gain (G)</p></td>
+							    <td><p>Normally the gain value is <b>G = 60</b>, although for proteins with lower fluorescence, it is recommended that G takes higher values.</p></td>
+							  </tr>
+							</table>
+						</p></li>
+						<li><p>
+							Samples are introduced into the 96 well plate and the experiment begins.
+						</p></li>
+						<li><p>
+							After the experiment, fluorescence and absorbance data obtained are exported to an Excel file.
+						</p></li>
+						<li><p>
+							We run the MATLAB <a href="https://static.igem.org/mediawiki/2018/1/13/T--Valencia_UPV--convert_dataUPV2018.txt" target="_blank">convert_data.m</a> script. This script uses several additional MATLAB files with which:
+							<ol style="list-style-type: lower-latin;">
+								<li><p>
+									We extract the fluorescence and absorbance matrices from the Excel file (see <a href="https://static.igem.org/mediawiki/2018/e/e4/T--Valencia_UPV--readExperiment3UPV2018.txt" target="_blank">readExperiment3.m</a> MATLAB function).
+								</p></li>
+								<li><p>
+									We apply the <a href="#esquema_ops" class="inner-link">corrections</a> to the data (see <a href="https://static.igem.org/mediawiki/2018/3/32/T--Valencia_UPV--plotExperiment_ybUPV2018.txt" target="_blank">plotExperiment_yb.m</a> MATLAB function).
+								</p></li>
+								<li><p>
+									We save the data of FOD, OD, molecules and time in a <b>.mat format</b> file.
+								</p></li>
+							</ol>
+						</p></li>
+					</ol>
+                    <a class="anchorOffset" id="optimization"></a>
+                    <h4>Multi-objective parameter optimization</h4>
+                    <p>
+                    <b>To what extent is our model valid? What values should the parameters take? To what extent do the theoretical results resemble the experimental ones? Is our model able to explain the behaviour of cells in reality?</b>
+                    </p>
+                    <p>
+                    One of the most important phases in the field of Printeria Modeling has been the <b>validation of theoretical models</b> with experimental data. In the Printeria team we wanted to confirm that the models designed are consistent with our experiments. However, we know that the <b>parameters will not take a fixed value</b>, but vary depending on the genetic construction with which we experiment. With this idea in mind, the need to apply an optimization process arises.
+                    </p>
+                    <p>
+                    The process of <b>multiobjective parameter optimization</b> consists in finding those parameters of the model that best fit a particular objective, i.e., a particular set of experimental data. From a dataset, and applying a mathematical algorithm, we can find the optimal solutions, also known as <b>Pareto-optimal or Pareto front solutions</b>. These solutions are such that there are no more solutions that do not improve one objective without worsening the rest.
+                    </p>
+                    <p>
+                    From the Printeria Modeling team we have established a protocol that allows us to obtain the optimal solutions of the model parameters for any experiment:
+                    </p>
+                    <ol>
+                    	<li><p>
+                    		We define the objectives to be optimized. These will be the values of Fluorescence F/Unit of Absorbance OD or FOD, and the cellular growth (OD) of a Transcriptional Unit (TU). Therefore, if we have different <i>n</i> TU in an experiment, we will optimize <i>2n</i> objectives.
+                    	</p></li>
+                    	<li><p>
+                    		We choose the parameters to optimize and establish the range of values they can take.
+                    	</p></li>
+                    	<li><p>
+                    		We select the identification and validation data. The identification data are those used in the optimization process. The validation data are used to compare the experiment with the simulated theoretical model and the optimized parameters.
+                    	</p></li>
+                    	<li><p>
+                    		We define the cost function, which specifies the error function to be minimized, simulates the model with different parameter values, and calculates the error between the simulation data and the experimental identification data.
+                    	</p></li>
+                    	<li><p>
+                    		The mathematical optimization algorithm is executed. In our case we use the <a href="https://es.mathworks.com/matlabcentral/fileexchange/39215-multi-objective-differential-evolution-algorithm-with-spherical-pruning" target="_blank" style="padding-right: 0">Multi-objective Differential Evolution Algorithm with Spherical Pruning</a>. This algorithm uses the cost function to test different parameter values and thus search for optimal parameter values.
+                    	</p></li>
+                    	<li><p>
+                    		The Pareto set and Pareto Front diagrams and the different solutions of our optimized parameters for the different objectives are obtained.
+                    	</p></li>
+                    	<li><p>
+							Decision Making is carried out by the designer. The optimal parameters are chosen for our <i>2n</i> objectives.
+                    	</p></li>
+                    	<li><p>
+                    		The model is simulated with the optimized parameters and compared with the validation data.
+                    	</p></li>
+                    </ol>
+                    <p>
+                    Following this protocol, and particularizing it for each experiment, we have been able to determine the validity of our models and obtain the best parameters for our experiments.
+                    </p>
+                    <a class="anchorOffset" id="const_models"></a>
+                    <h4>Constituve expression models</h4>
+                    <p>
+                    Constitutive expression models are those in which the expression of the <b>cell's messenger RNA (mRNA) is not regulated</b>. The gene does not depend on any inducing biochemical species, so it will produce RNA molecules continuously, and therefore the translation of the protein will also be constant.
+                    </p>
+                    <p><i>Model design</i></p>
+                    <p>
+                    A simple scheme that summarizes the cellular transcription and translation processes in a constitutive expression model is shown in the following Figure. We only consider as biochemical species involved in reactions to mRNA and our protein of interest (PoI).
+                    </p>
+                    <img src="https://static.igem.org/mediawiki/2018/5/54/T--Valencia_UPV--const_schemaUPV2018.png">
+                    <p>
+                    From the cell schema we can deduce the following biochemical reactions.
+                    </p>
+                    <img src="https://static.igem.org/mediawiki/2018/0/08/T--Valencia_UPV--react_constUPV2018.png">
+                    <table style="width:100%">
+						<tr>
+						  <th><p>Parameter</p></th>
+						  <th><p>Description</p></th>
+						  <th><p>Units</p></th>
+						  <th><p>Value</p></th>
+						</tr>
+						<tr>
+						  <td><p>C<sub>R</sub></p></td>
+						  <td><p>Constitutive transcription rate is calculated by multiplying the transcription rate per plasmid K<sub>R</sub> by the mean number of plasmids in cell c<sub>n</sub> . In our case, we are using the pMBI replication origin, so c<sub>n</sub> <meta charset="utf-8">&sim; 500, and C<sub>R</sub> = K<sub>R</sub>·c<sub>n</sub>.</p></td>
+						  <td><p>molecules.min<sup>-1</sup></p></td>
+						  <td><p><a href="">---</a></p></td>
+						 </tr>
+						 <tr>
+						  <td><p>p</p></td>
+						  <td><p>Translation rate</p></td>
+						  <td><p>min<sup>-1</sup></p></td>
+						  <td><p><a href="">---</a></p></td>
+						 </tr>
+						 <tr>
+						  <td><p>d<sub>m</sub></p></td>
+						  <td><p>mRNA degradation rate</p></td>
+						  <td><p>min<sup>-1</sup></p></td>
+						  <td><p><a href="">---</a></p></td>
+						 </tr>
+						 <tr>
+						  <td><p>d<sub>p</sub></p></td>
+						  <td><p>PoI degradation rate</p></td>
+						  <td><p>min<sup>-1</sup></p></td>
+						  <td><p><a href="">---</a></p></td>
+						 </tr>
+						 <tr>
+						  <td><p><meta charset="utf-8">&mu;</p></td>
+						  <td><p>Dilution rate</p></td>
+						  <td><p>min<sup>-1</sup></p></td>
+						  <td><p><a href="">---</a></p></td>
+						 </tr>
+						 <tr>
+						  <td><p>K<sub>max</sub></p></td>
+						  <td><p>Maximum growth capacity</p></td>
+						  <td><p>cells</p></td>
+						  <td><p><a href="">---</a></p></td>
+						 </tr>
+					</table>
+                    <p>
+                    With the reactions and approximations already established, we apply the <a href="">Law of Mass Action kynetics</a> (LMA) to deduce the differential equations. The LMA establishes that the variation of the species resulting from a reaction is proportional to the product of the reactants. In our case, to the set of differential equations deduced from the reactions we add a cell growth equation. Biochemical species and our model of equations obtained are shown below.
+                    </p>
+                    <table style="width:100%">
+						<tr>
+						  <th><p>Variable</p></th>
+						  <th><p>Biochemical species</p></th>
+						  <th><p>Units</p></th>
+						</tr>
+						 <tr>
+						  <td><p>n<sub>1</sub></p></td>
+						  <td><p>mRNA</p></td>
+						  <td><p>Molecules</p></td>
+						 </tr>
+						 <tr>
+						  <td><p>n<sub>2</sub></p></td>
+						  <td><p>PoI</p></td>
+						  <td><p>Molecules</p></td>
+						 </tr>
+						 <tr>
+						  <td><p>n<sub>3</sub></p></td>
+						  <td><p>Number of cells</p></td>
+						  <td><p>Cells</p></td>
+						 </tr>
+					</table>
+                    <img src="https://static.igem.org/mediawiki/2018/e/e1/T--Valencia_UPV--eq_constUPV2018.png">
+                    <p>
+                    By obtaining the equations, and seeking the simplicity of the model, a series of <a href="">quasi-static approximations</a> have been made. We have established that mRNA is generated with a constant C<sub>R</sub> constitutive transcription rate. Another consideration has been to assume that the copy number of plasmids in cell is constant. We have also proposed that other species such as ribosomes, which also participate in transcription, are constant over time.
+                    </p>
+                    <p>
+                    After developing the model and raising the reactions and equations, two files have been programmed in MATLAB: the <a href="https://static.igem.org/mediawiki/2018/e/e1/T--Valencia_UPV--mc_simpleUPV2018.txt" target="_blank">mc_simple.m</a> function, which describes the system of differential equations, and the script <a href="https://static.igem.org/mediawiki/2018/4/4e/T--Valencia_UPV--model_const_simpleUPV2018.txt" target="_blank">model_const_simple.m</a> in which the parameters are defined and the model is solved by the <a href="https://es.mathworks.com/help/matlab/ref/ode45.html?searchHighlight=ode45&s_tid=doc_srchtitle" style="padding-right: 0" target="_blank">ode45</a> MATLAB function. The result of the simulation of the model is shown in the following Figure.
+                	</p>
+                    <img src="https://static.igem.org/mediawiki/2018/f/f1/T--Valencia_UPV--model_const_sim_graphUPV2018.png">
+                    <p>
+                    In the simulation of mRNA and PoI we distinguish two phases in the temporal evolution of the species: the transient phase and the stationary phase. If we perform different simulations varying some parameters of the model, we can observe how the equilibrium value of the stationary phase varies. MATLAB script <a href="https://static.igem.org/mediawiki/2018/4/4d/T--Valencia_UPV--mc_sim_analysisUPV2018.txt" target="_blank">mc_sim_analysis.m</a> repeatedly simulates the model by varying the C<sub>R</sub> and p parameters, and graphs the results. The following Figure shows how an increase in both ratios implies a greater expression of the protein.
+                    </p>
+                    <img src="https://static.igem.org/mediawiki/2018/f/f9/T--Valencia_UPV--mc_sim_analysis_graphUPV2018.png">
+                    <p>
+                    At this point anyone could ask us: <b>Which are the advantages of a constitutive expression model like this? Why is a model like this really useful to us?</b> Faced with these questions, the Printeria Modeling team has looked for the answers in the great characteristics that the model offers us.
+                    </p>
+                    <ul>
+                    	<li><p>
+                    		A simple, <b>easy-to-understand model</b> that clearly explains the processes of cellular transcription and translation.
+                    	</p></li>
+                    	<li><p>
+                    		It is <b>valid for any Printeria construction</b> that has a constitutive expression promoter. The variation between different constructions occurs in the parameters values, not in the model equations!
+                    	</p></li>
+                    	<li><p>
+							It has <b>few parameters</b>, with a clear physical meaning, and <b>easily optimizable</b>.
+						</p></li>
+                    	<li><p>
+                    		As a compact model, <b>simulations and optimizations</b> are performed <b>at high speed</b>.
+                    	</p></li>
+                    </ul>
+                    <p><i>Experiments & optimization</i></p>
+                    <p>
+                    The simple model of constitutive expression is a model that characterizes a large number of Printeria Transcriptional Units (TU), specifically, all those constructions whose promoter is constitutive. We know that all these TU follow the <b>same model</b>, but each one has <b>different parameters</b>, and therefore different experimental values of fluorescence and absorbance. Faced with this situation, the Printeria Lab and Modeling teams have designed some experiments changing <a href="#RBS_exp" class="inner-link">RBS</a> and <a href="#prom_exp" class="inner-link">promoters</a>, in which we have applied the <a href="#optimization" class="inner-link">multi-objective parameter optimization</a>, and thus check whether the designed constitutive expression model responds to experimental reality.
+                	</p>
+                	<table style="width:100%">
+						  <tr>
+						    <th><p>Experiment parameters</p></th>
+						    <th><p>Description</p></th>
+						  </tr>
+						  <tr>
+						    <td><p>Time</p></td>
+						    <td><p>06:00:00 (HH:MM:SS). Measurement interval: 05:00 (MM:SS)</p></td>
+						  </tr>
+						  <tr>
+						    <td><p>Number of samples</p></td>
+						    <td><p>8 samples for each TU colony</p></td>
+						  </tr>
+						  <tr>
+						    <td><p>Number of medium samples</p></td>
+						    <td><p>8 samples of M9 medium</p></td>
+						  </tr>
+						  <tr>
+						    <td><p>Temperature</p></td>
+						    <td><p>37 ºC</p></td>
+						  </tr>
+						  <tr>
+						    <td><p>Shake</p></td>
+						    <td><p>Double Orbital. Continuously</p></td>
+						  </tr>
+							  <tr>
+						    <td><p>Absorbance. Optical Density (OD) measure</p></td>
+						    <td><p>Wavelenght at 600 nm emission</p></td>
+						  </tr>
+						  <tr>
+						    <td><p>Excitacion wavelength</p></td>
+						    <td><p>Wavelenght at 485 nm</p></td>
+						  </tr>
+						  <tr>
+						    <td><p>Emission wavelength</p></td>
+						    <td><p>Wavelenght at  528 nm</p></td>
+						  </tr>
+						  <tr>
+						    <td><p>Gain (G)</p></td>
+						    <td><p>60</p></td>
+						  </tr>
+					</table>
+                	<a class="anchorOffset" id="RBS_exp"></a>
+                	<h7><b>Experiments changing RBS</b></h7>
+                	<p>
+                	We have designed <b>two experiments</b> following the same <a href="#exp_protocol" class="inner-link">experimental protocol</a>. In them we have assembled different Printeria TU with the same promoters, CDS (<a href="http://parts.igem.org/Part:BBa_K2656013" target="_blank">sfGFP</a> reporter protein) and transcriptional terminator, but with <b>different RBS</b>. After obtaining the results, and following the <a href="#optimization" class="inner-link">optimization protocol</a>, we have obtained the parameters of the model and have validated our model.
+                	</p>
+                	<a class="anchorOffset" id="RBS_list"></a>
+                	<p><b>Printeria RBS</b>:</p>
+                	<ul>
+                		<li><p>
+                			<b>Strong expression</b>: <a href="http://parts.igem.org/Part:BBa_K2656009" target="_blank" style="padding-right: 0">BBa_K2656009</a>.
+                		</p></li>
+                		<li><p>
+                			<b>Medium expression</b>: <a href="http://parts.igem.org/Part:BBa_K2656011" target="_blank" style="padding-right: 0">BBa_K2656011</a>.
+                		</p></li>
+                		<li><p>
+                			<b>Low expression</b>: <a href="http://parts.igem.org/Part:BBa_K2656010" target="_blank" style="padding-right: 0">BBa_K2656010</a>.
+                		</p></li>
+                		<li><p>
+                			<b>Very low expression</b>: <a href="http://parts.igem.org/Part:BBa_K2656008" target="_blank" style="padding-right: 0">BBa_K2656008</a>, <a href="http://parts.igem.org/Part:BBa_K2656012" target="_blank" style="padding-right: 0">BBa_K2656012</a>.
+                		</p></li>
+                	</ul>
+					<table style="width:100%">
+						  <tr>
+						    <th><p>Optimization specifications</p></th>
+						    <th><p>Description</p></th>
+						  </tr>
+						  <tr>
+						    <td><p>Parameters</p></td>
+						    <td>
+						    	<p>Constitutive transcription rate C<sub>R</sub>: fixed</p>
+						    	<p><b>Translation rate p: to optimize</b></p>
+						    	<p><b>PoI degradation rate d<sub>p</sub>: to optimize</b></p>
+						    	<p>mRNA degradation rate d<sub>m</sub>: fixed</p>
+						    	<p><b>Dilution rate <meta charset="utf-8">&mu;: to optimize</b></p>
+						    	<p>Maximum growth capacity K<sub>max</sub>: experimental value</p>
+						    </td>
+						  </tr>
+						  <tr>
+						    <td><p>Objetives to optimize</p></td>
+						    <td>
+						    	<p>For each TU, we set 2 objetives: FOD (Fluorescence) and OD (Absorbance). In each experiment we have measured 3 TU, so we must set <b>6 objectives par experiment</b></p>
+						    </td>
+						  <tr>
+						    <td><p>Parameter ranges</p></td>
+						    <td>
+						    	<p>Translation rate p: [0.001 - 6] min<sup>-1</sup> </p>
+						    	<p>PoI degradation rate d<sub>p</sub>: [0.0058 - 0.0087] min<sup>-1</sup> </p>
+						    	<p>Dilution rate <meta charset="utf-8">&mu;: [0.0058 - 0.035] min<sup>-1</sup> </p>
+						    </td>
+						  </tr>
+						  <tr>
+						    <td><p>MATLAB files</p></td>
+						    <td>
+						    	<p><b>spMODEparam.m</b>:  it defines the parameters to be optimized, their value ranges, number of objectives, the cost function, the identification and validation experimental data and other spMODE algorithm parameters. When the script is executed, a spMODEDat structure variable is defined. This structure contains all the standardized optimization information that the spMODE algorithm will need for its execution.</p>
+						    	<ul>
+						    		<li><p>
+						    			Experiment 1: <a href="https://static.igem.org/mediawiki/2018/c/c9/T--Valencia_UPV--spMODEparam_exp1sfGFPUPV2018.txt" target="_blank">spMODEparam_exp1sfGFP.m</a>
+						    		</p></li>
+						    		<li><p>
+						    			Experiment 2: <a href="https://static.igem.org/mediawiki/2018/0/0d/T--Valencia_UPV--spMODEparam_exp2sfGFPUPV2018.txt" target="_blank">spMODEparam_exp2sfGFP.m</a>
+						    		</p></li>
+						    	</ul>
+						    	<p><a href="https://static.igem.org/mediawiki/2018/4/45/T--Valencia_UPV--spMODEUPV2018.txt"><b>spMODE.m</b></a>: it contains the <a href="https://es.mathworks.com/matlabcentral/fileexchange/39215-multi-objective-differential-evolution-algorithm-with-spherical-pruning" target="_blank" style="padding-right: 0">Multi-objective Differential Evolution Algorithm with Spherical Pruning</a>, which optimizes parameters for our experimental results (to execute spMODE.m file we need also <a href="https://static.igem.org/mediawiki/2018/9/99/T--Valencia_UPV--SphPruningUPV2018.txt">SphPruning.m</a> file).</p>
+						    	<p><a href="https://static.igem.org/mediawiki/2018/9/95/T--Valencia_UPV--CostFunction_RBS_RMS_2nUPV2018.txt" target="_blank"><b>CostFunction_RBS_RMS_2n.m</b></a>: it simulates the model with different vector parameters, and calculates the Root Mean Square Error with the identification dataset. Then, it returns the error for each objective and parameter vector.</p>
+						    	<p><a href="https://static.igem.org/mediawiki/2018/e/e4/T--Valencia_UPV--levelDiagramUPV2018.txt" target="_blank"><b>levelDiagram.m</b></a>: plots Pareto front and Pareto set that gives us back the algorithm.</p>
+						    	<p><a href="https://static.igem.org/mediawiki/2018/3/3f/T--Valencia_UPV--execute_RBS_2nUPV2018.txt" target="_blank"><b>execute_RBS_2n.m</b></a>: this script launches the optimization by executing spMODEparam.m, spMODEm algorithm and levelDiagram.m. Then, allows the user to enter the best parameters, and plots validation experimental data and simulation results in the same graph.</p>
+						    </td>
+						  </tr>
+					</table>
+					<img src="https://static.igem.org/mediawiki/2018/b/b1/T--Valencia_UPV--optimization_exp1_RBS_graphUPV2018.png">
+					<img src="https://static.igem.org/mediawiki/2018/b/b1/T--Valencia_UPV--optimization_exp1_RBS_graphUPV2018.png">
+					<img src="https://static.igem.org/mediawiki/2018/b/b1/T--Valencia_UPV--optimization_exp1_RBS_graphUPV2018.png">
+					<p>
+					The results of the graphs have been obtained after an optimization and decision-making process, in which the optimal parameters have been selected and the constituent model for each objective has been simulated and compared with the experimental data. The optimized parameters are summarized in the following table.
+					</p>
+					<table style="width:100%">
+						  <tr>
+						    <th><p>Optimized parameters</p></th>
+						    <th><p>Values</p></th>
+						  </tr>
+						  <tr>
+						  	<td><p>Translation rate <b>p</b></p></td>
+						  	<td>
+						  		<p><b>Experiment 1</b>:</p>
+						  		<ul>
+						  			<li><p>
+						  				For strong RBS <a href="http://parts.igem.org/Part:BBa_K2656009" target="_blank" style="padding-right: 0">BBa_K2656009</a>: p = 0.2059 min<sup>-1</sup>
+						  			</p></li>
+						  			<li><p>
+						  				For medium RBS <a href="http://parts.igem.org/Part:BBa_K2656011" target="_blank" style="padding-right: 0">BBa_K2656011</a>: p = 0.082 min<sup>-1</sup>
+						  			</p></li>
+						  			<li><p>
+						  				For low RBS <a href="http://parts.igem.org/Part:BBa_K2656010" target="_blank" style="padding-right: 0">BBa_K2656010</a>: p = 0.0072 min<sup>-1</sup>
+						  			</p></li>
+						  		</ul>
+						  		<p><b>Experiment 2</b>:</p>
+						  		<ul>
+						  			<li><p>
+						  				For strong RBS <a href="http://parts.igem.org/Part:BBa_K2656009" target="_blank" style="padding-right: 0">BBa_K2656009</a>: p = - min<sup>-1</sup>
+						  			</p></li>
+						  			<li><p>
+						  				For very low RBS <a href="http://parts.igem.org/Part:BBa_K2656008" target="_blank" style="padding-right: 0">BBa_K2656008</a>: p = - min<sup>-1</sup>
+						  			</p></li>
+						  			<li><p>
+						  				For very low RBS <a href="http://parts.igem.org/Part:BBa_K2656012" target="_blank" style="padding-right: 0">BBa_K2656012</a>: p = - min<sup>-1</sup>
+						  			</p></li>
+						  		</ul>
+						  	</td>
+						  </tr>
+						  <td><p>PoI degradation rate <b>d<sub>p</sub></b></p></td>
+						  	<td>
+						  		<p><b>Experiment 1</b>: d<sub>p</sub> = 0.0058 min<sup>-1</sup></p>
+						  		<p><b>Experiment 2</b>: d<sub>p</sub> = - min<sup>-1</sup></p>
+						  	</td>
+						  </tr>
+						  <tr>
+						  	<td><p>Dilution rate <b><meta charset="utf-8">&mu;</b></p></td>
+						  	<td>
+						  		<p><b>Experiment 1</b>:</p>
+						  		<ul>
+						  			<li><p>
+						  				For strong RBS <a href="http://parts.igem.org/Part:BBa_K2656009" target="_blank" style="padding-right: 0">BBa_K2656009</a>: <meta charset="utf-8">&mu; = 0.01492 min<sup>-1</sup>
+						  			</p></li>
+						  			<li><p>
+						  				For medium RBS <a href="http://parts.igem.org/Part:BBa_K2656009" target="_blank" style="padding-right: 0">BBa_K2656011</a>: <meta charset="utf-8">&mu; = 0.01631 min<sup>-1</sup>
+						  			</p></li>
+						  			<li><p>
+						  				For low RBS <a href="http://parts.igem.org/Part:BBa_K2656009" target="_blank" style="padding-right: 0">BBa_K2656010</a>: <meta charset="utf-8">&mu; = 0.0188 min<sup>-1</sup>
+						  			</p></li>
+						  		</ul>
+						  		<p><b>Experiment 2</b>:</p>
+						  		<ul>
+						  			<li><p>
+						  				For strong RBS <a href="http://parts.igem.org/Part:BBa_K2656009" target="_blank" style="padding-right: 0">BBa_K2656009</a>: <meta charset="utf-8">&mu; = - min<sup>-1</sup>
+						  			</p></li>
+						  			<li><p>
+						  				For very low RBS <a href="http://parts.igem.org/Part:BBa_K2656009" target="_blank" style="padding-right: 0">BBa_K2656008</a>: <meta charset="utf-8">&mu; = - min<sup>-1</sup>
+						  			</p></li>
+						  			<li><p>
+						  				For very low RBS <a href="http://parts.igem.org/Part:BBa_K2656009" target="_blank" style="padding-right: 0">BBa_K2656012</a>: <meta charset="utf-8">&mu; = - min<sup>-1</sup>
+						  			</p></li>
+						  		</ul>
+						  	</td>
+					</table>
+					<p>
+					In addition to simulating the optimized models, the Printeria Modeling team has also calculated the relative force between the different RBS, taking <a href="http://parts.igem.org/Part:BBa_K2656009" target="_blank" style="padding-right: 0">BBa_K2656009</a> strong RBS as a reference. The relative force has been defined as the quotient between the values of the protein in equilibrium of the results of the simulation of one RBS and another reference RBS. Likewise, a ratio between <i>p</i> parameters of the different <a href="#RBS_list" class="inner-link">RBS parts</a> and <i>p</i> parameter of the reference RBS has been calculated. The characterization of the RBS parts by their relative strength is shown below.
+					<img src="https://static.igem.org/mediawiki/2018/d/d7/T--Valencia_UPV--rel_force_eqUPV2018.png">
+					</p>
+					<table style="width:100%">
+						  <tr>
+						    <th><p>Printeria RBS part</p></th>
+						    <th><p>Relative strength</p></th>
+						    <th><p>p parameter ratio (p<sub>RBS</sub>/p<sub>ref</sub>)</p></th>
+						  </tr>
+						  <tr>
+						    <td><p><a href="http://parts.igem.org/Part:BBa_K2656009" target="_blank" style="padding-right: 0">BBa_K2656009</a> (<i>Reference</i>)</p></td>
+						    <td><p>1</p></td>
+						    <td><p>1</p></td>
+						  </tr>
+						  <tr>
+						    <td><p><a href="http://parts.igem.org/Part:BBa_K2656011" target="_blank" style="padding-right: 0">BBa_K2656011</a></p></td>
+						    <td><p>0.371</p></td>
+						    <td><p>0.398</p></td>
+						  </tr>
+						  <tr>
+						    <td><p><a href="http://parts.igem.org/Part:BBa_K2656010" target="_blank" style="padding-right: 0">BBa_K2656010</a></p></td>
+						    <td><p>0.029</p></td>
+						    <td><p>0.035</p></td>
+						  </tr>
+						  <tr>
+						    <td><p><a href="http://parts.igem.org/Part:BBa_K2656008" target="_blank" style="padding-right: 0">BBa_K2656008</a></p></td>
+						    <td><p>-</p></td>
+						    <td><p>-</p></td>
+						  </tr>
+						  <tr>
+						    <td><p><a href="http://parts.igem.org/Part:BBa_K2656012" target="_blank" style="padding-right: 0">BBa_K2656012</a></p></td>
+						    <td><p>-</p></td>
+						    <td><p>-</p></td>
+						  </tr>
+					</table>
+					<a class="anchorOffset" id="prom_exp"></a>
+                	<h7><b>Experiments changing promoters</b></h7>
+                	<p>
+                	We have designed <b>two experiments</b> following the same <a href="#exp_protocol" class="inner-link">experimental protocol</a>. In them we have assembled different Printeria TU with the same RBS, CDS (<a href="http://parts.igem.org/Part:BBa_K2656022" target="_blank">GFP</a> reporter protein) and transcriptional terminator, and with <b>different promoters</b>. After obtaining the results, and following the <a href="#optimization" class="inner-link">optimization protocol</a>, we have obtained the parameters of the model and have validated our model.
+                	</p>
+                	<a class="anchorOffset" id="promoter_list"></a>
+                	<p><b>Printeria promoters</b>:</p>
+                	<ul>
+                		<li><p>
+                			<b>Strong promoters</b>: <a href="http://parts.igem.org/Part:BBa_K2656005" target="_blank" style="padding-right: 0">BBa_K2656005</a>
+                		</p></li>
+                		<li><p>
+                			<b>Medium promoters</b>: <a href="http://parts.igem.org/Part:BBa_K2656007" target="_blank" style="padding-right: 0">BBa_K2656007</a>
+                		</p></li>
+                		<li><p>
+                			<b>Low promoters</b>: <a href="http://parts.igem.org/Part:BBa_K2656004" target="_blank" style="padding-right: 0">BBa_K2656004</a>
+                		</p></li>
+                	</ul>
+					<table style="width:100%">
+						  <tr>
+						    <th><p>Optimization specifications</p></th>
+						    <th><p>Description</p></th>
+						  </tr>
+						  <tr>
+						    <td><p>Parameters</p></td>
+						    <td>
+						    	<p><b>Constitutive transcription rate C<sub>R</sub>: to optimize</b></p>
+						    	<p>Translation rate p: fixed</p>
+						    	<p><b>PoI degradation rate d<sub>p</sub>: to optimize</b></p>
+						    	<p>mRNA degradation rate d<sub>m</sub>: fixed</p>
+						    	<p><b>Dilution rate <meta charset="utf-8">&mu;: to optimize</b></p>
+						    	<p>Maximum growth capacity K<sub>max</sub>: experimental value</p>
+						    </td>
+						  </tr>
+						  <tr>
+						    <td><p>Objetives to optimize</p></td>
+						    <td>
+						    	<p>For each TU, we set 2 objetives: FOD (Fluorescence) and OD (Absorbance). In each experiment we have measured 3 TU, so we must set <b>6 objectives par experiment</b></p>
+						    </td>
+						  <tr>
+						    <td><p>Parameter ranges</p></td>
+						    <td>
+						    	<p>Constitutive transcription rate C<sub>R</sub>: [ - ] min<sup>-1</sup> </p>
+						    	<p>PoI degradation rate d<sub>p</sub>: [0.0058 - 0.0087] min<sup>-1</sup> </p>
+						    	<p>Dilution rate <meta charset="utf-8">&mu;: [0.0058 - 0.035] min<sup>-1</sup> </p>
+						    </td>
+						  </tr>
+						  <tr>
+						    <td><p>MATLAB files</p></td>
+						    <td>
+						    	<p><b>spMODEparam.m</b>:  it defines the parameters to be optimized, their value ranges, number of objectives, the cost function, the identification and validation experimental data and other spMODE algorithm parameters. When the script is executed, a spMODEDat structure variable is defined. This structure contains all the standardized optimization information that the spMODE algorithm will need for its execution.</p>
+						    	<ul>
+						    		<li><p>
+						    			Experiment 1: <a href="" target="_blank">--</a>
+						    		</p></li>
+						    	</ul>
+						    	<p><a href="https://static.igem.org/mediawiki/2018/4/45/T--Valencia_UPV--spMODEUPV2018.txt"><b>spMODE.m</b></a>: it contains the <a href="https://es.mathworks.com/matlabcentral/fileexchange/39215-multi-objective-differential-evolution-algorithm-with-spherical-pruning" target="_blank" style="padding-right: 0">Multi-objective Differential Evolution Algorithm with Spherical Pruning</a>, which optimizes parameters for our experimental results (to execute spMODE.m file we need also <a href="https://static.igem.org/mediawiki/2018/9/99/T--Valencia_UPV--SphPruningUPV2018.txt">SphPruning.m</a> file).</p>
+						    	<p><a href="t" target="_blank"><b>CostFunction_Prom_RMS_2n.m</b></a>: it simulates the model with different vector parameters, and calculates the Root Mean Square Error with the identification dataset. Then, it returns the error for each objective and parameter vector.</p>
+						    	<p><a href="https://static.igem.org/mediawiki/2018/e/e4/T--Valencia_UPV--levelDiagramUPV2018.txt" target="_blank"><b>levelDiagram.m</b></a>: plots Pareto front and Pareto set that gives us back the algorithm.</p>
+						    	<p><a href="" target="_blank"><b>execute.m</b></a>: this script launches the optimization by executing spMODEparam.m, spMODEm algorithm and levelDiagram.m. Then, allows the user to enter the best parameters, and plots validation experimental data and simulation results in the same graph.</p>
+						    </td>
+						  </tr>
+					</table>
+					<img src=""> GRÁFICAS RESULTADOS 2 EXPERIMENTOS PROMOTORES
+                    <p>
+					The results of the graphs have been obtained after an optimization and decision-making process, in which the optimal parameters have been selected and the constituent model for each objective has been simulated and compared with the experimental data. The optimized parameters are summarized in the following table.
+					</p>
+					<table style="width:100%">
+						  <tr>
+						    <th><p>Optimized parameters</p></th>
+						    <th><p>Values</p></th>
+						  </tr>
+						  <tr>
+						  	<td><p>Constitutive transcription rate <b>C<sub>R</sub></b></p></td>
+						  	<td>
+						  		<ul>
+						  			<li><p>
+						  				For strong promoter <a href="http://parts.igem.org/Part:BBa_K2656005" target="_blank" style="padding-right: 0">BBa_K2656005</a>: C<sub>R</sub> = - min<sup>-1</sup>
+						  			</p></li>
+						  			<li><p>
+						  				For medium promoter <a href="http://parts.igem.org/Part:BBa_K2656007" target="_blank" style="padding-right: 0">BBa_K2656007</a>: C<sub>R</sub> = - min<sup>-1</sup>
+						  			</p></li>
+						  			<li><p>
+						  				For low promoter <a href="http://parts.igem.org/Part:BBa_K2656004" target="_blank" style="padding-right: 0">BBa_K2656004</a>: C<sub>R</sub> = - min<sup>-1</sup>
+						  			</p></li>
+						  		</ul>
+						  	</td>
+						  </tr>
+						  <td><p>PoI degradation rate <b>d<sub>p</sub></b></p></td>
+						  	<td>
+						  		<p><b>Experiment 1</b>: d<sub>p</sub> = - min<sup>-1</sup></p>
+						  	</td>
+						  </tr>
+						  <tr>
+						  	<td><p>Dilution rate <b><meta charset="utf-8">&mu;</b></p></td>
+						  	<td>
+						  		<ul>
+						  			<li><p>
+						  				For strong promoter <a href="http://parts.igem.org/Part:BBa_K2656005" target="_blank" style="padding-right: 0">BBa_K2656005</a>: <meta charset="utf-8">&mu; = - min<sup>-1</sup>
+						  			</p></li>
+						  			<li><p>
+						  				For medium promoter <a href="http://parts.igem.org/Part:BBa_K2656007" target="_blank" style="padding-right: 0">BBa_K2656007</a>: <meta charset="utf-8">&mu; = - min<sup>-1</sup>
+						  			</p></li>
+						  			<li><p>
+						  				For low promoter <a href="http://parts.igem.org/Part:BBa_K2656004" target="_blank" style="padding-right: 0">BBa_K2656004</a>: <meta charset="utf-8">&mu; = - min<sup>-1</sup>
+						  			</p></li>
+						  		</ul>
+						  	</td>
+					</table><p>
+					In addition to simulating the optimized models, the Printeria Modeling team has also calculated the relative force between the different promoters, taking <a href="http://parts.igem.org/Part:BBa_K2656005" target="_blank" style="padding-right: 0">BBa_K2656005</a> promoter as a reference. Likewise, a ratio between <i>C<sub>R</sub></i> parameters of the different <a href="#promoter_list" class="inner-link">promoter parts</a> and <i>C<sub>R</sub></i> parameter of the reference promoter has been calculated. The characterization of the promoter parts by their relative strength is shown below.
+					</p>
+					<table style="width:100%">
+						  <tr>
+						    <th><p>Printeria RBS part</p></th>
+						    <th><p>Relative strength</p></th>
+						    <th><p>p parameter ratio (p<sub>RBS</sub>/p<sub>ref</sub>)</p></th>
+						  </tr>
+						  <tr>
+						    <td><p><a href="http://parts.igem.org/Part:BBa_K2656009" target="_blank" style="padding-right: 0">BBa_K2656009</a> (<i>Reference</i>)</p></td>
+						    <td><p>1</p></td>
+						    <td><p>1</p></td>
+						  </tr>
+						  <tr>
+						    <td><p><a href="http://parts.igem.org/Part:BBa_K2656011" target="_blank" style="padding-right: 0">BBa_K2656011</a></p></td>
+						    <td><p>-</p></td>
+						    <td><p>-</p></td>
+						  </tr>
+						  <tr>
+						    <td><p><a href="http://parts.igem.org/Part:BBa_K2656010" target="_blank" style="padding-right: 0">BBa_K2656010</a></p></td>
+						    <td><p>-</p></td>
+						    <td><p>-</p></td>
+						  </tr>
+					</table>
+                    <a class="anchorOffset" id="ind_models"></a>
+                    <h4>Inducible expression models</h4>
+					<p>
+					Inducible expression models are those in which the expression of the <b>cell's messenger RNA (mRNA) depends on one or more inducing molecules</b>, which can act as <b>activators or repressors</b>. The concentration of the inducing species determines the level of expression of mRNA and protein in the cell.</p>
+					<p>
+					Our Printeria device gives us the possibility to create some inducible genetic constructs. Of all of them, we have managed to design and experiment with some models...
+					</p>
+					<ul>
+						<li><p>
+							<a href="#pBAD" class="inner-link">P<sub>BAD</sub>/araC inducible model</a>
+						</p></li>
+						<li><p>
+							<a href="#LuxR_LuxI" class="inner-link">LuxR/LuxI inducible model</a>
+						</p></li>
+					</ul>
+					<a class="anchorOffset" id="pBAD"></a>
+                    <h5>P<sub>BAD</sub>/araC inducible model</h5>
+					<p>
+					The inducible promoter <a href="" class="inner-link"><b>P<sub>BAD</sub>/araC</b></a> is a promoter whose expression depends on two inducing molecules: the <b>L-arabinose</b> (activator) and the <b>araC</b> (repressor). The dimeric protein araC binds to the DNA chain forming a DNA loop and thus prevents the binding between DNA and RNA polymerase. However, when two arabinose molecules bind to the araC dimer, the DNA loop is broken and the binding between the RNA polymerase to the promoting region is possible.
+					</p>
+					<p>
+					We have proposed a possible mathematical model describing the P<sub>BAD</sub>/araC mechanism. This <b>theoretical model</b> has been included in the <a href="#simtool" class="inner-link">Simulation Tool</a>, so that with Printeria it is possible to predict the behavior of P<sub>BAD</sub>/araC inducible promoter circuits.
+					</p>
+					<p><i>Model design</i></p>
+					<p>
+					After analyzing the biochemical functioning of the Pbad/araC promoter, we propose the following cellular scheme that would describe its behavior. It is worth mentioning that we have assumed a very high concentration of araC in the cell, as we have considered that the gene encoding for the araC protein is found in the cell genome (e.g. in DH5<meta charset="utf-8">&alpha; strains).
+					</p>
+					<img src="https://static.igem.org/mediawiki/2018/3/3c/T--Valencia_UPV--esq_pbadUPV2018.png">
+					<p>From the cell schema we can deduce the following biochemical reactions.</p>
+					<img src="https://static.igem.org/mediawiki/2018/d/d8/T--Valencia_UPV--react_pbadUPV2018.png">
+					<table style="width:100%">
+						<tr>
+						  <th><p>Parameter</p></th>
+						  <th><p>Description</p></th>
+						  <th><p>Units</p></th>
+						  <th><p>Value</p></th>
+						</tr>
+						 <tr>
+						  <td><p>D</p></td>
+						  <td><p>Diffusion coefficient</p>
+						  <td><p>min<sup>-1</sup></p></td>
+						  <td><p><a href="">---</a></p></td>
+						 </tr>
+						 <tr>
+						  <td><p>V<sub>c</sub></p></td>
+						  <td><p>Cell volume(<meta charset="utf-8">&mu;L)/External volume(<meta charset="utf-8">&mu;L)</p>
+						  <td><p>adimensional</sup></p></td>
+						  <td><p><a href="">---</a></p></td>
+						 </tr>
+						 <tr>
+						  <td><p>k<sub>u</sub></p></td>
+						  <td><p>Dissociation rate of (P<sub>BAD</sub>.araC.arabi) complex</p>
+						  <td><p>min<sup>-1</sup></p></td>
+						  <td><p><a href="">---</a></p></td>
+						 </tr>
+						 <tr>
+						  <td><p>k<sub>d</sub></p></td>
+						  <td><p>Dissociation constant of (P<sub>BAD</sub>.araC.arabi) complex</p>
+						  <td><p>molec<sup>2</sup></p></td>
+						  <td><p><a href="">---</a></p></td>
+						 </tr>
+						 <tr>
+						  <td><p>k<sub>b</sub></p></td>
+						  <td><p>Association rate of (P<sub>BAD</sub>.araC.arabi) complex</p>
+						  <td><p>molec<sup>-2</sup>min<sup>-1</sup></p></td>
+						  <td><p><a href="">---</a></p></td>
+						 </tr>
+						 <tr>
+						  <td><p>K</p></td>
+						  <td><p>Transcription rate is calculated by multiplying the transcription rate per plasmid K<sub>t</sub> by the mean number of plasmids in cell c<sub>n</sub> . In our case, we are using the pMBI replication origin, so c<sub>n</sub> <meta charset="utf-8">&sim; 500, and C<sub>R</sub> = K<sub>R</sub>·c<sub>n</sub>.</p></td>
+						  <td><p>min<sup>-1</sup></p></td>
+						  <td><p><a href="">---</a></p></td>
+						 </tr>
+						 <tr>
+						  <td><p><meta charset="utf-8">&alpha;</p></td>
+						  <td><p>(P<sub>BAD</sub>.araC) basal activity constant</p>
+						  <td><p>Admimensional</p></td>
+						  <td><p><a href="">---</a></p></td>
+						 </tr>
+						 <tr>
+						  <td><p>p</p></td>
+						  <td><p>Translation rate</p></td>
+						  <td><p>min<sup>-1</sup></p></td>
+						  <td><p><a href="">---</a></p></td>
+						 </tr>
+						 <tr>
+						  <td><p>d<sub>m</sub></p></td>
+						  <td><p>mRNA degradation rate</p></td>
+						  <td><p>min<sup>-1</sup></p></td>
+						  <td><p><a href="">---</a></p></td>
+						 </tr>
+						 <tr>
+						  <td><p>d<sub>p</sub></p></td>
+						  <td><p>PoI degradation rate</p></td>
+						  <td><p>min<sup>-1</sup></p></td>
+						  <td><p><a href="">---</a></p></td>
+						 </tr>
+						 <tr>
+						  <td><p><meta charset="utf-8">&mu;</p></td>
+						  <td><p>Dilution rate</p></td>
+						  <td><p>min<sup>-1</sup></p></td>
+						  <td><p><a href="">---</a></p></td>
+						 </tr>
+						 <tr>
+						  <td><p>K<sub>max</sub></p></td>
+						  <td><p>Maximum growth capacity</p></td>
+						  <td><p>cells</p></td>
+						  <td><p><a href="">---</a></p></td>
+						 </tr>
+					</table>
+                    <p>
+                    With the reactions and approximations already established, we apply the <a href="">Law of Mass Action kynetics</a> (LMA) to deduce the differential equations. The LMA establishes that the variation of the species resulting from a reaction is proportional to the product of the reactants. In our case, to the set of differential equations deduced from the reactions we add a cell growth equation. Biochemical species and our model of equations obtained are shown below.
+                    </p>
+                    <table style="width:100%">
+						<tr>
+						  <th><p>Variable</p></th>
+						  <th><p>Biochemical species</p></th>
+						  <th><p>Units</p></th>
+						</tr>
+						 <tr>
+						  <td><p>n<sub>1</sub></p></td>
+						  <td><p>mRNA</p></td>
+						  <td><p>Molecules</p></td>
+						 </tr>
+						 <tr>
+						  <td><p>n<sub>2</sub></p></td>
+						  <td><p>PoI</p></td>
+						  <td><p>Molecules</p></td>
+						 </tr>
+						 <tr>
+						  <td><p>n<sub>3</sub></p></td>
+						  <td><p>Number of cells</p></td>
+						  <td><p>Cells</p></td>
+						 </tr>
+					</table>
+                    <img src="https://static.igem.org/mediawiki/2018/e/e1/T--Valencia_UPV--eq_constUPV2018.png">
+					<a class="anchorOffset" id="LuxR_LuxI"></a>
+                    <h5>LuxR LuxI inducible model</h5>
+                  <a class="anchorOffset" id="charact"></a>
+                  <h3>Characterization procedures</h3>
+                    <p>
+					This year, Valencia UPV iGEM team has designed an extensive <a href="https://2018.igem.org/Team:Valencia_UPV/Part_Collection" target="_blank">Part Collection</a> in purpose of allowing the user to design multiple genetic constructions and experiments. One of our main objectives has been to <b>show the user clear and structured information</b> about the pieces that make up the Printeria kit. For this reason, we have considered the characterization of the parts as a priority when developing the project. In this way, we have elaborated some <b>procedures</b> which have allowed us to systematically obtain and structure information from the parts.
+					<ul>
+						<li><p>
+							<a href="#spectra" class="inner-link">Procedure for obtaining protein spectra</a>
+						</p></li>
+						<li><p>
+							<a href="#comparison" class="inner-link">Comparison between sfGFP and GFP relative fluorescence intensity</a>
+						</p></li>
+					</ul>
+					</p>
+					<a id="spectra" class="anchorOffset"></a>
+					<h4>Procedure for obtaining protein spectra </h4>
+					<p>
+					Obtaining <b>excitation and emission spectra</b> is a fundamental aspect in the process of characterization of a fluorescent protein. <b>Each protein has a characteristic spectrum</b>, which indicates the energy in which the molecule is excited or emits at a certain wavelength.
+					</p>
+					<p>
+					The characterization of the reporter proteins by excitation and emission spectra are of great importance in experimentation. When contrasting experimental information with the theoretical results of mathematical models, we experiment with numerous reporter proteins, such as fluorescence proteins or chromoproteins. However, the <b>fluorescence data obtained must be corrected applying diferent operations</b> in order to obtain representative fluorescence data:
+					</p>
+					<ol style="list-style-type: upper-roman;">
+						<li><p>
+							The subtraction of the medium fluorescence:
+							<img src="https://static.igem.org/mediawiki/2018/9/97/T--Valencia_UPV--fsubsUPV2018.png">
+						</li></p>
+						<li><p>
+							The quotient of the fluorescence with the gain:
+							<img src="https://static.igem.org/mediawiki/2018/b/b0/T--Valencia_UPV--fgainUPV2018.png">
+						</li></p>
+						<li><p>
+							The quotient of the fluorescence with excitation and emission efficiency:
+							<img src="https://static.igem.org/mediawiki/2018/2/21/T--Valencia_UPV--fcorrUPV2018.png">
+						</li></p>
+					</ol>
+					<a id="esquema_ops" class="anchorOffset"></a>
+					<img src="https://static.igem.org/mediawiki/2018/6/6e/T--Valencia_UPV--esquema_opsUPV2018.png">
+					<p>
+					Until now, the corrections applied to the experiments performed with these reporter proteins were the subtraction of the fluorescence of the medium, and the division by the gain factor of the measuring equipment. However, <b>with the protein spectra we can also normalize the fluorescence data to values that would have been obtained with maximum excitation and emission</b>.
+					</p>
+					<p>
+					Owing to this reason, <b>a protocol has been established</b> in the lab by Lab and Modeling team <b>to obtain the spectrum of any reporter protein</b>.
+					</p>
+					<p><i>Materials:</i></p>
+					<ul>
+						<li>
+							<p>Measuring equipment: <b>Biotek Cytation3</b></p>
+						</li>
+						<li>
+							<p><b>96 well plate</b></p>
+						</li>
+						<li>
+							<p><b>MATLAB 2018a software</b></p>
+						</li>
+					</ul>
+					<p><i>Procedure:</i></p>
+					<ol>
+						<li><p>
+							We look for the <a href="https://www.thermofisher.com/es/es/home/life-science/cell-analysis/labeling-chemistry/fluorescence-spectraviewer.html" target="_blank" style="padding-right: 0">theoretical spectra</a> of the protein to be measured or similar molecules in order to determine the wavelength at which the protein is excited or emitted at maximum energy, i.e. where the theoretical spectral peaks occurs.
+						</p></li>
+						<li><p>
+							We define the protocol of our equipment to get the absorbance and fluorescence dataset. In our protocol, the most important parameters to be established are summarized in the following Table.
+							<table style="width:100%">
+							  <tr>
+							    <th><p>Parameters</p></th>
+							    <th><p>Description</p></th>
+							  </tr>
+							  <tr>
+							    <td><p>Number of samples</p></td>
+							    <td><p>From <b>3 to 6 samples of reporter protein</b> and <b>3 samples of medium (LB or M9 medium)</b></p></td>
+							  </tr>
+							  <tr>
+							    <td><p>Temperature</p></td>
+							    <td><p>37 ºC</p></td>
+							  </tr>
+  							  <tr>
+							    <td><p>Shake</p></td>
+							    <td><p>Double Orbital. 01:00 (MM:SS)</p></td>
+							  </tr>
+  							  <tr>
+							    <td><p>Absorbance. Optical Density (OD)</p></td>
+							    <td><p>Wavelenght at <b>600 nm</b> emission</p></td>
+							  </tr>
+							  <tr>
+							    <td><p>Excitacion and emission scans</p></td>
+							    <td><p>The scans occur between <b>two wavelength limit values</b>. The established range will depend on the theoretical spectrum of the protein.</p></td>
+							  </tr>
+							  <tr>
+							    <td><p>Excitacion and emission wavelengths</p></td>
+							    <td><p>These values will depend on the range, and therefore on the spectrum. Values far from the theoretical peak lead to more attenuated fluorescence curves, and values very close to the peak can lead to overlap and error in reading the data. Therefore, <b>a compromise must be reached between curve resolution and reading overlap</b>.</p></td>
+							  </tr>
+							  <tr>
+							    <td><p>Gain (G)</p></td>
+							    <td><p>Normally the gain value is <b>G = 60</b>, although for proteins with lower fluorescence, it is recommended that G takes higher values.</p></td>
+							  </tr>
+							</table>
+						</p></li>
+					<li><p>
+						The experiment is introduced, and the experimental absorbance data and the fluorescence curves of the samples are obtained with the medium fluorescence correction applied.
+					</p></li>
+					<li><p>
+						The dataset is exported to an Excel file.
+					</p></li>
+					<li><p>
+						The MATLAB script <a href="https://static.igem.org/mediawiki/2018/f/f7/T--Valencia_UPV--spectrummUPV2018.txt" target="_blank">spectrum.m</a> for fluorescent proteins is executed:
+						<ol style="list-style-type: lower-latin;">
+							<li>
+								<p>Dataset is extracted from the Excel file. We discard readings that have suffered overlap, or that take negative values.</p>
+							</li>
+							<li>
+								<p>The fluorescence curves of all samples are averaged, and the result is normalised (from 0 to 100%).</p>
+							</li>
+							<li>
+								<p>Graphs of the normalized absorption and emission spectra are plotted. The X-axis represents the wavelength (nm), and the Y-axis represents the normalized fluorescence intensity (%).</p>
+							</li>
+						</ol>
+					</p></li>
+					</ol>
+					<img src="https://static.igem.org/mediawiki/2018/e/e8/T--Valencia_UPV--sfGFP_spectrumUPV2018.png">
+					<p>
+					In practical terms, the protocol has been applied to all of our reporter proteins: <a href="http://parts.igem.org/Part:BBa_K2656022" target="_blank">GFP</a>, <a href="http://parts.igem.org/Part:BBa_K2656013" target="_blank">sfGFP</a>, <a href="http://parts.igem.org/Part:BBa_K2656021" target="_blank">YFP</a> and <a href="http://parts.igem.org/Part:BBa_K2656014" target="_blank">mRFP</a>. All results can be found in our parts collection as well as in the <a href="http://parts.igem.org/Catalog" target="_blank">iGEM catalog</a>.
+					</p>
+					<p>
+					In the particular case of a reporter chromoprotein, such as <a href="http://parts.igem.org/Part:BBa_K2656018" target="_blank">amilCP</a>, we do not measure fluorescence, but absorbance. In this case, in rder to obtain the corrected absorbance curve, we must subtract from the cell absorbance data with the reporter protein the absorbance of a medium with cells without chromoprotein. Once the data have been corrected, we normalize them between values of 0 and 100 and with this we elaborate the graph. The protocol used can be found in the MATLAB script <a href="https://static.igem.org/mediawiki/2018/8/8b/T--Valencia_UPV--amilCP_spectrumUPV2018.txt" target="_blank">amilCP_spectrum.m</a>.
+					</p>
+					<p>
+					Finally, we have also obtained the spectra of the fluorescein molecule. These spectra have been used to correct the fluorescence data used in the <a href="" target="_blank">Interlab Study</a> to obtain the Relative Fluorescence Units (RFU) to Molecules of Equivalent Fluorochrome (MEFL) conversion factor. In addition, the <a href="#comparison" class="inner-link">comparison</a> between GFP and sfGFP proteins RFU has been possible thanks to fluorescein spectra.
+					</p>
+					<table style="width:100%">
+						  <tr>
+						    <th><p>Parameters</p></th>
+						    <th><p>Value</p></th>
+						  </tr>
+						  <tr>
+						    <td><p>Number of samples</p></td>
+						    <td><p>6 samples. 3 samples of medium (LB or M9 medium)</p></td>
+						  </tr>
+						  <tr>
+						    <td><p>Temperature</p></td>
+						    <td><p>37 ºC</p></td>
+						  </tr>
+							  <tr>
+						    <td><p>Shake</p></td>
+						    <td><p>Double Orbital. 01:00 (MM:SS)</p></td>
+						  </tr>
+							  <tr>
+						    <td><p>Absorbance. Optical Density (OD)</p></td>
+						    <td><p>Wavelenght at <b>600 nm</b> emission</p></td>
+						  </tr>
+						  <tr>
+						    <td><p>Emission range</p></td>
+						    <td><p>[495 - 580] nm</p></td>
+						  </tr>
+						  <tr>
+						    <td><p>Fixed excitation wavelength</p></td>
+						    <td><p>480 nm</p></td>
+						  </tr>
+						  <tr>
+						    <td><p>Excitation range</p></td>
+						    <td><p>[430 - 520] nm</p></td>
+						  </tr>
+						  <tr>
+						    <td><p>Fixed emission wavelength</p></td>
+						    <td><p>545 nm</p></td>
+						  </tr>
+						  <tr>
+						    <td><p>Gain (G)</p></td>
+						    <td><p>60</p></td>
+						  </tr>
+					</table>
+					<img src="https://static.igem.org/mediawiki/2018/1/11/T--Valencia_UPV--fluorescein_spectrumUPV2018.png">
+					<a id="comparison" class="anchorOffset"></a>
+					<h4>Comparison between sfGFP and GFP relative fluorescence intensity</h4>
+					<p>
+					<b>One of the main problems</b> we encounter when processing the results of experiments <b>in Synthetic Biology are the units of measurement of fluorescence data</b>. Unlike absorbance, where there is a simple conversion between Optical Density (OD) and cell number, <b>there is no direct relationship between the Relative Units of Fluorescence (RFU) and the number of protein molecules in the cell</b>. Moreover, RFUs vary among reporter fluoroproteins: for example, an RFU of the GFP protein does not have to be equivalent to an RFU of the mRFP protein, sfGFP, etc.
+					</p>
+					<p>
+					Thanks to initiatives such as the <a href="" target="_blank">Interlab Study</a>, we have been able to go a step further and obtain a <b>MEFL/cell<sub>GFP</sub> factor</b> of equivalence between the RFU of the <a href="http://parts.igem.org/Part:BBa_K2656022" target="_blank">GFP</a> protein and the Molecules of Equivalent Fluorochrome (MEFL). This relationship is an important breakthrough, as it can give us a more <b>accurate estimation</b> of the amount of GFP molecules in the cell.
+					</p>
+					<p>
+					However, <b>another reporter protein widely used in the experiments is the <a href="http://parts.igem.org/Part:BBa_K2656013" target="_blank">sfGFP</a></b>. This protein has a much faster folding than GFP, which translates into a <b>higher fluorescence intensity per molecule</b>. In order to obtain the MEFL/cell<sub>sfGFP</sub> factor from sfGFP, the Printeria Modeling and Lab teams have designed a <b>comparative experiment between both proteins</b>. The experiment consists of designing <b>two identical transcriptional units (TU), changing only the CDS sequence</b> so that each TU will produce GFP and sfGFP, respectively. It should be added that this experiment is based on <b>two fundamental assumptions</b>:
+					</p>
+					<p>
+					<ol style="list-style-type: upper-latin;">
+						<li><p>
+							The number of GFP molecules produced in the cells is equivalent to the number of MEFL.
+						</p></li>
+						<li id="cond_equiv_molec"><p>
+							Given two TU with identical promoters, RBS and terminators, but with different CDS, under the same experimental conditions, the number of molecules produced by each TU is the same.
+						</p></li>
+					</ol>
+					</p>
+					<p>
+					Taking these axioms into account, the materials and procedure followed to calculate the MEFL/cell<sub>sfGFP</sub> factor were as follows.
+					</p>
+					<p><i>Materials:</i></p>
+					<ul>
+						<li>
+							<p>Measuring equipment: <b>Biotek Cytation3</b></p>
+						</li>
+						<li>
+							<p><b>96 well plate</b></p>
+						</li>
+						<li>
+							<p><b>MATLAB 2018a software</b></p>
+						</li>
+					</ul>
+					<p><i>Procedure:</i></p>
+					<ol>
+						<li><p>
+							From the fluorescein spectrum, and the fluorescence data obtained from the <a href="" target="_blank">Interlab Study experiment</a>, we apply medium, gain and efficiency <a href="#esquema_ops" class="inner-link">corrections </a>.
+						</p></li>
+						<li><p>
+							The data are introduced in the <a href="" target="_blank">Excel file</a> of the Interlab Study. From the <i>Fluorescein standard curve</i> sheet we can obtain the MEFL/RFU factor, and then calculate the MEFL/cell factor.
+							<img src="https://static.igem.org/mediawiki/2018/a/af/T--Valencia_UPV--calc_MEFLUPV2018.png">
+							<table style="width:100%">
+							    <tr>
+							    	<th><p>Parameters</p></th>
+							    	<th><p>Value</p></th>
+							  	</tr>
+							  	<tr>
+							    	<td><p>Cuture volume (96 well plate)</p></td>
+							    	<td><p>200 <meta charset="utf-8">&mu;L</p></td>
+							  	</tr>
+							  	<tr>
+							    	<td><p>Cells per <meta charset="utf-8">&mu;L per OD unit</p></td>
+							    	<td><p>200 cells/OD <meta charset="utf-8">&mu; L</p></td>
+							  	</tr>
+							  	<tr>
+							    	<td><p>MEFL/RFU factor</p></td>
+							    	<td><p>4,38.10<sup>10</sup> MEFL/RFU</p></td>
+							  	</tr>
+							  	<tr>
+							    	<td><p><b>MEFL/cell<sub>GFP</sub> factor</b></p></td>
+							    	<td><p><b>273.7 MEFL OD/RFU cell</b></p></td>
+							  	</tr>
+							</table>
+						</p></li>
+						<li><p>
+							We establish the protocol of the experiment and the parameters of our equipment. The experiment consists in the measurement of the absorbance and fluorescence of <b>two TU with identical promoter, RBS and transcriptional terminator</b>, but whose CDS codifies for <a href="http://parts.igem.org/Part:BBa_K2656105" target="_blank">GFP</a> and <a href="http://parts.igem.org/Part:BBa_K2656101" target="_blank">sfGFP</a> proteins. The most relevant information of the experiment is detailed in the following Table:
+							<table style="width:100%">
+							  <tr>
+							    <th><p>Parameters</p></th>
+							    <th><p>Description</p></th>
+							  </tr>
+							  <tr>
+							    <td><p>Time</p></td>
+							    <td><p>06:00:00 (HH:MM:SS). Measurement interval: 05:00 (MM:SS)</p></td>
+							  </tr>
+							  <tr>
+							    <td><p>Number of samples</p></td>
+							    <td><p>8 samples for each TU. Total samples: 16</p></td>
+							  </tr>
+							  <tr>
+							    <td><p>Temperature</p></td>
+							    <td><p>37 ºC</p></td>
+							  </tr>
+  							  <tr>
+							    <td><p>Shake</p></td>
+							    <td><p>Double Orbital Continuous</p></td>
+							  </tr>
+  							  <tr>
+							    <td><p>Absorbance. Optical Density (OD)</p></td>
+							    <td><p>Wavelenght at <b>600 nm</b> emission</p></td>
+							  </tr>
+							  <tr>
+							    <td><p>Excitacion wavelength</p></td>
+							    <td><p>485 nm</p></td>
+							  </tr>
+							  <tr>
+							    <td><p>Emission wavelength</p></td>
+							    <td><p>528 nm</p></td>
+							  </tr>
+							  <tr>
+							    <td><p>Gain</p></td>
+							    <td><p>60</p></td>
+							  </tr>
+							</table>
+						</p></li>
+						<li><p>
+							The experiment is introduced into the equipment and the absorbance and fluorescence curves are obtained, as well as the curve of <b>Fluorescence F/Absorbance OD</b> or <b>FOD ratio</b> with the applied medium fluorescence correction.
+						</p></li>
+						<li><p>
+							The absorbance and fluorescence data are exported to an Excel file.
+						</p></li>
+						<li><p>
+							Running the MATLAB <a href="https://static.igem.org/mediawiki/2018/c/c9/T--Valencia_UPV--GPF_sfGFP_comparisonUPV2018.txt" target="_blank">GFP_sfGFP_comparison.m</a> script:
+							<ol style="list-style-type: lower-latin;">
+								<li>
+									<p>We extract the FOD data from the GFP and sfGFP from the Excel file.</p>
+								</li>
+								<li>
+									<p>We apply the gain and efficiency <a href="#esquema_ops">corrections </a> of the spectrum to the FOD.</p>
+								</li>
+								<li>
+									<p>We plot the FOD curves and look for a stationary equilibrium time interval in the expression of the GFP and sfGFP proteins. In our <a href="#comp_graph" class="inner-link">experiment</a>, we have decided to chose the interval [145,290] min</p>
+								</li>
+								<li>
+									<p>We obtain the average value of both regions.</p>
+								</li>
+								<li>
+									<p>We calculate the number of GFP MEFLs by multiplying the FOD data by the MEFL/cell<sub>GFP</sub> factor.</p>
+								</li>
+								<li>
+									<p>If we assume that the number of molecules expressed by both TU is the same, <a href="#cond_equiv_molec" class="inner-link">the number of MEFL for the calculated GFP is the same as the number of MEFL for the sfGFP</a>.</p>
+								</li>
+								<li>
+									<p>We calculate the MEFL/cell<sub>sfGFP</sub> factor by dividing the MEFL number by the average FOD value of the sfGFP in the stationary region. In our case, MEFL/cell<sub>sfGFP</sub> = 121.4 MEFL·OD/RFU·cell.</p>
+								</li>
+							</ol>
+						</p></li>
+					</ol>
+					<a id="comp_graph"></a>
+					<img src="https://static.igem.org/mediawiki/2018/9/95/T--Valencia_UPV--comp_graphUPV2018.png">
+					<table style="width:100%">
+							  <tr>
+							    <th><p>Reporter protein</p></th>
+							    <th><p>MEFL/cell factor (MEFL·OD/RFU·cell)</p></th>
+							  </tr>
+							  <tr>
+							    <td><p>GFP</p></td>
+							    <td><p>273.7</p></td>
+							  </tr>
+							  <tr>
+							    <td><p>sfGFP</p></td>
+							    <td><p>121.4</p></td>
+							  </tr>
+							  <tr>
+							</table>
+					<p>
+					With MEFL/cell factors, <b>FOD data</b> obtained in any experiment in which GFP or sfGFP was used as reporter proteins <b>can be transformed into equivalent fluorescein molecules</b> by applying the following ratio: <b>Molecules = (MEFL/cell factor) · FOD</b>.
+					</p>
+					<p>
+					The new data give us a more accurate estimation of the number of molecules in the cell. Consequently, by relating the experimental results with the theoretical mathematical models in the <a href="#optimization" class="inner-link">optimization process</a>, the <b>parameters of the model acquire values more consistent with their physical significance</b>, working in equivalent molecules and not in arbitrary units.
+					</p>
+					<a class="anchorOffset" id="simtool"></a>
+                  <h3>Simulation Tool</h3>
+                    <p>
+                    	<!--Contenido de Simulation Tool -->
+                    </p>
+               </div>
+            </div>
+         </div>
+      </div>
+   </div>
+</section>
+        <!--<div class="loader"></div>-->
+        <a class="back-to-top inner-link" href="#start" data-scroll-class="100vh:active" style="font-size: 82%;">
+            <i class="stack-interface stack-up-open-big"></i>
+        </a>
+    </body>
+        <script type="text/javascript"  src="https://2018.igem.org/Template:Valencia_UPV/scriptsJS?action=raw&ctype=text/javascript"></script>
+        <script type="text/javascript"  src="https://2018.igem.org/Template:Valencia_UPV/jquery311minJS?action=raw&ctype=text/javascript"></script>
+        <script type="text/javascript"  src="https://2018.igem.org/Template:Valencia_UPV/bootstrapMinJS?action=raw&ctype=text/javascript"></script>
+        <script type="text/javascript"  src="https://2018.igem.org/Template:Valencia_UPV/parallaxJS?action=raw&ctype=text/javascript"></script>
+       <script type="text/javascript"  src="https://2018.igem.org/Template:Valencia_UPV/smoothscrollminJS?action=raw&ctype=text/javascript"></script>
+       <script type="text/javascript"  src="https://2018.igem.org/Template:Valencia_UPV/jquerystepsminJS?action=raw&ctype=text/javascript"></script>
 </html>
+{{Valencia_UPV/footer}}

Difference between revisions of "Team:Valencia UPV/Experiments"

Revision as of 14:59, 6 October 2018

Index

Modeling

Models & Experiments

Modeling process

Experimental protocol

Multi-objective parameter optimization

Constituve expression models

Inducible expression models

PBAD/araC inducible model

LuxR LuxI inducible model

Characterization procedures

Procedure for obtaining protein spectra

Comparison between sfGFP and GFP relative fluorescence intensity

Simulation Tool

P_BAD/araC inducible model