Skip to main content
  • Research article
  • Open access
  • Published:

The tumor as an organ: comprehensive spatial and temporal modeling of the tumor and its microenvironment



Research related to cancer is vast, and continues in earnest in many directions. Due to the complexity of cancer, a better understanding of tumor growth dynamics can be gleaned from a dynamic computational model. We present a comprehensive, fully executable, spatial and temporal 3D computational model of the development of a cancerous tumor together with its environment.


The model was created using Statecharts, which were then connected to an interactive animation front-end that we developed especially for this work, making it possible to visualize on the fly the on-going events of the system’s execution, as well as the effect of various input parameters. We were thus able to gain a better understanding of, e.g., how different amounts or thresholds of oxygen and VEGF (vascular endothelial growth factor) affect the progression of the tumor. We found that the tumor has a critical turning point, where it either dies or recovers. If minimum conditions are met at that time, it eventually develops into a full, active, growing tumor, regardless of the actual amount; otherwise it dies.


This brings us to the conclusion that the tumor is in fact a very robust system: changing initial values of VEGF and oxygen can increase the time it takes to become fully developed, but will not necessarily completely eliminate it.


Cancer, a state of abnormal growth and regulation of cells, which proliferate in an uncontrolled way, is among the leading causes of mortality worldwide. Researching cancer is of high importance.

The dynamic microenvironment in which a tumor originates, also termed the stroma, plays a critical role in tumor initiation and progression, and may be an important factor in developing therapeutic approaches. For years, research has focused on understanding the transformation of normal cells into neoplastic or cancerous ones. However, it has become evident that the surrounding environment of the tumor cells is equally important. Cells that surround the tumor and can take part in its development include innate and adaptive immune cells, fibroblasts [1], cells that line the blood and lymphatic vessels, and the proteins that make up the structural component–the extra cellular matrix (ECM). The tumor cannot survive or progress on its own [2, 3]; it entirely depends on this dynamic microenvironment in which it originates and the bi-directional interactions with this surroundings [4, 5]. These include secretion of signals or cell-cell interactions and can act to either enhance or block tumor formation. The conditions within the tumor’s microenvironment differ considerably from those in normal tissue.

Angiogenesis, the process of new blood vessels growing from pre-existing ones, is an essential step in the transition of tumors from a dormant state to a malignant one. The new blood vessels grow towards the tumor and feed it with large supplies of oxygen and nutrients. Without angiogenesis, tumors cannot grow beyond the size of 1 mm3. The process leading to angiogenesis begins when the tumor cells lack oxygen (a state of hypoxia), and is followed by a series of events, orchestrated by a variety of activators and inhibitors.

There is a balance between the pro- and anti-angiogenic factors and when this balance tips in favor of the pro-angiogenesis, the onset of angiogenesis, or the angiogenic switch, is induced [6]. The tumor cells release (or cause nearby cells to produce) angiogenic factors that stimulate the formation of blood vessels and recruit them to the tumor’s area [7].

A major activator of the proliferation and migration of the endothelial-cell (the vessels lining) is the vascular endothelial growth factor (VEGF). The endothelial cell at the tip of the emerging vessel starts to migrate towards the angiogenic signals by sensing the concentration gradient [8, 9], while in the stem of the newly formed vessel the endothelial cells proliferate, adhere to each other and create a lumen (the inner space of the tubular vessel). Indeed, blocking angiogenesis by inserting anti-angiogenic factors has been a strategy for clinicians in their efforts to arrest tumor growth [10].

Computational approaches have become a big part of biological research. Computational modeling of biological systems provides a means to integrate a large amount of data and generate a comprehensive overview of the behavior of a system as a whole [11, 12]. Furthermore, being able to visualize animations of the model in operation can significantly aid in gaining a clear understanding of complex biological behavior, and together can serve as a beneficial way to analyze the system and make new discoveries.

Due to the complexity of cancer, a better understanding of tumor growth dynamics and insights into the cancer’s behavior can be gleaned from a dynamic computational visual model. We have therefore focused on modeling the cancerous tumor and its microenvironment and bringing it to life with an interactive animation tool.

Cancer, given its clinical importance, has been studied in detail, and is continuously under intensive investigation. A multitude of experimental, clinical and theoretical studies exist and have shed light on many aspects of cancer on all levels: the sub-cellular scale (e.g., DNA and proteins) [13, 14], cellular scale (activation, proliferation, interactions) [15, 16], and system scale (cell migration, diffusion, metastasis) [17, 18]. Given its complexity and multiscale nature, a better understanding of tumor growth dynamics can be expected from a suitable approach to computational modeling [1921].

Extensive attempts have been carried out to model and analyze cancer, or particular facets thereof [22, 23]; for the most part, this is done by traditional mathematical modeling [2430], using a top-down approach, whereby the behavior of the system is inserted into the model. [3141]. Modeling work on cancer using agent based methods has also been used ([42, 43]).

Our goal was to create a comprehensive model of the entire system, whereby we first model and only then ask the questions and not the other way around. This approach, sometimes termed ‘executable biology’ [4449], focuses on designing executable models that mimic complex biological phenomena. It is carried out in a bottom-up fashion, whereby the behavior of each of the elements of the system (e.g., cells) is described individually, and the system’s overall behavior emerges from those of its many elements. The main language used to build our model is the visual formalism of Statecharts.

We connected our Statechart model to SimuLife, an animation tool that we built especially in our group for viewing the behaviors of such biological models, and which thus serves as a sort of front end to them [50]. SimuLife, based on the technique of reactive animation (RA) [51, 52], is 3D, web-based, and easy to use via an intuitive interface. Visualization of the cancer model is important, as it enables one to see the development and morphology of the tumor and its surroundings based on its individual components. It can be used to fine-tune the model, visualize the effect of changing elements or parameter values, and to verify the behavior of the system.

Results and discussion

We created a comprehensive, fully executable, spatial and temporal 3D computational model that demonstrates the behavior of a typical cancerous solid tumor together with its microenvironment, treating it somewhat like a developing organ.

The model captures the ongoing bidirectional cross talk between the tumor and its surroundings, which plays a critical role in tumor initiation and progression, so that researching the dynamic behavior and morphology of this system via such a model should be of great interest.

The way we chose to build the model is in a bottom-up fashion, whereby the behavior of each of the system’s elements is described individually, using only the fundamental building blocks of that element. The system’s overall behavior emerges from that of the elements thereof, which is the essence of realistic modeling.

The tumor model in Statecharts

Using Statecharts, a generic program of behavior was created for each of the different types of the objects. During an execution of the model many instances of the objects are generated to represent each specific instance taking on its explicit states accordingly (see Fig. 1 for an example of a statechart). This resulted in a comprehensive and reactive computational model.

