Contact-inhibited chemotaxis in de novo and sprouting blood-vessel growth

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.

Agent-based model of angiogenesis simulates capillary sprout initiation in multicellular networks

Many biological processes are controlled by both deterministic and stochastic influences. However, efforts to model these systems often rely on either purely stochastic or purely rule-based methods. To better understand the balance between stochasticity and determinism in biological processes a computational approach that incorporates both influences may afford additional insight into underlying biological mechanisms that give rise to emergent system properties. We apply a combined approach to the simulation and study of angiogenesis, the growth of new blood vessels from existing networks. This complex multicellular process begins with selection of an initiating endothelial cell, or tip cell, which sprouts from the parent vessels in response to stimulation by exogenous cues. We have constructed an agent-based model of sprouting angiogenesis to evaluate endothelial cell sprout initiation frequency and location, and we have experimentally validated it using high-resolution time-lapse confocal microscopy. ABM simulations were then compared to a Monte Carlo model, revealing that purely stochastic simulations could not generate sprout locations as accurately as the rule-informed agent-based model. These findings support the use of rule-based approaches for modeling the complex mechanisms underlying sprouting angiogenesis over purely stochastic methods.

The Force at the Tip--Modelling Tension and Proliferation in Sprouting Angiogenesis

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.

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.

Dynamic mechanical finite element model of biological cells for studying cellular pattern formation

Understanding the geometric, topologic, and mechanical properties of cells and their interactions is critical for studying tissue pattern formation and organ development. Computational model and tools for simulating cell pattern formation have broad implications in studying embryogenesis, blood-vessel development, tissue regeneration, and tumor growth. Although a number of cell modeling methods exist, they do not simultaneously account for detailed cellular shapes as well as dynamic changes in cell geometry and topology. Here we describe a dynamic finite element cell model (dFEMC) for studying populations of cells and tissue development. By incorporating details of cell shape, cell growth and shrinkage, cell birth and death, cell division and fusion, our method can model realistically a variety problems of cell pattern formation. We give two examples of applying our method to the study of cell fusion and cell apoptosis. The dFEMC model developed here provides a general computational framework for studying dynamics pattern formation of tissue.

Emergent Stratification in Solid Tumors Selects for Reduced Cohesion of Tumor Cells: A Multi-Cell, Virtual-Tissue Model of Tumor Evolution Using CompuCell3D

Tumor cells and structure both evolve due to heritable variation of cell behaviors and selection over periods of weeks to years (somatic evolution). Micro-environmental factors exert selection pressures on tumor-cell behaviors, which influence both the rate and direction of evolution of specific behaviors, especially the development of tumor-cell aggression and resistance to chemotherapies. In this paper, we present, step-by-step, the development of a multi-cell, virtual-tissue model of tumor somatic evolution, simulated using the open-source CompuCell3D modeling environment. Our model includes essential cell behaviors, microenvironmental components and their interactions. Our model provides a platform for exploring selection pressures leading to the evolution of tumor-cell aggression, showing that emergent stratification into regions with different cell survival rates drives the evolution of less cohesive cells with lower levels of cadherins and higher levels of integrins. Such reduced cohesivity is a key hallmark in the progression of many types of solid tumors.

The Vital Role of Blood Flow-Induced Proliferation and Migration in Capillary Network Formation in a Multiscale Model of Angiogenesis

Sprouting angiogenesis and capillary network formation are tissue scale phenomena. There are also sub-scale phenomena involved in angiogenesis including at the cellular and intracellular (molecular) scales. In this work, a multiscale model of angiogenesis spanning intracellular, cellular, and tissue scales is developed in detail. The key events that are considered at the tissue scale are formation of closed flow path (that is called loop in this article) and blood flow initiation in the loop. At the cellular scale, growth, migration, and anastomosis of endothelial cells (ECs) are important. At the intracellular scale, cell phenotype determination as well as alteration due to blood flow is included, having pivotal roles in the model. The main feature of the model is to obtain the physical behavior of a closed loop at the tissue scale, relying on the events at the cellular and intracellular scales, and not by imposing physical behavior upon it. Results show that, when blood flow is considered in the loop, the anastomosed sprouts stabilize and elongate. By contrast, when the loop is modeled without consideration of blood flow, the loop collapses. The results obtained in this work show that proper determination of EC phenotype is the key for its survival.

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

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

Collective motion of cells: from experiments to models

Swarming or collective motion of living entities is one of the most common and spectacular manifestations of living systems that have been extensively studied in recent years. A number of general principles have been established. The interactions at the level of cells are quite different from those among individual animals, therefore the study of collective motion of cells is likely to reveal some specific important features which we plan to overview in this paper. In addition to presenting the most appealing results from the quickly growing related literature we also deliver a critical discussion of the emerging picture and summarize our present understanding of collective motion at the cellular level. Collective motion of cells plays an essential role in a number of experimental and real-life situations. In most cases the coordinated motion is a helpful aspect of the given phenomenon and results in making a related process more efficient (e.g., embryogenesis or wound healing), while in the case of tumor cell invasion it appears to speed up the progression of the disease. In these mechanisms cells both have to be motile and adhere to one another, the adherence feature being the most specific to this sort of collective behavior. One of the central aims of this review is to present the related experimental observations and treat them in light of a few basic computational models so as to make an interpretation of the phenomena at a quantitative level as well.

Systems biology of the microvasculature

