Team:UNAM-Genomics Mexico/Modeling/FBA

From 2011.igem.org

(Difference between revisions)
Line 8: Line 8:
Flux Balance Analysis (FBA) can serve to explore the fluxes of a given metabolic reconstruction. In this case, we wanted to determine the level and extent of interaction of our added pathway with the core metabolism. Since our chassis, R.etli, has two "flavors" (free living organism & plant symbiont), we were curious as to weather our transgenic system would remain functional under symbiont form. Moreover, since a key aspect of the project was Nitrogen fixation, we wanted to ensure said pathway was as functional as it could be.
Flux Balance Analysis (FBA) can serve to explore the fluxes of a given metabolic reconstruction. In this case, we wanted to determine the level and extent of interaction of our added pathway with the core metabolism. Since our chassis, R.etli, has two "flavors" (free living organism & plant symbiont), we were curious as to weather our transgenic system would remain functional under symbiont form. Moreover, since a key aspect of the project was Nitrogen fixation, we wanted to ensure said pathway was as functional as it could be.
-
Furthermore, key researchers at our University warned us against the project due to the widespread presence of endogenous hydrogenases (the word "impossible" was often used, which only got us rebelliously stuck on proving them wrong). These native enzymes play the opposite role of the activity we were interested in, they '''consume''' molecular hydrogen instead of '''producing''' it. Apparently, since Nitrogen fixation produces molecular hydrogen, some Rhizobial species developed ''capture hydrogenases'' to recycle the protons into the all-important proton gradient that feeds cellular machinery. However, our particular strain, CFN42, doesn't have any ''capture hydrogenase''. Nonetheless, since nature developed pathways to do the exact opposite of what we're doing, we were interested in finding out any toxic or deleterious effects our system might have. Finaly, out goals can be stated as being:
+
Furthermore, key researchers at our University warned us against the project due to the widespread presence of endogenous hydrogenases (the word "impossible" was often used, which only got us rebelliously stuck on proving them wrong). These native enzymes play the opposite role of the activity we were interested in, they '''consume''' molecular hydrogen instead of '''producing''' it. Apparently, since Nitrogen fixation produces molecular hydrogen, some Rhizobial species developed ''capture hydrogenases'' to recycle the protons into the all-important proton gradient that feeds cellular machinery. However, our particular strain, CFN42, doesn't have any ''capture hydrogenase''. Nonetheless, since nature developed pathways to do the exact opposite of what we're doing, we were interested in finding out any toxic or deleterious effects our system might have. Therefore, out goals can be stated as being:
# Under symbiont form, our transgenic chassis is capable of Hydrogen production?
# Under symbiont form, our transgenic chassis is capable of Hydrogen production?
# To what extent does the Hydrogen production affect core metabolism, and Nitrogen fixation?
# To what extent does the Hydrogen production affect core metabolism, and Nitrogen fixation?
Line 15: Line 15:
==Theoretical Background==
==Theoretical Background==
-
FBA is part of a larger set of models known as the Constraint Based Models (CBM). Since cells may adopt an astronomical level of possible solutions for a given set of conditions, CBM models create constrains on the solution space to limit the complexity of the solution and render it easier to compute, interpret, and/or understand.
+
FBA is part of a larger set of models known as the Constraint Based Analysis (CBA). Since cells may adopt an astronomical level of possible solutions for a given set of conditions, CBM models create constrains on the solution space to limit the complexity of the solution and render it easier to compute, interpret, and/or understand.
-
 
+
-
FBA starts at the assumption that pathways strive towards stability. Thus, it assumes the cellular chemoton will adjust whatever it has to do in order to maintain chemical stability. In other words, this model does not require pesky kinetic constants, something that made our lives easier. For some ''light reading'' on this, you can consult "Properties of Reconstructed Networks", by Bernhard Palsson '''[Ref1]'''.
+
 +
