Multicellular simulation predicts microvascular patterning and in silico tissue assembly

Multiscale computational model predicts how environmental changes and drug treatments affect microvascular remodeling in fibrotic disease

Investigating the molecular, cellular, and tissue-level changes caused by disease, and the effects of pharmacological treatments across these biological scales, necessitates the use of multiscale computational modeling in combination with experimentation. Many diseases dynamically alter the tissue microenvironment in ways that trigger microvascular network remodeling, which leads to the expansion or regression of microvessel networks. When microvessels undergo remodeling in idiopathic pulmonary fibrosis (IPF), functional gas exchange is impaired due to loss of alveolar structures and lung function declines. Here, we integrated a multiscale computational model with independent experiments to investigate how combinations of biomechanical and biochemical cues in IPF alter cell fate decisions leading to microvascular remodeling. Our computational model predicted that extracellular matrix (ECM) stiffening reduced microvessel area, which was accompanied by physical uncoupling of endothelial cell (ECs) and pericytes, the cells that comprise microvessels. Nintedanib, an FDA-approved drug for treating IPF, was predicted to further potentiate microvessel regression by decreasing the percentage of quiescent pericytes while increasing the percentage of pericytes undergoing pericyte-myofibroblast transition (PMT) in high ECM stiffnesses. Importantly, the model suggested that YAP/TAZ inhibition may overcome the deleterious effects of nintedanib by promoting EC-pericyte coupling and maintaining microvessel homeostasis. Overall, our combination of computational and experimental modeling can explain how cell decisions affect tissue changes during disease and in response to treatments.

Contributions of cell behavior to geometric order in embryonic cartilage

During early development, cartilage provides shape and stability to the embryo while serving as a precursor for the skeleton. Correct formation of embryonic cartilage is hence essential for healthy development. In vertebrate cranial cartilage, it has been observed that a flat and laterally extended macroscopic geometry is linked to regular microscopic structure consisting of tightly packed, short, transversal clonar columns. However, it remains an ongoing challenge to identify how individual cells coordinate to successfully shape the tissue, and more precisely which mechanical interactions and cell behaviors contribute to the generation and maintenance of this columnar cartilage geometry during embryogenesis. Here, we apply a three-dimensional cell-based computational model to investigate mechanical principles contributing to column formation. The model accounts for clonal expansion, anisotropic proliferation and the geometrical arrangement of progenitor cells in space. We confirm that oriented cell divisions and repulsive mechanical interactions between cells are key drivers of column formation. In addition, the model suggests that column formation benefits from the spatial gaps created by the extracellular matrix in the initial configuration, and that column maintenance is facilitated by sequential proliferative phases. Our model thus correctly predicts the dependence of local order on division orientation and tissue thickness. The present study presents the first cell-based simulations of cell mechanics during cranial cartilage formation and we anticipate that it will be useful in future studies on the formation and growth of other cartilage geometries.

Angiogenic Microvascular Wall Shear Stress Patterns Revealed Through Three-dimensional Red Blood Cell Resolved Modeling

The wall shear stress (WSS) exerted by blood flowing through microvascular capillaries is an established driver of new blood vessel growth, or angiogenesis. Such adaptations are central to many physiological processes in both health and disease, yet three-dimensional (3D) WSS characteristics in real angiogenic microvascular networks are largely unknown. This marks a major knowledge gap because angiogenesis, naturally, is a 3D process. To advance current understanding, we model 3D red blood cells (RBCs) flowing through rat angiogenic microvascular networks using state-of-the-art simulation. The high-resolution fluid dynamics reveal 3D WSS patterns occurring at sub-endothelial cell (EC) scales that derive from distinct angiogenic morphologies, including microvascular loops and vessel tortuosity. We identify the existence of WSS hot and cold spots caused by angiogenic surface shapes and RBCs, and notably enhancement of low WSS regions by RBCs. Spatiotemporal characteristics further reveal how fluctuations follow timescales of RBC "footprints." Altogether, this work provides a new conceptual framework for understanding how shear stress might regulate EC dynamics in vivo.

Mathematical models of developmental vascular remodelling: A review

Over the past 40 years, there has been a strong focus on the development of mathematical models of angiogenesis, while developmental remodelling has received little such attention from the mathematical community. Sprouting angiogenesis can be seen as a very crude way of laying out a primitive vessel network (the raw material), while remodelling (understood as pruning of redundant vessels, diameter control, and the establishment of vessel identity and hierarchy) is the key to turning that primitive network into a functional network. This multiscale problem is of prime importance in the development of a functional vasculature. In addition, defective remodelling (either during developmental remodelling or due to a reactivation of the remodelling programme caused by an injury) is associated with a significant number of diseases. In this review, we discuss existing mathematical models of developmental remodelling and explore the important contributions that these models have made to the field of vascular development. These mathematical models are effectively used to investigate and predict vascular development and are able to reproduce experimentally observable results. Moreover, these models provide a useful means of hypothesis generation and can explain the underlying mechanisms driving the observed structural and functional network development. However, developmental vascular remodelling is still a relatively new area in mathematical biology, and many biological questions remain unanswered. In this review, we present the existing modelling paradigms and define the key challenges for the field.

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.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Simulation of the process of angiogenesis: Quantification and assessment of vascular patterning in the chicken chorioallantoic membrane

