Hepatic toxicity is a key concern for novel pharmaceutical drugs since it is difficult to anticipate in preclinical models, and it can originate from pharmacologically unrelated drug effects, such as pathway interference, metabolism, and drug accumulation. Because liver toxicity still ranks among the top reasons for drug attrition, the reliable prediction of adverse hepatic effects is a substantial challenge in drug discovery and development. To this end, more effort needs to be focused on the development of improved predictive in-vitro and in-silico approaches. Current computational models often lack applicability to novel pharmaceutical candidates, typically due to insufficient coverage of the chemical space of interest, which is either imposed by size or diversity of the training data. Hence, there is an urgent need for better computational models to allow for the identification of safe drug candidates and to support experimental design. In this context, a large data set comprising 3712 compounds with liver related toxicity findings in humans and animals was collected from various sources. The complex pathology was clustered into 21 preclinical and human hepatotoxicity endpoints, which were organized into three levels of detail. Support vector machine models were trained for each endpoint, using optimized descriptor sets from chemometrics software. The optimized global human hepatotoxicity model has high sensitivity (68%) and excellent specificity (95%) in an internal validation set of 221 compounds. Models for preclinical endpoints performed similarly. To allow for reliable prediction of "truly external" novel compounds, all predictions are tagged with confidence parameters. These parameters are derived from a statistical analysis of the predictive probability densities. The whole approach was validated for an external validation set of 269 proprietary compounds. The models are fully integrated into our early safety in-silico workflow.
The electrophilic reactivity of Michael acceptors is an important determinant of their toxicity. For a set of 35 α,β-unsaturated aldehydes, ketones and esters with experimental rate constants of their reaction with glutathione (GSH), k(GSH), quantum chemical transition-state calculations of the corresponding Michael addition of the model nucleophile methane thiol (CH(3)SH) have been performed at the B3LYP/6-31G** level, focusing on the 1,2-olefin addition pathway without and with initial protonation. Inclusion of Boltzmann-weighting of conformational flexibility yields intrinsic reaction barriers ΔE(‡) that for the case of initial protonation correctly reflect the structural variation of k(GSH) across all three compound classes, except that they fail to account for a systematic (essentially incremental) decrease in reactivity upon α-substitution. By contrast, the reduction in k(GSH) through β-substitution is well captured by ΔE(‡). Empirical correction for the α-substitution effect yields a high squared correlation coefficient (r(2) = 0.96) for the quantum chemical prediction of log k(GSH), thus enabling an in silico screening of the toxicity-relevant electrophilicity of α,β-unsaturated carbonyls. The latter is demonstrated through application of the calculation scheme for a larger set of 46 Michael-acceptor aldehydes, ketones and esters with experimental values for their toxicity toward the ciliates Tetrahymena pyriformis in terms of 50% growth inhibition values after 48 h exposure (EC(50)). The developed approach may add in the predictive hazard evaluation of α,β-unsaturated carbonyls such as for the European REACH (Registration, Evaluation, Authorization and Restriction of Chemicals) Directive, enabling in particular an early identification of toxicity-relevant Michael-acceptor reactivity.
For a Michael-acceptor set of 45 α,β-unsaturated esters, the 2(nd) -order rate constant of reaction with glutathione, log kGSH , was modeled through the quantum chemical reaction barrier (ΔE(≠) ) employing methane thiol as model nucleophile. Regression of their 48-h toxicity toward the ciliates Tetrahymena pyriformis (log EC50 , 50 % growth inhibition) on log Kow (octanol/water partition coefficient) and log kGSH revealed a variation in the relative weights of hydrophobicity and electrophilic reactivity as determinants of the aquatic toxicity. The difference DKk =log Kow -log kGSH turned out as a suitable means for predictively discriminating between narcosis-level (DKk >3.0) and excess-toxic (DKk <2.0) compounds. In the intermediate DKk range (2.0≤DKk ≤3.0), both narcosis-level and reactive-toxicity models are applicable for predicting aquatic toxicity. As such, DKk represents the chemoavailability of Michael-acceptor esters, characterizing their likelihood for undertaking covalent reactions with thiol sites of endogenous peptides and proteins. At the same time, DKk introduces a straightforward way for characterizing the applicability domain of QSAR (quantitative structure-activity relationship) models for predicting the toxicity of Michael-acceptor esters. The resultant model suite comprising QSARs for reactive toxicity and baseline narcosis is triggered by the compounds' chemoavailability, and yields predictions superior to existing approaches.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.