Steady state approach to model gene regulatory networks--simulation of microarray experiments

Genetic regulatory networks (GRN) represent complex interactions between genes brought about through proteins that they code for. Quantification of expression levels in GRN either through experiments or theoretical modeling is a challenging task. Recently, microarray experiments have gained importance in evaluating GRN at the genome level. Microarray experiments yield log fold change in mRNA abundance which is helpful in deciphering connectivity in GRN. Current approaches such as data mining, Boolean or Bayesian modeling and combined use of expression and location data are useful in analyzing microarray data. However, these methodologies lack underlying mechanistic details present in GRN. We present here a steady state gene expression simulator (SSGES) which sets up steady state equations and simulates the response for a given network structure of a GRN. SSGES includes mechanistic details such as stoichiometry, protein-DNA and protein-protein interactions, translocation of regulatory proteins and autoregulation. SSGES can be used to simulate the response of a GRN in terms of fractional transcription and protein expression. SSGES can also be used to generate log fold change in mRNA abundance and protein expression implying that it is useful to simulate microarray type experiments. We have demonstrated these capabilities of SSGES by modeling the steady state response of GAL regulatory system in Saccharomyces cerevisiae. We have demonstrated that the predicted data qualitatively matched the microarray data obtained experimentally by Ideker et al. [Ideker, T., Thorsson, V., Ranish, J.A., Christmas, R., Buhler, J., Eng, J.K., Bumgarner, R., Goodlett, D.R., Aebersold, R., Hood, L., 2001. Integrated genomic and proteomic analyses of a systematically perturbed metabolic network. Science 292, 929-934]. SSGES is available from authors upon request.