Systems biology efforts are increasingly adopting quantitative, mechanistic modeling to study cellular signal transduction pathways and other networks. However, it is uncertain whether the particular set of kinetic parameter values of the model closely approximates the corresponding biological system. We propose that the parameters be assigned statistical distributions that reflect the degree of uncertainty for a comprehensive simulation analysis. From this analysis, we globally identify the key components and steps in signal transduction networks at a systems level. We investigated a recent mathematical model of interferon gamma induced Janus kinase-signal transducers and activators of transcription (JAK-STAT) signaling pathway by applying multi-parametric sensitivity analysis that is based on simultaneous variation of the parameter values. We find that suppressor of cytokine signaling-1, nuclear phosphatases, cytoplasmic STAT1, and the corresponding reaction steps are sensitive perturbation points of this pathway.