The vascular network carries blood throughout the body, delivering oxygen to tissues and providing a pathway for communication between distant organs. The network is hierarchical and structured, but also dynamic, especially at the smaller scales. Remodeling of the microvasculature occurs in response to local changes in oxygen, gene expression, cell-cell communication, and chemical and mechanical stimuli from the microenvironment. These local changes occur as a result of physiological processes such as growth and exercise, as well as acute and chronic diseases including stroke, cancer, and diabetes, and pharmacological intervention. While the vasculature is an important therapeutic target in many diseases, drugs designed to inhibit vascular growth have achieved only limited success, and no drug has yet been approved to promote therapeutic vascular remodeling. This highlights the challenges involved in identifying appropriate therapeutic targets in a system as complex as the vasculature. Systems biology approaches provide a means to bridge current understanding of the vascular system, from detailed signaling dynamics measured in vitro and pre-clinical animal models of vascular disease, to a more complete picture of vascular regulation in vivo. This will translate to an improved ability to identify multi-component biomarkers for diagnosis, prognosis, and monitoring of therapy that are easy to measure in vivo, as well as better drug targets for specific disease states. In this review, we summarize systems biology approaches that have advanced our understanding of vascular function and dysfunction in vivo, with a focus on computational modeling.

Collective migration exhibits greater sensitivity but slower dynamics of alignment to applied electric fields

During development and disease, cells migrate collectively in response to gradients in physical, chemical and electrical cues. Despite its physiological significance and potential therapeutic applications, electrotactic collective cell movement is relatively less well understood. Here, we analyze the combined effect of intercellular interactions and electric fields on the directional migration of non-transformed mammary epithelial cells, MCF-10A. Our data show that clustered cells exhibit greater sensitivity to applied electric fields but align more slowly than isolated cells. Clustered cells achieve half-maximal directedness with an electric field that is 50% weaker than that required by isolated cells; however, clustered cells take ∼2-4 fold longer to align. This trade-off in greater sensitivity and slower dynamics correlates with the slower speed and intrinsic directedness of collective movement even in the absence of an electric field. Whereas isolated cells exhibit a persistent random walk, the trajectories of clustered cells are more ballistic as evidenced by the superlinear dependence of their mean square displacement on time. Thus, intrinsically-directed, slower clustered cells take longer to redirect and align with an electric field. These findings help to define the operating space and the engineering trade-offs for using electric fields to affect cell movement in biomedical applications.

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

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

Large-scale parameter studies of cell-based models of tissue morphogenesis using CompuCell3D or VirtualLeaf

Computational, cell-based models, such as the cellular Potts model (CPM), have become a widely used tool to study tissue formation. Most cell-based models mimic the physical properties of cells and their dynamic behavior, and generate images of the tissue that the cells form due to their collective behavior. Due to these intuitive parameters and output, cell-based models are often evaluated visually and the parameters are fine-tuned by hand. To get better insight into how in a cell-based model the microscopic scale (e.g., cell behavior, secreted molecular signals, and cell-ECM interactions) determines the macroscopic scale, we need to generate morphospaces and perform parameter sweeps, involving large numbers of individual simulations. This chapter describes a protocol and presents a set of scripts for automatically setting up, running, and evaluating large-scale parameter sweeps of cell-based models. We demonstrate the use of the protocol using a recent cellular Potts model of blood vessel formation model implemented in CompuCell3D. We show the versatility of the protocol by adapting it to an alternative cell-based modeling framework, VirtualLeaf.

Cell-based computational modeling of vascular morphogenesis using Tissue Simulation Toolkit

Computational modeling has become a widely used tool for unraveling the mechanisms of higher level cooperative cell behavior during vascular morphogenesis. However, experimenting with published simulation models or adding new assumptions to those models can be daunting for novice and even for experienced computational scientists. Here, we present a step-by-step, practical tutorial for building cell-based simulations of vascular morphogenesis using the Tissue Simulation Toolkit (TST). The TST is a freely available, open-source C++ library for developing simulations with the two-dimensional cellular Potts model, a stochastic, agent-based framework to simulate collective cell behavior. We will show the basic use of the TST to simulate and experiment with published simulations of vascular network formation. Then, we will present step-by-step instructions and explanations for building a recent simulation model of tumor angiogenesis. Demonstrated mechanisms include cell-cell adhesion, chemotaxis, cell elongation, haptotaxis, and haptokinesis.

Direct measurement of cell protrusion force utilizing a robot-aided cell manipulation system with optical tweezers for cell migration control

Cell migration refers to the directional cell movement in response to a chemoattractant gradient, a key process that occurs in a wide variety of biological phenomena. Cell protrusion force is generated by the actin polymerization of a cell, which drives the cell to move toward the stimulus as induced by the chemoattractant gradient. This paper presents a new methodology for the direct measurement of cell protrusion force utilizing a robot-aided optical tweezer system. The functionalized beads that are robotically trapped and placed near the cell serve as both cell migration stimulators and protrusion force probes. The force generated by the actin polymerization of the cell propels the bead to move away from the trapping center when the cell comes in contact with the bead. Such a deviation can be determined and used to calculate the trapping force, which is equal to the protrusion force at a balanced position. With the quantitative measurement of the protrusion, we find that the protrusion force of a live cell in response to a chemoattractant within the range of hundreds of piconewtons. We further probe the protrusion force distribution at the cell leading edge and find that the highest protrusion force appears at the cell migration direction. These measurements can help us characterize the mechanism of cell migration and lay a solid foundation for further proactive control of cell movement.

Mechanical cell-matrix feedback explains pairwise and collective endothelial cell behavior in vitro

