Physical confinement and cell proximity increase cell migration rates and invasiveness: A mathematical model of cancer cell invasion through flexible channels.

Cancer cell migration between different body parts is the driving force behind cancer metastasis, which is the main cause of mortality of patients. Migration of cancer cells often proceeds by penetration through narrow cavities in locally stiff, yet flexible tissues. In our previous work, we developed a model for cell geometry evolution during invasion, which we extend here to investigate whether leader and follower (cancer) cells that only interact mechanically can benefit from sequential transmigration through narrow micro-channels and cavities. We consider two cases of cells sequentially migrating through a flexible channel: leader and follower cells being closely adjacent or distant. Using Wilcoxon's signed-rank test on the data collected from Monte Carlo simulations, we conclude that the modelled transmigration speed for the follower cell is significantly larger than for the leader cell when cells are distant, i.e. follower cells transmigrate after the leader has completed the crossing. Furthermore, it appears that there exists an optimum with respect to the width of the channel such that cell moves fastest. On the other hand, in the case of closely adjacent cells, effectively performing collective migration, the leader cell moves 12% faster since the follower cell pushes it. This work shows that mechanical interactions between cells can increase the net transmigration speed of cancer cells, resulting in increased invasiveness. In other words, interaction between cancer cells can accelerate metastatic invasion.

Stability of a two-dimensional biomorphoelastic model for post-burn contraction.

We consider the stability analysis of a two-dimensional model for post-burn contraction. The model is based on morphoelasticity for permanent deformations and combined with a chemical-biological model that incorporates cellular densities, collagen density, and the concentration of chemoattractants. We formulate stability conditions depending on the decay rate of signaling molecules for both the continuous partial differential equations-based problem and the (semi-)discrete representation. We analyze the difference and convergence between the resulting spatial eigenvalues from the continuous and semi-discrete problems.

Mathematical model of mechano-sensing and mechanically induced collective motility of cells on planar elastic substrates.

Cells mechanically interact with their environment to sense, for example, topography, elasticity and mechanical cues from other cells. Mechano-sensing has profound effects on cellular behaviour, including motility. The current study aims to develop a mathematical model of cellular mechano-sensing on planar elastic substrates and demonstrate the model's predictive capabilities for the motility of individual cells in a colony. In the model, a cell is assumed to transmit an adhesion force, derived from a dynamic focal adhesion integrin density, that locally deforms a substrate, and to sense substrate deformation originating from neighbouring cells. The substrate deformation from multiple cells is expressed as total strain energy density with a spatially varying gradient. The magnitude and direction of the gradient at the cell location define the cell motion. Cell-substrate friction, partial motion randomness, and cell death and division are included. The substrate deformation by a single cell and the motility of two cells are presented for several substrate elasticities and thicknesses. The collective motility of 25 cells on a uniform substrate mimicking the closure of a circular wound of 200 µm is predicted for deterministic and random motion. Cell motility on substrates with varying elasticity and thickness is explored for four cells and 15 cells, the latter again mimicking wound closure. Wound closure by 45 cells is used to demonstrate the simulation of cell death and division during migration. The mathematical model can adequately simulate the mechanically induced collective cell motility on planar elastic substrates. The model is suitable for extension to other cell and substrates shapes and the inclusion of chemotactic cues, offering the potential to complement in vitro and in vivo studies.

Sensitivity of a two-dimensional biomorphoelastic model for post-burn contraction.

We consider a two-dimensional biomorphoelastic model describing post-burn scar contraction. This model describes skin displacement and the development of the effective Eulerian strain in the tissue. Besides these mechanical components, signaling molecules, fibroblasts, myofibroblasts, and collagen also play a significant role in the model. We perform a sensitivity analysis for the independent parameters of the model and focus on the effects on features of the relative surface area and the total strain energy density. We conclude that the most sensitive parameters are the Poisson's ratio, the equilibrium collagen concentration, the contraction inhibitor constant, and the myofibroblast apoptosis rate. Next to these insights, we perform a sensitivity analysis where the proliferation rates of fibroblasts and myofibroblasts are not the same. The impact of this model adaptation is significant.

The Future of Burn Care From a Complexity Science Perspective.

