Phenotype analysis using network motifs derived from changes in regulatory network dynamics

Abstract

The intrinsic dynamic response of a transcriptional regulatory network depends directly on molecular interactions in the cellular transcription, translation, and degradation machineries. These interactions can be incorporated into dynamic mathematical models of the biochemical system using the biophysical relationship with the model parameters. Modifications of such interactions bring changes to the biological behavior of the cells, and therefore, many normal and pathological cellular states depend on them. It is important for analysis, prediction, diagnosis, and treatment of cellular function to have an experimentally derived model with parameters that adequately represent the molecular interactions of interest. Finding the model and parameters of a transcriptional regulatory network is a difficult task that has been approached at different levels and with different techniques. We develop here a new analysis method (based on previous work on network inference, modeling, and parameter identification) that finds the most changed parameters from yeast oligonucleotide microarray expression patterns in cases where a phenotype difference exists between two samples. We then relate and examine the changed parameters with their associated genes, corresponding genetic functional categories, and particular subnetworks and connectivities. The biophysical bases for these changes are also identified by studying the relationship of the changed parameters with the transcription, translation, and degradation mechanisms. The method is improved to cases where there are two or more transcription factors influencing transcription, and a statistical analysis is performed to give a measurement of the uniqueness and robustness of the parameter fit.