SkinSensDB: a curated database for skin sensitization assays

In Silico Prediction of Skin Sensitization for Compounds via Flexible Evidence Combination Based on Machine Learning and Dempster-Shafer Theory

Skin sensitization is increasingly becoming a significant concern in the development of drugs and cosmetics due to consumer safety and occupational health problems. In silico methods have emerged as alternatives to traditional in vivo animal testing due to ethical and economic considerations. In this study, machine learning methods were used to build quantitative structure-activity relationship (QSAR) models on five skin sensitization data sets (GPMT, LLNA, DPRA, KeratinoSens, and h-CLAT), achieving effective predictive accuracies (correct classification rates of 0.688-0.764 on test sets). To address the complex mechanisms of human skin sensitization, the Dempster-Shafer theory was applied to merge multiple QSAR models, resulting in an evidence-based integrated decision model. Various evidence combinations and combination rules were explored, with the self-defined Q3 rule showing superior balance. The combination of evidence such as GPMT and KeratinoSens and h-CLAT achieved a correct classification rate (CCR) of 0.880 and coverage of 0.893 while maintaining the competitiveness of other combinations. Additionally, the Shapley additive explanations (SHAP) method was used to interpret important features and substructures related to skin sensitization. A comparative analysis of an external human test set demonstrated the superior performance of the proposed method. Finally, to enhance accessibility, the workflow was implemented into a user-friendly software named HSkinSensDS.

Amino Chemoassay Profiling of Aromatic Aldehydes-Unraveling Drivers of Their Skin Sensitization Potency

Aromatic aldehydes are ubiquitous in humans' everyday life. As aldehydes, they can form imines (Schiff bases) with amino groups of skin proteins, leading to immune response-triggered allergic contact dermatitis. Many known aromatic aldehydes are considered as weak or nonsensitizers, but others like atranol and chloratranol, two components of the fragrance oak moss absolute, show strong sensitization potency. This large discrepancy in potency and, in particular, the underlying reaction mechanisms are only little understood so far. To reduce this knowledge gap, our chemoassay employing glycine-para-nitroanilide (Gly-pNA) as an amino model nucleophile was applied to 23 aromatic aldehydes. The determined Gly-pNA second-order rate constants for imine formation (k₁ ≤ 2.85 L·mol^-1·min^-1) and the imine stability constant (K ≤ 333 L·mol^-1) are on the lower end of the known amino reactivity scale for aldehydes, confirming many aromatic aldehydes as less potent sensitizers in line with animal and human data. The substantially higher sensitization potency of atranol and chloratranol, in turn, is reflected by their unique reaction chemistry profiles, inter alia, identifying them as cross-linkers able to form thermodynamically more stable epitopes with skin proteins (despite low formation kinetics, k₁). The discussion further includes a comparison of experimentally determined k₁ values with computed reactivity data (Taft σ*), the impact of the substitution pattern of the aryl ring on the reactivity with Gly-pNA, and analytically determined adduct patterns. Overall, this work provides new insights into the reaction of aromatic aldehydes with amino groups under aqueous conditions and fosters a better understanding of the chemistry underlying skin sensitization.

Reconfiguring the online tool of SkinSensPred for predicting skin sensitization of pesticides

Adverse outcome pathway (AOP)-based computational models provide state-of-the-art prediction for human skin sensitizers and are promising alternatives to animal testing. However, little is known about their applicability to pesticides due to scarce pesticide data for evaluation. Moreover, pesticides traditionally have been tested on animals without human data, making validation difficult. Direct application of AOP-based models to pesticides may be inappropriate since their original applicability domains were designed to maximize reliability for human response prediction on diverse chemicals but not pesticides. This study proposed to identify a consensus chemical space with concordant human responses predicted by the SkinSensPred online tool and animal testing data to reduce animal testing. The identified consensus chemical space for non-sensitizers achieved high concordance of 85% and 100% for the cross-validation and independent test, respectively. The reconfigured SkinSensPred can be applied as the first-tier tool for identifying non-sensitizers to reduce. animal testing for pesticides by 19.6%.

SkinSensPred as a Promising in Silico Tool for Integrated Testing Strategy on Skin Sensitization

Skin sensitization is an important regulatory endpoint associated with allergic contact dermatitis. Recently, several adverse outcome pathway (AOP)-based alternative methods were developed to replace animal testing for evaluating skin sensitizers. The AOP-based assays were further integrated as a two-out-of-three method with good predictivity. However, the acquisition of experimental data is resource-intensive. In contrast, an integrated testing strategy (ITS) capable of maximizing the usage of laboratory data from AOP-based and in silico methods was developed as defined approaches (DAs) to both hazard and potency assessment. There are currently two in silico models, namely Derek Nexus and OECD QSAR Toolbox, evaluated in the OECD Testing Guideline No. 497. Since more advanced machine learning algorithms have been proposed for skin sensitization prediction, it is therefore desirable to evaluate their performance under the ITS framework. This study evaluated the performance of a new ITS DA (ITS-SkinSensPred) adopting a transfer learning-based SkinSensPred model. Results showed that the ITS-SkinSensPred has similar or slightly better performance compared to the other ITS models. SkinSensPred-based ITS is expected to be a promising method for assessing skin sensitization.

A machine learning-driven approach for prioritizing food contact chemicals of carcinogenic concern based on complementary in silico methods

Carcinogenicity is one of the most critical endpoints for the risk assessment of food contact chemicals (FCCs). However, the carcinogenicity of FCCs remains insufficiently investigated. To fill the data gap, the application of standard experimental methods for identifying chemicals of carcinogenic concerns from a large set of FCCs is impractical due to their resource-intensive nature. In contrast, computational methods provide an efficient way to quickly screen chemicals with carcinogenic potential for subsequent experimental validation. Since every model was developed based on a limited number of training samples, the use of single models for carcinogenicity assessment may not cover the complex mechanisms of carcinogenesis. This study proposed a novel machine learning-based weight-of-evidence (WoE) model for prioritizing chemical carcinogenesis. The WoE model can nonlinearly integrate complementary computational methods of structural alerts, quantitative structure-activity relationship models and in silico toxicogenomics models into a WoE-score. Compared to the best single method, the WoE model gained 8% and 19.7% improvement in the area under the receiver operating characteristic curve (AUC) value and chemical coverage, respectively. The prioritization of 1623 FCCs concludes 44 chemicals of high carcinogenic concern. The machine learning-based WoE approach provides a fast and comprehensive way for prioritizing chemicals of carcinogenic concern.

Amino Reactivity of Glutardialdehyde and Monoaldehydes─Chemoassay Profile vs Skin Sensitization Potency

Chemoassay profiling of organic electrophiles through the direct peptide reactivity assay has become an OECD-accepted nonanimal component in the REACH evaluation of potential skin sensitizers. For aldehydes forming imines (Schiff bases), however, existing chemoassays yielded inconclusive results, indicating issues with their NH₂ sensitivity and the reversibility of the reaction. In the present study, a new kinetic chemoassay employing the N terminus of glycine-para-nitroanilide, Gly-pNA, as a model nucleophile for protein NH₂ groups is introduced and applied to nine aliphatic monoaldehydes and glutardialdehyde (1,5-pentanedial) that have log K_ow (octanol/water partition coefficient) values from 0.63 to 3.99. The Gly-pNA second-order rate constants k₁ range from 8.56 to 150 L·mol^-1·min^-1 for the monoaldehydes. Interestingly, glutardialdehyde with a k₁ of 17 731 L·mol^-1·min^-1 is 170-fold more reactive than its monoaldehyde counterpart pentanal. This can be rationalized by hydration or tautomerization of the dialdehyde to monoaldehydic forms, now facilitating Schiff base formation through an intramolecular H bond. Comparison with murine local lymph node assay data from the literature reveals that adduct stability in terms of reaction thermodynamics (K = k₁/k_-1^pseudo) rather than formation kinetics (k₁) governs the skin sensitization potency of Schiff-base-forming aldehydes. The discussion includes analytically determined adduct patterns, and the impact of α- and β-carbon substitution as well as hydrophobicity on aldehyde reactivity.

In silico Prediction of Skin Sensitization: Quo vadis?

