An immersed boundary framework for modelling the growth of individual cells: an application to the early tumour development

Mechanical deformation and death of circulating tumor cells in the bloodstream

The circulation of tumor cells through the bloodstream is a significant step in tumor metastasis. To better understand the metastatic process, circulating tumor cell (CTC) survival in the circulation must be explored. While immune interactions with CTCs in recent decades have been examined, research has yet to sufficiently explain some CTC behaviors in blood flow. Studies related to CTC mechanical responses in the bloodstream have recently been conducted to further study conditions under which CTCs might die. While experimental methods can assess the mechanical properties and death of CTCs, increasingly sophisticated computational models are being built to simulate the blood flow and CTC mechanical deformation under fluid shear stresses (FSS) in the bloodstream.Several factors contribute to the mechanical deformation and death of CTCs as they circulate. While FSS can damage CTC structure, diverse interactions between CTCs and blood components may either promote or hinder the next metastatic step-extravasation at a remote site. Overall understanding of how these factors influence the deformation and death of CTCs could serve as a basis for future experiments and simulations, enabling researchers to predict CTC death more accurately. Ultimately, these efforts can lead to improved metastasis-specific therapeutics and diagnostics specific in the future.

Selected aspects of avascular tumor growth reproduced by a hybrid model of cell dynamics and chemical kinetics

We here propose a hybrid computational framework to reproduce and analyze aspects of the avascular progression of a generic solid tumor. Our method first employs an individual-based approach to represent the population of tumor cells, which are distinguished in viable and necrotic agents. The active part of the disease is in turn differentiated according to a set of metabolic states. We then describe the spatio-temporal evolution of the concentration of oxygen and of tumor-secreted proteolytic enzymes using partial differential equations (PDEs). A differential equation finally governs the local degradation of the extracellular matrix (ECM) by the malignant mass. Numerical realizations of the model are run to reproduce tumor growth and invasion in a number scenarios that differ for cell properties (adhesiveness, duplication potential, proteolytic activity) and/or environmental conditions (level of tissue oxygenation and matrix density pattern). In particular, our simulations suggest that tumor aggressiveness, in terms of invasive depth and extension of necrotic tissue, can be reduced by (i) stable cell-cell contact interactions, (ii) poor tendency of malignant agents to chemotactically move upon oxygen gradients, and (iii) presence of an overdense matrix, if coupled by a disrupted proteolytic activity of the disease.

The Effect of Substrate Stiffness on Elastic Force Transmission in the Epithelial Monolayers over Short Timescales

Purpose

The importance of mechanical forces and microenvironment in guiding cellular behavior has been widely accepted. Together with the extracellular matrix (ECM), epithelial cells form a highly connected mechanical system subjected to various mechanical cues from their environment, such as ECM stiffness, and tensile and compressive forces. ECM stiffness has been linked to many pathologies, including tumor formation. However, our understanding of the effect of ECM stiffness and its heterogeneities on rapid force transduction in multicellular systems has not been fully addressed.

Methods

We used experimental and computational methods. Epithelial cells were cultured on elastic hydrogels with fluorescent nanoparticles. Single cells were moved by a micromanipulator, and epithelium and substrate deformation were recorded. We developed a computational model to replicate our experiments and quantify the force distribution in the epithelium. Our model further enabled simulations with local stiffness gradients.

Results

We found that substrate stiffness affects the force transduction and the cellular deformation following an external force. Also, our results indicate that the heterogeneities, e.g., gradients, in the stiffness can substantially influence the strain redistribution in the cell monolayers. Furthermore, we found that the cells' apico-basal elasticity provides a level of mechanical isolation between the apical cell-cell junctions and the basal focal adhesions.

Conclusions

Our simulation results show that increased ECM stiffness, e.g., due to a tumor, can mechanically isolate cells and modulate rapid mechanical signaling between cells over distances. Furthermore, the developed model has the potential to facilitate future studies on the interactions between epithelial monolayers and elastic substrates.

Supplementary Information

The online version of this article (10.1007/s12195-023-00772-0) contains supplementary material, which is available to authorized users.

Free and Interfacial Boundaries in Individual-Based Models of Multicellular Biological systems

Coordination of cell behaviour is key to a myriad of biological processes including tissue morphogenesis, wound healing, and tumour growth. As such, individual-based computational models, which explicitly describe inter-cellular interactions, are commonly used to model collective cell dynamics. However, when using individual-based models, it is unclear how descriptions of cell boundaries affect overall population dynamics. In order to investigate this we define three cell boundary descriptions of varying complexities for each of three widely used off-lattice individual-based models: overlapping spheres, Voronoi tessellation, and vertex models. We apply our models to multiple biological scenarios to investigate how cell boundary description can influence tissue-scale behaviour. We find that the Voronoi tessellation model is most sensitive to changes in the cell boundary description with basic models being inappropriate in many cases. The timescale of tissue evolution when using an overlapping spheres model is coupled to the boundary description. The vertex model is demonstrated to be the most stable to changes in boundary description, though still exhibits timescale sensitivity. When using individual-based computational models one should carefully consider how cell boundaries are defined. To inform future work, we provide an exploration of common individual-based models and cell boundary descriptions in frequently studied biological scenarios and discuss their benefits and disadvantages.

Computational approaches for simulating luminogenesis

Lumens, liquid-filled cavities surrounded by polarized tissue cells, are elementary units involved in the morphogenesis of organs. Theoretical modeling and computations, which can integrate various factors involved in biophysics of morphogenesis of cell assembly and lumens, may play significant roles to elucidate the mechanisms in formation of such complex tissue with lumens. However, up to present, it has not been documented well what computational approaches or frameworks can be applied for this purpose and how we can choose the appropriate approach for each problem. In this review, we report some typical lumen morphologies and basic mechanisms for the development of lumens, focusing on three keywords - mechanics, hydraulics and geometry - while outlining pros and cons of the current main computational strategies. We also describe brief guidance of readouts, i.e., what we should measure in experiments to make the comparison with the model's assumptions and predictions.

The Role of in silico Research in Developing Nanoparticle-Based Therapeutics

Nanoparticles (NPs) hold great potential as therapeutics, particularly in the realm of drug delivery. They are effective at functional cargo delivery and offer a great degree of amenability that can be used to offset toxic side effects or to target drugs to specific regions in the body. However, there are many challenges associated with the development of NP-based drug formulations that hamper their successful clinical translation. Arguably, the most significant barrier in the way of efficacious NP-based drug delivery systems is the tedious and time-consuming nature of NP formulation—a process that needs to account for downstream effects, such as the onset of potential toxicity or immunogenicity, in vivo biodistribution and overall pharmacokinetic profiles, all while maintaining desirable therapeutic outcomes. Computational and AI-based approaches have shown promise in alleviating some of these restrictions. Via predictive modeling and deep learning, in silico approaches have shown the ability to accurately model NP-membrane interactions and cellular uptake based on minimal data, such as the physicochemical characteristics of a given NP. More importantly, machine learning allows computational models to predict how specific changes could be made to the physicochemical characteristics of a NP to improve functional aspects, such as drug retention or endocytosis. On a larger scale, they are also able to predict the in vivo pharmacokinetics of NP-encapsulated drugs, predicting aspects such as circulatory half-life, toxicity, and biodistribution. However, the convergence of nanomedicine and computational approaches is still in its infancy and limited in its applicability. The interactions between NPs, the encapsulated drug and the body form an intricate network of interactions that cannot be modeled with absolute certainty. Despite this, rapid advancements in the area promise to deliver increasingly powerful tools capable of accelerating the development of advanced nanoscale therapeutics. Here, we describe computational approaches that have been utilized in the field of nanomedicine, focusing on approaches for NP design and engineering.