Angiogenesis, the formation of new blood vessels from pre-existing ones, begins during embryonic development and continues throughout life. Sprouting angiogenesis is a well-defined process, being mainly influenced by vascular endothelial growth factor (VEGF). In this study, we propose a meshless-based model capable of mimicking the angiogenic response to several VEGF concentrations. In this model, endothelial cells migrate according to a diffusion-reaction equation, following the VEGF gradient concentration. The chick chorioallantoic membrane (CAM) assay was used to model the branching process and to validate the obtained numerical results. To analyse the angiogenic response, the total vessel number and the total vessel length presented in the CAM images and in the simulations for all the VEGF concentrations tested were quantified. In both the CAM assay and simulation, the treatments with VEGF increased the total vessel number and the total vessel length. The obtained quantitative results were very similar between the two methodologies used. The proposed model accurately simulates the capillary network pattern concerning its structure and morphology, for the lowest VEGF concentration tested. For the highest VEGF concentration, the capillary network predicted by the model was less accurate when compared to the one presented in the CAM assay but this may be explained by changes in blood vessel width at higher VEGF concentrations. This remains to be tested.

Agent-based model provides insight into the mechanisms behind failed regeneration following volumetric muscle loss injury

Skeletal muscle possesses a remarkable capacity for repair and regeneration following a variety of injuries. When successful, this highly orchestrated regenerative process requires the contribution of several muscle resident cell populations including satellite stem cells (SSCs), fibroblasts, macrophages and vascular cells. However, volumetric muscle loss injuries (VML) involve simultaneous destruction of multiple tissue components (e.g., as a result of battlefield injuries or vehicular accidents) and are so extensive that they exceed the intrinsic capability for scarless wound healing and result in permanent cosmetic and functional deficits. In this scenario, the regenerative process fails and is dominated by an unproductive inflammatory response and accompanying fibrosis. The failure of current regenerative therapeutics to completely restore functional muscle tissue is not surprising considering the incomplete understanding of the cellular mechanisms that drive the regeneration response in the setting of VML injury. To begin to address this profound knowledge gap, we developed an agent-based model to predict the tissue remodeling response following surgical creation of a VML injury. Once the model was able to recapitulate key aspects of the tissue remodeling response in the absence of repair, we validated the model by simulating the tissue remodeling response to VML injury following implantation of either a decellularized extracellular matrix scaffold or a minced muscle graft. The model suggested that the SSC microenvironment and absence of pro-differentiation SSC signals were the most important aspects of failed muscle regeneration in VML injuries. The major implication of this work is that agent-based models may provide a much-needed predictive tool to optimize the design of new therapies, and thereby, accelerate the clinical translation of regenerative therapeutics for VML injuries.

An in-silico Investigation Into the Role of Strain and Structure on Vascular Smooth Muscle Cell Growth

The orientation of vascular cells can greatly influence the in vivo mechanical properties and functionality of soft vascular tissues. How cell orientation mediates the growth response of cells is of critical importance in understanding the response of soft tissues to mechanical stimuli or injury. To date, considerable evidence has shown that cells align with structural cues such as collagen fibers. However, in the presence of uniaxial cyclic strain on unstructured substrates, cells generally align themselves perpendicularly to the mechanical stimulus, such as strain, a phenomenon known as “strain avoidance.” The cellular response to this interplay between structural cues and a mechanical stimulus is poorly understood. A recent in vitro experimental study in our lab has investigated both the individual and collective response of rat aortic smooth muscle cells (RASMC) to structural (collagenous aligned constructs) and mechanical (cyclic strain) cues. In this study, a 2D agent-based model (ABM) is developed to simulate the collective response of RASMC to varying amplitudes of cyclic strain (0–10%, 2–8%, 4–6%) when seeded on unstructured (PDMS) and structured (decellularized collagenous tissue) constructs. An ABM is presented that is fit to the experimental outcomes in terms of cellular alignment and cell growth on PDMS substrates, under cyclic strain amplitudes of (4–6%, 2–8%, 0–10%) at 24 and 72 h timepoints. Furthermore, the ABM can predict RASMC alignment and change in cell number on a structured construct at a cyclic strain amplitude of 0–10% after 10 days. The ABM suggests that strain avoidance behavior observed in cells is dominated by selective cell proliferation and apoptosis at these early time points, as opposed to cell re-orientation, i.e., cells perpendicular to the strain increase their rate of proliferation, whilst the rate of apoptosis simultaneously increases in cells parallel to the strain direction. The development of in-silico modeling platforms, such as that presented here, allow for further understanding of the response of cells to changes in their mechanical environment. Such models offer an efficient and robust means to design and optimize the compliance and topological structure of implantable devices and could be used to aid the design of next-generation vascular grafts and stents.

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

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