Fig. 1
figure 1

Statechart of a cell. An example of a statechart that was used in the model

The tumor and its microenvironment act together as an organ, which defines the ‘world’ of the model. The model begins with a single cancerous cell, which proliferates to gradually form a primary tumor, consisting of cancerous cells that originated from the initial one. Each cell has its specific 3D position within this world and takes on its own behavior, depending on its current state and surroundings.

Controlled by the Statechart driving its behavior, the tumor cell constantly senses its immediate surroundings, and consumes available oxygen at a certain level. If it cannot consume a sufficient amount of oxygen, it will not be able to proliferate, will go into a state of hypoxia and will secrete angiogenic factors (VEGF). The VEGF is secreted in pulses, and continues to be secreted until the cell is able to consume enough oxygen so as not to be in a state of hypoxia. The VEGF diffuses from the specific tumor cell and moves in a random manner. If the cell continues to lack oxygen it enters a state of anoxia, becoming necrotic and no longer active.

A number of blood vessels, made up of individual endothelial cells, are initially located either around the border of the defined world, or at random positions therein. The endothelial cells are of roughly the same size as the tumor cells, and similarly sense their environment constantly. Once they bind to an amount of VEGF above a specified threshold within a certain amount of time, they become activated and begin the process of angiogenesis, where the vessels elongate in a direction that follows the VEGF gradient. This occurs by proliferation of those endothelial cells that have become activated and continues for as long as the endothelial cell continues to bind enough VEGF for its continuation. If the endothelial cell encounters another endothelial cell, it will join it and stop elongating. Activated endothelial cells can also split and branch out of their main vessel if they bind a high amount of VEGF in a short amount of time. Due to the delta-notch inhibition between adjacent cells, once an endothelial cell is activated, the neighboring cells cannot become activated too [53, 54]. The endothelial cells secrete oxygen at a constant rate, and like the VEGF molecule, the oxygen diffuses out of the cell in a random manner and is eliminated when consumed.

In this way, the newly produced blood vessels make their way to the tumor, in the process forming a unique spatial organization, and as a result the tumor continues to grow. Newly produced endothelial cells that do not continue to receive a minimum amount of VEGF cannot survive and eventually cause the death of the vessel of which they are part.

Fibroblast cells [1], which are the main components of the ECM, are initially placed randomly around the tumor. During the model’s execution, those that are close to the tumor cells and have enough oxygen have a greater chance of becoming activated into CAFs (cancer associated fibroblasts) [55]. CAFs secrete VEGF in correlation with their hypoxia state [56], which helps recruit the blood vessels. They also secrete HGF (Hepatocyte Growth Factor), which helps the tumor cells proliferate, as well as degrade the ECM around them, which, in turn, helps tumor cells move. CAFs themselves are motile and move towards the tumor [57, 58] by following the gradient of FGF (Fibroblast Growth Factor) that is secreted by the tumor cells. This type of movement can casue entry of CAFs into the tumor [55].

The SimuLife animation tool

Whilst the Statechart model holds the information of each of the individual objects, SimuLife makes it possible to visualize the information of all the objects together at once [50]. This helps greatly when developing the model; to be able to see the dynamics of the whole system, verify the behavior, fine-tune the model as a consequence, and examine the effect of the parameters on the system.

Results analysis

Many parameters were incorporated into in the model. Since we do not use every biological aspect whilst modeling, the parameter values do not represent exact measurement values but only values that are logical relative to each other. In addition, the values of the parameters in the model are unit-less (for a list of parameters, see Additional file 1). What guides us when adjusting parameter values is the behavior of the element in the model and that this matches the actual behavior of the real biological element it represents.

In order to check the correctness of the whole model and understand the behavior of the system, analysis was performed on the data retrieved from the model, and it was compared to existing biological data.

First, we wanted to verify that the behavior of the system modeled is consistent with that of a real cancerous tumor. This was done by looking at a number of points of comparison:

  1. 1)

    The overall fundamental behavior of a solid tumor, where the malignant cells proliferate into a primary tumor until they lack sufficient oxygen and nutrients, and then send out VEGF in order to recruit new blood vessels, which arrives at the endothelial cells of the original blood vessels [59, 60]. The VEGF molecules will have an affect only after enough molecules have arrived to initiate angiogenesis, a process that can take up to a few days. The vessels, in turn, grow towards the concentration gradient and supply the tumor with oxygen/nutrients, enabling it to continue growing. The vessels arrive at the tumor within anything from a few days to a few weeks [61], creating a new network [6, 10].

    This behavior is found in the model, and when visualized by using Simulife (Fig. 2) compares favorably to images found in [62, 63].

    Fig. 2
    figure 2

    SimuLife image of the tumor. A SimuLife image of the tumor surrounded by vessels and fibroblasts (green) as created by the Statecharts model. Shown from two different angles

    In addition, a necrotic core is formed in the inner part of the tumor where the cells do not receive enough oxygen and undergo necrosis [64, 65] (Figs. 3 and 4). This occurs as the tumor expands and the inner cells are no longer exposed to oxygen. The necrotic core occupies from a few percent of the tumor up to the majority of its volume, depending on the amount of available oxygen.

    Fig. 3
    figure 3

    Slicing of the tumor in SimuLife; the necrotic core (dark blue) can be seen

    Fig. 4
    figure 4

    Distribution of cells. Distribution of the different cells as a function of the distance from the center of the tumor (x = 0) at time step 788. Necrotic cells (light blue) occur at the center of the tumor. Vessel splits (red) occur more as they approach the tumor (blue)

    The blood vessels branch more as they approach the tumor (the brush-border effect, [6668]). This occurs as endothelial cells that are closer to the tumor sense more VEGF and are therefore activated more easily to sprout and form new vessels (Fig. 4).

    The diffusion behavior of the VEGF and oxygen molecules can be verified by their distribution over time. Molecules diffuse in space, forming a symmetric, decreasing gradient [69]. In the model, oxygen is initially present everywhere but is gradually consumed by the tumor and in parallel is secreted by vessels. VEGF is secreted by the tumor (and partly by fibroblasts) and spreads out while being consumed by the endothelial cells (Fig. 5). The molecules diffuse by moving in a random fashion from their point of secretion, gradually occupying available space.

    Fig. 5
    figure 5

    Distribution of molecules. Distribution of oxygen (blue) and VEGF molecules (green) as a function of the distance from the center of the tumor (x = 0) at different time steps throughout the simulation (left column followed by right column, from top to bottom). Oxygen is initially present everywhere but is gradually consumed by the tumor and in parallel secreted by vessels. VEGF is secreted by the tumor and spreads out while being consumed by the endothelial cells. The VEGF in the first two graphs is due to secretion by fibroblasts and not by the tumor

    It is important to note that the above anticipated behaviors of the system (necrotic core, increased branching near the tumor, and distribution of diffused molecules) emerged bottom-up from the behavior of the individual components and rules of the system, and not by inserting such behavior into the model explicitly.

    2) The dynamic behavior of the various cells and molecules as observed from running and analyzing the model, agrees with that found in literature. In (Fig. 6) we see that the emergent behavior of the system results in an initial linear growth of tumor cells, a halt in the growth (due to lack of oxygen), then growth of endothelial cells, and consequently an exponential growth in the number of tumor cells. This can be compared to results found in [7073], where the growth of tumor cells and blood vessels exhibit a similar pattern.

    Fig. 6
    figure 6

    Tumor cell and endothelial cell dynamics in the simulation. The behavior of the tumor cells and vessels in the model can be compared to biological results

  2. 3)

    Further verification included checking that “normal” cases occur in the model. This included testing the system under the following sets of circumstances:

    1. a)

      No angiogenesis, which resulted in a primary tumor that stopped growing at some point due to low oxygen/nutrient supply, and then started to die off [7477] (Fig. 7).

      Fig. 7
      figure 7

      No angiogenesis. Tumor does not develop. Left: tumor cell and endothelial cell dynamics. Right: image of the same run in SimuLife (blue cells are necrotic)

    2. b)

      Non-cancerous and non-proliferating cells, which are able to live on normally without angiogenesis (Fig. 8).

      Fig. 8
      figure 8

      No angiogenesis. Non-cancerous and non-proliferating cells live on. Left: tumor cell and endothelial cell dynamics. Right: image of the same run in SimuLife

    3. c)

      Blood vessels placed far from the tumor cells (at ~300 μm instead of ~200 μm, where each position in the model corresponds to 5 μm), which results in no angiogenesis, and hence not enough oxygen, leading to a non-active tumor [7779] (Fig. 9).

      Fig. 9
      figure 9

      Initial blood vessels are too far away. Tumor does not develop. Left: tumor cell and endothelial cell dynamics. Right: image of the same run in SimuLife (blue cells are necrotic)

