Difference between revisions of "Team:Korea U Seoul/Project/Biological Background/content"

Line 51: Line 51:
 
             <div class="content_">
 
             <div class="content_">
 
             <p class="ti2">
 
             <p class="ti2">
    The program ‘Gil’ is a project based on KEGG (Kyoto Encyclopedia of Genesand Genomes) database (Kanehisa M., et al. 2014; and Kanehisa M., Goto S., 2000). Its information were accumulated by genome sequencing and other biological research. It offers well-combined metabolic-related pathways.  
+
    The program ‘Gil’ is a project based on KEGG (Kyoto Encyclopedia of Genes and Genomes) database (KanehisaM., et al. 2014; and KanehisaM., GotoS., 2000). There are numerous advantages of KEGG (Adriaens M. E., et al. n.d.). Above all, KEGG has a nice interface. In an overview graphic, we can take all reactions of a pathway, and see the enzymes that are involved in the reaction by clicking individual reactions. In addition, there are most information about the gene linked with this enzyme. So we organized and parsed the databases that provide information of compounds, reactions, genes and enzymes.  
Also we organized and parsed these database that provide compound, reaction, gene and enzymatic information. Especially, we built RPair, a main network in four reaction pair groups. However, there are some shortcomings in KEGG database.</p><p class="ti2"> First, it doesn't provide thermodynamic information that is a crucial criterion to determine reaction direction. Though we cannot obtain thermodynamic information in KEGG database, we tried to compensate it by utilizing Equilibrater and anotherdatabase. Additionally KEGG database, however, show reactions in bidirectional way in which usually do not process as it is a thermodynamically unfavorable reaction. Lastly, every glycan in KEGG database is not connected with compound. By origin, glycans have G number (ex. G00000) but there are some cases that glycans have C number (ex. C00000) that kegg compounds have. Specially, they can react to compounds or glycans. However, many glycans don’t have C number, so they can not react with compounds and they even don’t have any reactions in glycans. For example, glycan named neoagarobiose in the path that has been being studying for experimental validation doesn’t have C number. So we added the information of the glycan found other databases such as “-------”* and -------** to data obtained from kegg compounds to complement kegg’s fault.
+
<p>
+
*  a database of chemical molecules maintained by the National Center for Biotechnology information (NCBI) (Wikipedia)
+
<br>** Eun Ju Yun et al. (2015), The novel catabolic pathway of 3,6-anhydro-Lgalactose,
+
the main component of red macroalgae, in a marine bacterium, Environmental Microbiology, 1677 – 1688p
+
</p>
+
 
+
We used K12 as our E. coli model strain. It is the most commonly used E. coli strain in biological lab.  
+
 
+
 
+
 
             </p>
 
             </p>
 +
<p class="ti2">However, KEGG database is definitely not perfect. First, it doesn't provide thermodynamic data. Additionally, KEGG database shows reactions in a bidirectional way, which lacks the information of thermodynamical feasibility. Also, KEGG is a pure gene database. While each enzyme has a large amount of gene information attached, it lacks a bit in the other departments. Therefore, users easily encounter empty pages that do not have any information. Lastly, every glycan in KEGG database is not connected with its corresponding compound. When glycan IDs are not connected to compound IDs, it is not possible to find a reaction with a glycan.</p>
 
</div>
 
</div>
 
     <div class="title">
 
     <div class="title">
Line 72: Line 63:
 
             <div class="content_">
 
             <div class="content_">
 
             <p class="ti2">
 
             <p class="ti2">
  Thermodynamic constrain is one of the most important factors of biological experiments. Therefore, the program ‘Gil’ provides Gibbs energy of each metabolic reaction.
+
  Thermodynamic constraint is one of the most important factors of biological experiments. Therefore, the program ‘Gil’ provides Gibbs energy changes of each metabolic reaction.  
The program ‘Gil’ obtained thermodynamic datafrom two separated databases even though they were obtained from the same program, eQuilibrator (Flamholz A., et al. 2012). The first database is composed of Gibbs energy from standard condition—pH 7.0 and 0.1 M ionic strength. These data are correct, but the size of the database is smaller than the other one due to the lack of wetlab experiments. On the other hand, the second version is made up of values from calculation. This database has much more values than the first one, however it is hard for researchers to apply this values directly to their experiments since they are not experimentally proved. Thus, additional attention is needed to utilize the second version of data.
+
</p>
 