Quantitative modeling identifies critical cell mechanics driving bile duct lumen formation

Biliary ducts collect bile from liver lobules, the smallest functional and anatomical units of liver, and carry it to the gallbladder. Disruptions in this process caused by defective embryonic development, or through ductal reaction in liver disease have a major impact on life quality and survival of patients. A deep understanding of the processes underlying bile duct lumen formation is crucial to identify intervention points to avoid or treat the appearance of defective bile ducts. Several hypotheses have been proposed to characterize the biophysical mechanisms driving initial bile duct lumen formation during embryogenesis. Here, guided by the quantification of morphological features and expression of genes in bile ducts from embryonic mouse liver, we sharpened these hypotheses and collected data to develop a high resolution individual cell-based computational model that enables to test alternative hypotheses in silico. This model permits realistic simulations of tissue and cell mechanics at sub-cellular scale. Our simulations suggest that successful bile duct lumen formation requires a simultaneous contribution of directed cell division of cholangiocytes, local osmotic effects generated by salt excretion in the lumen, and temporally-controlled differentiation of hepatoblasts to cholangiocytes, with apical constriction of cholangiocytes only moderately affecting luminal size.

The initial step in bile duct development is the formation of a biliary lumen, a process which involves several cellular mechanisms, such as cell division and polarization, and secretion of fluid. However, how these mechanisms are orchestrated in time and space is difficult to understand. Here, we built a computational model of biliary lumen formation which represents every cell and its function in detail. With the model we can simulate the effect of biophysical aspects that affect duct formation. We have tested the individual and combined effects of directed cell division, apical constriction, and osmotic effects on lumen expansion by varying the parameters that control their relative strength. Our simulations suggest that successful bile duct lumen formation requires the simultaneous contribution of directed cell division of cholangiocytes, local osmotic effects generated by salt excretion in the lumen, and temporally-controlled differentiation of hepatoblasts to cholangiocytes, with apical constriction of cholangiocytes only moderately affecting luminal size.

Optimal regulation of tumour-associated neutrophils in cancer progression

In a tumour microenvironment, tumour-associated neutrophils could display two opposing differential phenotypes: anti-tumour (N1) and pro-tumour (N2) effector cells. Converting N2 to N1 neutrophils provides innovative therapies for cancer treatment. In this study, a mathematical model for N1-N2 dynamics describing the cancer survival and immune inhibition in response to TGF-β and IFN-β is considered. The effects of exogenous intervention of TGF-β inhibitor and IFN-β are examined in order to enhance N1 recruitment to combat tumour progression. Our approach employs optimal control theory to determine drug infusion protocols that could minimize tumour volume with least administration cost possible. Four optimal control scenarios corresponding to different therapeutic strategies are explored, namely, TGF-β inhibitor control only, IFN-β control only, concomitant TGF-β inhibitor and IFN-β controls, and alternating TGF-β inhibitor and IFN-β controls. For each scheme, different initial conditions are varied to depict different pathophysiological condition of a cancer patient, leading to adaptive treatment schedule. TGF-β inhibitor and IFN-β drug dosages, total drug amount, infusion times and relative cost of drug administrations are obtained under various circumstances. The control strategies achieved could guide in designing individualized therapeutic protocols.

CBMOS: a GPU-enabled Python framework for the numerical study of center-based models

BACKGROUND

Cell-based models are becoming increasingly popular for applications in developmental biology. However, the impact of numerical choices on the accuracy and efficiency of the simulation of these models is rarely meticulously tested. Without concrete studies to differentiate between solid model conclusions and numerical artifacts, modelers are at risk of being misled by their experiments' results. Most cell-based modeling frameworks offer a feature-rich environment, providing a wide range of biological components, but are less suitable for numerical studies. There is thus a need for software specifically targeted at this use case.

RESULTS

We present CBMOS, a Python framework for the simulation of the center-based or cell-centered model. Contrary to other implementations, CBMOS' focus is on facilitating numerical study of center-based models by providing access to multiple ordinary differential equation solvers and force functions through a flexible, user-friendly interface and by enabling rapid testing through graphics processing unit (GPU) acceleration. We show-case its potential by illustrating two common workflows: (1) comparison of the numerical properties of two solvers within a Jupyter notebook and (2) measuring average wall times of both solvers on a high performance computing cluster. More specifically, we confirm that although for moderate accuracy levels the backward Euler method allows for larger time step sizes than the commonly used forward Euler method, its additional computational cost due to being an implicit method prohibits its use for practical test cases.

CONCLUSIONS

CBMOS is a flexible, easy-to-use Python implementation of the center-based model, exposing both basic model assumptions and numerical components to the user. It is available on GitHub and PyPI under an MIT license. CBMOS allows for fast prototyping on a central processing unit for small systems through the use of NumPy. Using CuPy on a GPU, cell populations of up to 10,000 cells can be simulated within a few seconds. As such, it will substantially lower the time investment for any modeler to check the crucial assumption that model conclusions are independent of numerical issues.

A computational model of glioma reveals opposing, stiffness-sensitive effects of leaky vasculature and tumor growth on tissue mechanical stress and porosity

A biphasic computational model of a growing, vascularized glioma within brain tissue was developed to account for unique features of gliomas, including soft surrounding brain tissue, their low stiffness relative to brain tissue, and a lack of draining lymphatics. This model is the first to couple nonlinear tissue deformation with porosity and tissue hydraulic conductivity to study the mechanical interaction of leaky vasculature and solid growth in an embedded glioma. The present model showed that leaky vasculature and elevated interstitial fluid pressure produce tensile stress within the tumor in opposition to the compressive stress produced by tumor growth. This tensile effect was more pronounced in softer tissue and resulted in a compressive stress concentration at the tumor rim that increased when tumor was softer than host. Aside from generating solid stress, fluid pressure-driven tissue deformation decreased the effective stiffness of the tumor while growth increased it, potentially leading to elevated stiffness in the tumor rim. A novel prediction of reduced porosity at the tumor rim was corroborated by direct comparison with estimates from our in vivo imaging studies. Antiangiogenic and radiation therapy were simulated by varying vascular leakiness and tissue hydraulic conductivity. These led to greater solid compression and interstitial pressure in the tumor, respectively, the former of which may promote tumor infiltration of the host. Our findings suggest that vascular leakiness has an important influence on in vivo solid stress, stiffness, and porosity fields in gliomas given their unique mechanical microenvironment.

A data-driven approach to modeling cancer cell mechanics during microcirculatory transport

