Team:Calgary/Project/OSCAR/FluxAnalysis

From 2012.igem.org

(Difference between revisions)
 
(78 intermediate revisions not shown)
Line 4: Line 4:
CONTENT=
CONTENT=
 +
[[Image:UCalgary2012_OSCAR_Flux_Analysis_Low-Res.png|250px|right|]]
<html>
<html>
-
<img src="https://static.igem.org/mediawiki/2012/6/61/UCalgary2012_OSCAR_Flux_Analysis_Low-Res.png" style="float: right; padding: 10px;"></img>
+
 
 +
<p>Flux-variability analysis (FVA) was applied to optimize the bioreactor system of OSCAR for the newly incorporated metabolic pathways
 +
(<a href="https://2012.igem.org/Team:Calgary/Project/OSCAR/Decarboxylation">Decarboxylation</a>,
 +
<a href="https://2012.igem.org/Team:Calgary/Project/OSCAR/CatecholDegradation">Decatecholization</a>,
 +
<a href="https://2012.igem.org/Team:Calgary/Project/OSCAR/Desulfurization">Desulfurization</a>, and
 +
<a href="https://2012.igem.org/Team:Calgary/Project/OSCAR/Denitrogenation">Denitrogenation</a>). FVA combines the framework of metabolic pathways with experimental enzymatic data to provide a computational platform for predicting what changes in metabolite levels will result in increased end-product. We developed a model using <a href="http://www.mathworks.com/products/matlab/">MATLAB</a> computer language for predicting what metabolites could be added to your growth media to increase production of hydrocarbons in the <i>E. coli</i> chassis. We also created a graphical user interface called the <i>OSCAR Optimizer</i> to make the the FVA program user-friendly and to allow current and future iGEM teams/scientists to input their own synthetic pathways into the FVA for analysis.  Finally, we validated this model in the wet-lab by optimizing the PetroBrick (<a href="http://partsregistry.org/wiki/index.php?title=Part:BBa_K590025">BBa_K590025</a>) system to increase hydrocarbon production, thereby saving time, and resources.</p>
 +
 
 +
<h3><center><b>Click <a href="https://static.igem.org/mediawiki/2012/5/50/UCalgary2012_OSCAR_optimizer_v1.zip"> here </a> to download the MATLAB package to test our model!</h3>
 +
<center><h3>Click <a href="https://static.igem.org/mediawiki/2012/d/d5/UCalgary2012_OSCAROptimizerManual.pdf"> here </a> to view our OSCAR Optimizer Manual.</b></center></h3>
 +
 
<h2>Background</h2>
<h2>Background</h2>
-
<p><b>What is Metabolic Flux Analysis?</b></p>
+
<h3>What is Flux Analysis?</h3>
-
<p>Metabolic Flux balance analysis (FBA) is an application of linear programming to metabolic network that uses the stoichiometric coefficients for each reaction in the system as the set of constraints for the optimization. Simply, it is a mathematical method for analyzing metabolism. This analysis requires the Steady State Assumption, which is, in chemistry, a steady state is a situation in which all state variables are constant in spite of ongoing processes that strive to change them. In addition, Flux Variability Analysis (FVA) is an extension of Flux balance analysis. FVA determines the ranges of fluxes that correspond to an optimal solution determined through FBA. Hence, FBA is able to calculate the full range of numerical values for each reaction flux within the network. </p>
+
<p>Flux Balance Analysis (FBA) uses linear programming of metabolic networks to convert each metabolite into a mathematical coefficient.  These coefficients can be connected to each other by changes associated with each enzymatic step of the pathway. FBA can then apply a mathematical method to examine how metabolites relate to each other in the network and allows the user to make generalized predictions for the growth of the organism, metabolite levels, and product output inside the cell. Flux Variability Analysis (FVA) is an extension of FBA that determines the range of reactions that result in an optimal flux through the metabolic network.  Unlike FBA, this allows the user to not only determine the most optimized network for increasing cell growth, but also other parameters like product output.</p>
-
<p><b>What are the constraints in model?</b></p>
+
<p><b>What are the constraints in the FVA model?</b></p>
-
<p>Networks can be encoded as stoichiometric matrices (S), in which each row represents a unique metabolite and each column represents a biochemical reaction. The entries in each column of this matrix are the stoichiometric coefficients of the metabolites in the reaction. Metabolites are consumed have a negative coefficient and metabolites that are produced have a positive coefficient. Since most reactions involve only a few metabolites, S is a sparse matrix. The size of S is m*n for a network with m metabolites and n reactions. </p>
+
<p>As illustrated in <b>Figure 1</b>, metabolic networks can be encoded as stoichiometric matrices, in which each row represents a specific metabolite and each column represents a biochemical reaction. The entries in each column represent the stoichiometric coefficients of the metabolites in the reaction. Metabolites that are consumed or produced have a negative or positive, respectively.</p>
 +
 
 +
</html>
 +
[[Image:Calgary FluxExample.jpg|700px|thumb|center|Figure 1: Flux Balance Analysis.  FBA uses metabolic network (left) and simulates the connections in the network as a linear algebra matrix (right). Each metabolite is listed vertically in the table and each reaction is listed horizontally. Based on the metabolites involved in each reaction, the state of the system is affected by each metabolite change.]]
 +