Sprouting Angiogenesis: A Numerical Approach with Experimental Validation

A functional vascular network is essential to the correct wound healing. In sprouting angiogenesis, vascular endothelial growth factor (VEGF) regulates the formation of new capillaries from pre-existing vessels. This is a very complex process and mathematical formulation permits to study angiogenesis using less time-consuming, reproducible and cheaper methodologies. This study aimed to mimic the chemoattractant effect of VEGF in stimulating sprouting angiogenesis. We developed a numerical model in which endothelial cells migrate according to a diffusion-reaction equation for VEGF. A chick chorioallantoic membrane (CAM) bioassay was used to obtain some important parameters to implement in the model and also to validate the numerical results. We verified that endothelial cells migrate following the highest VEGF concentration. We compared the parameters-total branching number, total vessel length and branching angle-that were obtained in the in silico and the in vivo methodologies and similar results were achieved (p-value smaller than 0.5; n = 6). For the difference between the total capillary volume fractions assessed using both methodologies values smaller than 15% were obtained. In this study we simulated, for the first time, the capillary network obtained during the CAM assay with a realistic morphology and structure.

A particle-based model for endothelial cell migration under flow conditions

Endothelial cells (ECs) play a major role in the healing process following angioplasty to inhibit excessive neointima. This makes the process of EC healing after injury, in particular EC migration in a stented vessel, important for recovery of normal vessel function. In that context, we present a novel particle-based model of EC migration and validate it against in vitro experimental data. We have developed a particle-based model of EC migration under flow conditions in an in vitro vessel with obstacles. Cell movement in the model is a combination of random walks and directed movement along the local flow velocity vector. For model calibration, a set of experimental data for cell migration in a similarly shaped channel has been used. We have calibrated the model for a baseline case of a channel with no obstacles and then applied it to the case of a channel with ridges on the bottom surface, representative of stent strut geometry. We were able to closely reproduce the cell migration speed and angular distribution of their movement relative to the flow direction reported in vitro. The model also reproduces qualitative aspects of EC migration, such as entrapment of cells downstream from the flow-disturbing ridge. The model has the potential, after more extensive in vitro validation, to study the effect of variation in strut spacing and shape, through modification of the local flow, on EC migration. The results of this study support the hypothesis that EC migration is strongly affected by the direction and magnitude of local wall shear stress.

A fully coupled framework for in silico investigation of in-stent restenosis

Finite element analysis (FEA) can be implemented along with Agent-based model (ABM) to investigate the biomechanical and mechanobiological mechanisms of pathophysiological processes. However, traditional ABM-FEA approaches are often partially coupled and lack the feedback responses from biological analysis. To overcome this problem, a fully coupled ABM-FEA framework is developed in this paper by linking the macro-scale and cell-scale modules bi-directionally. Numerical studies of the in-stent restenosis process are conducted using the proposed approach and comparisons are made between the two types of frameworks. A reduction in lumen loss rate, which is possibly caused by the time-varying stresses, is observed in the fully coupled simulations. The re-endothelialisation process is also simulated under different frameworks and the simulation results show strong inhibition of endothelial cells to vascular restenosis. The proposed method is proved to be effective to explain the biomechanical-mechanobiological coupling characteristics of the restenosis problem and can be utilized for stent design and optimization.

Modelling skin wound healing angiogenesis: A review

The occurrence of wounds is a main health concern in Western society due to their high frequency and treatment cost. During wound healing, the formation of a functional blood vessel network through angiogenesis is an essential process. Angiogenesis allows the reestablishment of the normal blood flow, the sufficient exchange of oxygen and nutrients and the removal of metabolic waste, necessary for cell proliferation and viability. Mathematical and computational models provide new tools to improve the healing process. In fact, over the last thirty years, in silico models have been continuously formulated to describe the effect of several biological and mechanical factors in angiogenesis during wound healing. Additionally, with different levels of complexity, these models allow coupling the human skin structure, to distinct cell types and growth factors, to study extracellular matrix composition and to understand its deformation. This paper discusses how in silico models, which are more economical and less time-consuming comparatively to laboratory methodologies, can help test new strategies to promote/optimize angiogenesis. The continuum, cell-based and hybrid mathematical models of wound healing angiogenesis are reviewed in the present paper, in order to identify possible improvements. Accordingly, the development of higher dimension models incorporating multiscale analysis at molecular, cellular and tissue level remains a challenge that future models should consider.