+
<p class="ti2">The program ‘Gil’ obtained thermodynamic data from two separated databases even though they were obtained from the same program, eQuilibrator (Flamholz A., et al. 2012). The first database is composed of Gibbs free energy from standard condition—pH 7.0 and 0.1 M ionic strength. These data are correct, but the size of the database is smaller than the other one due to the lack of wet lab experiments. On the other hand, the second version is made up of values from a calculation. This database has much more values than the first one however, it is hard for researchers to apply this values directly to their experiments since they are not experimentally proved. Thus, additional attention is needed to utilize the second version of data. </p>
            </p>
+
</div>
</div>
+
<p class="ti2">
+
Thermodynamic constrain is one of the most important factors of biological experiments. Therefore, the program ‘Gil’ provides Gibbs energy of each metabolic reaction.
+
</p><p class="ti2">The program ‘Gil’ obtained thermodynamic datafrom two separated databases even though they were obtained from the same program, eQuilibrator (Flamholz A., et al. 2012). The first database is composed of Gibbs energy from standard condition—pH 7.0 and 0.1 M ionic strength. These data are correct, but the size of the database is smaller than the other one due to the lack of wetlab experiments. On the other hand, the second version is made up of values from calculation. This database has much more values than the first one, however it is hard for researchers to apply this values directly to their experiments since they are not experimentally proved. Thus, additional attention is needed to utilize the second version of data.  
+
</p>
+
  
  
Line 98: Line 84:
 
<img src="https://static.igem.org/mediawiki/2015/4/43/KoreaUSeoul_BiologicalBackground_Figure1.jpg">
 
<img src="https://static.igem.org/mediawiki/2015/4/43/KoreaUSeoul_BiologicalBackground_Figure1.jpg">
  
  <p class="ti2"> The NAD + and NADH molecules are indistinguishable in structure to NADP + and NADPH, respectively, except the fact that the phosphate group is absent from the former. However, this smalldifference is crucial for their distinctive roles. Although the extra phosphate group on NADPH is located far away from the region involved in electron transfer (see Figure 1-b), it gives the molecule a slight difference in its shape from that of NADH, so NADPH and NADH bind to different sets of enzymes as substrates. Thus, thetwo types of carriers are used to transfer electrons (or hydride ions) between different sets of molecules. This role separationis crucial for organisms to regulate two sets of electron-transfer reactions independently. A NADPH operates mainly with enzymes that catalyze anabolic reactions, providing the high-energy electrons needed to synthesize energy-rich biological molecules. In contrast, a NADH has a special role as an intermediate in the catabolic system of reactions that generate ATP through the oxidation of compounds (Alberts, B., et al. 2002).
+
  <p class="ti2"> The NAD + and NADH molecules are indistinguishable in structure to NADP + and NADPH, respectively, except the fact that the phosphate group is absent from the former. However, this small difference is crucial for their distinctive roles. Although the extra phosphate group on NADPH is located far away from the region involved in electron transfer (see Figure 1-b), it gives the molecule a slight difference in its shape from that of NADH, so NADPH and NADH bind to different sets of enzymes as substrates. Thus, the two types of carriers are used to transfer electrons (or hydride ions) between different sets of molecules. This role separation is crucial for organisms to regulate two sets of electron-transfer reactions independently. An NADPH operates mainly with enzymes that catalyze anabolic reactions, providing the high-energy electrons needed to synthesize energy-rich biological molecules. In contrast, an NADH has a special role as an intermediate in the catabolic system of reactions that generate ATP through the oxidation of compounds (Alberts B., et al. 2002).
 
</p>
 
</p>
 
<p class="ti2">Then, why should we focus on the larger amount of NADH and NADPH as products? </p>
 
