Towards a kinetic model of the Entner-Doudoroff pathway in Zymomonas mobilis

Metabolic networks of cellular systems are complex, in that there are numerous components with multiple non-linear interactions. To understand how these networks work they are often split into manageable pieces and studied individually. However, an individual part is unable to account for the complex properties of systems. In order to study these interactions the field of systems biology has developed. Systems biology makes use of computers to construct models as a method to describe aspects of living systems. Using cellular pathways, kinetic models of metabolic pathways can be constructed and used as a tool to study the biological systems and provide a quantitative description. This thesis describes the quantitative analysis of a bacterium using a systems biology approach. Zymomonas mobilis is a rod shaped, Gram-negative, non-mobile facultative anaerobe and has one of the fastest observed fermentations, yet least energy efficient extractions found in nature. Furthermore it is the only known micro-organism to use the Entner-Doudoroff (2-keto-3-deoxy-6phosphogluconate) pathway anaerobically. The low energy yield of fermentation in Z. mobilis is a result of the usage of the Entner-Doudoroff glycolytic pathway, which has half the energy yield per mol substrate compared to the well known Embden-Meyerhof-Parnas glycolytic pathway. The work presented in this thesis forms part of a larger project to compare glycolytic regulation in different micro-organisms Z. mobilis, Escherichia coli, Saccharomyces cerevisiae and Lactococcus lactis. These organisms were chosen based on their usage of different glycolytic mechanisms. Kinetic models are suitable tools to draw a comparison between these organisms. The emphasis here is on the construction of a kinetic model of the Entner-Doudoroff glycolytic pathway as it occurs in Z. mobilis. The aim of this thesis was to characterise as many of the Entner-Doudoroff pathway enzymes as possible, under standard conditions. This was done using enzyme assays, to obtain the kinetic parameters of each of the enzymes. Microtitre plate assays were used to characterise most of the enzymes of the Entner-Doudoroff pathway. However, not all characterisations could be done using plate assay methods, as some intermediates were not commercially available to perform coupled assays. Nuclear magnetic resonance (NMR) spectroscopy was used to characterise these enzymes. These experimentally obtained parameters were then incorporated in a mathematical framework. Time simulations on the initial model were unable to reach a steady-state, with a build up of metabolic intermediates. A secondary model was constructed (using calculated maximal activities) which allowed us to identify discrepancies in the initial model. This showed that the experimentally determined maximal activities of three enzymes in lower glycolysis were unrealistically low, which might be due to protein denaturation by sonication. A final model was constructed which incorporated a correction factor for these three enzymes. The models’ predicted output (steady-state concentrations and flux) was compared to that of ii Stellenbosch University http://scholar.sun.ac.za