Simulating Properties of In Vitro Epithelial Cell Morphogenesis

How do individual epithelial cells (ECs) organize into multicellular structures? ECs are studied in vitro to help answer that question. Characteristic growth features include stable cyst formation in embedded culture, inverted cyst formation in suspension culture, and lumen formation in overlay culture. Formation of these characteristic structures is believed to be a consequence of an intrinsic program of differentiation and de-differentiation. To help discover how such a program may function, we developed an in silico analogue in which space, events, and time are discretized. Software agents and objects represent cells and components of the environment. “Cells” act independently. The “program” governing their behavior is embedded within each in the form of axioms and an inflexible decisional process. Relationships between the axioms and recognized cell functions are specified. Interactions between “cells” and environment components during simulation give rise to a complex in silico phenotype characterized by context-dependent structures that mimic counterparts observed in four different in vitro culture conditions: a targeted set of in vitro phenotypic attributes was matched by in silico attributes. However, for a particular growth condition, the analogue failed to exhibit behaviors characteristic of functionally polarized ECs. We solved this problem by following an iterative refinement method that improved the first analogue and led to a second: it exhibited characteristic differentiation and growth properties in all simulated growth conditions. It is the first model to simultaneously provide a representation of nonpolarized and structurally polarized cell types, and a mechanism for their interconversion. The second analogue also uses an inflexible axiomatic program. When specific axioms are relaxed, growths strikingly characteristic of cancerous and precancerous lesions are observed. In one case, the simulated cause is aberrant matrix production. Analogue design facilitates gaining deeper insight into such phenomena by making it easy to replace low-resolution components with increasingly detailed and realistic components

A language for modeling the morphogenesis process of biological system using MAS