In order to understand the effect of cellular level features on the transport of circulating cancer cells in the microcirculation, there has been an increasing reliance on high-resolution in silico models. Accurate simulation of cancer cells flowing with blood cells requires resolving cellular-scale interactions in 3D, which is a significant computational undertaking warranting a cancer cell model that is both computationally efficient yet sufficiently complex to capture relevant behavior. Given that the characteristics of metastatic spread are known to depend on cancer type, it is crucial to account for mechanistic behavior representative of a specific cancer's cells. To address this gap, in the present work we develop and validate a means by which an efficient and popular membrane model-based approach can be used to simulate deformable cancer cells and reproduce experimental data from specific cell lines. Here, cells are modeled using the immersed boundary method (IBM) within a lattice Boltzmann method (LBM) fluid solver, and the finite element method (FEM) is used to model cell membrane resistance to deformation. Through detailed comparisons with experiments, we (i) validate this model to represent cancer cells undergoing large deformation, (ii) outline a systematic approach to parameterize different cell lines to optimally fit experimental data over a range of deformations, and (iii) provide new insight into nucleated vs. non-nucleated cell models and their ability to match experiments. While many works have used the membrane-model based method employed here to model generic cancer cells, no quantitative comparisons with experiments exist in the literature for specific cell lines undergoing large deformation. Here, we describe a phenomenological, data-driven approach that can not only yield good agreement for large deformations, but explicitly detail how it can be used to represent different cancer cell lines. This model is readily incorporated into cell-resolved hemodynamic transport simulations, and thus offers significant potential to complement experiments towards providing new insights into various aspects of cancer progression.

Biophysical interactions between components of the tumor microenvironment promote metastasis

During metastasis, tumor cells need to adapt to their dynamic microenvironment and modify their mechanical properties in response to both chemical and mechanical stimulation. Physical interactions occur between cancer cells and the surrounding matrix including cell movements and cell shape alterations through the process of mechanotransduction. The latter describes the translation of external mechanical cues into intracellular biochemical signaling. Reorganization of both the cytoskeleton and the extracellular matrix (ECM) plays a critical role in these spreading steps. Migrating tumor cells show increased motility in order to cross the tumor microenvironment, migrate through ECM and reach the bloodstream to the metastatic site. There are specific factors affecting these processes, as well as the survival of circulating tumor cells (CTC) in the blood flow until they finally invade the secondary tissue to form metastasis. This review aims to study the mechanisms of metastasis from a biomechanical perspective and investigate cell migration, with a focus on the alterations in the cytoskeleton through this journey and the effect of biologic fluids on metastasis. Understanding of the biophysical mechanisms that promote tumor metastasis may contribute successful therapeutic approaches in the fight against cancer.

Computational models of cancer cell transport through the microcirculation

The transport of cancerous cells through the microcirculation during metastatic spread encompasses several interdependent steps that are not fully understood. Computational models which resolve the cellular-scale dynamics of complex microcirculatory flows offer considerable potential to yield needed insights into the spread of cancer as a result of the level of detail that can be captured. In recent years, in silico methods have been developed that can accurately and efficiently model the circulatory flows of cancer and other biological cells. These computational methods are capable of resolving detailed fluid flow fields which transport cells through tortuous physiological geometries, as well as the deformation and interactions between cells, cell-to-endothelium interactions, and tumor cell aggregates, all of which play important roles in metastatic spread. Such models can provide a powerful complement to experimental works, and a promising approach to recapitulating the endogenous setting while maintaining control over parameters such as shear rate, cell deformability, and the strength of adhesive binding to better understand tumor cell transport. In this review, we present an overview of computational models that have been developed for modeling cancer cells in the microcirculation, including insights they have provided into cell transport phenomena.

Immersed boundary method for multiphase transport phenomena

Multiphase flows with momentum, heat, and mass transfer exist widely in a variety of industrial applications. With the rapid development of numerical algorithms and computer capacity, advanced numerical simulation has become a promising tool in investigating multiphase transport problems. Immersed boundary (IB) method has recently emerged as such a popular interface capturing method for efficient simulations of multiphase flows, and significant achievements have been obtained. In this review, we attempt to give an overview of recent progresses on IB method for multiphase transport phenomena. Firstly, the governing equations, the basic ideas, and different boundary conditions for the IB methods are introduced. This is followed by numerical strategies, from which the IB methods are classified into two types, namely the artificial boundary method and the authentic boundary method. Discussions on the implementation of various boundary conditions at the interphase surface with momentum, heat, and mass transfer for different IB methods are then presented, together with a summary. Then, the state-of-the-art applications of IB methods to multiphase flows, including the isothermal flows, the heat transfer flows, and the mass transfer problems are outlined, with particular emphasis on the latter two topics. Finally, the conclusions and future challenges are identified.

Impact of Force Function Formulations on the Numerical Simulation of Centre-Based Models

Centre-based or cell-centre models are a framework for the computational study of multicellular systems with widespread use in cancer modelling and computational developmental biology. At the core of these models are the numerical method used to update cell positions and the force functions that encode the pairwise mechanical interactions of cells. For the latter, there are multiple choices that could potentially affect both the biological behaviour captured, and the robustness and efficiency of simulation. For example, available open-source software implementations of centre-based models rely on different force functions for their default behaviour and it is not straightforward for a modeller to know if these are interchangeable. Our study addresses this problem and contributes to the understanding of the potential and limitations of three popular force functions from a numerical perspective. We show empirically that choosing the force parameters such that the relaxation time for two cells after cell division is consistent between different force functions results in good agreement of the population radius of a two-dimensional monolayer relaxing mechanically after intense cell proliferation. Furthermore, we report that numerical stability is not sufficient to prevent unphysical cell trajectories following cell division, and consequently, that too large time steps can cause geometrical differences at the population level.

A Review of Cell-Based Computational Modeling in Cancer Biology

Cancer biology involves complex, dynamic interactions between cancer cells and their tissue microenvironments. Single-cell effects are critical drivers of clinical progression. Chemical and mechanical communication between tumor and stromal cells can co-opt normal physiologic processes to promote growth and invasion. Cancer cell heterogeneity increases cancer’s ability to test strategies to adapt to microenvironmental stresses. Hypoxia and treatment can select for cancer stem cells and drive invasion and resistance. Cell-based computational models (also known as discrete models, agent-based models, or individual-based models) simulate individual cells as they interact in virtual tissues, which allows us to explore how single-cell behaviors lead to the dynamics we observe and work to control in cancer systems. In this review, we introduce the broad range of techniques available for cell-based computational modeling. The approaches can range from highly detailed models of just a few cells and their morphologies to millions of simpler cells in three-dimensional tissues. Modeling individual cells allows us to directly translate biologic observations into simulation rules. In many cases, individual cell agents include molecular-scale models. Most models also simulate the transport of oxygen, drugs, and growth factors, which allow us to link cancer development to microenvironmental conditions. We illustrate these methods with examples drawn from cancer hypoxia, angiogenesis, invasion, stem cells, and immunosurveillance. An ecosystem of interoperable cell-based simulation tools is emerging at a time when cloud computing resources make software easier to access and supercomputing resources make large-scale simulation studies possible. As the field develops, we anticipate that high-throughput simulation studies will allow us to rapidly explore the space of biologic possibilities, prescreen new therapeutic strategies, and even re-engineer tumor and stromal cells to bring cancer systems under control.

A quantitative high-resolution computational mechanics cell model for growing and regenerating tissues

