Computational Model of Secondary Palate Fusion and Disruption

A computational dynamic systems model for in silico prediction of neural tube closure defects

Neural tube closure is a critical morphogenetic event during early vertebrate development. This complex process is susceptible to perturbation by genetic errors and chemical disruption, which can induce severe neural tube defects (NTDs) such as spina bifida. We built a computational agent-based model (ABM) of neural tube development based on the known biology of morphogenetic signals and cellular biomechanics underlying neural fold elevation, bending and fusion. The computer model functionalizes cell signals and responses to render a dynamic representation of neural tube closure. Perturbations in the control network can then be introduced synthetically or from biological data to yield quantitative simulation and probabilistic prediction of NTDs by incidence and degree of defect. Translational applications of the model include mechanistic understanding of how singular or combinatorial alterations in gene-environmental interactions and animal-free assessment of developmental toxicity for an important human birth defect (spina bifida) and potentially other neurological problems linked to development of the brain and spinal cord.

Engineering a computable epiblast for in silico modeling of developmental toxicity

Developmental hazard evaluation is an important part of assessing chemical risks during pregnancy. Toxicological outcomes from prenatal testing in pregnant animals result from complex chemical-biological interactions, and while New Approach Methods (NAMs) based on in vitro bioactivity profiles of human cells offer promising alternatives to animal testing, most of these assays lack cellular positional information, physical constraints, and regional organization of the intact embryo. Here, we engineered a fully computable model of the embryonic disc in the CompuCell3D.org modeling environment to simulate epithelial-mesenchymal transition (EMT) of epiblast cells and self-organization of mesodermal domains (chordamesoderm, paraxial, lateral plate, posterior/extraembryonic). Mesodermal fate is modeled by synthetic activity of the BMP4-NODAL-WNT signaling axis. Cell position in the epiblast determines timing with respect to EMT for 988 computational cells in the computer model. An autonomous homeobox (Hox) clock hidden in the epiblast is driven by WNT-FGF4-CDX signaling. Executing the model renders a quantitative cell-level computation of mesodermal fate and consequences of perturbation based on known biology. For example, synthetic perturbation of the control network rendered altered phenotypes (cybermorphs) mirroring some aspects of experimental mouse embryology, with electronic knockouts, under-activation (hypermorphs) or over-activation (hypermorphs) particularly affecting the size and specification of the posterior mesoderm. This foundational model is trained on embryology but capable of performing a wide variety of toxicological tasks conversing through anatomical simulation to integrate in vitro chemical bioactivity data with known embryology. It is amenable to quantitative simulation for probabilistic prediction of early developmental toxicity.

Design and mathematical analysis of activating transcriptional amplifiers that enable modular temporal control in synthetic juxtacrine circuits

The ability to control mammalian cells such that they self-organize or enact therapeutic effects as desired has incredible implications. Not only would it further our understanding of native processes such as development and the immune response, but it would also have powerful applications in medical fields such as regenerative medicine and immunotherapy. This control is typically obtained by synthetic circuits that use synthetic receptors, but control remains incomplete. The synthetic juxtacrine receptors (SJRs) are widely used as they are fully modular and enable spatial control, but they have limited gene expression amplification and temporal control. As these are integral facets to cell control, I therefore designed transcription factor based amplifiers that amplify gene expression and enable unidirectional temporal control by prolonging duration of target gene expression. Using a validated in silico framework for SJR signaling, I combined these amplifiers with SJRs and show that these SJR amplifier circuits can direct spatiotemporal patterning and improve the quality of self-organization. I then show that these circuits can improve chimeric antigen receptor (CAR) T cell tumor killing against various heterogenous antigen expression tumors. These amplifiers are flexible tools that improve control over SJR based circuits with both basic and therapeutic applications.

An improved multicellular human organoid model for the study of chemical effects on palatal fusion

BACKGROUND

Tissue fusion is a mechanism involved in the development of the heart, iris, genital tubercle, neural tube, and palate during embryogenesis. Failed fusion of the palatal shelves could result in cleft palate (CP), a common birth defect. Organotypic models constructed of human cells offer an opportunity to investigate developmental processes in the human. Previously, our laboratory developed an organoid model of the human palate that contains human mesenchyme and epithelial progenitor cells to study the effects of chemicals on fusion.

METHODS

Here, we developed an organoid model more representative of the embryonic palate that includes three cell types: mesenchyme, endothelial, and epithelial cells. We measured fusion by a decrease in epithelial cells at the contact point between the organoids and compared the effects of CP teratogens on fusion and toxicity in the previous and current organoid models. We further tested additional suspect teratogens in our new model.

RESULTS

We found that the three-cell-type model is more sensitive to fusion inhibition by valproic acid and inhibitors of FGF, BMP, and TGFβRI/II. In this new model, we tested other suspect CP teratogens and found that nocodazole, topiramate, and Y27632 inhibit fusion at concentrations that do not induce toxicity.

CONCLUSION

This sensitive human three-cell-type organotypic model accurately evaluates chemicals for cleft palate teratogenicity.

Human-Based New Approach Methodologies in Developmental Toxicity Testing: A Step Ahead from the State of the Art with a Feto-Placental Organ-on-Chip Platform

Developmental toxicity testing urgently requires the implementation of human-relevant new approach methodologies (NAMs) that better recapitulate the peculiar nature of human physiology during pregnancy, especially the placenta and the maternal/fetal interface, which represent a key stage for human lifelong health. Fit-for-purpose NAMs for the placental-fetal interface are desirable to improve the biological knowledge of environmental exposure at the molecular level and to reduce the high cost, time and ethical impact of animal studies. This article reviews the state of the art on the available in vitro (placental, fetal and amniotic cell-based systems) and in silico NAMs of human relevance for developmental toxicity testing purposes; in addition, we considered available Adverse Outcome Pathways related to developmental toxicity. The OECD TG 414 for the identification and assessment of deleterious effects of prenatal exposure to chemicals on developing organisms will be discussed to delineate the regulatory context and to better debate what is missing and needed in the context of the Developmental Origins of Health and Disease hypothesis to significantly improve this sector. Starting from this analysis, the development of a novel human feto-placental organ-on-chip platform will be introduced as an innovative future alternative tool for developmental toxicity testing, considering possible implementation and validation strategies to overcome the limitation of the current animal studies and NAMs available in regulatory toxicology and in the biomedical field.