Health care is undergoing a profound technological and digital transformation and has become increasingly complex. It is important for burns professionals and researchers to adapt to these developments which may require new ways of thinking and subsequent new strategies. As Einstein has put it: "We must learn to see the world anew." The relatively new scientific discipline "Complexity science" can give more direction to this and is the metaphorical open door that should not go unnoticed in view of the burn care of the future. Complexity science studies "why the whole is more than the sum of the parts." It studies how multiple separate components interact with each other and their environment and how these interactions lead to "behavior of the system." Biological systems are always part of smaller and larger systems and exhibit the behavior of adaptivity, hence the name complex adaptive systems. From the perspective of complexity science, a severe burn injury is an extreme disruption of the "human body system." But this disruption also applies to the systems at the organ and cellular levels. All these systems follow the principles of complex systems. Awareness of the scaling process at multilevel helps to understand and manage the complex situation when dealing with severe burn cases. This paper aims to create awareness of the concept of complexity and to demonstrate the value and possibilities of complexity science methods and tools for the future of burn care through examples from preclinical, clinical, and organizational perspectives in burn care.

A Bayesian finite-element trained machine learning approach for predicting post-burn contraction.

Burn injuries can decrease the quality of life of a patient tremendously, because of esthetic reasons and because of contractions that result from them. In severe case, skin contraction takes place at such a large extent that joint mobility of a patient is significantly inhibited. In these cases, one refers to a contracture. In order to predict the evolution of post-wounding skin, several mathematical model frameworks have been set up. These frameworks are based on complicated systems of partial differential equations that need finite element-like discretizations for the approximation of the solution. Since these computational frameworks can be expensive in terms of computation time and resources, we study the applicability of neural networks to reproduce the finite element results. Our neural network is able to simulate the evolution of skin in terms of contraction for over one year. The simulations are based on 25 input parameters that are characteristic for the patient and the injury. One of such input parameters is the stiffness of the skin. The neural network results have yielded an average goodness of fit ( R 2 ) of 0.9928 (± 0.0013). Further, a tremendous speed-up of 19354X was obtained with the neural network. We illustrate the applicability by an online medical App that takes into account the age of the patient and the length of the burn. Supplementary Information: The online version contains supplementary material available at 10.1007/s00521-021-06772-3.

Scar formation from the perspective of complexity science: a new look at the biological system as a whole.

A burn wound is a complex systemic disease at multiple levels. Current knowledge of scar formation after burn injury has come from traditional biological and clinical studies. These are normally focused on just a small part of the entire process, which has limited our ability to sufficiently understand the underlying mechanisms and to predict systems behaviour. Scar formation after burn injury is a result of a complex biological system-wound healing. It is a part of a larger whole. In this self-organising system, many components form networks of interactions with each other. These networks of interactions are typically non-linear and change their states dynamically, responding to the environment and showing emergent long-term behaviour. How molecular and cellular data relate to clinical phenomena, especially regarding effective therapies of burn wounds to achieve minimal scarring, is difficult to unravel and comprehend. Complexity science can help bridge this gap by integrating small parts into a larger whole, such that relevant biological mechanisms and data are combined in a computational model to better understand the complexity of the entire biological system. A better understanding of the complex biological system of post-burn scar formation could bring research and treatment regimens to the next level. The aim of this review/position paper is to create more awareness of complexity in scar formation after burn injury by describing the basic principles of complexity science and its potential for burn care professionals.

Stability of a one-dimensional morphoelastic model for post-burn contraction.

To deal with permanent deformations and residual stresses, we consider a morphoelastic model for the scar formation as the result of wound healing after a skin trauma. Next to the mechanical components such as strain and displacements, the model accounts for biological constituents such as the concentration of signaling molecules, the cellular densities of fibroblasts and myofibroblasts, and the density of collagen. Here we present stability constraints for the one-dimensional counterpart of this morphoelastic model, for both the continuous and (semi-) discrete problem. We show that the truncation error between these eigenvalues associated with the continuous and semi-discrete problem is of order [Formula: see text]. Next we perform numerical validation to these constraints and provide a biological interpretation of the (in)stability. For the mechanical part of the model, the results show the components reach equilibria in a (non) monotonic way, depending on the value of the viscosity. The results show that the parameters of the chemical part of the model need to meet the stability constraint, depending on the decay rate of the signaling molecules, to avoid unrealistic results.

Sensitivity and feasibility of a one-dimensional morphoelastic model for post-burn contraction.