Mathematical models are increasingly designed to guide experiments in biology, biotechnology, as well as to assist in medical decision making. They are in particular important to understand emergent collective cell behavior. For this purpose, the models, despite still abstractions of reality, need to be quantitative in all aspects relevant for the question of interest. This paper considers as showcase example the regeneration of liver after drug-induced depletion of hepatocytes, in which the surviving and dividing hepatocytes must squeeze in between the blood vessels of a network to refill the emerged lesions. Here, the cells' response to mechanical stress might significantly impact the regeneration process. We present a 3D high-resolution cell-based model integrating information from measurements in order to obtain a refined and quantitative understanding of the impact of cell-biomechanical effects on the closure of drug-induced lesions in liver. Our model represents each cell individually and is constructed by a discrete, physically scalable network of viscoelastic elements, capable of mimicking realistic cell deformation and supplying information at subcellular scales. The cells have the capability to migrate, grow, and divide, and the nature and parameters of their mechanical elements can be inferred from comparisons with optical stretcher experiments. Due to triangulation of the cell surface, interactions of cells with arbitrarily shaped (triangulated) structures such as blood vessels can be captured naturally. Comparing our simulations with those of so-called center-based models, in which cells have a largely rigid shape and forces are exerted between cell centers, we find that the migration forces a cell needs to exert on its environment to close a tissue lesion, is much smaller than predicted by center-based models. To stress generality of the approach, the liver simulations were complemented by monolayer and multicellular spheroid growth simulations. In summary, our model can give quantitative insight in many tissue organization processes, permits hypothesis testing in silico, and guide experiments in situations in which cell mechanics is considered important.

Morphophenotypic classification of tumor organoids as an indicator of drug exposure and penetration potential

The dynamics of tumor progression is driven by multiple factors, which can be exogenous to the tumor (microenvironment) or intrinsic (genetic, epigenetic or due to intercellular interactions). While tumor heterogeneity has been extensively studied on the level of cell genetic profiles or cellular composition, tumor morphological diversity has not been given as much attention. The limited analysis of tumor morphophenotypes may be attributed to the lack of accurate models, both experimental and computational, capable of capturing changes in tumor morphology with fine levels of spatial detail. Using a three-dimensional, agent-based, lattice-free computational model, we generated a library of multicellular tumor organoids, the experimental analogues of in vivo tumors. By varying three biologically relevant parameters-cell radius, cell division age and cell sensitivity to contact inhibition, we showed that tumor organoids with similar growth dynamics can express distinct morphologies and possess diverse cellular compositions. Taking advantage of the high-resolution of computational modeling, we applied the quantitative measures of compactness and accessible surface area, concepts that originated from the structural biology of proteins. Based on these analyses, we demonstrated that tumor organoids with similar sizes may differ in features associated with drug effectiveness, such as potential exposure to the drug or the extent of drug penetration. Both these characteristics might lead to major differences in tumor organoid's response to therapy. This indicates that therapeutic protocols should not be based solely on tumor size, but take into account additional tumor features, such as their morphology or cellular packing density.

Towards personalized computational oncology: from spatial models of tumour spheroids, to organoids, to tissues

A main goal of mathematical and computational oncology is to develop quantitative tools to determine the most effective therapies for each individual patient. This involves predicting the right drug to be administered at the right time and at the right dose. Such an approach is known as precision medicine. Mathematical modelling can play an invaluable role in the development of such therapeutic strategies, since it allows for relatively fast, efficient and inexpensive simulations of a large number of treatment schedules in order to find the most effective. This review is a survey of mathematical models that explicitly take into account the spatial architecture of three-dimensional tumours and address tumour development, progression and response to treatments. In particular, we discuss models of epithelial acini, multicellular spheroids, normal and tumour spheroids and organoids, and multi-component tissues. Our intent is to showcase how these in silico models can be applied to patient-specific data to assess which therapeutic strategies will be the most efficient. We also present the concept of virtual clinical trials that integrate standard-of-care patient data, medical imaging, organ-on-chip experiments and computational models to determine personalized medical treatment strategies.

Balance of mechanical forces drives endothelial gap formation and may facilitate cancer and immune-cell extravasation

The formation of gaps in the endothelium is a crucial process underlying both cancer and immune cell extravasation, contributing to the functioning of the immune system during infection, the unfavorable development of chronic inflammation and tumor metastasis. Here, we present a stochastic-mechanical multiscale model of an endothelial cell monolayer and show that the dynamic nature of the endothelium leads to spontaneous gap formation, even without intervention from the transmigrating cells. These gaps preferentially appear at the vertices between three endothelial cells, as opposed to the border between two cells. We quantify the frequency and lifetime of these gaps, and validate our predictions experimentally. Interestingly, we find experimentally that cancer cells also preferentially extravasate at vertices, even when they first arrest on borders. This suggests that extravasating cells, rather than initially signaling to the endothelium, might exploit the autonomously forming gaps in the endothelium to initiate transmigration.

Tumor growth and calcification in evolving microenvironmental geometries

In this paper, we apply the diffuse domain framework developed in Chen and Lowengrub (Tumor growth in complex, evolving microenvironmental geometries: A diffuse domain approach, J. Theor. Biol. 361 (2014) 14-30) to study the effects of a deformable basement membrane (BM) on the growth of a tumor in a confined, ductal geometry, such as ductal carcinoma in situ (DCIS). We use a continuum model of tumor microcalcification and investigate the tumor extent beyond the microcalcification. In order to solve the governing equations efficiently, we develop a stable nonlinear multigrid finite difference method. Two dimensional simulations are performed where the adhesion between tumor cells and the basement membrane is varied. Additional simulations considering the variation of duct radius and membrane stiffness are also conducted. The results demonstrate that enhanced membrane deformability promotes tumor growth and tumor calcification. When the duct radius is small, the cell-BM adhesion is weak or when the membrane is slightly deformed, the mammographic and pathologic tumor extents are linearly correlated, as predicted by Macklin et al. (J. Theor. Biol. 301 (2012) 122-140) using an agent-based model that does not account for the deformability of the basement membrane and the active forces that the membrane imparts on the tumor cells. Interestingly, we predict that when the duct radius is large, there is strong cell-BM adhesion or the membrane is highly deformed, the extents of the mammographic and pathologic tumors are instead quadratically correlated. The simulations can help surgeons to measure DCIS surgical margins while removing less non-cancerous tissue, and can improve targeting of intra- and post-operative radiotherapy.

Scalable population-level modelling of biological cells incorporating mechanics and kinetics in continuous time

The processes taking place inside the living cell are now understood to the point where predictive computational models can be used to gain detailed understanding of important biological phenomena. A key challenge is to extrapolate this detailed knowledge of the individual cell to be able to explain at the population level how cells interact and respond with each other and their environment. In particular, the goal is to understand how organisms develop, maintain and repair functional tissues and organs. In this paper, we propose a novel computational framework for modelling populations of interacting cells. Our framework incorporates mechanistic, constitutive descriptions of biomechanical properties of the cell population, and uses a coarse-graining approach to derive individual rate laws that enable propagation of the population through time. Thanks to its multiscale nature, the resulting simulation algorithm is extremely scalable and highly efficient. As highlighted in our computational examples, the framework is also very flexible and may straightforwardly be coupled with continuous-time descriptions of biochemical signalling within, and between, individual cells.

A general framework dedicated to computational morphogenesis Part I - Constitutive equations