The values of the parameters used in the model can be changed easily. Playing with the values of these parameters, as well as with other parts of the model, helped to gain a better understanding of their effect on the behavior of the system.

High vs. low oxygen secretion

Playing with the parameter OxygenSecretionAmount, which represents the number of oxygen units secreted every pulse in the simulation, reveals that very low oxygen secretion slows down the growth of the tumor until it is almost completely eliminated. However, as long as some cells are still active, they manage to recover from the situation and end up growing into a large tumor, as large as when OxygenSecretionAmount is high, but later in time (Fig. 10). Therefore, above a certain threshold of minimum oxygen secretion, the tumor and vessels will eventually develop to their full competence. Below this threshold, a tumor will not develop in the model. Although in reality tumor cells deprived of oxygen can produce energy by fermenting sugars (Warburg effect, [80]), and become even more aggressive, the oxygen in our model represents both actual oxygen and nutrients.

Fig. 10
figure 10

High vs. low oxygen secretion simulations. Top: graph of tumor cells showing that at a low oxygen level the tumor cells drop to almost zero at ~800ts but recover and reach 16,000 cells at ~1100ts, whereas at a high oxygen level they do not demonstrate a drop and reach 16,000 cells shortly after ~800ts. Bottom: images of these runs in SimuLife (left for low oxygen, right for high oxygen), presenting amounts on their left tab. Both images are presented at approximately the same time step – at low oxygen only few active tumor cells are present, at high oxygen the tumor consists of many cells where almost all are active

High vs. low VEGF secretion

VEGF secretion by the tumor cells is what recruits the vessels towards the tumor in order to supply it with oxygen/nutrients. Playing with the parameter VEGFSecretionAmount, the amount of VEGF units secreted by the tumor cells at every pulse in the simulation, reveals that, here too, very low VEGF secretion results in tumor and vessels similar to those resulting from high VEGF, although it takes longer, as vessels take more time to become activated by the VEGF (Fig. 11). This indicates that the tumor system is insensitive to the amount of VEGF, as long as there is some VEGF secreted by the tumor cells. Since there are many secreting cells in a tumor, the VEGF will eventually reach and activate the nearby vessels.

Fig. 11
figure 11

High vs. low VEGF secretion simulations. Top: graph of endothelial cells showing that at low VEGF the endothelial cells are activated at ~600ts and reach ~5000 cells at ~900ts, much later than at high VEGF secretion, which begins at ~400ts and reaches ~5000 cells at ~700ts. Bottom: images of these runs in SimuLife (left for low VEGF, right for high VEGF), presenting amounts on their left tab. Both images are presented at approximately the same time step–at low VEGF angiogenesis has only begun, whereas at high VEGF there are many activated and branched vessels

High vs. low angiogenic switch threshold

In order for the endothelial cell to be activated for angiogenesis, it needs to meet a minimum amount of VEGF in a certain amount of time [76]. Playing with the parameter AngiogenicSwitchThreshold, which represents the minimum amount of VEGF needed to initiate angiogenesis, reveals once again that there is a threshold for this parameter; any value above it means that it will take the vessels too long to become activated, and they won’t arrive at the tumor in time before it dies. Below this value the tumor survives and arrives at the same end point, but the time it takes it takes the tumor to develop depends on the value: the lower the threshold the faster the occurrence of angiogenesis (Fig. 12).

Fig. 12
figure 12

High vs. low angiogenic switch threshold simulations. Graph of tumor cells and endothelial cells at low angiogenic switch threshold (light colors), and high angiogenic switch thresholds (dark colors). They reach approximately the same levels, but at a time shift of ~200ts. Also, at high angiogenic thresholds we see that the tumor almost dies at ~900ts, but then recovers

High vs. low hypoxia level

Hypoxia (insufficient oxygen supply) in the model occurs after a cell has not consumed a sufficient amount of oxygen for the duration of HypoxiaLevel time steps. At this point the cell continues to live but cannot proliferate. If it is a tumor cell, it also begins to secrete VEGF in order to recruit blood vessels. At high levels of HypoxiaLevel the tumor grows more before needing the help of angiogenesis. When the tumor cells do eventually secrete VEGF molecules, there are many cells that do so, hence angiogenesis occurs fast and the tumor continues to grow. When HypoxiaLevel is too low, the tumor cells lose their ability to proliferate very fast and so the VEGF they secrete is not enough to recruit the vessels in time before the tumor dies. Once again, values higher than the minimum threshold will eventually result in a developed tumor, but at different time durations–the lower the value the more time it will take.

