We demonstrate a computational network model that integrates 18 in vitro, high-throughput screening assays measuring estrogen receptor (ER) binding, dimerization, chromatin binding, transcriptional activation, and ER-dependent cell proliferation. The network model uses activity patterns across the in vitro assays to predict whether a chemical is an ER agonist or antagonist, or is otherwise influencing the assays through a manner dependent on the physics and chemistry of the technology platform ("assay interference"). The method is applied to a library of 1812 commercial and environmental chemicals, including 45 ER positive and negative reference chemicals. Among the reference chemicals, the network model correctly identified the agonists and antagonists with the exception of very weak compounds whose activity was outside the concentration range tested. The model agonist score also correlated with the expected potency class of the active reference chemicals. Of the 1812 chemicals evaluated, 111 (6.1%) were predicted to be strongly ER active in agonist or antagonist mode. This dataset and model were also used to begin a systematic investigation of assay interference. The most prominent cause of false-positive activity (activity in an assay that is likely not due to interaction of the chemical with ER) is cytotoxicity. The model provides the ability to prioritize a large set of important environmental chemicals with human exposure potential for additional in vivo endocrine testing. Finally, this model is generalizable to any molecular pathway for which there are multiple upstream and downstream assays available.
Background:Humans are exposed to thousands of man-made chemicals in the environment. Some chemicals mimic natural endocrine hormones and, thus, have the potential to be endocrine disruptors. Most of these chemicals have never been tested for their ability to interact with the estrogen receptor (ER). Risk assessors need tools to prioritize chemicals for evaluation in costly in vivo tests, for instance, within the U.S. EPA Endocrine Disruptor Screening Program.Objectives:We describe a large-scale modeling project called CERAPP (Collaborative Estrogen Receptor Activity Prediction Project) and demonstrate the efficacy of using predictive computational models trained on high-throughput screening data to evaluate thousands of chemicals for ER-related activity and prioritize them for further testing.Methods:CERAPP combined multiple models developed in collaboration with 17 groups in the United States and Europe to predict ER activity of a common set of 32,464 chemical structures. Quantitative structure–activity relationship models and docking approaches were employed, mostly using a common training set of 1,677 chemical structures provided by the U.S. EPA, to build a total of 40 categorical and 8 continuous models for binding, agonist, and antagonist ER activity. All predictions were evaluated on a set of 7,522 chemicals curated from the literature. To overcome the limitations of single models, a consensus was built by weighting models on scores based on their evaluated accuracies.Results:Individual model scores ranged from 0.69 to 0.85, showing high prediction reliabilities. Out of the 32,464 chemicals, the consensus model predicted 4,001 chemicals (12.3%) as high priority actives and 6,742 potential actives (20.8%) to be considered for further testing.Conclusion:This project demonstrated the possibility to screen large libraries of chemicals using a consensus of different in silico approaches. This concept will be applied in future projects related to other end points.Citation:Mansouri K, Abdelaziz A, Rybacka A, Roncaglioni A, Tropsha A, Varnek A, Zakharov A, Worth A, Richard AM, Grulke CM, Trisciuzzi D, Fourches D, Horvath D, Benfenati E, Muratov E, Wedebye EB, Grisoni F, Mangiatordi GF, Incisivo GM, Hong H, Ng HW, Tetko IV, Balabin I, Kancherla J, Shen J, Burton J, Nicklaus M, Cassotti M, Nikolov NG, Nicolotti O, Andersson PL, Zang Q, Politi R, Beger RD, Todeschini R, Huang R, Farag S, Rosenberg SA, Slavov S, Hu X, Judson RS. 2016. CERAPP: Collaborative Estrogen Receptor Activity Prediction Project. Environ Health Perspect 124:1023–1033; http://dx.doi.org/10.1289/ehp.1510267
The reaction exoergicity of K + HgBr\ensuremath{\rightarrow}KBr + HgBr ($\ensuremath{\Delta}{D}_{0}=\ensuremath{-}0.86\ifmmode\pm\else\textpm\fi{}0.09$ eV) plus the energy of a photon of $\ensuremath{\lambda}<606\ifmmode\pm\else\textpm\fi{}26$ nm provides sufficient energy to open the channel to HgBr* ($\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{B}^{2}\ensuremath{\Sigma}^{+}$). Although neither reagents nor products absorb near 595 nm, luminescence at 500 nm is observed from the crossing region of beams of K, Hg${\mathrm{Br}}_{2}$, and photons ($\ensuremath{\lambda}\ensuremath{\sim}595$ nm). This emission is attributed to HgBr* formed by light absorption during reactive collision. The cross section is \ensuremath{\sim} ${10}^{\ensuremath{-}17}$ ${\mathrm{cm}}^{2}$ at \ensuremath{\sim} 2 MW/${\mathrm{cm}}^{2}$.
High-content imaging (HCI) allows simultaneous measurement of multiple cellular phenotypic changes and is an important tool for evaluating the biological activity of chemicals.Our goal was to analyze dynamic cellular changes using HCI to identify the "tipping point" at which the cells did not show recovery towards a normal phenotypic state.HCI was used to evaluate the effects of 967 chemicals (in concentrations ranging from 0.4 to 200 μM) on HepG2 cells over a 72-hr exposure period. The HCI end points included p53, c-Jun, histone H2A.x, α-tubulin, histone H3, alpha tubulin, mitochondrial membrane potential, mitochondrial mass, cell cycle arrest, nuclear size, and cell number. A computational model was developed to interpret HCI responses as cell-state trajectories.Analysis of cell-state trajectories showed that 336 chemicals produced tipping points and that HepG2 cells were resilient to the effects of 334 chemicals up to the highest concentration (200 μM) and duration (72 hr) tested. Tipping points were identified as concentration-dependent transitions in system recovery, and the corresponding critical concentrations were generally between 5 and 15 times (25th and 75th percentiles, respectively) lower than the concentration that produced any significant effect on HepG2 cells. The remaining 297 chemicals require more data before they can be placed in either of these categories.These findings show the utility of HCI data for reconstructing cell state trajectories and provide insight into the adaptation and resilience of in vitro cellular systems based on tipping points. Cellular tipping points could be used to define a point of departure for risk-based prioritization of environmental chemicals.Shah I, Setzer RW, Jack J, Houck KA, Judson RS, Knudsen TB, Liu J, Martin MT, Reif DM, Richard AM, Thomas RS, Crofton KM, Dix DJ, Kavlock RJ. 2016. Using ToxCast™ data to reconstruct dynamic cell state trajectories and estimate toxicological points of departure. Environ Health Perspect 124:910-919; http://dx.doi.org/10.1289/ehp.1409029.
Many parameters treated as constants in traditional physiologically based pharmacokinetic models must be formulated as time-varying quantities when modeling pregnancy and gestation due to the dramatic physiological and anatomical changes that occur during this period. While several collections of empirical models for such parameters have been published, each has shortcomings. We sought to create a repository of empirical models for tissue volumes, blood flow rates, and other quantities that undergo substantial changes in a human mother and her fetus during the time between conception and birth, and to address deficiencies with similar, previously published repositories. We used maximum likelihood estimation to calibrate various models for the time-varying quantities of interest, and then used the Akaike information criterion to select an optimal model for each quantity. For quantities of interest for which time-course data were not available, we constructed composite models using percentages and/or models describing related quantities. In this way, we developed a comprehensive collection of formulae describing parameters essential for constructing a PBPK model of a human mother and her fetus throughout the approximately 40 weeks of pregnancy and gestation. We included models describing blood flow rates through various fetal blood routes that have no counterparts in adults. Our repository of mathematical models for anatomical and physiological quantities of interest provides a basis for PBPK models of human pregnancy and gestation, and as such, it can ultimately be used to support decision-making with respect to optimal pharmacological dosing and risk assessment for pregnant women and their developing fetuses. The views expressed in this article are those of the authors and do not necessarily represent the views or policies of the U.S. Environmental Protection Agency.
CCEMD (Center for Computational Engineering Molecular Dynamics) is a general purpose molecular dynamics program written in the C language, built on top of the MD program of Windemuth and Schulten. CCEMD can perform molecular dynamics, gradient minimization, genetic algorithm-based conformation searching, and ligand-protein docking. This report documents the algorithms used and provides a users` guide for the code.
In order to compare between in vivo toxicity studies, dosimetry is needed to translate study-specific dose regimens into dose metrics such as tissue concentration. These tissue concentrations may then be compared with in vitro bioactivity assays to perhaps identify mechanisms relevant to the lowest observed effect level (LOEL) dose group and the onset of the observed in vivo toxicity. Here, we examine the perfluorinated compounds (PFCs) perfluorooctanoate (PFOA) and perfluorooctanesulfonate (PFOS). We analyzed 9 in vivo toxicity studies for PFOA and 13 in vivo toxicity studies for PFOS. Both PFCs caused multiple effects in various test species, strains, and genders. We used a Bayesian pharmacokinetic (PK) modeling framework to incorporate data from 6 PFOA PK studies and 2 PFOS PK studies (conducted in 3 species) to predict dose metrics for the in vivo LOELs and no observed effect levels (NOELs). We estimated PK parameters for 11 combinations of chemical, species, strain, and gender. Despite divergent study designs and species-specific PK, for a given effect, we found that the predicted dose metrics corresponding to the LOELs (and NOELs where available) occur at similar concentrations. In vitro assay results for PFOA and PFOS from EPA's ToxCast project were then examined. We found that most in vitro bioactivity occurs at concentrations lower than the predicted concentrations for the in vivo LOELs and higher than the predicted concentrations for the in vivo NOELs (where available), for a variety of nonimmunological effects. These results indicate that given sufficient PK data, the in vivo LOELs dose regimens, but not necessarily the effects, could have been predicted from in vitro studies for these 2 PFCs.
Presentation for an invited talk at the Endocrine Platform Session at the Society of Environmental Toxicology and Chemistry 40th annual meeting held November 2019 in Toronto Canada.
