Search

The dynamics of tumor growth and associated events cover multiple time and spatial scales, generally including extracellular, cellular and intracellular modifications. The main goal of this study is to model the biological and physical behavior of tumor evolution in presence of normal healthy tissue, considering a variety of events involved in the process. These include hyper and hypoactivation of signaling pathways during tumor growth, vessels’ growth, intratumoral vascularization and competition of cancer cells with healthy host tissue. The work addresses two distinctive phases in tumor development—the avascular and vascular phases—and in each stage two cases are considered—with and without normal healthy cells. The tumor growth rate increases considerably as closed vessel loops (anastomoses) form around the tumor cells resulting from tumor induced vascularization. When taking into account the host tissue around the tumor, the results show that competition between normal cells and cancer cells leads to the formation of a hypoxic tumor core within a relatively short period of time. Moreover, a dense intratumoral vascular network is formed throughout the entire lesion as a sign of a high malignancy grade, which is consistent with reported experimental data for several types of solid carcinomas. In comparison with other mathematical models of tumor development, in this work we introduce a multiscale simulation that models the cellular interactions and cell behavior as a consequence of the activation of oncogenes and deactivation of gene signaling pathways within each cell. Simulating a therapy that blocks relevant signaling pathways results in the prevention of further tumor growth and leads to an expressive decrease in its size (82% in the simulation)., Author summary Mathematical modeling and simulation of cancer across different biological scales is becoming increasingly important in the development of therapeutic strategies. In the current work, a multiscale model is presented to study the growth and progression of tumor and angiogenesis based on tumor-host interactions which allows investigating the effects of tumor-targeted therapy. Considering the signal-transduction networks involved in various types of cancers, we proposed a cascade that encompasses some significant signaling pathways. A Boolean network model is employed to describe the receptors cross-talk. As a result of the activation of oncogenes and deactivation of pertinent gene signaling pathways within each cell, the cellular interactions and cell behavior are modeled. By linking cells state with environmental cues, the tumor morphology is determined. Consistent with the experimental observations, the intratumoral vascularization density resulting from the simulation reports malignancy grade as a prognostic parameter. Moreover, our model permits to explore possible novel therapeutic procedures, including therapies targeting specific pathways. It captures cellular apoptosis by receptor inhibition in tumor development as a new area of mathematical modeling of targeted therapy.

We examined keratin aggregate formation and the possible mechanisms involved. With this aim, we observed the effect that different ratios between mutant and wild-type keratins expressed in cultured keratinocytes may have on aggregate formation in vitro, as well as how keratin aggregate formation affects the mechanical properties of cells at the cell cortex. To this end we prepared clones with expression rates as close as possible to 25%, 50% and 100% of the EGFP-K14 proteins (either WT or R125P and V270M mutants). Our results showed that only in the case of the 25% EGFP-K14 R125P mutant significant differences could be seen. Namely, we observed in this case the largest accumulation of keratin aggregates and a significant reduction in cell stiffness. To gain insight into the possible mechanisms behind this observation, we extended our previous mathematical model of keratin dynamics by implementing a more complex reaction network that considers the coexistence of wild-type and mutant keratins in the cell. The new model, consisting of a set of coupled, non-linear, ordinary differential equations, allowed us to draw conclusions regarding the relative amounts of intermediate filaments and aggregates in cells, and suggested that aggregate formation by asymmetric binding between wild-type and mutant keratins could explain the data obtained on cells grown in culture.

Keratin intermediate filaments are the principal structural element of epithelial cells. Their importance in providing bulk cellular stiffness is well recognized, but their role in the mechanics of cell cortex is less understood. In this study, we therefore compared the cortical stiffness of three keratinocyte lines: primary wild type cells (NHEK2), immortalized wild type cells (NEB1) and immortalized mutant cells (KEB7). The cortical stiffness was measured by lateral indentation of cells with AOD-steered optical tweezers without employing any moving mechanical elements. The method was validated on fixed cells and Cytochalasin-D treated cells to ensure that the observed variations in stiffness within a single cell line were not a consequence of low measurement precision. The measurements of the cortical stiffness showed that primary wild type cells were significantly stiffer than immortalized wild type cells, which was also detected in previous studies of bulk elasticity. In addition, a small difference between the mutant and the wild type cells was detected, showing that mutation of keratin impacts also the cell cortex. Thus, our results indicate that the role of keratins in cortical stiffness is not negligible and call for further investigation of the mechanical interactions between keratins and elements of the cell cortex.