Skin direct contact with chemical or physical substances is predisposed to allergic contact dermatitis (ACD), producing various allergic reactions, namely rash, blister, or itchy, in the contacted skin area. ACD can be triggered by various extremely complicated adverse outcome pathways (AOPs) remains to be causal for biosafety warrant. As such, commercial products such as ointments or cosmetics can fulfill the topically safe requirements in animal and non-animal models including allergy. Europe, nevertheless, has banned animal tests for the safety evaluations of cosmetic ingredients since 2013, followed by other countries. A variety of non-animal in vitro tests addressing different key events of the AOP, the direct peptide reactivity assay (DPRA), KeratinoSens™, LuSens and human cell line activation test h-CLAT and U-SENS™ have been developed and were adopted in OECD test guideline to identify the skin sensitizers. Other methods, such as the SENS-IS are not yet fully validated and regulatorily accepted. A broad spectrum of in silico models, alternatively, to predict skin sensitization have emerged based on various animal and non-animal data using assorted modeling schemes. In this article, we extensively summarize a number of skin sensitization predictive models that can be used in the biopharmaceutics and cosmeceuticals industries as well as their future perspectives, and the underlined challenges are also discussed.

Advancing computer-aided drug discovery (CADD) by big data and data-driven machine learning modeling

Advancing a new drug to market requires substantial investments in time as well as financial resources. Crucial bioactivities for drug candidates, including their efficacy, pharmacokinetics (PK), and adverse effects, need to be investigated during drug development. With advancements in chemical synthesis and biological screening technologies over the past decade, a large amount of biological data points for millions of small molecules have been generated and are stored in various databases. These accumulated data, combined with new machine learning (ML) approaches, such as deep learning, have shown great potential to provide insights into relevant chemical structures to predict in vitro, in vivo, and clinical outcomes, thereby advancing drug discovery and development in the big data era.

Skin sensitization in silico protocol

The assessment of skin sensitization has evolved over the past few years to include in vitro assessments of key events along the adverse outcome pathway and opportunistically capitalize on the strengths of in silico methods to support a weight of evidence assessment without conducting a test in animals. While in silico methods vary greatly in their purpose and format; there is a need to standardize the underlying principles on which such models are developed and to make transparent the implications for the uncertainty in the overall assessment. In this contribution, the relationship between skin sensitization relevant effects, mechanisms, and endpoints are built into a hazard assessment framework. Based on the relevance of the mechanisms and effects as well as the strengths and limitations of the experimental systems used to identify them, rules and principles are defined for deriving skin sensitization in silico assessments. Further, the assignments of reliability and confidence scores that reflect the overall strength of the assessment are discussed. This skin sensitization protocol supports the implementation and acceptance of in silico approaches for the prediction of skin sensitization.

Transfer learning for predicting human skin sensitizers

Computational prioritization of chemicals for potential skin sensitization risks plays essential roles in the risk assessment of environmental chemicals and drug development. Given the huge number of chemicals for testing, computational methods enable the fast identification of high-risk chemicals for experimental validation and design of safer alternatives. However, the development of robust prediction model requires a large dataset of tested chemicals that is usually not available for most toxicological endpoints, especially for human data. A small training dataset makes the development of effective models difficult with insufficient coverage and accuracy. In this study, an ensemble tree-based multitask learning method was developed incorporating three relevant tasks in the well-defined adverse outcome pathway (AOP) of skin sensitization to transfer shared knowledge to the major task of human sensitizers. The results show both largely improved coverage and accuracy compared with three state-of-the-art methods. A user-friendly prediction server was available at https://cwtung.kmu.edu.tw/skinsensdb/predict . As AOPs for various toxicity endpoints are being actively developed, the proposed method can be applied to develop prediction models for other endpoints.

Adverse outcome pathways: a concise introduction for toxicologists

Adverse outcome pathways (AOPs) are designed to provide a clear-cut mechanistic representation of critical toxicological effects that propagate over different layers of biological organization from the initial interaction of a chemical with a molecular target to an adverse outcome at the individual or population level. Adverse outcome pathways are currently gaining momentum, especially in view of their many potential applications as pragmatic tools in the fields of human toxicology, ecotoxicology, and risk assessment. A number of guidance documents, issued by the Organization for Economic Cooperation and Development, as well as landmark papers, outlining best practices to develop, assess and use AOPs, have been published in the last few years. The present paper provides a synopsis of the main principles related to the AOP framework for the toxicologist less familiar with this area, followed by two case studies relevant for human toxicology and ecotoxicology.