Team:Valencia/Modeling

From 2012.igem.org



Modeling


Light, sucrose and AHL yield in Synechococcus elongatus:


In first place we adapted a metabolic model of Synechococcus with the functions we included with our constructs. The model is an algorithm applied to a metabolic network of reactions, originally designed for Synechocystis sp. by Montagud et al. 2010 and 2011. and recently adapted for Synechococcus sp., where optimization of determined fluxes with constraints on certain reactions which make it behave like a living cell.
We transformed a wildtype S. elongatus with the luciferase genes (luxAB) regulated by the promoter psbA. The reactions of bioluminescence expressed by the LuxAB construct were introduced into the network as follows:

  • LuxE8 : octanoic acid + ATP + NADPH -> octanal + AMP + NADP+ + diphosphate

  • LuxAB8 : octanal + FMNH2 + oxygen -> FMN + H2O + octanoic acid + light

  • LuxE14 : tetradecanoic acid + ATP -> tetradecanal + AMP + NADP+ + diphosphate

  • LuxAB14 : tetradecanal + FMNH2 + oxygen -> FMN + H2O + tetradecanoic acid + light

  • LuxE16 : hexadecanoic acid + ATP -> hexadecanal + AMP + NADP+ + diphosphate

  • LuxAB16 : hexadecanal + FMNH2 + O2 -> FMN + H2O + hexadecanoic acid + luz

  • FMN_red : FMNH2 + NADP+ <-> FMN + NADPH


May be noted that the luciferin (FMNH2) is reduced (regenerated) by the reaction FMN_red, which is not naturally present in Synechococcus, but which we added to our model to ensure light production in a persistent way. As a future aim, we plan to include luciferin regeneration gene cassettes to this wildtype. If we didn’t insert this adjustment, the model would calculate flow 0, as luminescence would extinguish in a few seconds with the scarce FMNH2, so the model has no solutions for infinite time flow. This helped us to notice (and then assured from bibliographic support) this handicap of cyanobacteria as light producers. This partly triggered the reorientation of our research towards the bispecific coculture idea.

Moreover, we programmed bioluminescence functions with 3 types of fatty acid (8, 14 and 16 carbon chain length) which are most abundant in the cyanobacterial cell so could be used as ‘fuel’ with greatest probability, and lied between the operational chain-lenght limits for the luciferase.

This yielded results as we optimized the model’s algorithm to maximize light production, by adjusting the flow rates of the other reaction of the metabolism under certain restrictions which make it a biologically reasonable maximum (principle of operation for every result we obtained from the algorithm)value of 22509mmol/g DCW/h, equivalent to 6.48x108photons/cell/s.

In second place, we adapted the wildtype model with the sucrose export induced by the expression of the cscB transporter gene included in the transformed strain from Harvard. Here we modeled the maximum values of sucrose export for 2 different growth restrictions, 0.09mmol Biomass/gDCW*/h (a value near to the maximum growth, similar to an averaged exponential phase culture) and 1x10-6mmol Biomass/gDCW/h (a negligible growth value which will keep all the vital reactions working but not divert a significant amount of fixed carbon to increasing biomass, as in an averaged dynamic equilibrium of a stationary phase culture).

Maximum export values yielded:
With minimal growth constraint: 0.3075mmol/g DCW/h
With maximal growth constraint: 0.027mmol/g DCW/h
As we planned to transform this cscB strain to express luxI, the protein that synthesizes AHL (autoinducer molecule), we introduced the reaction in the model as follows:

AHL: Hexanoyl-(acyl carrier protein) + S-adenosyl-L-methionine -> AHL + S-methyl-5'-thioadenosine + acyl-carrier protein

Our biomachine has 2 synthetic functions inserted, for which we have an interest of flow maximization. Unfortunately, export of sucrose and export of AHL, both carbon based molecules, sets a biochemical competition for carbon redirection. Therefore we came out with 3 modeled outputs: Parabolic curve of maximization of AHL+Sucrose (fig a), and the linear functions of AHL vs. sucrose export in both growth scenarios (fig b, c). All units are expressed in mmol/g DCW/h.


Fig a. Parabolic curve representing the product of AHL·Sucrose export flows
versus net sucrose export constraint, at the maximum growth constrain.
The maximum value, of near 0.016mmol/g DCW/h, corresponds
to the 0.013mmol/g DCW/h sucrose export constraint.



Fig b. Linear equation representing AHL export versus sucrose export
constraint, at the minimum growth constraint.



Fig c. Linear equation representing AHL export versus sucrose export constraint,
at the maximum growth constraint. In this scenario, the value of AHL export
for a sucrose export of 0.049 mmol/g DCW/h (units adapted from Ducat et
al, 2012, at the experimental conditions of our setting) is 0.221mmol/g
DCW/h, which is a very high value, near to the the maximum
(0.2633mmol/g DCW/h).



We just require 1 or 2 molecules per cell of Aliivibrio fischeri to induce biolumminescence (Kaplan & Greenberg, 1985). In Fig c we can see the how small sucrose export values, such as our real fluxes affect very little the export of AHL. This was tested on a maximal growth restriction, as the forward idea is to keep the bioreactor in a continuous exponential phase, which corresponds to a flow of 0.00974 mmol/g DCW/h . This imposes a constant flux of sucrose and AHL, which is optimum for our system. The answer of the model to this constant export of sucrose is a flow of 0.015mmol AHL/g DCW/h, which is a pretty high value. As you will see later, this value of AHL export is fundamental for the development of our next model.

