A computational framework for gene regulatory network inference that combines multiple methods and datasets
© Gupta et al; licensee BioMed Central Ltd. 2011
Received: 7 January 2011
Accepted: 13 April 2011
Published: 13 April 2011
Reverse engineering in systems biology entails inference of gene regulatory networks from observational data. This data typically include gene expression measurements of wild type and mutant cells in response to a given stimulus. It has been shown that when more than one type of experiment is used in the network inference process the accuracy is higher. Therefore the development of generally applicable and effective methodologies that embed multiple sources of information in a single computational framework is a worthwhile objective.
This paper presents a new method for network inference, which uses multi-objective optimisation (MOO) to integrate multiple inference methods and experiments. We illustrate the potential of the methodology by combining ODE and correlation-based network inference procedures as well as time course and gene inactivation experiments. Here we show that our methodology is effective for a wide spectrum of data sets and method integration strategies.
The approach we present in this paper is flexible and can be used in any scenario that benefits from integration of multiple sources of information and modelling procedures in the inference process. Moreover, the application of this method to two case studies representative of bacteria and vertebrate systems has shown potential in identifying key regulators of important biological processes.
In the last ten years the development of functional genomics technologies has provided us with the ability to generate quantitative data representing the molecular state of cells and tissues at a genome level [1, 2]. These datasets can be in the form of a time series representing the dynamics of gene expression profiles (e.g. mRNA, proteins and metabolites) in response to a given stimulus, such as an environmental perturbation, the effect of a growth factor or an experimentally induced gene deletion. Despite the relatively large amount of information, predicting underlying regulatory networks from observational data is still not trivial and is a matter of intense research .
A number of reverse-engineering approaches have been proposed. Some of these are designed to infer networks from a compendium of perturbation experiments while others are able to use time course data to develop dynamical models of gene interaction. Bayesian networks have been among the first to be applied to biological problems . They work by inferring probabilistic relationships between variables, can use either time course or steady state data and allow integration of prior knowledge in the model. Correlation-based methods [5, 6] compute correlation coefficients between variables to infer the underlying network topology. State-space models (SSMs) [7, 8], and ODE-based methods [9, 10], on the other hand use time-course data to develop dynamic models of gene regulatory networks (GRN). For an extensive overview of these methodologies see: [11, 12].
The general validity of the principal of integrating multiple data sources in the reverse-engineering process is exemplified by the observation that the best performing methods utilize some degree of integration between different experiments . For example, the top performing method in the third edition of the "Dialogue for Reverse Engineering Assessments and Methods" (DREAM), developed by Yip et al. , was based on a combination of a statistical error-model and ODE modeling to integrate gene knock-out (KO) and time-course experiments. Interestingly, Yip et al.  also noted that a relatively simple differential gene-expression analysis, comparing wild-type and mutant strains, was in itself a very good representation of the underlying gene regulatory network. However, not all KO experiments are likely to be equally informative and identifying a priori the most relevant genes is not a trivial task. Moreover, large-scale gene-inactivation experiments are not a viable option for many non-model species.
Therefore, there is the need to expand the repertoire of available network inference tools by developing more methods that allow integration of multiple data sources and have the flexibility to use a wide range of datasets and information. In order to achieve this objective, we set out to develop a computational framework that has the potential to combine different inference methodologies, multiple datasets, as well as any pre-existing biological knowledge. We based this approach on an ODE framework combined to a multi-objective optimization (MOO) procedure for parameter estimation. We named this method "Network-Inference with Multi Objective Optimization" (NIMOO).
The basic network inference framework: Model Equations and parameter estimation of a single objective optimization procedure
The gene regulatory network (GRN) is then inferred from the optimized parameter matrix w. The matrix element |wij| indicates the strength of the interaction between genes i and j (with gene j regulating gene i), and, sign (wij) indicates whether the effect is excitatory (wij > 0) or inhibitory (wij < 0)
In our implementation of single objective optimisation (SOO), minimisation of E SQE was achieved using the trust-region method based on the interior-reflective Newton method [15, 16]. In this method the minimisation process involves defining a trust region where the objective function SQE can be approximated with a simpler function q. For successive iterations, function q, in conjunction with the Preconditioned Conjugate Gradient Method , is used to find a new trust region where the function SQE is lower. The process is terminated when the change in function value is less than a pre-determined tolerance (10-6).
Parameter estimation using a multi-objective optimization procedure
The term θ k γ introduces flexibility in the degree of goals attained. Also, the weight factor θk can be used to assign relative importance to the objectives: Thus, from Equation 4, θ k = 0 implies hard goal for the corresponding objective function F k . For all results presented in this paper, unless mentioned otherwise, the goal and weight corresponding to the objective SQE were set to 0.
Overall strategy for the development of a MOO-based inference method
The principle behind NIMOO, as detailed in the above sections, is to infer the gene regulatory matrix (GRM) w by minimizing the E SQE for the ODE system in conjunction with additional objective functions, E object , which represent the distance between the parameter matrix wij and objectives Oij constructed from additional information.
In principal, objectives Oij can be constructed from any data available on the underlying regulatory network. In this paper, we focused on two possible scenarios.
In the first case, we considered the possibility that MOO might be used to integrate two different network-inference procedures, for example applied to independent replicates of a time-course experiment. In this implementation we used time-delayed Spearman rank-correlation  to develop a matrix Oij (Equation 3) representing the degree of statistical correlation from any pair of genes within a set time delay interval (Figure 1, objective DSp).
Alternatively, O ij can be built from the results of gene inactivation experiments. We reasoned that these experiments might fall in at least two categories. In the first case the gene is deleted at some stage of the lifecycle of the organism so that gene expression measurements can only be acquired after the new steady state has been reached. This can be easily achieved, by a plethora of gene knockouts (KOs) methodologies in model systems ranging from E. coli to mice. Alternatively, gene inactivation could be achieved by using biochemical inhibitors or RNA interference. In this scenario mRNA expression can be monitored at different time points following intervention. In the context of MOO both scenarios lead to a gene-expression matrices where rows are represented by genes and columns by gene KO experiments. From each of these matrices an objective Oij can be computed to represent the relationship between every gene pair (Figure 1, objectives Tc, Tr are derived from expression data at a given time-point shortly after gene inactivation whereas objectives Sc and Sr derive from expression data at a single time point at steady state; for further details on how to compute Oij see sections below).
Different procedures may be used in combination using an ensemble approach; in this paper we describe the results of combining MOO-Tr with MOO-Tc (Figure 1, MOO-Tens) and MOO-Sr with MOO-Sc (Figure 1, MOO-Sens).
All MOO procedures developed within NIMOO have been compared with the accuracy of an ODE model developed by minimizing ESQE, a procedure that we called single-objective optimization (Figure 1, SOO).
The paragraphs below describe in detail how the different objectives were computed.
Construction of a time-delay correlation matrix (objective DSp)
To test the potential of MOO to combine different network inference approaches we choose to build an objective based on time-delayed Spearman Rank-correlation  (Figure 1, objective DSp). DSp was computed as follows: For each gene pair (i,j), the expression profile of gene i is shifted along the time axis with respect to that of gene j. The Spearman Rank-correlation coefficient is calculated for varying time delays and the largest coefficient from this list forms the (i,j)th element of the delayed Spearman Rank-correlation matrix d-SRC. We also construct a time delay matrix dt from the corresponding values. The objective DSp is then obtained from d-SRC by equating all d-SRC(i,j) = 0 for which dt(i,j) < t o , so that only gene pairs with delay of to or more are considered.
Construction of a gene KO matrix: a ratios-based procedure
The objectives Tr and Sr (Figure 1) were constructed by computing the ratios between the expressions of each gene i in the mutant j and the expression of gene i in the wild type. The expression of gene i is taken either at a given time point tp after inactivation (Tr) or at the steady state (Sr). We selected tp as the time point where the largest numbers of genes have the highest derivative in absolute value. We found that this heuristic rule allowed us to identify a value of tp, which often (8 out of 9 networks tested) corresponded to the highest AUC values within a tolerance of 5% (Figure S1).
Construction of a gene KO matrix: a correlation-based procedure
The objectives Tc and Sc (Figure 1) were computed by calculating the correlation between the expressions of every pair of genes (gene i, gene j) across all KO samples. Similarly to the ratio procedure, the Tc matrix was built using the measure of gene i expression at time t p , where t p was chosen as detailed above.
Combining MOO procedures using an ensemble approach
The ratio and correlation methods were integrated to produce a single model by using an ensemble approach. Within this procedure, a GRM w a was constructed so that |wa(i,j)| = |wr (i,j)| and sign(wa(i,j)) = sign(wc (i,j); Where w r and w c represent two GRMs obtained from the ratio and correlation procedures, respectively. As exemplified in Figure 1, MOO-Sens represent the result of combining the MOO-Sr and MOO-Sc procedures whereas MOO-Tens, is the result of combining the MOO-Tr and MOO-Tc procedures.
The validation study has been performed using the java application GeneNetWeaver (GNW) http://gnw.sourceforge.net. This network generator has been used as part of the DREAM Initiative . It builds synthetic networks by specifying a biologically relevant topology and implementing an ODE model to generate synthetic data. GNW grows the initial topology from a seed node (selected randomly) in a Source Gene Network (E. Coli in this application) by progressively adding a randomly selected neighbouring node till the desired size is reached. Each model can be used to generate simulated time course gene expression data either with the intact network or following deletion of one of the nodes.
We tested the performance of MOO in conjunction with the objectives D-Sp (MOO-Sp), Tc (MOO-Tc), Tr (MOO-Tr), Sc (MOO-Sc) and Sr (MOO-Sr). Each of these procedures was applied to ten independent networks of size 20, 35 and 50 genes. The gene KOs dataset associated with every network was build by generating synthetic data after the stepwise deletion of each gene in the network.
All GNW-generated network-models were used to simulate time-series datasets (26 time points, t_max = 200) as well as steady-state data for all KOs.
Data processing and optimisation procedure
Noise was added to the simulated data (5% of the signal) to mimic variability in experimental replication. Polynomial fitting was used for interpolation of the time-series data  after averaging three experimental replicates. 200 interpolated, equally spaced time-points were then computed and used as input of the MOO procedures. Optimisation of the matrix w was initiated from a randomly generated matrix. In order to test the reproducibility of the optimization methods, fifty independent runs of optimization from each MOO procedure were carried out for a subset of the GNW networks. We found that the AUCs values were always within 0.2%.
Network inference accuracy
In order to evaluate various MOO methods we compared the inferred gene-regulatory matrix w with the true network topologies. The accuracy of the inference process for undirected networks was quantified by using the area under curve (AUC) of a ROC plot (False Positive Rate (FPR) versus True Positive Rate (TPA)). For direct-signed networks the AUC was computed by plotting TNR (True Negative Rate) versus TPR as described in . The distribution of AUC values for boxplots and these represented each batch of networks were compared when appropriate using a Wilcoxon's non-parametric rank sum test .
Modelling in vivo tumour development
In order to assess the potential of NIMOO to model true biological systems we have used two microarray datasets generated in our laboratory.
We first used an in vivo model of glioblastoma development to test the MOO-Sp procedure. In this experimental model  U87 human glioma cells (ATCC, USA) were maintained in DMEM with 10% FBS, antibiotics, and l-glutamine. Fertilized chicken eggs (Gallus gallus; E.A.R.L. Morizeau, Dangers, France) were incubated at 37°C and 80% humidified atmosphere. On day 4 of development, a window was made in the eggshell after punctuating the air chamber and sealed with Durapore tape. On embryonic day 10, a plastic ring was placed on the embryo chorioallantoic membrane (CAM), and 3 million to 5 million U87 cells in 20 μl of medium were deposited after gentle laceration of surface. Implantation experiments were performed in triplicate, and, from day 11 to day 15, mRNA from the growing tumour was extracted every 12 hours using the standard protocol provided in the Qiagen RNeasy kit. Labelling was performed using protocol V5.7 provided in Agilent's Quick Amp One-Colour labelling kit and hybridized onto human Agilent microarrays (AMADID:014850). Data were normalized using quantile normalization and genes differentially expressed over time were identified using the statistical methodology SAM . 58 genes were selected among the most differentially expressed across the time course (Table S3) and used as input of the modelling procedure.
Modelling E. coli acid stress
In order to fully test the potential of MOO methodology we have applied the MOO-Sens procedure to model the E. coli response to mild acid conditions, a stress signal relevant to pathogenesis in diarrheagenic E. coli strains . The procedure was used to integrate two microarray datasets representing the dynamic response of the E. coli MG1655 strain to acid exposure (pH = 5.5) and a gene KO experiment performed in the related strain E. coli BW 25113, representing the transcriptional state of strains mutated in the 26 most differentially expressed genes. In this analysis we aim to reverse engineer the interactions between these 26 genes. The time-course analysis of the response of the E. coli strain MG1655 to acid exposure was performed maintaining a constant cell number (OD600 nm = 2) using a media replenishment procedure. Samples were collected every 5 minutes for 1 hour in three replicated experiments. Mutant strains representing 26 of the most differentially regulated genes over time were selected from the BW25113 KEIO mutant collection  and analysed using microarrays as described below. Experiments were performed exactly in the same conditions as the MG1655 strain but only control and 15 minutes in acid were processed for microarray analysis.
Microarray analysis was performed as follows. 10 ml of cultures were samples at the different time points and stabilized by adding a solution of phenol-ethanol (final concentration of 19% phenol and 1% ethanol). Cell pellets were recovered by centrifugation and stored at -80°C. mRNA was extracted using the standard protocol provided in the Quiagen RNEasy kit (QUIAGEN, USA) and labelled with Cy5 labelled dCTP (Amersham Biosciences, USA) using the CyScribe Post-Labelling Kit (Amersham Biosciences, USA). Probes were then purified using CyScribe purification Kit (Amersham Biosciences, USA) according to the manufacturer's instructions. Labelled probes (80 pmol) were then hybridized on Operon E. coli Ultra GAPS microarray slides (Corning, USA). Slides were washed in AdvaWash automated washing station (Adavlytix, USA) and scanned with the ScanArray® GX (PerkinElmer®, USA), using the ScanArray® software. Data were normalized using quantile normalization and genes differentially expressed over time were identified using the statistical methodology SAM . We modelled the E. coli datasets by using the ensemble approach integrating both correlation and ratio procedures as described above. In order to generate comparable sparse networks we thresholded the connectivity matrix w to obtain networks with same number of genes in the networks (25).
Method implementation and datasets availability
Combining different inference methodologies within the MOO framework improves the accuracy of network reconstruction
The first scenario we considered involved combining two network inference methods to model replicated time course experiments. To achieve this, we used delayed Spearman Rank-correlation  to build the objective Oij (Equation 3) for MOO.
Accuracy of GRN inference with MOO-dSp
Overall, we can conclude that in the event that only replicated time-course experiments are available, a situation which is not uncommon, the integration between two methodologies can lead to a dynamical model with better accuracy than one solely based on a SOO procedure.
Integrating time-course and gene inactivation experiments within the MOO framework: A ratio-based procedure
Having shown that MOO is an effective approach to combine different network-inference methodologies we set to test whether it may also provide a solution to integrate time-course and gene-inactivation experiments.
Accuracy of GRN inference by integrating gene KO datasets in the MOO framework
Ratio methods MOO-Sr/MOO-Tr
Corr. methods MOO-Sc/MOO-Tc
Ensemble MOO Sens/MOO Tens
Integrating time course and gene inactivation experiments within the MOO framework: A correlation-based procedure
In this section, we describe the results of the correlation-based procedure to construct MOO objectives from mutant gene expression data. As detailed in the methods section, this approach works by computing the correlation between the expression profiles of every pair of genes across the mutant samples.
Accuracy of GRN inference with partial coverage gene KO datasets
Combining ratio and correlation-based procedures further improve inference acuracy
Since we have shown that correlation and ratio-based methods provide complementary information, we decided to test whether combining them using an ensemble approach could result in an even higher accuracy of the network inference process.
Overall, MOO-Sens was the best performing procedure in inferring directed-signed networks. Therefore, we concluded that if both time-course and KO data are available for a sub-set of genes of interest, MOO-Sens may be the procedure of choice.
Modelling biological systems with NIMOO
In order to test the validity of NIMOO to model real biological systems, we have analysed two datasets generated in our own laboratory. The first was a replicated gene-expression-profiling time-course experiment representing a model of in vivo glioblastoma development. A sub-set of these data were modelled with the MOO-Sp procedure. The second dataset included a time course representing the transcriptional response of E. coli during acid adaptation and the expression profiling of a compendium of 26 mutants exposed to acid. Because of the availability of both time-course and mutant steady-state data we applied the MOO-Sens procedure.
Modelling in vivo tumour development
The other network hubs are also known important signalling factors. ERBB2 is a gene that encode for a member of the epidermal growth factor (EGF) receptor family of receptor tyrosine kinases. This protein has no ligand-binding domain but it does bind tightly to other ligand-bound EGF receptor family members enhancing kinase-mediated activation of downstream signaling pathways. HSPB1has a cytoprotective function and support of cell survival under stress conditions. This gene is also involved in the apoptotic signalling pathway and interacts with actin and intermediate filaments to protect actin filaments from fragmentation. It also preserves the focal contacts fixed at the cell membrane. These observations support the hypothesis that Nrf2 sits high in the hierarchy of events leading to the development of a fully vascularized tumour.
Reverse engineering an E coli acid response network
Both single objective (SOO) and multiobjective (MOO) optimization methods were applied to investigate regulatory networks representative of E. coli acid adaptation. In order to test the full potential of the NIMOO methodology, we used both time-course and gene-inactivation experiments.
The networks identified using either the time course data (SOO procedure) or the combination of time course and gene KO profiles (MOO procedure) are represented in Figure 7C and 7F respectively. In order to validate the model, we have compared our results with the gene interactions known in literature or extracted from the REGULON DB database .
The SOO method identified a number of gene connections that were known to play a role in acid adaptation. These included the interaction between two of the glutamate-dependent acid-stress response genes gadW and gadX . However, in this model the directions of the gene interactions are mostly incorrect and not representative of the known E. coli acid response mechanisms. For example, the coding glutamate decarboxylase gene gadB is unlikely to be involved in the modulation of the two-component system PhoP/PhoQ.
On the contrary, the gene regulatory network derived from the application of the MOO procedure (Figure 7D) includes several gene interactions known to be important in acid adaptation.
A key interaction involved the two-component system PhoP/PhoQ . This complex is a known upstream regulator of acid adaptation. The model we developed (Figure 7C) reflects the upstream regulatory role of this complex and correctly predicted its influence in the regulation of the acid resistance genes gadW and hdeA . The network also shows the known negative interaction between gadX and gadW  and the inhibition of the crp gene by fis [35, 36]. Another validated interaction found by the MOO procedure is represented by the link between the histone-like protein hns and cadA . Our model shows that hns may activate the expression of cadA. The connection is consistent with the literature, however, in GNB7145K hns mutants Shi et al.  have shown that hns acts as a repressor.
Some of the interactions in the network represent potentially novel regulatory mechanisms in E. coli response to acid stress. An example is the hypothesis that phoP may be involved in the activation of narX, a nitrite/nitrate sensor protein. This is a gene that is part of a two-component system regulating a component of anaerobic metabolism, which is a function known to be regulated during acid response .
In this paper we presented the gene regulatory network inference method "Network Inference with Multi Objective Optimization" (NIMOO).
When tested on simulated and "real world" datasets, NIMOO performs well even if incomplete data are available. The main feature of this methodology is that it can be used to develop dynamical models of gene regulatory networks integrating multiple data sources and combining any existing network inference methodology to identifying the network topology.
Although other methods have the potential to include prior knowledge in the inference process the ability to plug-in different inference methods in the same modelling procedure is so far a unique feature of NIMOO. In this paper we tested this concept and demonstrated that the approach can be successful even if a relatively simple procedure is integrated in the ODE model parameter estimation. However, a more comprehensive testing may be required to explore the full potential of this approach, for example combining more advanced methods in the MOO optimization procedure.
In terms of data integration, we have mainly focused on gene KO experiments. However, some of the procedures we have tested (e.g. MOO-Tc and MOO-Sc) are directly applicable to other types of experimental data. For example, a compendium of environmental and growth factor-induced perturbations could be employed to develop an objective compatible with these procedures. Such objectives could be for example computed by using the information theoretical approach ARACNE .
Moreover, additional information, for instance the confidence level in transcription factor binding consensus sequences in a gene's promoter region could also be incorporated within the optimisation process. More generally, in the event that multiple objectives are used within a MOO procedure, each function's relative importance could be weighted by adjusting the optimization parameters, such as weights θk (Equation 4). Additionally, any definite qualitative knowledge of the presence or absence of gene connections may be used as a constraint on the inferred gene-regulatory matrix (hard prior).
Because of the ability to integrate different methods the user can very easily customize NIMOO. In this respect, NIMOO is an integration framework rather than a specific method. Comparing its performance with existing methods is therefore not necessarily consequential. However, we have performed an initial assessment comparing some implementations of NIMOO to other methods. For example, all NIMOO procedures outperformed TSNI  in inferring undirected networks (Table S1 in Additional File 1) and the MOO-Sens and MOO-Tens performed better with both undirected and direct-signed networks (Table S1 in Additional File 1). Moreover, NIMOO performed in a comparable manner to the method developed by Yip et al. , which won the DREAM3 competition http://wiki.c2b2.columbia.edu/dream/index.php/ (Table S2 in Additional File 1).
So far the application of multi-objective optimization methods to inferring biological networks has been limited: For example, van Someren et al.  and Fomekong-Nanfack et al.  used MOO to incorporate multiple constraints arising from the requirement of stability and robustness of gene networks, and, Liu and Wang  have used MOO to infer biochemical networks by simultaneously minimizing for the concentration error and the slope error. However, in all these cases a single data set and a single reverse engineering criterion were used.
The network-inference framework NIMOO is flexible and can be used in many different scenarios, even when available information is incomplete. The application of NIMOO to biological datasets representing two different "real world" scenarios produced very interesting results. The analysis of the experimental datasets illustrated that inclusion of additional objectives from the same dataset could significantly improve our ability to identify key regulators of relevant biological processes.
Acknowledgements and Funding
The work described in this paper was funded by the CRUK grant C8504/A9488 and partially funded by the BBSRC grant BBC5151041. AS is a recipient of a Darwin Trust PhD fellowship and PA is a recipient of a BBSRC PhD studentship. We would like to thank Victor Trevino for useful comments on the manuscript.
- Fodor SP, Rava RP, Huang XC, Pease AC, Holmes CP, Adams CL: Multiplexed biochemical assays with biological chips. Nature. 1993, 364: 555-556. 10.1038/364555a0View ArticlePubMedGoogle Scholar
- Schena M, Shalon D, Davis RW, Brown PO: Quantitative monitoring of gene expression patterns with a complementary DNA microarray. Science. 1995, 270: 467-470. 10.1126/science.270.5235.467View ArticlePubMedGoogle Scholar
- Stolovitzky G, Monroe D, Califano A: Dialogue on reverse-engineering assessment and methods: the DREAM of high-throughput pathway inference. Ann N Y Acad Sci. 2007, 1115: 1-22. 10.1196/annals.1407.021View ArticlePubMedGoogle Scholar
- Yu J, Smith VA, Wang PP, Hartemink AJ, Jarvis ED: Advances to Bayesian network inference for generating causal networks from observational biological data. Bioinformatics. 2004, 20: 3594-3603. 10.1093/bioinformatics/bth448View ArticlePubMedGoogle Scholar
- Basso K, Margolin AA, Stolovitzky G, Klein U, Dalla-Favera R, Califano A: Reverse engineering of regulatory networks in human B cells. Nature Genetics. 2005, 37: 382-390. 10.1038/ng1532View ArticlePubMedGoogle Scholar
- Schmitt WA, Raab RM, Stephanopoulos G: Elucidation of gene interaction networks through time-lagged correlation analysis of transcriptional data. Genome Research. 2004, 14: 1654-63. 10.1101/gr.2439804PubMed CentralView ArticlePubMedGoogle Scholar
- Rangel C, Angus J, Ghahramani Z, Lioumi M, Sotheran E, Gaiba A, Wild DL, Falciani F: Modeling T-cell activation using gene expression profiling and state-space models. Bioinformatics. 2004, 20: 1361-72. 10.1093/bioinformatics/bth093View ArticlePubMedGoogle Scholar
- Hirose O, Yoshida R, Imoto S, Yamaguchi R, Higuchi T, Charnock-Jones DS, Print C, Miyano S: Statistical inference of transcriptional module-based gene networks from time course gene expression profiles by using state space models. Bioinformatics. 2008, 24: 932-942. 10.1093/bioinformatics/btm639View ArticlePubMedGoogle Scholar
- Guthke R, Möller U, Hoffmann M, Thies F, Töpfer S: Dynamic network reconstruction from gene expression data applied to immune response during bacterial infection. Bioinformatics. 2005, 21: 1626-1634. 10.1093/bioinformatics/bti226View ArticlePubMedGoogle Scholar
- Bansal M, Della Gatta G, 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. Molecular System Biology. 2007, 3: 78-View ArticleGoogle Scholar
- Hecker M, Lambeck S, Toepfer S, van Someren E, Guthke R: Gene regulatory network inference: Data integration in dynamic models: A review. Biosystems. 2009, 96: 86-103. 10.1016/j.biosystems.2008.12.004View ArticlePubMedGoogle Scholar
- Marbach D, Prill RJ, Schaffter T, Mattiussia C, Floreanoa D, Stolovitzkyc G: Revealing strengths and weaknesses of methods for gene network inference. PNAS. 2010, 107 (14): 6286-91. 10.1073/pnas.0913357107PubMed CentralView ArticlePubMedGoogle Scholar
- Yip KY, Alexander RP, Yan K-K, Gerstein M: Improved Reconstruction of In Silico Gene Regulatory Networks by Integrating Knockout and Perturbation Data. PLoS ONE. 2010, 5 (1): e8121- 10.1371/journal.pone.0008121PubMed CentralView ArticlePubMedGoogle Scholar
- Coleman TF, Li Y: On the Convergence of Reflective Newton Methods for Large- Scale Nonlinear Minimization Subject to Bounds. Mathematical Programming. 1994, 67 (2): 189-224. 10.1007/BF01582221.View ArticleGoogle Scholar
- Coleman TF, Li Y: An Interior, Trust Region Approach for Nonlinear Minimization Subject to Bounds. SIAM Journal on Optimization. 1996, 6: 418-445. 10.1137/0806023.View ArticleGoogle Scholar
- Gembicki FW: Vector Optimization for Control with Performance and Parameter Sensitivity Indices. Ph.D. Dissertation. 1974, Case Western Reserve Univ., Cleveland, OH,Google Scholar
- Gembicki F, Haimes Y: Approach to Performance and Sensitivity Multiobjective Optimization: The Goal Attainment Method. IEEE Transactions on Automatic Control. 1975, 20: 769-10.1109/TAC.1975.1101105.View ArticleGoogle Scholar
- Marbach D, Schaffter T, Mattiussi C, Floreano D: Generating Realistic In Silico Gene Networks for Performance Assessment of Reverse Engineering Methods. Journal of Computational Biology. 2009, 16: 229- 10.1089/cmb.2008.09TTView ArticlePubMedGoogle Scholar
- Prill RJ, Marbach D, Saez-Rodriguez J, Sorger PK, Alexopoulos LG, Xue X, Clarke ND, Altan-Bonnet G, Stolovitzky G: Towards a rigorous assessment of systems biology models: the DREAM3 challenges. PLoS One. 2010, 5 (2): e9202- 10.1371/journal.pone.0009202PubMed CentralView ArticlePubMedGoogle Scholar
- Atkinson Kendell: An Introduction to Numerical Analysis. 1988, Chapter 3: John Wiley and Sons, 2,Google Scholar
- Wilcoxon F: Individual comparisons by ranking methods. Biometrics. 1945, 1: 80-83. 10.2307/3001968.View ArticleGoogle Scholar
- Hagedorn M, Javerzat S, Gilges D, Meyre A, de Lafarge B, Eichmann A, Bikfalvi A: Accessing key steps of human tumor progression in vivo by using an avian embryo model. PNAS. 2005, 102 (5): 1643-1648. 10.1073/pnas.0408622102PubMed CentralView ArticlePubMedGoogle Scholar
- Tusher VG, Tibshirani R, Chu G: Significance analysis of microarrays applied to the ionizing radiation response. PNAS. 2001, 98 (9): 5116-21. 10.1073/pnas.091062498PubMed CentralView ArticlePubMedGoogle Scholar
- Foster JW: Escherichia coli acid resistance: tales of an amateur acidophile. Nat Rev Microbiol. 2004, 2 (11): 898-907. 10.1038/nrmicro1021View ArticlePubMedGoogle Scholar
- Baba T, Ara T, Hasegawa M, Takai Y, Okumura Y, Baba M, Datsenko KA, Tomita M, Wanner BL, Mori H: Construction of Escherichia coli K-12 in-frame, single-gene knockout mutants: the Keio collection. Mol Syst Biol. 2006, 2: 10.1038/msb4100050. 2006 0008,Google Scholar
- MATLAB version 7.10.0 (R2010a). The MathWorks Inc. Natick, Massachusetts,
- Kensler TW, Wakabayashi N, Biswal S: Cell survival responses to environmental stresses via the Keap1-Nrf2-ARE pathway. Annu Rev Pharmacol Toxicol. 2007, 47: 89-116. 10.1146/annurev.pharmtox.46.120604.141046View ArticlePubMedGoogle Scholar
- Giudice A, Arra C, Turco MC: Review of molecular mechanisms involved in the activation of the Nrf2-ARE signaling pathway by chemopreventive agents. Methods Mol Biol. 2010, 647: 37-74. full_text full_textView ArticlePubMedGoogle Scholar
- Kensler TW, Wakabayashi N: Nrf2: friend or foe for chemoprevention?. Carcinogenesis. 31 (1): 90-99.PubMed CentralView ArticlePubMedGoogle Scholar
- Gama-Castro S, Salgado H, Peralta-Gil M, Santos-Zavaleta A, Muñiz-Rascado L, Solano- Lira H, Jimenez-Jacinto V, Weiss V, García-Sotelo JS, López-Fuentes A, Porrón-Sotelo L, Alquicira-Hernández S, Medina-Rivera A, Martínez-Flores I, Alquicira-Hernández K, Martínez-Adame R, Bonavides-Martínez C, Miranda-Ríos J, Huerta AM, Mendoza-Vargas A, Collado-Torres L, Taboada B, Vega-Alvarado L, Olvera M, Olvera L, Grande R, Morett E, Collado-Vides J: RegulonDB version 7: transcriptional regulation of Escherichia coli K-12 integrated within genetic response units (Gensor Units). Nuclear Acid Research. 2010, D98-105. 9 Database,Google Scholar
- Ma Z, Richard H, Tucker DL, Conway T, Foster JW: Collaborative regulation of Escherichia coli glutamate-dependent acid resistance by two AraC-like regulators, GadX and GadW (YhiW). J Bacteriol. 2002, 184 (24): 7001-12. 10.1128/JB.184.24.7001-7012.2002PubMed CentralView ArticlePubMedGoogle Scholar
- Groisman EA: The pleiotropic two-component regulatory system PhoP-PhoQ. J Bacteriol. 2001, 183 (6): 1835-42. 10.1128/JB.183.6.1835-1842.2001PubMed CentralView ArticlePubMedGoogle Scholar
- Eguchi Y, Itou J, Yamane M, Demizu R, Yamato F, Okada A, Mori H, Kato A, Utsumi R: B1500, a small membrane protein, connects the two-component systems EvgS/EvgA and PhoQ/PhoP in Escherichia coli. PNAS. 2007, 104 (47): 18712-7. 10.1073/pnas.0705768104PubMed CentralView ArticlePubMedGoogle Scholar
- Kolb A, Igarashi K, Ishihama A, Lavigne M, Buckle M, Buc H: E. coli RNA polymerase, deleted in the C-terminal part of its alpha-subunit, interacts differently with the cAMP-CRP complex at the lacP1 and at the galP1 promoter. Nucleic Acids Res. 1993, 21 (2): 319-26. 10.1093/nar/21.2.319PubMed CentralView ArticlePubMedGoogle Scholar
- González-Gil G, Kahmann R, Muskhelishvili G: Regulation of crp transcription by oscillation between distinct nucleoprotein complexes. EMBO J. 1998, 17 (10): 2877-85.PubMed CentralView ArticlePubMedGoogle Scholar
- Shi X, Waasdorp BC, Bennett GN: Modulation of acid-induced amino acid decarboxylase gene expression by hns in Escherichia coli. J Bacteriol. 1993, 175 (4): 1182-6.PubMed CentralPubMedGoogle Scholar
- Stewart V: Nitrate regulation of anaerobic respiratory gene expression in Escherichia coli. Molecular Microbiology. 1993, 9: 425-434. 10.1111/j.1365-2958.1993.tb01704.xView ArticlePubMedGoogle Scholar
- Van Someren EP, Wessels LFA, Backer E, Reinders MJT: Multi-criterion optimization for genetic network modeling. Signal Processing. 2003, 83: 763-775. 10.1016/S0165-1684(02)00473-5.View ArticleGoogle Scholar
- Fomekong-Nanfack Y, Postma M, Kaandorp JA: Inferring Drosophila gap gene regulatory network: a parameter sensitivity and perturbation analysis. BMC Systems Biology. 2009, 3: 94- 10.1186/1752-0509-3-94PubMed CentralView ArticlePubMedGoogle Scholar
- Liu , Wang : Inference of biochemical models in S-system using multi-objective optimization approach. Bioinformatics. 2008, 24: 1085-1092. 10.1093/bioinformatics/btn075View ArticlePubMedGoogle Scholar