A computer-implemented method for determining an optimal treatment protocol for a disease related to angiogenesis, comprising creating an angiogenesis model including pro-angiogenesis and anti-angiogenesis factors. Effective vessel density (EVD) is incorporated as a factor regulating switching on and switching off of at least one component in the angiogenesis model. Effects of vasculature maturation and mature vessels destabilization are incorporated. Pro-angiogenesis and anti-angiogenesis factors, which can influence changes in state of a tissue are selected. Effects of drugs in the pro-angiogenesis and anti-angiogenesis factors are incorporated. A plurality of treatment protocols in a protocol space is generated. A best treatment protocol based on a pre-determined criteria.