In order to understand living organisms, considerable experimental efforts and resources have been devoted to correlate genes and their expressions with cell, tissue, organ and whole organisms' phenotypes. This data driven approach to knowledge discovery has led to many breakthrough in our understanding of healthy and diseased states, and is paving the way to improve the diagnosis and treatment of diseases. Complementary to this data-driven approach, computational models of biological systems based on first principles have been developed in order to deepen our understanding of the multi-scale dynamics that drives normal and pathological biological functions. In this paper we describe the biological, physical and mathematical concepts that led to the design of a Computational Morphogenesis (CM) platform baptized Generic Modeling and Simulating Platform (GMSP). Its role is to generate realistic 3D multi-scale biological tissues from virtual stem cells and the intended target applications include in virtuo studies of normal and abnormal tissue (re)generation as well as the development of complex diseases such as carcinogenesis. At all space-scales of interest, biological agents interact with each other via biochemical, bioelectrical, and mechanical fields that operate in concert during embryogenesis, growth and adult life. The spatio-temporal dependencies of these fields can be modeled by physics-based constitutive equations that we propose to examine in relation to the landmark biological events that occur during embryogenesis.

Cell signaling heterogeneity is modulated by both cell-intrinsic and -extrinsic mechanisms: An integrated approach to understanding targeted therapy

During the last decade, our understanding of cancer cell signaling networks has significantly improved, leading to the development of various targeted therapies that have elicited profound but, unfortunately, short-lived responses. This is, in part, due to the fact that these targeted therapies ignore context and average out heterogeneity. Here, we present a mathematical framework that addresses the impact of signaling heterogeneity on targeted therapy outcomes. We employ a simplified oncogenic rat sarcoma (RAS)-driven mitogen-activated protein kinase (MAPK) and phosphoinositide 3-kinase-protein kinase B (PI3K-AKT) signaling pathway in lung cancer as an experimental model system and develop a network model of the pathway. We measure how inhibition of the pathway modulates protein phosphorylation as well as cell viability under different microenvironmental conditions. Training the model on this data using Monte Carlo simulation results in a suite of in silico cells whose relative protein activities and cell viability match experimental observation. The calibrated model predicts distributional responses to kinase inhibitors and suggests drug resistance mechanisms that can be exploited in drug combination strategies. The suggested combination strategies are validated using in vitro experimental data. The validated in silico cells are further interrogated through an unsupervised clustering analysis and then integrated into a mathematical model of tumor growth in a homogeneous and resource-limited microenvironment. We assess posttreatment heterogeneity and predict vast differences across treatments with similar efficacy, further emphasizing that heterogeneity should modulate treatment strategies. The signaling model is also integrated into a hybrid cellular automata (HCA) model of tumor growth in a spatially heterogeneous microenvironment. As a proof of concept, we simulate tumor responses to targeted therapies in a spatially segregated tissue structure containing tumor and stroma (derived from patient tissue) and predict complex cell signaling responses that suggest a novel combination treatment strategy.

A signaling pathway is a network of molecules in a cell that is typically initiated by stimuli (e.g., microenvironmental cues) acting on receptors and internal signaling molecules to determine cell fate. Signaling pathways in cancer cells are different from those in normal cells, and this difference helps cancer cells to grow and thrive indefinitely. Drugs that target the aberrant signaling pathways in cancer cells (often referred to as targeted therapy) are promising for improving treatment outcomes of many different cancers in patients. However, most patients eventually develop resistance to these drugs. Resistance may already be present in the tumor or may emerge via mutation or via microenvironmental mediation. Tumor heterogeneity, which is characterized by subtle or dramatic differences among tumor cells, plays a key role in the development of drug resistance. Some tumor cells respond well to therapy, while others may adapt to the stress induced by the drug within the microenvironment. Moreover, removal of drug-sensitive cells may result in the competitive release of drug-resistant cells. Here, we present mathematical models to assess the impact of heterogeneity in signaling pathways within tumor cells on the outcomes of targeted therapy. We consider a simplified version of two well-known signaling pathways that modulate the growth of lung cancer cells. By using different targeted therapies, we quantify the effect of pathway inhibition on protein activity and cell viability and developed a mathematical model of the network, which is trained to reproduce these data and to develop a panel of heterogeneous in silico cells. The model predicts potential mechanisms of drug resistance and proposes combination therapies that are effective across the panel. We validate these combination therapies experimentally using the lung cancer cells and integrated the in silico cells into a computational lung tissue model that explicitly captures the microenvironment of lung cancer. Our results suggest that heterogeneity in both the tumor and microenvironment impacts treatment response in different ways and suggest a novel combination therapy for a better response.

Hybrid cell-centred/vertex model for multicellular systems with equilibrium-preserving remodelling

We present a hybrid cell-centred/vertex model for mechanically simulating planar cellular monolayers undergoing cell reorganisation. Cell centres are represented by a triangular nodal network, while the cell boundaries are formed by an associated vertex network. The two networks are coupled through a kinematic constraint which we allow to relax progressively. Special attention is paid to the change of cell-cell connectivity due to cell reorganisation or remodelling events. We handle these situations by using a variable resting length and applying an Equilibrium-Preserving Mapping on the new connectivity, which computes a new set of resting lengths that preserve nodal and vertex equilibrium. We illustrate the properties of the model by simulating monolayers subjected to imposed extension and during a wound healing process. The evolution of forces and the Equilibrium-Preserving Mapping are analysed during the remodelling events. As a by-product, the proposed technique enables to recover fully vertex or fully cell-centred models in a seamless manner by modifying a numerical parameter of the model.

Agent-Based Modeling of Cancer Stem Cell Driven Solid Tumor Growth

Computational modeling of tumor growth has become an invaluable tool to simulate complex cell-cell interactions and emerging population-level dynamics. Agent-based models are commonly used to describe the behavior and interaction of individual cells in different environments. Behavioral rules can be informed and calibrated by in vitro assays, and emerging population-level dynamics may be validated with both in vitro and in vivo experiments. Here, we describe the design and implementation of a lattice-based agent-based model of cancer stem cell driven tumor growth.

Simulation of Morphogen and Tissue Dynamics

Morphogenesis, the process by which an adult organism emerges from a single cell, has fascinated humans for a long time. Modeling this process can provide novel insights into development and the principles that orchestrate the developmental processes. This chapter focuses on the mathematical description and numerical simulation of developmental processes. In particular, we discuss the mathematical representation of morphogen and tissue dynamics on static and growing domains, as well as the corresponding tissue mechanics. In addition, we give an overview of numerical methods that are routinely used to solve the resulting systems of partial differential equations. These include the finite element method and the Lattice Boltzmann method for the discretization as well as the arbitrary Lagrangian-Eulerian method and the Diffuse-Domain method to numerically treat deforming domains.

Dynamic cellular finite-element method for modelling large-scale cell migration and proliferation under the control of mechanical and biochemical cues: a study of re-epithelialization