<html>
-
<p><b>Why use Flux Variability Analysis?</b></p>
+
<p><b>Why use FVA as opposed to FBA?</b></p>
-
<p>Biological systems often contain redundancies that contribute to their robustness. However, flux balance analysis only returns a single flux distribution that corresponds to maximal growth under given growth conditions regardless alternate optimal solutions may exist. FVA is capable to exanimate these redundancies by calculating the full range of numerical values for each reaction flux in a network. Consequently, FVA can be employed to study the entire range of achievable cellular functions as well as the redundancy in optimal phenotypes.</p>
+
<p>Biological systems often contain redundancies that contribute to their robustness. However, FBA only returns a single flux distribution that corresponds to maximal bacterial growth under the given growth conditions, even if other optimal parameters may exist. FVA is capable of examining these redundancies by calculating the full range of numerical values for each reaction flux in a network. Thus, FVA can be employed to study the entire range of achievable cellular functions, as well as the redundancy in optimal phenotypes.  FVA can also examine different ranges of bacterial growth vs. product output which is valuable in assessing validity of models in the wet-lab.</p>
<br>
<br>
 +
<h2>Using FVA to optimize OSCAR</h2>
 +
<h3>What Are We Trying To Model?</h3>
 +
<p>FVA provides a method to modulate the inputs of endogenous metabolic pathways, but could also investigate what chemicals are needed in the growth media to upregulate an introduced pathway in <i>E. coli</i>.  In the latter case, we need to specifically model the flux rate of metabolic pathways responding to different growth media conditions and generate an optimal set of metabolites that should be added to improve production rate. Development of a tool to model the optimal flux would also benefit numerous iGEM teams and research labs who engineer bacteria with novel pathways. </p>
 +
<h3>How could systems like OSCAR benefit from this model?</h3>
 +
<p>Much like chemical reactions which need optimal environmental conditions to achieve maximum production rate, microbes require certain conditions to accomplish cellular tasks at maximum speed. During industrial scale up, the optimal conditions for production needs to be maximized to reduce cost of production to a minimum. Conditions in the microbial bioreactor systems are more crucial than chemical synthesis reactors ⎯ simply due to the increase in simultaneous reactions occurring. Furthermore, the selection and concentration of compounds in the growth media can severely affect the growth rate. If a computational model can predict an optimal set of metabolites for media composition, this will save time, resources, and funds. </p>
-
<h2>Introduction</h2>
+
</html>
-
<p><b>What is it modeling?</b></p>
+
[[Image:Ucalgary2012_OptimalPoint.png|right|thumb|300px|Figure 2: Diagram of the relationship between growth rate and production rate, and the computed optimal growth rate.]]
-
<p>It models flux rate of metabolic pathways responding to different growth media conditions. The model is expected to generate an optimal set of metabolites that should be added to growth media in order to improve production rate. </p>
+
<html><h3>How does the program work?</h3>
-
<p><b>Why needs this kind of model?</b></p>
+
<p>Our program is built upon flux variability analysis applied by the functions in the Cobra Toolbox. It uses the published <i>E.coli</i> iAF1260 and <i>E.coli</i> core models provided from the  
-
<p>Same as chemical reactions need optimal environmental conditions to achieve maximum production rate, microbes also require optimized growth conditions to accomplish their tasks in maximum speed. In industrial scale, the optimal conditions for chemical reactions could increase the production rate in hundreds times which means millions of money. This principle also works even if taking biology pathways over chemical reactions. The difference is that in chemical way, the conditions of enzyme, pH, temperature and pressure are more important; however, in microbiological method, the conditions of growth media is more crucial. Further, the selection of media compounds is one of the most significant conditions for growth media. If a model can predict an optimal set of metabolites that need to be added into media, there will be tons of time and lots of money got saved. </p>
+
<a href="http://gcrg.ucsd.edu/Home">
-
 
