Model evaluation for glycolytic oscillations in yeast biotransformations of xenobiotics

On the role of lateral stabilization during early patterning in the pancreas

The cell fate decision of multi-potent pancreatic progenitor cells between the exocrine and endocrine lineages is regulated by Notch signalling, mediated by cell-cell interactions. However, canonical models of Notch-mediated lateral inhibition cannot explain the scattered spatial distribution of endocrine cells and the cell-type ratio in the developing pancreas. Based on evidence from acinar-to-islet cell transdifferentiation in vitro, we propose that lateral stabilization, i.e. positive feedback between adjacent progenitor cells, acts in parallel with lateral inhibition to regulate pattern formation in the pancreas. A simple mathematical model of transcriptional regulation and cell-cell interaction reveals the existence of multi-stability of spatial patterns whose simultaneous occurrence causes scattering of endocrine cells in the presence of noise. The scattering pattern allows for control of the endocrine-to-exocrine cell-type ratio by modulation of lateral stabilization strength. These theoretical results suggest a previously unrecognized role for lateral stabilization in lineage specification, spatial patterning and cell-type ratio control in organ development.

Exploring the genetic control of glycolytic oscillations in Saccharomyces cerevisiae

BACKGROUND

A 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.

RESULTS

We 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.

CONCLUSIONS

The 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.

Biosimulation of drug metabolism--a yeast based model

Computationally predicting the metabolic fates of drugs is a very complex task which is owed not only to the huge and diverse biochemical network in the living cell, but also to the majority of in vivo transformations that occur through the action of hepatocytes and gastro-intestinal micro-flora. Thus, xenobiotics are metabolised by more than a single cell type. However, the prediction of metabolic fates is definitely a problem worth solving since it would allow facilitate the development of drugs in a way less relying on animal testing. As a first step in this direction, PharmBiosim is being developed, a biosimulation tool which is based on substantial data reduction and on attributing metabolic fates of drug molecules to functional groups and substituents. This approach works with yeast as a model organism and is restricted to drugs that are mainly transformed by enzymes of the central metabolism, especially sugar metabolism. The reason for the latter is that the qualitative functioning of the involved biochemistry is very similar in diverse cell types involved in drug metabolism. Further it allows for using glycolytic oscillations as a tool to quantify interactions of a drug with this metabolic pathway.

Reduction of a biochemical model with preservation of its basic dynamic properties

The complexity of full-scale metabolic models is a major obstacle for their effective use in computational systems biology. The aim of model reduction is to circumvent this problem by eliminating parts of a model that are unimportant for the properties of interest. The choice of reduction method is influenced both by the type of model complexity and by the objective of the reduction; therefore, no single method is superior in all cases. In this study we present a comparative study of two different methods applied to a 20D model of yeast glycolytic oscillations. Our objective is to obtain biochemically meaningful reduced models, which reproduce the dynamic properties of the 20D model. The first method uses lumping and subsequent constrained parameter optimization. The second method is a novel approach that eliminates variables not essential for the dynamics. The applications of the two methods result in models of eight (lumping), six (elimination) and three (lumping followed by elimination) dimensions. All models have similar dynamic properties and pin-point the same interactions as being crucial for generation of the oscillations. The advantage of the novel method is that it is algorithmic, and does not require input in the form of biochemical knowledge. The lumping approach, however, is better at preserving biochemical properties, as we show through extensive analyses of the models.