In vitro cultures of endothelial cells are a widely used model system of the collective behavior of endothelial cells during vasculogenesis and angiogenesis. When seeded in an extracellular matrix, endothelial cells can form blood vessel-like structures, including vascular networks and sprouts. Endothelial morphogenesis depends on a large number of chemical and mechanical factors, including the compliancy of the extracellular matrix, the available growth factors, the adhesion of cells to the extracellular matrix, cell-cell signaling, etc. Although various computational models have been proposed to explain the role of each of these biochemical and biomechanical effects, the understanding of the mechanisms underlying in vitro angiogenesis is still incomplete. Most explanations focus on predicting the whole vascular network or sprout from the underlying cell behavior, and do not check if the same model also correctly captures the intermediate scale: the pairwise cell-cell interactions or single cell responses to ECM mechanics. Here we show, using a hybrid cellular Potts and finite element computational model, that a single set of biologically plausible rules describing (a) the contractile forces that endothelial cells exert on the ECM, (b) the resulting strains in the extracellular matrix, and (c) the cellular response to the strains, suffices for reproducing the behavior of individual endothelial cells and the interactions of endothelial cell pairs in compliant matrices. With the same set of rules, the model also reproduces network formation from scattered cells, and sprouting from endothelial spheroids. Combining the present mechanical model with aspects of previously proposed mechanical and chemical models may lead to a more complete understanding of in vitro angiogenesis.

In Vitro Microvessel Growth and Remodeling within a Three-dimensional Microfluidic Environment

This paper presents in vitro microvascular network formation within 3D gel scaffolds made from different concentrations of type-I collagen, fibrin, or a mixture of collagen and fibrin, using a simple microfluidic platform. Initially, microvascular network formation of human umbilical vein endothelial cells was examined using live time-lapse confocal microscopy every 90 min from 3 h to 12 h after seeding within three different concentrations of collagen gel scaffolds. Among the three conditions of collagen gel scaffolds (2.0 mg/ml, 2.5 mg/ml, and 3.0 mg/ml), the number of skeleton within collagen gel scaffolds was consistently the highest (3.0 mg/ml), followed by those of collagen gel scaffolds (2.5 mg/ml and 2.0 mg/ml). Results demonstrated that concentration of collagen gel scaffolds, which influences matrix stiffness and ligand density, may affect microvascular network formation during the early stages of vasculogenesis. In addition, the maturation of microvascular networks in monoculture under different gel compositions within gel scaffolds (2.5 mg/ml) was examined for 7 d using live confocal microscopy. It was confirmed that pure fibrin gel scaffolds are preferable to collagen gel or collagen/fibrin combinations, significantly reducing matrix retractions during maturation of microvascular networks for 7 d. Finally, early steps in the maturation process of microvascular networks for 14 d were characterized by demonstrating sequential steps of branching, expanding, remodeling, pruning, and clear delineation of lumens within fibrin gel scaffolds. Our findings demonstrate an in vitro model for generating mature microvascular networks within 3D microfluidic fibrin gel scaffolds (2.5 mg/ml), and furthermore suggest the importance of gel concentration and composition in promoting the maturation of microvascular networks.

Cell-based multi-parametric model of cleft progression during submandibular salivary gland branching morphogenesis

Cleft formation during submandibular salivary gland branching morphogenesis is the critical step initiating the growth and development of the complex adult organ. Previous experimental studies indicated requirements for several epithelial cellular processes, such as proliferation, migration, cell-cell adhesion, cell-extracellular matrix (matrix) adhesion, and cellular contraction in cleft formation; however, the relative contribution of each of these processes is not fully understood since it is not possible to experimentally manipulate each factor independently. We present here a comprehensive analysis of several cellular parameters regulating cleft progression during branching morphogenesis in the epithelial tissue of an early embryonic salivary gland at a local scale using an on lattice Monte-Carlo simulation model, the Glazier-Graner-Hogeweg model. We utilized measurements from time-lapse images of mouse submandibular gland organ explants to construct a temporally and spatially relevant cell-based 2D model. Our model simulates the effect of cellular proliferation, actomyosin contractility, cell-cell and cell-matrix adhesions on cleft progression, and it was used to test specific hypotheses regarding the function of these parameters in branching morphogenesis. We use innovative features capturing several aspects of cleft morphology and quantitatively analyze clefts formed during functional modification of the cellular parameters. Our simulations predict that a low epithelial mitosis rate and moderate level of actomyosin contractility in the cleft cells promote cleft progression. Raising or lowering levels of contractility and mitosis rate resulted in non-progressive clefts. We also show that lowered cell-cell adhesion in the cleft region and increased cleft cell-matrix adhesions are required for cleft progression. Using a classifier-based analysis, the relative importance of these four contributing cellular factors for effective cleft progression was determined as follows: cleft cell contractility, cleft region cell-cell adhesion strength, epithelial cell mitosis rate, and cell-matrix adhesion strength.

Branching morphogenesis is a complex and dynamic embryonic process that creates the structure of many adult organs, including the salivary gland. During this process, many cellular changes occur in the epithelial cells, including changes in cell-cell adhesions, cell-extracellular matrix (matrix) adhesions, cell proliferation, and cellular contraction, resulting in formation of clefts in the epithelial cells of the organ. A comprehensive understanding of the relative contributions of these cellular processes has crucial therapeutic implications for organ regeneration and functional restoration of organ structure in diseased salivary glands. Here, we have developed a cell-based model of cleft progression and simulated cleft progression under conditions of altered cell-cell adhesions, cellular contractility, cell-matrix adhesion and cell proliferation to identify the optimum cellular conditions that cause clefts to progress. The model predicts that cleft progression requires a moderate level of cleft cell contractility, a low epithelial proliferation rate, reduced cell-cell adhesion strength in the cleft and high cell-matrix adhesion strength also in the cleft region. The results of our classification analysis demonstrate that cellular contractility in the cleft cells has a significant effect on cleft progression, followed by cell-cell adhesion strength, rate of cell proliferation, and strength of cell-matrix adhesion energies.