We consider a one-dimensional morphoelastic model describing post-burn scar contraction. Contraction can lead to a limited range of motion (contracture). Reported prevalence of burn scar contractures are 58.6% at 3-6 weeks and 20.9% at 12 months post-reconstructive surgery after burns. This model describes the displacement of the dermal layer of the skin and the development of the effective Eulerian strain in the tissue. Besides these components, the model also contains components that play a major role in the skin repair after trauma. These components are signaling molecules, fibroblasts, myofibroblasts, and collagen. We perform a sensitivity analysis for many parameters of the model and use the results for a feasibility study. In this study, we test whether the model is suitable for predicting the extent of contraction in different age groups. To this end, we conduct an extensive literature review to find parameter values. From the sensitivity analysis, we conclude that the most sensitive parameters are the equilibrium collagen concentration in the dermal layer, the apoptosis rate of fibroblasts and myofibroblasts, and the secretion rate of signaling molecules. Further, although we can use the model to simulate significant distinct contraction densities in different age groups, our results differ from what is seen in the clinic. This particularly concerns children and elderly patients. In children we see more intense contractures if the burn injury occurs near a joint, because the growth induces extra forces on the tissue. Elderly patients seem to suffer less from contractures, possibly because of excess skin.

Agent-based modelling and parameter sensitivity analysis with a finite-element method for skin contraction.

In this paper, we extend the model of wound healing by Boon et al. (J Biomech 49(8):1388-1401, 2016). In addition to explaining the model explicitly regarding every component, namely cells, signalling molecules and tissue bundles, we categorized fibroblasts as regular fibroblasts and myofibroblasts. We do so since it is widely documented that myofibroblasts play a significant role during wound healing and skin contraction and that they are the main phenotype of cells that is responsible for the permanent deformations. Furthermore, we carried out some sensitivity tests of the model by modifying certain parameter values, and we observe that the model shows some consistency with several biological phenomena. Using Monte Carlo simulations, we found that there is a significant strong positive correlation between the final wound area and the minimal wound area. The high correlation between the wound area after 4 days and the final/minimal wound area makes it possible for physicians to predict the most probable time evolution of the wound of the patient. However, the collagen density ratio at the time when the wound area reaches its equilibrium and minimum, cannot indicate the degree of wound contractions, whereas at the 4th day post-wounding, when the collagen is accumulating from null, there is a strong negative correlation between the area and the collagen density ratio. Further, under the circumstances that we modelled, the probability that patients will end up with 5% contraction is about 0.627.

Uncertainty quantification on a spatial Markov-chain model for the progression of skin cancer.

A spatial Markov-chain model is formulated for the progression of skin cancer. The model is based on the division of the computational domain into nodal points, that can be in a binary state: either in 'cancer state' or in 'non-cancer state'. The model assigns probabilities for the non-reversible transition from 'non-cancer' state to the 'cancer state' that depend on the states of the neighbouring nodes. The likelihood of transition further depends on the life burden intensity of the UV-rays that the skin is exposed to. The probabilistic nature of the process and the uncertainty in the input data is assessed by the use of Monte Carlo simulations. A good fit between experiments on mice and our model has been obtained.

Computational modeling of therapy on pancreatic cancer in its early stages.

More than eighty percent of pancreatic cancer involves ductal adenocarcinoma with an abundant desmoplastic extracellular matrix surrounding the solid tumor entity. This aberrant tumor microenvironment facilitates a strong resistance of pancreatic cancer to medication. Although various therapeutic strategies have been reported to be effective in mice with pancreatic cancer, they still need to be tested quantitatively in wider animal-based experiments before being applied as therapies. To aid the design of experiments, we develop a cell-based mathematical model to describe cancer progression under therapy with a specific application to pancreatic cancer. The displacement of cells is simulated by solving a large system of stochastic differential equations with the Euler-Maruyama method. We consider treatment with the PEGylated drug PEGPH20 that breaks down hyaluronan in desmoplastic stroma followed by administration of the chemotherapy drug gemcitabine to inhibit the proliferation of cancer cells. Modeling the effects of PEGPH20 + gemcitabine concentrations is based on Green's fundamental solutions of the reaction-diffusion equation. Moreover, Monte Carlo simulations are performed to quantitatively investigate uncertainties in the input parameters as well as predictions for the likelihood of success of cancer therapy. Our simplified model is able to simulate cancer progression and evaluate treatments to inhibit the progression of cancer.

Can Mathematics and Computational Modeling Help Treat Deep Tissue Injuries?