+
Systems Biology Research Group (Dr. B. Palsson)
-
<p><b>How dose the program work?</b></p>
+
</a>
-
<p>This program is built upon constraint-based reconstruction analysis and flux variability analysis. It uses published E.coli (iAF1260) and E.coli core models as base chassis. Upon those models, we reconstruct new models of interests. Specifically, new reactions corresponding to the genes we try to engineer into E.coli from other organisms will be added into base chassis. By running flux variability analysis, program will give different sets of flux rates based on distinct constraints. Finally, the program will analysis those data with designed algorithm to generate a set of media compounds that is expected to accelerate production rate. </p>
+
at the University of San Diego. Using this as a base, we constructed reactions and metabolites for the hydrocarbon production components of OSCAR. Specifically, new reactions corresponding to the PetroBrick as well as the upgrading (desulfurization and denitrogenation) pathways were engineered into the <i>E.coli</i> base chassis. Running FVA will give an output correlated to a optimized growth rate. Moreover, this FVA algorithm running at optimal biomass flux rates can provide a range of outputs for your product of interest (<b>Figure 2</b>). This allows the user to select an optimum point where compound production is maximized.  Finally, the program will analysis the data with an algorithm to generate a set of media compounds that is expected to accelerate production rate. </p>
-
<br>
+
<h2>Algorithm</h2>
<h2>Algorithm</h2>
-
<p><b>Conceptual</b></p>
+
<h3>Conceptualization</h3>
-
<p>How to improve the products flux rates through data from FVA? The answer was remained unknown. There were only two things being noticed. One was that FVA could determine full range of numerical values for each reaction flux within the network and its output were able to use for analyze, and the other one was the biomass rate normally had trade-off relation with production rate. Since biomass rate reflects the growth condition, cell must have positive value of biomass flux rate in order to producing. On the other hand, the production flux rate should be higher than zero as well. This implied among all possible set of fluxes, the optimal flux set should locate a place where growth rate multiplies production rate is maximum.</p>
+
</html>
-
<p>Once the optimal flux rate of biomass was obtained, the value would be set as a new constraint of biomass. Then flux variability analysis would find out the full range of numerical values for each reaction flux within the network that was restricted to the new biological objective. </p>
+
[[Image:UCalgary2012_MaxAndMin.png|thumb|right|300px|Figure 3: Diagram of maximum and minimum production rates computed by flux variability analysis based on optimal growth rate.]]
-
<p>The differences of values for each reaction in a set of flux that maximized production rate and a set of flux that minimized production rate became interesting. By comparing two sets of fluxes based on visual maps, the results showed some reactions had higher flux rates in production maximum set than production minimum set, some were higher in production minimum set than production maximum set and some had opposite flux directions as most of biological reactions were reversible. In chemical, adding the amount of reactants would force the reactions equilibrium to move forwards, and adding the amount of products could drive the reactions equilibrium to go backwards. Consequently, the question becomes how to find out metabolites that need additional amount to improve the production rate.</p>
+
<html>
-
<p>One of the possible solutions could be comparing two sets of fluxes, determining differences of each reaction between two sets and changing constraints according to reaction needs. For example, if a metabolite needs more in production maximum set than production minimum set, then add more amount of this metabolite by change constraints to improve the production. However, in reality, cell could only uptake limited kinds of metabolites. Some metabolites were able to be produced by cell but not able to be absorbed from growth media. Hence, only the metabolites that had natural transporters in cell would count.</p>
+
<p>FVA can determine the full range of numerical values for each reaction flux within the network.  Additionally, it allows for a broader analysis of growth and production rates compared to FBA. Since biomass rate reflects the growth condition, cells must have a positive biomass flux rate to survive and proliferate. This positive value is indicative of a real cellular system since cells prefer to have increases in growth rather than increases in product output; however, our goal is to increase the production flux rate above a zero value. Therefore, among all the possible set of fluxes, the optimal flux is located where growth and production rate are at the optimum point for cell survival and compound output.</p>
-
<p>Last but not least, to improve production by adding more metabolites to growth media, the analysis should start from a model that was built upon glucose minimum growth media. </p>
+
-
<p><b>Concrete</b></p>
+
<p> Differences between the <i>maximum</i> and <i>minimum</i> production rates are compared for each reaction in a set of fluxes. Reactions that demonstrated higher flux in the maximum production rate versus the minimum production rate represent metabolites that could be increased to increase overall product output (<b>Figure 3</b>). This analysis generates a list of compounds that should increase product output if increased in concentration in the growth media.</p>
-
<p>Precondition: The original model is built with glucose minimum media.</p>
+
 
 +
 
 +
 
 +
<p>Note that not all compounds can be imported into the cell. Hence, only the metabolites that have native transporters in <i>E. coli</i> were considered for wet-lab verification. We used a model based on glucose minimal media; however, this could be applied to any media type.</p>
 +
 
 +
 
 +
<h3>Flowcart of the Model</h3>
 +
 
 +
</html>
 +
[[Image:Calgary_FluxModelMethod.jpg|center|thumb|600px|Figure 4: Methodology the program uses to identify metabolite hits which should be supplemented to the media of your organism to increase output.]]
 +
<html>
 +
 
 +
<p>Disclaimer: Original model built using minimum media with glucose.</p>
<p>1. Define relationship between growth rate and production rate.</p>
<p>1. Define relationship between growth rate and production rate.</p>
<p>2. Find out the optimal growth rate that can maximize the production.</p>
<p>2. Find out the optimal growth rate that can maximize the production.</p>
-
<p>3. Get the difference percentage of flux rate for each reaction between production maximum set and production minimum set.</p>
+
<p>3. Get the differences in percentage of flux rate for each reaction between production maximum and minimum set.</p>
-
<p>4. Collect all reactions have difference percentage between two sets that exceed threshold.</p>
+
<p>4. Collect all reactions with different percentages between two sets that exceed the user input threshold. (The input threshold determines the size of the difference in flux that the user is interested in.  We used a 100% difference in our model.).</p>
-
<p>5. Score each compound in all collected reactions (Initial score is zero for each compound). </p>
+
<p>5. Score each compound in all collected reactions, with the initial score of zero for each compound (Scoring is determined by the difference in the value of a particular compound when the flux sets of the desired compound production are maximized and minimized as determined by FVA). This process is additive if the particular compound is found in more than one reaction (i.e., ATP) and only includes the reactions identified from step 4.</p>
-
<p><span style="padding-left:20px">5.1 The difference of flux rates of one reaction from production maximum set to production minimum set is added to the score for all reactants of this reaction. </span></p>
+
-
        <p><span style="padding-left:20px">5.2 The difference of flux rates of one reaction from production minimum set to production maximum set is added to the score for all products of this reaction.</span></p>
