A cell-based model of extracellular-matrix-guided endothelial cell migration during angiogenesis

Modeling the extracellular matrix in cell migration and morphogenesis: a guide for the curious biologist

The extracellular matrix (ECM) is a highly complex structure through which biochemical and mechanical signals are transmitted. In processes of cell migration, the ECM also acts as a scaffold, providing structural support to cells as well as points of potential attachment. Although the ECM is a well-studied structure, its role in many biological processes remains difficult to investigate comprehensively due to its complexity and structural variation within an organism. In tandem with experiments, mathematical models are helpful in refining and testing hypotheses, generating predictions, and exploring conditions outside the scope of experiments. Such models can be combined and calibrated with in vivo and in vitro data to identify critical cell-ECM interactions that drive developmental and homeostatic processes, or the progression of diseases. In this review, we focus on mathematical and computational models of the ECM in processes such as cell migration including cancer metastasis, and in tissue structure and morphogenesis. By highlighting the predictive power of these models, we aim to help bridge the gap between experimental and computational approaches to studying the ECM and to provide guidance on selecting an appropriate model framework to complement corresponding experimental studies.

Computational modeling of angiogenesis: The importance of cell rearrangements during vascular growth

Angiogenesis is the process wherein endothelial cells (ECs) form sprouts that elongate from the pre-existing vasculature to create new vascular networks. In addition to its essential role in normal development, angiogenesis plays a vital role in pathologies such as cancer, diabetes and atherosclerosis. Mathematical and computational modeling has contributed to unraveling its complexity. Many existing theoretical models of angiogenic sprouting are based on the "snail-trail" hypothesis. This framework assumes that leading ECs positioned at sprout tips migrate toward low-oxygen regions while other ECs in the sprout passively follow the leaders' trails and proliferate to maintain sprout integrity. However, experimental results indicate that, contrary to the snail-trail assumption, ECs exchange positions within developing vessels, and the elongation of sprouts is primarily driven by directed migration of ECs. The functional role of cell rearrangements remains unclear. This review of the theoretical modeling of angiogenesis is the first to focus on the phenomenon of cell mixing during early sprouting. We start by describing the biological processes that occur during early angiogenesis, such as phenotype specification, cell rearrangements and cell interactions with the microenvironment. Next, we provide an overview of various theoretical approaches that have been employed to model angiogenesis, with particular emphasis on recent in silico models that account for the phenomenon of cell mixing. Finally, we discuss when cell mixing should be incorporated into theoretical models and what essential modeling components such models should include in order to investigate its functional role. This article is categorized under: Cardiovascular Diseases > Computational Models Cancer > Computational Models.

Long-term day-by-day tracking of microvascular networks sprouting in fibrin gels: From detailed morphological analyses to general growth rules

Understanding and controlling of the evolution of sprouting vascular networks remains one of the basic challenges in tissue engineering. Previous studies on the vascularization dynamics have typically focused only on the phase of intense growth and often lacked spatial control over the initial cell arrangement. Here, we perform long-term day-by-day analysis of tens of isolated microvasculatures sprouting from endothelial cell-coated spherical beads embedded in an external fibrin gel. We systematically study the topological evolution of the sprouting networks over their whole lifespan, i.e., for at least 14 days. We develop a custom image analysis toolkit and quantify (i) the overall length and area of the sprouts, (ii) the distributions of segment lengths and branching angles, and (iii) the average number of branch generations-a measure of network complexity. We show that higher concentrations of vascular endothelial growth factor (VEGF) lead to earlier sprouting and more branched networks, yet without significantly affecting the speed of growth of individual sprouts. We find that the mean branching angle is weakly dependent on VEGF and typically in the range of 60°-75°, suggesting that, by comparison with the available diffusion-limited growth models, the bifurcating tips tend to follow local VEGF gradients. At high VEGF concentrations, we observe exponential distributions of segment lengths, which signify purely stochastic branching. Our results-due to their high statistical relevance-may serve as a benchmark for predictive models, while our new image analysis toolkit, offering unique features and high speed of operation, could be exploited in future angiogenic drug tests.

Three-component contour dynamics model to simulate and analyze amoeboid cell motility in two dimensions

Amoeboid cell motility is relevant in a wide variety of biomedical processes such as wound healing, cancer metastasis, and embryonic morphogenesis. It is characterized by pronounced changes of the cell shape associated with expansions and retractions of the cell membrane, which result in a crawling kind of locomotion. Despite existing computational models of amoeboid motion, the inference of expansion and retraction components of individual cells, the corresponding classification of cells, and the a priori specification of the parameter regime to achieve a specific motility behavior remain challenging open problems. We propose a novel model of the spatio-temporal evolution of two-dimensional cell contours comprising three biophysiologically motivated components: a stochastic term accounting for membrane protrusions and two deterministic terms accounting for membrane retractions by regularizing the shape and area of the contour. Mathematically, these correspond to the intensity of a self-exciting Poisson point process, the area-preserving curve-shortening flow, and an area adjustment flow. The model is used to generate contour data for a variety of qualitatively different, e.g., polarized and non-polarized, cell tracks that visually resemble experimental data very closely. In application to experimental cell tracks, we inferred the protrusion component and examined its correlation to common biomarkers: the F-actin density close to the membrane and its local motion. Due to the low model complexity, parameter estimation is fast, straightforward, and offers a simple way to classify contour dynamics based on two locomotion types: the amoeboid and a so-called fan-shaped type. For both types, we use cell tracks segmented from fluorescence imaging data of the model organism Dictyostelium discoideum. An implementation of the model is provided within the open-source software package AmoePy, a Python-based toolbox for analyzing and simulating amoeboid cell motility.

Self-assembled and perfusable microvasculature-on-chip for modeling leukocyte trafficking

Leukocyte recruitment from blood to tissue is a process that occurs at the level of capillary vessels during both physiological and pathological conditions. This process is also relevant for evaluating novel adoptive cell therapies, in which the trafficking of therapeutic cells such as chimeric antigen receptor (CAR)-T cells throughout the capillaries of solid tumors is important. Local variations in blood flow, mural cell concentration, and tissue stiffness contribute to the regulation of capillary vascular permeability and leukocyte trafficking throughout the capillary microvasculature. We developed a platform to mimic a biologically functional human arteriole-venule microcirculation system consisting of pericytes (PCs) and arterial and venous primary endothelial cells (ECs) embedded within a hydrogel, which self-assembles into a perfusable, heterogeneous microvasculature. Our device shows a preferential association of PCs with arterial ECs that drives the flow-dependent formation of microvasculature networks. We show that PCs stimulate basement membrane matrix synthesis, which affects both vessel diameter and permeability in a manner correlating with the ratio of ECs to PCs. Moreover, we demonstrate that hydrogel concentration can affect capillary morphology but has no observed effect on vascular permeability. The biological function of our capillary network was demonstrated using an inflammation model, where significantly higher expression of cytokines, chemokines, and adhesion molecules was observed after tumor necrosis factor-alpha (TNF-α) treatment. Accordingly, T cell adherence and transendothelial migration were significantly increased in the immune-activated state. Taken together, our platform allows the generation of a perfusable microvasculature that recapitulates the structure and function of an in vivo capillary bed that can be used as a model for developing potential immunotherapies.

A three-dimensional model with two-body interactions for endothelial cells in angiogenesis

We introduce a three-dimensional mathematical model for the dynamics of vascular endothelial cells during sprouting angiogenesis. Angiogenesis is the biological process by which new blood vessels form from existing ones. It has been the subject of numerous theoretical models. These models have successfully replicated various aspects of angiogenesis. Recent studies using particle-based models have highlighted the significant influence of cell shape on network formation, with elongated cells contributing to the formation of branching structures. While most mathematical models are two-dimensional, we aim to investigate whether ellipsoids also form branch-like structures and how their shape affects the pattern. In our model, the shape of a vascular endothelial cell is represented as a spheroid, and a discrete dynamical system is constructed based on the simple assumption of two-body interactions. Numerical simulations demonstrate that our model reproduces the patterns of elongation and branching observed in the early stages of angiogenesis. We show that the pattern formation of the cell population is strongly dependent on the cell shape. Finally, we demonstrate that our current mathematical model reproduces the cell behaviours, specifically cell-mixing, observed in sprouts.

Approaches to vascular network, blood flow, and metabolite distribution modeling in brain tissue

The cardiovascular system plays a key role in the transport of nutrients, ensuring a continuous supply of all cells of the body with the metabolites necessary for life. The blood supply to the brain is carried out by the large arteries located on its surface, which branch into smaller arterioles that penetrate the cerebral cortex and feed the capillary bed, thereby forming an extensive branching network. The formation of blood vessels is carried out via vasculogenesis and angiogenesis, which play an important role in both embryo and adult life. The review presents approaches to modeling various aspects of both the formation of vascular networks and the construction of the formed arterial tree. In addition, a brief description of models that allows one to study the blood flow in various parts of the circulatory system and the spatiotemporal metabolite distribution in brain tissues is given. Experimental study of these issues is not always possible due to both the complexity of the cardiovascular system and the mechanisms through which the perfusion of all body cells is carried out. In this regard, mathematical models are a good tool for studying hemodynamics and can be used in clinical practice to diagnose vascular diseases and assess the need for treatment.

Hybrid cellular Potts and bead-spring modeling of cells in fibrous extracellular matrix

The mechanical interaction between cells and the extracellular matrix (ECM) is fundamental to coordinate collective cell behavior in tissues. Relating individual cell-level mechanics to tissue-scale collective behavior is a challenge that cell-based models such as the cellular Potts model (CPM) are well-positioned to address. These models generally represent the ECM with mean-field approaches, which assume substrate homogeneity. This assumption breaks down with fibrous ECM, which has nontrivial structure and mechanics. Here, we extend the CPM with a bead-spring model of ECM fiber networks modeled using molecular dynamics. We model a contractile cell pulling with discrete focal adhesion-like sites on the fiber network and demonstrate agreement with experimental spatiotemporal fiber densification and displacement. We show that at high network cross-linking, contractile cell forces propagate over at least eight cell diameters, decaying with distance with power law exponent n= 0.35 - 0.65 typical of viscoelastic ECMs. Further, we use in silico atomic force microscopy to measure local cell-induced network stiffening consistent with experiments. Our model lays the foundation for investigating how local and long-ranged cell-ECM mechanobiology contributes to multicellular morphogenesis.

Semi-autonomous wound invasion via matrix-deposited, haptotactic cues

Proper wound healing relies on invasion of fibroblasts via directed migration. While the related experimental and mathematical modeling literature has mainly focused on cell migration directed by soluble cues (chemotaxis), there is ample evidence that fibroblast migration is also directed by insoluble, matrix-bound cues (haptotaxis). Furthermore, numerous studies indicate that fibronectin (FN), a haptotactic ligand for fibroblasts, is present and dynamic in the provisional matrix throughout the proliferative phase of wound healing. In the present work, we show the plausibility of a hypothesis that fibroblasts themselves form and maintain haptotactic gradients in a semi-autonomous fashion. As a precursor to this, we examine the positive control scenario where FN is pre-deposited in the wound matrix, and fibroblasts maintain haptotaxis by removing FN at an appropriate rate. After developing conceptual and quantitative understanding of this scenario, we consider two cases in which fibroblasts activate the latent form of a matrix-loaded cytokine, TGFβ, which upregulates the fibroblasts' own secretion of FN. In the first of these, the latent cytokine is pre-patterned and released by the fibroblasts. In the second, fibroblasts in the wound produce the latent TGFβ, with the presence of the wound providing the only instruction. In all cases, wound invasion is more effective than a negative control model with haptotaxis disabled; however, there is a trade-off between the degree of fibroblast autonomy and the rate of invasion.

The Matrix Reloaded—The Role of the Extracellular Matrix in Cancer

As the core component of all organs, the extracellular matrix (ECM) is an interlocking macromolecular meshwork of proteins, glycoproteins, and proteoglycans that provides mechanical support to cells and tissues. In cancer, the ECM can be remodelled in response to environmental cues, and it controls a plethora of cellular functions, including metabolism, cell polarity, migration, and proliferation, to sustain and support oncogenesis. The biophysical and biochemical properties of the ECM, such as its structural arrangement and being a reservoir for bioactive molecules, control several intra- and intercellular signalling pathways and induce cytoskeletal changes that alter cell shapes, behaviour, and viability. Desmoplasia is a major component of solid tumours. The abnormal deposition and composition of the tumour matrix lead to biochemical and biomechanical alterations that determine disease development and resistance to treatment. This review summarises the complex roles of ECM in cancer and highlights the possible therapeutic targets and how to potentially remodel the dysregulated ECM in the future. Furthering our understanding of the ECM in cancer is important as the modification of the ECM will probably become an important tool in the characterisation of individual tumours and personalised treatment options.

Data-driven spatio-temporal modelling of glioblastoma

Mathematical oncology provides unique and invaluable insights into tumour growth on both the microscopic and macroscopic levels. This review presents state-of-the-art modelling techniques and focuses on their role in understanding glioblastoma, a malignant form of brain cancer. For each approach, we summarize the scope, drawbacks and assets. We highlight the potential clinical applications of each modelling technique and discuss the connections between the mathematical models and the molecular and imaging data used to inform them. By doing so, we aim to prime cancer researchers with current and emerging computational tools for understanding tumour progression. By providing an in-depth picture of the different modelling techniques, we also aim to assist researchers who seek to build and develop their own models and the associated inference frameworks. Our article thus strikes a unique balance. On the one hand, we provide a comprehensive overview of the available modelling techniques and their applications, including key mathematical expressions. On the other hand, the content is accessible to mathematicians and biomedical scientists alike to accommodate the interdisciplinary nature of cancer research.

Pattern formation of elliptic particles by two-body interactions: A model for dynamics of endothelial cells in angiogenesis

A two-dimensional mathematical model for dynamics of endothelial cells in angiogenesis is investigated. Angiogenesis is a morphogenic process in which new blood vessels emerge from an existing vascular network. Recently a one-dimensional discrete dynamical model has been proposed to reproduce elongation, bifurcation, and cell motility such as cell-mixing during angiogenesis on the assumption of a simple two-body interaction between endothelial cells. The present model is its two-dimensional extension, where endothelial cells are represented as the ellipses with the two-body interactions: repulsive interaction due to excluded volume effect, attractive interaction through pseudopodia and rotation by contact. We show that the oblateness of ellipses and the magnitude of contact rotation significantly affect the shape of created vascular patterns and elongation of branches.

Through history to growth dynamics: deciphering the evolution of spatial networks

Many ramified, network-like patterns in nature, such as river networks or blood vessels, form as a result of unstable growth of moving boundaries in an external diffusive field. Here, we pose the inverse problem for the network growth-can the growth dynamics be inferred from the analysis of the final pattern? We show that by evolving the network backward in time one can not only reconstruct the growth rules but also get an insight into the conditions under which branch splitting occurs. Determining the growth rules from a single snapshot in time is particularly important for growth processes so slow that they cannot be directly observed, such as growth of river networks and deltas or cave passages. We apply this approach to analyze the growth of a real river network in Vermont, USA. We determine its growth rule and argue that branch splitting events are triggered by an increase in the tip growth velocity.

On the role of mechanical signals on sprouting angiogenesis through computer modeling approaches

Sprouting angiogenesis, the formation of new vessels from preexisting vasculature, is an essential process in the regeneration of new tissues as well as in the development of some diseases like cancer. Although early studies identified chemical signaling as the main driver of this process, many recent studies have shown a strong role of mechanical signals in the formation of new capillaries. Different types of mechanical signals (e.g., external forces, cell traction forces, and blood flow-induced shear forces) have been shown to play distinct roles in the process; however, their interplay remains still largely unknown. During the last decades, mathematical and computational modeling approaches have been developed to investigate and better understand the mechanisms behind mechanically driven angiogenesis. In this manuscript, we review computational models of angiogenesis with a focus on models investigating the role of mechanics on the process. Our aim is not to provide a detailed review on model methodology but to describe what we have learnt from these models. We classify models according to the mechanical signals being investigated and describe how models have looked into their role on the angiogenic process. We show that a better understanding of the mechanobiology of the angiogenic process will require the development of computer models that incorporate the interactions between the multiple mechanical signals and their effect on cellular responses, since they all seem to play a key in sprout patterning. In the end, we describe some of the remaining challenges of computational modeling of angiogenesis and discuss potential avenues for future research.

Detachable Microneedle Patches Deliver Mesenchymal Stromal Cell Factor-Loaded Nanoparticles for Cardiac Repair

Intramyocardial injection is a direct and efficient approach to deliver therapeutics to the heart. However, the injected volume must be very limited, and there is injury to the injection site and leakage issues during heart beating. Herein, we developed a detachable therapeutic microneedle (MN) patch, which is comprised of mesenchymal stromal cell-secreted factors (MSCF)-loaded poly(lactic-co-glycolic acid) nanoparticles (NP) in MN tips made of elastin-like polypeptide gel, with a resolvable non-cross-linked hyaluronic acid (HA) gel as the MN base. The tips can be firmly inserted into the infarcted myocardium after base removal, and no suture is needed. In isolated neonatal rat cardiac cells, we found that the cellular uptake of MSCF-NP in the cardiomyocytes was higher than in cardiac fibroblasts. MSCF-NP promoted the proliferation of injured cardiomyocytes. In a rat model of myocardial infarction, MN-MSCF-NP treatment reduced cardiomyocyte apoptosis, restored myocardium volume, and reduced fibrosis during the cardiac remodeling process. Our work demonstrated the therapeutic potential of MN to deliver MSCF directly into the myocardium and provides a promising treatment approach for cardiac diseases.

A mechanistic protrusive-based model for 3D cell migration

Cell migration is essential for a variety of biological processes, such as embryogenesis, wound healing, and the immune response. After more than a century of research-mainly on flat surfaces-, there are still many unknowns about cell motility. In particular, regarding how cells migrate within 3D matrices, which more accurately replicate in vivo conditions. We present a novel in silico model of 3D mesenchymal cell migration regulated by the chemical and mechanical profile of the surrounding environment. This in silico model considers cell's adhesive and nuclear phenotypes, the effects of the steric hindrance of the matrix, and cells ability to degradate the ECM. These factors are crucial when investigating the increasing difficulty that migrating cells find to squeeze their nuclei through dense matrices, which may act as physical barriers. Our results agree with previous in vitro observations where fibroblasts cultured in collagen-based hydrogels did not durotax toward regions with higher collagen concentrations. Instead, they exhibited an adurotactic behavior, following a more random trajectory. Overall, cell's migratory response in 3D domains depends on its phenotype, and the properties of the surrounding environment, that is, 3D cell motion is strongly dependent on the context.

A novel nonlocal partial differential equation model of endothelial progenitor cell cluster formation during the early stages of vasculogenesis

The formation of new vascular networks is essential for tissue development and regeneration, in addition to playing a key role in pathological settings such as ischemia and tumour development. Experimental findings in the past two decades have led to the identification of a new mechanism of neovascularisation, known as cluster-based vasculogenesis, during which endothelial progenitor cells (EPCs) mobilised from the bone marrow are capable of bridging distant vascular beds in a variety of hypoxic settings in vivo. This process is characterised by the formation of EPC clusters during its early stages and, while much progress has been made in identifying various mechanisms underlying cluster formation, we are still far from a comprehensive description of such spatio-temporal dynamics. In order to achieve this, we propose a novel mathematical model of the early stages of cluster-based vasculogenesis, comprising of a system of nonlocal partial differential equations including key mechanisms such as endogenous chemotaxis, matrix degradation, cell proliferation and cell-to-cell adhesion. We conduct a linear stability analysis on the system and solve the equations numerically. We then conduct a parametric analysis of the numerical solutions of the one-dimensional problem to investigate the role of underlying dynamics on the speed of cluster formation and the size of clusters, measured via appropriate metrics for the cluster width and compactness. We verify the key results of the parametric analysis with simulations of the two-dimensional problem. Our results, which qualitatively compare with data from in vitro experiments, elucidate the complementary role played by endogenous chemotaxis and matrix degradation in the formation of clusters, suggesting chemotaxis is responsible for the cluster topology while matrix degradation is responsible for the speed of cluster formation. Our results also indicate that the nonlocal cell-to-cell adhesion term in our model, even though it initially causes cells to aggregate, is not sufficient to ensure clusters are stable over long time periods. Consequently, new modelling strategies for cell-to-cell adhesion are required to stabilise in silico clusters. We end the paper with a thorough discussion of promising, fruitful future modelling and experimental research perspectives.

Capillary network formation and structure in a modified discrete mathematical model of angiogenesis

Angiogenesis, as part of cancer development, involves hierarchical complicated events and processes. Multiple studies have revealed the significance of the formation and structure of tumor-induced capillary networks. In this study, a discrete mathematical model of angiogenesis is studied and modified to capture the realistic physics of capillary network formation. Modifications are performed on the mathematical foundations of an existing discrete model of angiogenesis. The main modifications are the imposition of the matrix density effect, implementation of realistic boundary and initial conditions, and improvement of the method of governing equations based on physical observation. Results show that endothelial cells accelerate angiogenesis and capillary formation as they migrate toward the tumor and clearly exhibit the physical concept of haptotactic movement. On the other hand, consideration of blood flow-induced stress leads to a dynamic adaptive vascular network of capillaries which intelligibly reflects the brush border effect . The present modified model of capillary network formation is based on the physical rationale that defines a clear mathematical and physical interpretation of angiogenesis, which is likely to be used in cancer development modeling and anti-angiogenic therapies.

Mathematical simulation and prediction of tumor volume using RBF artificial neural network at different circumstances in the tumor microenvironment

Uncontrolled proliferation of cells in a tissue caused by genetic mutations inside a cell is referred to as a tumor. A tumor which grows rapidly encounters a barrier when it grows to a certain size in presence of preexisting vasculature. This is the time when it has to find a way to go on the growth. The tumor starts to secrete tumor angiogenic factors (TAFs) and stimulate preexisting vessels to grow new sprouts. These new sprouts will find their way to the tumor in the extracellular matrix (ECM) by the gradient of TAF. As these new capillaries anastomose and reach tumor, fresh oxygen is available for the tumor and it will reinitiate the growth. Number of initial sprouts, distance of initial tumor cells from the vessel(s) and initial density of the tumor at the time of sprout formation are questions which are to be investigated. In the present study, the aim is to find the response of tumor cells and vessels to the reciprocal effects of each other in different circumstances in the tissue. Together with a mathematical formulation, a radial basis function (RBF) neural network is established to predict the number of tumor cells at different circumstances including size and distance of initial tumors from the parent vessel. A final formulation is given for the final number of tumor cells as a function of initial tumor size and distance between a parent vessel and a tumor. Results of this simulation demonstrate that, increasing the distance between a tumor and a parent vessel decreases the number of final tumor cells. Specially, this decrement becomes faster beyond a certain distance. Moreover, initial tumors in bigger domains must become much bigger before inducing angiogenesis which makes it harder for them to survive.

A hybrid model of intercellular tension and cell-matrix mechanical interactions in a multicellular geometry

Epithelial cells form continuous sheets of cells that exist in tensional homeostasis. Homeostasis is maintained through cell-to-cell junctions that distribute tension and balance forces between cells and their underlying matrix. Disruption of tensional homeostasis can lead to epithelial-mesenchymal transition (EMT), a transdifferentiation process in which epithelial cells adopt a mesenchymal phenotype, losing cell-cell adhesion and enhancing cellular motility. This process is critical during embryogenesis and wound healing, but is also dysregulated in many disease states. To further understand the role of intercellular tension in spatial patterning of epithelial cell monolayers, we developed a multicellular computational model of cell-cell and cell-substrate forces. This work builds on a hybrid cellular Potts model (CPM)-finite element model to evaluate cell-matrix mechanical feedback of an adherent multicellular cluster. Cellular movement is governed by thermodynamic constraints from cell volume, cell-cell and cell-matrix contacts, and durotaxis, which arises from cell-generated traction forces on a finite element substrate. Junction forces at cell-cell contacts balance these traction forces, thereby producing a mechanically stable epithelial monolayer. Simulations were compared to in vitro experiments using fluorescence-based junction force sensors in clusters of cells undergoing EMT. Results indicate that the multicellular CPM model can reproduce many aspects of EMT, including epithelial monolayer formation dynamics, changes in cell geometry, and spatial patterning of cell-cell forces in an epithelial tissue.

A multiscale hybrid mathematical model of epidermal-dermal interactions during skin wound healing

Following injury, skin activates a complex wound healing programme. While cellular and signalling mechanisms of wound repair have been extensively studied, the principles of epidermal-dermal interactions and their effects on wound healing outcomes are only partially understood. To gain new insight into the effects of epidermal-dermal interactions, we developed a multiscale, hybrid mathematical model of skin wound healing. The model takes into consideration interactions between epidermis and dermis across the basement membrane via diffusible signals, defined as activator and inhibitor. Simulations revealed that epidermal-dermal interactions are critical for proper extracellular matrix deposition in the dermis, suggesting these signals may influence how wound scars form. Our model makes several theoretical predictions. First, basal levels of epidermal activator and inhibitor help to maintain dermis in a steady state, whereas their absence results in a raised, scar-like dermal phenotype. Second, wound-triggered increase in activator and inhibitor production by basal epidermal cells, coupled with fast re-epithelialization kinetics, reduces dermal scar size. Third, high-density fibrin clot leads to a raised, hypertrophic scar phenotype, whereas low-density fibrin clot leads to a hypotrophic phenotype. Fourth, shallow wounds, compared to deep wounds, result in overall reduced scarring. Taken together, our model predicts the important role of signalling across dermal-epidermal interface and the effect of fibrin clot density and wound geometry on scar formation. This hybrid modelling approach may be also applicable to other complex tissue systems, enabling the simulation of dynamic processes, otherwise computationally prohibitive with fully discrete models due to a large number of variables.

All Roads Lead to Directional Cell Migration

Directional cell migration normally relies on a variety of external signals, such as chemical, mechanical, or electrical, which instruct cells in which direction to move. Many of the major molecular and physical effects derived from these cues are now understood, leading to questions about whether directional cell migration is alike or distinct under these different signals, and how cells might be directed by multiple simultaneous cues, which would be expected in complex in vivo environments. In this review, we compare how different stimuli are spatially distributed, often as gradients, to direct cell movement and the mechanisms by which they steer cells. A comparison of the downstream effectors of directional cues suggests that different external signals regulate a common set of components: small GTPases and the actin cytoskeleton, which implies that the mechanisms downstream of different signals are likely to be closely related and underlies the idea that cell migration operates by a common set of physical principles, irrespective of the input.

Simulating flow induced migration in vascular remodelling

Shear stress induces directed endothelial cell (EC) migration in blood vessels leading to vessel diameter increase and induction of vascular maturation. Other factors, such as EC elongation and interaction between ECs and non-vascular areas are also important. Computational models have previously been used to study collective cell migration. These models can be used to predict EC migration and its effect on vascular remodelling during embryogenesis. We combined live time-lapse imaging of the remodelling vasculature of the quail embryo yolk sac with flow quantification using a combination of micro-Particle Image Velocimetry and computational fluid dynamics. We then used the flow and remodelling data to inform a model of EC migration during remodelling. To obtain the relation between shear stress and velocity in vitro for EC cells, we developed a flow chamber to assess how confluent sheets of ECs migrate in response to shear stress. Using these data as an input, we developed a multiphase, self-propelled particles (SPP) model where individual agents are driven to migrate based on the level of shear stress while maintaining appropriate spatial relationship to nearby agents. These agents elongate, interact with each other, and with avascular agents at each time-step of the model. We compared predicted vascular shape to real vascular shape after 4 hours from our time-lapse movies and performed sensitivity analysis on the various model parameters. Our model shows that shear stress has the largest effect on the remodelling process. Importantly, however, elongation played an especially important part in remodelling. This model provides a powerful tool to study the input of different biological processes on remodelling.

Shear stress is known to play a leading role in endothelial cell (EC) migration and hence, vascular remodelling. Vascular remodelling is, however, more complicated than only EC migration. To achieve a better understanding of this process, we developed a computational model in which, shear stress mediated EC migration has the leading role and other factors, such as avascular regions and EC elongation, are also accounted for. We have tested this model for different vessel shapes during remodelling and could study the role that each of these factors play in remodelling. This model gives us the possibility of addition of other factors such as biochemical signals and angiogenesis which will help us in the study of vascular remodelling in both development and disease.

Umbilical Cord-Derived Mesenchymal Stem Cell-Derived Exosomes Combined Pluronic F127 Hydrogel Promote Chronic Diabetic Wound Healing and Complete Skin Regeneration

Purpose

Chronic refractory wounds are a multifactorial comorbidity of diabetes mellitus with the characteristic of impaired vascular networks. Currently, there is a lack of effective treatments for such wounds. Various types of mesenchymal stem cell-derived exosomes (MSC-exos) have been shown to exert multiple therapeutic effects on skin regeneration. We aimed to determine whether a constructed combination of human umbilical cord MSC (hUCMSC)-derived exosomes (hUCMSC-exos) and Pluronic F-127 (PF-127) hydrogel could improve wound healing.

Materials and Methods

We topically applied human umbilical cord-derived MSC (hUCMSC)-derived exosomes (hUCMSC-exos) encapsulated in a thermosensitive PF-127 hydrogel to a full-thickness cutaneous wound in a streptozotocin-induced diabetic rat model. The material properties and wound healing ability of the hydrogel and cellular responses were analyzed.

Results

Compared with hUCMSC-exos, PF-127-only or control treatment, the combination of PF-127 and hUCMSC-exos resulted in a significantly accelerated wound closure rate, increased expression of CD31 and Ki67, enhanced regeneration of granulation tissue and upregulated expression of vascular endothelial growth factor (VEGF) and factor transforming growth factor beta-1 (TGFβ-1).

Conclusion

The efficient delivery of hUCMSC-exos in PF-127 gel and improved exosome ability could promote diabetic wound healing. Thus, this biomaterial-based exosome therapy may represent a new therapeutic approach for cutaneous regeneration of chronic wounds.

Transcription box‑3 protects human umbilical vein endothelial cells in a high‑glucose environment through sirtuin 1/AKT signaling

The increasing burden of diabetes in low and middle-income countries is attributable to both genetic and epigenetic factors. Environmental- and lifestyle-associated changes are also considered to be important contributors to this disease. The resultant co-morbidities arising from micro-and macrovascular changes in diabetes are difficult to manage and are an economic burden. However, very little is known about the molecular mechanisms that drive this phenotype. The present study aimed to investigate the role of sirtuin 1 (SIRT1)- and transcription box-3 (TBX-3)-mediated regulation of endothelial dysfunction, given the significance of SIRT1 in glucose metabolism and the role of TBX-3 in the maintenance of cellular proliferation, senescence and apoptosis. Following the recruitment of adult patients with and without diabetes, both SIRT1 and TBX-3 expression was confirmed to be present in the sera of the patients with diabetes and the patients without diabetes; however, both SIRT1 and TBX-3 expression levels were higher in the sera of the patients with diabetes. Human umbilical vein endothelial cells (HUVECs) were further used for in vitro studies. Using TBX-3 and SIRT1 knockdown models, the cellular responses to proliferation, migration, invasion and tube formation were investigated using an MTS, cell cycle analysis, wound healing, Transwell and tube formation assay, respectively. Western blotting was also used to determine the downstream signaling pathways involved. The genetic knockdown of TBX-3 in hyperglycemic conditions significantly decreased the cellular proliferation, migration, invasion and angiogenesis of HUVECs. It was subsequently identified that TBX-3 mediated its effects through the activation of AKT and vascular endothelial growth factor (VEGF) signaling. However, the genetic knockdown of SIRT1 in the presence of TBX-3 overexpression and glucose failed to activate the AKT and VEGF signaling pathways. In conclusion, the results of the present study suggested that SIRT1 may positively regulate TBX-3 in endothelial cells, therefore, SIRT1 and/or TBX-3 may serve as potential novel biomarkers for disease progression.

Angiogenesis is Differentially Modulated by Anthocyanin and Phenolic Acid Extracts from Wild Blueberry (V. angustifolium) Through PI3K Pathway

Wild blueberries rank very high in anthocyanins (ACNs) and phenolic acids (PA) among other berries. Previous work from our group has documented their beneficial effects. In this study, human umbilical vein endothelial cells were used, and ACNs, PAs, and their combination (ACNs:PAs) at concentrations of 0.002, 8, 15, 60, and 120 μg/mL were tested for endothelial tube formation. Treatment with ACNs decreased, while treatment with PAs and ACNs:PAs increased overall endothelial cell tube formation compared to control. Endothelial cells exposed to ACNs downregulated gene expression of AKT1 and endothelial nitric oxide synthase (eNOS), while PAs upregulated AKT1 and vascular endothelial growth factor (VEGF) gene expression. Combination of ACNs:PAs decreased gene expression of AKT1 and eNOS, while protein levels of AKT1 increased. In summary, based on the type of wild blueberry extract, angiogenesis is differentially modulated and is concentration dependent. Further experiments will delineate the mechanism(s) of the differential action of the aforementioned extracts on angiogenesis.

Functional evaluation of prevascularization in one-stage versus two-stage tissue engineering approach of human bio-artificial muscle

A common shortcoming of current tissue engineered constructs is the lack of a functional vasculature, limiting their size and functionality. Prevascularization is a possible strategy to introduce vascular networks in these constructs. It includes among others co-culturing target cells with endothelial (precursor) cells that are able to form endothelial networks through vasculogenesis. In this paper, we compared two different prevascularization approaches of bio-artificial skeletal muscle tissue (BAM) in vitro and in vivo. In a one-stage approach, human muscle cells were directly co-cultured with endothelial cells in 3D. In a two-stage approach, a one week old BAM containing differentiated myotubes was coated with a fibrin hydrogel containing endothelial cells. The obtained endothelial networks were longer and better interconnected with the two-stage approach. We evaluated whether prevascularization had a beneficial effect on in vivo perfusion of the BAM and improved myotube survival by implantation on the fascia of the latissimus dorsi muscle of NOD/SCID mice for 5 or 14 d. Also in vivo, the two-stage approach displayed the highest vascular density. At day 14, anastomosis of implanted endothelial networks with the host vasculature was apparent. BAMs without endothelial networks contained longer and thicker myotubes in vitro, but their morphology degraded in vivo. In contrast, maintenance of myotube morphology was well supported in the two-stage prevascularized BAMs. To conclude, a two-stage prevascularization approach for muscle engineering improved the vascular density in the construct and supported myotube maintenance in vivo.

Cell-Scale Degradation of Peritumoural Extracellular Matrix Fibre Network and Its Role Within Tissue-Scale Cancer Invasion

Local cancer invasion of tissue is a complex, multiscale process which plays an essential role in tumour progression. During the complex interaction between cancer cell population and the extracellular matrix (ECM), of key importance is the role played by both bulk two-scale dynamics of ECM fibres within collective movement of the tumour cells and the multiscale leading edge dynamics driven by proteolytic activity of the matrix-degrading enzymes (MDEs) that are secreted by the cancer cells. As these two multiscale subsystems share and contribute to the same tumour macro-dynamics, in this work we develop further the model introduced in Shuttleworth and Trucu (Bull Math Biol 81:2176–2219, 2019. 10.1007/s11538-019-00598-w) by exploring a new aspect of their interaction that occurs at the cell scale. Specifically, here we will focus on understanding the cell-scale cross talk between the micro-scale parts of these two multiscale subsystems which get to interact directly in the peritumoural region, with immediate consequences both for MDE micro-dynamics occurring at the leading edge of the tumour and for the cell-scale rearrangement of the naturally oriented ECM fibres in the peritumoural region, ultimately influencing the way tumour progresses in the surrounding tissue. To that end, we will propose a new modelling that captures the ECM fibres degradation not only at macro-scale in the bulk of the tumour but also explicitly in the micro-scale neighbourhood of the tumour interface as a consequence of the interactions with molecular fluxes of MDEs that exercise their spatial dynamics at the invasive edge of the tumour.

Quasi-periodic migration of single cells on short microlanes

Cell migration on microlanes represents a suitable and simple platform for the exploration of the molecular mechanisms underlying cell cytoskeleton dynamics. Here, we report on the quasi-periodic movement of cells confined in stripe-shaped microlanes. We observe persistent polarized cell shapes and directed pole-to-pole motion within the microlanes. Cells depolarize at one end of a given microlane, followed by delayed repolarization towards the opposite end. We analyze cell motility via the spatial velocity distribution, the velocity frequency spectrum and the reversal time as a measure for depolarization and spontaneous repolarization of cells at the microlane ends. The frequent encounters of a boundary in the stripe geometry provides a robust framework for quantitative investigations of the cytoskeleton protrusion and repolarization dynamics. In a first advance to rigorously test physical models of cell migration, we find that the statistics of the cell migration is recapitulated by a Cellular Potts model with a minimal description of cytoskeleton dynamics. Using LifeAct-GFP transfected cells and microlanes with differently shaped ends, we show that the local deformation of the leading cell edge in response to the tip geometry can locally either amplify or quench actin polymerization, while leaving the average reversal times unaffected.

Fluidity of Poly (ε-Caprolactone)-Based Material Induces Epithelial-to-Mesenchymal Transition

BACKGROUND

We propose the potential studies on material fluidity to induce epithelial to mesenchymal transition (EMT) in MCF-7 cells. In this study, we examined for the first time the effect of material fluidity on EMT using poly(ε-caprolactone-co-D,L-lactide) (P(CL-co-DLLA)) with tunable elasticity and fluidity.

METHODS

The fluidity was altered by chemically crosslinking the polymer networks. The crosslinked P(CL-co-DLLA) substrate showed a solid-like property with a stiffness of 261 kPa, while the non-crosslinked P(CL-co-DLLA) substrate of 100 units (high fluidity) and 500 units (low fluidity) existed in a quasi-liquid state with loss modulus of 33 kPa and 30.8 kPa, respectively, and storage modulus of 10.8 kPa and 20.1 kPa, respectively.

RESULTS

We observed that MCF-7 cells on low fluidic substrates decreased the expression of E-cadherin, an epithelial marker, and increased expression of vimentin, a mesenchymal marker. This showed that the cells lose their epithelial phenotype and gain a mesenchymal property. On the other hand, MCF-7 cells on high fluidic substrates maintained their epithelial phenotype, suggesting that the cells did not undergo EMT.

CONCLUSION

Considering these results as the fundamental information for material fluidity induced EMT, our system could be used to regulate the degree of EMT by turning the fluidity of the material.

Roles of Remote and Contact Forces in Epithelial Cell Structure Formation

Cancer cells collectively form a large-scale structure for their growth. In this article, we report that HeLa cells, epithelial-like human cervical cancer cells, aggressively migrate on Matrigel and form a large-scale structure in a cell-density-dependent manner. To explain the experimental results, we develop a simple model in which cells interact and migrate using the two fundamentally different types of force, remote and contact forces, and show how cells form a large-scale structure. We demonstrate that the simple model reproduces experimental observations, suggesting that the remote and contact forces considered in this work play a major role in large-scale structure formation of HeLa cells. This article provides important evidence that cancer cells form a large-scale structure and develops an understanding into the poorly understood mechanisms of their structure formation.

Mechanochemical Coupling and Junctional Forces during Collective Cell Migration

Cell migration, a fundamental physiological process in which cells sense and move through their surrounding physical environment, plays a critical role in development and tissue formation, as well as pathological processes, such as cancer metastasis and wound healing. During cell migration, dynamics are governed by the bidirectional interplay between cell-generated mechanical forces and the activity of Rho GTPases, a family of small GTP-binding proteins that regulate actin cytoskeleton assembly and cellular contractility. These interactions are inherently more complex during the collective migration of mechanically coupled cells because of the additional regulation of cell-cell junctional forces. In this study, we adapted a recent minimal modeling framework to simulate the interactions between mechanochemical signaling in individual cells and interactions with cell-cell junctional forces during collective cell migration. We find that migration of individual cells depends on the feedback between mechanical tension and Rho GTPase activity in a biphasic manner. During collective cell migration, waves of Rho GTPase activity mediate mechanical contraction/extension and thus synchronization throughout the tissue. Further, cell-cell junctional forces exhibit distinct spatial patterns during collective cell migration, with larger forces near the leading edge. Larger junctional force magnitudes are associated with faster collective cell migration and larger tissue size. Simulations of heterogeneous tissue migration exhibit a complex dependence on the properties of both leading and trailing cells. Computational predictions demonstrate that collective cell migration depends on both the emergent dynamics and interactions between cellular-level Rho GTPase activity and contractility and multicellular-level junctional forces.

Role of Berry Anthocyanins and Phenolic Acids on Cell Migration and Angiogenesis: An Updated Overview

Cell migration is a critical process that is highly involved with normal and pathological conditions such as angiogenesis and wound healing. Important members of the RHO GTPase family are capable of controlling cytoskeleton conformation and altering motility characteristics of cells. There is a well-known relationship between small GTPases and the PI3K/AKT pathway. Endothelial cell migration can lead to angiogenesis, which is highly linked to wound healing processes. Phenolics, flavonoids, and anthocyanins are major groups of phytochemicals and are abundant in many natural products. Their antioxidant, antimicrobial, anti-inflammatory, antidiabetic, angiogenenic, neuroprotective, hepatoprotective, and cardioprotective properties have been extensively documented. This comprehensive review focuses on the in vitro and in vivo role of berry extracts and single anthocyanin and phenolic acid compounds on cell migration and angiogenesis. We aim to summarize the most recent published studies focusing on the experimental model, type of berry extract, source, dose/concentration and overall effect(s) of berry extracts, anthocyanins, and phenolic acids on the above processes.

Integration of in vitro and in silico Models Using Bayesian Optimization With an Application to Stochastic Modeling of Mesenchymal 3D Cell Migration

Cellular migration plays a crucial role in many aspects of life and development. In this paper, we propose a computational model of 3D migration that is solved by means of the tau-leaping algorithm and whose parameters have been calibrated using Bayesian optimization. Our main focus is two-fold: to optimize the numerical performance of the mechano-chemical model as well as to automate the calibration process of in silico models using Bayesian optimization. The presented mechano-chemical model allows us to simulate the stochastic behavior of our chemically reacting system in combination with mechanical constraints due to the surrounding collagen-based matrix. This numerical model has been used to simulate fibroblast migration. Moreover, we have performed in vitro analysis of migrating fibroblasts embedded in 3D collagen-based fibrous matrices (2 mg/ml). These in vitro experiments have been performed with the main objective of calibrating our model. Nine model parameters have been calibrated testing 300 different parametrizations using a completely automatic approach. Two competing evaluation metrics based on the Bhattacharyya coefficient have been defined in order to fit the model parameters. These metrics evaluate how accurately the in silico model is replicating in vitro measurements regarding the two main variables quantified in the experimental data (number of protrusions and the length of the longest protrusion). The selection of an optimal parametrization is based on the balance between the defined evaluation metrics. Results show how the calibrated model is able to predict the main features observed in the in vitro experiments.

A local uPAR-plasmin-TGFβ1 positive feedback loop in a qualitative computational model of angiogenic sprouting explains the in vitro effect of fibrinogen variants

In experimental assays of angiogenesis in three-dimensional fibrin matrices, a temporary scaffold formed during wound healing, the type and composition of fibrin impacts the level of sprouting. More sprouts form on high molecular weight (HMW) than on low molecular weight (LMW) fibrin. It is unclear what mechanisms regulate the number and the positions of the vascular-like structures in cell cultures. To address this question, we propose a mechanistic simulation model of endothelial cell migration and fibrin proteolysis by the plasmin system. The model is a hybrid, cell-based and continuum, computational model based on the cellular Potts model and sets of partial-differential equations. Based on the model results, we propose that a positive feedback mechanism between uPAR, plasmin and transforming growth factor β1 (TGFβ1) selects cells in the monolayer for matrix invasion. Invading cells releases TGFβ1 from the extracellular matrix through plasmin-mediated fibrin degradation. The activated TGFβ1 further stimulates fibrin degradation and keeps proteolysis active as the sprout invades the fibrin matrix. The binding capacity for TGFβ1 of LMW is reduced relative to that of HMW. This leads to reduced activation of proteolysis and, consequently, reduced cell ingrowth in LMW fibrin compared to HMW fibrin. Thus our model predicts that endothelial cells in LMW fibrin matrices compared to HMW matrices show reduced sprouting due to a lower bio-availability of TGFβ1.

Therapies for a range of medical conditions, including cancer, wound healing and diabetic retinopathy can benefit from a better control over the growth of blood vessels. The chemical properties of fibrin, the material that forms scabs in wounds and can also occur in large concentrations in tumors, can regulate the degree of blood vessel growth (angiogenesis). Angiogenesis can be mimicked in cell cultures. These allow us to modulate the chemical properties of fibrin and study the effect on angiogenesis. Fibrin occurs in high molecular weight (HMW) and in low molecular weight (LMW) forms. Interestingly, there is more ingrowth of angiogenic-like structures into HMW than in LMW fibrin, but the mechanisms are poorly understood. To get more insight into these, we constructed a computational model. Using the model, we propose and analyse a hypothetical mechanism for sprouting that could explain the differences in endothelial cell sprouting in LMW and HMW fibrin matrices. Our model suggests that cells digest fibrin, thus creating space for ingrowth. At the same time, digestion frees growth factors bound to fibrin, that activates further secretion of digestive enzymes by the cells. We propose that the resulting positive feedback loop spontaneously selects cells in the endothelial monolayer for ingrowth and helps the blood vessel sprout move deeper into the fibrin. This could be a complementary mechanism to lateral-inhibition by Delta-Notch for the selection of leader cells, also called ‘tip cells’. Our model predicts that endothelial cells in LMW fibrin compared to HMW fibrin show reduced sprouting due to a lower bio-availability of TGFβ1.

Human lung fibroblast-derived matrix facilitates vascular morphogenesis in 3D environment and enhances skin wound healing

Extracellular matrix (ECM) is crucial to many aspects of vascular morphogenesis and maintenance of vasculature function. Currently the recapitulation of angiogenic ECM microenvironment is still challenging, due mainly to its diverse components and complex organization. Here we investigate the angiogenic potential of human lung fibroblast-derived matrix (hFDM) in creating a three-dimensional (3D) vascular construct. hFDM was obtained via decellularization of in vitro cultured human lung fibroblasts and analyzed via immunofluorescence staining and ELISA, which detect multiple ECM macromolecules and angiogenic growth factors (GFs). Human umbilical vein endothelial cells (HUVECs) morphology was more elongated and better proliferative on hFDM than on gelatin-coated substrate. To prepare 3D construct, hFDM is collected, quantitatively analyzed, and incorporated in collagen hydrogel (Col) with HUVECs. Capillary-like structure (CLS) formation at 7day was significantly better with the groups containing higher doses of hFDM compared to the Col group (control). Moreover, the group (Col/hFDM/GFs) with both hFDM and angiogenic GFs (VEGF, bFGF, SDF-1) showed the synergistic activity on CLS formation and found much larger capillary lumen diameters with time. Further analysis of hFDM via angiogenesis antibody array kit reveals abundant biochemical cues, such as angiogenesis-related cytokines, GFs, and proteolytic enzymes. Significantly up-regulated expression of VE-cadherin and ECM-specific integrin subunits was also noticed in Col/hFDM/GFs. In addition, transplantation of Col/hFMD/GFs with HUVECs in skin wound model presents more effective re-epithelialization, many regenerated hair follicles, better transplanted cells viability, and advanced neovascularization. We believe that current system is a very promising platform for 3D vasculature construction in vitro and for cell delivery toward therapeutic applications in vivo.

STATEMENT OF SIGNIFICANCE

Functional 3D vasculature construction in vitro is still challenging due to the difficulty of recapitulating the complex angiogenic extracellular matrix (ECM) environment. Herein, we present a simple and practical method to create an angiogenic 3D environment via incorporation of human lung fibroblast-derived matrix (hFDM) into collagen hydrogel. We found that hFDM offers a significantly improved angiogenic microenvironment for HUVECs on 2D substrates and in 3D construct. A synergistic effect of hFDM and angiogenic growth factors has been well confirmed in 3D condition. The prevascularized 3D collagen constructs also facilitate skin wound healing. We believe that current system should be a convenient and powerful platform in engineering 3D vasculature in vitro, and in delivering cells for therapeutic purposes in vivo.

Spatiotemporal Dynamics of Insulitis in Human Type 1 Diabetes

Type 1 diabetes (T1D) is an auto-immune disease characterized by the selective destruction of the insulin secreting beta cells in the pancreas during an inflammatory phase known as insulitis. Patients with T1D are typically dependent on the administration of externally provided insulin in order to manage blood glucose levels. Whilst technological developments have significantly improved both the life expectancy and quality of life of these patients, an understanding of the mechanisms of the disease remains elusive. Animal models, such as the NOD mouse model, have been widely used to probe the process of insulitis, but there exist very few data from humans studied at disease onset. In this manuscript, we employ data from human pancreases collected close to the onset of T1D and propose a spatio-temporal computational model for the progression of insulitis in human T1D, with particular focus on the mechanisms underlying the development of insulitis in pancreatic islets. This framework allows us to investigate how the time-course of insulitis progression is affected by altering key parameters, such as the number of the CD20+ B cells present in the inflammatory infiltrate, which has recently been proposed to influence the aggressiveness of the disease. Through the analysis of repeated simulations of our stochastic model, which track the number of beta cells within an islet, we find that increased numbers of B cells in the peri-islet space lead to faster destruction of the beta cells. We also find that the balance between the degradation and repair of the basement membrane surrounding the islet is a critical component in governing the overall destruction rate of the beta cells and their remaining number. Our model provides a framework for continued and improved spatio-temporal modeling of human T1D.

Effects of endothelial cell proliferation and migration rates in a computational model of sprouting angiogenesis

Angiogenesis, the recruitment of new blood vessels, is a critical process for the growth, expansion, and metastatic dissemination of developing tumors. Three types of cells make up the new vasculature: tip cells, which migrate in response to gradients of vascular endothelial growth factor (VEGF), stalk cells, which proliferate and extend the vessels, and phalanx cells, which are quiescent and support the sprout. In this study we examine the contribution of tip cell migration rate and stalk cell proliferation rate on the formation of new vasculature. We calculate several vascular metrics, such as the number of vascular bifurcations per unit volume, vascular segment length per unit volume, and vascular tortuosity. These measurements predict that proliferation rate has a greater effect on the spread and extent of vascular growth compared to migration rate. Together, these findings provide strong implications for designing anti-angiogenic therapies that may differentially target endothelial cell proliferation and migration. Computational models can be used to predict optimal anti-angiogenic therapies in combination with other therapeutics to improve outcome.

Decellularized myocardial matrix hydrogels: In basic research and preclinical studies

A variety of decellularized materials have been developed that have demonstrated potential for treating cardiovascular diseases and improving our understanding of cardiac development. Of these biomaterials, decellularized myocardial matrix hydrogels have shown great promise for creating cellular microenvironments representative of the native cardiac tissue and treating the heart after a myocardial infarction. Decellularized myocardial matrix hydrogels derived from porcine cardiac tissue form a nanofibrous hydrogel once thermally induced at physiological temperatures. Use of isolated cardiac extracellular matrix in 2D and 3D in vitro platforms has demonstrated the capability to provide tissue specific cues for cardiac cell growth and differentiation. Testing of the myocardial matrix hydrogel as a therapy after myocardial infarction in both small and large animal models has demonstrated improved left ventricular function, increased cardiac muscle, and cellular recruitment into the treated infarct. Based on these results, steps are currently being taken to translate these hydrogels into a clinically used injectable biomaterial therapy. In this review, we will focus on the basic science and preclinical studies that have accelerated the development of decellularized myocardial matrix hydrogels into an emerging novel therapy for treating the heart after a myocardial infarction.

Tip cell overtaking occurs as a side effect of sprouting in computational models of angiogenesis

BACKGROUND

During angiogenesis, the formation of new blood vessels from existing ones, endothelial cells differentiate into tip and stalk cells, after which one tip cell leads the sprout. More recently, this picture has changed. It has become clear that endothelial cells compete for the tip position during angiogenesis: a phenomenon named tip cell overtaking. The biological function of tip cell overtaking is not yet known. From experimental observations, it is unclear to what extent tip cell overtaking is a side effect of sprouting or to what extent it is regulated through a VEGF-Dll4-Notch signaling network and thus might have a biological function. To address this question, we studied tip cell overtaking in computational models of angiogenic sprouting in absence and in presence of VEGF-Dll4-Notch signaling.

RESULTS

We looked for tip cell overtaking in two existing Cellular Potts models of angiogenesis. In these simulation models angiogenic sprouting-like behavior emerges from a small set of plausible cell behaviors. In the first model, cells aggregate through contact-inhibited chemotaxis. In the second model the endothelial cells assume an elongated shape and aggregate through (non-inhibited) chemotaxis. In both these sprouting models the endothelial cells spontaneously migrate forwards and backwards within sprouts, suggesting that tip cell overtaking might occur as a side effect of sprouting. In accordance with other experimental observations, in our simulations the cells' tendency to occupy the tip position can be regulated when two cell lines with different levels of Vegfr2 expression are contributing to sprouting (mosaic sprouting assay), where cell behavior is regulated by a simple VEGF-Dll4-Notch signaling network.

CONCLUSIONS

Our modeling results suggest that tip cell overtaking can occur spontaneously due to the stochastic motion of cells during sprouting. Thus, tip cell overtaking and sprouting dynamics may be interdependent and should be studied and interpreted in combination. VEGF-Dll4-Notch can regulate the ability of cells to occupy the tip cell position in our simulations. We propose that the function of VEGF-Dll4-Notch signaling might not be to regulate which cell ends up at the tip, but to assure that the cell that randomly ends up at the tip position acquires the tip cell phenotype.

Autonomy and Non-autonomy of Angiogenic Cell Movements Revealed by Experiment-Driven Mathematical Modeling

Angiogenesis is a multicellular phenomenon driven by morphogenetic cell movements. We recently reported morphogenetic vascular endothelial cell (EC) behaviors to be dynamic and complex. However, the principal mechanisms orchestrating individual EC movements in angiogenic morphogenesis remain largely unknown. Here we present an experiment-driven mathematical model that enables us to systematically dissect cellular mechanisms in branch elongation. We found that cell-autonomous and coordinated actions governed these multicellular behaviors, and a cell-autonomous process sufficiently illustrated essential features of the morphogenetic EC dynamics at both the single-cell and cell-population levels. Through refining our model and experimental verification, we further identified a coordinated mode of tip EC behaviors regulated via a spatial relationship between tip and follower ECs, which facilitates the forward motility of tip ECs. These findings provide insights that enhance our mechanistic understanding of not only angiogenic morphogenesis, but also other types of multicellular phenomenon.