Objective: Improving the treatment of deep tissue injuries, such as burns, by the use of computational modeling, instead of by animal experiments. Approach: Development of mathematical relations between various parameters and processes. Furthermore, solving the resulting problems through the use of numerical methods, such as finite-element methods. Results: Using our framework, we are able to simulate wound contraction in two dimensions, in which the wound area is followed over time. Our studies indicate that the degree of contraction can be reduced if the appearance of myofibroblasts is inhibited and if their apoptosis is enhanced. Furthermore, after skin grafting, splinting procedures are to be continued as long as TG-beta like growth factor levels are significant. Innovation: A morphoelasticity-based and computational-probabilistic framework for studying the evolution of burn injuries. Conclusion: The current framework is able to reproduce the time evolution of the wound area as observed in clinical results for skin grafts.

A phenomenological model for cell and nucleus deformation during cancer metastasis.

Cell migration plays an essential role in cancer metastasis. In cancer invasion through confined spaces, cells must undergo extensive deformation, which is a capability related to their metastatic potentials. Here, we simulate the deformation of the cell and nucleus during invasion through a dense, physiological microenvironment by developing a phenomenological computational model. In our work, cells are attracted by a generic emitting source (e.g., a chemokine or stiffness signal), which is treated by using Green's Fundamental solutions. We use an IMEX integration method where the linear parts and the nonlinear parts are treated by using an Euler backward scheme and an Euler forward method, respectively. We develop the numerical model for an obstacle-induced deformation in 2D or/and 3D. Considering the uncertainty in cell mobility, stochastic processes are incorporated and uncertainties in the input variables are evaluated using Monte Carlo simulations. This quantitative study aims at estimating the likelihood for invasion and the length of the time interval in which the cell invades the tissue through an obstacle. Subsequently, the two-dimensional cell deformation model is applied to simplified cancer metastasis processes to serve as a model for in vivo or in vitro biomedical experiments.

A model for cell migration in non-isotropic fibrin networks with an application to pancreatic tumor islets.

Cell migration, known as an orchestrated movement of cells, is crucially important for wound healing, tumor growth, immune response as well as other biomedical processes. This paper presents a cell-based model to describe cell migration in non-isotropic fibrin networks around pancreatic tumor islets. This migration is determined by the mechanical strain energy density as well as cytokines-driven chemotaxis. Cell displacement is modeled by solving a large system of ordinary stochastic differential equations where the stochastic parts result from random walk. The stochastic differential equations are solved by the use of the classical Euler-Maruyama method. In this paper, the influence of anisotropic stromal extracellular matrix in pancreatic tumor islets on T-lymphocytes migration in different immune systems is investigated. As a result, tumor peripheral stromal extracellular matrix impedes the immune response of T-lymphocytes through changing direction of their migration.

Biomedical implications from a morphoelastic continuum model for the simulation of contracture formation in skin grafts that cover excised burns.

A continuum hypothesis-based model is developed for the simulation of the (long term) contraction of skin grafts that cover excised burns in order to obtain suggestions regarding the ideal length of splinting therapy and when to start with this therapy such that the therapy is effective optimally. Tissue is modeled as an isotropic, heterogeneous, morphoelastic solid. With respect to the constituents of the tissue, we selected the following constituents as primary model components: fibroblasts, myofibroblasts, collagen molecules, and a generic signaling molecule. Good agreement is demonstrated with respect to the evolution over time of the surface area of unmeshed skin grafts that cover excised burns between outcomes of computer simulations obtained in this study and scar assessment data gathered previously in a clinical study. Based on the simulation results, we suggest that the optimal point in time to start with splinting therapy is directly after placement of the skin graft on its recipient bed. Furthermore, we suggest that it is desirable to continue with splinting therapy until the concentration of the signaling molecules in the grafted area has become negligible such that the formation of contractures can be prevented. We conclude this study with a presentation of some alternative ideas on how to diminish the degree of contracture formation that are not based on a mechanical intervention, and a discussion about how the presented model can be adjusted.

A biomechanical mathematical model for the collagen bundle distribution-dependent contraction and subsequent retraction of healing dermal wounds.