FBA starts at the assumption that pathways strive towards homeostasis. Thus, it assumes the cellular chemoton will adjust whatever it has to do in order to maintain chemical stability. In other words, this model does not require pesky kinetic constants, something that made our lives easier. For some ''light'' reading on this, you can consult [http://dx.doi.org/10.2277/0521859034 "Systems Biology: Properties of Reconstructed Networks", by Bernhard Palsson]. Or you can always go ask the All-Knowing-Oracle [http://en.wikipedia.org/wiki/Flux_balance_analysis here].
==The Simulation==
==The Simulation==
-
We performed the FBA calculations using a particular toolbox in iGEM's favorite program: MATLAB.
+
We performed the FBA calculations using a particular toolbox in iGEM's favorite program: MATLAB. In particular, we used the [http://opencobra.sourceforge.net/openCOBRA/Welcome.html COnstrained Based Reconstruction & Analysis toolbox], originally by the Palsson Lab. I have to confess, the poor documentation made our life difficult, but eventually we got it working. It appears the toolbox has a lot of functions for various CBA. At first we decided to parse the data we had into The Stochiometric Matrix, and we did it. It was then we discovered the elusive '''parseModel''' function, after which we curated manually our data and set to re-analyzing and tweaking the metabolic reconstruction to our heart's content.
-
 
+
-
===The Tool===
+
-
We used the COBRA toolbox '''[Ref2]''' to simulate our system.
+
===The Different Models===
===The Different Models===
Line 42: Line 38:
==References==
==References==
-
# Properties of Reconstructed Networks, Bernhard Palsson
+
<references />
-
# [http://opencobra.sourceforge.net/openCOBRA/Welcome.html openCOBRA Project]
+
 
}}
}}

Revision as of 23:58, 25 September 2011

UNAM-Genomics_Mexico

Flux Balance Analysis

Our synthetic pathway produces Hydrogen. Our chassis fixes Nitrogen, which is influenced by Hydrogen availability. Therefore, we were keenly interested in finding out how our synthetic pathway would interact with the host metabolism. After some exploratory consultation, we determined Flux Balance Analysis was an effective tool for this.


Introduction

Flux Balance Analysis (FBA) can serve to explore the fluxes of a given metabolic reconstruction. In this case, we wanted to determine the level and extent of interaction of our added pathway with the core metabolism. Since our chassis, R.etli, has two "flavors" (free living organism & plant symbiont), we were curious as to weather our transgenic system would remain functional under symbiont form. Moreover, since a key aspect of the project was Nitrogen fixation, we wanted to ensure said pathway was as functional as it could be.

Furthermore, key researchers at our University warned us against the project due to the widespread presence of endogenous hydrogenases (the word "impossible" was often used, which only got us rebelliously stuck on proving them wrong). These native enzymes play the opposite role of the activity we were interested in, they consume molecular hydrogen instead of producing it. Apparently, since Nitrogen fixation produces molecular hydrogen, some Rhizobial species developed capture hydrogenases to recycle the protons into the all-important proton gradient that feeds cellular machinery. However, our particular strain, CFN42, doesn't have any capture hydrogenase. Nonetheless, since nature developed pathways to do the exact opposite of what we're doing, we were interested in finding out any toxic or deleterious effects our system might have. Therefore, out goals can be stated as being:

  1. Under symbiont form, our transgenic chassis is capable of Hydrogen production?
  2. To what extent does the Hydrogen production affect core metabolism, and Nitrogen fixation?
  3. Is there a toxic effect derived from Hydrogen production?


Theoretical Background

FBA is part of a larger set of models known as the Constraint Based Analysis (CBA). Since cells may adopt an astronomical level of possible solutions for a given set of conditions, CBM models create constrains on the solution space to limit the complexity of the solution and render it easier to compute, interpret, and/or understand.

FBA starts at the assumption that pathways strive towards homeostasis. Thus, it assumes the cellular chemoton will adjust whatever it has to do in order to maintain chemical stability. In other words, this model does not require pesky kinetic constants, something that made our lives easier. For some light reading on this, you can consult [http://dx.doi.org/10.2277/0521859034 "Systems Biology: Properties of Reconstructed Networks", by Bernhard Palsson]. Or you can always go ask the All-Knowing-Oracle [http://en.wikipedia.org/wiki/Flux_balance_analysis here].

The Simulation

We performed the FBA calculations using a particular toolbox in iGEM's favorite program: MATLAB. In particular, we used the [http://opencobra.sourceforge.net/openCOBRA/Welcome.html COnstrained Based Reconstruction & Analysis toolbox], originally by the Palsson Lab. I have to confess, the poor documentation made our life difficult, but eventually we got it working. It appears the toolbox has a lot of functions for various CBA. At first we decided to parse the data we had into The Stochiometric Matrix, and we did it. It was then we discovered the elusive parseModel function, after which we curated manually our data and set to re-analyzing and tweaking the metabolic reconstruction to our heart's content.

The Different Models

We also made different models to check how they each behaved

Wildtype

The standard and everyday chassis, our Institute's pet.

Transgenic full

Transgenic without H2 sink

Transgenic without H2 sink nor H2 production

Transgenic without H2 sink nor H2 production nor PFOR activity

Transgenic full without PHB carbon storage pathway

The Results

0x00000

References

<references />