1. A method, comprising executing computer executable instructions by one or more processors for predicting gene expression profile changes upon inhibiting proteins or genes of drug targets on treating a disease, the method further comprising: constructing a genetic network comprising a structure of a dynamic Bayesian network based at least in part on knowledge of drug inhibiting effects on a disease; associating a set of parameters with the constructed dynamic Bayesian network; determining values of a joint probability distribution of the dynamic Bayesian network via an automatic procedure; deriving a mean dynamic Bayesian network with one or more averaged parameters based at least in part on the joint probability values; and calculating a quantitative prediction based at least in part on the mean dynamic Bayesian network, wherein said method searches for an optimal combination of drug targets whose perturbed gene expression profiles are most similar to healthy cells.