Modeling persistence of motion in a crowded environment: The diffusive limit of excluding velocity-jump processes

Persistence of motion is the tendency of an object to maintain motion in a direction for short time scales without necessarily being biased in any direction in the long term. One of the most appropriate mathematical tools to study this behavior is an agent-based velocity-jump process. In the absence of agent-agent interaction, the mean-field continuum limit of the agent-based model (ABM) gives rise to the well known hyperbolic telegraph equation. When agent-agent interaction is included in the ABM, a strictly advective system of partial differential equations (PDEs) can be derived at the population level. However, no diffusive limit of the ABM has been obtained from such a model. Connecting the microscopic behavior of the ABM to a diffusive macroscopic description is desirable, since it allows the exploration of a wider range of scenarios and establishes a direct connection with commonly used statistical tools of movement analysis. In order to connect the ABM at the population level to a diffusive PDE at the population level, we consider a generalization of the agent-based velocity-jump process on a two-dimensional lattice with three forms of agent interaction. This generalization allows us to take a diffusive limit and obtain a faithful population-level description. We investigate the properties of the model at both the individual and population levels and we elucidate some of the models' key characteristic features. In particular, we show an intrinsic anisotropy inherent to the models and we find evidence of a spontaneous form of aggregation at both the micro- and macroscales.

VEGFR1 promotes cell migration and proliferation through PLCγ and PI3K pathways

The ability to control vascular endothelial growth factor (VEGF) signaling offers promising therapeutic potential for vascular diseases and cancer. Despite this promise, VEGF-targeted therapies are not clinically effective for many pathologies, such as breast cancer. VEGFR1 has recently emerged as a predictive biomarker for anti-VEGF efficacy, implying a functional VEGFR1 role beyond its classically defined decoy receptor status. Here we introduce a computational approach that accurately predicts cellular responses elicited via VEGFR1 signaling. Aligned with our model prediction, we show empirically that VEGFR1 promotes macrophage migration through PLCγ and PI3K pathways and promotes macrophage proliferation through a PLCγ pathway. These results provide new insight into the basic function of VEGFR1 signaling while offering a computational platform to quantify signaling of any receptor.

Mechanobiological model of arterial growth and remodeling

A coupled agent-based model (ABM) and finite element analysis (FEA) computational framework is developed to study the interplay of bio-chemo-mechanical factors in blood vessels and their role in maintaining homeostasis. The agent-based model implements the power of REPAST Simphony libraries and adapts its environment for biological simulations. Coupling a continuum-level model (FEA) to a cellular-level model (ABM) has enabled this computational framework to capture the response of blood vessels to increased or decreased levels of growth factors, proteases and other signaling molecules (on the micro scale) as well as altered blood pressure. Performance of the model is assessed by simulating porcine left anterior descending artery under normotensive conditions and transient increases in blood pressure and by analyzing sensitivity of the model to variations in the rule parameters of the ABM. These simulations proved that the model is stable under normotensive conditions and can recover from transient increases in blood pressure. Sensitivity studies revealed that the model is most sensitive to variations in the concentration of growth factors that affect cellular proliferation and regulate extracellular matrix composition (mainly collagen).

Computational Model-Based Analysis of Strategies to Enhance Scaffold Vascularization

Stable and extensive blood vessel networks are required for cell function and survival in engineered tissues. A number of different strategies are currently being investigated to enhance biomaterial vascularization with screening primarily through extensive in vitro and in vivo experiments. In this article, we describe an agent-based model (ABM) developed to evaluate various strategies in silico, including design of optimal biomaterial structure, delivery of angiogenic factors, and application of prevascularized biomaterials. The model predictions are evaluated using experimental data. The ABM developed provides insight into different strategies currently applied for scaffold vascularization and will enable researchers to rapidly screen new hypotheses and explore alternative strategies for enhancing vascularization.

Computer Simulations of the Tumor Vasculature: Applications to Interstitial Fluid Flow, Drug Delivery, and Oxygen Supply

