No abstract

Vascular development is a complex process regulated by dynamic biological networks that vary in topology and state across different tissues and developmental stages. Signals regulating de novo blood vessel formation (vasculogenesis) and remodeling (angiogenesis) come from a variety of biological pathways linked to endothelial cell (EC) behavior, extracellular matrix (ECM) remodeling and the local generation of chemokines and growth factors. Simulating these interactions at a systems level requires sufficient biological detail about the relevant molecular pathways and associated cellular behaviors, and tractable computational models that offset mathematical and biological complexity. Here, we describe a novel multicellular agent-based model of vasculogenesis using the CompuCell3D (http://www.compucell3d.org/) modeling environment supplemented with semi-automatic knowledgebase creation. The model incorporates vascular endothelial growth factor signals, pro- and anti-angiogenic inflammatory chemokine signals, and the plasminogen activating system of enzymes and proteases linked to ECM interactions, to simulate nascent EC organization, growth and remodeling. The model was shown to recapitulate stereotypical capillary plexus formation and structural emergence of non-coded cellular behaviors, such as a heterologous bridging phenomenon linking endothelial tip cells together during formation of polygonal endothelial cords. Molecular targets in the computational model were mapped to signatures of vascular disruption derived from in vitro chemical profiling using the EPA's ToxCast high-throughput screening (HTS) dataset. Simulating the HTS data with the cell-agent based model of vascular development predicted adverse effects of a reference anti-angiogenic thalidomide analog, 5HPP-33, on in vitro angiogenesis with respect to both concentration-response and morphological consequences. These findings support the utility of cell agent-based models for simulating a morphogenetic series of events and for the first time demonstrate the applicability of these models for predictive toxicology.

Morphogenetic events are driven by cell-generated physical forces and complex cellular dynamics. To improve our capacity to predict developmental effects from chemical-induced cellular alterations, we built a multicellular agent-based model in CompuCell3D that recapitulates the cellular networks and collective cell behavior underlying growth and fusion of the mammalian secondary palate. The model incorporated multiple signaling pathways (TGFβ, BMP, FGF, EGF, and SHH) in a biological framework to recapitulate morphogenetic events from palatal outgrowth through midline fusion. It effectively simulated higher-level phenotypes (e.g., midline contact, medial edge seam (MES) breakdown, mesenchymal confluence, and fusion defects) in response to genetic or environmental perturbations. Perturbation analysis of various control features revealed model functionality with respect to cell signaling systems and feedback loops for growth and fusion, diverse individual cell behaviors and collective cellular behavior leading to physical contact and midline fusion, and quantitative analysis of the TGF/EGF switch that controls MES breakdown-a key event in morphogenetic fusion. The virtual palate model was then executed with theoretical chemical perturbation scenarios to simulate switch behavior leading to a disruption of fusion following chronic (e.g., dioxin) and acute (e.g., retinoic acid) chemical exposures. This computer model adds to similar systems models toward an integrative "virtual embryo" for simulation and quantitative prediction of adverse developmental outcomes following genetic perturbation and/or environmental disruption.