Structure-based virtual screening of CYP1A1 inhibitors: towards rapid tier-one assessment of potential developmental toxicants.

Cytochrome P450 1A1 (CYP1A1) metabolizes estrogens, melatonin, and other key endogenous signaling molecules critical for embryonic/fetal development. The enzyme has increasing expression during pregnancy, and its inhibition or knockout increases embryonic/fetal lethality and/or developmental problems. Here, we present a virtual screening model for CYP1A1 inhibitors based on the orthosteric and predicted allosteric sites of the enzyme. Using 1001 reference compounds with CYP1A1 activity data, we optimized the decision thresholds of our model and classified the training compounds with 68.3% balanced accuracy (91.0% sensitivity and 45.7% specificity). We applied our final model to 11 known CYP1A1 orthosteric binders and related compounds, and found that our ranking of the known orthosteric binders generally agrees with the relative activity of CYP1A1 in metabolizing these compounds. We also applied the model to 22 new test compounds with unknown/unclear CYP1A1 inhibitory activity, and predicted 16 of them are CYP1A1 inhibitors. The CYP1A1 potency and modes of inhibition of these 22 compounds were experimentally determined. We confirmed that most predicted inhibitors, including drugs contraindicated during pregnancy (amiodarone, bicalutamide, cyproterone acetate, ketoconazole, and tamoxifen) and environmental agents suspected to be endocrine disruptors (bisphenol A, diethyl and dibutyl phthalates, and zearalenone), are indeed potent inhibitors of CYP1A1. Our results suggest that virtual screening may be used as a rapid tier-one method to screen for potential CYP1A1 inhibitors, and flag them out for further experimental evaluations.

Optimum concentration-response curve metrics for supervised selection of discriminative cellular phenotypic endpoints for chemical hazard assessment.

High-content imaging (HCI) provides quantitative and information-rich measurements of chemical effects on human in vitro cell models. Identification of discriminative phenotypic endpoints from cellular features obtained from HCI is required for accurate assessments of potential chemical hazards. However, the use of suboptimal metrics to quantify the concentration-response curves (CRC) of chemicals based on these features may obscure discriminative features, and lead to non-predictive endpoints and poor chemical classifications or hazard assessments. Here, we present a systematic and data-driven study on the performances of different CRC metrics in identifying image-based phenotypic features that can accurately classify the effects of reference chemicals with known in vivo toxicities. We studied four previous HCI in vitro nephro- or pulmono-toxicity datasets, which contain phenotypic feature measurements from different cell and feature types. Within a feature type, we found that efficacy metrics at higher chemical concentrations tend to give higher classification accuracy, whereas potency metrics do not have obvious trends across different response levels. Across different cell and feature types, efficacy metrics generally gave higher classification accuracy than potency metrics and area under the curve (AUC). Our results suggest that efficacy metrics, especially at higher concentrations, are more likely to help us to identify discriminative phenotypic endpoints. Therefore, HCI experiments for toxicological applications should include measurements at sufficiently high chemical concentrations, and efficacy metrics should always be analyzed. The identified features may be used as specific toxicity endpoints for further chemical hazard assessment.

A Case Study with Triazole Fungicides to Explore Practical Application of Next-Generation Hazard Assessment Methods for Human Health.

The ongoing developments in chemical risk assessment have led to new concepts building on integration of sophisticated nonanimal models for hazard characterization. Here we explore a pragmatic approach for implementing such concepts, using a case study of three triazole fungicides, namely, flusilazole, propiconazole, and cyproconazole. The strategy applied starts with evaluating the overall level of concern by comparing exposure estimates to toxicological potential, followed by a combination of in silico tools and literature-derived high-throughput screening assays and computational elaborations to obtain insight into potential toxicological mechanisms and targets in the organism. Additionally, some targeted in vitro tests were evaluated for their utility to confirm suspected mechanisms of toxicity and to generate points of departure. Toxicological mechanisms instead of the current "end point-by-end point" approach should guide the selection of methods and assays that constitute a toolbox for next-generation risk assessment. Comparison of the obtained in silico and in vitro results with data from traditional in vivo testing revealed that, overall, nonanimal methods for hazard identification can produce adequate qualitative hazard information for risk assessment. Follow-up studies are needed to further refine the proposed approach, including the composition of the toolbox, toxicokinetics models, and models for exposure assessment.