<p class="ti2">Then, why should we focus on the larger amount of NADH and NADPH as products? </p>
<p class="ti2">First, in the case of NADPH, it is the key for determining the efficiency of certain biosynthetic pathways. The oxidation during a pentose phosphorylation produce most NADPH which is needed in the biosynthesis of certain metabolites. Therefore, thecarbon flow of the pentose phosphorylation path is considered to be the targets of metabolic engineering. For example, microorganisms found in the Corynebacterium genus, having enhanced L-ornithine production, not only activate pentose phosphate pathway but also increase the specific activities of G6PDH and 6PGD, which are enzymes related to the pentose phosphate pathway, by reducing or inactivating intrinsic GDH (glucose dehydrogenase) (see Figure 2). Consequently, the improved NADPH production makes L-ornithine yield to increase, which may have the possibility of being widely utilized for efficient and economical applications (Cho J., Ko Y. 2012). To summarize, if we find a path with improved NADPH yield, we could get more products in other relevant paths.
+
<p class="ti2">First, in the case of NADPH, it is the key for determining the efficiency of certain biosynthetic pathways. The oxidation during a pentose phosphorylation produces most NADPH which is needed in the biosynthesis of certain metabolites. Therefore, the carbon flow of the pentose phosphorylation path is considered to be the targets of metabolic engineering. For example, microorganisms found in the Corynebacterium genus, having enhanced L-ornithine production, not only activate pentose phosphate pathway but also increase the specific activities of G6PDH and 6PGD, which are enzymes related to the pentose phosphate pathway, by reducing or inactivating intrinsic GDH (glucose dehydrogenase) (see Figure 2). Consequently, the improved NADPH production makes L-ornithine yield to increase, which may have the possibility of being widely utilized for efficient and economical applications (Cho J., Ko Y. 2012). To summarize, if we find a path with improved NADPH yield, we could get more products in other relevant paths.
 
</p>
 
</p>
 
<img src="https://static.igem.org/mediawiki/2015/f/f6/KoreaUSeoul_BiologicalBackground_Figure2.png" width="700" height="800">
 
<img src="https://static.igem.org/mediawiki/2015/f/f6/KoreaUSeoul_BiologicalBackground_Figure2.png" width="700" height="800">
 
<p>Figure 2</p>
 
<p>Figure 2</p>
<p class="ti2">In the same context, if theNADH yield is increased, the efficiency of the path using the NADH will also be enhanced. In addition, cells grown aerobically may possibly use the extra NADH for ATP generation through the electron transport system. According to the research on the availability of intracellular NADH (Berríos-Rivera SJ., Bennett GN, San KY 2002), the higher NADH availability increased dramatically the ethanol to acetate (Et/Ac) ratio by 27-fold. The cause of this amazing increase is proven by an experiment comparing BS1 (pSBF2), which is inactivated by the native FDH (formate dehydrogenase) that regenerates NADH, with GJT001 (pDHK29), the control, under anaerobic conditions.
+
<p class="ti2">In the same context, if the NADH yield is increased, the efficiency of the path using the NADH will also be enhanced. In addition, cells grown aerobically may possibly use the extra NADH for ATP generation through the electron transport system. According to the research on the availability of intracellular NADH (Berríos-Rivera SJ., Bennett GN, San KY 2002), the higher NADH availability dramatically increased the ethanol to acetate (Et/Ac) ratio by 27-fold. The cause of this amazing increase is proven by an experiment comparing BS1 (pSBF2), which is inactivated by the native FDH (formate dehydrogenase) that regenerates NADH, with GJT001 (pDHK29), the control, under anaerobic conditions.
 
</p>
 
</p>
 
<p class="ti2">For these reasons, our program suggests the paths in the order of the amount of NADPH and NADH production. Users will be able to get efficient and productive paths with the information provided by our program, with less effort.
 
<p class="ti2">For these reasons, our program suggests the paths in the order of the amount of NADPH and NADH production. Users will be able to get efficient and productive paths with the information provided by our program, with less effort.
Line 111: Line 97:
  
 
<h2>② ATP</h2>
 
