- Research article
- Open Access
Integrated cellular network of transcription regulations and protein-protein interactions
© Wang and Chen; licensee BioMed Central Ltd. 2010
- Received: 6 August 2009
- Accepted: 8 March 2010
- Published: 8 March 2010
With the accumulation of increasing omics data, a key goal of systems biology is to construct networks at different cellular levels to investigate cellular machinery of the cell. However, there is currently no satisfactory method to construct an integrated cellular network that combines the gene regulatory network and the signaling regulatory pathway.
In this study, we integrated different kinds of omics data and developed a systematic method to construct the integrated cellular network based on coupling dynamic models and statistical assessments. The proposed method was applied to S. cerevisiae stress responses, elucidating the stress response mechanism of the yeast. From the resulting integrated cellular network under hyperosmotic stress, the highly connected hubs which are functionally relevant to the stress response were identified. Beyond hyperosmotic stress, the integrated network under heat shock and oxidative stress were also constructed and the crosstalks of these networks were analyzed, specifying the significance of some transcription factors to serve as the decision-making devices at the center of the bow-tie structure and the crucial role for rapid adaptation scheme to respond to stress. In addition, the predictive power of the proposed method was also demonstrated.
We successfully construct the integrated cellular network which is validated by literature evidences. The integration of transcription regulations and protein-protein interactions gives more insight into the actual biological network and is more predictive than those without integration. The method is shown to be powerful and flexible and can be used under different conditions and for different species. The coupling dynamic models of the whole integrated cellular network are very useful for theoretical analyses and for further experiments in the fields of network biology and synthetic biology.
- Gene Regulatory Network
- Hyperosmotic Stress
- Saccharomyces Genome Database
- Signaling Regulatory Pathway
- Stress Response Mechanism
With the advance of technologies (whole genome sequencing, expression profiling, and other high-throughput experiments), vast amounts of data, which cover different aspects of cellular physiology, have emerged [1, 2]. These kinds of 'omics' data, which include the genome sequencing data (genomics), microarray-based genome-wide expression profiles (transcriptomics), protein abundances data (proteomics), etc., provide unprecedented views of cellular components in the biological systems . However, the challenge for current researchers lies in how to interpret these large-scale data sets and extract true information to understand biological systems more thoroughly. With the amount of data accumulated, it is appropriate to understand living organisms from a system point of view. Computational techniques, which are able to combine these large and heterogeneous data sets, will provide us useful tools to gain more biological insights.
Though many studies have focused on gene regulatory networks [3–6] and on signaling transduction pathways [7–9], few works have combined these two kinds of networks. Yeger-Lotem and Margalit  used classical graph algorithms to integrate the cellular networks of protein-protein and protein-DNA interactions. They successfully identified known simple and complex regulatory circuits and discovered many putative circuits. Yeger-Lotem et al.  searched for composite network motifs consisting of both transcription-regulation and protein-protein interactions that recur significantly more often than in random networks. Mazurie et al.  investigated the integrated network comprising transcriptional and protein-protein interaction data. They compared network motifs for different species and showed that the network motifs are not subject to any particular evolutionary pressure to preserve the corresponding interaction patterns. Recently, researchers also have integrated protein-protein and protein-DNA interactions to infer signaling-regulatory pathways, but they focused only on the pathways that explain gene expression changes in response to the knockout genes . Although these works successfully extracted some characteristics of the integrated cellular network, they did not propose a method to construct the integrated cellular network of transcription regulations and protein-protein interactions under all kinds of conditions.
Many compelling biological questions center on how interactions among genes and proteins give rise to specific cellular functions. Genetic and high throughput methods have successfully identified many circuit components and their interactions. However, to explain the cellular function, there still are deficiencies, such as the incompleteness of the interactions among genes/proteins and the lack of consideration of different cellular conditions [14, 15]. To address this problem, we here propose a method to construct an integrated cellular network, which can be presented under all kinds of cellular conditions, based on coupling dynamic models. A dynamic model, which is typically represented by difference equations or differential equations, is often used in many fields to describe a complex and kinetic system. Recently, systems biology and computational biology methods have widely employed the dynamic model to describe the biological functions from the dynamic system point of view [3, 16, 17]. The advantage of using dynamic models to construct an integrated cellular network is not only that it can provide quantitative descriptions of the network, but it also can predict the behavior of the network under different conditions, i.e., gene knockout, treatment with an external agent, etc .
Two coupling dynamic models were given in this study to describe the interplay of gene regulatory network and signaling regulatory pathway. The Akaike Information Criterion (AIC) and the p-value statistical method [25–27] were used to determine the significant interactions via microarray data. By the interplay of the transcription regulations and protein-protein interactions, the cellular machinery can be elucidated comprehensively. The stress response mechanism is the most adequate example to demonstrate the integrated cellular network since signaling regulatory pathway sense and transmit the stress information quickly to the corresponding TFs to activate gene regulatory network, encoding protection proteins to respond to the stresses. Therefore, the proposed method was applied to yeast Saccharomyces cerevisiae to investigate the integrated mechanisms of stress responses.
Living organisms have evolved complex mechanisms to respond to changes in different environmental conditions, even for the unicellular organism yeast Saccharomyces cerevisiae[19, 28]. Such environmental changes, commonly termed as "stress", are harmful or even lethal to the survival of the cells, especially in microorganisms whose environment is highly variable. The responses to these stresses require complex networks of sensing and signal transduction pathways leading to adaptations of cell growth and proliferation as well as to adjustments of the gene expression program, metabolic activities, and other features of the cell . Several regulatory systems have been implicated in modulating these responses, but the complete network of regulators of stress responses and the details of their actions, including the signals that activate them and the downstream targets they regulated, remain to be elucidated . Consequently, these regulatory systems in response to environmental stresses are very suitable for the topic of integrated cellular network of transcription regulations and protein-protein interactions. In the microarray contributed by Gasch et al. , they identified the yeast stress responsive genes including stress-specific scheme and the common response to all of the stress conditions, termed the "environmental stress response" (ESR). However, they failed to determine the regulatory interactions between these genes and the dynamic characteristic of the system. In this paper, from the network or system point of view, we focused on yeast stress response and applied the proposed method to different kinds of stress (hyperosmotic stress, heat shock stress, oxidative stress), illustrating the integrated cellular network of yeast stress responses in which transcription regulations (gene regulatory network) and protein-protein interactions (signaling regulatory pathway) are integrated. Furthermore, the crosstalks of yeast stress responses were analyzed, specifying the significance of some transcription factors to serve as the decision-making devices at the center of the bow-tie structure and the crucial role for rapid adaptation scheme to respond to stress.
Construction of integrated cellular network
where f j denotes the sigmoid function, μ j and σ j represent the mean and deviation of protein activity level of the TF j. The biological meaning of the equation (1) is that the mRNA expression x i [t+1] of the target gene i at the next time t+1 is determined by the present mRNA expression x i [t] plus the mRNA expression due to the present transcription regulation of N TFs binding to the target gene at time t, minus the mRNA due to the degradation of the present time, plus the basal level of mRNA expression, and some stochastic noises due to measurement noises and some random fluctuations. Because of computational simplicity, the cooperative interaction between TFs is not included in the model .
where y n [t] represents the protein activity level at time t of the target protein n, b nm denotes the interaction ability of the m-th interactive protein to n-th target protein, y m [t] represents the protein activity level of m-th protein interacting with the target protein n, α n denotes the translation effect from mRNA to protein, x n [t] represents the mRNA expression level of the corresponding target protein n, β n indicates the degradation effect of the protein, h n represents the basal activity level, and ω n [t] is the stochastic noise. The rate of protein-protein interaction is proportional to the product of two proteins' concentrations , i.e., proportional to the probability of molecular collisions between two proteins, thus the interaction was modeled as the nonlinear multiplication scheme. For example, in the signaling transduction pathways, the phosphorylation of y n [t] by kinase y m [t] is proportional to the concentration of kinase y m [t] times the concentration of its substrate y n [t] . The biological meaning of the equation (3) is that the protein activity of the target protein n at the next time t+1 is contributed by the present protein activity, plus the regulatory interactions with M interactive proteins, plus the translation effect from mRNA, minus the protein degradation of the present state, plus basal protein level from other sources beyond the M interactive proteins in the cell, and some stochastic noises. Because of the undirected nature of protein interactions, there is no direction between interacting proteins in the protein-protein interaction subnetwork.
The interactions or couplings among genes and proteins are described in the following. Some TFs y j [t] at the end of signaling regulatory pathway regulate their target genes through the regulation function z j [t] = f j (y j [t]) in equation (1) and then the regulated genes influence their corresponding proteins through translation effect α n x n [t] in equation(3). The interplay between genes and proteins can be seen from their coupling dynamic equations (1), (3) and Figure 1, in which TFs serve as the interface between the signaling regulatory pathway and gene regulatory network. In other words, the interplay of transcription regulations and protein-protein interactions constitutes the framework of the integrated cellular network.
Based on the above dynamic models in (1) and (3), the candidate gene regulatory network can be linked through the regulatory parameters a ij in (1) between genes and their possible regulatory TFs, and the candidate signaling regulatory pathway can be linked through the interaction parameters b nm in (3) between possible interactive proteins. Since the candidate gene regulatory network and candidate signaling regulatory pathway only propose possible TF-gene regulations and possible protein interactions in omics data, their reality should be confirmed by microarray data, i.e., the values of a ij in (1) and b nm in (3) should be identified and validated by microarray data. The regulatory parameters were identified via microarray data by solving the constrained least square parameter estimation problem (see Methods). The strategy was to identify the regulatory parameters a ij gene by gene (interaction parameters b nm protein by protein), so that the network identification process was not limited by the number of target genes (proteins). In other words, when identifying the regulatory parameters in candidate gene regulatory network, the regulatory parameters in equation (1) were first identified for target gene i and then for target genes i+1, i+2, etc. By means of model selection method such as Akaike Information Criterion (AIC) and statistical assessment like student's t-test [25–27], the significant regulations and interactions between genes and proteins were detected (see Methods). In this way, the candidate gene regulatory network and the candidate signaling regulatory pathway were pruned by microarray data, leading to construction of the gene regulatory network and signaling regulatory pathway, respectively. Based on the interface between the gene regulatory network and the signaling regulatory pathway, i.e., the transcription factors, these two networks were coupled to become the integrated cellular network, which is the output of the proposed method (see Methods). The details of each step of the algorithm are shown in Methods.
Validation of the feasibility of the proposed method
Before applying the proposed method to construct the global view of the integrated cellular network under stress condition, we first validated the feasibility of the proposed method. While construction of gene regulatory networks using dynamic model has been successfully shown to characterize the biological system from a system point of view [3, 16, 17, 31, 32, 34], to the best of our knowledge, no research has used dynamic model to infer a signaling regulatory pathway. Thus, we focused on validating the protein interaction subnetwork of the integrated cellular network.
From the constructed subnetwork shown in Figure 3, we can find that there are two independent branches originating from Sho1p and Sln1p and converging to Pbs2p. This conforms to the observation made by Maeda et al. , who demonstrated that Pbs2p is activated by two independent signals that emanate from distinct cell-surface osmosensors by mutation analysis. In the Sln1p branch, Posas et al. showed that HOG MAPK cascade is regulated by a multistep phosphorelay mechanism in the Sln1p-Ypd1p-Ssk1p two component osmosensor . The phosphate group is transferred sequentially from Sln1p-His576 to Sln1p-Asp1144, then to Ypd1p-His64, and finally to Ssk1p-Asp554 . Posas and Saito also indicated that the Ssk2p/Ssk22p MAPKKKs are activated by the same two component osmosensor upon hyperosmotic treatment . Phosphorylated Ssk2p then activates the MAPKK Pbs2p [38, 40, 41]. Activated Pbs2p can phosphorylate and activate Hog1p . Among all these interactions observed by experiments, all were uncovered in the constructed subnetwork except for the Ypd1p-Ssk1p and Ssk1p-Ssk2p, which can be regarded as false negatives of the constructed subnetwork.
In the Sho1p branch, under hyperosmotic shock, Sho1p binds to Pbs2p and thereby recruits it to the cell surface . This event may mark the generation of the protein signaling complex that recruits the MAPKKK Ste11p together with proteins required for Ste11p activation such as Cdc42p [43, 44], Ste20p , Ste50p . The interaction of Sho1p and Pbs2p and the interactions among the protein complex were identified in the constructed subnetwork. The assembly of the appropriate protein signaling complex then leads to activation of Ste20p, phosphorylation of Ste11p and subsequently phosphorylation of Pbs2p . The interactions Ste20p-Ste11p, Ste11p-Pbs2p were also recognized. In yeast, the GTPase Cdc42p, together with its GTP exchange factor Cdc24p and its target Ste20p, is required to establish cell polarity during the cell cycle and is involved in cellular responses to mating pheromone and to nutritional limitation . The interactions between Cdc24p, Cdc42p, and Ste20p were uncovered in the constructed subnetwork, indicating that these proteins may also be required under hyperosmotic stress.
In addition to the well-described HOG MAPK pathway osmosensors Sho1p and Sln1p, the existence of other inputs into the HOG pathway has been suggested. O'Rourke and Herskowitz identified the Msb2p protein as a weak but physiologically relevant osmosensing components in S. cerevisiae which functions in parallel with the Shop1 branch and activate Ste11p . Recently, Tatebayashi et al. showed that the mucin-like transmembrane proteins Msb2p and Hkr1p are the potential osmosensors for the Shop1 branch . They demonstrated that hyperactive forms of Msb2p and Hkr1p can activate the HOG pathway only in the presence of Sho1p, whereas a hyperactive Sho1p mutant activates the HOG pathway in the absence of both Msb2p and Hkr1p, indicating that Msb2p and Hkr1p are the most upstream elements in the Sho1p branch . Although the interaction between Hkr1p and Sho1p was not uncovered, Msb2p was identified to interact with Sho1p in the constructed subnetwork. Besides, Msb2p can affect MAPKKK Ste11p by the interactions with Cdc42p, Ste20p, and Ste50p without the involvement of Sho1p, which is consistent with the experimental observation that Msb2p, but not Hkr1p can generate an intracellular signal in a Sho1p-independent manner .
Pbs2p is activated by phosphorylation by any of the three MAPKKKs Ssk2p, Ssk22p, and Ste11p and then phosphorylates the MAPK Hog1p. Phosphorylation promotes a rapid nuclear concentration of Hog1p, indicating that Hog1p is the factor that establishes the connection between the cytoplasmic part of signal pathway and its nuclear response system [48, 49]. In addition to nuclear translocation of the activated Hog1p, it also mediates regulatory effects outside the nucleus. For example, Hog1p regulates the activation of protein kinase Rck2p, which controls the translational efficiency. These two proteins are necessary for attenuation of protein synthesis in response to osmotic stress , and the interaction among them was observed in the constructed subnetwork. Several proteins influence nuclear residence of Hog1p, such as the tyrosine phosphatases Ptp2p and Ptp3p . Mattison and Ota showed that Ptp2p is a nuclear tether for Hog1p and Ptp3p is a cytoplasmic anchor, thus modulating the subcellular localization of Hog1p . The Hog1p-Ptp3p interaction was indicated in the constructed subnetwork whereas Hog1p-Ptp2p was absent. Besides the absence of Hog1p-Ptp2p interaction, Ste20p-Ptp2p interaction was recognized. Although Ptp2p has been shown to be a substrate of kinase Ste20p , no literature evidence correlates the interaction with HOG MAPK pathway until now. Thus, further investigation is needed to clarify the Ste20p-Ptp2p interaction. Three protein phosphatases Ptc1p, Ptc2p, and Ptc3p, inactivate the HOG MAPK pathway by acting on Hog1p [53, 54], however, only Ptc2p was detected to interact with Hog1p.
In summary, the constructed protein-protein interaction subnetwork was compared with the well-studied HOG MAPK pathway under osmotic stress to validate the feasibility of the proposed method. The comparison demonstrated that the well-studied MAPK pathway was mostly uncovered by the proposed method. However, there still were some interactions validated by literature evidences absent in the constructed subnetwork. The misidentification of the interaction may be because the protein activity levels were substituted by the gene expression levels (see Methods). The sensitivity of the HOG MAPK pathway is 78.57%, which shows that the signaling regulatory pathway can be successfully reconstructed by the proposed method, thus validating the feasibility of the proposed method. From the feasibility of the protein interaction network reconstruction along with the previous success of gene regulatory network characterization, we believe that the proposed method can be used to construct the integrated cellular network from the network or system point of view.
Global view of the yeast integrated cellular network under hyperosmotic stress
The proposed method was applied to yeast Saccharomyces cerevisiae hyperosmotic stress, illustrating the global view of yeast integrated cellular network under hyperosmotic stress in which transcription regulations (gene regulatory network) and protein-protein interactions (signaling regulatory pathway) are integrated. Based on the microarray data from Gasch et al. , 331 hyperosmotic stress responsive genes whose transcription levels change by threefold in at least one of the time courses were identified. Among the documented TFs obtained from the YEASTRACT database and ChIP-chip data from Harbison et al. , we identified 195 TFs which have expression profiles under hyperosmotic condition. The 195 TFs, 331 genes and 4333 TF-gene regulations retrieved were regarded as candidate gene network from which we constructed the gene regulatory network of the integrated cellular network. For the protein interaction subnetwork, 102 hyperosmotic stress-related proteins were selected from GO and SGD databases. We also retrieved 714 protein-protein interactions among the 102 proteins and 195 TFs from BioGRID database to be the candidate signaling pathway to construct the signaling regulatory pathway.
Expression profile reconstruction and prediction
Crosstalks of the yeast stress responses
Comparison of the three yeast integrated cellular networks under hyperosmotic stress, heat shock, and oxidative stress.
# TF-gene interactions
# Protein-protein interactions
Osmotic & Heat
Heat & Oxidative
Osmotic & Oxidative
Since the network reconstruction process is organism specific and is based on different datasets such as annotated genome sequence, high-throughput network-wide datasets and bibliomic data on the detailed properties of individual network components , it is not easy for biologists to efficiently construct the integrated cellular network by traditional genetic/molecular manipulation. By integrating different kinds of omics data including microarray gene expression data, regulatory associations between TFs and genes, and protein-protein interaction data, we have provided an efficient method based on coupling dynamic models and statistical assessments to construct an integrated cellular network in which transcription regulations (gene regulatory network) and protein-protein interactions (signaling regulatory pathway) are integrated. The stress response mechanism is the most adequate example to demonstrate the integrated cellular network because the signaling regulatory pathway is highly expressed to sense and transmit the stress information and the gene regulatory network is expressed to encode protection proteins to respond to the stress. Therefore the proposed method was applied to construct the integrated cellular network of S. cerevisiae under different stress conditions. Our results showed that the constructed integrated network can be validated by literature evidences, indicating the feasibility of the proposed method. The highly connected hubs which are functionally relevant to the hyperosmotic stress response were also identified and the predictive power of the proposed method was also demonstrated. Beyond hyperosmotic stress, the integrated network under heat shock and oxidative stress were constructed and the crosstalks of these networks were analyzed, which is not easily accomplished by traditional biological experiments.
In this work, a coupling dynamic model was used to characterize the yeast stress response mechanism from the network or system point of view. We provide a systematic and efficient way to move forwards studying the stress response networks rather than the stress response genes, which can be done by performing microarray experiments. However, one inconvenience of using this method is the lack of protein activity profiles. Thus there may be some inaccuracies when modeling the integrated cellular network with the substitution of gene expression levels. Even so, the well-studied HOG MAPK pathway was mostly uncovered by the proposed method and the results were validated by literature evidences. Once the genome-wide protein activity levels have been obtained, the accuracy can be improved. The advantage of the proposed method over others is its flexibility. It can be employed under all kinds of experimental conditions such as cell cycle or different kinds of diseases and for all kinds of organisms. Also, our method can be applied to identify and integrate the cellular networks of all sizes (see Results and Methods). In other words, the use of the proposed method is not limited by the number of genes/proteins. This benefit allows us to explore both a global view of the integrated network and a more detailed local view of the gene/protein regulation of interest. Since the global view of the integrated cellular network is not easily obtained by time- and labor-consuming biological experiments, the proposed method can be used by biologists to extract some useful information such as highly connected hubs which is possibly for further investigation. With slight modification, for example, by adding one clustering step before dynamic modeling, the proposed method can also be used to construct the module networks, in which sets of genes are co-regulated to respond to different conditions . The major advantage of dynamic modeling is the supply of quantitative dynamics of the network rather than qualitative descriptions. Because of the gene/protein dynamic system we constructed, computational methods can easily be employed to investigate the integrated gene/protein network. For example, we can make use of the dynamic systems to simulate any kind of perturbation such as gene knockout or the interference of external disturbances to see the system response for the perturbation. This can be very useful since the integrated gene/protein system can be modified arbitrarily to see the predicted consequences. In addition, the dynamics of the system can be used for theoretical analysis of the network such as robustness analysis and the signal transduction ability analysis. The system analysis and the computational simulation can be served as the pre-screen and the guidance for further experiments in the field of synthetic biology . Furthermore, the integrated cellular network can be mapped to the metabolic network to see how the cells make use of the metabolism to respond to different stresses and to unravel where and when certain parts of the network are active, elucidating the cellular machinery from a more comprehensive perspective [75, 76].
At present, there is no satisfactory method to construct the integrated cellular network, combining a gene regulatory network and a signaling regulatory pathway. We here provide a systematic and efficient method to integrate different kinds of omics data to construct the integrated cellular network via microarray data based on coupling dynamic models and statistical assessments. The integration of transcription regulations and protein-protein interactions gives more insight into actual biological networks and it is more predictive than those without integration . The integrated network construction represents an initial step for further network comparison and analysis. In our analyses of the S. cerevisiae stress response mechanism, the significance of some transcription factors to serve as the decision-making devices at the center of the bow-tie structure and the crucial role for rapid adaptation scheme to respond to stress are indicated. Furthermore, we also identify some genes/proteins which are relevant to the stress responses or are attractive targets for potential disease therapies . The method is shown to be powerful and flexible and can be used under different conditions and for different species. The dynamic models of the overall integrated cellular network provide a very useful tool for theoretical analyses and for further experiments in the fields of network biology and synthetic biology.
Data selection and preprocessing
There are three kinds of data used in the method--microarray gene expression data, regulatory associations between TFs and genes, and protein-protein interaction data. We used the genome-wide microarray data from Gasch et al.  as the stress responses data. They explored genomic expression patterns in the yeast Saccharomyces cerevisiae responding to diverse environmental transitions for almost every yeast gene, including temperature shocks, hyper- and hypo-osmotic shock, amino acid starvation, hydrogen peroxide shock, etc . The genes whose mRNA levels changed by threefold in at least one of the time courses were chosen as the significantly responsive genes. The responsive genes with more than 30% of the time points missing, which were regarded as unreliable experimental data, were excluded from further analyses. The transcription factors and the regulatory associations between TFs and target genes were obtained from YEASRTACT database http://www.yeastract.com/ and genome-wide location (ChIP-chip) data from Harbison et al. , in which the genomic occupancy of 203 DNA-binding TFs was determined in yeast. YEASTRACT (Yeast Search for Transcriptional Regulators And Consensus Tracking) is a curated repository of more than 34469 regulatory associations between TFs and target genes in Saccharomyces cerevisiae, based on more than 1000 bibliographic references . The significant binding of Harbison et al.  was selected as p < 0.001, as indicated in their paper. By the stress-responsive genes and the TF-gene regulatory associations obtained, we have the candidate gene regulatory network, which is the rough TF-gene regulation pool for building the gene regulatory network. The last dataset, protein-protein interaction data, were extracted from BioGRID database http://www.thebiogrid.org/. The Biological General Repository for Interaction Datasets (BioGRID) database was developed to house and distribute collections of protein and genetic interactions from major model organism species. BioGRID currently contains over 198000 interactions from six different species, as derived from both high-throughput studies and conventional focused studies . The stress-related proteins of the candidate integrated cellular network were selected from Gene Ontology (GO) http://www.geneontology.org/ and Saccharomyces Genome Database (SGD) http://www.yeastgenome.org/. Again, based on the stress-related proteins and the protein-protein interactions extracted, we have the candidate signaling regulatory pathway, which provides the rough protein interaction pool for constructing the signaling regulatory pathway.
In the microarray data from Gasch et al. , there are some missing values. The cubic spline interpolation method, which employs piecewise third-order polynomials to fit data points, was used to complement the missing values [78–80]. In order to fit the dynamic model, the microarray gene expression data was transformed from log2 scale to linear scale.
Dynamic model of the integrated cellular network
The candidate gene regulatory network and the candidate signaling regulatory pathway were constructed respectively (see Figure 2), and the stochastic discrete coupling dynamic models of these candidate networks are shown in equation (1) and equation (3). The interplay of the transcription regulations and protein-protein interactions constitutes the framework of the candidate integrated cellular network.
Identification of the regulatory parameters a ij and interaction parameters b nm
After constructing the coupling dynamic models of the candidate integrated cellular network, the regulatory/interaction parameters in the models have to be identified using the microarray data we have. The strategy is to identify the integrated cellular work gene by gene (protein by protein). Before determining the identification method, we first examine the dynamic models carefully. In equation (1), the basal expression level k i should be always non-negative, since the microarray expression of the genes are always non-negative. Because the parameters in equation (1) have certain constraints, the regulatory parameters were identified by solving the constrained least square problems.
where A = [0 ··· 0 0 -1], b = 0 give the constraints to force the basal level k i in equation (1) to be always non-negative, i.e., k i ≥ 0. The constrained least square problem can be solved using the active set method for quadratic programming [81, 82].
where C = diag [0 ··· 0 -1 0 -1] and d = [0 ··· 0] T , indicating that the translation effect α n and the basal activity level h n are non-negative. Since large-scale measurement of protein activities has yet to be realized and it was demonstrated that 73% of the variance in yeast protein abundance can be explained by mRNA abundance , mRNA expression profiles were used to substitute for the protein activity levels when identifying the interaction parameters. Due to the lack of protein activity data and the undirected nature of the protein interactions in the candidate signaling regulatory pathway, the sign of the regulatory parameters b nm 's does not imply activation or repression and there is no direction between interacting proteins.
Determination of significant interaction pairs
When the regulatory parameters were identified, Akaike Information Criterion (AIC) [25, 26] and student's t-test , which is used to calculate the p-values of the regulatory/interaction abilities, were employed for both model selection and determination of significant interactions in the integrated network. The AIC, which attempts to include both the estimated residual variance and model complexity in one statistics, decreases as the residual variance decreases and increases as the number of parameters increases. As the expected residual variance decreases with increasing parameter numbers for non-adequate model complexities, there should be a minimum around the correct parameter number [25, 26]. Therefore, AIC can be used to select model structure based on the regulatory abilities and the interaction abilities (a ij 's, b nm 's) identified above. Due to computational efficiency, it is impractical to compute the AIC statistics for all possible regression models. Stepwise methods such as forward selection method and backward elimination method are developed to avoid the complexity of exhaustive search [84, 85]. However, in the case of backward selection method, a variable once eliminated can never be reintroduced into the model, and in the case of forward selection, once included can never be removed . Thus, the stepwise regression method which combines forward selection method and backward elimination method was applied to compute the AIC statistics. Once the estimated regulatory/interaction parameters were examined using the AIC model selection criteria, student's t-test was employed to calculate the p-values  to determine the significant regulations/interactions. The p-values computed were then adjusted by Bonferroni correction to avoid a lot of spurious positives . The regulations/interactions which adjusted p-value ≦ 0.05 were determined as significant regulations/interactions and preserved in the integrated cellular network.
Combination of gene regulatory network with signaling regulatory pathway
which described one part of the integrated cellular network in Figure 1.
The study is supported by the National Science Council of Taiwan under grants NSC 97-2627-B-007-001 and NSC 98-2627-B-007-014.
- Joyce AR, Palsson BO: The model organism as a system: integrating 'omics' data sets. Nat Rev Mol Cell Biol. 2006, 7: 198-210. 10.1038/nrm1857View ArticlePubMedGoogle Scholar
- Ideker T, Bafna V, Lemberger T: Integrating scientific cultures. Mol Syst Biol. 2007, 3: 105- 10.1038/msb4100145PubMed CentralView ArticlePubMedGoogle Scholar
- Chen HC, Lee HC, Lin TY, Li WH, Chen BS: Quantitative characterization of the transcriptional regulatory network in the yeast cell cycle. Bioinformatics. 2004, 20: 1914-1927. 10.1093/bioinformatics/bth178View ArticlePubMedGoogle Scholar
- Workman CT, Mak HC, McCuine S, Tagne JB, Agarwal M, Ozier O, Begley TJ, Samson LD, Ideker T: A systems approach to mapping DNA damage response pathways. Science. 2006, 312: 1054-1059. 10.1126/science.1122088PubMed CentralView ArticlePubMedGoogle Scholar
- Xing B, Laan van der MJ: A causal inference approach for constructing transcriptional regulatory networks. Bioinformatics. 2005, 21: 4007-4013. 10.1093/bioinformatics/bti648View ArticlePubMedGoogle Scholar
- Veber P, Guziolowski C, Le Borgne M, Radulescu O, Siegel A: Inferring the role of transcription factors in regulatory networks. BMC Bioinformatics. 2008, 9: 228- 10.1186/1471-2105-9-228PubMed CentralView ArticlePubMedGoogle Scholar
- Draghici S, Khatri P, Tarca AL, Amin K, Done A, Voichita C, Georgescu C, Romero R: A systems biology approach for pathway level analysis. Genome Res. 2007, 17: 1537-1545. 10.1101/gr.6202607PubMed CentralView ArticlePubMedGoogle Scholar
- Visvanathan M, Breit M, Pfeifer B, Baumgartner C, Modre-Osprian R, Tilg B: Systematic analysis of signaling pathways using an integrative environment. Methods Inf Med. 2007, 46: 386-391.PubMedGoogle Scholar
- Kuepfer L, Peter M, Sauer U, Stelling J: Ensemble modeling for analysis of cell signaling dynamics. Nat Biotechnol. 2007, 25: 1001-1006. 10.1038/nbt1330View ArticlePubMedGoogle Scholar
- Yeger-Lotem E, Margalit H: Detection of regulatory circuits by integrating the cellular networks of protein-protein interactions and transcription regulation. Nucleic Acids Res. 2003, 31: 6053-6061. 10.1093/nar/gkg787PubMed CentralView ArticlePubMedGoogle Scholar
- Yeger-Lotem E, Sattath S, Kashtan N, Itzkovitz S, Milo R, Pinter RY, Alon U, Margalit H: Network motifs in integrated cellular networks of transcription-regulation and protein-protein interaction. Proc Natl Acad Sci USA. 2004, 101: 5934-5939. 10.1073/pnas.0306752101PubMed CentralView ArticlePubMedGoogle Scholar
- Mazurie A, Bottani S, Vergassola M: An evolutionary and functional assessment of regulatory network motifs. Genome Biol. 2005, 6: R35- 10.1186/gb-2005-6-4-r35PubMed CentralView ArticlePubMedGoogle Scholar
- Ourfali O, Shlomi T, Ideker T, Ruppin E, Sharan R: SPINE: a framework for signaling-regulatory pathway inference from cause-effect experiments. Bioinformatics. 2007, 23: i359-366. 10.1093/bioinformatics/btm170View ArticlePubMedGoogle Scholar
- Sprinzak D, Elowitz MB: Reconstruction of genetic circuits. Nature. 2005, 438: 443-448. 10.1038/nature04335View ArticlePubMedGoogle Scholar
- Feist AM, Herrgard MJ, Thiele I, Reed JL, Palsson BO: Reconstruction of biochemical networks in microorganisms. Nat Rev Microbiol. 2009, 7: 129-143.PubMed CentralView ArticlePubMedGoogle Scholar
- Chang YH, Wang YC, Chen BS: Identification of transcription factor cooperativity via stochastic system model. Bioinformatics. 2006, 22: 2276-2282. 10.1093/bioinformatics/btl380View ArticlePubMedGoogle Scholar
- Bansal M, Gatta GD, di Bernardo D: Inference of gene regulatory networks and compound mode of action from time course gene expression profiles. Bioinformatics. 2006, 22: 815-822. 10.1093/bioinformatics/btl003View ArticlePubMedGoogle Scholar
- Bansal M, Belcastro V, Ambesi-Impiombato A, di Bernardo D: How to infer gene networks from expression profiles. Mol Syst Biol. 2007, 3:Google Scholar
- Gasch AP, Spellman PT, Kao CM, Carmel-Harel O, Eisen MB, Storz G, Botstein D, Brown PO: Genomic expression programs in the response of yeast cells to environmental changes. Mol Biol Cell. 2000, 11: 4241-4257.PubMed CentralView ArticlePubMedGoogle Scholar
- Harbison CT, Gordon DB, Lee TI, Rinaldi NJ, Macisaac KD, Danford TW, Hannett NM, Tagne JB, Reynolds DB, Yoo J, et al.: Transcriptional regulatory code of a eukaryotic genome. Nature. 2004, 431: 99-104. 10.1038/nature02800PubMed CentralView ArticlePubMedGoogle Scholar
- Cherry JM, Adler C, Ball C, Chervitz SA, Dwight SS, Hester ET, Jia Y, Juvik G, Roe T, Schroeder M, et al.: SGD: Saccharomyces Genome Database. Nucl Acids Res. 1998, 26: 73-79. 10.1093/nar/26.1.73PubMed CentralView ArticlePubMedGoogle Scholar
- Teixeira MC, Monteiro P, Jain P, Tenreiro S, Fernandes AR, Mira NP, Alenquer M, Freitas AT, Oliveira AL, Sa-Correia I: The YEASTRACT database: a tool for the analysis of transcription regulatory associations in Saccharomyces cerevisiae. Nucleic Acids Res. 2006, 34: D446-451. 10.1093/nar/gkj013PubMed CentralView ArticlePubMedGoogle Scholar
- Stark C, Breitkreutz B-J, Reguly T, Boucher L, Breitkreutz A, Tyers M: BioGRID: a general repository for interaction datasets. Nucl Acids Res. 2006, 34: D535-539. 10.1093/nar/gkj109PubMed CentralView ArticlePubMedGoogle Scholar
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al.: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000, 25: 25-29. 10.1038/75556PubMed CentralView ArticlePubMedGoogle Scholar
- Akaike H: A new look at the statistical model identification. IEEE Transactions on Automatic Control. 1974, 19: 716-723. 10.1109/TAC.1974.1100705.View ArticleGoogle Scholar
- Johansson R: System modeling and identification. 1993, Englewood Cliffs, N.J.: Prentice HallGoogle Scholar
- Mendenhall W, Sincich T: Statistics for engineering and the sciences. 1995, Englewood Cliffs, N.J.: Prentice-Hall, 4Google Scholar
- Hohmann S, Mager WH: Yeast stress responses. 2003, Berlin; New York: SpringerView ArticleGoogle Scholar
- Wu WS, Li WH, Chen BS: Identifying regulatory targets of cell cycle transcription factors using gene expression and ChIP-chip data. BMC Bioinformatics. 2007, 8: 188- 10.1186/1471-2105-8-188PubMed CentralView ArticlePubMedGoogle Scholar
- Alon U: An introduction to systems biology: design principles of biological circuits. 2007, Boca Raton, FL: Chapman & Hall/CRCGoogle Scholar
- Chen KC, Wang TY, Tseng HH, Huang CY, Kao CY: A stochastic differential equation model for quantifying transcriptional regulatory network in Saccharomyces cerevisiae. Bioinformatics. 2005, 21: 2883-2890. 10.1093/bioinformatics/bti415View ArticlePubMedGoogle Scholar
- Wu WS, Li WH, Chen BS: Reconstructing a network of stress-response regulators via dynamic system modeling of gene regulation. Gene Regulation and Systems Biology. 2008, 2: 53-62.PubMed CentralPubMedGoogle Scholar
- de Jong H: Modeling and simulation of genetic regulatory systems: a literature review. J Comput Biol. 2002, 9: 67-103. 10.1089/10665270252833208View ArticlePubMedGoogle Scholar
- Tegner J, Yeung MK, Hasty J, Collins JJ: Reverse engineering gene networks: integrating genetic perturbations with dynamical modeling. Proc Natl Acad Sci USA. 2003, 100: 5944-5949. 10.1073/pnas.0933416100PubMed CentralView ArticlePubMedGoogle Scholar
- Roberts CJ, Nelson B, Marton MJ, Stoughton R, Meyer MR, Bennett HA, He YD, Dai H, Walker WL, Hughes TR, et al.: Signaling and circuitry of multiple MAPK pathways revealed by a matrix of global gene expression profiles. Science. 2000, 287: 873-880. 10.1126/science.287.5454.873View ArticlePubMedGoogle Scholar
- Steffen M, Petti A, Aach J, D'Haeseleer P, Church G: Automated modelling of signal transduction networks. BMC Bioinformatics. 2002, 3: 34- 10.1186/1471-2105-3-34PubMed CentralView ArticlePubMedGoogle Scholar
- Hersen P, McClean MN, Mahadevan L, Ramanathan S: Signal processing by the HOG MAP kinase pathway. Proc Natl Acad Sci USA. 2008, 105: 7165-7170. 10.1073/pnas.0710770105PubMed CentralView ArticlePubMedGoogle Scholar
- Maeda T, Takekawa M, Saito H: Activation of yeast PBS2 MAPKK by MAPKKKs or by binding of an SH3-containing osmosensor. Science. 1995, 269: 554-558. 10.1126/science.7624781View ArticlePubMedGoogle Scholar
- Posas F, Wurgler-Murphy SM, Maeda T, Witten EA, Thai TC, Saito H: Yeast HOG1 MAP kinase cascade is regulated by a multistep phosphorelay mechanism in the SLN1-YPD1-SSK1 "two-component" osmosensor. Cell. 1996, 86: 865-875. 10.1016/S0092-8674(00)80162-2View ArticlePubMedGoogle Scholar
- Posas F, Saito H: Activation of the yeast SSK2 MAP kinase kinase kinase by the SSK1 two-component response regulator. EMBO J. 1998, 17: 1385-1394. 10.1093/emboj/17.5.1385PubMed CentralView ArticlePubMedGoogle Scholar
- Maeda T, Wurgler-Murphy SM, Saito H: A two-component system that regulates an osmosensing MAP kinase cascade in yeast. Nature. 1994, 369: 242-245. 10.1038/369242a0View ArticlePubMedGoogle Scholar
- Brewster JL, de Valoir T, Dwyer ND, Winter E, Gustin MC: An osmosensing signal transduction pathway in yeast. Science. 1993, 259: 1760-1763. 10.1126/science.7681220View ArticlePubMedGoogle Scholar
- Reiser V, Salah SM, Ammerer G: Polarized localization of yeast Pbs2 depends on osmostress, the membrane protein Sho1 and Cdc42. Nat Cell Biol. 2000, 2: 620-627. 10.1038/35023568View ArticlePubMedGoogle Scholar
- Raitt DC, Posas F, Saito H: Yeast Cdc42 GTPase and Ste20 PAK-like kinase regulate Sho1-dependent activation of the Hog1 MAPK pathway. EMBO J. 2000, 19: 4623-4631. 10.1093/emboj/19.17.4623PubMed CentralView ArticlePubMedGoogle Scholar
- Posas F, Witten EA, Saito H: Requirement of STE50 for osmostress-induced activation of the STE11 mitogen-activated protein kinase kinase kinase in the high-osmolarity glycerol response pathway. Mol Cell Biol. 1998, 18: 5788-5796.PubMed CentralView ArticlePubMedGoogle Scholar
- O'Rourke SM, Herskowitz I: A third osmosensing branch in Saccharomyces cerevisiae requires the Msb2 protein and functions in parallel with the Sho1 branch. Mol Cell Biol. 2002, 22: 4739-4749. 10.1128/MCB.22.13.4739-4749.2002PubMed CentralView ArticlePubMedGoogle Scholar
- Tatebayashi K, Tanaka K, Yang HY, Yamamoto K, Matsushita Y, Tomida T, Imai M, Saito H: Transmembrane mucins Hkr1 and Msb2 are putative osmosensors in the SHO1 branch of yeast HOG pathway. EMBO J. 2007, 26: 3521-3533. 10.1038/sj.emboj.7601796PubMed CentralView ArticlePubMedGoogle Scholar
- Reiser V, Ruis H, Ammerer G: Kinase activity-dependent nuclear export opposes stress-induced nuclear accumulation and retention of Hog1 mitogen-activated protein kinase in the budding yeast Saccharomyces cerevisiae. Mol Biol Cell. 1999, 10: 1147-1161.PubMed CentralView ArticlePubMedGoogle Scholar
- Ferrigno P, Posas F, Koepp D, Saito H, Silver PA: Regulated nucleo/cytoplasmic exchange of HOG1 MAPK requires the importin beta homologs NMD5 and XPO1. EMBO J. 1998, 17: 5606-5614. 10.1093/emboj/17.19.5606PubMed CentralView ArticlePubMedGoogle Scholar
- Teige M, Scheikl E, Reiser V, Ruis H, Ammerer G: Rck2, a member of the calmodulin-protein kinase family, links protein synthesis to high osmolarity MAP kinase signaling in budding yeast. Proc Natl Acad Sci USA. 2001, 98: 5625-5630. 10.1073/pnas.091610798PubMed CentralView ArticlePubMedGoogle Scholar
- Mattison CP, Ota IM: Two protein tyrosine phosphatases, Ptp2 and Ptp3, modulate the subcellular localization of the Hog1 MAP kinase in yeast. Genes Dev. 2000, 14: 1229-1235.PubMed CentralPubMedGoogle Scholar
- Ptacek J, Devgan G, Michaud G, Zhu H, Zhu X, Fasolo J, Guo H, Jona G, Breitkreutz A, Sopko R, et al.: Global analysis of protein phosphorylation in yeast. Nature. 2005, 438: 679-684. 10.1038/nature04187View ArticlePubMedGoogle Scholar
- Young C, Mapes J, Hanneman J, Al-Zarban S, Ota I: Role of Ptc2 type 2C Ser/Thr phosphatase in yeast high-osmolarity glycerol pathway inactivation. Eukaryot Cell. 2002, 1: 1032-1040. 10.1128/EC.1.6.1032-1040.2002PubMed CentralView ArticlePubMedGoogle Scholar
- Warmka J, Hanneman J, Lee J, Amin D, Ota I: Ptc1, a type 2C Ser/Thr phosphatase, inactivates the HOG pathway by dephosphorylating the mitogen-activated protein kinase Hog1. Mol Cell Biol. 2001, 21: 51-60. 10.1128/MCB.21.1.51-60.2001PubMed CentralView ArticlePubMedGoogle Scholar
- Barabasi AL, Oltvai ZN: Network biology: understanding the cell's functional organization. Nat Rev Genet. 2004, 5: 101-113. 10.1038/nrg1272View ArticlePubMedGoogle Scholar
- Jeong H, Mason SP, Barabasi AL, Oltvai ZN: Lethality and centrality in protein networks. Nature. 2001, 411: 41-42. 10.1038/35075138View ArticlePubMedGoogle Scholar
- Chen BS, Yang SK, Lan CY, Chuang YJ: A systems biology approach to construct the gene regulatory network of systemic inflammation via microarray and databases mining. BMC Med Genomics. 2008, 1: 46- 10.1186/1755-8794-1-46PubMed CentralView ArticlePubMedGoogle Scholar
- Yang XX, Maurer KC, Molanus M, Mager WH, Siderius M, Vies van der SM: The molecular chaperone Hsp90 is required for high osmotic stress response in Saccharomyces cerevisiae. FEMS Yeast Res. 2006, 6: 195-204. 10.1111/j.1567-1364.2006.00026.xView ArticlePubMedGoogle Scholar
- Pearl LH: Hsp90 and Cdc37 -- a chaperone cancer conspiracy. Curr Opin Genet Dev. 2005, 15: 55-61. 10.1016/j.gde.2004.12.011View ArticlePubMedGoogle Scholar
- Varela J, Praekelt U, Meacock P, Planta R, Mager W: The Saccharomyces cerevisiae HSP12 gene is activated by the high-osmolarity glycerol pathway and negatively regulated by protein kinase A. Mol Cell Biol. 1995, 15: 6232-6245.PubMed CentralView ArticlePubMedGoogle Scholar
- Karreman RJ, Dague E, Gaboriaud F, Quiles F, Duval JF, Lindsey GG: The stress response protein Hsp12p increases the flexibility of the yeast Saccharomyces cerevisiae cell wall. Biochim Biophys Acta. 2007, 1774: 131-137.View ArticlePubMedGoogle Scholar
- Schuller C, Brewster JL, Alexander MR, Gustin MC, Ruis H: The HOG pathway controls osmotic regulation of transcription via the stress response element (STRE) of the Saccharomyces cerevisiae CTT1 gene. EMBO J. 1994, 13: 4382-4389.PubMed CentralPubMedGoogle Scholar
- Marchler G, Schuller C, Adam G, Ruis H: A Saccharomyces cerevisiae UAS element controlled by protein kinase A activates transcription in response to a variety of stress conditions. EMBO J. 1993, 12: 1997-2003.PubMed CentralPubMedGoogle Scholar
- Chen BS, Li CW: Analysing microarray data in drug discovery using systems biology. Expert Opinion on Drug Discovery. 2007, 2: 755-768. 10.1517/17460418.104.22.1685.View ArticlePubMedGoogle Scholar
- Diderich JA, Schuurmans JM, Van Gaalen MC, Kruckeberg AL, Van Dam K: Functional analysis of the hexose transporter homologue HXT5 in Saccharomyces cerevisiae. Yeast. 2001, 18: 1515-1524. 10.1002/yea.779View ArticlePubMedGoogle Scholar
- Verwaal R, Paalman JW, Hogenkamp A, Verkleij AJ, Verrips CT, Boonstra J: HXT5 expression is determined by growth rates in Saccharomyces cerevisiae. Yeast. 2002, 19: 1029-1038. 10.1002/yea.895View ArticlePubMedGoogle Scholar
- Borkovich KA, Farrelly FW, Finkelstein DB, Taulien J, Lindquist S: hsp82 is an essential protein that is required in higher concentrations for growth of cells at higher temperatures. Mol Cell Biol. 1989, 9: 3919-3930.PubMed CentralView ArticlePubMedGoogle Scholar
- Picard D: Heat-shock protein 90, a chaperone for folding and regulation. Cell Mol Life Sci. 2002, 59: 1640-1648. 10.1007/PL00012491View ArticlePubMedGoogle Scholar
- Csete M, Doyle J: Bow ties, metabolism and disease. Trends Biotechnol. 2004, 22: 446-450. 10.1016/j.tibtech.2004.07.007View ArticlePubMedGoogle Scholar
- Li R, Brawley S: Improved survival under heat stress in intertidal embryos (Fucus spp.) simultaneously exposed to hypersalinity and the effect of parental thermal history. Marine Biology. 2004, 144: 205-213. 10.1007/s00227-003-1190-9.View ArticleGoogle Scholar
- Xiong L, Ishitani M, Zhu JK: Interaction of osmotic stress, temperature, and abscisic acid in the regulation of gene expression in Arabidopsis. Plant Physiol. 1999, 119: 205-212. 10.1104/pp.119.1.205PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang S, Jin G, Zhang XS, Chen L: Discovering functions and revealing mechanisms at molecular level from biological networks. Proteomics. 2007, 7: 2856-2869. 10.1002/pmic.200700095View ArticlePubMedGoogle Scholar
- Segal E, Shapira M, Regev A, Pe'er D, Botstein D, Koller D, Friedman N: Module networks: identifying regulatory modules and their condition-specific regulators from gene expression data. NatGenet. 2003, 34: 166-176.Google Scholar
- Cuccato G, Della Gatta G, di Bernardo D: Systems and Synthetic biology: tackling genetic networks and complex diseases. Heredity. 2009, 102: 527-532. 10.1038/hdy.2009.18View ArticlePubMedGoogle Scholar
- Chechik G, Oh E, Rando O, Weissman J, Regev A, Koller D: Activity motifs reveal principles of timing in transcriptional control of the yeast metabolic network. Nat Biotechnol. 2008, 26: 1251-1259. 10.1038/nbt.1499PubMed CentralView ArticlePubMedGoogle Scholar
- Walhout AJ: Unraveling transcription regulatory networks by protein-DNA and protein-protein interaction mapping. Genome Res. 2006, 16: 1445-1454. 10.1101/gr.5321506View ArticlePubMedGoogle Scholar
- Zhu J, Zhang B, Smith EN, Drees B, Brem RB, Kruglyak L, Bumgarner RE, Schadt EE: Integrating large-scale functional genomic data to dissect the complexity of yeast regulatory networks. Nat Genet. 2008, 40: 854-861. 10.1038/ng.167PubMed CentralView ArticlePubMedGoogle Scholar
- Bar-Joseph Z, Gerber G, Gifford DK, Jaakkola TS, Simon I: A new approach to analyzing gene expression time series data. Proceedings of the sixth annual international conference on Computational biology. 2002, Washington, DC, USA: ACMGoogle Scholar
- Faires JD, Burden RL: Numerical methods. 1998, Pacific Grove, CA: Brooks/Cole Pub. Co., 2Google Scholar
- de Boor C: A practical guide to splines Rev edn. 2001, New York: SpringerGoogle Scholar
- Coleman TF, Hulbert LA: A direct active set algorithm for large sparse quadratic programs with simple bounds. Mathematical Programming. 1989, 45: 373-406. 10.1007/BF01589112.View ArticleGoogle Scholar
- Gill PE, Murray W, Wright MH: Practical optimization. 1981, London; New York: Academic PressGoogle Scholar
- Lu P, Vogel C, Wang R, Yao X, Marcotte EM: Absolute protein expression profiling estimates the relative contributions of transcriptional and translational regulation. Nat Biotechnol. 2007, 25: 117-124. 10.1038/nbt1270View ArticlePubMedGoogle Scholar
- Hocking RR: A Biometrics invited paper. The analysis and selection of variables in linear regression. Biometrics. 1976, 32: 1-49. 10.2307/2529336.View ArticleGoogle Scholar
- Seber GAF, Lee AJ: Linear regression analysis. 2003, Hoboken, N.J.: Wiley-Interscience, 2View ArticleGoogle Scholar
- Barsky A, Gardy JL, Hancock RE, Munzner T: Cerebral: a Cytoscape plugin for layout of and interaction with biological networks using subcellular localization annotation. Bioinformatics. 2007, 23: 1040-1042. 10.1093/bioinformatics/btm057View ArticlePubMedGoogle Scholar
- Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T: Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13: 2498-2504. 10.1101/gr.1239303PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.