From steady-state to synchronized yeast glycolytic oscillations I: model construction

An existing detailed kinetic model for the steady-state behavior of yeast glycolysis was tested for its ability to simulate dynamic behavior. Using a small subset of experimental data, the original model was adapted by adjusting its parameter values in three optimization steps. Only small adaptations to the original model were required for realistic simulation of experimental data for limit-cycle oscillations. The greatest changes were required for parameter values for the phosphofructokinase reaction. The importance of ATP for the oscillatory mechanism and NAD(H) for inter-and intra-cellular communications and synchronization was evident in the optimization steps and simulation experiments. In an accompanying paper [du Preez F et al. (2012) FEBS J279, 2823-2836], we validate the model for a wide variety of experiments on oscillatory yeast cells. The results are important for re-use of detailed kinetic models in modular modeling approaches and for approaches such as that used in the Silicon Cell initiative.

The mathematical models described here have been submitted to the JWS Online Cellular Systems Modelling Database and can be accessed at http://jjj.biochem.sun.ac.za/database/dupreez/index.html.