<h2>② ATP</h2>
<p class="ti2">Another scoring factor which the program ‘Gil’ provides is the ATP quantity. This function is expected to be useful for synthetic biologists to design reaction pathway necessary for their experiments which is designed to have higher efficiency when there is incread ATP level. For instance, according to the research that used Saccharomyces cerevisiae to generate ATP from adenosinein the process of glucose breakdown into ethanol and carbon dioxide (Murata K, et al.. 1981), theATP-generating reaction enhances the biosynthesis of Escherichia coli which synthesizes glutathione and Brevibacterium ammoniagenes which produces NADP.</p><p class="ti2"> Therefore, substituting the biological reaction with the ones that have higher ATP yield, designed by the program ‘Gil’, will possibly increase the productivity of the both organisms. The following KEGG compound number order is the examples of the reaction pathway: C00031 (glucose), C00394, C00035, C00144, C00044, C04895, C00084, C00469 (ethanol).
+
<p class="ti2">Another scoring factor which the program ‘Gil’ provides is the ATP quantity. This function is expected to be useful for synthetic biologists to design reaction pathway necessary for their experiments which are designed to have higher efficiency when there is increased ATP level. For instance, according to the research that used Saccharomyces cerevisiae to generate ATP from adenosine in the process of glucose breakdown into ethanol and carbon dioxide (Murata K, et al.. 1981), the ATP-generating reaction enhances the biosynthesis of Escherichia coli which synthesizes glutathione and Brevibacterium ammoniagenes which produces NADP. Therefore, substituting the biological reaction with the ones that have higher ATP yield, designed by the program ‘Gil’, will possibly increase the productivity of both organisms. The following KEGG compound number order is the examples of the reaction pathway: C00031 (glucose), C00394, C00035, C00144, C00044, C04895, C00084, C00469 (ethanol).
 
</p>
 
</p>
  
 
<h2>③ CO2</h2>
 
<h2>③ CO2</h2>
<p class="ti2">Most chemical compounds with carbon produce a product through a synthetic-dissociation pathway. The decarboxylation, thedissociation of carbon, will lose carbon equivalent and limit the product’s theoretical carbon yield. The Formula 1 show you how to calculate the efficiency of a pathway.</p>
+
<p class="ti2">Most chemical compounds with carbon produce a product through a synthetic-dissociation pathway. The decarboxylation, the dissociation of carbon, will lose carbon equivalent and limit the product’s theoretical carbon yield. The Formula 1 show you how to calculate the efficiency of a pathway.</p>
 
<img src="https://static.igem.org/mediawiki/2015/1/14/KoreaUSeoul_BiologicalBackground_Formula1.gif">
 
<img src="https://static.igem.org/mediawiki/2015/1/14/KoreaUSeoul_BiologicalBackground_Formula1.gif">
<p class="ti2">We assume that non-oxidative glycolysis (NOG) as model mechanism to develop our software program. NOG is one of the complete carbon yield pathwaysthat starts with sugar to acetyl-CoA. It needs conjunction with carbon dioxide or assimilation with other one-carbon (C1) source to produce desirable fuels and chemicals (Bogorad I. W., Lin T. S., Liao J. C. 2013). </p>
+
<p class="ti2">We assume that non-oxidative glycolysis (NOG) as a model mechanism to develop our software program. NOG is one of the complete carbon yield pathways that starts with sugar to acetyl-CoA. It needs conjunction with carbon dioxide or assimilation with other one-carbon (C1) source to produce desirable fuels and chemicals (Bogorad I. W., Lin T. S., Liao J. C. 2013). </p>
 
<p class="ti2">While simulating the NOG pathway, we achieved complete carbon yield. So without losing any of carbons, we are able to obtain desirable results. During the reaction, there are lots of chemicals that result in carbon loss. The program ‘Gil’ calculates the CO2 reduction to calculate the carbon loss. This is due to the fact that CO2 is less likely to be utilized unless it is fixed compared to other chemicals.
 
<p class="ti2">While simulating the NOG pathway, we achieved complete carbon yield. So without losing any of carbons, we are able to obtain desirable results. During the reaction, there are lots of chemicals that result in carbon loss. The program ‘Gil’ calculates the CO2 reduction to calculate the carbon loss. This is due to the fact that CO2 is less likely to be utilized unless it is fixed compared to other chemicals.
 
</p>
 
</p>
Line 137: Line 123:
 
<p>BogoradI. W., LinT. S., LiaoJ. C. “Synthetic non-oxidative glycolysis enables complete carbon conservation.” “Nature”, 2013: 693-697.
 
<p>BogoradI. W., LinT. S., LiaoJ. C. “Synthetic non-oxidative glycolysis enables complete carbon conservation.” “Nature”, 2013: 693-697.
 
</p>
 
</p>
<p>ChoJ., KoY. Microorganism found in the corynebacterium genus having enhanced l-ornithine production and method for preparing l-ornithine using same. South Korea 특허권 WO 2012161522 A2. 2012년 5월 23일.
+
<p>Cho, J., and Y. Ko. Microorganism found in the corynebacterium genus having enhanced l-ornithine production and method for preparing l-ornithine using same. South Korea Patent WO 2012161522 A2. 5 23, 2012.
 
</p>
 