+
-
        <p><span style="padding-left:20px">5.3 Repeat 3.1 to 3.2 till all collected reactions are analyzed.</span></p>
+
<p>6. Determine whether compounds with positive scores have natural transporters in cell. If so, mark the compound as candidate.</p>
<p>6. Determine whether compounds with positive scores have natural transporters in cell. If so, mark the compound as candidate.</p>
-
<p>7. Add each candidate to growth media, and run FVA under optimal growth rate computed in Step 2. Compare the production rate from novel model to that from raw model, if the rate is improved, mark as effector.</p>
+
<p>7. Add each candidate to the growth media and the run FVA under optimal growth rate computed in Step 2. Compare the production rate from novel model (step 7) to that from the raw model (step 2), if the rate is improved, mark the compound as an effector.</p>
<br>
<br>
 +
<h2>Graphical user interface (GUI) development - Creating OSCAR Optimizer</h2>
 +
<p>The GUI allows for easier use of our program by everyone in iGEM and beyond. This will enable individuals to use our program who are inexperienced with programming programs such as MATLAB and/or the Cobra Toolbox.  We developed the GUI using the GUIDE program designed in MATLAB.  For more information please see our manual.</p>
 +
<br>
 +
 +
<h2>Code</h2>
 +
<p>This application is a Matlab extension that runs on top of the Cobra Toolbox and SBML Toolbox. To run the application, one must have the Cobra Toolbox and SBML Toolbox installed.  SBML Toolbox can download from <a href="http://sbml.org/Software/SBMLToolbox">SBML.org</a> or <a href="http://sourceforge.net/projects/sbml/files/SBMLToolbox/4.1.0/SBMLToolbox-4.1.0.zip/download">here</a>. Cobra Toolbox can download from <a href="http://opencobra.sourceforge.net/openCOBRA/Welcome.html">openCOBRA</a> or <a href="http://sourceforge.net/projects/opencobra/files/cobra/cobra_2.0.5.zip/download">here</a>. Application Package and source code<a href="https://static.igem.org/mediawiki/2012/5/50/UCalgary2012_OSCAR_optimizer_v1.zip"> download </a>.</p>
 +
<br>
<h2>Demo</h2>
<h2>Demo</h2>
 +
<p>Below is an uploaded tutorial to show anyone how to use the GUI of the OSCAR Optimizer.  The GUI interface allows for easy building of different synthetic constructs into the <i>E. coli</i> network but this could based on any model from any organism that is available in SBML format.</p>
 +
<div align="center">
<div align="center">
-
<iframe width="600" height="480" align="center" src="http://www.youtube.com/watch?v=9oHJhQs5wQk" frameborder="0" allowfullscreen></iframe>
+
<iframe width="420" height="315" src="http://www.youtube.com/embed/9oHJhQs5wQk" frameborder="0" allowfullscreen></iframe>
</div>
</div>
-
<br>
 
 +
</html>
 +
<h4><center>Screen Shots of Our Application</center></h4>
 +
[[Image:UCalgary2012_RunTimeBuild.png|350px|left|]]
 +
[[Image:UCalgary2012_DemoOutput.png|350px|right|]]
-
<h2>Drawbacks</h2>
+
[[Image:UCalgary2012_BuildOutput.png|400px|center|]]
-
<p>This application is built upon Cobra Toolbox, and Cobra Toolbox is an application of SBML Toolbox. As consequence, any flaws in Cobra Toolbox and SBML Toolbox will affect this application.</P>
+
[[Image:UCalgary2012_AnalysisOutput.png|400px|center|]]
-
<p>In this program, pathways added to base chassis model (E. Coli iAF1260) contain constraints only relied on Stoichiometric Matrix such as Stoichiometric coefficients, lower bounds and upper bounds of reactions. They are lack of genetic and enzymatic regulation, which makes the connections between reactions in the network much weaker than those in real. The missing rules could lead to program outputs inaccuracy results.</P>
+
<html>
-
<p>At current stage, the algorithm can only pick metabolites with natural transporters in cell. Many other intermediate metabolites are ignored. The algorithm has no power to trace the intermediate metabolites back to initial metabolites and take those initial metabolites into account. This lack of power could again make the results weak.</P>
+
<br>
<br>
 +
<a name="Flux"></a><h2>Wet-lab validation of OSCAR Optimizer</h2>
 +
<p>We tested our OSCAR Optimizer FVA Model in the wet-lab by assessing if a suite of identified compounds  would increase the output of the Petrobrick. Therefore we ran the model with the <i>AAR</i> and <i>ADC</i> gene components of the Petrobrick system and looked at the predicted metabolites for growth media supplementation (<b>Figure 5</b>).</p>
-
<h2>Validation</h2>
+
</html>[[File:Calgary ModelingPetroOutput.jpg|center|thumb|600px|Figure 5: OSCAR Optimizer output for increase hydrocarbon production in OSCAR. This program identified a series of metabolites that are predicted to increase hydrocarbon output, including fructose and AMP with the highest outputs suggesting that their presence should contribute to increased product output the most.]]<html>
-
<br>
+
 +
