A cell-centered, agent-based framework that enables flexible environment granularities
© Kennedy et al. 2016
Received: 25 October 2015
Accepted: 20 January 2016
Published: 2 February 2016
Mechanistic explanations of cell-level phenomena typically adopt an observer perspective. Explanations developed from a cell’s perspective may offer new insights. Agent-based models lend themselves to model from an individual perspective, and existing agent-based models generally utilize a regular lattice-based environment. A framework which utilizes a cell’s perspective in an off-lattice environment could improve the overall understanding of biological phenomena.
We present an agent-based, discrete event framework, with a demonstrative focus on biomimetic agents. The framework was first developed in 2-dimensions and then extended, with a subset of behaviors, to 3-dimensions. The framework is expected to facilitate studies of more complex biological phenomena through exploitation of a dynamic Delaunay and Voronoi off-lattice environment. We used the framework to model biological cells and to specifically demonstrate basic biological cell behaviors in two- and three-dimensional space. Potential use cases are highlighted, suggesting the utility of the framework in various scenarios.
The framework presented in this manuscript expands on existing cell- and agent-centered methods by offering a new perspective in an off-lattice environment. As the demand for biomimetic models grows, the demand for new methods, such as the presented Delaunay and Voronoi framework, is expected to increase.
There is a growing need for new methods to cater to increasingly complex biological models, which aim to provide better mechanistic explanations of biological phenomena. We describe and make available an early stage simulation framework that enables a dynamic Delaunay and Voronoi (D/V) off-lattice environment to be created and used by biomimetic agents. This framework can accommodate a variety of uses, among which are those that adopt a cell- or entity-centered perspective. It is intended to expand the repertoire available to modelers by helping to make irregular grids part of their standard toolkit.
Three scientifically important needs provided motivation for this simulation framework. When one is engaged in improving mechanistic explanations of cell level phenomena, we aim to make it easier to: 1) modify mechanistic granularity within a simulation where and when that is needed to improve explanatory insight ; 2) alter the focus or perspective of a simulation in much the same way wet-lab biologists adjust the focus of their experiments; 3) acknowledge, identify, represent, and begin explaining multilevel uncertainties within and between comparable observations made using different wet-lab systems. An important biological focus for us has been improving mechanistic explanations of particular phenomena at cell and multicell levels, primarily in vitro, during such fundamental processes as wound healing, formation and maintenance of a single central lumen during early cystogenesis , and multicellular collective invasion that is characteristic of carcinomas [3, 4].
Research presented in the cited papers employ methods representative of the field. The focus of experiments necessarily shifts from one aspect to another as phenomena change and evolve. For example, early in cystogenesis , attention may focus on pre-luminal events occurring at the apical interface of two or more cells. Multiple visualization methods are employed; examples include differentially staining particular proteins or using cells capable of expressing fluorescent versions of particular proteins. Events elsewhere in the multicell structures are deemed less critical, and thus may not be measured or observed. Later in the cystogenesis process, attention may shift to characteristics of whole cysts; for example, fine grain details at cell-cell interfaces may be deemed less influential to evolving system-level phenotype. Fluorescent staining of cell nuclei may enable measuring the relative arrangement of cells in a cyst, yet data identifying locations of cell boundaries and/or cell-cell interfaces may not be available because they were not visualized or measured. When studying cell behavior in tumor organoids , an experiment’s scale of focus may shift between numbers of invasive multicell leader structures emanating from particular tumor organoids to behavior of individual leader cells within a single leader structure.
The preceding examples illustrate that in such experiments there is always uncertainty about aspects and features not measured. Ideally, we would like to acknowledge, even represent, such uncertainty within our simulation models by avoiding over-committing to the simulation of particular details when and where there is little or no data against which to validate those commitments. However, doing so is challenging  and even problematic if one begins a modeling and simulation project by specifying in advance how and at what level of detail a mechanism’s spaces, entities, and activities will be represented. Those commitments begin when the modeler follows dominant practices and, for example, simply chooses in advance to use a regular grid. We envision the D/V grid providing the capability to fine or coarsen a model’s local scale of focus. An irregular grid allows us to selectively fine or coarsen the scale in one region more or less than in other regions, thereby adjusting granularity.
Explanations of, and theories about, cell level and multicell phenomena, such as those cited above, are offered based on inferences drawn from observing many images, yet only a few images, thought by the said research authors to be typical or representative, are included in publications. A natural inclination for a modeler is to focus first on the images provided, begin creating computational descriptions of those images, and then hypothesize and describe possible event sequences to fill the gaps separating the created depictions. So doing is made straightforward when one commits upfront to using a regular grid. Once that commitment is made, the modeler is at risk of restricting attention to those mechanistic entities and activities that are easily instantiated using the selected grid.
Early commitment to a regular grid also supports favoring the observer’s perspective rather than the biology’s perspective. Longer term, we anticipate that “decisions” made by each cell will contribute to a greater phenomenon based only on information available in its local environment—the cell’s perspective—and the resolution or granularity with which the cell samples or measures those features, thereby enabling observation and validation from cell-and system-levels. Further, we can anticipate that location within the larger system will lead to location-dependent differences in mechanisms. We envision the D/V grid being used by an individual agent representing an individual cell (cell cluster or other entity) to fine or coarsen its local focus. Again, because our framework utilizes an irregular grid, the cell agent will be able to selectively fine or coarsen granularity at one location more or less than elsewhere.
We expect researchers interested in our framework will have familiarity with agent-based modeling (ABM) techniques, and we encourage those interested in further details to seek Macal and North’s tutorial . Agent-based methods are adept at modeling natural phenomena, and an ABM approach helps facilitate the needs of the cell-centered framework. ABMs can represent complex entities and behaviors at an agent-level; each entity, or agent, implements a set of rules and has its own set of properties. Numerous agent classes can exist within a model, and the interactions between agents can lead to global, or emergent, properties. Because mechanisms in ABM are often informally specified, with few or no theoretical constraints, natural phenomena can be simulated more directly, in contrast to those methods that require arching principles which may bias simulation results. ABM has been applied to a wide range of disciplines, from social networking to biology. D/V grids are a natural fit with ABM, as D/V grids facilitate the emergence of a spatial structure rather than one that is based upon prior assumptions.
Lattice and off-lattice spatial representations
Spatial representations in cell modeling can be categorized as lattice (discrete) or off-lattice (continuous); we focus the framework on an off-lattice spatial representation, D/V grids, which give us the capability to utilize neighborhood-based information at an individual agent’s level.
Lattice-based models limit cells to defined spatial boundaries or grids and are most frequently used within the cellular Potts model (CPM). The CPM is an extensible technique to simulate a variety of cell behaviors, such as morphogenesis, sorting, adhesion, as well as tissue simulations and virtual tumors [6–14]. In the CPM, cells exist within a lattice and are connected with bonds represented by various equations; cell behaviors are governed largely by “energies.” They have the ability to model cells in a variety of shapes and sizes and can utilize a number of grid representations. Conventional agent-based models utilize lattices and the correlating neighborhoods for most agent-to-agent interactions. Hexagonal grids are increasingly used  and have some inherent advantages.
There has been only limited work incorporating off-lattice D/V grids into ABM methods. Most notably, the Meyer-Hermann group has developed dynamic force-based models that utilize D/V grids , which have been used in studies of multicellular tumor spheroids [35, 36]. Drasdo et al.  built upon the models from the Meyer-Hermann group to model three-dimensional (3D) tumor growth. An additional example of off-lattice modeling is provided by Dillon et al.  in their model of cell growth and division. Pathmanathan et al.  describe tissue models that utilize D/V grids, but not in a direct ABM syntax. Honda used D/V grids to represent the wing epidermis in butterflies . The Chaste software library  supports models with D/V grids, but largely focuses on cardiac modeling. Virtual Cell  provides a framework for off-lattice representations and provides mathematical solvers but without an agent-based core.
D/V grids are most appropriate for modeling use cases where the spatial medium (or organization) is part of the experimental target, an aspect of the model to be studied, as opposed to an assumption programmed by the developer. Similar modeling objectives are sought by network-based models where the organization and interactions of model components arise as a consequence of the model. It is in this manner that we aim to differentiate our framework from existing approaches.
Framework and Methods
The framework is implemented in MASON  with heavy use of the Visualization Toolkit (VTK) ; it is available at simtk.org . MASON is a powerful Java-based simulation framework, and notable advantages of MASON include speed, a distinct separation of display objects from the core simulation, and the built-in ability to checkpoint simulations . VTK is an advanced system for image processing and 3D graphics. As an open source toolkit, VTK has been incorporated into many systems, has a rich set of algorithms, and has been extensively used for medical imaging. While it is written in C++, wrappers exist that make it accessible to other programming languages. VTK supports Delaunay triangulations, both in two-dimensional (2D) and 3D space, but it does not have built in mechanisms to compute Voronoi tessellations. With MASON, we utilize (Java) VTK and build upon its capabilities to allow for Voronoi tessellations. MASON and VTK allow for the 2- and 3D display of data.
Improving explanatory insight into the generation of biological phenomena, for example those observed during wound healing and morphogenesis, is critical to advancing biomedical science. Cell-centered models can represent the behavior of systems of individual interacting cells and tissue functional units to study their collective behaviors. While tracking the behavior of an individual cell during a wet-lab experiment can at best be challenging, system-level measurements are routine. We engineered the framework to facilitate study of biological cells and have demonstrated use by modeling basic cell behavior, from a cell’s perspective, utilizing 2D and 3D D/V grids as a basis for V-cells.
The framework can support an arbitrary number of V-cell agent types, representing many entities. Each V-cell type is recognized by specific behaviors in the method and operates from a V-cell-based perspective. V-cells agents rely on their surrounding neighborhoods to make decisions. Partly because V-cell agents are not on a grid, decisions are more isolatable from the surrounding environment. The use of a D/V grid means that agents have no regular neighborhood shape. Instead, local computations are made for an individual V-cell. Neighbors are considered as V-cells that are a single hop away from the observing V-cell, meaning other V-cell centers to which the V-cell center point is directly connected (there exists a edge between the V-cells). Since V-cell center points make up the nodes for the triangulation, neighbors can be identified (calculated) from the triangulation. When a V-cell agent moves, its center point moves and the triangulation must be recalculated. The resulting V-cells are calculated and change location and shape automatically depending on the agent’s rules and agenda; many V-cells are impacted when a single V-cell moves. Similarly, when a given V-cell grows or shrinks, surrounding V-cells may be pushed or pulled.
Use case: biological cells
Agent behaviors are governed by axioms developed by domain experts. We have implemented mechanisms to mimic movement, growth, division, and death in 2D and 3D space. Lec and mec can move randomly or in a general direction, based on the underlying Cartesian space. For growth, the agent increases its area or volume, while pushing out neighboring agents. In division, the agent creates a copy of itself. The copy, or daughter agent, inherits traits from its parent. Additionally, agents can be removed from the system. Removal can map to a cell disintegration following death. When an agent is removed, surrounding agents can compete to fill the void. V-cell absorption is also demonstrated; specifically, a mec cell is shown to be absorbed by a cyst. We have implemented an additional behavior in 2D: V-cell rotation, in which a set of agents rotates around another agent. We demonstrate lec agents rotating around a lmn agent in Additional file 4. It is important to note that we have not attempted to implement a fully falsifiable hypothesis for an actual referent. Rather, we have developed a framework whereby the individual agents make decisions rather than utilizing observer requirements, and we have specified several examples through the form of basic use cases. We do not, at this stage, believe it necessary to have a fully falsifiable referent to demonstrate how this method could be used to glean new explanatory insights from an agent’s perspective.
Additional files 4, 5, 6, 7, 8, 9 and 10 are illustrative videos depicting the aforementioned behaviors. Additional file 4 demonstrates 2D agent rotation. Additional file 5 shows 2D agent movement, division, and death, while Additional file 6 demonstrates 2D agent growth. Additional file 7 shows 2D absorption. We visualize 3D agent movement, division, and death in Additional file 8; we show 3D agent growth in Additional file 9. Additional file 10 shows absorption in 3D. Finally, Additional file 11 shows the underlying tetrahedral mesh similar to what is used in Additional file 10. The speed of each video is decreased for clarity.
We provide a new framework for use within the realm of cell- or entity-centered modeling. Although we demonstrate only a few basic cell behaviors, the framework is designed to be extensible to accommodate a variety of use cases. Examples include aspects of morphogenesis and tissue repair in vitro and in vivo. Objects that map to particular mobile molecules, including one or more drugs, can be added. Doing so is expected to facilitate using simulations to explore competing explanations of disease progression in the presence and absence of a therapeutic intervention.
Networks naturally lend themselves to study via D/V grids. Nodes can map to V-cell centers, with the range of “awareness” or capacity denoted by the size of the V-cell. Agents at nodes can communicate with other agents at other nodes based on the number of nodes between them. In this manner, irregular grid structures can be used to model, for example, neuronal tissue, where mobile objects representing transmitters can traverse the graph.
The framework’s design enables study of interfacial regions between neighboring entities. In the demonstrations, an interface between agents need not have a specific biological counterpart. Rather, it can be mapped to cell-cell and/or cell-environment interfaces. Dynamic interfaces between agents during simulations are necessarily precise. The mapping to biology, however, should include appropriate uncertainties.
We demonstrate and make available a new framework that expands upon traditional cell- and agent-centered methods. Our framework provides opportunities for several use cases. The methods can be applied across domains as a framework. V-cells can be used to represent any entity, at various granularities. We foresee these methods facilitating new ways to think about and better explain biological phenomena.
Availability and Requirements
Our framework is available as source code, and we also include the VTK Java bindings and MASON Java archive (JAR) files in the download. Installation instructions on Ubuntu are included as documentation within the download.
Project Name: DVFramework
Project home page: https://simtk.org/home/dvframework
Operating system: Platform independent
Development Environment: Java 1.6
Other requirements: MASON 15, VTK 5.8
License: GNU General Public License v3
Any restrictions to use by non-academics: none
We thank Brenden Petersen and Scott Christley for their helpful discussion and commentary.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Kirschner DE, Hunt CA, Marino S, Fallahi-Sichani M, Linderman JJ. Tuneable resolution as a systems biology approach for multi-scale, multi-compartment computational models. Wiley Interdiscip Rev Syst Biol and Med. 2014;6(4):289–309.View ArticleGoogle Scholar
- Gálvez-Santisteban M, Rodriguez-Fraticelli AE, Bryant DM, Vergarajauregui S, Yasuda T, Bañón-Rodríguez I, et al. Synaptotagmin-like proteins control the formation of a single apical membrane domain in epithelial cells. Nat Cell Biol. 2012;14(8):838–49.View ArticlePubMedPubMed CentralGoogle Scholar
- Ewald AJ, Huebner RJ, Palsdottir H, Lee JK, Perez MJ, Jorgens DM, et al. Mammary collective cell migration involves transient loss of epithelial features and individual cell migration within the epithelium. J Cell Sci. 2012;125(Pt 11):2638–54.View ArticlePubMedPubMed CentralGoogle Scholar
- Cheung KJ, Gabrielson E, Werb Z, Ewald AJ. Collective Invasion in Breast Cancer Requires a Conserved Basal Epithelial Program. Cell. 2013;155(7):1639–51.View ArticlePubMedPubMed CentralGoogle Scholar
- Macal CM, North MJ. Tutorial on agent-based modelling and simulation. J Simulat. 2010;4(3):151–62.View ArticleGoogle Scholar
- Graner F, Glazier JA. Simulation of biological cell sorting using a two-dimensional extended Potts model. Phys Rev Lett. 1992;69:2013–6.View ArticlePubMedGoogle Scholar
- Merks RMH, Glazier JA. A cell-centered approach to developmental biology. Physica A. 2005;352(1):113–30.View ArticleGoogle Scholar
- Ouchi NB, Glazier JA, Rieu JP, Upadhyaya A, Sawada Y. Improving the realism of the cellular Potts model in simulations of biological cells. Physica A. 2003;329(3–4):451–8.View ArticleGoogle Scholar
- Glazier JA, Graner F. Simulation of the differential adhesion driven rearrangement of biological cells. Phys Rev E. 1993;47:2128–54.View ArticleGoogle Scholar
- Savill NJ, Hogeweg P. Modelling Morphogenesis: From Single Cells to Crawling Slugs. J Theor Biol. 1997;184(3):229–35.View ArticleGoogle Scholar
- Izaguirre JA, Chaturvedi R, Huang C, Cickovski T, Coffland J, Thomas G, et al. CompuCell, a multi-model framework for simulation of morphogenesis. Bioinformatics. 2004;20(7):1129–37.View ArticlePubMedGoogle Scholar
- Hogeweg P. Evolving Mechanisms of Morphogenesis: on the Interplay between Differential Adhesion and Cell Differentiation. J Theor Biol. 2000;203(4):317–33.View ArticlePubMedGoogle Scholar
- Zhang L, Athale CA, Deisboeck TS. Development of a three-dimensional multiscale agent-based tumor model: Simulating gene-protein interaction profiles, cell phenotypes and multicellular patterns in brain cancer. J Theor Biol. 2007;244(1):96–107.View ArticlePubMedGoogle Scholar
- Voss-Böhme A. Multi-Scale Modeling in Morphogenesis: A Critical Analysis of the Cellular Potts Model. PLoS One. 2012;7(9), e42852.View ArticlePubMedPubMed CentralGoogle Scholar
- Rejniak KA, Anderson ARA. Hybrid models of tumor growth. Wiley Interdiscip Rev Syst Biol and Med. 2011;3(1):115–25.View ArticleGoogle Scholar
- Farhadifar R, Röper JC, Aigouy B, Eaton S, Jülicher F. The Influence of Cell Mechanics, Cell-Cell Interactions, and Proliferation on Epithelial Packing. Curr Biol. 2007;17(24):2095–104.View ArticlePubMedGoogle Scholar
- Landsberg KP, Farhadifar R, Ranft J, Umetsu D, Widmann TJ, Bittig T, et al. Increased Cell Bond Tension Governs Cell Sorting at the Drosophila Anteroposterior Compartment Boundary. Curr Biol. 2009;19(22):1950–5.View ArticlePubMedGoogle Scholar
- Weliky M, Oster G. The mechanical basis of cell rearrangement. I. Epithelial morphogenesis during Fundulus epiboly. Development. 1990;109(2):373–86.PubMedGoogle Scholar
- Palsson E. A three-dimensional model of cell movement in multicellular systems. Future Gener Comp Sy. 2001;17(7):835–52.View ArticleGoogle Scholar
- Newman TJ. Modeling multicellular systems using subcellular elements. Math Biosci Eng. 2005;2(3):613–24.View ArticlePubMedGoogle Scholar
- Galle J, Sittig D, Hanisch I, Wobus M, Wandel E, Loeffler M, et al. Individual Cell-Based Models of Tumor-Environment Interactions. Am J Pathol. 2006;169(5):1802–11.View ArticlePubMedPubMed CentralGoogle Scholar
- Christley S, Lee B, Dai X, Nie Q. Integrative multicellular biological modeling: a case study of 3D epidermal development using GPU algorithms. BMC Syst Biol. 2010;4:107.View ArticlePubMedPubMed CentralGoogle Scholar
- Honda H. Description of cellular patterns by Dirichlet domains: The two-dimensional case. J Theor Biol. 1978;72(3):523–43.View ArticlePubMedGoogle Scholar
- Jeon J, Quaranta V, Cummings PT. An Off-Lattice Hybrid Discrete-Continuum Model of Tumor Growth and Invasion. Biophys J. 2010;98(1):37–47.View ArticlePubMedPubMed CentralGoogle Scholar
- Dyson L, Maini PK, Baker RE. Macroscopic limits of individual-based models for motile cell populations with volume exclusion. Phys Rev E. 2012;86:031903.View ArticleGoogle Scholar
- Drasdo D. On Selected Individual-based Approaches to the Dynamics in Multicellular Systems. In: Alt W, Chaplain M, Griebel M, Lenz J, editors. Polymer and Cell Dynamics. Mathematics and Biosciences in Interaction. Switzerland: Birkhäuser Verlag Basel; 2003. p. 169–203.View ArticleGoogle Scholar
- Meineke FA, Potten CS, Loeffler M. Cell migration and organization in the intestinal crypt using a lattice-free model. Cell Proliferat. 2001;34(4):253–66.View ArticleGoogle Scholar
- Galle J, Hoffmann M, Aust G. From single cells to tissue architecture-a bottom-up approach to modelling the spatio-temporal organisation of complex multi-cellular systems. J Math Biol. 2009;58(1–2):261–83.View ArticlePubMedGoogle Scholar
- Murray PJ, Edwards CM, Tindall MJ, Maini PK. From a discrete to a continuum model of cell dynamics in one dimension. Phys Rev E. 2009;80:031912.View ArticleGoogle Scholar
- Byrne H, Drasdo D. Individual-based and continuum models of growing cell populations: a comparison. J Math Biol. 2009;58(4–5):657–87.View ArticlePubMedGoogle Scholar
- Van Leeuwen IMM, Mirams GR, Walter A, Fletcher A, Murray P, Osborne J, et al. An integrative computational model for intestinal tissue renewal. Cell Proliferat. 2009;42(5):617–36.View ArticleGoogle Scholar
- Bowyer A. Computing Dirichlet tessellations. Comput J. 1981;24(2):162–6.View ArticleGoogle Scholar
- Watson DF. Computing the n-dimensional Delaunay tessellation with application to Voronoi polytopes. Comput J. 1981;24(2):167–72.View ArticleGoogle Scholar
- Schaller G, Meyer-Hermann M. Kinetic and dynamic Delaunay tetrahedralizations in three dimensions. Comput Phys Commun. 2004;162(1):9–23.View ArticleGoogle Scholar
- Schaller G, Meyer-Hermann M. Multicellular tumor spheroid in an off-lattice Voronoi-Delaunay cell model. Phys Rev E. 2005;71:051910.View ArticleGoogle Scholar
- Meyer-Hermann M. Delaunay-Object-Dynamics: Cell Mechanics with a 3D Kinetic and Dynamic Weighted Delaunay-Triangulation. In: Schnell S, Maini PK, Newman SA, Newman TJ, editors. Multiscale Modeling of Developmental Systems, volume 81 of Current Topics in Developmental Biology. London, United Kingdom: Academic; 2008. p. 373–99.View ArticleGoogle Scholar
- Drasdo D, Jagiella N, Ramis-Conde I, Vignon-Clementel I, Weens W. Modeling steps from a benign tumor to an invasive cancer: examples of intrinsically multi-scale problems. In: Chauviere A, Preziosi L, Verdier C, editors. Cell Mechanics: From Single Scale-Based Models to Multiscale Modeling. Boca Raton, Florida: Chapman & Hall CRC; 2010. p. 379–417.View ArticleGoogle Scholar
- Dillon R, Owen M, Painter K. A single-cell based model of multicellular growth using the immersed boundary method. In: Khoo BC, Li Z, Lin P, editors. Contemporary Mathematics: Moving Interface Problems and Applications in Fluid Dynamics. Providence, Rhode Island: American Mathematical Society; 2008. p. 1–16.View ArticleGoogle Scholar
- Pathmanathan P, Cooper J, Fletcher A, Mirams G, Murray P, Osborne J, et al. A computational study ofdiscrete mechanical tissue models. Phys Biol. 2009;6(3):036001.View ArticlePubMedGoogle Scholar
- Honda H, Tanemura M, Yoshida A. Differentiation of wing epidermal scale cells in a butterfly under the lateral inhibition model–appearance of large cells in a polygonal pattern. Acta Biotheor. 2000;48(2):121–36.View ArticlePubMedGoogle Scholar
- Pitt-Francis J, Pathmanathan P, Bernabeu MO, Bordas R, Cooper J, Fletcher AG, et al. Chaste: A test-driven approach to software development for biological modelling. Comput Phys Commun. 2009;180(12):2452–71.View ArticleGoogle Scholar
- Resasco DC, Gao F, Morgan F, Novak IL, Schaff JC, Slepchenko BM. Virtual Cell: computational tools for modeling in cell biology. Wiley Interdiscip Rev Syst Biol and Med. 2012;4(2):129–40.View ArticleGoogle Scholar
- Luke S, Cioffi-Revilla C, Panait L, Sullivan K, Balan G. MASON: A Multiagent Simulation Environment. Simulation. 2005;81(7):517–27.View ArticleGoogle Scholar
- Schroeder W, Martin K, Lorensen B. The Visualization Toolkit: An Object-Oriented Approach to 3D Graphics. 4th ed. Clifton Park: New York, Kitware, Inc.; 2006.Google Scholar
- A Cell-centered, Agent-based Framework with flexible Granularity, Available from: https://simtk.org/home/dvframework.
- Railsback SF, Lytinen SL, Jackson SK. Agent-based Simulation Platforms: Review and Development Recommendations. Simulation. 2006;82(9):609–23.View ArticleGoogle Scholar
- Engelberg JA, Datta A, Mostov KE, Hunt CA. MDCK Cystogenesis Driven by Cell Stabilization within Computational Analogues. PLoS Comput Biol. 2011;7(4), e1002030.View ArticlePubMedPubMed CentralGoogle Scholar
- Engelberg JA, Datta A, Mostov KE, Hunt CA. A Robust in Silico Analogue of MDCK Cystogenesis Mimics Growth in Multiple Culture Conditions. In: Proceedings of the 2011 Workshop on Agent-Directed Simulation. ADS ’11. San Diego, CA, USA: Society for Computer Simulation International; 2011. p. 52–6.Google Scholar
- Kim SH, Yu W, Mostov K, Matthay MA, Hunt CA. A computational approach to understand in vitro alveolar morphogenesis. PLoS One. 2009;4(3), e4819.View ArticlePubMedPubMed CentralGoogle Scholar