Genome Scale Modeling
One of our main objectives for this project is to build a comprehensive, genome-scale model that details the function of the lux system and its interaction with the host organism, Saccharomyces Cerevisiae (s288c).
One of our main objectives for this project is to build a comprehensive, genome-scale model that details the function of the lux system and its interaction with the host organism, Saccharomyces Cerevisiae (s288c).
There are currently many M-Models available in literature, mostly detailing the networks of different strains of bacteria and the simpler eukaryotes [2]. The yeast strain of Saccharomyces Cerevisiae (s288c) was one of the first to be converted into a network model. For this and for reasons of experimental simplicity, we chose this strain as our base model.
The UC system-wide biochemically, genetically, and genomically (BiGG) database features two different networks for the genome of Saccharomyces Cerevisiae. The latest model reconstruction, iND750, consists of 1059 metabolites, 1266 reactions, and 750 genes [2]. In early 2007, a community of researchers came together to standardize the metabolic map of Saccharomyces Cerevisiae. They collaborated to periodically update the model and make it available in the public domain (yeast.sf.net) for the years that followed. The latest version to date, version 7, includes reactions published by Aung et al. 2013 [3].
The UC system-wide biochemically, genetically, and genomically (BiGG) database features two different networks for the genome of Saccharomyces Cerevisiae. The latest model reconstruction, iND750, consists of 1059 metabolites, 1266 reactions, and 750 genes [2]. In early 2007, a community of researchers came together to standardize the metabolic map of Saccharomyces Cerevisiae. They collaborated to periodically update the model and make it available in the public domain (yeast.sf.net) for the years that followed. The latest version to date, version 7, includes reactions published by Aung et al. 2013 [3].
Metabolic Control Analysis and Flux Balance Analysis
The UC system-wide biochemically, genetically, and genomically (BiGG) database features two different networks for the genome of Saccharomyces Cerevisiae. The latest model reconstruction, iND750, consists of 1059 metabolites, 1266 reactions, and 750 genes [2]. In early 2007, a community of researchers came together to standardize the metabolic map of Saccharomyces Cerevisiae. They collaborated to periodically update the model and make it available in the public domain (yeast.sf.net) for the years that followed. The latest version to date, version 7, includes reactions published by Aung et al. 2013 [3].
However, FBA has many limitations. FBA cannot predict metabolic concentrations, because it does not utilize kinetic parameters. Also, it cannot account for regulatory agents such as those involved in transcription or translation of the genome [7,9]].
Constraint-Based Modeling
We used the COBRA Toolbox plugin in Matlab to edit and analyze our SBML model [1]. The steady-state assumption that is employed in flux balance analysis is one example of a constraint. Constraints define the allowable solution space in flux balance analysis and accounts for many reaction properties that cannot be encoded by simple stoichiometry. Metabolic constraints, such as flux coefficients, define the physico-chemical boundaries of a reaction. The lower bound of the flux coefficient also tells us if the reaction is reversible or irreversible [6].
Space constraints, such as the use of compartments, define the physical boundaries on which metabolites can react with one another. There are many other types of constraints that can be embedded into a genome-scale model such as regulatory constraints, thermodynamic constraints, and temporal constraints. All parameters seek to eliminate infeasible solutions and define the solution space as closely as possible [6,8].
The most important constraint is the objective function. The COBRA Toolbox allows us to maximize or minimize the flux through certain reactions of interest [1]. For instance, setting the objective to biomass production allows us to predict the growth rate of a yeast cell once it reaches steady-state. For our purposes, maximizing the flux through the light reaction was set as the main objective.
Results
In both models, we first defined a minimum media to mimic experimental conditions. The components of Yeast Extract Peptone Dextrose (YPD) were listed and linked to exchange reactions in both the v7 and the BIGG model. Exchange reactions are defined to describe the metabolite composition in the immediate surroundings of the cell. The upper and lower bounds for all other exchange reactions were then set to ‘0’.
Once the lux reactions were inserted into both models, we analyzed the metabolic burden on the cell using a two-step process. First, we ran the built-in COBRA function robustnessAnalysis. This function plots the flux through a pathway against the flux through the optimized pathway [1]. In this case, the flux through the light pathway was negatively correlated with the flux through the growth pathway. We concluded that this result could only have been obtained if the light pathway was redirecting a metabolite that would normally have lent itself to growth.
Lorem ipsum dolor sit amet, consectetur adipiscing elit. Mauris eu suscipit lectus. Pellentesque ut risus rhoncus, congue tortor at, aliquam augue. Vestibulum ex mi, varius quis sollicitudin at, blandit ac lorem. Vivamus mattis sapien turpis, in fringilla nulla cursus id. Vestibulum vestibulum velit et accumsan aliquet. Aenean nulla justo, scelerisque id pulvinar eu, fringilla et nisl. Cras sapien magna, tincidunt in sapien et, sagittis sodales lacus. Praesent a ex ut augue fringilla interdum.
Lorem ipsum dolor sit amet, consectetur adipiscing elit. Mauris eu suscipit lectus. Pellentesque ut risus rhoncus, congue tortor at, aliquam augue. Vestibulum ex mi, varius quis sollicitudin at, blandit ac lorem. Vivamus mattis sapien turpis, in fringilla nulla cursus id. Vestibulum vestibulum velit et accumsan aliquet. Aenean nulla justo, scelerisque id pulvinar eu, fringilla et nisl. Cras sapien magna, tincidunt in sapien et, sagittis sodales lacus. Praesent a ex ut augue fringilla interdum.
Lorem ipsum dolor sit amet, consectetur adipiscing elit. Mauris eu suscipit lectus. Pellentesque ut risus rhoncus, congue tortor at, aliquam augue. Vestibulum ex mi, varius quis sollicitudin at, blandit ac lorem. Vivamus mattis sapien turpis, in fringilla nulla cursus id. Vestibulum vestibulum velit et accumsan aliquet. Aenean nulla justo, scelerisque id pulvinar eu, fringilla et nisl. Cras sapien magna, tincidunt in sapien et, sagittis sodales lacus. Praesent a ex ut augue fringilla interdum.
Future Directions
Lorem ipsum dolor sit amet, consectetur adipiscing elit. Mauris eu suscipit lectus. Pellentesque ut risus rhoncus, congue tortor at, aliquam augue. Vestibulum ex mi, varius quis sollicitudin at, blandit ac lorem. Vivamus mattis sapien turpis, in fringilla nulla cursus id. Vestibulum vestibulum velit et accumsan aliquet. Aenean nulla justo, scelerisque id pulvinar eu, fringilla et nisl. Cras sapien magna, tincidunt in sapien et, sagittis sodales lacus. Praesent a ex ut augue fringilla interdum.
Lorem ipsum dolor sit amet, consectetur adipiscing elit. Mauris eu suscipit lectus. Pellentesque ut risus rhoncus, congue tortor at, aliquam augue. Vestibulum ex mi, varius quis sollicitudin at, blandit ac lorem. Vivamus mattis sapien turpis, in fringilla nulla cursus id. Vestibulum vestibulum velit et accumsan aliquet. Aenean nulla justo, scelerisque id pulvinar eu, fringilla et nisl. Cras sapien magna, tincidunt in sapien et, sagittis sodales lacus. Praesent a ex ut augue fringilla interdum.
Scripts
Link to Constraint-Based Modeling Scripts
References
REFERENCES