A systems biology view of blood vessel growth and remodelling

Blood travels throughout the body in an extensive network of vessels – arteries, veins and capillaries. This vascular network is not static, but instead dynamically remodels in response to stimuli from cells in the nearby tissue. In particular, the smallest vessels – arterioles, venules and capillaries – can be extended, expanded or pruned, in response to exercise, ischaemic events, pharmacological interventions, or other physiological and pathophysiological events. In this review, we describe the multi-step morphogenic process of angiogenesis – the sprouting of new blood vessels – and the stability of vascular networks in vivo. In particular, we review the known interactions between endothelial cells and the various blood cells and plasma components they convey. We describe progress that has been made in applying computational modelling, quantitative biology and high-throughput experimentation to the angiogenesis process.

Simulating cortical development as a self constructing process: a novel multi-scale approach combining molecular and physical aspects

Current models of embryological development focus on intracellular processes such as gene expression and protein networks, rather than on the complex relationship between subcellular processes and the collective cellular organization these processes support. We have explored this collective behavior in the context of neocortical development, by modeling the expansion of a small number of progenitor cells into a laminated cortex with layer and cell type specific projections. The developmental process is steered by a formal language analogous to genomic instructions, and takes place in a physically realistic three-dimensional environment. A common genome inserted into individual cells control their individual behaviors, and thereby gives rise to collective developmental sequences in a biologically plausible manner. The simulation begins with a single progenitor cell containing the artificial genome. This progenitor then gives rise through a lineage of offspring to distinct populations of neuronal precursors that migrate to form the cortical laminae. The precursors differentiate by extending dendrites and axons, which reproduce the experimentally determined branching patterns of a number of different neuronal cell types observed in the cat visual cortex. This result is the first comprehensive demonstration of the principles of self-construction whereby the cortical architecture develops. In addition, our model makes several testable predictions concerning cell migration and branching mechanisms.

The proper operation of the brain depends on the correct developmental wiring of billions of neurons. Understanding this process of living self-construction is crucial not only for biological explanation and medical therapy, but could also provide an entirely new approach to industrial fabrication. We are approaching this problem through detailed simulation of cortical development. We have previously presented a software package that allows for simulation of cellular growth in a 3D space that respects physical forces and diffusion of substances, as well as an instruction language for specifying biologically plausible ‘genetic codes’. Here we apply this novel formalism to understanding the principles of cortical development in the context of multiple, spatially distributed agents that communicate only by local metabolic messages.

Endothelial cell motility, coordination and pattern formation during vasculogenesis

How vascular networks assemble is a fundamental problem of developmental biology that also has medical importance. To explain the organizational principles behind vascular patterning, we must understand how can tissue level structures be controlled through cell behavior patterns like motility and adhesion that, in turn, are determined by biochemical signal transduction processes? We discuss the various ideas that have been proposed as mechanisms for vascular network assembly: cell motility guided by extracellular matrix alignment (contact guidance), chemotaxis guided by paracrine and autocrine morphogens, and multicellular sprouting guided by cell-cell contacts. All of these processes yield emergent patterns, thus endothelial cells can form an interconnected structure autonomously, without guidance from an external pre-pattern.

Collective cell streams in epithelial monolayers depend on cell adhesion

We report a spontaneously emerging, randomly oriented, collective streaming behavior within a monolayer culture of a human keratinocyte cell line, and explore the effect of modulating cell adhesions by perturbing the function of calcium-dependent cell adhesion molecules. We demonstrate that decreasing cell adhesion induces narrower and more anisotropic cell streams, reminiscent of decreasing the Taylor scale of turbulent liquids. To explain our empirical findings, we propose a cell-based model that represents the dual nature of cell-cell adhesions. Spring-like connections provide mechanical stability, while a cellular Potts model formalism represents surface-tension driven attachment. By changing the relevance and persistence of mechanical links between cells, we are able to explain the experimentally observed changes in emergent flow patterns.

Role of differential adhesion in cell cluster evolution: from vasculogenesis to cancer metastasis

Cell-cell and cell-matrix adhesions are fundamental to numerous physiological processes, including angiogenesis, tumourigenesis, metastatic spreading and wound healing. We use cellular potts model to computationally predict the organisation of cells within a 3D matrix. The energy potentials regulating cell-cell (JCC) and cell-matrix (JMC) adhesive interactions are systematically varied to represent different, biologically relevant adhesive conditions. Chemotactically induced cell migration is also addressed. Starting from a cluster of cells, variations in relative cell adhesion alone lead to different cellular patterns such as spreading of metastatic tumours and angiogenesis. The combination of low cell-cell adhesion (high JCC) and high heterotypic adhesion (low JMC) favours the fragmentation of the original cluster into multiple, smaller cell clusters (metastasis). Conversely, cellular systems exhibiting high-homotypic affinity (low JCC) preserve their original configuration, avoiding fragmentation (organogenesis). For intermediate values of JCC and JMC (i.e. JCC/JMC ∼ 1), tubular and corrugated structures form. Fully developed vascular trees are assembled only in systems in which contact-inhibited chemotaxis is activated upon cell contact. Also, the rate of secretion, diffusion and sequestration of chemotactic factors, cell deformability and motility do not significantly affect these trends. Further developments of this computational model will predict the efficacy of therapeutic interventions to modulate the diseased microenvironment by directly altering cell cohesion.

Cellular potts modeling of tumor growth, tumor invasion, and tumor evolution

