Rule-Based Simulation of Multi-Cellular Biological Systems-A Review of Modeling Techniques

Fibre crosslinking drives the emergence of order in a three-dimensional dynamical network model

The extracellular-matrix (ECM) is a complex interconnected three-dimensional network that provides structural support for the cells and tissues and defines organ architecture as key for their healthy functioning. However, the intimate mechanisms by which ECM acquire their three-dimensional architecture are still largely unknown. In this paper, we study this question by means of a simple three-dimensional individual based model of interacting fibres able to spontaneously crosslink or unlink to each other and align at the crosslinks. We show that such systems are able to spontaneously generate different types of architectures. We provide a thorough analysis of the emerging structures by an exhaustive parametric analysis and the use of appropriate visualization tools and quantifiers in three dimensions. The most striking result is that the emergence of ordered structures can be fully explained by a single emerging variable: the number of links per fibre in the network. If validated on real tissues, this simple variable could become an important putative target to control and predict the structuring of biological tissues, to suggest possible new therapeutic strategies to restore tissue functions after disruption, and to help in the development of collagen-based scaffolds for tissue engineering. Moreover, the model reveals that the emergence of architecture is a spatially homogeneous process following a unique evolutionary path, and highlights the essential role of dynamical crosslinking in tissue structuring.

A feedback mechanism controls rDNA copy number evolution in yeast independently of natural selection

Ribosomal DNA (rDNA) is the genetic loci that encodes rRNA in eukaryotes. It is typically arranged as tandem repeats that vary in copy number within the same species. We have recently shown that rDNA repeats copy number in the yeast Saccharomyces cerevisiae is controlled by cell volume via a feedback circuit that senses cell volume by means of the concentration of the free upstream activator factor (UAF). The UAF strongly binds the rDNA gene promoter, but is also able to repress SIR2 deacetylase gene transcription that, in turn, represses rDNA amplification. In this way, the cells with a smaller DNA copy number than what is optimal evolve to increase that copy number until they reach a number that sequestrates free UAF and provokes SIR2 derepression that, in turn, blocks rDNA amplification. Here we propose a mathematical model to show that this evolutionary process can amplify rDNA repeats independently of the selective advantage of yeast cells having bigger or smaller rDNA copy numbers. We test several variants of this process and show that it can explain the observed experimental results independently of natural selection. These results predict that an autoregulated feedback circuit may, in some instances, drive to non Darwinian deterministic evolution for a limited time period.

BioSANS: A software package for symbolic and numeric biological simulation

Modeling biochemical systems can provide insights into behaviors that are difficult to observe or understand. It requires software, programming, and understanding of the system to build a model and study it. Softwares exist for systems biology modeling, but most support only certain types of modeling tasks. Desirable features including ease in preparing input, symbolic or analytical computation, parameter estimation, graphical user interface, and systems biology markup language (SBML) support are not seen concurrently in one software package. In this study, we developed a python-based software that supports these features, with both deterministic and stochastic propagations. The software can be used by graphical user interface, command line, or as a python import. We also developed a semi-programmable and intuitively easy topology input method for the biochemical reactions. We tested the software with semantic and stochastic SBML test cases. Tests on symbolic solution and parameter estimation were also included. The software we developed is reliable, well performing, convenient to use, and compliant with most of the SBML tests. So far it is the only systems biology software that supports symbolic, deterministic, and stochastic modeling in one package that also features parameter estimation and SBML support. This work offers a comprehensive set of tools and allows for better availability and accessibility for studying kinetics and dynamics in biochemical systems.

Data-driven learning how oncogenic gene expression locally alters heterocellular networks

Developing drugs increasingly relies on mechanistic modeling and simulation. Models that capture causal relations among genetic drivers of oncogenesis, functional plasticity, and host immunity complement wet experiments. Unfortunately, formulating such mechanistic cell-level models currently relies on hand curation, which can bias how data is interpreted or the priority of drug targets. In modeling molecular-level networks, rules and algorithms are employed to limit a priori biases in formulating mechanistic models. Here we combine digital cytometry with Bayesian network inference to generate causal models of cell-level networks linking an increase in gene expression associated with oncogenesis with alterations in stromal and immune cell subsets from bulk transcriptomic datasets. We predict how increased Cell Communication Network factor 4, a secreted matricellular protein, alters the tumor microenvironment using data from patients diagnosed with breast cancer and melanoma. Predictions are then tested using two immunocompetent mouse models for melanoma, which provide consistent experimental results.

While mechanistic models play increasing roles in immuno-oncology, hand network curation is current practice. Here the authors use a Bayesian data-driven approach to infer how expression of a secreted oncogene alters the cellular landscape within the tumor.

Multiscale Computational Modeling of Vascular Adaptation: A Systems Biology Approach Using Agent-Based Models

The widespread incidence of cardiovascular diseases and associated mortality and morbidity, along with the advent of powerful computational resources, have fostered an extensive research in computational modeling of vascular pathophysiology field and promoted in-silico models as a support for biomedical research. Given the multiscale nature of biological systems, the integration of phenomena at different spatial and temporal scales has emerged to be essential in capturing mechanobiological mechanisms underlying vascular adaptation processes. In this regard, agent-based models have demonstrated to successfully embed the systems biology principles and capture the emergent behavior of cellular systems under different pathophysiological conditions. Furthermore, through their modular structure, agent-based models are suitable to be integrated with continuum-based models within a multiscale framework that can link the molecular pathways to the cell and tissue levels. This can allow improving existing therapies and/or developing new therapeutic strategies. The present review examines the multiscale computational frameworks of vascular adaptation with an emphasis on the integration of agent-based approaches with continuum models to describe vascular pathophysiology in a systems biology perspective. The state-of-the-art highlights the current gaps and limitations in the field, thus shedding light on new areas to be explored that may become the future research focus. The inclusion of molecular intracellular pathways (e.g., genomics or proteomics) within the multiscale agent-based modeling frameworks will certainly provide a great contribution to the promising personalized medicine. Efforts will be also needed to address the challenges encountered for the verification, uncertainty quantification, calibration and validation of these multiscale frameworks.

A rigid body framework for multicellular modeling

Off-lattice models are a well-established approach in multicellular modeling, where cells are represented as points that are free to move in space. The representation of cells as point objects is useful in a wide range of settings, particularly when large populations are involved; however, a purely point-based representation is not naturally equipped to deal with objects that have length, such as cell boundaries or external membranes. Here we introduce an off-lattice modeling framework that exploits rigid body mechanics to represent objects using a collection of conjoined one-dimensional edges in a viscosity-dominated system. This framework can be used to represent cells as free moving polygons, to allow epithelial layers to smoothly interact with themselves, to model rod-shaped cells such as bacteria and to robustly represent membranes. We demonstrate that this approach offers solutions to the problems that limit the scope of current off-lattice multicellular models.

Shaping the Organ: A Biologist Guide to Quantitative Models of Plant Morphogenesis

Organ morphogenesis is the process of shape acquisition initiated with a small reservoir of undifferentiated cells. In plants, morphogenesis is a complex endeavor that comprises a large number of interacting elements, including mechanical stimuli, biochemical signaling, and genetic prerequisites. Because of the large body of data being produced by modern laboratories, solving this complexity requires the application of computational techniques and analyses. In the last two decades, computational models combined with wet-lab experiments have advanced our understanding of plant organ morphogenesis. Here, we provide a comprehensive review of the most important achievements in the field of computational plant morphodynamics. We present a brief history from the earliest attempts to describe plant forms using algorithmic pattern generation to the evolution of quantitative cell-based models fueled by increasing computational power. We then provide an overview of the most common types of "digital plant" paradigms, and demonstrate how models benefit from diverse techniques used to describe cell growth mechanics. Finally, we highlight the development of computational frameworks designed to resolve organ shape complexity through integration of mechanical, biochemical, and genetic cues into a quantitative standardized and user-friendly environment.

Morphological switch is associated with increase in cell-cell contacts, ALP, and confluence above a minimum island area to perimeter ratio

During osteogenic differentiation in vitro, stem-like cells seeded at a low-density spread and are isolated. As the cells proliferate and mature, they become more cuboidal in shape with more cell-cell contacts. However, the coordination of this switch in cell morphology from elongated to cuboidal, cell-cell contacts, and differentiation is not known. In this work, we present results from experiments and a simulation of cell proliferation on protein-micropatterned islands that, independent of island size (25-1,000 μm) or shape (circles, squares, and hollow squares), shows a distinct morphological switch that is better described as a function of island confluence than time in culture, the standard measure in cell culture experiments. The simulation and experiments show cell morphology and island cell density versus confluence collapse to a single curve for all islands if the island area to perimeter ratio is ≥25 μm. Cell-cell contacts in the simulation and alkaline phosphatase (ALP) expression in experiments, a common marker for osteogenic differentiation, show exponential growth with confluence, rapidly increasing after the switch at ≈0.5 confluence. Furthermore, cell morphology, density, contacts, and ALP are better predicted by confluence than time in culture. The variability with time in culture leads to challenges in not only interpreting data but also in comparing data across research laboratories. This simulation can be used to predict cell behavior on different size and shape islands and to plan and optimize experiments that explore cell behavior as a function of a wide range of island geometries.

How and why to build a mathematical model: A case study using prion aggregation

Biological systems are inherently complex, and the increasing level of detail with which we are able to experimentally probe such systems continually reveals new complexity. Fortunately, mathematical models are uniquely positioned to provide a tool suitable for rigorous analysis, hypothesis generation, and connecting results from isolated in vitro experiments with results from in vivo and whole-organism studies. However, developing useful mathematical models is challenging because of the often different domains of knowledge required in both math and biology. In this work, we endeavor to provide a useful guide for researchers interested in incorporating mathematical modeling into their scientific process. We advocate for the use of conceptual diagrams as a starting place to anchor researchers from both domains. These diagrams are useful for simplifying the biological process in question and distinguishing the essential components. Not only do they serve as the basis for developing a variety of mathematical models, but they ensure that any mathematical formulation of the biological system is led primarily by scientific questions. We provide a specific example of this process from our own work in studying prion aggregation to show the power of mathematical models to synergistically interact with experiments and push forward biological understanding. Choosing the most suitable model also depends on many different factors, and we consider how to make these choices based on different scales of biological organization and available data. We close by discussing the many opportunities that abound for both experimentalists and modelers to take advantage of collaborative work in this field.

MpFEW RHIZOIDS1 miRNA-Mediated Lateral Inhibition Controls Rhizoid Cell Patterning in Marchantia polymorpha

Lateral inhibition patterns differentiated cell types among equivalent cells during development in bacteria, metazoans, and plants. Tip-growing rhizoid cells develop among flat epidermal cells in the epidermis of the early-diverging land plant Marchantia polymorpha. We show that the majority of rhizoid cells develop individually, but some develop in linear, one-dimensional groups (chains) of between 2 and 7 rhizoid cells in wild-type plants. The distribution of rhizoid cells can be accounted for within a simple cellular automata model of lateral inhibition. The model predicted that in the absence of lateral inhibition, two-dimensional rhizoid cell groups (clusters) form. These can be larger than those formed with lateral inhibition. M. polymorpha rhizoid differentiation is positively regulated by the ROOT HAIR DEFECTIVE SIX-LIKE1 (MpRSL1) basic-helix-loop-helix (bHLH) transcription factor, which is directly repressed by the FEW RHIZOIDS1 (MpFRH1) microRNA (miRNA). To test if MpFRH1 miRNA acts during lateral inhibition, we generated loss-of-function (lof) mutants without the MpFRH1 miRNA. Two-dimensional clusters of rhizoids develop in Mpfrh1^lof mutants as predicted by the model for plants that lack lateral inhibition. Furthermore, two-dimensional clusters of up to 9 rhizoid cells developed in the Mpfrh1^lof mutants compared to a maximum number of 7 observed in wild-type groups. The higher steady-state levels of MpRSL1 mRNA in Mpfrh1^lof mutants indicate that MpFRH1-mediated lateral inhibition involves the repression of MpRSL1 activity. Together, the modeling and genetic data indicate that MpFRH1 miRNA mediates lateral inhibition by repressing MpRSL1 during pattern formation in the M. polymorpha epidermis.

Computational modelling of magnesium stent mechanical performance in a remodelling artery: Effects of multiple remodelling stimuli

Significant research has been conducted in the area of coronary stents/scaffolds made from resorbable metallic and polymeric biomaterials. These next-generation bioabsorbable stents have the potential to completely revolutionise the treatment of coronary artery disease. The primary advantage of resorbable devices over permanent stents is their temporary presence which, from a theoretical point of view, means only a healed coronary artery will be left behind following degradation of the stent potentially eliminating long-term clinical problems associated with permanent stents. The healing of the artery following coronary stent/scaffold implantation is crucial for the long-term safety of these devices. Computational modelling can be used to evaluate the performance of complex stent devices in silico and assist in the design and development and understanding of the next-generation resorbable stents. What is lacking in computational modelling literature is the representation of the active response of the arterial tissue in the weeks and months following stent implantation, ie, neointimal remodelling, in particular for the case of biodegradable stents. In this paper, a computational modelling framework is developed, which accounts for two major physiological stimuli responsible for neointimal remodelling and combined with a magnesium corrosion model that is capable of simulating localised pitting (realistic) stent corrosion. The framework is used to simulate different neointimal growth patterns and to explore the effects the neointimal remodelling has on the mechanical performance (scaffolding support) of the bioabsorbable magnesium stent.

Formal description of plant morphogenesis

Plant morphogenesis may be characterized by complex feedback mechanisms between signals specifying growth and by the growth of the plant body itself. Comprehension of such feedback mechanisms is an ongoing research task and can be aided with formal descriptions of morphogenesis. In this review, we present a number of established mathematical paradigms that are useful to the formal representation of plant shape, and of biomechanical and biochemical signaling. Specifically, we discuss work from a range of research areas including plant biology, material sciences, fluid dynamics, and computer graphics. Treating plants as organized systems of information processing allows us to compare these different mathematical methods in terms of their expressive power of biological hypotheses. This is an attempt to bring together a large number of computational modeling concepts and make them accessible to the analytical as well as empirical student of plant morphogenesis.

Homeostasis and systematic ageing as non-equilibrium phase transitions in computational multicellular organizations

Being a fatal threat to life, the breakdown of homeostasis in tissues is believed to involve multiscale factors ranging from the accumulation of genetic damages to the deregulation of metabolic processes. Here, we present a prototypical multicellular homeostasis model in the form of a two-dimensional stochastic cellular automaton with three cellular states, cell division, cell death and cell cycle arrest, of which the state-updating rules are based on fundamental cell biology. Despite the simplicity, this model illustrates how multicellular organizations can develop into diverse homeostatic patterns with distinct morphologies, turnover rates and lifespans without considering genetic, metabolic or other exogenous variations. Through mean-field analysis and Monte-Carlo simulations, those homeostatic states are found to be classified into extinctive, proliferative and degenerative phases, whereas healthy multicellular organizations evolve from proliferative to degenerative phases over a long time, undergoing a systematic ageing akin to a transition into an absorbing state in non-equilibrium physical systems. It is suggested that the collapse of homeostasis at the multicellular level may originate from the fundamental nature of cell biology regarding the physics of some non-equilibrium processes instead of subcellular details.

Metabolic Engineering and Synthetic Biology: Synergies, Future, and Challenges

The “-omics” era has brought a new set of tools and methods that have created a significant impact on the development of Metabolic Engineering and Synthetic Biology. These fields, rather than working separately, depend on each other to prosper and achieve their individual goals. Synthetic Biology aims to design libraries of genetic components (promoters, coding sequences, terminators, transcriptional factors and their binding sequences, and more), the assembly of devices, genetic circuits and even organism; in addition to obtaining quantitative information for the creation of models that can predict the behavior of biological systems (Cameron et al., 2014). Metabolic engineering seeks for the optimization of cellular processes, endemic to a specific organism, to produce a compound of interest from a substrate, preferably cheap and simple. It uses different databases, libraries of components and conditions to generate the maximum production rate of a desired chemical compound and avoiding inhibitors and conditions that affect the growth rate and other vital functions in the specific organism to achieve these goals; metabolic fluxes manipulation represents an important alternative (Stephanopoulos, 2012).

Mathematical models of tumor cell proliferation: A review of the literature

INTRODUCTION

A defining hallmark of cancer is aberrant cell proliferation. Efforts to understand the generative properties of cancer cells span all biological scales: from genetic deviations and alterations of metabolic pathways to physical stresses due to overcrowding, as well as the effects of therapeutics and the immune system. While these factors have long been studied in the laboratory, mathematical and computational techniques are being increasingly applied to help understand and forecast tumor growth and treatment response. Advantages of mathematical modeling of proliferation include the ability to simulate and predict the spatiotemporal development of tumors across multiple experimental scales. Central to proliferation modeling is the incorporation of available biological data and validation with experimental data. Areas covered: We present an overview of past and current mathematical strategies directed at understanding tumor cell proliferation. We identify areas for mathematical development as motivated by available experimental and clinical evidence, with a particular emphasis on emerging, non-invasive imaging technologies. Expert commentary: The data required to legitimize mathematical models are often difficult or (currently) impossible to obtain. We suggest areas for further investigation to establish mathematical models that more effectively utilize available data to make informed predictions on tumor cell proliferation.

A versatile hybrid agent-based, particle and partial differential equations method to analyze vascular adaptation

Peripheral arterial occlusive disease is a chronic pathology affecting at least 8–12 million people in the USA, typically treated with a vein graft bypass or through the deployment of a stent in order to restore the physiological circulation. Failure of peripheral endovascular interventions occurs at the intersection of vascular biology, biomechanics, and clinical decision making. It is our hypothesis that the majority of endovascular treatment approaches share the same driving mechanisms and that a deep understanding of the adaptation process is pivotal in order to improve the current outcome of the procedure. The postsurgical adaptation of vein graft bypasses offers the perfect example of how the balance between intimal hyperplasia and wall remodeling determines the failure or the success of the intervention. Accordingly, this work presents a versatile computational model able to capture the feedback loop that describes the interaction between events at cellular/tissue level and mechano-environmental conditions. The work here presented is a generalization and an improvement of a previous work by our group of investigators, where an agent-based model uses a cellular automata principle on a fixed hexagonal grid to reproduce the leading events of the graft’s restenosis. The new hybrid model here presented allows a more realistic simulation both of the biological laws that drive the cellular behavior and of the active role of the membranes that separate the various layers of the vein. The novel feature is to use an immersed boundary implementation of a highly viscous flow to represent SMC motility and matrix reorganization in response to graft adaptation. Our implementation is modular, and this makes us able to choose the right compromise between closeness to the physiological reality and complexity of the model. The focus of this paper is to offer a new modular implementation that combines the best features of an agent-based model, continuum mechanics, and particle-tracking methods to cope with the multiscale nature of the adaptation phenomena. This hybrid method allows us to quickly test various hypotheses with a particular attention to cellular motility, a process that we demonstrated should be driven by mechanical homeostasis in order to maintain the right balance between cells and extracellular matrix in order to reproduce a distribution similar to histological experimental data from vein grafts.

Simple mechanical cues could explain adipose tissue morphology

The mechanisms by which organs acquire their functional structure and realize its maintenance (or homeostasis) over time are still largely unknown. In this paper, we investigate this question on adipose tissue. Adipose tissue can represent 20 to 50% of the body weight. Its investigation is key to overcome a large array of metabolic disorders that heavily strike populations worldwide. Adipose tissue consists of lobular clusters of adipocytes surrounded by an organized collagen fiber network. By supplying substrates needed for adipogenesis, vasculature was believed to induce the regroupment of adipocytes near capillary extremities. This paper shows that the emergence of these structures could be explained by simple mechanical interactions between the adipocytes and the collagen fibers. Our assumption is that the fiber network resists the pressure induced by the growing adipocytes and forces them to regroup into clusters. Reciprocally, cell clusters force the fibers to merge into a well-organized network. We validate this hypothesis by means of a two-dimensional Individual Based Model (IBM) of interacting adipocytes and extra-cellular-matrix fiber elements. The model produces structures that compare quantitatively well to the experimental observations. Our model seems to indicate that cell clusters could spontaneously emerge as a result of simple mechanical interactions between cells and fibers and surprisingly, vasculature is not directly needed for these structures to emerge.

Computational Modeling of the Mechanical Performance of a Magnesium Stent Undergoing Uniform and Pitting Corrosion in a Remodeling Artery

Coronary stents made from degradable biomaterials such as magnesium alloy are an emerging technology in the treatment of coronary artery disease. Biodegradable stents provide mechanical support to the artery during the initial scaffolding period after which the artery will have remodeled. The subsequent resorption of the stent biomaterial by the body has potential to reduce the risk associated with long-term placement of these devices, such as in-stent restenosis, late stent thrombosis, and fatigue fracture. Computational modeling such as finite-element analysis has proven to be an extremely useful tool in the continued design and development of these medical devices. What is lacking in computational modeling literature is the representation of the active response of the arterial tissue in the weeks and months following stent implantation, i.e., neointimal remodeling. The phenomenon of neointimal remodeling is particularly interesting and significant in the case of biodegradable stents, when both stent degradation and neointimal remodeling can occur simultaneously, presenting the possibility of a mechanical interaction and transfer of load between the degrading stent and the remodeling artery. In this paper, a computational modeling framework is developed that combines magnesium alloy degradation and neointimal remodeling, which is capable of simulating both uniform (best case) and localized pitting (realistic) stent corrosion in a remodeling artery. The framework is used to evaluate the effects of the neointima on the mechanics of the stent, when the stent is undergoing uniform or pitting corrosion, and to assess the effects of the neointimal formation rate relative to the overall stent degradation rate (for both uniform and pitting conditions).

Hemodynamic Influence on Smooth Muscle Cell Kinetics and Phenotype During Early Vein Graft Adaptation

Pathologic vascular adaptation following local injury is the primary driver for accelerated intimal hyperplasia and an occlusive phenotype. Smooth muscle cell (SMC) proliferation within the wall, and migration into the developing intima, is a major component of this remodeling response. The primary objective in the current study was to investigate the effect of the local biomechanical forces on early vein graft adaptation, specifically focusing on the spatial and temporal response of SMC proliferation and conversion from a contractile to synthetic architecture. Taking advantage of the differential adaptation that occurs during exposure to divergent flow environments, vein grafts were implanted in rabbits to create two distinct flow environments and harvested at times ranging from 2 h to 28 days. Using an algorithm for the virtual reconstruction of unfixed, histologic specimens, immunohistochemical tracking of DNA synthesis, and high-throughput transcriptional analysis, the spatial and temporal changes in graft morphology, cell proliferation, and SMC phenotype were catalogued. Notable findings include a burst of cell proliferation at 7 days post-implantation, which was significantly augmented by exposure to a reduced flow environment. Compared to the adjacent media, proliferation rates were 3-fold greater in the intima, and a specific spatial distribution of these proliferating cells was identified, with a major peak in the sub-endothelial region and a second peak centering on the internal elastic lamina. Genomic markers of a contractile SMC phenotype were reduced as early as 2 h post-implantation and reached a nadir at 7 days. Network analysis of upstream regulatory pathways identified GATA6 and KLF5 as important transcription factors that regulate this shift in SMC phenotype.

A Six-Stage Workflow for Robust Application of Systems Pharmacology

Quantitative and systems pharmacology (QSP) is increasingly being applied in pharmaceutical research and development. One factor critical to the ultimate success of QSP is the establishment of commonly accepted language, technical criteria, and workflows. We propose an integrated workflow that bridges conceptual objectives with underlying technical detail to support the execution, communication, and evaluation of QSP projects.