<p>Once these compounds were identified, we set-up an assay where we supplemented minimal M9 media and glucose with each of the compounds alone, or in combination. Compounds were added at concentrations of 50 mM except for Ethanol (2.5% v/v), AMP (100mg/L), and L-aspartate (100mg/L).  As a positive control, instead of using minimal media a solution of 50:50 LB:Washington Production Media (see
 +
<a href="https://2011.igem.org/Team:Washington/alkanebiosynthesis">
 +
2011 iGEM University of Washington Protocols
 +
</a>
 +
) was used to represent normal production.  Additionally, we found that the compound formate was predicted to NOT increase hydrocarbons production and this represented a chemical negative control.  All compounds were added to <i>E. coli</i> culture containing the PetroBrick at an OD<sub>600</sub> of ~0.05 in the first experiment, and a higher OD<sub>600</sub> of 0.400 for the second experiment and grown for 72 hours at 37<sup>o</sup>C.  Cultures were then sonicated after an OD<sub>600</sub> measurement and any produced hydrocarbons were extracted by 1 mL ethyl acetate. This was quantitated based on the peak area of a C15 hydrocarbon product as described in the
 +
<a href="https://2012.igem.org/Team:Calgary/Project/OSCAR/Decarboxylation">Decarboxylation</a> section (also see
 +
<a href="https://2012.igem.org/Team:Calgary/Notebook/Protocols">Protocols</a>
 +
for relevant procedures).  The tabulated results of the quantitated yields are shown in <b>Figure 6</b>. </p>
-
<h2>Source Code</h2>
+
</html>[[File:Calgary ModellingWetlabDataFinal.png|center|thumb|600px|Figure 6: Model validation, showing the relative abundance of C15 hydrocarbons of <i>E. coli</i> containing the Petrobrick in combination with different compounds as growth media supplements.  Data is normalized to the positive control (50:50 LB:Washington Production Media) and the relative production of hydrocarbons is displayed relative to optical density measurements.  This data suggests hydrocarbon production can be increased or decreased if the Petrobrick is incubated with specific compounds.  Dashed line represents the average output of hydrocarbons in the positive control. PMFA (Pyruvate, Malate, Fumarate, Aspartic Acid) and G+P (Glycine and Pyruvate). ]]<html>
-
<br>
+
-
<h2>Documents</h2>
+
<p> These hydrocarbon production results indicate there is variability in the output of hydrocarbons from the PetroBrick with different compounds (<b>Figure 6</b>). Interestingly, five of the compounds demonstrated production levels higher than that of the minimal media control, suggesting our model could be used for predicting compounds to increase metabolite production.  Although, AMP and fructose both decreased output despite having the highest predicted increases in hydrcarbon output, suggesting that our model may produce false positives and there is some error in assessing these predictions.  Most excitingly, addition of two compounds, pyruvate and glycine, had increased the relative number of hydrocarbons higher than the complex media of the positive control.  As noted by the multiple experiments, there are some differences in the relative output of hydrocarbons in each experiment, however, the trend for each compound is relatively consistent for the different samples.  Some of the runs were not able to be repeated in the second assay, which is the reason for some of the lack of data in Figure 6.  <b>These results suggest our OSCAR Optimizer model can optimize OSCAR (along with other metabolic systems); however, these predicted candidates need to be tested in the wet-lab to ensure the model is accurate.</b>
 +
 
 +
<h2> Modelling informing our Killswitch </h2>
 +
 
 +
<p> Although we had initially rejected the idea of a killswitch based on <b>auxotrophy</b> due to the issue of cost, the results of our model made us think that this might be more feasible than we had initially thought.  As glycine alone in minimal media was able to give us higher alkane production than Washington’s production media, supplementing our system with glycine, achieving higher output may actually be feasible.  We saw auxotrophy as a method of creating <b>layers</b> in our killswitch system, where we could use it in tandem with the inducible kill systems that we have already been developing.  We obtained a glycine knockout stain of <i>E. coli</i> and have begun some initial <a href="https://2012.igem.org/Team:Calgary/Project/HumanPractices/Killswitch/Regulation">characterization data</a>.  This was a great example of how our modelling project was actually able to <b>inform</b> other parts of our project.  In order to validate that this system was feasible for use we performed three key experiments:
 +
 
 +
<ul>
 +
<li>Perform a glycine requirement assay to determine the survivability of the knockout strain in our minimal media's used with OSCAR <.</li>
 +
<li>Determine if the Petrobrick can be used with a similar activity as in the wild type strain</li>
 +
<li>Use the auxotrophic strain in conjunction with our killswitch to show that these can act together.</li>
 +
</ul>
 +
<p><a href="https://2012.igem.org/Team:Calgary/Project/Synergy"> Click here</a> to learn more about these assays. </p>
 +
 
 +