</p>
 
<p>FlamholzA., NoorE., Bar-EvenA., MiloR. “eQuilibrator—the biochemical thermodynamics calculator.” “Nucleic Acids Res.”, 2012: D770-D775.
 
<p>FlamholzA., NoorE., Bar-EvenA., MiloR. “eQuilibrator—the biochemical thermodynamics calculator.” “Nucleic Acids Res.”, 2012: D770-D775.

Revision as of 23:21, 18 September 2015

1) KEGG database

The program ‘Gil’ is a project based on KEGG (Kyoto Encyclopedia of Genes and Genomes) database (KanehisaM., et al. 2014; and KanehisaM., GotoS., 2000). There are numerous advantages of KEGG (Adriaens M. E., et al. n.d.). Above all, KEGG has a nice interface. In an overview graphic, we can take all reactions of a pathway, and see the enzymes that are involved in the reaction by clicking individual reactions. In addition, there are most information about the gene linked with this enzyme. So we organized and parsed the databases that provide information of compounds, reactions, genes and enzymes.

However, KEGG database is definitely not perfect. First, it doesn't provide thermodynamic data. Additionally, KEGG database shows reactions in a bidirectional way, which lacks the information of thermodynamical feasibility. Also, KEGG is a pure gene database. While each enzyme has a large amount of gene information attached, it lacks a bit in the other departments. Therefore, users easily encounter empty pages that do not have any information. Lastly, every glycan in KEGG database is not connected with its corresponding compound. When glycan IDs are not connected to compound IDs, it is not possible to find a reaction with a glycan.

2) Thermodynamic information

Thermodynamic constraint is one of the most important factors of biological experiments. Therefore, the program ‘Gil’ provides Gibbs energy changes of each metabolic reaction.

The program ‘Gil’ obtained thermodynamic data from two separated databases even though they were obtained from the same program, eQuilibrator (Flamholz A., et al. 2012). The first database is composed of Gibbs free energy from standard condition—pH 7.0 and 0.1 M ionic strength. These data are correct, but the size of the database is smaller than the other one due to the lack of wet lab experiments. On the other hand, the second version is made up of values from a calculation. This database has much more values than the first one however, it is hard for researchers to apply this values directly to their experiments since they are not experimentally proved. Thus, additional attention is needed to utilize the second version of data.

3) Four path-scoring criteria

① NADH and NADPH

Like an ATP (adenosine triphosphate) which serves as a convenient and versatile storage in cells, NADH (nicotinamide adenine dinucleotide) and NADPH (nicotinamide adenine dinucleotide phosphate) serve as activated carriers specialized in carrying high-energy electrons and hydrogen atoms (see Figure 1).

The NAD + and NADH molecules are indistinguishable in structure to NADP + and NADPH, respectively, except the fact that the phosphate group is absent from the former. However, this small difference is crucial for their distinctive roles. Although the extra phosphate group on NADPH is located far away from the region involved in electron transfer (see Figure 1-b), it gives the molecule a slight difference in its shape from that of NADH, so NADPH and NADH bind to different sets of enzymes as substrates. Thus, the two types of carriers are used to transfer electrons (or hydride ions) between different sets of molecules. This role separation is crucial for organisms to regulate two sets of electron-transfer reactions independently. An NADPH operates mainly with enzymes that catalyze anabolic reactions, providing the high-energy electrons needed to synthesize energy-rich biological molecules. In contrast, an NADH has a special role as an intermediate in the catabolic system of reactions that generate ATP through the oxidation of compounds (Alberts B., et al. 2002).

Then, why should we focus on the larger amount of NADH and NADPH as products?

First, in the case of NADPH, it is the key for determining the efficiency of certain biosynthetic pathways. The oxidation during a pentose phosphorylation produces most NADPH which is needed in the biosynthesis of certain metabolites. Therefore, the carbon flow of the pentose phosphorylation path is considered to be the targets of metabolic engineering. For example, microorganisms found in the Corynebacterium genus, having enhanced L-ornithine production, not only activate pentose phosphate pathway but also increase the specific activities of G6PDH and 6PGD, which are enzymes related to the pentose phosphate pathway, by reducing or inactivating intrinsic GDH (glucose dehydrogenase) (see Figure 2). Consequently, the improved NADPH production makes L-ornithine yield to increase, which may have the possibility of being widely utilized for efficient and economical applications (Cho J., Ko Y. 2012). To summarize, if we find a path with improved NADPH yield, we could get more products in other relevant paths.