Computational model for fetal skeletal defects potentially linked to disruption of retinoic acid signaling

All-trans retinoic acid (ATRA) gradients determine skeletal patterning morphogenesis and can be disrupted by diverse genetic or environmental factors during pregnancy, leading to fetal skeleton defects. Adverse Outcome Pathway (AOP) frameworks for ATRA metabolism, signaling, and homeostasis allow for the development of new approach methods (NAMs) for predictive toxicology with less reliance on animal testing. Here, a data-driven model was constructed to identify chemicals associated with both ATRA pathway bioactivity and prenatal skeletal defects. The phenotype data was culled from ToxRefDB prenatal developmental toxicity studies and produced a list of 363 ToxRefDB chemicals with altered skeletal observations. Defects were classified regionally as cranial, post-cranial axial, appendicular, and other (unspecified) features based on ToxRefDB descriptors. To build a multivariate statistical model, high-throughput screening bioactivity data from >8,070 chemicals in ToxCast/Tox21 across 10 in vitro assays relevant to the retinoid signaling system were evaluated and compared to literature-based candidate reference chemicals in the dataset. There were 48 chemicals identified for effects on both in vivo skeletal defects and in vitro ATRA pathway targets for computational modeling. The list included 28 chemicals with prior evidence of skeletal defects linked to retinoid toxicity and 20 chemicals without prior evidence. The combination of thoracic cage defects and DR5 (direct repeats of 5 nucleotides for RAR/RXR transactivation) disruption was the most frequently occurring phenotypic and target disturbance, respectively. This data model provides valuable AOP elucidation and validates current mechanistic understanding. These findings also shed light on potential avenues for new mechanistic discoveries related to ATRA pathway disruption and associated skeletal dysmorphogenesis due to environmental exposures.

Pluripotent stem cell assays: Modalities and applications for predictive developmental toxicity

This manuscript provides a review focused on embryonic stem cell-based models and their place within the landscape of alternative developmental toxicity assays. Against the background of the principles of developmental toxicology, the wide diversity of alternative methods using pluripotent stem cells developed in this area over the past half century is reviewed. In order to provide an overview of available models, a systematic scoping review was conducted following a published protocol with inclusion criteria, which were applied to select the assays. Critical aspects including biological domain, readout endpoint, availability of standardized protocols, chemical domain, reproducibility and predictive power of each assay are described in detail, in order to review the applicability and limitations of the platform in general and progress moving forward to implementation. The horizon of innovative routes of promoting regulatory implementation of alternative methods is scanned, and recommendations for further work are given.

Parameterized Computational Framework for the Description and Design of Genetic Circuits of Morphogenesis Based on Contact-Dependent Signaling and Changes in Cell-Cell Adhesion

Synthetic development is a nascent field of research that uses the tools of synthetic biology to design genetic programs directing cellular patterning and morphogenesis in higher eukaryotic cells, such as mammalian cells. One specific example of such synthetic genetic programs was based on cell-cell contact-dependent signaling using synthetic Notch pathways and was shown to drive the formation of multilayered spheroids by modulating cell-cell adhesion via differential expression of cadherin family proteins in a mouse fibroblast cell line (L929). The design method for these genetic programs relied on trial and error, which limited the number of possible circuits and parameter ranges that could be explored. Here, we build a parameterized computational framework that, given a cell-cell communication network driving changes in cell adhesion and initial conditions as inputs, predicts developmental trajectories. We first built a general computational framework where contact-dependent cell-cell signaling networks and changes in cell-cell adhesion could be designed in a modular fashion. We then used a set of available in vitro results (that we call the "training set" in analogy to similar pipelines in the machine learning field) to parameterize the computational model with values for adhesion and signaling. We then show that this parameterized model can qualitatively predict experimental results from a "testing set" of available in vitro data that varied the genetic network in terms of adhesion combinations, initial number of cells, and even changes to the network architecture. Finally, this parameterized model is used to recommend novel network implementation for the formation of a four-layered structure that has not been reported previously. The framework that we develop here could function as a testing ground to identify the reachable space of morphologies that can be obtained by controlling contact-dependent cell-cell communications and adhesion with these molecular tools and in this cellular system. Additionally, we discuss how the model could be expanded to include other forms of communication or effectors for the computational design of the next generation of synthetic developmental trajectories.

Dose Addition in the Induction of Craniofacial Malformations in Zebrafish Embryos Exposed to a Complex Mixture of Food-Relevant Chemicals with Dissimilar Modes of Action

BACKGROUND

Humans are exposed to combinations of chemicals. In cumulative risk assessment (CRA), regulatory bodies such as the European Food Safety Authority consider dose addition as a default and sufficiently conservative approach. The principle of dose addition was confirmed previously for inducing craniofacial malformations in zebrafish embryos in binary mixtures of chemicals with either similar or dissimilar modes of action (MOAs).

OBJECTIVES

In this study, we explored a workflow to select and experimentally test multiple compounds as a complex mixture with each of the compounds at or below its no observed adverse effect level (NOAEL), in the same zebrafish embryo model.

METHODS

Selection of candidate compounds that potentially induce craniofacial malformations was done using in silico methods-structural similarity, molecular docking, and quantitative structure-activity relationships-applied to a database of chemicals relevant for oral exposure in humans via food (EuroMix inventory, n = 1,598 ). A final subselection was made manually to represent different regulatory fields (e.g., food additives, industrial chemicals, plant protection products), different chemical families, and different MOAs.

RESULTS