During angiogenesis, new blood vessels sprout and grow from existing ones. This process plays a crucial role in organ development and repair, in wound healing and in numerous pathological processes such as cancer progression or diabetes. Here, we present a mathematical model of early stage angiogenesis that permits exploration of the relative importance of mechanical, chemical and cellular cues. Endothelial cells proliferate and move over an extracellular matrix by following external gradients of Vessel Endothelial Growth Factor, adhesion and stiffness, which are incorporated to a Cellular Potts model with a finite element description of elasticity. The dynamics of Notch signaling involving Delta-4 and Jagged-1 ligands determines tip cell selection and vessel branching. Through their production rates, competing Jagged-Notch and Delta-Notch dynamics determine the influence of lateral inhibition and lateral induction on the selection of cellular phenotypes, branching of blood vessels, anastomosis (fusion of blood vessels) and angiogenesis velocity. Anastomosis may be favored or impeded depending on the mechanical configuration of strain vectors in the ECM near tip cells. Numerical simulations demonstrate that increasing Jagged production results in pathological vasculatures with thinner and more abundant vessels, which can be compensated by augmenting the production of Delta ligands., Author summary Angiogenesis is the process by which new blood vessels grow from existing ones. This process plays a crucial role in organ development, in wound healing and in numerous pathological processes such as cancer growth or in diabetes. Angiogenesis is a complex, multi-step and well regulated process where biochemistry and physics are intertwined. The process entails signaling in vessel cells being driven by both chemical and mechanical mechanisms that result in vascular cell movement, deformation and proliferation. Mathematical models have the ability to bring together these mechanisms in order to explore their relative relevance in vessel growth. Here, we present a mathematical model of early stage angiogenesis that is able to explore the role of biochemical signaling and tissue mechanics. We use this model to unravel the regulating role of Jagged, Notch and Delta dynamics in vascular cells. These membrane proteins have an important part in determining the leading cell in each neo-vascular sprout. Numerical simulations demonstrate that increasing Jagged production results in pathological vasculatures with thinner and more abundant vessels, which can be compensated by augmenting the production of Delta ligands.

We relate vascular network structure to hemodynamics after vessel obstructions. We consider tree-like networks with a viscoelastic fluid with the rheological characteristics of blood. We analyze the network hemodynamic response, which is a function of the frequencies involved in the driving, and a measurement of the resistance to flow. This response function allows the study of the hemodynamics of the system, without the knowledge of a particular pressure gradient. We find analytical expressions for the network response, which explicitly show the roles played by the network structure, the degree of obstruction, and the geometrical place in which obstructions occur. Notably, we find that the sequence of resistances of the network without occlusions strongly determines the tendencies that the response function has with the anatomical place where obstructions are located. We identify anatomical sites in a network that are critical for its overall capacity to supply blood to a tissue after obstructions. We demonstrate that relatively small obstructions in such critical sites are able to cause a much larger decrease on flow than larger obstructions placed in non-critical sites. Our results indicate that, to a large extent, the response of the network is determined locally. That is, it depends on the structure that the vasculature has around the place where occlusions are found. This result is manifest in a network that follows Murray's law, which is in reasonable agreement with several mammalian vasculatures. For this one, occlusions in early generation vessels have a radically different effect than occlusions in late generation vessels occluding the same percentage of area available to flow. This locality implies that whenever there is a tissue irrigated by a tree-like in vivo vasculature, our model is able to interpret how important obstructions are for the irrigation of such tissue.