Computational modelling of cells can reveal insight into the mechanisms of the important processes of tissue development. However, current cell models have limitations and are challenged to model detailed changes in cellular shapes and physical mechanics when thousands of migrating and interacting cells need to be modelled. Here we describe a novel dynamic cellular finite-element model (DyCelFEM), which accounts for changes in cellular shapes and mechanics. It also models the full range of cell motion, from movements of individual cells to collective cell migrations. The transmission of mechanical forces regulated by intercellular adhesions and their ruptures are also accounted for. Intra-cellular protein signalling networks controlling cell behaviours are embedded in individual cells. We employ DyCelFEM to examine specific effects of biochemical and mechanical cues in regulating cell migration and proliferation, and in controlling tissue patterning using a simplified re-epithelialization model of wound tissue. Our results suggest that biochemical cues are better at guiding cell migration with improved directionality and persistence, while mechanical cues are better at coordinating collective cell migration. Overall, DyCelFEM can be used to study developmental processes when a large population of migrating cells under mechanical and biochemical controls experience complex changes in cell shapes and mechanics.

Mechanocellular models of epithelial morphogenesis

Embryonic epithelia achieve complex morphogenetic movements, including in-plane reshaping, bending and folding, through the coordinated action and rearrangement of individual cells. Technical advances in molecular and live-imaging studies of epithelial dynamics provide a very real opportunity to understand how cell-level processes facilitate these large-scale tissue rearrangements. However, the large datasets that we are now able to generate require careful interpretation. In combination with experimental approaches, computational modelling allows us to challenge and refine our current understanding of epithelial morphogenesis and to explore experimentally intractable questions. To this end, a variety of cell-based modelling approaches have been developed to describe cell-cell mechanical interactions, ranging from vertex and 'finite-element' models that approximate each cell geometrically by a polygon representing the cell's membrane, to immersed boundary and subcellular element models that allow for more arbitrary cell shapes. Here, we review how these models have been used to provide insights into epithelial morphogenesis and describe how such models could help future efforts to decipher the forces and mechanical and biochemical feedbacks that guide cell and tissue-level behaviour. In addition, we discuss current challenges associated with using computational models of morphogenetic processes in a quantitative and predictive way.This article is part of the themed issue 'Systems morphodynamics: understanding the development of tissue hardware'.

Simulation framework for generating intratumor heterogeneity patterns in a cancer cell population

As cancer cell populations evolve, they accumulate a number of somatic mutations, resulting in heterogeneous subclones in the final tumor. Understanding the mechanisms that produce intratumor heterogeneity is important for selecting the best treatment. Although some studies have involved intratumor heterogeneity simulations, their model settings differed substantially. Thus, only limited conditions were explored in each. Herein, we developed a general framework for simulating intratumor heterogeneity patterns and a simulator (tumopp). Tumopp offers many setting options so that simulations can be carried out under various settings. Setting options include how the cell division rate is determined, how daughter cells are placed, and how driver mutations are treated. Furthermore, to account for the cell cycle, we introduced a gamma function for the waiting time involved in cell division. Tumopp also allows simulations in a hexagonal lattice, in addition to a regular lattice that has been used in previous simulation studies. A hexagonal lattice produces a more biologically reasonable space than a regular lattice. Using tumopp, we investigated how model settings affect the growth curve and intratumor heterogeneity pattern. It was found that, even under neutrality (with no driver mutations), tumopp produced dramatically variable patterns of intratumor heterogeneity and tumor morphology, from tumors in which cells with different genetic background are well intermixed to irregular shapes of tumors with a cluster of closely related cells. This result suggests a caveat in analyzing intratumor heterogeneity with simulations with limited settings, and tumopp will be useful to explore intratumor heterogeneity patterns in various conditions.

Are Tumor Cell Lineages Solely Shaped by Mechanical Forces?

This paper investigates cell proliferation dynamics in small tumor cell aggregates using an individual-based model (IBM). The simulation model is designed to study the morphology of the cell population and of the cell lineages as well as the impact of the orientation of the division plane on this morphology. Our IBM model is based on the hypothesis that cells are incompressible objects that grow in size and divide once a threshold size is reached, and that newly born cell adhere to the existing cell cluster. We performed comparisons between the simulation model and experimental data by using several statistical indicators. The results suggest that the emergence of particular morphologies can be explained by simple mechanical interactions.

An agent-based model of triple-negative breast cancer: the interplay between chemokine receptor CCR5 expression, cancer stem cells, and hypoxia

Background

Triple-negative breast cancer lacks estrogen, progesterone, and HER2 receptors and is thus not possible to treat with targeted therapies for these receptors. Therefore, a greater understanding of triple-negative breast cancer is necessary for the treatment of this cancer type. In previous work from our laboratory, we found that chemokine ligand-receptor CCL5-CCR5 axis is important for the metastasis of human triple-negative breast cancer cell MDA-MB-231 to the lymph nodes and lungs, in a mouse xenograft model. We collected relevant experimental data from our and other laboratories for numbers of cancer stem cells, numbers of CCR5+ cells, and cell migration rates for different breast cancer cell lines and different experimental conditions.

Results

Using these experimental data we developed an in silico agent-based model of triple-negative breast cancer that considers surface receptor CCR5-high and CCR5-low cells and breast cancer stem cells, to predict the tumor growth rate and spatio-temporal distribution of cells in primary tumors. We find that high cancer stem cell percentages greatly increase tumor growth. We find that anti-stem cell treatment decreases tumor growth but may not lead to dormancy unless all stem cells get eliminated. We further find that hypoxia increases overall tumor growth and treatment with a CCR5 inhibitor maraviroc slightly decreases overall tumor growth. We also characterize 3D shapes of solid and invasive tumors using several shape metrics.

Conclusions

Breast cancer stem cells and CCR5+ cells affect the overall growth and morphology of breast tumors. In silico drug treatments demonstrate limited efficacy of incomplete inhibition of cancer stem cells after which tumor growth recurs, and CCR5 inhibition causes only a slight reduction in tumor growth.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-017-0445-x) contains supplementary material, which is available to authorized users.

Homogeneous pancreatic cancer spheroids mimic growth pattern of circulating tumor cell clusters and macrometastases: displaying heterogeneity and crater-like structure on inner layer

PURPOSE

Pancreatic cancer 3D in vitro models including multicellular tumor spheroid (MCTS), single cell-derived tumor spheroid (SCTS), tissue-derived tumor spheroid, and organotypic models provided powerful platforms to mimic in vivo tumor. Recent work supports that circulating tumor cell (CTC) clusters are more efficient in metastasis seeding than single CTCs. The purpose of this study is to establish 3D culture models which can mimic single CTC, monoclonal CTC clusters, and the expansion of macrometastases.

METHODS

Seven pancreatic ductal adenocarcinoma cell lines were used to establish MCTS and SCTS using hanging drop and ultra-low attachment plates. Spheroid immunofluorescence staining, spheroid formation assay, immunoblotting, and literature review were performed to investigate molecular biomarkers and the morphological characteristics of pancreatic tumor spheroids.

RESULTS

Single cells experienced different growth patterns to form SCTS, like signet ring-like cells, blastula-like structures, and solid core spheroids. However, golf ball-like hollow spheroids could also be detected, especially when DanG and Capan-1 cells were cultivated with fibroblast-conditioned medium (p < 0.05). The size of golf ball-like hollow spheroids hardly grew after getting matured. Only DanG and Capan-1 could establish SCTS- and MCTS-derived hollow spheroids using hanging drop plates and ultra-low attachment plates. Other PDA cell lines could also establish tumor spheroid with hanging drop plates by adding methylated cellulose. Tumor spheroids derived from pancreatic cancer cell line DanG possessed asymmetrically distributed proliferation center, immune-checkpoint properties. ß-catenin, Ki-67, and F-actin were active surrounding the crater-like structure distributing on the inner layer of viable rim cover of the spheroids, which was relevant to well-differentiated tumor cells.

