Computational modeling with forward and reverse engineering links signaling network and genomic regulatory responses:...

Signal transduction is the major mechanism through which cells transmit external stimuli to evoke intracellular biochemical responses. Diverse cellular stimuli create a wide variety of transcription factor activities through signal transduction pathways, resulting in different gene expression patterns. Understanding the relationship between external stimuli and the corresponding cellular responses, as well as the subsequent effects on downstream genes, is a major challenge in systems biology. Thus, a systematic approach is needed to integrate experimental data and theoretical hypotheses to identify the physiological consequences of environmental stimuli.

We proposed a systematic approach that combines forward and reverse engineering to link the signal transduction cascade with the gene responses. To demonstrate the feasibility of our strategy, we focused on linking the NF-kappaB signaling pathway with the inflammatory gene regulatory responses because NF-kappaB has long been recognized to play a crucial role in inflammation. We first utilized forward engineering (Hybrid Functional Petri Nets) to construct the NF-kappaB signaling pathway and reverse engineering (Network Components Analysis) to build a gene regulatory network (GRN). Then, we demonstrated that the corresponding IKK profiles can be identified in the GRN and are consistent with the experimental validation of the IKK kinase assay. We found that the time-lapse gene expression of several cytokines and chemokines (TNF-alpha, IL-1, IL-6, CXCL1, CXCL2 and CCL3) is concordant with the NF-kappaB activity profile, and these genes have stronger influence strength within the GRN. Such regulatory effects have highlighted the crucial roles of NF-kappaB signaling in the acute inflammatory response and enhance our understanding of the systemic inflammatory response syndrome.

We successfully identified and distinguished the corresponding signaling profiles among three microarray datasets with different stimuli strengths. In our model, the crucial genes of the NF-kappaB regulatory network were also identified to reflect the biological consequences of inflammation. With the experimental validation, our strategy is thus an effective solution to decipher cross-talk effects when attempting to integrate new kinetic parameters from other signal transduction pathways. The strategy also provides new insight for systems biology modeling to link any signal transduction pathways with the responses of downstream genes of interest.