High vs. low anoxia level

Anoxia (total depletion in the level of oxygen) refers to when the cell is completely out of oxygen and cannot continue to live. At this stage it will become necrotic (an un-programmed, unnaturally occurring cell death–as opposed to apoptosis). This occurs at a later stage than hypoxia, after the cell has had insufficient oxygen for a duration of AnoxiaLevel time steps. When this level is very high no necrotic core occurs, since although all cells may not be proliferating, they live on. At very low levels the entire tumor dies out very fast. Values in between form necrotic cores of different sizes, depending on the value of the parameter, but since the outer cells continue to receive oxygen, the tumor as a whole continues to develop.

A list of the VEGF and Oxygen parameters that were used for analysis, as well as the range of values thereof that ensures tumor recovery and development can be found in Additional file 1.

In each of the conditions described above, if the tumor recovered from just a few surviving cells, those cells are considered the stronger, more aggressive ones (survival of the fittest [81]).

This analysis indicates that the tumor is in fact a very robust system. On the one hand, each of the key parameters has a threshold, which, if crossed, the tumor does not develop. However, at the same time, any value within the range allows the tumor to eventually become fully developed, even if it takes longer and the tumor has to overcome tough conditions. (A short clip showing these results in SimuLife can be found in Additional file 2.) The issue of robustness of the tumor was raised in the past [82], where it was suggested that this fact calls for new therapeutics.

Additional file 2 Video clip of the cancer model in SimuLife. A recording of the process of the model from Statecharts to the SimuLife animation presenting some of the biological results.

Moreover, what we saw in all runs of the model, especially in those where the tumor recovered from tough conditions, was the phase transition in the tumor cell dynamics. At first, the tumor grows linearly, consuming oxygen from its surroundings. The growth then slows down, and sometimes even decreases, due to low oxygen. When enough blood vessels finally arrive at the tumor and feed it with a large quantity of oxygen, the tumor cells start growing exponentially (Fig. 13). At this point, the amount of cells increases, requiring more oxygen, and the cells therefore secrete more VEGF, and in turn receive more oxygen. In this way, the tumor and its surroundings maintain themselves. This is a critical turning point for the tumor, since if it does not manage to pass this point, it will simply die. This finding supports the need to treat tumors as early as possible [83, 84].

Fig. 13
figure 13

Tumor cell dynamics. Just before 800ts there is a turning point for the tumor and a phase transition


Cancer affects many people. It is researched in numerous labs around the world, in order to understand its function and behavior better and to develop possible treatments. Cancer is a multi-scale and complex system, where the tumor and its microenvironment work together in a way similar to an organ. Applying advanced computational techniques to this system is complementary to classical ongoing research. It enables the integration of a variety of data and helps one see the big picture.

In our research we used the language of Statecharts with the Rhapsody tool in order to create a comprehensive 3D model of a cancerous solid tumor, together with its microenvironment. The model was constructed in a bottom up method, where the behavior of each kind of participating component was built individually, and the system’s overall dynamic behavior emerged from those of the elements thereof.

In order to better understand biological models there is a need for detailed and realistic visualization. For this purpose the SimuLife tool was developed, where an animation is dynamically constructed at real time, producing an interactive visualization of the system’s runs. It shows the tumor cells at their precise 3D locations, together with the blood vessels that consist of the individual endothelial cells. The blood vessels elongate towards the tumor by following the VEGF gradient, and in turn secrete oxygen. Cell proliferation or death (in the case of tumor cells this is necrosis), and molecule movement can be observed too. SimuLife allows one to easily play with the animation, send commands back to the model during runtime and observe the immediate resulting output. A more detailed description of SimuLife and its abilities can be found in [50].

Using the SimuLife tool, we were able to see that the model matched the behavior of a solid tumor in a number of ways: It developed an inner necrotic core, the branching of the blood vessels occurred more often as they approached the tumor, the VEGF molecules that were secreted from the tumor cells diffused and were distributed throughout space until reaching the blood vessels, the oxygen secreted from the endothelial cells was also diffused and finally reached and was consumed by the tumor cells. In addition, the dynamic behavior of the various cells showed patterns similar to those found in real tumors. Further verification was carried out to test known cases within the model. This included verifying that the tumor cannot continue to grow without angiogenesis, that a non-cancerous cell can continue to live on without angiogenesis, and that when placing blood vessels far from the tumor angiogenesis does not take place.

From analyzing the components of the model by looking at their qualitative and quantitative dynamics, playing with the different parameters, and inspecting and analyzing the resulting animations, we conclude that the tumor has a turning point, which depends on thresholds of the key parameters that effect amounts of VEGF and oxygen. At this point, the tumor either dies, or recovers and continues to develop to become a full, actively growing tumor. Thus, the tumor’s growth may be halted or declined while it waits for the blood vessels to deliver oxygen. If the oxygen arrives in time, while there are still surviving tumor cells, a phase transition will occur in the tumor’s growth, and from a linear growth rate it will suddenly start growing exponentially. The VEGF and oxygen parameters (amount secreted, consumption, activation threshold….) are what affect the fate of the tumor in this case. If the thresholds are met and the fate of the tumor is to continue, the values of these parameters do not make a big difference. A tumor that survives will develop similarly within the range of allowed values, the main difference being in the time it takes.

It is known that a tumor cannot continue to grow above a certain size (~1 mm3) without angiogenesis [6, 85, 86]. The angiogenic switch point is when the blood vessels have bound enough VEGF to begin angiogenesis [87] and deliver oxygen and nutrients to the tumor. Evidence has shown that tumor growth is slow and linear before vascularization and rapid and nearly exponential after vascularization [88, 89]. Therefore, the essential need of VEGF to recruit blood vessels, and oxygen and nutrients for the tumor to continue growing is confirmed. Here, we suggest that not only is there an angiogenic switch turning point that causes the tumor to enter exponential growth, but there is also a recovery turning point. This means that: 1) if the tumor cells secrete VEGF in any amount, the blood vessels will eventually arrive at the tumor, and 2) if at least some active tumor cells are alive when the blood vessels arrive at the tumor, and if they deliver a sufficient amount of oxygen for those cells to continue their activity, the tumor will eventually become fully developed, even if it takes longer. This also suggests that the tumor, together with its microenvironment, is a robust system, reaching its maximum outcome if the minimum conditions are met, regardless of the actual amounts. It therefore seems like the tumor does not economize resources, and sends well over the needed angiogenic factors.