Tumor vasculature, the blood vessel network supplying a growing tumor with nutrients such as oxygen or glucose, is in many respects different from the hierarchically organized arterio-venous blood vessel network in normal tissues. Angiogenesis (the formation of new blood vessels), vessel cooption (the integration of existing blood vessels into the tumor vasculature), and vessel regression remodel the healthy vascular network into a tumor-specific vasculature. Integrative models, based on detailed experimental data and physical laws, implement, in silico, the complex interplay of molecular pathways, cell proliferation, migration, and death, tissue microenvironment, mechanical and hydrodynamic forces, and the fine structure of the host tissue vasculature. With the help of computer simulations high-precision information about blood flow patterns, interstitial fluid flow, drug distribution, oxygen and nutrient distribution can be obtained and a plethora of therapeutic protocols can be tested before clinical trials. This chapter provides an overview over the current status of computer simulations of vascular remodeling during tumor growth including interstitial fluid flow, drug delivery, and oxygen supply within the tumor. The model predictions are compared with experimental and clinical data and a number of longstanding physiological paradigms about tumor vasculature and intratumoral solute transport are critically scrutinized.

Agent-based modeling of osteogenic differentiation of mesenchymal stem cells in porous biomaterials

Mesenchymal stem cells (MSC) have shown promise in tissue engineering applications due to their potential for differentiating into mesenchymal tissues such as osteocytes, chondrocytes, and adipocytes and releasing proteins to promote tissue regeneration. One application involves seeding MSCs in biomaterial scaffolds to promote osteogenesis in the repair of bone defects following implantation. However, predicting in vivo survival and differentiation of MSCs in biomaterials is challenging. Rapid and stable vascularization of scaffolds is required to supply nutrients and oxygen that MSCs need to survive as well as to go through osteogenic differentiation. The objective of this study is to develop an agent-based model and simulator that can be used to investigate the effects of using gradient growth factors on survival and differentiation of MSCs seeded in scaffolds. An agent-based model is developed to simulate the MSC behavior. The effect of vascular endothelial growth factor (VEGF) and bone morphogenic protein-2 (BMP-2) on both survival and osteogenic differentiation is studied. Results showed that the survival ratio of MSCs can be enhanced by increasing VEGF concentration. BMP-2 caused a slight increase on survival ratio. Osteogenesis strongly depends on the VEGF concentration as well because of its effect on vascularization. BMP-2 increased the osteogenic differentiation of MSCs.

PROKARYO: an illustrative and interactive computational model of the lactose operon in the bacterium Escherichia coli

Background

We are creating software for agent-based simulation and visualization of bio-molecular processes in bacterial and eukaryotic cells. As a first example, we have built a 3-dimensional, interactive computer model of an Escherichia coli bacterium and its associated biomolecular processes. Our illustrative model focuses on the gene regulatory processes that control the expression of genes involved in the lactose operon. Prokaryo, our agent-based cell simulator, incorporates cellular structures, such as plasma membranes and cytoplasm, as well as elements of the molecular machinery, including RNA polymerase, messenger RNA, lactose permease, and ribosomes.

Results

The dynamics of cellular ’agents’ are defined by their rules of interaction, implemented as finite state machines. The agents are embedded within a 3-dimensional virtual environment with simulated physical and electrochemical properties. The hybrid model is driven by a combination of (1) mathematical equations (DEQs) to capture higher-scale phenomena and (2) agent-based rules to implement localized interactions among a small number of molecular elements. Consequently, our model is able to capture phenomena across multiple spatial scales, from changing concentration gradients to one-on-one molecular interactions.

We use the classic gene regulatory mechanism of the lactose operon to demonstrate our model’s resolution, visual presentation, and real-time interactivity. Our agent-based model expands on a sophisticated mathematical E. coli metabolism model, through which we highlight our model’s scientific validity.

Conclusion

We believe that through illustration and interactive exploratory learning a model system like Prokaryo can enhance the general understanding and perception of biomolecular processes. Our agent-DEQ hybrid modeling approach can also be of value to conceptualize, illustrate, and—eventually—validate cell experiments in the wet lab.

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.

Applications of computational models to better understand microvascular remodelling: a focus on biomechanical integration across scales

Microvascular network remodelling is a common denominator for multiple pathologies and involves both angiogenesis, defined as the sprouting of new capillaries, and network patterning associated with the organization and connectivity of existing vessels. Much of what we know about microvascular remodelling at the network, cellular and molecular scales has been derived from reductionist biological experiments, yet what happens when the experiments provide incomplete (or only qualitative) information? This review will emphasize the value of applying computational approaches to advance our understanding of the underlying mechanisms and effects of microvascular remodelling. Examples of individual computational models applied to each of the scales will highlight the potential of answering specific questions that cannot be answered using typical biological experimentation alone. Looking into the future, we will also identify the needs and challenges associated with integrating computational models across scales.

Genome scale modeling in systems biology: algorithms and resources

