"
Team:Alberta/Project/Modeling
From 2009.igem.org
 |
|
|
Modeling
Literature information is insufficient for making good predictions about what's required for a minimal E.coli metabolism
Many of the literature essential gene lists are either based on single gene deletions or the metabolic pathways of mycoplasma rather than E.coli. Due to pathway redundancy, many metabolic genes that would be essential in a minimal genome would not appear essential in single gene knock out experiments. As E.coli possesses structures such as the cell wall that are not present in mycoplasma, minimal mycoplasma genome propositions are not sufficient. To test which metabolic genes are essential for E.coli, we developed a computer model.The E.coli MG1655 genome was modeled by the Palsson group at the University of San Diego and was used for our modeling experiments. Additionally, the Cobra Toolbox developed by the System`s Biology Research Group was used to interface the model with the Matlab program. A series of multiple gene deletions were performed using the model in order to determine the essential metabolic genes. These were compared to the literature genes selected and based on the individual gene`s function and involvement in various metabolic pathways, the gene was either added to our master essential gene list or removed. Additionally, media conditions were altered for the cells environment allowing for predictions for the conditions which should be applied to the minimal cell once it is developed. Once completed our MatLab model contributed 116 genes to our master gene list. Metabolic modeling allows for computational analysis of entire genomes which would be impossible to accomplish any other way. The various sources and methods used to collect data has allowed for an unique gene list which has the best possible chance of producing a minimal genome. The MatLab protocols demonstrated in the modeling section can be used to identify any organism’s essential genes provided a model is available. |
Constraint Based Flux Analysis – Cobra Toolbox and SBMLConstraint Based Flux Analysis is a mathematical way of representing biological system information and allows for easy manipulation of this data. It is based on a stoichiometric matrix of reactions which correspond to individual enzymatic or transport reactions which have been characterized inside of the organism. In other words, it computationally represents each reaction using a linear array of numbers (see Figure 1). The flux can be defined as the amount of substrate moving to product for each individual reaction. The model assumes that the system is at steady state therefore, the overall flux is zero since each product becomes a substrate of another reaction. All substrates entering the system will have the same amount leaving the system. The products leaving the system can be removed by changing the boundary condition of the compound (that is making it unavailable to the system) or by using it to produce growth of the organism. A master growth equation determines which products are required for the cell to grow and this represented in units of DW/unit time. Systems Biology Markup Language (SBML) and the Cobra Toolbox (both produced from System’s Biology Research Group) allows for flux analysis to be performed in the MATLAB program. |
Finding Essential GenesWe can perform single gene knockouts in the e. coli metabolic gene model and determine the resulting growth rate. Knockouts that result in a growth rate of approximately zero are put in a list for the researcher to study later. Deletion of any one of these genes will result in no growth. These genes are clearly essential. The Cobra Toolbox makes this single gene deletion analysis for every gene very easy. We found 187 of the model's genes to be essential with this code. After enriching the in silico media with various substances (MediaChange.m), the code was run again, and a list of 112 genes resulted. Clearly, enriching the media relieves some genes of their function, perhaps by allowing the bypassing of a reaction catalyzed by a protein coded for by that gene. Comparing this list 112 genes to our previously complied list of researched literature essential genes, we found many differences which we found to be very helpful. Many of these missing genes were added to our master list of essential genes. Obviously, since the code only analyses single gene deletions, the e. coli cell would not survive with just these 112 genes (these 112 genes are not the minimal genome). Multiple gene deletions were not taken into account. For instance, consider 2 genes that at first appear unessential. When one is knocked out, the cell lives. When the other is knocked out, the cell lives. But lets say that once both in the pair of genes is deleted, the cell dies. The single gene deletion analysis would not detect this. Another code was written to find out how many pair of genes are essentially entangled in this way. We found that there are 81 of these pairs in the model (assuming glucose minimal media). To validate this model's claim of paired essentiality, we performed single and double gene knockouts in the lab for these essential pairs. See the Model Validation section for details. |
Finding the Model's Minimal GenomeIt was our goal to find the model's minimal genome. It is clear that single gene knockout analysis would not suffice. Since performing every possible combination of multiple gene knockouts in the model would take most computers many years to accomplish, a different strategy was taken. Before the code is explained, let it be known that there are a few factors that will affect the size of the minimal genome. The first factor is the order of gene deletions. For example, once you delete a particular gene, it may cause 3 other genes that were previously inessential to become essential. So depending on which gene was deleted first, the size of the resulting minimal genome after further deletions will change. Another factor to consider is the final growth rate. If you want a growth rate that is 90% what it was originally, your minimal genome may be quite large. If you want a final growth rate of 10% of what it was originally, then your minimal genome can afford to become much smaller. The growth medium will also affect the size of the minimal genome. For these reasons, we wrote a code that will accommodate these factors, and allow the researcher to find the minimal genome of any in silico model. For the curious, here is the code.
For something that is easier on the eyes, here is a simplified flow chart depicting an algorithm that part of the code goes through.
Throughout the processes described above, we were constantly crosschecking our previously compiled list of researched genes from the literature with these model genes. Each gene in the model's minimal genome were researched further to judge whether or not it truly needed to be included in a minimal genome. As a result, most of these model genes were consequently included in our master list of essential genes needed for the minimal e. coli genome. Modeling in MATLAB has proven to be an essential aspect of our project. And in the process, we have created tools for other researchers to use to help them find their organism's essential genes and minimal genome. |
Our Codes & How We Used ThemStart.m - this will set the solver to glpk (a way of solving the vast amounts of matrix equations), and load in the model that the researcher will be working with. FindEssentials.m - Uses single gene deletion analyses to find a list of essential genes. A deletion of just one of these genes will cause cell death. CheckForRepeats.m - Checks for repeats in a gene list. We use this to check our imported list (our researched literature essential gene list). InvertList.m - "Inverts" the input List1 with respect to the input List2. The output is a list of genes in List2 that did not appear in List1. We use this to find the genes in our imported literature list that we would be deleting from the model. This helps us analyze the effects of deleting genes that we previously thought to be inessential. FindMatches.m - Finds the matches between two lists of genes. We use this to compare our researched literature essential genes with our model's predicted essential genes. It outputs the number of matches and the list of matches. MediaChange.m - We could reduce the amount of genes that were essential by changing the in silico media. We added to the glucose minimal media: 20 Standard Amino Acids Adenosine, Guanosine, Uridine, Thymidine, Cytidine Thiamine Pantothenic Acid Xylose Shikimate NMN Meso Diaminopimelic Acid Biotin FindMinimalGenome.m - This is the code that finds the model's minimal genome(s) through iterative semi-random multiple gene deletions. Because of the fact that many genes code for the same reaction, many reactions perform the same function, and that one gene may control a few reactions, there are easily thousands of potential minimal genomes. This code finds many of them. We are not necessarily interested in the smallest minimal genome; we want to use the resulting lists to gain a better understanding of what genes should be included in an actual minimal genome that could be created in a real lab environment. You can specify the number of minimal genomes you want (the more you specify, the longer it will take). You can also specify the final growth rate wanted, and 2 other factors relating to how the algorithm is carried out which will affect the time it takes, and how small the minimal genomes will be. MaxMinAvg.m - This is a simple code which lets the user know the largest minimal genome, the smallest minimal genome, and the average size of the minimal genomes created in the FindMinimalGenome.m code. It can be used on other gene lists (for instance, on an array of gene lists that represent what genes are missing from the researched literature list with respect to each minimal genome found). InvertList.m - "Inverts" input1 with respect to input2. (Output is a list of genes that includes all genes in input2 that are not present in input1. InvertListArray.m - "Inverts" input1 with respect to every gene list in input2. input2 may include any number of gene lists. This code is typically used when comparing your researched literature gene list with each of the model's minimal genome lists (produced by FindMinimalGenome). FindMatchesArray.m - Typically only used with the array of minimal genomes and the researched literature essential genes. Outputs a list of matches for every minimal genome creating the output MatchListArray. The code also outputs the number of matches for each list and the position of the smallest minimal genome. This is very helpful when trying to find what genes you are missing in your researched literature essential genes list when comparing them with the many minimal genomes produced by our code. CellToMatrix.m - Code that makes it more convenient for you to copy and paste the array of lists into Excel. Without this step, you'd have to copy and paste the hundreds of unique minimal genome lists one at a time. |