These conclusions may provide further evidence as to why inhibiting VEGF, or reducing the supply of oxygen and nutrients to the tumor, does not always result in its complete elimination. This is especially relevant to anti-cancer treatments other than chemotherapy, such as VEGF inhibitors [90], or destruction of the tumor’s surrounding blood vessels and hence elimination of its oxygen supply, such as is done in photodynamic therapy [9195]. These treatments, according to the model’s results, can extend the time it takes the tumor to become fully developed.


The Statecharts modeling language

Our model was designed using the visual language of Statecharts [96, 97], which was invented as a system engineering tool to aid in the design of complex reactive systems [96, 98], where the components react to each other and to the environment. The language makes possible the dynamic and visual specification and execution of reactive behavior via the use of intuitive, yet completely formal and fully executable, diagrams. Statecharts describe discrete behavior using states, and events that cause transitions between the states. Orthogonal/concurrent states may also be specified, such that the system or parts thereof may be in several different states simultaneously, in accordance with the different stages of the simulation. In addition, the language is hierarchical, so that states may contain substates, which enable description at multiple levels, as well as level-rich transitions. The object-oriented version of the Statecharts language [97] is based on an intra-object philosophy; i.e., on supplying the full description of the internal behavior of each of the participating objects.

Statecharts and their execution are supported by several appropriate tools, such as Rhapsody, which was co-designed by the second-listed author, and is available from IBM (

The language has become very useful in modeling a variety of biological systems (see past [4447, 99]), since these are in fact complex reactive systems; they interact with, and respond to, both the environment and other parts within the system [48, 98]. We have used Statecharts to specify the behavior of the individual entities that take part in the cancer process, in order to capture the dynamic behavior and morphology of the system.

SimuLife, a 3D animation tool

SimuLife is an interactive animation tool we have built for visualizing models of cellular biology [50]. It can receive inputs from the model, as well as send information to the model (both sent as XML files) and draws the appropriate animated graphics based on the changes. This makes it possible to see the ongoing events of the entire model on the screen in real time [51, 52]. SimuLife is based on WebGL (Web Graphics Library), a JavaScript API (THREE.js framework in our case) for rendering interactive 3D graphics within any compatible web browser without the use of plug-ins. The client side is Chrome and the communication with external engines is done via sockets. The realistic images used within SimuLife (e.g., of a cell) are in the COLLADA format.



Cancer associated fibroblast


Extra cellular matrix


Fibroblast growth factor


Hepatocyte growth factor


Reactive animation


Vascular endothelial growth factor


  1. Kalluri R, Zeisberg M. Fibroblasts in cancer. Nat Rev Cancer. 2006;6(5):392–401.

    Article  CAS  PubMed  Google Scholar 

  2. Hanahan D, Weinberg RA. The hallmarks of cancer. Cell. 2000;100(1):57–70.

    Article  CAS  PubMed  Google Scholar 

  3. Hanahan D, Weinberg RA. Hallmarks of cancer: the next generation. Cell. 2011;144(5):646–74.

    Article  CAS  PubMed  Google Scholar 

  4. Witz IP, Levy-Nissenbaum O. The tumor microenvironment in the post-PAGET era. Cancer Lett. 2006;242(1):1–10.

    Article  CAS  PubMed  Google Scholar 

  5. Tlsty TD, Coussens LM. Tumor stroma and regulation of cancer development. Annu Rev Pathol. 2006;1:119–50.

    Article  CAS  PubMed  Google Scholar 

  6. Bergers G, Benjamin LE. Tumorigenesis and the angiogenic switch. Nat Rev Cancer. 2003;3(6):401–10.

    Article  CAS  PubMed  Google Scholar 

  7. Ehrbar M, Djonov VG, Schnell C, Tschanz SA, Martiny-Baron G, Schenk U, Wood J, Burri PH, Hubbell JA, Zisch AH. Cell-demanded liberation of VEGF121 from fibrin implants induces local and controlled blood vessel growth. Circ Res. 2004;94(8):1124–32.

    Article  CAS  PubMed  Google Scholar 

  8. Jakobsson L, Franco CA, Bentley K, Collins RT, Ponsioen B, Aspalter IM, Rosewell I, Busse M, Thurston G, Medvinsky A, et al. Endothelial cells dynamically compete for the tip cell position during angiogenic sprouting. Nat Cell Biol. 2010;12(10):943–53.

    Article  CAS  PubMed  Google Scholar 

  9. Benedito R, Roca C, Sorensen I, Adams S, Gossler A, Fruttiger M, Adams RH. The notch ligands Dll4 and Jagged1 have opposing effects on angiogenesis. Cell. 2009;137(6):1124–35.

    Article  CAS  PubMed  Google Scholar 

  10. Carmeliet P, Jain RK. Angiogenesis in cancer and other diseases. Nature. 2000;407(6801):249–57.

    Article  CAS  PubMed  Google Scholar 

  11. Aebersold R, Auffray C, Baney E, Barillot E, Brazma A, Brett C, Brunak S, Butte A, Califano A, Celis J, et al. Report on EU-USA workshop: how systems biology can advance cancer research (27 October 2008). Mol Oncol. 2009;3(1):9–17.

    Article  PubMed  Google Scholar 

  12. Fisher J, Harel D, Henzinger TA. Biology as reactivity. Comm Assoc Comput Mach. 2011;54(10):72–82.

    Google Scholar 

  13. Howk CL, Voller Z, Beck BB, Dai D. Genetic diversity in normal cell populations is the earliest stage of oncogenesis leading to intra-tumor heterogeneity. Front Oncol. 2013;3:61.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Jain HV, Jackson TL. A hybrid model of the role of VEGF binding in endothelial cell migration and capillary formation. Front Oncol. 2013;3:102.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Tang L, van de Ven AL, Guo D, Andasari V, Cristini V, Li KC, Zhou X. Computational modeling of 3D tumor growth and angiogenesis for chemotherapy evaluation. PloS one. 2014;9(1):e83962.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Enderling H, Hlatky L, Hahnfeldt P. Cancer stem cells: a minor cancer subpopulation that redefines global cancer features. Front Oncol. 2013;3:76.

    PubMed  Google Scholar 

  17. Szabo A, Merks RM. Cellular potts modeling of tumor growth, tumor invasion, and tumor evolution. Front Oncol. 2013;3:87.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Cheung KJ, Gabrielson E, Werb Z, Ewald AJ. Collective invasion in breast cancer requires a conserved basal epithelial program. Cell. 2013;155(7):1639–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Materi W, Wishart DS. Computational systems biology in cancer: modeling methods and applications. Gene Regul Syst Bio. 2007;1:91–110.

    PubMed  PubMed Central  Google Scholar 

  20. Wang Z, Butner JD, Kerketta R, Cristini V, Deisboeck TS. Simulating cancer growth with multiscale agent-based modeling. Semin Cancer Biol. 2015;30:70–8.

    Article  PubMed  Google Scholar 

  21. Kourou K, Fotiadis DI. Computational modelling in cancer: methods and applications. Biomed Data Journal. 2015;1(1):15–25.

    Article  Google Scholar 

  22. Araujo RP, McElwain DL. A history of the study of solid tumour growth: the contribution of mathematical modelling. Bull Math Biol. 2004;66(5):1039–91.

    Article  CAS  PubMed  Google Scholar 

  23. Cobb K. Modeling cancer biology. Biomedical computation review. 2007.

    Google Scholar 

  24. Shirinifard A, Gens JS, Zaitlen BL, Poplawski NJ, Swat M, Glazier JA. 3D multi-cell simulation of tumor growth and angiogenesis. PloS one. 2009;4(10):e7190.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Vipin Narang SYW, Shiang Rong Leong, Jean-Pierre Abastado, Alexandre Gouaillard: Comparing mathematical models of cell adhesion in tumors. Defense Science Research Conference and Expo (DSR) 2011. p. 1–4.

  26. Qutub AA, Popel AS. Elongation, proliferation & migration differentiate endothelial cell phenotypes and determine capillary sprouting. BMC Syst Biol. 2009;3:13.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Mahoney AW, Podgorski GJ, Flann NS. A multi-objective optimization based-approach for discovering novel cancer therapies. IEEE/ACM transactions on computational biology and bioinformatics/IEEE, ACM. 2010.

    Google Scholar 

  28. Stamper IJ, Byrne HM, Owen MR, Maini PK. Modelling the role of angiogenesis and vasculogenesis in solid tumour growth. Bull Math Biol. 2007;69(8):2737–72.

    Article  CAS  PubMed  Google Scholar 

  29. Swat MH, Thomas GL, Belmonte JM, Shirinifard A, Hmeljak D, Glazier JA. Multi-scale modeling of tissues using CompuCell3D. Methods Cell Biol. 2012;110:325–66.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Mahoney AW, Smith BG, Flann NS, Podgorski GJ: Discovering novel cancer therapies: A computational modeling and search approach. In: Computational Intelligence in Bioinformatics and Computational Biology, 2008 CIBCB ‘08 IEEE Symposium on: 15–17 Sept. 2008 2008. p. 233–240.

  31. Anderson AR. A hybrid mathematical model of solid tumour invasion: the importance of cell adhesion. Math Med Biol. 2005;22(2):163–86.

    Article  PubMed  Google Scholar 

  32. Anderson AR, Chaplain MA. Continuous and discrete mathematical models of tumor-induced angiogenesis. Bull Math Biol. 1998;60(5):857–99.

    Article  CAS  PubMed  Google Scholar 

  33. Anderson AR, Quaranta V. Integrative mathematical oncology. Nat Rev Cancer. 2008;8(3):227–34.

    Article  CAS  PubMed  Google Scholar 

  34. Anderson AR, Rejniak KA, Gerlee P, Quaranta V. Microenvironment driven invasion: a multiscale multimodel investigation. J Math Biol. 2009;58(4–5):579–624.

    Article  PubMed  Google Scholar 

  35. Chaplain MA, McDougall SR, Anderson AR. Mathematical modeling of tumor-induced angiogenesis. Annu Rev Biomed Eng. 2006;8:233–57.

    Article  CAS  PubMed  Google Scholar 

  36. Gevertz JL, Gillies GT, Torquato S. Simulating tumor growth in confined heterogeneous environments. Phys Biol. 2008;5(3):36010.

    Article  Google Scholar 

  37. Quaranta V, Weaver AM, Cummings PT, Anderson AR. Mathematical modeling of cancer: the future of prognosis and treatment. Clin Chim Acta. 2005;357(2):173–9.

    Article  CAS  PubMed  Google Scholar 

  38. Hinow P, Gerlee P, McCawley LJ, Quaranta V, Ciobanu M, Wang S, Graham JM, Ayati BP, Claridge J, Swanson KR, et al. A spatial model of tumor-host interaction: application of chemotherapy. Math Biosci Eng. 2009;6(3):521–46.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Kansal AR, Torquato S, Harsh GI, Chiocca EA, Deisboeck TS. Simulated brain tumor growth dynamics using a three-dimensional cellular automaton. J Theor Biol. 2000;203(4):367–82.

    Article  CAS  PubMed  Google Scholar 

  40. Macklin P, McDougall S, Anderson AR, Chaplain MA, Cristini V, Lowengrub J. Multiscale modelling and nonlinear simulation of vascular tumour growth. J Math Biol. 2009;58(4–5):765–98.

    Article  PubMed  Google Scholar 

  41. McDougall SR, Anderson AR, Chaplain MA. Mathematical modelling of dynamic adaptive tumour-induced angiogenesis: clinical implications and therapeutic targeting strategies. J Theor Biol. 2006;241(3):564–89.

    Article  PubMed  Google Scholar 

  42. Athale C, Mansury Y, Deisboeck TS. Simulating the impact of a molecular ‘decision-process’ on cellular phenotype and multicellular patterns in brain tumors. J Theor Biol. 2005;233(4):469–81.

    Article  CAS  PubMed  Google Scholar 

  43. Wang Z, Zhang L, Sagotsky J, Deisboeck TS. Simulating non-small cell lung cancer with a multiscale agent-based model. Theor Biol Med Model. 2007;4:50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Kam N, Cohen IR, Harel D: The Immune System as a Reactive System: Modeling T Cell Activation with Statecharts. Proc Visual Languages and Formal Methods. IEEE Symp on Human-Centric Computing. 2001. p. 15–22.

  45. Efroni S, Harel D, Cohen IR. Toward rigorous comprehension of biological complexity: modeling, execution, and visualization of thymic T-cell maturation. Genome Res. 2003;13(11):2485–97.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Setty Y, Cohen IR, Dor Y, Harel D. Four-dimensional realistic modeling of pancreatic organogenesis. Proc Natl Acad Sci U S A. 2008;105(51):20374–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Swerdlin N, Cohen IR, Harel D. The Lymph Node B Cell Immune Response: Dynamic Analysis in-silico. In: Proceedings of the IEEE (special issue on Computational System Biology), 8, vol. 96. 2008. p. 1421–43.

    Google Scholar 

  48. Harel D. A grand challenge for computing: full reactive modeling of a multi-cellular animal. In: Bulletin of the EATCS, European Association for Theoretical Computer Science. 81st ed. 2003. p. 226–35.

    Google Scholar 

  49. Fisher J, Henzinger TA. Executable cell biology. Nat Biotechnol. 2007;25(11):1239–49.

    Article  CAS  PubMed  Google Scholar 

  50. Bloch N, Weiss G, Szekely S, Harel D. An interactive tool for animating biology, and its use in spatial and temporal modeling of a cancerous tumor and its microenvironment. PloS one. 2015;10(7):e0133484.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Harel D, Efroni S, Cohen IR. Reactive animation. Lec Notes Comp Sci. 2003;2852:136–53.

    Article  Google Scholar 

  52. Efroni S, Harel D, Cohen IR. Reactive animation: realistic modeling of complex dynamic systems. IEEE. 2005;38(1):38–47.

    Google Scholar 

  53. Appel B, Givan LA, Eisen JS. Delta-notch signaling and lateral inhibition in zebrafish spinal cord development. BMC Dev Biol. 2001;1:13.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Harvey NL. To sprout or “notch” to sprout? Blood. 2011;118(4):836–7.

    Article  CAS  Google Scholar 

  55. Elenbaas B, Weinberg RA. Heterotypic signaling between epithelial tumor cells and fibroblasts in carcinoma formation. Exp Cell Res. 2001;264(1):169–84.

    Article  CAS  PubMed  Google Scholar 

  56. Giaccia AJ, Schipani E. Role of carcinoma-associated fibroblasts and hypoxia in tumor progression. Curr Top Microbiol Immunol. 2010;345:31–45.

    CAS  PubMed  Google Scholar 

  57. Beacham DA, Cukierman E. Stromagenesis: the changing face of fibroblastic microenvironments during tumor progression. Semin Cancer Biol. 2005;15(5):329–41.

    Article  PubMed  Google Scholar 

  58. Bhowmick NA, Neilson EG, Moses HL. Stromal fibroblasts in cancer initiation and progression. Nature. 2004;432(7015):332–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Chary SR, Jain RK. Direct measurement of interstitial convection and diffusion of albumin in normal and neoplastic tissues by fluorescence photobleaching. Proc Natl Acad Sci U S A. 1989;86(14):5385–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Chen RR, Silva EA, Yuen WW, Mooney DJ. Spatio-temporal VEGF and PDGF delivery patterns blood vessel formation and maturation. Pharm Res. 2007;24(2):258–64.

    Article  PubMed  Google Scholar 

  61. Ausprunk DH, Folkman J. Migration and proliferation of endothelial cells in preformed and newly formed blood vessels during tumor angiogenesis. Microvasc Res. 1977;14(1):53–65.

    Article  CAS  PubMed  Google Scholar 

  62. Gibby K, You WK, Kadoya K, Helgadottir H, Young LJ, Ellies LG, Chang Y, Cardiff RD, Stallcup WB. Early vascular deficits are correlated with delayed mammary tumorigenesis in the MMTV-PyMT transgenic mouse following genetic ablation of the NG2 proteoglycan. Breast Cancer Res. 2012;14(2):R67.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Laufer J, Johnson P, Zhang E, Treeby B, Cox B, Pedley B, Beard P. In vivo preclinical photoacoustic imaging of tumor vasculature development and therapy. J Biomed Opt. 2012;17(5):056016.

    Article  PubMed  Google Scholar 

  64. Freyer JP, Sutherland RM. Regulation of growth saturation and development of necrosis in EMT6/Ro multicellular spheroids by the glucose and oxygen supply. Cancer Res. 1986;46(7):3504–12.

    CAS  PubMed  Google Scholar 

  65. Folkman J, Hochberg M. Self-regulation of growth in three dimensions. J Exp Med. 1973;138(4):745–53.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  66. Bauer AL, Jackson TL, Jiang Y. A cell-based model exhibiting branching and anastomosis during tumor-induced angiogenesis. Biophys J. 2007;92(9):3105–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  67. Byrne HM, Chaplain MAJ. Explicit solutions of a simplified model of capillary sprout growth during tumor angiogenesis. Appl Math Lett. 1996;9(1):69–74.

    Article  Google Scholar 

  68. Muthukkaruppan VR, Kubai L, Auerbach R. Tumor-induced neovascularization in the mouse eye. J Natl Cancer Inst. 1982;69(3):699–708.

    CAS  PubMed  Google Scholar 

  69. Helm CL, Fleury ME, Zisch AH, Boschetti F, Swartz MA. Synergy between interstitial flow and VEGF directs capillary morphogenesis in vitro through a gradient amplification mechanism. Proc Natl Acad Sci U S A. 2005;102(44):15779–84.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  70. Tomayko MM, Reynolds CP. Determination of subcutaneous tumor size in athymic (nude) mice. Cancer Chemother Pharmacol. 1989;24(3):148–54.

    Article  CAS  PubMed  Google Scholar 

  71. Laird AK. Dynamics of tumor growth. Br J Cancer. 1964;13:490–502.

    Article  CAS  PubMed  Google Scholar 

  72. Zhang HT, Craft P, Scott PA, Ziche M, Weich HA, Harris AL, Bicknell R. Enhancement of tumor growth and vascular density by transfection of vascular endothelial cell growth factor into MCF-7 human breast carcinoma cells. J Natl Cancer Inst. 1995;87(3):213–9.

    Article  CAS  PubMed  Google Scholar 

  73. Kanwar JR, Berg RW, Yang Y, Kanwar RK, Ching LM, Sun X, Krissansen GW. Requirements for ICAM-1 immunogene therapy of lymphoma. Cancer Gene Ther. 2003;10(6):468–76.

    Article  CAS  PubMed  Google Scholar 

  74. Avni R, Cohen B, Neeman M. Hypoxic stress and cancer: imaging the axis of evil in tumor metastasis. NMR Biomed. 2011;24(6):569–81.

    CAS  PubMed  PubMed Central  Google Scholar 

  75. Fukumura D, Jain RK. Imaging angiogenesis and the microenvironment. APMIS. 2008;116(7–8):695–715.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  76. Ribatti D, Nico B, Crivellato E, Roccaro AM, Vacca A. The history of the angiogenic switch concept. Leukemia. 2007;21(1):44–52.

    Article  CAS  PubMed  Google Scholar 

  77. Folkman U, Kalluri R. Beginning of angiogenesis research. In: Kufe D, Pollock R, Weichselbaum R, editors. Holland-frei cancer medicine. 6th ed. 2003.

    Google Scholar 

  78. Bullitt E, Ewend MG, Aylward S, Lin W, Gerig G, Joshi S, Jung I, Muller K, Smith JK. Abnormal vessel tortuosity as a marker of treatment response of malignant gliomas: preliminary report. Technol Cancer Res Treat. 2004;3(6):577–84.

    Article  PubMed  PubMed Central  Google Scholar 

  79. Chen Y, Cairns R, Papandreou I, Koong A, Denko NC. Oxygen consumption can regulate the growth of tumors, a new perspective on the Warburg effect. PLoS One. 2009;4(9):e7033.

    Article  PubMed  PubMed Central  Google Scholar 

  80. Warburg O. On the origin of cancer cells. Science. 1956;123(3191):309–14.

    Article  CAS  PubMed  Google Scholar 

  81. Giaccia AJ. Hypoxic stress proteins: survival of the fittest. Semin Radiat Oncol. 1996;6(1):46–58.

    Article  PubMed  Google Scholar 

  82. Kitano H. Cancer as a robust system: implications for anticancer therapy. Nat Rev Cancer. 2004;4(3):227–35.

    Article  CAS  PubMed  Google Scholar 

  83. Brown PO, Palmer C. The preclinical natural history of serous ovarian cancer: defining the target for early detection. PLoS Med. 2009;6(7):e1000114.

    Article  PubMed  PubMed Central  Google Scholar 

  84. Dome B, Paku S, Somlai B, Timar J. Vascularization of cutaneous melanoma involves vessel co-option and has clinical significance. J Pathol. 2002;197(3):355–62.

    Article  PubMed  Google Scholar 

  85. Ribatti D. The history of angiogenesis inhibitors. Leukemia. 2007;21(8):1606–9.

    Article  CAS  PubMed  Google Scholar 

  86. Kerbel RS. Tumor angiogenesis: past, present and the near future. Carcinogenesis. 2000;21(3):505–15.

    Article  CAS  PubMed  Google Scholar 

  87. Hanrahan V, Currie MJ, Gunningham SP, Morrin HR, Scott PA, Robinson BA, Fox SB. The angiogenic switch for vascular endothelial growth factor (VEGF)-A, VEGF-B, VEGF-C, and VEGF-D in the adenoma-carcinoma sequence during colorectal cancer progression. J Pathol. 2003;200(2):183–94.

    Article  CAS  PubMed  Google Scholar 

  88. Ribatti D, Vacca A, Dammacco F. The role of the vascular phase in solid tumor growth: a historical review. Neoplasia. 1999;1(4):293–302.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  89. Folkman J. What is the evidence that tumors are angiogenesis dependent? J Natl Cancer Inst. 1990;82(1):4–6.

    Article  CAS  PubMed  Google Scholar 

  90. Koukourakis GV, Sotiropoulou-Lontou A. Targeted therapy with bevacizumab (Avastin) for metastatic colorectal cancer. Clinical Translational Oncol. 2011;13(10):710–4.

    Article  CAS  Google Scholar 

  91. Fleshker S, Preise D, Kalchenko V, Scherz A, Salomon Y. Prompt assessment of WST11-VTP outcome using luciferase transfected tumors enables second treatment and increase in overall therapeutic rate. Photochem Photobiol. 2008;84(5):1231–7.

    Article  CAS  PubMed  Google Scholar 

  92. Koudinova NV, Pinthus JH, Brandis A, Brenner O, Bendel P, Ramon J, Eshhar Z, Scherz A, Salomon Y. Photodynamic therapy with Pd-bacteriopheophorbide (TOOKAD): successful in vivo treatment of human prostatic small cell carcinoma xenografts. Int J Cancer. 2003;104(6):782–9.

    Article  CAS  PubMed  Google Scholar 

  93. Zilberstein J, Schreiber S, Bloemers MC, Bendel P, Neeman M, Schechtman E, Kohen F, Scherz A, Salomon Y. Antivascular treatment of solid melanoma tumors with bacteriochlorophyll-serine-based photodynamic therapy. Photochem Photobiol. 2001;73(3):257–66.

    Article  CAS  PubMed  Google Scholar 

  94. Azzouzi AR, Barret E, Moore CM, Villers A, Allen C, Scherz A, Muir G, de Wildt M, Barber NJ, Lebdai S, et al. TOOKAD ((R)) Soluble vascular-targeted photodynamic (VTP) therapy: determination of optimal treatment conditions and assessment of effects in patients with localised prostate cancer. BJU Int. 2013;112(6):766–74.

    Article  CAS  PubMed  Google Scholar 

  95. Ben-Nun Y, Merquiol E, Brandis A, Turk B, Scherz A, Blum G. Photodynamic quenched cathepsin activity based probes for cancer detection and macrophage targeted therapy. Theranostics. 2015;5(8):847–62.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  96. Harel D. Statecharts: a visual formalism for complex systems. Sci Comput Programming. 1987;8:231–74.

    Article  Google Scholar 

  97. Harel D, Gery E: Executable object modeling with Statecharts. Computer, IEEE computer society press 1997;30(7):31-42.

  98. Harel D, Pnueli A: On the development of reactive systems. In: Logics and models of concurrent systems. Edited by Krzysztof RA: Springer-Verlag New York, Inc.; 1985: 477-498.

  99. Fisher J, Piterman N, Hubbard EJ, Stern MJ, Harel D. Computational insights into Caenorhabditis elegans vulval development. Proc Natl Acad Sci U S A. 2005;102(6):1951–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references


We are deeply grateful to Irun Cohen for many helpful suggestions throughout the research period, and also for originally sowing in our minds the seed of the idea that cancer could perhaps be modeled as analogous to organogenesis. We also thank Raviv Kribus for his help in creating automatic analysis for simulation results, and Guy Aidelberg for useful information on Fibroblasts that was used in the model.


This research was supported in part by an Advanced Research Grant from the European Union’s ERC Programme (FP7/2007-2013), by the Israel Science Foundation, and by the Israeli Center of Research Excellence in Algorithms (I-Core ALGO).

Availability of data and materials

The data sets supporting the results of this article are included within the article and its additional file.

Authors’ contributions

DH and NB conceived the research project. NB developed the model, analyzed the data, and wrote the manuscript. DH supervised the entire research and manuscript writing. All authors read and approved the final manuscript.

Authors’ information

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Not applicable.


Not applicable.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Naamah Bloch.

Additional files

Additional file 1.

Further details of the model. Details regarding the specification of the model, together with a list of the parameters used in the model and their default values, and the range of values for VEGF and oxygen parameters that ensure tumor recovery. (PDF 244 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Bloch, N., Harel, D. The tumor as an organ: comprehensive spatial and temporal modeling of the tumor and its microenvironment. BMC Bioinformatics 17, 317 (2016).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: