Toward a Multi-Scale Computational Model of Arterial Adaptation in Hypertension: Verification of a Multi-Cell Agent Based Model

Agent-based models (ABMs) represent a novel approach to study and simulate complex mechano chemo-biological responses at the cellular level. Such models have been used to simulate a variety of emergent responses in the vasculature, including angiogenesis and vasculogenesis. Although not used previously to study large vessel adaptations, we submit that ABMs will prove equally useful in such studies when combined with well-established continuum models to form multi-scale models of tissue-level phenomena. In order to couple agent-based and continuum models, however, there is a need to ensure that each model faithfully represents the best data available at the relevant scale and that there is consistency between models under baseline conditions. Toward this end, we describe the development and verification of an ABM of endothelial and smooth muscle cell responses to mechanical stimuli in a large artery. A refined rule-set is proposed based on a broad literature search, a new scoring system for assigning confidence in the rules, and a parameter sensitivity study. To illustrate the utility of these new methods for rule selection, as well as the consistency achieved with continuum-level models, we simulate the behavior of a mouse aorta during homeostasis and in response to both transient and sustained increases in pressure. The simulated responses depend on the altered cellular production of seven key mitogenic, synthetic, and proteolytic biomolecules, which in turn control the turnover of intramural cells and extracellular matrix. These events are responsible for gross changes in vessel wall morphology. This new ABM is shown to be appropriately stable under homeostatic conditions, insensitive to transient elevations in blood pressure, and responsive to increased intramural wall stress in hypertension

Integration of experimental and computational approaches to sprouting angiogenesis

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

Multiscale computational analysis of Xenopus laevis morphogenesis reveals key insights of systems-level behavior

<p>Abstract</p> <p>Background</p> <p>Tissue morphogenesis is a complex process whereby tissue structures self-assemble by the aggregate behaviors of independently acting cells responding to both intracellular and extracellular cues in their environment. During embryonic development, morphogenesis is particularly important for organizing cells into tissues, and although key regulatory events of this process are well studied in isolation, a number of important systems-level questions remain unanswered. This is due, in part, to a lack of integrative tools that enable the coupling of biological phenomena across spatial and temporal scales. Here, we present a new computational framework that integrates intracellular signaling information with multi-cell behaviors in the context of a spatially heterogeneous tissue environment.</p> <p>Results</p> <p>We have developed a computational simulation of mesendoderm migration in the <it>Xenopus laevis </it>explant model, which is a well studied biological model of tissue morphogenesis that recapitulates many features of this process during development in humans. The simulation couples, via a JAVA interface, an ordinary differential equation-based mass action kinetics model to compute intracellular Wnt/β-catenin signaling with an agent-based model of mesendoderm migration across a fibronectin extracellular matrix substrate. The emergent cell behaviors in the simulation suggest the following properties of the system: maintaining the integrity of cell-to-cell contact signals is necessary for preventing fractionation of cells as they move, contact with the Fn substrate and the existence of a Fn gradient provides an extracellular feedback loop that governs migration speed, the incorporation of polarity signals is required for cells to migrate in the same direction, and a delicate balance of integrin and cadherin interactions is needed to reproduce experimentally observed migratory behaviors.</p> <p>Conclusion</p> <p>Our computational framework couples two different spatial scales in biology: intracellular with multicellular. In our simulation, events at one scale have quantitative and dynamic impact on events at the other scale. This integration enables the testing and identification of key systems-level hypotheses regarding how signaling proteins affect overall tissue-level behavior during morphogenesis in an experimentally verifiable system. Applications of this approach extend to the study of tissue patterning processes that occur during adulthood and disease, such as tumorgenesis and atherogenesis.</p

Flt-1 (VEGFR-1) coordinates discrete stages of blood vessel formation

In developing blood vessel networks, the overall level of vessel branching often correlates with angiogenic sprout initiations, but in some pathological situations, increased sprout initiations paradoxically lead to reduced vessel branching and impaired vascular function. We examine the hypothesis that defects in the discrete stages of angiogenesis can uniquely contribute to vessel branching outcomes

Polarized localization of phosphatidylserine in the endothelium regulates Kir2.1