Despite a growing wealth of available molecular data, the growth of tumors, invasion of tumors into healthy tissue, and response of tumors to therapies are still poorly understood. Although genetic mutations are in general the first step in the development of a cancer, for the mutated cell to persist in a tissue, it must compete against the other, healthy or diseased cells, for example by becoming more motile, adhesive, or multiplying faster. Thus, the cellular phenotype determines the success of a cancer cell in competition with its neighbors, irrespective of the genetic mutations or physiological alterations that gave rise to the altered phenotype. What phenotypes can make a cell "successful" in an environment of healthy and cancerous cells, and how? A widely used tool for getting more insight into that question is cell-based modeling. Cell-based models constitute a class of computational, agent-based models that mimic biophysical and molecular interactions between cells. One of the most widely used cell-based modeling formalisms is the cellular Potts model (CPM), a lattice-based, multi particle cell-based modeling approach. The CPM has become a popular and accessible method for modeling mechanisms of multicellular processes including cell sorting, gastrulation, or angiogenesis. The CPM accounts for biophysical cellular properties, including cell proliferation, cell motility, and cell adhesion, which play a key role in cancer. Multiscale models are constructed by extending the agents with intracellular processes including metabolism, growth, and signaling. Here we review the use of the CPM for modeling tumor growth, tumor invasion, and tumor progression. We argue that the accessibility and flexibility of the CPM, and its accurate, yet coarse-grained and computationally efficient representation of cell and tissue biophysics, make the CPM the method of choice for modeling cellular processes in tumor development.

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

Angiogenesis, the formation of new blood vessels sprouting from existing ones, occurs in several situations like wound healing, tissue remodeling, and near growing tumors. Under hypoxic conditions, tumor cells secrete growth factors, including VEGF. VEGF activates endothelial cells (ECs) in nearby vessels, leading to the migration of ECs out of the vessel and the formation of growing sprouts. A key process in angiogenesis is cellular self-organization, and previous modeling studies have identified mechanisms for producing networks and sprouts. Most theoretical studies of cellular self-organization during angiogenesis have ignored the interactions of ECs with the extra-cellular matrix (ECM), the jelly or hard materials that cells live in. Apart from providing structural support to cells, the ECM may play a key role in the coordination of cellular motility during angiogenesis. For example, by modifying the ECM, ECs can affect the motility of other ECs, long after they have left. Here, we present an explorative study of the cellular self-organization resulting from such ECM-coordinated cell migration. We show that a set of biologically-motivated, cell behavioral rules, including chemotaxis, haptotaxis, haptokinesis, and ECM-guided proliferation suffice for forming sprouts and branching vascular trees.

A Cellular Potts Model simulating cell migration on and in matrix environments

Cell migration on and through extracellular matrix is fundamental in a wide variety of physiological and pathological phenomena, and is exploited in scaffold-based tissue engineering. Migration is regulated by a number of extracellular matrix- or cell-derived biophysical parameters, such as matrix fiber orientation, pore size, and elasticity, or cell deformation, proteolysis, and adhesion. We here present an extended Cellular Potts Model (CPM) able to qualitatively and quantitatively describe cell migration efficiencies and phenotypes both on two-dimensional substrates and within three-dimensional matrices, close to experimental evidence. As distinct features of our approach, cells are modeled as compartmentalized discrete objects, differentiated into nucleus and cytosolic region, while the extracellular matrix is composed of a fibrous mesh and a homogeneous fluid. Our model provides a strong correlation of the directionality of migration with the topological extracellular matrix distribution and a biphasic dependence of migration on the matrix structure, density, adhesion, and stiffness, and, moreover, simulates that cell locomotion in highly constrained fibrillar obstacles requires the deformation of the cell's nucleus and/or the activity of cell-derived proteolysis. In conclusion, we here propose a mathematical modeling approach that serves to characterize cell migration as a biological phenomenon in healthy and diseased tissues and in engineering applications.

Vascular networks due to dynamically arrested crystalline ordering of elongated cells

Recent experimental and theoretical studies suggest that crystallization and glass-like solidification are useful analogies for understanding cell ordering in confluent biological tissues. It remains unexplored how cellular ordering contributes to pattern formation during morphogenesis. With a computational model we show that a system of elongated, cohering biological cells can get dynamically arrested in a network pattern. Our model provides an explanation for the formation of cellular networks in culture systems that exclude intercellular interaction via chemotaxis or mechanical traction.

Comparison of mixed and lamellar coculture spatial arrangements for tissue engineering capillary networks in vitro

Coculture of endothelial cells (ECs) and smooth muscle cells (SMCs) in vitro can yield confluent monolayers or EC networks. Factors influencing this transition are not known. In this study, we examined whether the spatial arrangement of EC-SMC cocultures affected EC migration, network morphology, and angiogenic protein secretion. Human umbilical cord blood-derived ECs (hCB-ECs) were grown in coculture with human aortic SMCs in either a mixed or lamellar spatial geometry and analyzed over a culture period of 12 days. The hCB-ECs cultured on SMCs in a mixed system had higher cell speeds, shorter persistence times, and lower random motility coefficients than ECs in a lamellar system. By day 12 of coculture, mixed systems demonstrated greater anastomoses and capillary loop formation than lamellar systems as evidenced by a higher number of branch points, angle of curvature between branch points, and percentage of imaged area covered by networks. The network morphology was more uniform in the mixed systems than the lamellar systems with fewer EC clusters present after several days in culture. Proliferation of hCB-ECs was higher for mixed cocultures during the first 24 h of coculture, and then declined dramatically suggesting that proliferation only contributed to network formation during the early stages of coculture. Proteome assay results show reduced solution levels, but no change in intracellular levels of angiogenic proteins in lamellar systems compared to mixed systems. These data suggest that mixing ECs and SMCs together favors the formation of EC networks to a greater extent than a lamellar arrangement in which ECs form a cell layer above a confluent, quiescent layer of SMCs.