<h2>Drawbacks</h2>
 +
<p>This application is built upon the Cobra Toolbox, and the SBML Toolbox. As a consequence, any flaws in the Cobra Toolbox and SBML Toolbox will affect this application.</P>
 +
<p>In this program, pathways added to base chassis model (<i>E. Coli</i> iAF1260) which contains constraints that rely on the Stoichiometric Matrix (such as Stoichiometric Coefficients), lower bounds and upper bounds of reactions. These pathways lack genetic and enzymatic regulation, making the connections between reactions in the network much weaker than those in the wet-lab. These missing components may lead to inaccurate results.</p>
 +
 
 +
<p>At the current stage, our algorithm can only pick metabolites with natural transporters in cell. Therefore, many other intermediate metabolites are ignored. The algorithm has no power to trace the intermediate metabolites back to initial metabolites and take those initial metabolites into account. <b>This model does not take into consideration rate limiting steps of reactions, toxic intermediate accumulation, or any form of regulation of the enzymes (e.g., negative feedback).  The OSCAR Optimizer should act as a starting point to identify any limitation of available endogenous metabolites in <i>E. coli</i> to increase the starting compound levels of your synthetic circuit for optimal production.</b></P>
<br>
<br>
 +
 +
</html>
 +
[[Image:UCalgary2012_Manual.png|left|200px|link=https://static.igem.org/mediawiki/2012/d/d5/UCalgary2012_OSCAROptimizerManual.pdf]]
 +
<html><h2>Documents</h2>
 +
<p><a href="https://static.igem.org/mediawiki/2012/d/d5/UCalgary2012_OSCAROptimizerManual.pdf">Download</a> the OSCAR Optimizer Manual for more information on how to use the program!</p>
</html>
</html>
}}
}}

Latest revision as of 03:43, 27 October 2012

Hello! iGEM Calgary's wiki functions best with Javascript enabled, especially for mobile devices. We recommend that you enable Javascript on your device for the best wiki-viewing experience. Thanks!

Flux-Variability Analysis for Optimization

UCalgary2012 OSCAR Flux Analysis Low-Res.png

Flux-variability analysis (FVA) was applied to optimize the bioreactor system of OSCAR for the newly incorporated metabolic pathways (Decarboxylation, Decatecholization, Desulfurization, and Denitrogenation). FVA combines the framework of metabolic pathways with experimental enzymatic data to provide a computational platform for predicting what changes in metabolite levels will result in increased end-product. We developed a model using MATLAB computer language for predicting what metabolites could be added to your growth media to increase production of hydrocarbons in the E. coli chassis. We also created a graphical user interface called the OSCAR Optimizer to make the the FVA program user-friendly and to allow current and future iGEM teams/scientists to input their own synthetic pathways into the FVA for analysis. Finally, we validated this model in the wet-lab by optimizing the PetroBrick (BBa_K590025) system to increase hydrocarbon production, thereby saving time, and resources.

Click here to download the MATLAB package to test our model!

Click here to view our OSCAR Optimizer Manual.

Background

What is Flux Analysis?

Flux Balance Analysis (FBA) uses linear programming of metabolic networks to convert each metabolite into a mathematical coefficient. These coefficients can be connected to each other by changes associated with each enzymatic step of the pathway. FBA can then apply a mathematical method to examine how metabolites relate to each other in the network and allows the user to make generalized predictions for the growth of the organism, metabolite levels, and product output inside the cell. Flux Variability Analysis (FVA) is an extension of FBA that determines the range of reactions that result in an optimal flux through the metabolic network. Unlike FBA, this allows the user to not only determine the most optimized network for increasing cell growth, but also other parameters like product output.

What are the constraints in the FVA model?

As illustrated in Figure 1, metabolic networks can be encoded as stoichiometric matrices, in which each row represents a specific metabolite and each column represents a biochemical reaction. The entries in each column represent the stoichiometric coefficients of the metabolites in the reaction. Metabolites that are consumed or produced have a negative or positive, respectively.

Figure 1: Flux Balance Analysis. FBA uses metabolic network (left) and simulates the connections in the network as a linear algebra matrix (right). Each metabolite is listed vertically in the table and each reaction is listed horizontally. Based on the metabolites involved in each reaction, the state of the system is affected by each metabolite change.

Why use FVA as opposed to FBA?

Biological systems often contain redundancies that contribute to their robustness. However, FBA only returns a single flux distribution that corresponds to maximal bacterial growth under the given growth conditions, even if other optimal parameters may exist. FVA is capable of examining these redundancies by calculating the full range of numerical values for each reaction flux in a network. Thus, FVA can be employed to study the entire range of achievable cellular functions, as well as the redundancy in optimal phenotypes. FVA can also examine different ranges of bacterial growth vs. product output which is valuable in assessing validity of models in the wet-lab.


Using FVA to optimize OSCAR

What Are We Trying To Model?

FVA provides a method to modulate the inputs of endogenous metabolic pathways, but could also investigate what chemicals are needed in the growth media to upregulate an introduced pathway in E. coli. In the latter case, we need to specifically model the flux rate of metabolic pathways responding to different growth media conditions and generate an optimal set of metabolites that should be added to improve production rate. Development of a tool to model the optimal flux would also benefit numerous iGEM teams and research labs who engineer bacteria with novel pathways.

