Can yeast glycolysis be understood in terms of in vitro kinetics of the constituent enzymes? Testing biochemistry
Catherine
Lloyd
Auckland Bioengineering Institute, The University of Auckland
Model Status
This CellML runs in both OpenCell and COR to recreate the published results. The units have been checked and they are consistent. This particular CellML model is based on a Matlab file which, in turn, was exported from the working SBML file in the BioModels database. There are several subtle differences between the published paper and the working Matlab code which have been previously documented by the BioModels curators. In brief these are:
in the kinetic law for ADH the species a should denote NAD and b Ethanol the last term in the equation should read bpq/(Kib Kiq Kp) and the rate equation for the PFK contains two errors: the R term should read: R = 1 + lambda1 + lambda2 + gR lambda1 lambda2 and the last term in the L equation should be squared.
Also there are slight differences in the parameterisation of the model depending on whether it comes from the JWS or the BioModels database and, for Vmax values, please note that there is a conversion factor of approxinately 270 to convert from U/mg-protein as shown in Table 1 of the paper to mmol/(min*L_cytosol).
Model Structure
ABSTRACT: This paper examines whether the in vivo behavior of yeast glycolysis can be understood in terms of the in vitro kinetic properties of the constituent enzymes. In nongrowing, anaerobic, compressed Saccharomyces cerevisiae the values of the kinetic parameters of most glycolytic enzymes were determined. For the other enzymes appropriate literature values were collected. By inserting these values into a kinetic model for glycolysis, fluxes and metabolites were calculated. Under the same conditions fluxes and metabolite levels were measured. In our first model, branch reactions were ignored. This model failed to reach the stable steady state that was observed in the experimental flux measurements. Introduction of branches towards trehalose, glycogen, glycerol and succinate did allow such a steady state. The predictions of this branched model were compared with the empirical behavior. Half of the enzymes matched their predicted flux in vivo within a factor of 2. For the other enzymes it was calculated what deviation between in vivo and in vitro kinetic characteristics could explain the discrepancy between in vitro rate and in vivo flux.
The original paper reference is cited below:
Can yeast glycolysis be understood in terms of in vitro kinetics of the constituent enzymes? Testing biochemistry, Bas Teusink, Jutta Passarge, Corinne A. Reijenga, Eugenia Esgalhado, Coen C. van der Weijden, Mike Schepper, Michael C. Walsh, Barbara M. Bakker, Karel van Dam, Hans V. Westerhoff, and Jacky L. Snoep, 2000, European Journal of Biochemistry, 267, 5313-5329. PubMed ID: 10951190
reaction diagram
Schematic diagram of the glycolysis pathway model described by Teusink et al. 2000. Reactions in boxes show the branches introduced in the extended model.
GLCo
extracellular glucose
GLCi
intracellular glucose
G6P
glucose 6-phosphate
F6P
fructose 6-phosphate
F16P
fructose 1,6-bisphosphate
F26P
fructose 2,6-bisphosphate
TRIO
triosphosphate
BPG
1,3-bisphosphoglycerate
P3G
3-phosphoglycerate
P2G
2-phosphoglycerate
PEP
phosphoenolpyruvate
PYR
pyruvate
ACE
acetaldehyde
P
high energy phosphates
NADH
nicotinamide adenine dinucleotide
NAD
charged nicotinamide adenine dinucleotide
GAP
D-glyceraldehyde-3-phosphate
ADP
adenosine diphoshate
ATP
adenosine triphosphate
AMP
adenosine monophoshate
ETOH
ethanol
GLY
glycerol
This is the CellML description of Teusink et al.'s 2000 model for glycolysis in yeast.
Auckland Bioengineering Institute, The University of Auckland
Hans
Westerhoff
V
2000-09-01
Michael
Walsh
C
Catherine Lloyd
keyword
glycoloysis
yeast
metabolism
enzyme kinetics
Catherine Lloyd
Eugenia
Esgalhado
Jacky
Snoep
L
2008-01-07T07:39:20+13:00
European Journal of Biochemistry
10951190
2003-05-12T00:00:00+00:00
Coen
van der Weijden
C
This CellML version of the model has been recoded to remove the reaction element. The model has been opened in both PCEnv and COR and the units have been checked and are consistent. However, unfortunately the model can not be run in either PCEnv or COR at the moment due to the presence of 'unsuitable constraints' (PCEnv) or 'circular arguments' (COR). These coupled equations need to be solved iteratively at each time step. We intend this feature to become a future development of our simulation tools.
Barbara
Bakker
M
Teusink et al.'s 2000 model for glycolysis in yeast.
Saccharomyces cerevisae
Karel
van Dam
Catherine
Lloyd
May
Bas
Teusink
Corrine
Reijenga
A
Jutta
Passarge
Catherine
Lloyd
May
Recoded the model to remove the reaction element.
Mike
Schepper
Can yeast glycolysis be understood in terms of in vitro kinetics of the constituent enzymes? Testing biochemistry.
267
5313
5329
The University of Auckland
Auckland Bioengineering Institute
c.lloyd@auckland.ac.nz