  Understanding tumor invasion and metastasis is of crucial importance for both
fundamental cancer research and clinical practice. In vitro experiments have
established that the invasive growth of malignant tumors is characterized by
the dendritic invasive branches composed of chains of tumor cells emanating
from the primary tumor mass. The preponderance of previous tumor simulations
focused on non-invasive (or proliferative) growth. The formation of the
invasive cell chains and their interactions with the primary tumor mass and
host microenvironment are not well understood. Here, we present a novel
cellular automaton (CA) model that enables one to efficiently simulate invasive
tumor growth in a heterogeneous host microenvironment. By taking into account a
variety of microscopic-scale tumor-host interactions, including the short-range
mechanical interactions between tumor cells and tumor stroma, degradation of
extracellular matrix by the invasive cells and oxygen/nutrient gradient driven
cell motions, our CA model predicts a rich spectrum of growth dynamics and
emergent behaviors of invasive tumors. Besides robustly reproducing the salient
features of dendritic invasive growth, such as least resistance and intrabranch
homotype attraction, we also predict nontrivial coupling of the growth dynamics
of the primary tumor mass and the invasive cells. In addition, we show that the
properties of the host microenvironment can significantly affect tumor
morphology and growth dynamics, emphasizing the importance of understanding the
tumor-host interaction. The capability of our CA model suggests that
well-developed in silico tools could eventually be utilized in clinical
situations to predict neoplastic progression and propose individualized optimal
treatment strategies.