A continuum hypothesis-based, biomechanical model is presented for the simulation of the collagen bundle distribution-dependent contraction and subsequent retraction of healing dermal wounds that cover a large surface area. Since wound contraction mainly takes place in the dermal layer of the skin, solely a portion of this layer is included explicitly into the model. This portion of dermal layer is modeled as a heterogeneous, orthotropic continuous solid with bulk mechanical properties that are locally dependent on both the local concentration and the local geometrical arrangement of the collagen bundles. With respect to the dynamic regulation of the geometrical arrangement of the collagen bundles, it is assumed that a portion of the collagen molecules are deposited and reoriented in the direction of movement of (myo)fibroblasts. The remainder of the newly secreted collagen molecules are deposited by ratio in the direction of the present collagen bundles. Simulation results show that the distribution of the collagen bundles influences the evolution over time of both the shape of the wounded area and the degree of overall contraction of the wounded area. Interestingly, these effects are solely a consequence of alterations in the initial overall distribution of the collagen bundles, and not a consequence of alterations in the evolution over time of the different cell densities and concentrations of the modeled constituents. In accordance with experimental observations, simulation results show furthermore that ultimately the majority of the collagen molecules ends up permanently oriented toward the center of the wound and in the plane that runs parallel to the surface of the skin.

A mathematical model for the simulation of the contraction of burns.

A continuum hypothesis-based model is developed for the simulation of the contraction of burns in order to gain new insights into which elements of the healing response might have a substantial influence on this process. Tissue is modeled as a neo-Hookean solid. Furthermore, (myo)fibroblasts, collagen molecules, and a generic signaling molecule are selected as model components. An overview of the custom-made numerical algorithm is presented. Subsequently, good agreement is demonstrated with respect to variability in the evolution of the surface area of burns over time between the outcomes of computer simulations and measurements obtained in an experimental study. In the model this variability is caused by varying the values for some of its parameters simultaneously. A factorial design combined with a regression analysis are used to quantify the individual contributions of these parameter value variations to the dispersion in the surface area of healing burns. The analysis shows that almost all variability in the surface area can be explained by variability in the value for the myofibroblast apoptosis rate and, to a lesser extent, the value for the collagen molecule secretion rate. This suggests that most of the variability in the evolution of the surface area of burns over time in the experimental study might be attributed to variability in these two rates. Finally, a probabilistic analysis is used in order to investigate in more detail the effect of variability in the values for the two rates on the healing process. Results of this analysis are presented and discussed.

A mathematical model for the simulation of the formation and the subsequent regression of hypertrophic scar tissue after dermal wounding.

A continuum hypothesis-based model is presented for the simulation of the formation and the subsequent regression of hypertrophic scar tissue after dermal wounding. Solely the dermal layer of the skin is modeled explicitly and it is modeled as a heterogeneous, isotropic and compressible neo-Hookean solid. With respect to the constituents of the dermal layer, the following components are selected as primary model components: fibroblasts, myofibroblasts, a generic signaling molecule and collagen molecules. A good match with respect to the evolution of the thickness of the dermal layer of scars between the outcomes of simulations and clinical measurements on hypertrophic scars at different time points after injury in human subjects is demonstrated. Interestingly, the comparison between the outcomes of the simulations and the clinical measurements demonstrates that a relatively high apoptosis rate of myofibroblasts results in scar tissue that behaves more like normal scar tissue with respect to the evolution of the thickness of the tissue over time, while a relatively low apoptosis rate results in scar tissue that behaves like hypertrophic scar tissue with respect to the evolution of the thickness of the tissue over time. Our ultimate goal is to construct models with which the properties of newly generated tissues that form during wound healing can be predicted with a high degree of certainty. The development of the presented model is considered by us as a step toward their construction.

Review on experiment-based two- and three-dimensional models for wound healing.

Traumatic and chronic wounds are a considerable medical challenge that affects many populations and their treatment is a monetary and time-consuming burden in an ageing society to the medical systems. Because wounds are very common and their treatment is so costly, approaches to reveal the responses of a specific wound type to different medical procedures and treatments could accelerate healing and reduce patient suffering. The effects of treatments can be forecast using mathematical modelling that has the predictive power to quantify the effects of induced changes to the wound-healing process. Wound healing involves a diverse and complex combination of biophysical and biomechanical processes. We review a wide variety of contemporary approaches of mathematical modelling of gap closure and wound-healing-related processes, such as angiogenesis. We provide examples of the understanding and insights that may be garnered using those models, and how those relate to experimental evidence. Mathematical modelling-based simulations can provide an important visualization tool that can be used for illustrational purposes for physicians, patients and researchers.