In recent years, in silico studies and trial simulations have complemented experimental procedures. A model is a description of a system, and a system is any collection of interrelated objects; an object, moreover, is some elemental unit upon which observations can be made but whose internal structure either does not exist or is ignored. Therefore, any network analysis approach is critical for successful quantitative modeling of biological systems. This review highlights some of most popular and important modeling algorithms, tools, and emerging standards for representing, simulating and analyzing cellular networks in five sections. Also, we try to show these concepts by means of simple example and proper images and graphs. Overall, systems biology aims for a holistic description and understanding of biological processes by an integration of analytical experimental approaches along with synthetic computational models. In fact, biological networks have been developed as a platform for integrating information from high to low-throughput experiments for the analysis of biological systems. We provide an overview of all processes used in modeling and simulating biological networks in such a way that they can become easily understandable for researchers with both biological and mathematical backgrounds. Consequently, given the complexity of generated experimental data and cellular networks, it is no surprise that researchers have turned to computer simulation and the development of more theory-based approaches to augment and assist in the development of a fully quantitative understanding of cellular dynamics.

The effects of cell compressibility, motility and contact inhibition on the growth of tumor cell clusters using the Cellular Potts Model

There are numerous biological examples where genes associated with migratory ability of cells also confer the cells with an increased fitness even though these genes may not have any known effect on the cell mitosis rates. Here, we provide insight into these observations by analyzing the effects of cell migration, compression, and contact inhibition on the growth of tumor cell clusters using the Cellular Potts Model (CPM) in a monolayer geometry. This is a follow-up of a previous study (Thalhauser et al. 2010) in which a Moran-type model was used to study the interaction of cell proliferation, migratory potential and death on the emergence of invasive phenotypes. Here, we extend the study to include the effects of cell size and shape. In particular, we investigate the interplay between cell motility and compressibility within the CPM and find that the CPM predicts that increased cell motility leads to smaller cells. This is an artifact in the CPM. An analysis of the CPM reveals an explicit inverse-relationship between the cell stiffness and motility parameters. We use this relationship to compensate for motility-induced changes in cell size in the CPM so that in the corrected CPM, cell size is independent of the cell motility. We find that subject to comparable levels of compression, clusters of motile cells grow faster than clusters of less motile cells, in qualitative agreement with biological observations and our previous study. Increasing compression tends to reduce growth rates. Contact inhibition penalizes clumped cells by halting their growth and gives motile cells an even greater advantage. Finally, our model predicts cell size distributions that are consistent with those observed in clusters of neuroblastoma cells cultured in low and high density conditions.

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.

Emergence of spatial structure in the tumor microenvironment due to the Warburg effect

Drastic metabolic alterations, such as the Warburg effect, are found in most if not all types of malignant tumors. Emerging evidence shows that cancer cells benefit from these alterations, but little is known about how they affect noncancerous stromal cells within the tumor microenvironment. Here we show that cancer cells are better adapted to metabolic changes in the microenvironment, leading to the emergence of spatial structure. A clear example of tumor spatial structure is the localization of tumor-associated macrophages (TAMs), one of the most common stromal cell types found in tumors. TAMs are enriched in well-perfused areas, such as perivascular and cortical regions, where they are known to potentiate tumor growth and invasion. However, the mechanisms of TAM localization are not completely understood. Computational modeling predicts that gradients--of nutrients, gases, and metabolic by-products such as lactate--emerge due to altered cell metabolism within poorly perfused tumors, creating ischemic regions of the tumor microenvironment where TAMs struggle to survive. We tested our modeling prediction in a coculture system that mimics the tumor microenvironment. Using this experimental approach, we showed that a combination of metabolite gradients and differential sensitivity to lactic acid is sufficient for the emergence of macrophage localization patterns in vitro. This suggests that cancer metabolic changes create a microenvironment where tumor cells thrive over other cells. Understanding differences in tumor-stroma sensitivity to these alterations may open therapeutic avenues against cancer.

A hybrid model of the role of VEGF binding in endothelial cell migration and capillary formation

Vascular endothelial growth factor (VEGF) is the most studied family of soluble, secreted mediators of endothelial cell migration, survival, and proliferation. VEGF exerts its function by binding to specific tyrosine kinase receptors on the cell surface and transducing the effect through downstream signaling. In order to study the influence of VEGF binding on endothelial cell motion, we develop a hybrid model of VEGF-induced angiogenesis, based on the theory of reinforced random walks. The model includes the chemotactic response of endothelial cells to angiogenic factors bound to cell-surface receptors, rather than approximating this as a function of extracellular chemical concentrations. This allows us to capture biologically observed phenomena such as activation and polarization of endothelial cells in response to VEGF gradients across their lengths, as opposed to extracellular gradients throughout the tissue. We also propose a novel and more biologically reasonable functional form for the chemotactic sensitivity of endothelial cells, which is also governed by activated cell-surface receptors. This model is able to predict the threshold level of VEGF required to activate a cell to move in a directed fashion as well as an optimal VEGF concentration for motion. Model validation is achieved by comparison of simulation results directly with experimental data.