Lipid regulation of ion channels is largely explored using in silico modeling with minimal experimentation in intact tissue; thus, the functional consequences of these predicted lipid-channel interactions within native cellular environments remain elusive. The goal of this study is to investigate how lipid regulation of endothelial Kir2.1 - an inwardly rectifying potassium channel that regulates membrane hyperpolarization - contributes to vasodilation in resistance arteries. First, we show that phosphatidylserine (PS) localizes to a specific subpopulation of myoendothelial junctions (MEJs), crucial signaling microdomains that regulate vasodilation in resistance arteries, and in silico data have implied that PS may compete with phosphatidylinositol 4,5-bisphosphate (PIP2) binding on Kir2.1. We found that Kir2.1-MEJs also contained PS, possibly indicating an interaction where PS regulates Kir2.1. Electrophysiology experiments on HEK cells demonstrate that PS blocks PIP2 activation of Kir2.1 and that addition of exogenous PS blocks PIP2-mediated Kir2.1 vasodilation in resistance arteries. Using a mouse model lacking canonical MEJs in resistance arteries (Elnfl/fl/Cdh5-Cre), PS localization in endothelium was disrupted and PIP2 activation of Kir2.1 was significantly increased. Taken together, our data suggest that PS enrichment to MEJs inhibits PIP2-mediated activation of Kir2.1 to tightly regulate changes in arterial diameter, and they demonstrate that the intracellular lipid localization within the endothelium is an important determinant of vascular function

Rapid Analysis of Vessel Elements (RAVE): A Tool for Studying Physiologic, Pathologic and Tumor Angiogenesis

Quantification of microvascular network structure is important in a myriad of emerging research fields including microvessel remodeling in response to ischemia and drug therapy, tumor angiogenesis, and retinopathy. To mitigate analyst-specific variation in measurements and to ensure that measurements represent actual changes in vessel network structure and morphology, a reliable and automatic tool for quantifying microvascular network architecture is needed. Moreover, an analysis tool capable of acquiring and processing large data sets will facilitate advanced computational analysis and simulation of microvascular growth and remodeling processes and enable more high throughput discovery. To this end, we have produced an automatic and rapid vessel detection and quantification system using a MATLAB graphical user interface (GUI) that vastly reduces time spent on analysis and greatly increases repeatability. Analysis yields numerical measures of vessel volume fraction, vessel length density, fractal dimension (a measure of tortuosity), and radii of murine vascular networks. Because our GUI is open sourced to all, it can be easily modified to measure parameters such as percent coverage of non-endothelial cells, number of loops in a vascular bed, amount of perfusion and two-dimensional branch angle. Importantly, the GUI is compatible with standard fluorescent staining and imaging protocols, but also has utility analyzing brightfield vascular images, obtained, for example, in dorsal skinfold chambers. A manually measured image can be typically completed in 20 minutes to 1 hour. In stark comparison, using our GUI, image analysis time is reduced to around 1 minute. This drastic reduction in analysis time coupled with increased repeatability makes this tool valuable for all vessel research especially those requiring rapid and reproducible results, such as anti-angiogenic drug screening

Agent-Based Model of Therapeutic Adipose-Derived Stromal Cell Trafficking during Ischemia Predicts Ability To Roll on P-Selectin

Intravenous delivery of human adipose-derived stromal cells (hASCs) is a promising option for the treatment of ischemia. After delivery, hASCs that reside and persist in the injured extravascular space have been shown to aid recovery of tissue perfusion and function, although low rates of incorporation currently limit the safety and efficacy of these therapies. We submit that a better understanding of the trafficking of therapeutic hASCs through the microcirculation is needed to address this and that selective control over their homing (organ- and injury-specific) may be possible by targeting bottlenecks in the homing process. This process, however, is incredibly complex, which merited the use of computational techniques to speed the rate of discovery. We developed a multicell agent-based model (ABM) of hASC trafficking during acute skeletal muscle ischemia, based on over 150 literature-based rules instituted in Netlogo and MatLab software programs. In silico, trafficking phenomena within cell populations emerged as a result of the dynamic interactions between adhesion molecule expression, chemokine secretion, integrin affinity states, hemodynamics and microvascular network architectures. As verification, the model reasonably reproduced key aspects of ischemia and trafficking behavior including increases in wall shear stress, upregulation of key cellular adhesion molecules expressed on injured endothelium, increased secretion of inflammatory chemokines and cytokines, quantified levels of monocyte extravasation in selectin knockouts, and circulating monocyte rolling distances. Successful ABM verification prompted us to conduct a series of systematic knockouts in silico aimed at identifying the most critical parameters mediating hASC trafficking. Simulations predicted the necessity of an unknown selectin-binding molecule to achieve hASC extravasation, in addition to any rolling behavior mediated by hASC surface expression of CD15s, CD34, CD62e, CD62p, or CD65. In vitro experiments confirmed this prediction; a subpopulation of hASCs slowly rolled on immobilized P-selectin at speeds as low as 2 µm/s. Thus, our work led to a fundamentally new understanding of hASC biology, which may have important therapeutic implications