Chemotactic migration of T cells towards dendritic cells promotes the detection of rare antigens

In many immunological processes chemoattraction is thought to play a role in guiding cells to their sites of action. However, based on in vivo two-photon microscopy experiments in the absence of cognate antigen, T cell migration in lymph nodes (LNs) has been roughly described as a random walk. Although it has been shown that dendritic cells (DCs) carrying cognate antigen in some circumstances attract T cells chemotactically, it is currently still unclear whether chemoattraction of T cells towards DCs helps or hampers scanning. Chemoattraction towards DCs could on the one hand help T cells to rapidly find DCs. On the other hand, it could be deleterious if DCs become shielded by a multitude of attracted yet non-specific T cells. Results from a recent simulation study suggested that the deleterious effect dominates. We re-addressed the question whether T cell chemoattraction towards DCs is expected to promote or hamper the detection of rare antigens using the Cellular Potts Model, a formalism that allows for dynamic, flexible cellular shapes and cell migration. Our simulations show that chemoattraction of T cells enhances the DC scanning efficiency, leading to an increased probability that rare antigen-specific T cells find DCs carrying cognate antigen. Desensitization of T cells after contact with a DC further improves the scanning efficiency, yielding an almost threefold enhancement compared to random migration. Moreover, the chemotaxis-driven migration still roughly appears as a random walk, hence fine-tuned analysis of cell tracks will be required to detect chemotaxis within microscopy data.

Modeling the influence of nucleus elasticity on cell invasion in fiber networks and microchannels

Cell migration in highly constrained extracellular matrices is exploited in scaffold-based tissue engineering and is fundamental in a wide variety of physiological and pathological phenomena, among others in cancer invasion and development. Research into the critical processes involved in cell migration has mainly focused on cell adhesion and proteolytic degradation of the external environment. However, rising evidence has recently shown that a number of cell-derived biophysical and mechanical parameters, among others nucleus stiffness and cell deformability, plays a major role in cell motility, especially in the ameboid-like migration mode in 3D confined tissue structures. We here present an extended cellular Potts model (CPM) first used to simulate a micro-fabricated migration chip, which tests the active invasive behavior of cancer cells into narrow channels. As distinct features of our approach, cells are modeled as compartmentalized discrete objects, differentiated in the nucleus and in the cytosolic region, while the migration chamber is composed of channels of different widths. We find that cell motile phenotype and velocity in open spaces (i.e., 2D flat surfaces or large channels) are not significantly influenced by cell elastic properties. On the contrary, the migratory behavior of cells within subcellular and subnuclear structures strongly relies on the deformability of the cytosol and of the nuclear cluster, respectively. Further, we characterize two migration dynamics: a stepwise way, characterized by fluctuations in cell length, within channels smaller than nucleus dimensions and a smooth sliding (i.e., maintaining constant cell length) behavior within channels larger than the nuclear cluster. These resulting observations are then extended looking at cell migration in an artificial fiber network, which mimics cell invasion in a 3D extracellular matrix. In particular, in this case, we analyze the effect of variations in elasticity of the nucleus on cell movement. In order to summarize, with our simulated migration assays, we demonstrate that the dimensionality of the environment strongly affects the migration phenotype and we suggest that the cytoskeletal and nuclear elastic characteristics correlate with the tumor cell's invasive potential.

Multi-scale modeling in morphogenesis: a critical analysis of the cellular Potts model

Cellular Potts models (CPMs) are used as a modeling framework to elucidate mechanisms of biological development. They allow a spatial resolution below the cellular scale and are applied particularly when problems are studied where multiple spatial and temporal scales are involved. Despite the increasing usage of CPMs in theoretical biology, this model class has received little attention from mathematical theory. To narrow this gap, the CPMs are subjected to a theoretical study here. It is asked to which extent the updating rules establish an appropriate dynamical model of intercellular interactions and what the principal behavior at different time scales characterizes. It is shown that the longtime behavior of a CPM is degenerate in the sense that the cells consecutively die out, independent of the specific interdependence structure that characterizes the model. While CPMs are naturally defined on finite, spatially bounded lattices, possible extensions to spatially unbounded systems are explored to assess to which extent spatio-temporal limit procedures can be applied to describe the emergent behavior at the tissue scale. To elucidate the mechanistic structure of CPMs, the model class is integrated into a general multiscale framework. It is shown that the central role of the surface fluctuations, which subsume several cellular and intercellular factors, entails substantial limitations for a CPM's exploitation both as a mechanistic and as a phenomenological model.

Multilevel complexity of calcium signaling: Modeling angiogenesis

Intracellular calcium signaling is a universal, evolutionary conserved and versatile regulator of cell biochemistry. The complexity of calcium signaling and related cell machinery can be investigated by the use of experimental strategies, as well as by computational approaches. Vascular endothelium is a fascinating model to study the specific properties and roles of calcium signals at multiple biological levels. During the past 20 years, live cell imaging, patch clamp and other techniques have allowed us to detect and interfere with calcium signaling in endothelial cells (ECs), providing a huge amount of information on the regulation of vascularization (angiogenesis) in normal and tumoral tissues. These data range from the spatiotemporal dynamics of calcium within different cell microcompartments to those in entire multicellular and organized EC networks. Beside experimental strategies, in silico endothelial models, specifically designed for simulating calcium signaling, are contributing to our knowledge of vascular physiology and pathology. They help to investigate and predict the quantitative features of proangiogenic events moving through subcellular, cellular and supracellular levels. This review focuses on some recent developments of computational approaches for proangiogenic endothelial calcium signaling. In particular, we discuss the creation of hybrid simulation environments, which combine and integrate discrete Cellular Potts Models. They are able to capture the phenomenological mechanisms of cell morphological reorganization, migration, and intercellular adhesion, with single-cell spatiotemporal models, based on reaction-diffusion equations that describe the agonist-induced intracellular calcium events.