-*: DCW=Dry Cell Weight, where an average Synechococcus cell has a dry weight of 3.87ng (Rosales et al. 2005).


Synergic model:


We achieved our main goal of integrating and connecting all the system into a single model capable of predicting the luminescence of the consortium system in different bioreactor setting scenarios, assuming fast diffusion mechanisms of substances in the medium with a pump system and a pump-filter-return control of population density in both S. elongatus and A. fischeri modules.

The main part is based on the model of bioluminescence regulation developed by Belta et al. 2001 (source of the following figures), a hybrid model of 9 differential equations which predicts the behavior of the whole quorum sensing regulatory mechanism of luciferase expression (including cAMP, AHL, LuxR and LuxI expression) in A. fischeri.










We modified this model to restrict the growth of the A. fischeri population into a constrained volume (biolamp compartment) meanwhile we adopted a greater volume for the dilution of the autoinducer (AHL), as the population cannot cross the membrane to colonize the whole of the medium.
Such 'dilution volume' represents the annexing of the photobioreactor of S. elongatus and the extra volume from tubing and pumping systems where the medium flows (50ml).

We gave the system a new input of AHL derived from the S. elongatus population, which has a Ks (constant rate of AHL production/cell resultant from our modified Synechococcus metabolic network model) of 2.38x10-7 nmol AHL /cell/s, multiplying the number of cells (culture density x culture volume – of the photobioreactor compartment). The model still counts with the AHL produced by the A. fischeri population, and the degradation half-life in the cells. As we said above, the diffusion speed of AHL is considered instantaneous (as the rate’s scale is very high compared with other variables, such as cell growth or gene expression).

The application of this model is focused in the optimization of the design of the future bioreactor which may contain the system. We are looking for the values of compartment volumes, relative volume, population densities and total volume which can assure a diel control of luminescence (in 12h cycles, considering export rates and the AHL degradation half-life of 10h) by the light induced activity of our S. elongatus biomachine; and, within those values, the settings with maximum luminescence values (values higher than 1000nM luciferase in the A. fischeri module volume).

A second filter is be applied to the group of solutions obtained from these premises, with the result of a further model on export and consumption rates of sucrose. This filter discards the settings where insufficient sucrose is exported to the culture to sustain the growth of the A. fischeri population. We managed to adapt some experimental values of sucrose export and consumption from the literature:

The export of sucrose from our cyanobacteria (Ke) is: 2.07x10-11 nmol/cell/s (units adapted from Ducat et al, 2012)
The consumption of sucrose by A. fischeri (Kc) is: 0.00000150 nmol/cell/s (units adapted from oxygen consumption rates at maximum bioluminescent activity, at normal glycolytic route assumptions in aerobic conditions, from Makemson 1985).

dSucrose=Ke·(S. elongatus cell density).(photobioreactor module volume) - Kc·(A. fischeri cell density)·(biolammp module volume) ; where only dSucrose (> or =) 0nM/s are accepted.

This amendment model sets the total volume of dilution, the compartmental volumes of occupation (in a range between 0.001-10l) and the cell densities for S. elongatus cscB and for A. fischeri (in a range between 0 and 109 cell/ml) as controlled variables.
The output gives the multivariate scenarios where the rate of change of sucrose concentration in the common broth is 0 or positive, so that the energy budget is not a shortfall.

The results of our model were non-unimodal. This means that we cannot rationally present in short the outcome of multivariate solutions we obtained as min/max limits for each variable. Our high resolution discrete outcome is an enormous multidimensional permutation. If we lower the resolution for a presentable report, the interpolation error is too high. We look forward to rationalize this by the means of analytic resolution of the differential equations system into a simpler function applicable to an easy assessment of the design parameters for the bioreactor.


References


Belta, C., Schug, J., Dang, T., Kumar, V., Pappas, G.J., Rubin, H., Dunlap, P. 2001. Stability and reachability analysis of a hybrid model of luminescence in the marine bacterium Vibrio fischeri. Decis. In Cont., 1, 869-874.

Ducat, D. C., Avelar-Rivas, A. J., Way, J. C. & Silver, P. A. (2012) Rerouting carbon flux to enhance photosynthetic productivity. Applied and Environmental Microbiology, 78(8):2660–2668.

Kaplan H., Greenberg, E.B. 2004. Diffusion of autoinducer is involved in regulation of the Vibrio fischeri luminescence system. Journ. Of Bacteriol. 1210-1214.

Makemson J., 1985. Luciferase-dependent oxygen consumption by bioluminescent Vibrios. Amer. Soc. for Microb.165(2): 465-466.

Montagud, A., Navarro E., Fernández de Córdoba, P., Urchueguía, J.F., Patil, K.R. (2010). Reconstruction and analysis of genome-scale metabolic model of a photosynthetic bacterium. BMC Syst. Biol. 4:156.

Montagud, A., Zelezniak, A., Navarro E., Fernández de Córdoba, P., Urchueguía, J.F., Patil, K.R. (2011). Flux coupling and transcriptional regulation within the metabolic network of the photosynthetic bacterium Synechocystis sp. PCC6803. Biotechnol. J. 6(3):330-42.

Rosales, N., Ortega, J., Mora, R., Morales, E. 2005. Influence of salinity on the growth and biochemical composition of the cyanobacterium Synechococcus sp. Ciencias Marinas, 31(2): 349–355.