How could systems like OSCAR benefit from this model?

Much like chemical reactions which need optimal environmental conditions to achieve maximum production rate, microbes require certain conditions to accomplish cellular tasks at maximum speed. During industrial scale up, the optimal conditions for production needs to be maximized to reduce cost of production to a minimum. Conditions in the microbial bioreactor systems are more crucial than chemical synthesis reactors ⎯ simply due to the increase in simultaneous reactions occurring. Furthermore, the selection and concentration of compounds in the growth media can severely affect the growth rate. If a computational model can predict an optimal set of metabolites for media composition, this will save time, resources, and funds.

Figure 2: Diagram of the relationship between growth rate and production rate, and the computed optimal growth rate.

How does the program work?

Our program is built upon flux variability analysis applied by the functions in the Cobra Toolbox. It uses the published E.coli iAF1260 and E.coli core models provided from the Systems Biology Research Group (Dr. B. Palsson) at the University of San Diego. Using this as a base, we constructed reactions and metabolites for the hydrocarbon production components of OSCAR. Specifically, new reactions corresponding to the PetroBrick as well as the upgrading (desulfurization and denitrogenation) pathways were engineered into the E.coli base chassis. Running FVA will give an output correlated to a optimized growth rate. Moreover, this FVA algorithm running at optimal biomass flux rates can provide a range of outputs for your product of interest (Figure 2). This allows the user to select an optimum point where compound production is maximized. Finally, the program will analysis the data with an algorithm to generate a set of media compounds that is expected to accelerate production rate.

Algorithm

Conceptualization

Figure 3: Diagram of maximum and minimum production rates computed by flux variability analysis based on optimal growth rate.

FVA can determine the full range of numerical values for each reaction flux within the network. Additionally, it allows for a broader analysis of growth and production rates compared to FBA. Since biomass rate reflects the growth condition, cells must have a positive biomass flux rate to survive and proliferate. This positive value is indicative of a real cellular system since cells prefer to have increases in growth rather than increases in product output; however, our goal is to increase the production flux rate above a zero value. Therefore, among all the possible set of fluxes, the optimal flux is located where growth and production rate are at the optimum point for cell survival and compound output.

Differences between the maximum and minimum production rates are compared for each reaction in a set of fluxes. Reactions that demonstrated higher flux in the maximum production rate versus the minimum production rate represent metabolites that could be increased to increase overall product output (Figure 3). This analysis generates a list of compounds that should increase product output if increased in concentration in the growth media.

Note that not all compounds can be imported into the cell. Hence, only the metabolites that have native transporters in E. coli were considered for wet-lab verification. We used a model based on glucose minimal media; however, this could be applied to any media type.

Flowcart of the Model

Figure 4: Methodology the program uses to identify metabolite hits which should be supplemented to the media of your organism to increase output.

Disclaimer: Original model built using minimum media with glucose.

1. Define relationship between growth rate and production rate.

2. Find out the optimal growth rate that can maximize the production.

3. Get the differences in percentage of flux rate for each reaction between production maximum and minimum set.

4. Collect all reactions with different percentages between two sets that exceed the user input threshold. (The input threshold determines the size of the difference in flux that the user is interested in. We used a 100% difference in our model.).

5. Score each compound in all collected reactions, with the initial score of zero for each compound (Scoring is determined by the difference in the value of a particular compound when the flux sets of the desired compound production are maximized and minimized as determined by FVA). This process is additive if the particular compound is found in more than one reaction (i.e., ATP) and only includes the reactions identified from step 4.

6. Determine whether compounds with positive scores have natural transporters in cell. If so, mark the compound as candidate.

7. Add each candidate to the growth media and the run FVA under optimal growth rate computed in Step 2. Compare the production rate from novel model (step 7) to that from the raw model (step 2), if the rate is improved, mark the compound as an effector.


Graphical user interface (GUI) development - Creating OSCAR Optimizer

The GUI allows for easier use of our program by everyone in iGEM and beyond. This will enable individuals to use our program who are inexperienced with programming programs such as MATLAB and/or the Cobra Toolbox. We developed the GUI using the GUIDE program designed in MATLAB. For more information please see our manual.


Code

This application is a Matlab extension that runs on top of the Cobra Toolbox and SBML Toolbox. To run the application, one must have the Cobra Toolbox and SBML Toolbox installed. SBML Toolbox can download from SBML.org or here. Cobra Toolbox can download from openCOBRA or here. Application Package and source code download .


Demo

Below is an uploaded tutorial to show anyone how to use the GUI of the OSCAR Optimizer. The GUI interface allows for easy building of different synthetic constructs into the E. coli network but this could based on any model from any organism that is available in SBML format.

Screen Shots of Our Application

UCalgary2012 RunTimeBuild.png
UCalgary2012 DemoOutput.png
UCalgary2012 BuildOutput.png
UCalgary2012 AnalysisOutput.png

Wet-lab validation of OSCAR Optimizer