A Modeling Insight into Adipose-Derived Stem Cell Myogenesis

Adipose-derived stem cells (ASCs) are clinically important in regenerative medicine as they are relatively easy to obtain, are characterized by low morbidity, and can differentiate into myogenic progenitor cells. Although studies have elucidated the principal markers, PAX7, Desmin, MyoD, and MHC, the underlying mechanisms are not completely understood. This motivates the application of computational methods to facilitate greater understanding of such processes. In the following, we present a multi-stage kinetic model comprising a system of ordinary differential equations (ODEs). We sought to model ASC differentiation using data from a static culture, where no strain is applied, and a dynamic culture, where 10% strain is applied. The coefficients of the equations have been modulated by those experimental data points. To correctly represent the trajectories, various switches and a feedback factor based on total cell number have been introduced to better represent the biology of ASC differentiation. Furthermore, the model has then been applied to predict ASC fate for strains different from those used in the experimental conditions and for times longer than the duration of the experiment. Analysis of the results reveals unique characteristics of ASC myogenesis under dynamic conditions of the applied strain.

A Multiscale Computational Framework to Understand Vascular Adaptation

The failure rate for vascular interventions (vein bypass grafting, arterial angioplasty/stenting) remains unacceptably high. Over the past two decades, researchers have applied a wide variety of approaches to investigate the primary failure mechanisms, neointimal hyperplasia and aberrant remodeling of the wall, in an effort to identify novel therapeutic strategies. Despite incremental progress, specific cause/effect linkages among the primary drivers of the pathology, (hemodynamic factors, inflammatory biochemical mediators, cellular effectors) and vascular occlusive phenotype remain lacking. We propose a multiscale computational framework of vascular adaptation to develop a bridge between theory and experimental observation and to provide a method for the systematic testing of relevant clinical hypotheses. Cornerstone to our model is a feedback mechanism between environmental conditions and dynamic tissue plasticity described at the cellular level with an agent based model. Our implementation (i) is modular, (ii) starts from basic mechano-biology principle at the cell level and (iii) facilitates the agile development of the model.

From pairwise to group interactions in games of cyclic dominance

We study the rock-paper-scissors game in structured populations, where the invasion rates determine individual payoffs that govern the process of strategy change. The traditional version of the game is recovered if the payoffs for each potential invasion stem from a single pairwise interaction. However, the transformation of invasion rates to payoffs also allows the usage of larger interaction ranges. In addition to the traditional pairwise interaction, we therefore consider simultaneous interactions with all nearest neighbors, as well as with all nearest and next-nearest neighbors, thus effectively going from single pair to group interactions in games of cyclic dominance. We show that differences in the interaction range affect not only the stationary fractions of strategies but also their relations of dominance. The transition from pairwise to group interactions can thus decelerate and even revert the direction of the invasion between the competing strategies. Like in evolutionary social dilemmas, in games of cyclic dominance, too, the indirect multipoint interactions that are due to group interactions hence play a pivotal role. Our results indicate that, in addition to the invasion rates, the interaction range is at least as important for the maintenance of biodiversity among cyclically competing strategies.

Dynamics between cancer cell subpopulations reveals a model coordinating with both hierarchical and stochastic concepts

Tumors are often heterogeneous in which tumor cells of different phenotypes have distinct properties. For scientific and clinical interests, it is of fundamental importance to understand their properties and the dynamic variations among different phenotypes, specifically under radio- and/or chemo-therapy. Currently there are two controversial models describing tumor heterogeneity, the cancer stem cell (CSC) model and the stochastic model. To clarify the controversy, we measured probabilities of different division types and transitions of cells via in situ immunofluorescence. Based on the experiment data, we constructed a model that combines the CSC with the stochastic concepts, showing the existence of both distinctive CSC subpopulations and the stochastic transitions from NSCCs to CSCs. The results showed that the dynamic variations between CSCs and non-stem cancer cells (NSCCs) can be simulated with the model. Further studies also showed that the model can be used to describe the dynamics of the two subpopulations after radiation treatment. More importantly, analysis demonstrated that the experimental detectable equilibrium CSC proportion can be achieved only when the stochastic transitions from NSCCs to CSCs occur, indicating that tumor heterogeneity may exist in a model coordinating with both the CSC and the stochastic concepts. The mathematic model based on experimental parameters may contribute to a better understanding of the tumor heterogeneity, and provide references on the dynamics of CSC subpopulation during radiotherapy.