Figure 2

In the same context, if the NADH yield is increased, the efficiency of the path using the NADH will also be enhanced. In addition, cells grown aerobically may possibly use the extra NADH for ATP generation through the electron transport system. According to the research on the availability of intracellular NADH (Berríos-Rivera SJ., Bennett GN, San KY 2002), the higher NADH availability dramatically increased the ethanol to acetate (Et/Ac) ratio by 27-fold. The cause of this amazing increase is proven by an experiment comparing BS1 (pSBF2), which is inactivated by the native FDH (formate dehydrogenase) that regenerates NADH, with GJT001 (pDHK29), the control, under anaerobic conditions.

For these reasons, our program suggests the paths in the order of the amount of NADPH and NADH production. Users will be able to get efficient and productive paths with the information provided by our program, with less effort.

② ATP

Another scoring factor which the program ‘Gil’ provides is the ATP quantity. This function is expected to be useful for synthetic biologists to design reaction pathway necessary for their experiments which are designed to have higher efficiency when there is increased ATP level. For instance, according to the research that used Saccharomyces cerevisiae to generate ATP from adenosine in the process of glucose breakdown into ethanol and carbon dioxide (Murata K, et al.. 1981), the ATP-generating reaction enhances the biosynthesis of Escherichia coli which synthesizes glutathione and Brevibacterium ammoniagenes which produces NADP. Therefore, substituting the biological reaction with the ones that have higher ATP yield, designed by the program ‘Gil’, will possibly increase the productivity of both organisms. The following KEGG compound number order is the examples of the reaction pathway: C00031 (glucose), C00394, C00035, C00144, C00044, C04895, C00084, C00469 (ethanol).

③ CO2

Most chemical compounds with carbon produce a product through a synthetic-dissociation pathway. The decarboxylation, the dissociation of carbon, will lose carbon equivalent and limit the product’s theoretical carbon yield. The Formula 1 show you how to calculate the efficiency of a pathway.

We assume that non-oxidative glycolysis (NOG) as a model mechanism to develop our software program. NOG is one of the complete carbon yield pathways that starts with sugar to acetyl-CoA. It needs conjunction with carbon dioxide or assimilation with other one-carbon (C1) source to produce desirable fuels and chemicals (Bogorad I. W., Lin T. S., Liao J. C. 2013).

While simulating the NOG pathway, we achieved complete carbon yield. So without losing any of carbons, we are able to obtain desirable results. During the reaction, there are lots of chemicals that result in carbon loss. The program ‘Gil’ calculates the CO2 reduction to calculate the carbon loss. This is due to the fact that CO2 is less likely to be utilized unless it is fixed compared to other chemicals.

References

Alberts, B., A. Johnson, J. Lewis, and et al. Molecular Biology of the Cell. 4th edition. New York: Garland Science, 2002.

Berríos-RiveraSJ., BennettGN, SanKY. “Metabolic Engineering of Escherichia coli : Increase of NADH Availability by Overexpressing an NAD+-Dependent Formate Dehydrogenase.” “Metabolic Engineering”, 2002: 217-229.

BogoradI. W., LinT. S., LiaoJ. C. “Synthetic non-oxidative glycolysis enables complete carbon conservation.” “Nature”, 2013: 693-697.

Cho, J., and Y. Ko. Microorganism found in the corynebacterium genus having enhanced l-ornithine production and method for preparing l-ornithine using same. South Korea Patent WO 2012161522 A2. 5 23, 2012.

FlamholzA., NoorE., Bar-EvenA., MiloR. “eQuilibrator—the biochemical thermodynamics calculator.” “Nucleic Acids Res.”, 2012: D770-D775.

KanehisaM., GotoS. “KEGG: Kyoto Encyclopedia of Genes and Genomes.” “Nucleic Acids Res”, 2000: 27-30.

KanehisaM., GotoS., SatoY., KawashimaM., FurumichiM., TanabeM. “Data, information, knowledge and principle: back to metabolism in KEGG.” “Nucleic Acids Res”, 2014: D199-D205.

MurataK, TaniK., KatoJ., ChibataI. “Glycolytic pathway as an ATP generation system and its application to the production of glutathione and NADP.” “Enzyme and Microbial Technology”, 1981: 233-242.