We tested our OSCAR Optimizer FVA Model in the wet-lab by assessing if a suite of identified compounds would increase the output of the Petrobrick. Therefore we ran the model with the AAR and ADC gene components of the Petrobrick system and looked at the predicted metabolites for growth media supplementation (Figure 5).

Figure 5: OSCAR Optimizer output for increase hydrocarbon production in OSCAR. This program identified a series of metabolites that are predicted to increase hydrocarbon output, including fructose and AMP with the highest outputs suggesting that their presence should contribute to increased product output the most.

Once these compounds were identified, we set-up an assay where we supplemented minimal M9 media and glucose with each of the compounds alone, or in combination. Compounds were added at concentrations of 50 mM except for Ethanol (2.5% v/v), AMP (100mg/L), and L-aspartate (100mg/L). As a positive control, instead of using minimal media a solution of 50:50 LB:Washington Production Media (see 2011 iGEM University of Washington Protocols ) was used to represent normal production. Additionally, we found that the compound formate was predicted to NOT increase hydrocarbons production and this represented a chemical negative control. All compounds were added to E. coli culture containing the PetroBrick at an OD600 of ~0.05 in the first experiment, and a higher OD600 of 0.400 for the second experiment and grown for 72 hours at 37oC. Cultures were then sonicated after an OD600 measurement and any produced hydrocarbons were extracted by 1 mL ethyl acetate. This was quantitated based on the peak area of a C15 hydrocarbon product as described in the Decarboxylation section (also see Protocols for relevant procedures). The tabulated results of the quantitated yields are shown in Figure 6.

Figure 6: Model validation, showing the relative abundance of C15 hydrocarbons of E. coli containing the Petrobrick in combination with different compounds as growth media supplements. Data is normalized to the positive control (50:50 LB:Washington Production Media) and the relative production of hydrocarbons is displayed relative to optical density measurements. This data suggests hydrocarbon production can be increased or decreased if the Petrobrick is incubated with specific compounds. Dashed line represents the average output of hydrocarbons in the positive control. PMFA (Pyruvate, Malate, Fumarate, Aspartic Acid) and G+P (Glycine and Pyruvate).

These hydrocarbon production results indicate there is variability in the output of hydrocarbons from the PetroBrick with different compounds (Figure 6). Interestingly, five of the compounds demonstrated production levels higher than that of the minimal media control, suggesting our model could be used for predicting compounds to increase metabolite production. Although, AMP and fructose both decreased output despite having the highest predicted increases in hydrcarbon output, suggesting that our model may produce false positives and there is some error in assessing these predictions. Most excitingly, addition of two compounds, pyruvate and glycine, had increased the relative number of hydrocarbons higher than the complex media of the positive control. As noted by the multiple experiments, there are some differences in the relative output of hydrocarbons in each experiment, however, the trend for each compound is relatively consistent for the different samples. Some of the runs were not able to be repeated in the second assay, which is the reason for some of the lack of data in Figure 6. These results suggest our OSCAR Optimizer model can optimize OSCAR (along with other metabolic systems); however, these predicted candidates need to be tested in the wet-lab to ensure the model is accurate.

Modelling informing our Killswitch

Although we had initially rejected the idea of a killswitch based on auxotrophy due to the issue of cost, the results of our model made us think that this might be more feasible than we had initially thought. As glycine alone in minimal media was able to give us higher alkane production than Washington’s production media, supplementing our system with glycine, achieving higher output may actually be feasible. We saw auxotrophy as a method of creating layers in our killswitch system, where we could use it in tandem with the inducible kill systems that we have already been developing. We obtained a glycine knockout stain of E. coli and have begun some initial characterization data. This was a great example of how our modelling project was actually able to inform other parts of our project. In order to validate that this system was feasible for use we performed three key experiments:

  • Perform a glycine requirement assay to determine the survivability of the knockout strain in our minimal media's used with OSCAR <.
  • Determine if the Petrobrick can be used with a similar activity as in the wild type strain
  • Use the auxotrophic strain in conjunction with our killswitch to show that these can act together.

Click here to learn more about these assays.

Drawbacks

This application is built upon the Cobra Toolbox, and the SBML Toolbox. As a consequence, any flaws in the Cobra Toolbox and SBML Toolbox will affect this application.

In this program, pathways added to base chassis model (E. Coli iAF1260) which contains constraints that rely on the Stoichiometric Matrix (such as Stoichiometric Coefficients), lower bounds and upper bounds of reactions. These pathways lack genetic and enzymatic regulation, making the connections between reactions in the network much weaker than those in the wet-lab. These missing components may lead to inaccurate results.

At the current stage, our algorithm can only pick metabolites with natural transporters in cell. Therefore, many other intermediate metabolites are ignored. The algorithm has no power to trace the intermediate metabolites back to initial metabolites and take those initial metabolites into account. This model does not take into consideration rate limiting steps of reactions, toxic intermediate accumulation, or any form of regulation of the enzymes (e.g., negative feedback). The OSCAR Optimizer should act as a starting point to identify any limitation of available endogenous metabolites in E. coli to increase the starting compound levels of your synthetic circuit for optimal production.


UCalgary2012 Manual.png

Documents

Download the OSCAR Optimizer Manual for more information on how to use the program!