A dynamical system that describes vein graft adaptation and failure

Adaptation of vein bypass grafts to the mechanical stresses imposed by the arterial circulation is thought to be the primary determinant for lesion development, yet an understanding of how the various forces dictate local wall remodeling is lacking. We develop a dynamical system that summarizes the complex interplay between the mechanical environment and cell/matrix kinetics, ultimately dictating changes in the vein graft architecture. Based on a systematic mapping of the parameter space, three general remodeling response patterns are observed: (1) shear stabilized intimal thickening, (2) tension induced wall thinning and lumen expansion, and (3) tension stabilized wall thickening. Notable is our observation that the integration of multiple feedback mechanisms leads to a variety of non-linear responses that would be unanticipated by an analysis of each system component independently. This dynamic analysis supports the clinical observation that the majority of vein grafts proceed along an adaptive trajectory, where grafts dilate and mildly thicken in response to the increased tension and shear, but a small portion of the grafts demonstrate a maladaptive phenotype, where progressive inward remodeling and accentuated wall thickening lead to graft failure.

Rule-based model of vein graft remodeling

When vein segments are implanted into the arterial system for use in arterial bypass grafting, adaptation to the higher pressure and flow of the arterial system is accomplished thorough wall thickening and expansion. These early remodeling events have been found to be closely coupled to the local hemodynamic forces, such as shear stress and wall tension, and are believed to be the foundation for later vein graft failure. To further our mechanistic understanding of the cellular and extracellular interactions that lead to global changes in tissue architecture, a rule-based modeling method is developed through the application of basic rules of behaviors for these molecular and cellular activities. In the current method, smooth muscle cell (SMC), extracellular matrix (ECM), and monocytes are selected as the three components that occupy the elements of a grid system that comprise the developing vein graft intima. The probabilities of the cellular behaviors are developed based on data extracted from in vivo experiments. At each time step, the various probabilities are computed and applied to the SMC and ECM elements to determine their next physical state and behavior. One- and two-dimensional models are developed to test and validate the computational approach. The importance of monocyte infiltration, and the associated effect in augmenting extracellular matrix deposition, was evaluated and found to be an important component in model development. Final model validation is performed using an independent set of experiments, where model predictions of intimal growth are evaluated against experimental data obtained from the complex geometry and shear stress patterns offered by a mid-graft focal stenosis, where simulation results show good agreements with the experimental data.

The dynamics of vein graft remodeling induced by hemodynamic forces: a mathematical model

Although vein bypass grafting is one of the primary options for the treatment of arterial occlusive disease and provides satisfactory results at an early stage of the treatment, the patency is limited to a few months in many patients. When the vein is implanted in the arterial system, it adapts to the high flow rate and high pressure of the arterial environment by changing the sizes of its layers, and this remodeling is believed to be a precursor of future graft failure. Hemodynamic forces, such as wall shear stress (WSS) and wall tension, have been recognized as major factors impacting vein graft remodeling. Although a wide range of experimental evidence relating hemodynamic forces to vein graft remodeling has been reported, a comprehensive mathematical model describing the relationship among WSS, wall tension, and the structural adaptation of each individual layer of the vein graft wall is lacking. The current manuscript presents a comprehensive and robust framework for treating the complex interaction between the WSS, wall tension, and the structural adaptation of each individual layer of the vein graft wall. We modeled the intimal and medial area and the radius of external elastic lamina, which in combination dictate luminal narrowing and the propensity for graft occlusion. Central to our model is a logistic relationship between independent and dependent variables to describe the initial increase and later decrease in the growth rate. The detailed understanding of the temporal changes in vein graft morphology that can be extracted from the current model is critical in identifying the dominant contributions to vein graft failure and the further development of strategies to improve their longevity.