Pattern formation during vasculogenesis

Vasculogenesis, the assembly of the first vascular network, is an intriguing developmental process that yields the first functional organ system of the embryo. In addition to being a fundamental part of embryonic development, vasculogenic processes also have medical importance. To explain the organizational principles behind vascular patterning, we must understand how morphogenesis of tissue level structures can be controlled through cell behavior patterns that, in turn, are determined by biochemical signal transduction processes. Mathematical analyses and computer simulations can help conceptualize how to bridge organizational levels and thus help in evaluating hypotheses regarding the formation of vascular networks. Here, we discuss the ideas that have been proposed to explain the formation of the first vascular pattern: cell motility guided by extracellular matrix alignment (contact guidance), chemotaxis guided by paracrine and autocrine morphogens, and sprouting guided by cell-cell contacts.

Adhesion failures determine the pattern of choroidal neovascularization in the eye: a computer simulation study

Choroidal neovascularization (CNV) of the macular area of the retina is the major cause of severe vision loss in adults. In CNV, after choriocapillaries initially penetrate Bruch's membrane (BrM), invading vessels may regress or expand (CNV initiation). Next, during Early and Late CNV, the expanding vasculature usually spreads in one of three distinct patterns: in a layer between BrM and the retinal pigment epithelium (sub-RPE or Type 1 CNV), in a layer between the RPE and the photoreceptors (sub-retinal or Type 2 CNV) or in both loci simultaneously (combined pattern or Type 3 CNV). While most studies hypothesize that CNV primarily results from growth-factor effects or holes in BrM, our three-dimensional simulations of multi-cell model of the normal and pathological maculae recapitulate the three growth patterns, under the hypothesis that CNV results from combinations of impairment of: 1) RPE-RPE epithelial junctional adhesion, 2) Adhesion of the RPE basement membrane complex to BrM (RPE-BrM adhesion), and 3) Adhesion of the RPE to the photoreceptor outer segments (RPE-POS adhesion). Our key findings are that when an endothelial tip cell penetrates BrM: 1) RPE with normal epithelial junctions, basal attachment to BrM and apical attachment to POS resists CNV. 2) Small holes in BrM do not, by themselves, initiate CNV. 3) RPE with normal epithelial junctions and normal apical RPE-POS adhesion, but weak adhesion to BrM (e.g. due to lipid accumulation in BrM) results in Early sub-RPE CNV. 4) Normal adhesion of RBaM to BrM, but reduced apical RPE-POS or epithelial RPE-RPE adhesion (e.g. due to inflammation) results in Early sub-retinal CNV. 5) Simultaneous reduction in RPE-RPE epithelial binding and RPE-BrM adhesion results in either sub-RPE or sub-retinal CNV which often progresses to combined pattern CNV. These findings suggest that defects in adhesion dominate CNV initiation and progression.

Integration of experimental and computational approaches to sprouting angiogenesis

PURPOSE OF REVIEW

We summarize recent experimental and computational studies that investigate molecular and cellular mechanisms of sprouting angiogenesis. We discuss how experimental tools have unveiled new opportunities for computational modeling by providing detailed phenomenological descriptions and conceptual models of cell-level behaviors underpinned by high-quality molecular data. Using recent examples, we show how new understanding results from bridging computational and experimental approaches.

RECENT FINDINGS

Experimental data extends beyond the tip cell vs. stalk cell paradigm, and involves numerous molecular inputs such as vascular endothelial growth factor and Notch. This data is being used to generate and validate computational models, which can then be used to predict the results of hypothetical experiments that are difficult to perform in the laboratory, and to generate new hypotheses that account for system-wide interactions. As a result of this integration, descriptions of critical gradients of growth factor-receptor complexes have been generated, and new modulators of cell behavior have been described.

SUMMARY

We suggest that the recent emphasis on the different stages of sprouting angiogenesis, and integration of experimental and computational approaches, should provide a way to manage the complexity of this process and help identify new regulatory paradigms and therapeutic targets.

On the evolution of morphogenetic models: mechano-chemical interactions and an integrated view of cell differentiation, growth, pattern formation and morphogenesis

In the 1950s, embryology was conceptualized as four relatively independent problems: cell differentiation, growth, pattern formation and morphogenesis. The mechanisms underlying the first three traditionally have been viewed as being chemical in nature, whereas those underlying morphogenesis have usually been discussed in terms of mechanics. Often, morphogenesis and its mechanical processes have been regarded as subordinate to chemical ones. However, a growing body of evidence indicates that the biomechanics of cells and tissues affect in striking ways those phenomena often thought of as mainly under the control of cell-cell signalling. This accumulation of data has led to a revival of the mechano-transduction concept in particular, and of complexity in general, causing us now to consider whether we should retain the traditional conceptualization of development. The researchers' semantic preferences for the terms 'patterning', 'pattern formation' or 'morphogenesis' can be used to describe three main 'schools of thought' which emerged in the late 1970s. In the 'molecular school', the term patterning is deeply tied to the positional information concept. In the 'chemical school', the term 'pattern formation' regularly implies reaction-diffusion models. In the 'mechanical school', the term 'morphogenesis' is more frequently used in relation to mechanical instabilities. Major differences among these three schools pertain to the concept of self-organization, and models can be classified as morphostatic or morphodynamic. Various examples illustrate the distorted picture that arises from the distinction among differentiation, growth, pattern formation and morphogenesis, based on the idea that the underlying mechanisms are respectively chemical or mechanical. Emerging quantitative approaches integrate the concepts and methods of complex sciences and emphasize the interplay between hierarchical levels of organization via mechano-chemical interactions. They draw upon recent improvements in mathematical and numerical morphogenetic models and upon considerable progress in collecting new quantitative data. This review highlights a variety of such models, which exhibit important advances, such as hybrid, stochastic and multiscale simulations.

