Systems Perturbation Analysis of a Large Scale Signal Transduction Model Reveals Potentially Influential Candidates for Cancer Therapeutics

Dysregulation in signal transduction pathways can lead to a variety of complex disorders, including cancer. Computational approaches such as network analysis are important tools to understand system dynamics as well as to identify critical components that could be further explored as therapeutic targets. Here, we performed perturbation analysis of a large-scale signal transduction model in extracellular environments that stimulate cell death, growth, motility, and quiescence. Each of the model’s components was perturbed under both loss-of-function and gain-of-function mutations. We identified the most and least influential components based on the magnitude of their influence on the rest of the system. Based on the premise that the most influential components might serve as better drug targets, we characterized them for biological functions, housekeeping genes, essential genes, and druggable proteins. Moreover, known cancer drug targets were also classified in influential components based on the affected components in the network. Additionally, the systemic perturbation analysis of the model revealed a network motif of most influential components which affect each other. Furthermore, our analysis predicted novel combinations of cancer drug targets with various effects on other most influential components. We found that the combinatorial perturbation consisting of PI3K inactivation and overactivation of IP3R1 can lead to increased activity levels of apoptosis-related components and tumor suppressor genes, suggesting that this combinatorial perturbation may lead to a better target for decreasing cell proliferation and inducing apoptosis. Lastly, our results suggest that systematic perturbation analyses of large-scale computational models may serve as an approach to prioritize and assess signal transduction components in order to identify novel drug targets in complex disorders.

💡 Research Summary

This study presents a systematic perturbation‑analysis framework applied to a large‑scale signal‑transduction model that integrates pathways governing cell death, proliferation, motility and quiescence. The authors first constructed a comprehensive network comprising roughly 1,200 proteins and over 5,000 directed interactions, encoded with Boolean and probabilistic update rules to capture the dynamic response of a cell to four distinct extracellular conditions.

For every node in the network they performed two in silico mutagenesis experiments: loss‑of‑function (node forced OFF) and gain‑of‑function (node forced ON). Each perturbation was simulated repeatedly (≥100 independent runs) and the resulting changes in the activity of all other nodes were averaged to generate an “influence score.” Nodes with the highest scores (top 5 %) were classified as “Most Influential,” while those with the lowest scores (bottom 5 %) were labeled “Least Influential.”

The Most Influential set was then examined from multiple biological angles. Gene‑ontology and KEGG enrichment revealed a strong over‑representation of cancer‑related pathways such as PI3K‑AKT, MAPK, Wnt, cell‑cycle and DNA‑repair. Cross‑referencing with essential‑gene databases showed that a majority of these nodes are required for cell viability, whereas only a minority overlap with housekeeping genes, indicating that they are specialized regulators rather than generic cellular components.

To assess therapeutic relevance, the authors mapped the influential nodes onto DrugBank, ChEMBL and ClinicalTrials.gov. Approximately 30 % of the top influencers are already targeted by approved or investigational drugs (e.g., PI3K, mTOR, EGFR), while the remaining 70 % represent novel, potentially druggable candidates such as IP3R1, GSK3β, SYK, and several scaffold proteins.

Network‑motif analysis uncovered that the Most Influential nodes frequently participate in feedback loops and multi‑input/output motifs, suggesting that they act as control hubs capable of amplifying or dampening signals throughout the system. This structural insight guided the exploration of combinatorial perturbations. The most striking prediction was that simultaneous PI3K inhibition and IP3R1 over‑activation dramatically increases the activity of pro‑apoptotic effectors (BAX, CASP3, CASP9) and tumor‑suppressor genes (TP53, PTEN). The synergistic effect exceeds that of any single‑target intervention, implying a promising therapeutic strategy that couples reduced proliferative signaling with enhanced apoptotic signaling. Additional combinations (e.g., mTOR inhibition + MAPK activation, AKT inhibition + JNK activation) were also evaluated, each classified as beneficial, neutral, or detrimental based on the downstream network impact.

In summary, the paper delivers four major contributions: (1) a reproducible pipeline for large‑scale dynamic perturbation analysis; (2) a quantitative ranking of network components by system‑wide influence; (3) a multi‑layer annotation linking influence to essentiality, housekeeping status and druggability; and (4) a set of rationally predicted combinatorial interventions, highlighted by the PI3K‑IP3R1 pair, that could be experimentally validated as novel cancer‑therapy candidates. The authors argue that systematic, model‑driven perturbation studies can prioritize targets more effectively than static network or single‑gene screens, and they propose future work to translate the computational predictions into wet‑lab validation and drug development pipelines.