Computational medicine: translating models to clinical care

Because of the inherent complexity of coupled nonlinear biological systems, the development of computational models is necessary for achieving a quantitative understanding of their structure and function in health and disease. Statistical learning is applied to high-dimensional biomolecular data to create models that describe relationships between molecules and networks. Multiscale modeling links networks to cells, organs, and organ systems. Computational approaches are used to characterize anatomic shape and its variations in health and disease. In each case, the purposes of modeling are to capture all that we know about disease and to develop improved therapies tailored to the needs of individuals. We discuss advances in computational medicine, with specific examples in the fields of cancer, diabetes, cardiology, and neurology. Advances in translating these computational methods to the clinic are described, as well as challenges in applying models for improving patient health.

Systems-based approaches toward wound healing

Wound healing in the pediatric patient is of utmost clinical and social importance because hypertrophic scarring can have aesthetic and psychological sequelae, from early childhood to late adolescence. Wound healing is a well-orchestrated reparative response affecting the damaged tissue at the cellular, tissue, organ, and system scales. Although tremendous progress has been made toward understanding wound healing at the individual temporal and spatial scales, its effects across the scales remain severely understudied and poorly understood. Here, we discuss the critical need for systems-based computational modeling of wound healing across the scales, from short-term to long-term and from small to large. We illustrate the state of the art in systems modeling by means of three key signaling mechanisms: oxygen tension-regulating angiogenesis and revascularization; transforming growth factor-β (TGF-β) kinetics controlling collagen deposition; and mechanical stretch stimulating cellular mitosis and extracellular matrix (ECM) remodeling. The complex network of biochemical and biomechanical signaling mechanisms and the multiscale character of the healing process make systems modeling an integral tool in exploring personalized strategies for wound repair. A better mechanistic understanding of wound healing in the pediatric patient could open new avenues in treating children with skin disorders such as birth defects, skin cancer, wounds, and burn injuries.

Three-dimensional modeling of angiogenesis in porous biomaterial scaffolds

Vascularization of biomaterial scaffolds is essential for the successful clinical application of engineered tissues. Experimental studies are often performed to investigate the role of scaffold architecture on vascularized tissue formation. However, experiments are expensive and time-consuming and synthesis protocols often do not allow for independent investigation of specific scaffold properties. Computational models allow for rapid screening of potential material designs with control over scaffold properties that is difficult in laboratory settings. We have developed and tested a three-dimensional agent-based framework for investigating the effect of scaffold pore architecture on angiogenesis. Software agents represent endothelial cells, interacting together and with their micro-environment, leading to the invasion of blood vessels into the scaffold. A rule base, driven by experimental findings, governs the behavior of individual agents. 3D scaffold models with well-defined homogeneous and heterogeneous pore architectures were simulated to investigate the impact of various design parameters. Simulation results indicate that pores of larger size with higher interconnectivity and porosity support rapid and extensive angiogenesis. The developed framework can be used to screen biomaterial scaffold designs for optimal vascularization and investigate complex interactions among invading blood vessels and their micro-environment.

Multiscale models of breast cancer progression

Breast cancer initiation, invasion and metastasis span multiple length and time scales. Molecular events at short length scales lead to an initial tumorigenic population, which left unchecked by immune action, acts at increasingly longer length scales until eventually the cancer cells escape from the primary tumor site. This series of events is highly complex, involving multiple cell types interacting with (and shaping) the microenvironment. Multiscale mathematical models have emerged as a powerful tool to quantitatively integrate the convective-diffusion-reaction processes occurring on the systemic scale, with the molecular signaling processes occurring on the cellular and subcellular scales. In this study, we reviewed the current state of the art in cancer modeling across multiple length scales, with an emphasis on the integration of intracellular signal transduction models with pro-tumorigenic chemical and mechanical microenvironmental cues. First, we reviewed the underlying biomolecular origin of breast cancer, with a special emphasis on angiogenesis. Then, we summarized the development of tissue engineering platforms which could provide high-fidelity ex vivo experimental models to identify and validate multiscale simulations. Lastly, we reviewed top-down and bottom-up multiscale strategies that integrate subcellular networks with the microenvironment. We present models of a variety of cancers, in addition to breast cancer specific models. Taken together, we expect as the sophistication of the simulations increase, that multiscale modeling and bottom-up agent-based models in particular will become an increasingly important platform technology for basic scientific discovery, as well as the identification and validation of potentially novel therapeutic targets.

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.

Multiscale models of cell signaling

