Embodiments of the invention encompass methods for modeling of complex systems, which include, but are not limited to gene regulatory networks, biological systems, and the like. Other embodiments of the invention include the development of computational tools for the identification and discovery of potential targets for therapeutic intervention in diseases such as cancer. The embodiments discussed utilize methods that model the potential effect of individual genes on the global dynamical network behavior, both from the view of random gene mutation as well as intervention in order to elicit desired network behavior.