No Arabic abstract
The Cellular Potts Model (CPM) is a lattice based modeling technique which is widely used for simulating cellular patterns such as foams or biological tissues. Despite its realism and generality, the standard Monte Carlo algorithm used in the scientific literature to evolve this model preserves connectivity of cells on a limited range of simulation temperature only. We present a new algorithm in which cell fragmentation is forbidden for all simulation temperatures. This allows to significantly enhance realism of the simulated patterns. It also increases the computational efficiency compared with the standard CPM algorithm even at same simulation temperature, thanks to the time spared in not doing unrealistic moves. Moreover, our algorithm restores the detailed balance equation, ensuring that the long-term stage is independent of the chosen acceptance rate and chosen path in the temperature space.
Soft cellular systems, such as foams or biological tissues, exhibit highly complex rheological properties, even in the quasistatic regime, that numerical modeling can help to apprehend. We present a numerical implementation of quasistatic strain within the widely used cellular Potts model. The accuracy of the method is tested by simulating the quasistatic strain 2D dry foams, both ordered and disordered. The implementation of quasistatic strain in CPM allows the investigation of sophisticated interplays between stress-strain relationship and structural changes that take place in cellular systems.
In the absence of environmental cues, a migrating cell performs an isotropic random motion. Recently, the breaking of this isotropy has been observed when cells move in the presence of asymmetric adhesive patterns. However, up to now the mechanisms at work to direct cell migration in such environments remain unknown. Here, we show that a non-adhesive surface with asymmetric micro-geometry consisting of dense arrays of tilted micro-pillars can direct cell motion. Our analysis reveals that most features of cell trajectories, including the bias, can be reproduced by a simple model of active Brownian particle in a ratchet potential, which we suggest originates from a generic elastic interaction of the cell body with the environment. The observed guiding effect, independent of adhesion, is therefore robust and could be used to direct cell migration both in vitro and in vivo.
Heterogeneity is a hallmark of all cancers. Tumor heterogeneity is found at different levels -- interpatient, intrapatient, and intratumor heterogeneity. All of them pose challenges for clinical treatments. The latter two scenarios can also increase the risk of developing drug resistance. Although the existence of tumor heterogeneity has been known for two centuries, a clear understanding of its origin is still elusive, especially at the level of intratumor heterogeneity (ITH). The coexistence of different subpopulations within a single tumor has been shown to play crucial roles during all stages of carcinogenesis. Here, using concepts from evolutionary game theory and public goods game, often invoked in the context of the tragedy of commons, we explore how the interactions among subclone populations influence the establishment of ITH. By using an evolutionary model, which unifies several experimental results in distinct cancer types, we develop quantitative theoretical models for explaining data from {it in vitro} experiments involving pancreatic cancer as well as {it vivo} data in glioblastoma multiforme. Such physical and mathematical models complement experimental studies, and could optimistically also provide new ideas for the design of efficacious therapies for cancer patients.
We introduce and analyze several aspects of a new model for cell differentiation. It assumes that differentiation of progenitor cells is a continuous process. From the mathematical point of view, it is based on partial differential equations of transport type. Specifically, it consists of a structured population equation with a nonlinear feedback loop. This models the signaling process due to cytokines, which regulate the differentiation and proliferation process. We compare the continuous model to its discrete counterpart, a multi-compartmental model of a discrete collection of cell subpopulations recently proposed by Marciniak-Czochra et al. in 2009 to investigate the dynamics of the hematopoietic system. We obtain uniform bounds for the solutions, characterize steady state solutions, and analyze their linearized stability. We show how persistence or extinction might occur according to values of parameters that characterize the stem cells self-renewal. We also perform numerical simulations and discuss the qualitative behavior of the continuous model vis a vis the discrete one.
We study a mathematical model describing the dynamics of a pluripotent stem cell population involved in the blood production process in the bone marrow. This model is a differential equation with a time delay. The delay describes the cell cycle duration and is uniformly distributed on an interval. We obtain stability conditions independent of the delay. We also show that the distributed delay can destabilize the entire system. In particularly, it is shown that Hopf bifurcations can occur.