A Boolean dynamical system integrating the main signaling pathways involved in cancer is constructed based on the currently known protein-protein interaction network. This system exhibits stationary protein activation patterns – attractors – dependent on the cell's microenvironment. These dynamical attractors were determined through simulations and their stabilities against mutations were tested. In a higher hierarchical level, it was possible to group the network attractors into distinct cell phenotypes and determine driver mutations that promote phenotypic transitions. We find that driver nodes are not necessarily central in the network topology, but at least they are direct regulators of central components towards which converge or through which crosstalk distinct cancer signaling pathways. The predicted drivers are in agreement with those pointed out by diverse census of cancer genes recently performed for several human cancers. Furthermore, our results demonstrate that cell phenotypes can evolve towards full malignancy through distinct sequences of accumulated mutations. In particular, the network model supports routes of carcinogenesis known for some tumor types. Finally, the Boolean network model is employed to evaluate the outcome of molecularly targeted cancer therapies. The major find is that monotherapies were additive in their effects and that the association of targeted drugs is necessary for cancer eradication.

Recently, we have proposed a nutrient-limited model for the avascular growth of tumors including cell proliferation, motility, and death [S. C. Ferreira, Jr., M. L. Martins, and M. J. Vilela, Phys. Rev. E 65, 021907 (2002)], which qualitatively reproduces commonly observed morphologies for carcinomas in situ. In the present work, we analyze the effects of distinct chemotherapeutic strategies on the patterns, scaling, and growth laws obtained for the nutrient-limited model. Two kinds of chemotherapeutic strategies were considered, namely, those that kill cancer cells and those that block cell mitosis but allow the cell to survive for some time. Depending on the chemotherapeutic schedule used, the tumors are completely eliminated, reach a stationary size, or grow following power laws. The model suggests that the scaling properties of the tumors are not affected by the mild cytotoxic treatments, although a reduction in growth rates and an increase in invasiveness are observed. For the strategies based on antimitotic drugs, a morphological transition in which compact tumors become more fractal under aggressive treatments was seen.

We review some computer algorithms for the simulation of off-lattice clusters grown from a seed, with emphasis on the diffusion-limited aggregation, ballistic aggregation and Eden models. Only those methods which can be immediately extended to distinct off-lattice aggregation processes are discussed. The computer efficiencies of the distinct algorithms are compared.

In multicellular organisms, cell motility is central in all morphogenetic processes, tissue maintenance, wound healing and immune surveillance. Hence, the control of cell motion is a major demand in the creation of artificial tissues and organs. Here, cell migration assays on plastic 2D surfaces involving normal (MDCK) and tumoral (B16F10) epithelial cell lines were performed varying the initial density of plated cells. Through time-lapse microscopy quantities such as speed distributions, velocity autocorrelations and spatial correlations, as well as the scaling of mean-squared displacements were determined. We find that these cells exhibit anomalous diffusion with q-Weibull speed distributions that evolves non-monotonically to a Maxwellian distribution as the initial density of plated cells increases. Although short-ranged spatial velocity correlations mark the formation of small cell clusters, the emergence of collective motion was not observed. Finally, simulational results from a correlated random walk and the Vicsek model of collective dynamics evidence that fluctuations in cell velocity orientations are sufficient to produce q-Weibull speed distributions seen in our migration assays.

Microorganisms and plants very commonly release toxic secondary chemical compounds (allelochemicals) that inhibit or kill sensitive strains or individuals from their own or other species. In this work we study a model that describes two species interacting through allelopathic suppression and competing for resources. Employing linear stability analysis, the conditions for coexistence or extinction of species in spatially homogeneous systems were determined. We found that the borders between the regimes of bistability, coexistence, and the extinction of the weaker by the stronger competitor, are altered by allelopathic interactions. In addition, traveling wave solutions for one species invasion were obtained considering the spatially explicit nature of the model. Our findings indicate that the minimum speed of the invasion wavefronts depends primarily on the competition coefficients and the parameters characterizing the species' functional responses to their allelochemicals. As a general rule, the species provided with the most effective chemical weapons dominates the population dynamics. Finally, we found a tristability at the coexistence region due to the combination of allelopathy and patchy population distributions in space. So, our model provides a distinct mechanism, independent of social behaviors, that produces such unexpected tristability impossible in classical competition models involving one-to-one individual interactions.