Invasion from a cell aggregate--the roles of active cell motion and mechanical equilibrium

Cell invasion from an aggregate into a surrounding extracellular matrix (ECM) is an important process during development disease, e.g., vascular network assembly or tumor progression. To describe the behavior emerging from autonomous cell motility, cell-cell adhesion and contact guidance by ECM filaments, we propose a suitably modified cellular Potts model. We consider an active cell motility process in which internal polarity is governed by a positive feedback from cell displacements, a mechanism that can result in highly persistent motion when constrained by an oriented ECM structure. The model allows us to explore the interplay between haptotaxis, matrix degradation and active cell movement. We show that for certain conditions the cells are able to both invade the ECM and follow the ECM tracks. Furthermore, we argue that enforcing mechanical equilibrium within a bulk cell mass is of key importance in multicellular simulations.

Multi-scale modeling of tissues using CompuCell3D

The study of how cells interact to produce tissue development, homeostasis, or diseases was, until recently, almost purely experimental. Now, multi-cell computer simulation methods, ranging from relatively simple cellular automata to complex immersed-boundary and finite-element mechanistic models, allow in silico study of multi-cell phenomena at the tissue scale based on biologically observed cell behaviors and interactions such as movement, adhesion, growth, death, mitosis, secretion of chemicals, chemotaxis, etc. This tutorial introduces the lattice-based Glazier-Graner-Hogeweg (GGH) Monte Carlo multi-cell modeling and the open-source GGH-based CompuCell3D simulation environment that allows rapid and intuitive modeling and simulation of cellular and multi-cellular behaviors in the context of tissue formation and subsequent dynamics. We also present a walkthrough of four biological models and their associated simulations that demonstrate the capabilities of the GGH and CompuCell3D.

Multiobjective optimization based-approach for discovering novel cancer therapies

Solid tumors must recruit new blood vessels for growth and maintenance. Discovering drugs that block tumor-induced development of new blood vessels (angiogenesis) is an important approach in cancer treatment. The complexity of angiogenesis presents both challenges and opportunities for cancer therapies. Intuitive approaches, such as blocking VegF activity, have yielded important therapies. But there maybe opportunities to alter nonintuitive targets either alone or in combination. This paper describes the development of a high-fidelity simulation of angiogenesis and uses this as the basis for a parallel search-based approach for the discovery of novel potential cancer treatments that inhibit blood vessel growth. Discovering new therapies is viewed as a multiobjective combinatorial optimization over two competing objectives: minimizing the estimated cost of practically developing the intervention while minimizing the simulated oxygen provided to the tumor by angiogenesis. Results show the effectiveness of the search process by finding interventions that are currently in use, and more interestingly, discovering potential new approaches that are nonintuitive yet effective.

Vascular Network Formation in Expanding versus Static Tissues: Embryos and Tumors

In this perspectives article, we review scientific literature regarding de novo formation of vascular networks within tissues undergoing a significant degree of motion. Next, we contrast dynamic pattern formation in embryos to the vascularization of relatively static tissues, such as the retina. We argue that formation of primary polygonal vascular networks is an emergent process, which is regulated by biophysical mechanisms. Dynamic empirical data, derived from quail embryos, show that vascular beds readily form within a moving extracellular matrix (ECM) microenvironment-which we analogize to the de novo vascularization of small rapidly growing tumors. Our perspective is that the biophysical rules, which govern cell motion during vasculogenesis, may hold important clues to understanding how the first vessels form in certain malignancies.

Proteolytic Activity Attenuates the Response of Endothelial Cells to Fluid Shear Stress

Recent evidence indicates that several experimental pathophysiological conditions are associated with elevated protease activity in plasma, which impacts endothelial function. We hypothesize that extracellular structures bound to the endothelial cell (EC) membrane may be degraded by proteolytic activity and cause the cells to respond abnormally to physiological shear stress (12 dyn/cm(2)). To test this hypothesis, cultured bovine aortic endothelial cells (BAECs) were exposed to low levels of a serine protease, trypsin. Extracellular mechanosensor densities of the glycocalyx and vascular endothelial growth factor receptor 2 (VEGFR-2) were determined. Metabolic dysfunction was tested by examining insulin receptor and glucose uptake levels. Protease treatment impaired the cells' ability to align in the direction of fluid flow after 12 hours of shear stress; however, cells realigned after an additional 12 hours of shear stress with protease inhibition. Proteases caused reduction in the densities of glycocalyx, VEGFR-2, and insulin receptor in static and shear conditions, except for static VEGFR-2 cells. Under static conditions, protease-treated endothelial cells had reduced glucose uptake compared to untreated controls. Under shear, however, glucose uptake for protease-treated BAECs was greater than untreated controls. In conclusion, protease activity in plasma alters the exofacial membrane components of ECs and may interfere with mechanotransduction.