CONCLUSIONS

It is possible to establish 3D CTC cluster models from homogenous PDA cell lines using hanging drop and ultra-low attachment plates. PDA cell line displays its own intrinsic properties or heterogeneity. The mechanism of formation of the crater-like structure as well as golf ball-like structure needs further exploration.

Multi-scale computational study of the mechanical regulation of cell mitotic rounding in epithelia

Mitotic rounding during cell division is critical for preventing daughter cells from inheriting an abnormal number of chromosomes, a condition that occurs frequently in cancer cells. Cells must significantly expand their apical area and transition from a polygonal to circular apical shape to achieve robust mitotic rounding in epithelial tissues, which is where most cancers initiate. However, how cells mechanically regulate robust mitotic rounding within packed tissues is unknown. Here, we analyze mitotic rounding using a newly developed multi-scale subcellular element computational model that is calibrated using experimental data. Novel biologically relevant features of the model include separate representations of the sub-cellular components including the apical membrane and cytoplasm of the cell at the tissue scale level as well as detailed description of cell properties during mitotic rounding. Regression analysis of predictive model simulation results reveals the relative contributions of osmotic pressure, cell-cell adhesion and cortical stiffness to mitotic rounding. Mitotic area expansion is largely driven by regulation of cytoplasmic pressure. Surprisingly, mitotic shape roundness within physiological ranges is most sensitive to variation in cell-cell adhesivity and stiffness. An understanding of how perturbed mechanical properties impact mitotic rounding has important potential implications on, amongst others, how tumors progressively become more genetically unstable due to increased chromosomal aneuploidy and more aggressive.

Dynamics of cancerous tissue correlates with invasiveness

Two of the classical hallmarks of cancer are uncontrolled cell division and tissue invasion, which turn the disease into a systemic, life-threatening condition. Although both processes are studied, a clear correlation between cell division and motility of cancer cells has not been described previously. Here, we experimentally characterize the dynamics of invasive and non-invasive breast cancer tissues using human and murine model systems. The intrinsic tissue velocities, as well as the divergence and vorticity around a dividing cell correlate strongly with the invasive potential of the tissue, thus showing a distinct correlation between tissue dynamics and aggressiveness. We formulate a model which treats the tissue as a visco-elastic continuum. This model provides a valid reproduction of the cancerous tissue dynamics, thus, biological signaling is not needed to explain the observed tissue dynamics. The model returns the characteristic force exerted by an invading cell and reveals a strong correlation between force and invasiveness of breast cancer cells, thus pinpointing the importance of mechanics for cancer invasion.

The role of myosin II in glioma invasion: A mathematical model

Gliomas are malignant tumors that are commonly observed in primary brain cancer. Glioma cells migrate through a dense network of normal cells in microenvironment and spread long distances within brain. In this paper we present a two-dimensional multiscale model in which a glioma cell is surrounded by normal cells and its migration is controlled by cell-mechanical components in the microenvironment via the regulation of myosin II in response to chemoattractants. Our simulation results show that the myosin II plays a key role in the deformation of the cell nucleus as the glioma cell passes through the narrow intercellular space smaller than its nuclear diameter. We also demonstrate that the coordination of biochemical and mechanical components within the cell enables a glioma cell to take the mode of amoeboid migration. This study sheds lights on the understanding of glioma infiltration through the narrow intercellular spaces and may provide a potential approach for the development of anti-invasion strategies via the injection of chemoattractants for localization.

Comparing individual-based approaches to modelling the self-organization of multicellular tissues

