BackgroundA well known example of oscillatory phenomena is the transient oscillations of glycolytic intermediates in Saccharomyces cerevisiae, their regulation being predominantly investigated by mathematical modeling. To our knowledge there has not been a genetic approach to elucidate the regulatory role of the different enzymes of the glycolytic pathway.ResultsWe report that the laboratory strain BY4743 could also be used to investigate this oscillatory phenomenon, which traditionally has been studied using S. cerevisiae X2180. This has enabled us to employ existing isogenic deletion mutants and dissect the roles of isoforms, or subunits of key glycolytic enzymes in glycolytic oscillations. We demonstrate that deletion of TDH3 but not TDH2 and TDH1 (encoding glyceraldehyde-3-phosphate dehydrogenase: GAPDH) abolishes NADH oscillations. While deletion of each of the hexokinase (HK) encoding genes (HXK1 and HXK2) leads to oscillations that are longer lasting with lower amplitude, the effect of HXK2 deletion on the duration of the oscillations is stronger than that of HXK1. Most importantly our results show that the presence of beta (Pfk2) but not that of alpha subunits (Pfk1) of the hetero-octameric enzyme phosphofructokinase (PFK) is necessary to achieve these oscillations. Furthermore, we report that the cAMP-mediated PKA pathway (via some of its components responsible for feedback down-regulation) modulates the activity of glycoytic enzymes thus affecting oscillations. Deletion of both PDE2 (encoding a high affinity cAMP-phosphodiesterase) and IRA2 (encoding a GTPase activating protein- Ras-GAP, responsible for inactivating Ras-GTP) abolished glycolytic oscillations.ConclusionsThe genetic approach to characterising the glycolytic oscillations in yeast has demonstrated differential roles of the two types of subunits of PFK, and the isoforms of GAPDH and HK. Furthermore, it has shown that PDE2 and IRA2, encoding components of the cAMP pathway responsible for negative feedback regulation of PKA, are required for glycolytic oscillations, suggesting an enticing link between these cAMP pathway components and the glycolysis pathway enzymes shown to have the greatest role in glycolytic oscillation. This study suggests that a systematic genetic approach combined with mathematical modelling can advance the study of oscillatory phenomena.

Supplementary data are available at Bioinformatics online.

The study of biological systems at the genome scale helps us understand fundamental biological processes that govern the activity of living organisms and regulate their interactions with the environment. Genome-scale metabolic models are usually analysed using constraint-based methods, since detailed rate equations and kinetic parameters are often missing. However, constraint-based analysis is limited in capturing the dynamics of cellular processes. In this paper, we present an approach to build a genome-scale kinetic model of Mycobacterium tuberculosis metabolism using generic rate equations. M. tuberculosis causes tuberculosis which remains one of the largest killer infectious diseases. Using a genetic algorithm, we estimated kinetic parameters for a genome-scale metabolic model of M. tuberculosis based on flux distributions derived from Flux Balance Analysis. Our results show that an excellent agreement with flux values is obtained under several growth conditions, although kinetic parameters may vary in different conditions. Parameter variability analysis indicates that a high degree of redundancy remains present in model parameters, which suggests that the integration of other types of high-throughput datasets will enable the development of better constrained models accounting for a variety of in vivo phenotypes.