A final selection of eight compounds was examined in the zebrafish embryo model, and craniofacial malformations were observed in embryos exposed to each of the compounds, thus confirming the developmental toxicity as predicted by the in silico methods. When exposed to a mixture of the eight compounds, each at its NOAEL, substantial craniofacial malformations were observed; according to a dose-response analysis, even embryos exposed to a 7-fold dilution of this mixture still exhibited a slight abnormal phenotype. The cumulative effect of the compounds in the mixture was in accordance with dose addition (added doses of the individual compounds after adjustment for relative potencies), despite different MOAs of the compounds involved.

DISCUSSION

This case study of a complex mixture inducing craniofacial malformations in zebrafish embryos shows that dose addition can adequately predicted the cumulative effect of a mixture of multiple substances at low doses, irrespective of the (expected) MOA. The applied workflow may be useful as an approach for CRA in general. https://doi.org/10.1289/EHP9888.

Histological and Immunohistochemical Studies to Determine the Mechanism of Cleft Palate Induction after Palatal Fusion in Mice Exposed to TCDD

Rupture of the basement membrane in fused palate tissue can cause the palate to separate after fusion in mice, leading to the development of cleft palate. Here, we further elucidate the mechanism of palatal separation after palatal fusion in 8–10-week-old ICR female mice. On day 12 of gestation, 40 μg/kg of 2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD), sufficient to cause cleft palate in 100% of mice, was dissolved in 0.4 mL of olive oil containing toluene and administered as a single dose via a gastric tube. Fetal palatine frontal sections were observed by H&E staining, and epithelial cell adhesion factors, apoptosis, and cell proliferation were observed from the anterior to posterior palate. TUNEL-positive cells and Ki67-positive cells were observed around the posterior palatal dissection area of the TCDD-treated group. Moreover, in fetal mice exposed to TCDD, some fetuses exhibited cleft palate dehiscence during fusion. The results suggest that palatal dehiscence may be caused by abnormal cell proliferation in epithelial tissues, decreased intercellular adhesion, and inhibition of mesenchymal cell proliferation. By elucidating the mechanism of cleavage after palatal fusion, this research can contribute to establishing methods for the prevention of cleft palate development.

Endoderm and mesoderm derivatives in embryonic stem cell differentiation and their use in developmental toxicity testing

Embryonic stem cell differentiation models have increasingly been applied in non-animal test systems for developmental toxicity. After the initial focus on cardiac differentiation, attention has also included an array of neuro-ectodermal differentiation routes. Alternative differentiation routes in the mesodermal and endodermal germ lines have received less attention. This review provides an inventory of achievements in the latter areas of embryonic stem cell differentiation, with a view to possibilities for their use in non-animal test systems in developmental toxicology. This includes murine and human stem cell differentiation models, and also gains information from the field of stem cell use in regenerative medicine. Endodermal stem cell derivatives produced in vitro include hepatocytes, pancreatic cells, lung epithelium, and intestinal epithelium, and mesodermal derivatives include cardiac muscle, osteogenic, vascular and hemopoietic cells. This inventory provides an overview of studies on the different cell types together with biomarkers and culture conditions that stimulate these differentiation routes from embryonic stem cells. These models may be used to expand the spectrum of embryonic stem cell based new approach methodologies in non-animal developmental toxicity testing.

An overview of current practices for regulatory risk assessment with lessons learnt from cosmetics in the European Union

Risk assessments of various types of chemical compounds are carried out in the European Union (EU) foremost to comply with legislation and to support regulatory decision-making with respect to their safety. Historically, risk assessment has relied heavily on animal experiments. However, the EU is committed to reduce animal experimentation and has implemented several legislative changes, which have triggered a paradigm shift towards human-relevant animal-free testing in the field of toxicology, in particular for risk assessment. For some specific endpoints, such as skin corrosion and irritation, validated alternatives are available whilst for other endpoints, including repeated dose systemic toxicity, the use of animal data is still central to meet the information requirements stipulated in the different legislations. The present review aims to provide an overview of established and more recently introduced methods for hazard assessment and risk characterisation for human health, in particular in the context of the EU Cosmetics Regulation (EC No 1223/2009) as well as the Registration, Evaluation, Authorisation and Restriction of Chemicals (REACH) Regulation (EC 1907/2006).

FutureTox IV Workshop Summary: Predictive Toxicology for Healthy Children