The coordinated behaviour of populations of cells plays a central role in tissue growth and renewal. Cells react to their microenvironment by modulating processes such as movement, growth and proliferation, and signalling. Alongside experimental studies, computational models offer a useful means by which to investigate these processes. To this end a variety of cell-based modelling approaches have been developed, ranging from lattice-based cellular automata to lattice-free models that treat cells as point-like particles or extended shapes. However, it remains unclear how these approaches compare when applied to the same biological problem, and what differences in behaviour are due to different model assumptions and abstractions. Here, we exploit the availability of an implementation of five popular cell-based modelling approaches within a consistent computational framework, Chaste (http://www.cs.ox.ac.uk/chaste). This framework allows one to easily change constitutive assumptions within these models. In each case we provide full details of all technical aspects of our model implementations. We compare model implementations using four case studies, chosen to reflect the key cellular processes of proliferation, adhesion, and short- and long-range signalling. These case studies demonstrate the applicability of each model and provide a guide for model usage.

Simulating the multicellular homeostasis with a cell-based discrete receptor dynamics model: The non-mutational origin of cancer and aging

The purpose of the study is to investigate the multicellular homeostasis in epithelial tissues over very large timescales. Inspired by the receptor dynamics of IBCell model proposed by Rejniak et al. an on-grid agent-based model for multicellular system is constructed. Instead of observing the multicellular architectural morphologies, the diversity of homeostatic states is quantitatively analyzed through a substantial number of simulations by measuring three new order parameters, the phenotypic population structure, the average proliferation age and the relaxation time to stable homeostasis. Nearby the interfaces of distinct homeostatic phases in 3D phase diagrams of the three order parameters, intermediate quasi-stable phases of slow dynamics that features quasi-stability with a large spectrum of relaxation timescales are found. A further exploration on the static and dynamic correlations among the three order parameters reveals that the quasi-stable phases evolve towards two terminations, tumorigenesis and degeneration, which are respectively accompanied by rejuvenation and aging. With the exclusion of the environmental impact and the mutational strategies, the results imply that cancer and aging may share the non-mutational origin in the intrinsic slow dynamics of the multicellular systems.

Prediction of Growth Factor-Dependent Cleft Formation During Branching Morphogenesis Using A Dynamic Graph-Based Growth Model

This study considers the problem of describing and predicting cleft formation during the early stages of branching morphogenesis in mouse submandibular salivary glands (SMG) under the influence of varied concentrations of epidermal growth factors (EGF). Given a time-lapse video of a growing SMG, first we build a descriptive model that captures the underlying biological process and quantifies the ground truth. Tissue-scale (global) and morphological features related to regions of interest (local features) are used to characterize the biological ground truth. Second, we devise a predictive growth model that simulates EGF-modulated branching morphogenesis using a dynamic graph algorithm, which is driven by biological parameters such as EGF concentration, mitosis rate, and cleft progression rate. Given the initial configuration of the SMG, the evolution of the dynamic graph predicts the cleft formation, while maintaining the local structural characteristics of the SMG. We determined that higher EGF concentrations cause the formation of higher number of buds and comparatively shallow cleft depths. Third, we compared the prediction accuracy of our model to the Glazier-Graner-Hogeweg (GGH) model, an on-lattice Monte-Carlo simulation model, under a specific energy function parameter set that allows new rounds of de novo cleft formation. The results demonstrate that the dynamic graph model yields comparable simulations of gland growth to that of the GGH model with a significantly lower computational complexity. Fourth, we enhanced this model to predict the SMG morphology for an EGF concentration without the assistance of a ground truth time-lapse biological video data; this is a substantial benefit of our model over other similar models that are guided and terminated by information regarding the final SMG morphology. Hence, our model is suitable for testing the impact of different biological parameters involved with the process of branching morphogenesis in silico, while reducing the requirement of in vivo experiments.

Circulating Tumor Cells: When a Solid Tumor Meets a Fluid Microenvironment

Solid tumor dissemination from the primary site to the sites of metastasis involves tumor cell transport through the blood or lymph circulation systems. Once the tumor cells enter the bloodstream, they encounter a new hostile microenvironment. The cells must withstand hemodynamic forces and overcome the effects of fluid shear. The cells are exposed to immunological signaling insults from leukocytes, to collisions with erythrocytes, and to interactions with platelets or macrophages. Finally, the cells need to attach to the blood vessel walls and extravasate to the surrounding stroma to form tumor metastases. Although only a small fraction of invasive cells is able to complete the metastatic process, most cancer-related deaths are the result of tumor metastasis. Thus, investigating the intracellular properties of circulating tumor cells and the extracellular conditions that allow the tumor cells to survive and thrive in this microenvironment is of vital interest. In this chapter, we discuss the intravascular microenvironment that the circulating tumor cells must endure. We summarize the current experimental and computational literature on tumor cells in the circulation system. We also illustrate various aspects of the intravascular transport of circulating tumor cells using a mathematical model based on immersed boundary principles.

Three-dimensional simulations of the cell growth and cytokinesis using the immersed boundary method

In this paper, we present a three-dimensional immersed boundary method to simulate the eukaryotic cell growth and cytokinesis. The proposed model and numerical method are a non-trivial three-dimensional extension of the previous work (Li et al., 2012). Unstructured triangular meshes are employed to discretize the cell membrane. The nodes of the surface mesh constitute a set of Lagrangian control points used to track the motion of the cell. A surface remeshing algorithm is applied to prevent mesh distortion during evolution. We also use a volume-conserving algorithm to maintain the mass of cells in cytokinesis. The ability of the proposed method to simulate cell growth and division processes is numerically demonstrated.

Mechanical model of geometric cell and topological algorithm for cell dynamics from single-cell to formation of monolayered tissues with pattern

Geometric and mechanical properties of individual cells and interactions among neighboring cells are the basis of formation of tissue patterns. Understanding the complex interplay of cells is essential for gaining insight into embryogenesis, tissue development, and other emerging behavior. Here we describe a cell model and an efficient geometric algorithm for studying the dynamic process of tissue formation in 2D (e.g. epithelial tissues). Our approach improves upon previous methods by incorporating properties of individual cells as well as detailed description of the dynamic growth process, with all topological changes accounted for. Cell size, shape, and division plane orientation are modeled realistically. In addition, cell birth, cell growth, cell shrinkage, cell death, cell division, cell collision, and cell rearrangements are now fully accounted for. Different models of cell-cell interactions, such as lateral inhibition during the process of growth, can be studied in detail. Cellular pattern formation for monolayered tissues from arbitrary initial conditions, including that of a single cell, can also be studied in detail. Computational efficiency is achieved through the employment of a special data structure that ensures access to neighboring cells in constant time, without additional space requirement. We have successfully generated tissues consisting of more than 20,000 cells starting from 2 cells within 1 hour. We show that our model can be used to study embryogenesis, tissue fusion, and cell apoptosis. We give detailed study of the classical developmental process of bristle formation on the epidermis of D. melanogaster and the fundamental problem of homeostatic size control in epithelial tissues. Simulation results reveal significant roles of solubility of secreted factors in both the bristle formation and the homeostatic control of tissue size. Our method can be used to study broad problems in monolayered tissue formation. Our software is publicly available.

LBIBCell: a cell-based simulation environment for morphogenetic problems

MOTIVATION

The simulation of morphogenetic problems requires the simultaneous and coupled simulation of signalling and tissue dynamics. A cellular resolution of the tissue domain is important to adequately describe the impact of cell-based events, such as cell division, cell-cell interactions and spatially restricted signalling events. A tightly coupled cell-based mechano-regulatory simulation tool is therefore required.

RESULTS

We developed an open-source software framework for morphogenetic problems. The environment offers core functionalities for the tissue and signalling models. In addition, the software offers great flexibility to add custom extensions and biologically motivated processes. Cells are represented as highly resolved, massless elastic polygons; the viscous properties of the tissue are modelled by a Newtonian fluid. The Immersed Boundary method is used to model the interaction between the viscous and elastic properties of the cells, thus extending on the IBCell model. The fluid and signalling processes are solved using the Lattice Boltzmann method. As application examples we simulate signalling-dependent tissue dynamics.

AVAILABILITY AND IMPLEMENTATION

The documentation and source code are available on http://tanakas.bitbucket.org/lbibcell/index.html

Multiscale Modeling of Cellular Epigenetic States: Stochasticity in Molecular Networks, Chromatin Folding in Cell Nuclei, and Tissue Pattern Formation of Cells

Genome sequences provide the overall genetic blueprint of cells, but cells possessing the same genome can exhibit diverse phenotypes. There is a multitude of mechanisms controlling cellular epigenetic states and that dictate the behavior of cells. Among these, networks of interacting molecules, often under stochastic control, depending on the specific wirings of molecular components and the physiological conditions, can have a different landscape of cellular states. In addition, chromosome folding in three-dimensional space provides another important control mechanism for selective activation and repression of gene expression. Fully differentiated cells with different properties grow, divide, and interact through mechanical forces and communicate through signal transduction, resulting in the formation of complex tissue patterns. Developing quantitative models to study these multi-scale phenomena and to identify opportunities for improving human health requires development of theoretical models, algorithms, and computational tools. Here we review recent progress made in these important directions.

Simulating tissue morphogenesis and signaling

During embryonic development tissue morphogenesis and signaling are tightly coupled. It is therefore important to simulate both tissue morphogenesis and signaling simultaneously in in silico models of developmental processes. The resolution of the processes depends on the questions of interest. As part of this chapter we introduce different descriptions of tissue morphogenesi s. In the simplest approximation tissue is a continuous domain and tissue expansion is described according to a predefined function of time (and possibly space). In a slightly more advanced version the expansion speed and direction of the tissue may depend on a signaling variable that evolves on the domain. Both versions will be referred to as "prescribed growth." Alternatively tissue can be regarded as incompressible fluid and can be described with Navier-Stokes equations. Local cell expansion, proliferation, and death are then incorporated by a source term. In other applications the cell boundaries may be important and cell-based models must be introduced. Finally, cells may move within the tissue, a process best described by agent-based models.

A new model for cell division and migration with spontaneous topology changes

Tissue topology, in particular proliferating epithelium topology, is remarkably similar between various species. Understanding the mechanisms that result in the observed topologies is needed for better insight into the processes governing tissue formation. We present a two-dimensional single-cell based model for cell divisions and tissue growth. The model accounts for cell mechanics and allows cell migration. Cells do not have pre-existing shapes or topologies. Shape changes and local rearrangements occur naturally as a response to the evolving cellular environment and cell-cell interactions. We show that the commonly observed tissue topologies arise spontaneously from this model. We consider different cellular rearrangements that accompany tissue growth and study their effects on tissue topology.