Computational models of signal transduction face challenges of scale below the resolution of a single cell. Here, we organize these challenges around three key interfaces for multiscale models of cell signaling: molecules to pathways, pathways to networks, and networks to outcomes. Each interface requires its own set of computational approaches and systems-level data, and no single approach or dataset can effectively bridge all three interfaces. This suggests that realistic "whole-cell" models of signaling will need to agglomerate different model types that span critical intracellular scales. Future multiscale models will be valuable for understanding the impact of signaling mutations or population variants that lead to cellular diseases such as cancer.

Evaluation of angiogenesis and osteogenesis

Bone regeneration has long been a major focus for tissue engineers and the importance of vascularization to the bone regeneration process has been well documented. Over the past decade, technological advances in the areas of stem cell biology, scaffold fabrication, and protein engineering have significantly enhanced our understanding of the interplay between vascularization and bone growth. This review, therefore, describes the commonly used models for investigating the complex interactions between osteoblastic cells and endothelial cells, evaluates the different tools utilized to investigate the relationship between vascularization and bone growth in vivo, and finally, summarizes possible areas of research related to therapeutic development.

Computational modeling of interacting VEGF and soluble VEGF receptor concentration gradients

Experimental data indicates that soluble vascular endothelial growth factor (VEGF) receptor 1 (sFlt-1) modulates the guidance cues provided to sprouting blood vessels by VEGF-A. To better delineate the role of sFlt-1 in VEGF signaling, we have developed an experimentally based computational model. This model describes dynamic spatial transport of VEGF, and its binding to receptors Flt-1 and Flk-1, in a mouse embryonic stem cell model of vessel morphogenesis. The model represents the local environment of a single blood vessel. Our simulations predict that blood vessel secretion of sFlt-1 and increased local sFlt-1 sequestration of VEGF results in decreased VEGF–Flk-1 levels on the sprout surface. In addition, the model predicts that sFlt-1 secretion increases the relative gradient of VEGF–Flk-1 along the sprout surface, which could alter endothelial cell perception of directionality cues. We also show that the proximity of neighboring sprouts may alter VEGF gradients, VEGF receptor binding, and the directionality of sprout growth. As sprout distances decrease, the probability that the sprouts will move in divergent directions increases. This model is a useful tool for determining how local sFlt-1 and VEGF gradients contribute to the spatial distribution of VEGF receptor binding, and can be used in conjunction with experimental data to explore how multi-cellular interactions and relationships between local growth factor gradients drive angiogenesis.

An agent-based model for the investigation of neovascularization within porous scaffolds

The ability to control blood vessel assembly in polymer scaffolds is important for clinical success in tissue engineering. A mathematical and computational representation of the relationship between scaffold properties and neovascularization may provide a better understanding of the fundamental process itself and help guide the design of new therapeutic approaches. This article proposes a multilayered, multiagent framework to model sprouting angiogenesis in porous scaffolds and examines the impact of pore structure on vessel invasion and network structure. We have defined the speed of vessel sprouting in the agent-based model based on in vivo results in the absence of a polymer scaffold. A number of cases were run to investigate the effect of scaffold pore size on angiogenesis. The simulation results indicate that the rate of scaffold vascularization increases with pore size. Pores of larger size (160-270 μm) support rapid and extensive angiogenesis throughout the scaffold. Model predictions were compared to experimental results of vascularization in porous poly(ethylene glycol) hydrogels to validate the results. This model can be used to provide insight into optimal scaffold properties that support vascularization of engineered tissues.

Emergence of self-reproduction in cooperative chemical evolution of prebiological molecules

The paper presents a model of coevolution of short peptides (P) and short oligonucleotides (N) at an early stage of chemical evolution leading to the origin of life. The model describes polymerization of both P and N types of molecules on mineral surfaces in aqueous solution at moderate temperatures. It is assumed that amino acid and nucleotide monomers were available in a prebiotic milieu, that periodic variation in environmental conditions between dry/warm and wet/cool took place and that energy sources were available for the polymerization. An artificial chemistry approach in combination with agent-based modeling was used to explore chemical evolution from an initially random mixture of monomers. It was assumed that the oligonucleotides could serve as templates for self-replication and for translation of peptide compositional sequences, and that certain peptides could serve as weak catalysts. Important features of the model are the short lengths of the peptide and oligonucleotide molecules that prevent an error catastrophe caused by copying errors and a finite diffusion rate of the molecules on a mineral surface that prevents excessive development of parasitism. The result of the simulation was the emergence of self-replicating molecular systems consisting of peptide catalysts and oligonucleotide templates. In addition, a smaller but significant number of molecules with alternative compositions also survived due to imprecise reproduction and translation of templates providing variability for further evolution. In a more general context, the model describes not only peptide-oligonucleotide molecular systems, but any molecular system containing two types of polymer molecules: one of which serves as templates and the other as catalysts.The presented coevolutionary system suggests a possible direction towards finding the origin of molecular functionality in a prebiotic environment.