EFFECTIVE MODELS FOR GENE NETWORKS AND THEIR APPLICATIONS

Abstract

Networks are used to abstract large, highly-coupled sets of objects. Their analyses have included network classification into a few broad classes and selection of small substructures that perform simple yet common tasks. One issue that has received little attention is how the state of a network can be moved according to a pre-specified set of guidelines. In this paper, we address this question in the context of gene networks. In general, neither the full membership of the gene network associated with a biological process nor the precise form of interactions between nodes is known. What is available, through microarrays or sequencing, are gene expression profiles of an organism or its viable mutants. Our approach relies only on these expression profiles, and not on the availability of an accurate model for the network. The first step is to select a small set of core- or master- nodes, such as transcription factors or microRNAs, that can be used to alter the levels of many of the remaining genes in the network. We ask how the state — or solution — of the gene network changes as the levels of these master nodes are altered externally. The object of our study is, not the network, but the surface of these solutions. We argue that it can be approximated using gene expression profiles of the organism and single manipulation of master node activity. This is done through an "effective model." The effective model as well as error estimates for its predictions can be derived from experimental data. The method is validated using synthetic gene networks that have stationary solutions and those that are periodically driven, e.g., circadian networks. An effective model for the oxygen-deprivation network of E.coli is constructed using previously published gene expression profiles, and used to predict the expression levels in a double knockout mutant. Less that 30% of the predictions lie outside the 5% confidence level. We propose the use of the effective model methodology to compute how Drosophila melanogaster in the normal state can be genetically altered into a pre-defined sleep deprived-like state.