Understanding tumor induced angiogenesis is a challenging problem with important consequences for diagnosis and treatment of cancer. Recently, strong evidences suggest the dual role of endothelial cells on the migrating tips and on the proliferating body of blood vessels, in consonance with further events behind lumen formation and vascular patterning. In this paper we present a multi-scale phase-field model that combines the benefits of continuum physics description and the capability of tracking individual cells. The model allows us to discuss the role of the endothelial cells' chemotactic response and proliferation rate as key factors that tailor the neovascular network. Importantly, we also test the predictions of our theoretical model against relevant experimental approaches in mice that displayed distinctive vascular patterns. The model reproduces the in vivo patterns of newly formed vascular networks, providing quantitative and qualitative results for branch density and vessel diameter on the order of the ones measured experimentally in mouse retinas. Our results highlight the ability of mathematical models to suggest relevant hypotheses with respect to the role of different parameters in this process, hence underlining the necessary collaboration between mathematical modeling, in vivo imaging and molecular biology techniques to improve current diagnostic and therapeutic tools.

Sprouting angiogenesis, where new blood vessels grow from pre-existing ones, is a complex process where biochemical and mechanical signals regulate endothelial cell proliferation and movement. Therefore, a mathematical description of sprouting angiogenesis has to take into consideration biological signals as well as relevant physical processes, in particular the mechanical interplay between adjacent endothelial cells and the extracellular microenvironment. In this work, we introduce the first phase-field continuous model of sprouting angiogenesis capable of predicting sprout morphology as a function of the elastic properties of the tissues and the traction forces exerted by the cells. The model is very compact, only consisting of three coupled partial differential equations, and has the clear advantage of a reduced number of parameters. This model allows us to describe sprout growth as a function of the cell-cell adhesion forces and the traction force exerted by the sprout tip cell. In the absence of proliferation, we observe that the sprout either achieves a maximum length or, when the traction and adhesion are very large, it breaks. Endothelial cell proliferation alters significantly sprout morphology, and we explore how different types of endothelial cell proliferation regulation are able to determine the shape of the growing sprout. The largest region in parameter space with well formed long and straight sprouts is obtained always when the proliferation is triggered by endothelial cell strain and its rate grows with angiogenic factor concentration. We conclude that in this scenario the tip cell has the role of creating a tension in the cells that follow its lead. On those first stalk cells, this tension produces strain and/or empty spaces, inevitably triggering cell proliferation. The new cells occupy the space behind the tip, the tension decreases, and the process restarts. Our results highlight the ability of mathematical models to suggest relevant hypotheses with respect to the role of forces in sprouting, hence underlining the necessary collaboration between modelling and molecular biology techniques to improve the current state-of-the-art., Author Summary Sprouting angiogenesis—a process by which new blood vessels grow from existing ones—is an ubiquitous phenomenon in health and disease of higher organisms, playing a crucial role in organogenesis, wound healing, inflammation, as well as on the onset and progression of over 50 different diseases such as cancer, rheumatoid arthritis and diabetes. Mathematical models have the ability to suggest relevant hypotheses with respect to the mechanisms of cell movement and rearrangement within growing vessel sprouts. The inclusion of both biochemical and mechanical processes in a mathematical model of sprouting angiogenesis permits to describe sprout extension as a function of the forces exerted by the cells in the tissue. It also allows to question the regulation of biochemical processes by mechanical forces and vice-versa. In this work we present a compact model of sprouting angiogenesis that includes the mechanical characteristics of the vessel and the tissue. We use this model to suggest the mechanism for the regulation of proliferation within sprout formation. We conclude that the tip cell has the role of creating a tension in the cells that follow its lead. On those first cells of the stalk, this tension produces strain and/or empty spaces, inevitably triggering cell proliferation. The new cells occupy the space behind the tip, the tension decreases, and the process restarts. The modelling strategy used, deemed phase-field, permits to describe the evolution of the shape of different domains in complex systems. It is focused on the movement of the interfaces between the domains, and not on an exhaustive description of the transport properties within each domain. For this reason, it requires a reduced number of parameters, and has been used extensively in modelling other biological phenomena such as tumor growth. The coupling of mechanical and biochemical processes in a compact mathematical model of angiogenesis will enable the study of lumen formation and aneurisms in the near future. Also, this framework will allow the study of the action of flow in vessel remodelling, since local forces can readily be coupled with cell movement to obtain the final vessel morphology.