The biological systems are Complex Adaptive Systems that can build a complex structure using interaction between cell and its environment. The process is called morphogenesis, it is study by developmental biology and applied by tissue engineering and regenerative medicine. Computational models are used by researchers to organize and communicate ideas, test hypotheses and abstractly represent behaviors, etc. Although, Agent Based Models are used in biological systems, it is necessary appropriate domain-oriented language. This work shows the progress made in the definition of a language for modeling and simulation of such systems, in order to hide details of implementation and execution of a simulation. The language concepts are demonstrated using, as study case, the morphogenesis of a heart valve in vitro.Sociedad Argentina de Informática e Investigación Operativa (SADIO

A language for modeling the morphogenesis process of biological system using MAS

The biological systems are Complex Adaptive Systems that can build a complex structure using interaction between cell and its environment. The process is called morphogenesis, it is study by developmental biology and applied by tissue engineering and regenerative medicine. Computational models are used by researchers to organize and communicate ideas, test hypotheses and abstractly represent behaviors, etc. Although, Agent Based Models are used in biological systems, it is necessary appropriate domain-oriented language. This work shows the progress made in the definition of a language for modeling and simulation of such systems, in order to hide details of implementation and execution of a simulation. The language concepts are demonstrated using, as study case, the morphogenesis of a heart valve in vitro.Sociedad Argentina de Informática e Investigación Operativa (SADIO

A language for modeling the morphogenesis process of biological system using MAS

The biological systems are Complex Adaptive Systems that can build a complex structure using interaction between cell and its environment. The process is called morphogenesis, it is study by developmental biology and applied by tissue engineering and regenerative medicine. Computational models are used by researchers to organize and communicate ideas, test hypotheses and abstractly represent behaviors, etc. Although, Agent Based Models are used in biological systems, it is necessary appropriate domain-oriented language. This work shows the progress made in the definition of a language for modeling and simulation of such systems, in order to hide details of implementation and execution of a simulation. The language concepts are demonstrated using, as study case, the morphogenesis of a heart valve in vitro.Sociedad Argentina de Informática e Investigación Operativa (SADIO

Multi-agent model of hepatitis C virus infection

Objectives: The objective of this study is to design a method for modeling hepatitis C virus (HCV) infection using multi-agent simulation and to verify it in practice. Methods and materials: In this paper, first, the modeling of HCV infection using a multi-agent system is compared with the most commonly used model type, which is based on differential equations. Then, the implementation and results of the model using a multi-agent simulation is presented. To find the values of the parameters used in the model, a method using inverted simulation flow and genetic algorithm is proposed. All of the data regarding HCV infection are taken from the paper describing the model based on the differential equation to which the proposed method is compared. Results: Important advantages of the proposed method are noted and demonstrated; these include flexibility, clarity, re-usability and the possibility to model more complex dependencies. Then, the simulation framework that uses the proposed approach is successfully implemented in C++ and is verified by comparing it to the approach based on differential equations. The verification proves that an objective function that performs the best is the function that minimizes the maximal differences in the data. Finally, an analysis of one of the already known models is performed, and it is proved that it incorrectly models a decay in the hepatocytes number by 40%. Conclusions: The proposed method has many advantages in comparison to the currently used model types and can be used successfully for analyzing HCV infection. With almost no modifications, it can also be used for other types of viral infections

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

Computational investigation of epithelial cell dynamic phenotype in vitro

<p>Abstract</p> <p>Background</p> <p>When grown in three-dimensional (3D) cultures, epithelial cells typically form cystic organoids that recapitulate cardinal features of in vivo epithelial structures. Characterizing essential cell actions and their roles, which constitute the system's dynamic phenotype, is critical to gaining deeper insight into the cystogenesis phenomena.</p> <p>Methods</p> <p>Starting with an earlier in silico epithelial analogue (ISEA1) that validated for several Madin-Darby canine kidney (MDCK) epithelial cell culture attributes, we built a revised analogue (ISEA2) to increase overlap between analogue and cell culture traits. Both analogues used agent-based, discrete event methods. A set of axioms determined ISEA behaviors; together, they specified the analogue's operating principles. A new experimentation framework enabled tracking relative axiom use and roles during simulated cystogenesis along with establishment of the consequences of their disruption.</p> <p>Results</p> <p>ISEA2 consistently produced convex cystic structures in a simulated embedded culture. Axiom use measures provided detailed descriptions of the analogue's dynamic phenotype. Dysregulating key cell death and division axioms led to disorganized structures. Adhering to either axiom less than 80% of the time caused ISEA1 to form easily identified morphological changes. ISEA2 was more robust to identical dysregulation. Both dysregulated analogues exhibited characteristics that resembled those associated with an in vitro model of early glandular epithelial cancer.</p> <p>Conclusion</p> <p>We documented the causal chains of events, and their relative roles, responsible for simulated cystogenesis. The results stand as an early hypothesis–a theory–of how individual MDCK cell actions give rise to consistently roundish, cystic organoids.</p

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

A computational approach to resolve cell level contributions to early glandular epithelial cancer progression

<p>Abstract</p> <p>Background</p> <p>Three-dimensional (3D) embedded cell cultures provide an appropriate physiological environment to reconstruct features of early glandular epithelial cancer. Although these are orders of magnitude simpler than tissues, they too are complex systems that have proven challenging to understand. We used agent-based, discrete event simulation modeling methods to build working hypotheses of mechanisms of epithelial 3D culture phenotype and early cancer progression. Starting with an earlier software analogue, we validated an improved in silico epithelial analogue (ISEA) for cardinal features of a normally developed MDCK cyst. A set of axiomatic operating principles defined simulated cell actions. We explored selective disruption of individual simulated cell actions. New framework features enabled recording detailed measures of ISEA cell activities and morphology.</p> <p>Results</p> <p>Enabled by a small set of cell operating principles, ISEA cells multiplied and self-organized into cyst-like structures that mimicked those of MDCK cells in a 3D embedded cell culture. Selective disruption of "anoikis" or directional cell division caused the ISEA to develop phenotypic features resembling those of in vitro tumor reconstruction models and cancerous tissues in vivo. Disrupting either process, or both, altered cell activity patterns that resulted in morphologically similar outcomes. Increased disruption led to a prolonged presence of intraluminal cells.</p> <p>Conclusions</p> <p>ISEA mechanisms, behaviors, and morphological properties may have biological counterparts. To the extent that in silico-to-in vitro mappings are valid, the results suggest plausible, additional mechanisms of in vitro cancer reconstruction or reversion, and raise potentially significant implications for early cancer diagnosis based on histology. Further ISEA development and use are expected to provide a viable platform to complement in vitro methods for unraveling the mechanistic basis of epithelial morphogenesis and cancer progression.</p

Identifying the Rules of Engagement Enabling Leukocyte Rolling, Activation, and Adhesion

The LFA-1 integrin plays a pivotal role in sustained leukocyte adhesion to the endothelial surface, which is a precondition for leukocyte recruitment into inflammation sites. Strong correlative evidence implicates LFA-1 clustering as being essential for sustained adhesion, and it may also facilitate rebinding events with its ligand ICAM-1. We cannot challenge those hypotheses directly because it is infeasible to measure either process during leukocyte adhesion following rolling. The alternative approach undertaken was to challenge the hypothesized mechanisms by experimenting on validated, working counterparts: simulations in which diffusible, LFA1 objects on the surfaces of quasi-autonomous leukocytes interact with simulated, diffusible, ICAM1 objects on endothelial surfaces during simulated adhesion following rolling. We used object-oriented, agent-based methods to build and execute multi-level, multi-attribute analogues of leukocytes and endothelial surfaces. Validation was achieved across different experimental conditions, in vitro, ex vivo, and in vivo, at both the individual cell and population levels. Because those mechanisms exhibit all of the characteristics of biological mechanisms, they can stand as a concrete, working theory about detailed events occurring at the leukocyte–surface interface during leukocyte rolling and adhesion experiments. We challenged mechanistic hypotheses by conducting experiments in which the consequences of multiple mechanistic events were tracked. We quantified rebinding events between individual components under different conditions, and the role of LFA1 clustering in sustaining leukocyte–surface adhesion and in improving adhesion efficiency. Early during simulations ICAM1 rebinding (to LFA1) but not LFA1 rebinding (to ICAM1) was enhanced by clustering. Later, clustering caused both types of rebinding events to increase. We discovered that clustering was not necessary to achieve adhesion as long as LFA1 and ICAM1 object densities were above a critical level. Importantly, at low densities LFA1 clustering enabled improved efficiency: adhesion exhibited measurable, cell level positive cooperativity