FutureTox IV, a Society of Toxicology Contemporary Concepts in Toxicology workshop, was held in November 2018. Building upon FutureTox I, II, and III, this conference focused on the latest science and technology for in vitro profiling and in silico modeling as it relates to predictive developmental and reproductive toxicity (DART). Publicly available high-throughput screening data sets are now available for broad in vitro profiling of bioactivities across large inventories of chemicals. Coupling this vast amount of mechanistic data with a deeper understanding of molecular embryology and post-natal development lays the groundwork for using new approach methodologies (NAMs) to evaluate chemical toxicity, drug efficacy, and safety assessment for embryo-fetal development. NAM is a term recently adopted in reference to any technology, methodology, approach, or combination thereof that can be used to provide information on chemical hazard and risk assessment to avoid the use of intact animals (U.S. Environmental Protection Agency [EPA], Strategic plan to promote the development and implementation of alternative test methods within the tsca program, 2018, https://www.epa.gov/sites/production/files/2018-06/documents/epa_alt_strat_plan_6-20-18_clean_final.pdf). There are challenges to implementing NAMs to evaluate chemicals for developmental toxicity compared with adult toxicity. This forum article reviews the 2018 workshop activities, highlighting challenges and opportunities for applying NAMs for adverse pregnancy outcomes (eg, preterm labor, malformations, low birth weight) as well as disorders manifesting postnatally (eg, neurodevelopmental impairment, breast cancer, cardiovascular disease, fertility). DART is an important concern for different regulatory statutes and test guidelines. Leveraging advancements in such approaches and the accompanying efficiencies to detecting potential hazards to human development are the unifying concepts toward implementing NAMs in DART testing. Although use of NAMs for higher level regulatory decision making is still on the horizon, the conference highlighted novel testing platforms and computational models that cover multiple levels of biological organization, with the unique temporal dynamics of embryonic development, and novel approaches for estimating toxicokinetic parameters essential in supporting in vitro to in vivo extrapolation.

A Microphysiological Approach to Evaluate Effectors of Intercellular Hedgehog Signaling in Development

Paracrine signaling in the tissue microenvironment is a central mediator of morphogenesis, and modeling this dynamic intercellular activity in vitro is critical to understanding normal and abnormal development. For example, Sonic Hedgehog (Shh) signaling is a conserved mechanism involved in multiple developmental processes and strongly linked to human birth defects including orofacial clefts of the lip and palate. SHH ligand produced, processed, and secreted from the epithelial ectoderm is shuttled through the extracellular matrix where it binds mesenchymal receptors, establishing a gradient of transcriptional response that drives orofacial morphogenesis. In humans, complex interactions of genetic predispositions and environmental insults acting on diverse molecular targets are thought to underlie orofacial cleft etiology. Consequently, there is a need for tractable in vitro approaches that model this complex cellular and environmental interplay and are sensitive to disruption across the multistep signaling cascade. We developed a microplate-based device that supports an epithelium directly overlaid onto an extracellular matrix-embedded mesenchyme, mimicking the basic tissue architecture of developing orofacial tissues. SHH ligand produced from the epithelium generated a gradient of SHH-driven transcription in the adjacent mesenchyme, recapitulating the gradient of pathway activity observed in vivo. Shh pathway activation was antagonized by small molecule inhibitors of epithelial secretory, extracellular matrix transport, and mesenchymal sensing targets, supporting the use of this approach in high-content chemical screening of the complete Shh pathway. Together, these findings demonstrate a novel and practical microphysiological model with broad utility for investigating epithelial-mesenchymal interactions and environmental signaling disruptions in development.

Meeting report on the NIDDK/AUA Workshop on Congenital Anomalies of External Genitalia: challenges and opportunities for translational research

Congenital anomalies of the external genitalia (CAEG) are a prevalent and serious public health concern with lifelong impacts on the urinary function, sexual health, fertility, tumor development, and psychosocial wellbeing of affected individuals. Complications of treatment are frequent, and data reflecting long-term outcomes in adulthood are limited. To identify a path forward to improve treatments and realize the possibility of preventing CAEG, the National Institute of Diabetes and Digestive and Kidney Diseases and the American Urological Association convened researchers from a range of disciplines to coordinate research efforts to fully understand the different etiologies of these common conditions, subsequent variation in clinical phenotypes, and best practices for long term surgical success. Meeting participants concluded that a central data hub for clinical evaluations, including collection of DNA samples from patients and their parents, and short interviews to determine familial penetrance (small pedigrees), would accelerate research in this field. Such a centralized datahub will advance efforts to develop detailed multi-dimensional phenotyping and will enable access to genome sequence analyses and associated metadata to define the genetic bases for these conditions. Inclusion of tissue samples and integration of clinical studies with basic research using human cells and animal models will advance efforts to identify the developmental mechanisms that are disrupted during development and will add cellular and molecular granularity to phenotyping CAEG. While the discussion focuses heavily on hypospadias, this can be seen as a potential template for other conditions in the realm of CAEG, including cryptorchidism or the exstrophy-epispadias complex. Taken together with long-term clinical follow-up, these data could inform surgical choices and improve likelihood for long-term success.

Computational Biology and in silico Toxicodynamics

New approach methodologies (NAMs) refer to any non-animal technology, methodology, approach, or combination thereof that can be used to provide information on chemical hazard and risk assessment that avoids the use of intact animals. A spectrum of in silico models is needed for the integrated analysis of various domains in toxicology to improve predictivity and reduce animal testing. This review focuses on in silico approaches, computer models, and computational intelligence for developmental and reproductive toxicity (predictive DART), providing a means to measure toxicodynamics in simulated systems for quantitative prediction of adverse outcomes phenotypes.

Culture Conditions Affect Chemical-Induced Developmental Toxicity In Vitro: The Case of Folic Acid, Methionine and Methotrexate in the Neural Embryonic Stem Cell Test

In vitro tests are increasingly applied in chemical hazard assessment. Basic culture conditions may affect the outcome of in vitro tests and should be optimised to reduce false predictions. The neural embryonic stem cell test (ESTn) can predict early neurodevelopmental effects of chemicals, as it mimics the differentiation of stem cells towards the neuroectodermal lineage. Normal early neural differentiation depends crucially on folic acid (FA) and methionine (MET), both elements of the one-carbon (1C) cycle. The aim of this study was to assess the concentration-dependent influence of FA and MET on neural differentiation in the ESTn, and its consequences for assay sensitivity to methotrexate (MTX), a compound that interferes with the 1C cycle. Neural differentiation was inhibited below 0.007 mM and above 0.22 mM FA, while both stem cell viability (< 0.097 mM, > 1.52 mM) and neural differentiation (< 0.181 mM, > 1.35 mM) were affected when changing MET concentrations. A 10-day exposure to 13 nM MTX inhibited neural differentiation, especially in FA- and MET-deficient conditions. However, a 24-hour exposure to 39 nM MTX decreased neural cell and neural crest cell differentiation markers only when the concentration of FA in the medium was three times the standard concentration, which was expected to have a protective effect against MTX. These results show the importance of nutrient concentrations, exposure scenarios and timing of read-outs for cell differentiation and compound sensitivity in the ESTn. Caution should be taken when interpreting results from a single in vitro test, especially when extrapolating to effects on complex morphogenetic processes, like neural tube development.

An ontology for developmental processes and toxicities of neural tube closure

In recent years, the development and implementation of animal-free approaches to chemical and pharmaceutical hazard and risk assessment has taken off. Alternative approaches are being developed starting from the perspective of human biology and physiology. Neural tube closure is a vital step that occurs early in human development. Correct closure of the neural tube depends on a complex interplay between proteins along a number of protein concentration gradients. The sensitivity of neural tube closure to chemical disturbance of signalling pathways such as the retinoid pathway, is well known. To map the pathways underlying neural tube closure, literature data on the molecular regulation of neural tube closure were collected. As the process of neural tube closure is highly conserved in vertebrates, the extensive literature available for the mouse was used whilst considering its relevance for humans. Thus, important cell compartments, regulatory pathways, and protein interactions essential for neural tube closure under physiological circumstances were identified and mapped. An understanding of aberrant processes leading to neural tube defects (NTDs) requires detailed maps of neural tube embryology, including the complex genetic signals and responses underlying critical cellular dynamical and biomechanical processes. The retinoid signaling pathway serves as a case study for this ontology because of well-defined crosstalk with the genetic control of neural tube patterning and morphogenesis. It is a known target for mechanistically-diverse chemical structures that disrupt neural tube closure The data presented in this manuscript will set the stage for constructing mathematical models and computer simulation of neural tube closure for human-relevant AOPs and predictive toxicology.

Potential frameworks to support evaluation of mechanistic data for developmental neurotoxicity outcomes: A symposium report

A key challenge in systematically incorporating mechanistic data into human health assessments is that, compared to studies of apical health endpoints, these data are both more abundant (mechanistic studies routinely outnumber other studies by several orders of magnitude) and more heterogeneous (e.g. different species, test system, tissue, cell type, exposure paradigm, or specific assays performed). A structured decision-making process for organizing, integrating, and weighing mechanistic DNT data for use in human health risk assessments will improve the consistency and efficiency of such evaluations. At the Developmental Neurotoxicology Society (DNTS) 2016 annual meeting, a symposium was held to address the application of existing organizing principles and frameworks for evaluation of mechanistic data relevant to interpreting neurotoxicology data. Speakers identified considerations with potential to advance the use of mechanistic DNT data in risk assessment, including considering the context of each exposure, since epigenetics, tissue type, sex, stress, nutrition and other factors can modify toxicity responses in organisms. It was also suggested that, because behavior is a manifestation of complex nervous system function, the presence and absence of behavioral change itself could be used to organize the interpretation of multiple complex simultaneous mechanistic changes. Several challenges were identified with frameworks and their implementation, and ongoing research to develop these approaches represents an early step toward full evaluation of mechanistic DNT data for assessments.

Predicting the safety of medicines in pregnancy: A workshop report

The framework for developmental toxicity testing has remained largely unchanged for over 50 years and although it remains invaluable in assessing potential risks in pregnancy, knowledge gaps exist, and some outcomes do not necessarily correlate with clinical experience. Advances in omics, in silico approaches and alternative assays are providing opportunities to enhance our understanding of embryo-fetal development and the prediction of potential risks associated with the use of medicines in pregnancy. A workshop organised by the Medicines and Healthcare products Regulatory Agency (MHRA), "Predicting the Safety of Medicines in Pregnancy - a New Era?", was attended by delegates representing regulatory authorities, academia, industry, patients, funding bodies and software developers to consider how to improve the quality of and access to nonclinical developmental toxicity data and how to use this data to better predict the safety of medicines in human pregnancy. The workshop delegates concluded that based on comparative data to date alternative methodologies are currently no more predictive than conventional methods and not qualified for use in regulatory submissions. To advance the development and qualification of alternative methodologies, there is a requirement for better coordinated multidisciplinary cross-sector interactions coupled with data sharing. Furthermore, a better understanding of human developmental biology and the incorporation of this knowledge into the development of alternative methodologies is essential to enhance the prediction of adverse outcomes for human development. The output of the workshop was a series of recommendations aimed at supporting multidisciplinary efforts to develop and validate these alternative methodologies.

A multiscale in silico model of endothelial to mesenchymal transformation in a tumor microenvironment

Endothelial to mesenchymal transformation (EndMT) is a process in which endothelial cells gain a mesenchymal-like phenotype in response to mechanobiological signals that results in the remodeling or repair of underlying tissue. While initially associated with embryonic development, this process has since been shown to occur in adult tissue remodeling including wound healing, fibrosis, and cancer. In an attempt to understand the role of EndMT in cancer progression and metastasis, we present a multiscale, three-dimensional, in silico model. The model couples tissue level phenomena such as extracellular matrix remodeling, cellular level phenomena such as migration and proliferation, and chemical transport in the tumor microenvironment to mimic in vitro tissue models of the cancer microenvironment. The model is used to study the presence of EndMT-derived activated fibroblasts (EDAFs) and varying substrate stiffness on tumor cell migration and proliferation. The simulations accurately model the behavior of tumor cells under given conditions. The presence of EDAFs and/or an increase in substrate stiffness resulted in an increase in tumor cell activity. This model lays the foundation of further studies of EDAFs in a tumor microenvironment on a cellular and subcellular physiological level.

Synthetic development: building mammalian multicellular structures with artificial genetic programs

Synthetic biology efforts began in simple single-cell systems, which were relatively easy to manipulate genetically (Cameron et al., 2014). The field grew exponentially in the last two decades, and one of the latest frontiers are synthetic developmental programs for multicellular mammalian systems (Black et al., 2017; Wieland and Fussenegger, 2012) to genetically control features such as patterning or morphogenesis. These programs rely on engineered cell-cell communications, multicellular gene regulatory networks and effector genes. Here, we contextualize the first of these synthetic developmental programs, examine molecular and computational tools that can be used to generate next generation versions, and present the general logic that underpins these approaches. These advances are exciting as they represent a novel way to address both control and understanding in the field of developmental biology and tissue development (Elowitz and Lim, 2010; Velazquez et al., 2018; White et al., 2018; Morsut, 2017). This field is just at the beginning, and it promises to be of major interest in the upcoming years of biomedical research.

Characterizing cleft palate toxicants using ToxCast data, chemical structure, and the biomedical literature

Cleft palate has been linked to both genetic and environmental factors that perturb key events during palatal morphogenesis. As a developmental outcome, it presents a challenging, mechanistically complex endpoint for predictive modeling. A data set of 500 chemicals evaluated for their ability to induce cleft palate in animal prenatal developmental studies was compiled from Toxicity Reference Database and the biomedical literature, which included 63 cleft palate active and 437 inactive chemicals. To characterize the potential molecular targets for chemical-induced cleft palate, we mined the ToxCast high-throughput screening database for patterns and linkages in bioactivity profiles and chemical structural descriptors. ToxCast assay results were filtered for cytotoxicity and grouped by target gene activity to produce a "gene score." Following unsuccessful attempts to derive a global prediction model using structural and gene score descriptors, hierarchical clustering was applied to the set of 63 cleft palate positives to extract local structure-bioactivity clusters for follow-up study. Patterns of enrichment were confirmed on the complete data set, that is, including cleft palate inactives, and putative molecular initiating events identified. The clusters corresponded to ToxCast assays for cytochrome P450s, G-protein coupled receptors, retinoic acid receptors, the glucocorticoid receptor, and tyrosine kinases/phosphatases. These patterns and linkages were organized into preliminary decision trees and the resulting inferences were mapped to a putative adverse outcome pathway framework for cleft palate supported by literature evidence of current mechanistic understanding. This general data-driven approach offers a promising avenue for mining chemical-bioassay drivers of complex developmental endpoints where data are often limited.

The Next Generation Blueprint of Computational Toxicology at the U.S. Environmental Protection Agency

The U.S. Environmental Protection Agency (EPA) is faced with the challenge of efficiently and credibly evaluating chemical safety often with limited or no available toxicity data. The expanding number of chemicals found in commerce and the environment, coupled with time and resource requirements for traditional toxicity testing and exposure characterization, continue to underscore the need for new approaches. In 2005, EPA charted a new course to address this challenge by embracing computational toxicology (CompTox) and investing in the technologies and capabilities to push the field forward. The return on this investment has been demonstrated through results and applications across a range of human and environmental health problems, as well as initial application to regulatory decision-making within programs such as the EPA's Endocrine Disruptor Screening Program. The CompTox initiative at EPA is more than a decade old. This manuscript presents a blueprint to guide the strategic and operational direction over the next 5 years. The primary goal is to obtain broader acceptance of the CompTox approaches for application to higher tier regulatory decisions, such as chemical assessments. To achieve this goal, the blueprint expands and refines the use of high-throughput and computational modeling approaches to transform the components in chemical risk assessment, while systematically addressing key challenges that have hindered progress. In addition, the blueprint outlines additional investments in cross-cutting efforts to characterize uncertainty and variability, develop software and information technology tools, provide outreach and training, and establish scientific confidence for application to different public health and environmental regulatory decisions.

Development of an organotypic stem cell model for the study of human embryonic palatal fusion

BACKGROUND

Cleft palate (CP) is a common birth defect, occurring in an estimated 1 in 1,000 births worldwide. The secondary palate is formed by paired palatal shelves, consisting of a mesenchymal core with an outer layer of epithelial cells that grow toward each other, attach, and fuse. One of the mechanisms that can cause CP is failure of fusion, that is, failure to remove the epithelial seam between the palatal shelves to allow the mesenchyme confluence. Epidermal growth factor (EGF) plays an important role in palate growth and differentiation, while it may impede fusion.

METHODS

We developed a 3D organotypic model using human mesenchymal and epithelial stem cells to mimic human embryonic palatal shelves, and tested the effects of human EGF (hEGF) on proliferation and fusion. Spheroids were generated from human umbilical-derived mesenchymal stem cells (hMSCs) directed down an osteogenic lineage. Heterotypic spheroids, or organoids, were constructed by coating hMSC spheroids with extracellular matrix solution followed by a layer of human progenitor epithelial keratinocytes (hPEKs). Organoids were incubated in co-culture medium with or without hEGF and assessed for cell proliferation and time to fusion.

RESULTS

Osteogenic differentiation in hMSC spheroids was highest by Day 13. hEGF delayed fusion of organoids after 12 and 18 hr of contact. hEGF increased proliferation in organoids at 4 ng/ml, and proliferation was detected in hPEKs alone.

CONCLUSION

Our results show that this model of human palatal fusion appropriately mimics the morphology of the developing human palate and responds to hEGF as expected.

Possibilities and limitations of three-dimensional reconstruction and simulation techniques to identify patterns, rhythms and functions of apoptosis in the early developing neural tube

The now classical idea that programmed cell death (apoptosis) contributes to a plethora of developmental processes still has lost nothing of its impact. It is, therefore, important to establish effective three-dimensional (3D) reconstruction as well as simulation techniques to decipher the exact patterns and functions of such apoptotic events. The present study focuses on the question whether and how apoptosis promotes neurulation-associated processes in the spinal cord of Tupaia belangeri (Tupaiidae, Scandentia, Mammalia). Our 3D reconstructions demonstrate that at least two craniocaudal waves of apoptosis consecutively pass through the dorsal spinal cord. The first wave appears to be involved in neural fold fusion and/or in selection processes among premigratory neural crest cells. The second one seems to assist in establishing the dorsal signaling center known as the roof plate. In the hindbrain, in contrast, apoptosis among premigratory neural crest cells progresses craniocaudally but discontinuously, in a segment-specific manner. Unlike apoptosis in the spinal cord, these segment-specific apoptotic events, however, precede later ones that seemingly support neural fold fusion and/or postfusion remodeling. Arguing with Whitehead that biological patterns and rhythms differ in that biological rhythms depend "upon the differences involved in each exhibition of the pattern" (Whitehead in An enquiry concerning the principles of natural knowledge. Cambridge University Press, London, 1919, p. 198) we show that 3D reconstruction and simulation techniques can contribute to distinguish between (static) patterns and (dynamic) rhythms of apoptosis. By deciphering novel patterns and rhythms of developmental apoptosis, our reconstructions help to reconcile seemingly inconsistent earlier findings in chick and mouse embryos, and to create rules for computer simulations.

Blood-brain barrier development: Systems modeling and predictive toxicology

The blood-brain barrier (BBB) serves as a gateway for passage of drugs, chemicals, nutrients, metabolites, and hormones between vascular and neural compartments in the brain. Here, we review BBB development with regard to the microphysiology of the neurovascular unit (NVU) and the impact of BBB disruption on brain development. Our focus is on modeling these complex systems. Extant in silico models are available as tools to predict the probability of drug/chemical passage across the BBB; in vitro platforms for high-throughput screening and high-content imaging provide novel data streams for profiling chemical-biological interactions; and engineered human cell-based microphysiological systems provide empirical models with which to investigate the dynamics of NVU function. Computational models are needed that bring together kinetic and dynamic aspects of NVU function across gestation and under various physiological and toxicological scenarios. This integration will inform adverse outcome pathways to reduce uncertainty in translating in vitro data and in silico models for use in risk assessments that aim to protect neurodevelopmental health.

Design and validation of an ontology-driven animal-free testing strategy for developmental neurotoxicity testing

Developmental neurotoxicity entails one of the most complex areas in toxicology. Animal studies provide only limited information as to human relevance. A multitude of alternative models have been developed over the years, providing insights into mechanisms of action. We give an overview of fundamental processes in neural tube formation, brain development and neural specification, aiming at illustrating complexity rather than comprehensiveness. We also give a flavor of the wealth of alternative methods in this area. Given the impressive progress in mechanistic knowledge of human biology and toxicology, the time is right for a conceptual approach for designing testing strategies that cover the integral mechanistic landscape of developmental neurotoxicity. The ontology approach provides a framework for defining this landscape, upon which an integral in silico model for predicting toxicity can be built. It subsequently directs the selection of in vitro assays for rate-limiting events in the biological network, to feed parameter tuning in the model, leading to prediction of the toxicological outcome. Validation of such models requires primary attention to coverage of the biological domain, rather than classical predictive value of individual tests. Proofs of concept for such an approach are already available. The challenge is in mining modern biology, toxicology and chemical information to feed intelligent designs, which will define testing strategies for neurodevelopmental toxicity testing.

Recommendation on test readiness criteria for new approach methods in toxicology: Exemplified for developmental neurotoxicity

Multiple non-animal-based test methods have never been formally validated. In order to use such new approach methods (NAMs) in a regulatory context, criteria to define their readiness are necessary. The field of developmental neurotoxicity (DNT) testing is used to exemplify the application of readiness criteria. The costs and number of untested chemicals are overwhelming for in vivo DNT testing. Thus, there is a need for inexpensive, high-throughput NAMs, to obtain initial information on potential hazards, and to allow prioritization for further testing. A background on the regulatory and scientific status of DNT testing is provided showing different types of test readiness levels, depending on the intended use of data from NAMs. Readiness criteria, compiled during a stakeholder workshop, uniting scientists from academia, industry and regulatory authorities are presented. An important step beyond the listing of criteria, was the suggestion for a preliminary scoring scheme. On this basis a (semi)-quantitative analysis process was assembled on test readiness of 17 NAMs with respect to various uses (e.g. prioritization/screening, risk assessment). The scoring results suggest that several assays are currently at high readiness levels. Therefore, suggestions are made on how DNT NAMs may be assembled into an integrated approach to testing and assessment (IATA). In parallel, the testing state in these assays was compiled for more than 1000 compounds. Finally, a vision is presented on how further NAM development may be guided by knowledge of signaling pathways necessary for brain development, DNT pathophysiology, and relevant adverse outcome pathways (AOP).

Rethinking developmental toxicity testing: Evolution or revolution?

BACKGROUND

Current developmental toxicity testing adheres largely to protocols suggested in 1966 involving the administration of test compound to pregnant laboratory animals. After more than 50 years of embryo-fetal development testing, are we ready to consider a different approach to human developmental toxicity testing?

METHODS

A workshop was held under the auspices of the Developmental and Reproductive Toxicology Technical Committee of the ILSI Health and Environmental Sciences Institute to consider how we might design developmental toxicity testing if we started over with 21st century knowledge and techniques (revolution). We first consider what changes to the current protocols might be recommended to make them more predictive for human risk (evolution).

RESULTS

The evolutionary approach includes modifications of existing protocols and can include humanized models, disease models, more accurate assessment and testing of metabolites, and informed approaches to dose selection. The revolution could start with hypothesis-driven testing where we take what we know about a compound or close analog and answer specific questions using targeted experimental techniques rather than a one-protocol-fits-all approach. Central to the idea of hypothesis-driven testing is the concept that testing can be done at the level of mode of action. It might be feasible to identify a small number of key events at a molecular or cellular level that predict an adverse outcome and for which testing could be performed in vitro or in silico or, rarely, using limited in vivo models. Techniques for evaluating these key events exist today or are in development.

DISCUSSION

Opportunities exist for refining and then replacing current developmental toxicity testing protocols using techniques that have already been developed or are within reach.

3S - Systematic, systemic, and systems biology and toxicology

A biological system is more than the sum of its parts - it accomplishes many functions via synergy. Deconstructing the system down to the molecular mechanism level necessitates the complement of reconstructing functions on all levels, i.e., in our conceptualization of biology and its perturbations, our experimental models and computer modelling. Toxicology contains the somewhat arbitrary subclass "systemic toxicities"; however, there is no relevant toxic insult or general disease that is not systemic. At least inflammation and repair are involved that require coordinated signaling mechanisms across the organism. However, the more body components involved, the greater the challenge to reca-pitulate such toxicities using non-animal models. Here, the shortcomings of current systemic testing and the development of alternative approaches are summarized. We argue that we need a systematic approach to integrating existing knowledge as exemplified by systematic reviews and other evidence-based approaches. Such knowledge can guide us in modelling these systems using bioengineering and virtual computer models, i.e., via systems biology or systems toxicology approaches. Experimental multi-organ-on-chip and microphysiological systems (MPS) provide a more physiological view of the organism, facilitating more comprehensive coverage of systemic toxicities, i.e., the perturbation on organism level, without using substitute organisms (animals). The next challenge is to establish disease models, i.e., micropathophysiological systems (MPPS), to expand their utility to encompass biomedicine. Combining computational and experimental systems approaches and the chal-lenges of validating them are discussed. The suggested 3S approach promises to leverage 21st century technology and systematic thinking to achieve a paradigm change in studying systemic effects.

Engineering human cell spheroids to model embryonic tissue fusion in vitro

Epithelial-mesenchymal interactions drive embryonic fusion events during development, and perturbations of these interactions can result in birth defects. Cleft palate and neural tube defects can result from genetic defects or environmental exposures during development, yet very little is known about the effect of chemical exposures on fusion events during human development because of a lack of relevant and robust human in vitro assays of developmental fusion behavior. Given the etiology and prevalence of cleft palate and the relatively simple architecture and composition of the embryonic palate, we sought to develop a three-dimensional culture system that mimics the embryonic palate and could be used to study fusion behavior in vitro using human cells. We engineered size-controlled human Wharton’s Jelly stromal cell (HWJSC) spheroids and established that 7 days of culture in osteogenesis differentiation medium was sufficient to promote an osteogenic phenotype consistent with embryonic palatal mesenchyme. HWJSC spheroids supported the attachment of human epidermal keratinocyte progenitor cells (HPEKp) on the outer spheroid surface likely through deposition of collagens I and IV, fibronectin, and laminin by mesenchymal spheroids. HWJSC spheroids coated in HPEKp cells exhibited fusion behavior in culture, as indicated by the removal of epithelial cells from the seams between spheroids, that was dependent on epidermal growth factor signaling and fibroblast growth factor signaling in agreement with palate fusion literature. The method described here may broadly apply to the generation of three-dimensional epithelial-mesenchymal co-cultures to study developmental fusion events in a format that is amenable to predictive toxicology applications.

Programming microphysiological systems for children's health protection

Microphysiological systems (MPS) and computer simulation models that recapitulate the underlying biology and toxicology of critical developmental transitions are emerging tools for developmental effects assessment of drugs/chemicals. Opportunities and challenges exist for their application to alternative, more public health relevant and efficient chemical toxicity testing methods. This is especially pertinent to children's health research and the evaluation of complex embryological and reproductive impacts of drug/chemical exposure. Scaling these technologies to higher throughput is a key challenge and drives the need for in silico models for quantitative prediction of developmental toxicity to inform safety assessments. One example is cellular agent-based models, constructed from extant embryology, that produce data useful to simulate critical developmental transitions and thereby predict phenotypic consequences of disruption in silico. Biologically inspired MPS models built from human induced pluripotent stem (iPS)-derived cells and synthetic matrices that recapitulate organ-specific physiologies and native tissue architectures are providing exciting new research opportunities to advance the assessment of developmental toxicity and offer the possibility of deriving a full 'human on a chip' system, or a 'Homunculus.' Impact statement This 'commentary' summarizes research needs and opportunities for engineered MPS models for developmental and reproductive toxicity testing. Emerging concepts can be taken forward to a virtual tissue modeling framework for assessing chemical (and non-chemical) stressors on human development. These models will advance children's health research, both basic and translational and new ways to evaluate complex embryological and reproductive impacts of drug and chemical exposures to inform safety assessments.

Retinoic acid in developmental toxicology: Teratogen, morphogen and biomarker

This review explores the usefulness retinoic acid (RA) related physiological factors as possible biomarkers of embryotoxicity. RA is involved in the morphogenesis of the early embryo as well as in the development and maturation of a wide variety of organ anlagen. The region-specific homeostasis of RA in the embryo is in many ways the driving force determining developmental cell proliferation versus differentiation. As a consequence, RA concentrations are carefully controlled in time and space in the developing embryo. RA deficiency and overdosing both result in characteristic patterns of malformations that may involve many different organ systems. The central role of RA in embryo development provides us with a set of sensitive biomarkers that may be employed in developmental toxicity testing. This includes the synthesizing and metabolizing enzymes of RA, but also a myriad of related morphogenetic factors and their genes, of which the expression may be affected by changes in RA balance. Several examples of embryotoxicants interfering with the homeostasis of RA and related parameters have been described. A preliminary adverse outcome pathway framework for RA mediated malformations has been published. Expansion of this framework and its application in developmental toxicity testing may allow the detection of a large variety of embryotoxicants with diverse modes of action. RA homeostasis therefore provides a promising set of molecular tools that may be employed in the advancement of mode of action driven animal-free developmental toxicity testing.

Systems Toxicology: Real World Applications and Opportunities

Systems Toxicology aims to change the basis of how adverse biological effects of xenobiotics are characterized from empirical end points to describing modes of action as adverse outcome pathways and perturbed networks. Toward this aim, Systems Toxicology entails the integration of in vitro and in vivo toxicity data with computational modeling. This evolving approach depends critically on data reliability and relevance, which in turn depends on the quality of experimental models and bioanalysis techniques used to generate toxicological data. Systems Toxicology involves the use of large-scale data streams ("big data"), such as those derived from omics measurements that require computational means for obtaining informative results. Thus, integrative analysis of multiple molecular measurements, particularly acquired by omics strategies, is a key approach in Systems Toxicology. In recent years, there have been significant advances centered on in vitro test systems and bioanalytical strategies, yet a frontier challenge concerns linking observed network perturbations to phenotypes, which will require understanding pathways and networks that give rise to adverse responses. This summary perspective from a 2016 Systems Toxicology meeting, an international conference held in the Alps of Switzerland, describes the limitations and opportunities of selected emerging applications in this rapidly advancing field. Systems Toxicology aims to change the basis of how adverse biological effects of xenobiotics are characterized, from empirical end points to pathways of toxicity. This requires the integration of in vitro and in vivo data with computational modeling. Test systems and bioanalytical technologies have made significant advances, but ensuring data reliability and relevance is an ongoing concern. The major challenge facing the new pathway approach is determining how to link observed network perturbations to phenotypic toxicity.