Quantitative structure–property relationship modeling of skin sensitization: A quantitative prediction

PreS/MD: Predictor of Sensitization Hazard for Chemical Substances Released From Medical Devices

In the United States, a pre-market regulatory submission for any medical device that comes into contact with either a patient or the clinical practitioner must include an adequate toxicity evaluation of chemical substances that can be released from the device during its intended use. These substances, also referred to as extractables and leachables, must be evaluated for their potential to induce sensitization/allergenicity, which traditionally has been done in animal assays such as the guinea pig maximization test (GPMT). However, advances in basic and applied science are continuously presenting opportunities to employ new approach methodologies, including computational methods which, when qualified, could replace animal testing methods to support regulatory submissions. Herein, we developed a new computational tool for rapid and accurate prediction of the GPMT outcome that we have named PreS/MD (predictor of sensitization for medical devices). To enable model development, we (1) collected, curated, and integrated the largest publicly available dataset for GPMT results; (2) succeeded in developing externally predictive (balanced accuracy of 70%-74% as evaluated by both 5-fold external cross-validation and testing of novel compounds) quantitative structure-activity relationships (QSAR) models for GPMT using machine learning algorithms, including deep learning; and (3) developed a publicly accessible web portal integrating PreS/MD models that can predict GPMT outcomes for any molecule of interest. We expect that PreS/MD will be used by both industry and regulatory scientists in medical device safety assessments and help replace, reduce, or refine the use of animals in toxicity testing. PreS/MD is freely available at https://presmd.mml.unc.edu/.

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.

SkinSensDB: a curated database for skin sensitization assays

Skin sensitization is an important toxicological endpoint for chemical hazard determination and safety assessment. Prediction of chemical skin sensitizer had traditionally relied on data from rodent models. The development of the adverse outcome pathway (AOP) and associated alternative in vitro assays have reshaped the assessment of skin sensitizers. The integration of multiple assays as key events in the AOP has been shown to have improved prediction performance. Current computational models to predict skin sensitization mainly based on in vivo assays without incorporating alternative in vitro assays. However, there are few freely available databases integrating both the in vivo and the in vitro skin sensitization assays for development of AOP-based skin sensitization prediction models. To facilitate the development of AOP-based prediction models, a skin sensitization database named SkinSensDB has been constructed by curating data from published AOP-related assays. In addition to providing datasets for developing computational models, SkinSensDB is equipped with browsing and search tools which enable the assessment of new compounds for their skin sensitization potentials based on data from structurally similar compounds. SkinSensDB is publicly available at http://cwtung.kmu.edu.tw/skinsensdb.

QSAR models of human data can enrich or replace LLNA testing for human skin sensitization

Skin sensitization is a major environmental and occupational health hazard. Although many chemicals have been evaluated in humans, there have been no efforts to model these data to date. We have compiled, curated, analyzed, and compared the available human and LLNA data. Using these data, we have developed reliable computational models and applied them for virtual screening of chemical libraries to identify putative skin sensitizers. The overall concordance between murine LLNA and human skin sensitization responses for a set of 135 unique chemicals was low (R = 28-43%), although several chemical classes had high concordance. We have succeeded to develop predictive QSAR models of all available human data with the external correct classification rate of 71%. A consensus model integrating concordant QSAR predictions and LLNA results afforded a higher CCR of 82% but at the expense of the reduced external dataset coverage (52%). We used the developed QSAR models for virtual screening of CosIng database and identified 1061 putative skin sensitizers; for seventeen of these compounds, we found published evidence of their skin sensitization effects. Models reported herein provide more accurate alternative to LLNA testing for human skin sensitization assessment across diverse chemical data. In addition, they can also be used to guide the structural optimization of toxic compounds to reduce their skin sensitization potential.

Predicting chemically-induced skin reactions. Part I: QSAR models of skin sensitization and their application to identify potentially hazardous compounds

Repetitive exposure to a chemical agent can induce an immune reaction in inherently susceptible individuals that leads to skin sensitization. Although many chemicals have been reported as skin sensitizers, there have been very few rigorously validated QSAR models with defined applicability domains (AD) that were developed using a large group of chemically diverse compounds. In this study, we have aimed to compile, curate, and integrate the largest publicly available dataset related to chemically-induced skin sensitization, use this data to generate rigorously validated and QSAR models for skin sensitization, and employ these models as a virtual screening tool for identifying putative sensitizers among environmental chemicals. We followed best practices for model building and validation implemented with our predictive QSAR workflow using Random Forest modeling technique in combination with SiRMS and Dragon descriptors. The Correct Classification Rate (CCR) for QSAR models discriminating sensitizers from non-sensitizers was 71-88% when evaluated on several external validation sets, within a broad AD, with positive (for sensitizers) and negative (for non-sensitizers) predicted rates of 85% and 79% respectively. When compared to the skin sensitization module included in the OECD QSAR Toolbox as well as to the skin sensitization model in publicly available VEGA software, our models showed a significantly higher prediction accuracy for the same sets of external compounds as evaluated by Positive Predicted Rate, Negative Predicted Rate, and CCR. These models were applied to identify putative chemical hazards in the Scorecard database of possible skin or sense organ toxicants as primary candidates for experimental validation.

Recombinant human hyaluronidase PH20 (rHuPH20) facilitates subcutaneous infusions of large volumes of immunoglobulin in a swine model

Many patients with primary immunodeficiency disease (PIDD) require lifelong immunoglobulin (Ig) replacement therapy. Home-based subcutaneous (SC) infusion provides advantages to patients with PIDD compared to hospital-based intravenous infusion. One limitation of current practice with SCIg infusion is the need for small-volume infusions at multiple injection sites on a frequent basis. A method was developed for large-volume SC infusion that uses preinfusion of recombinant human hyaluronidase (rHuPH20) to facilitate fluid dispersion. Miniature swine was used as a preclinical model to assess the effects of rHuPH20-facilitated infusions, of a single monthly dose, on fluid dispersion, infusion-related pressure, swelling, induration, and tissue damage. Preinfusion of vehicle (control) or rHuPH20 (75 U/g Ig) was performed simultaneously on contralateral abdominal sites on each animal, followed by infusion of 300 mL 10 % Ig (30 g) at each site. Compared to control infusions, rHuPH20 significantly reduced infusion pressure and induration (p < 0.05) and accelerated postinfusion Ig dispersion. Histological evaluation of infusion site tissue showed moderate to severe swelling for the control. Swelling after rHuPH20-facilitated infusion was mild on day 1 and had completely resolved shortly thereafter. Laser Doppler imaging of control infusion sites revealed local cutaneous hypoperfusion during Ig infusion, which was reduced almost 7-fold (p < 0.05) with the use of rHuPH20. These results demonstrate that rHuPH20-facilitated Ig infusion is associated with improved dispersion of Ig, resulting in reduced tissue pressure, induration, and reduced risk of tissue damage from mechanical trauma or local ischemia, thus enabling SC administration of large volumes of Ig at a single site.

Use of artificial neural networks in the QSAR prediction of physicochemical properties and toxicities for REACH legislation

With the introduction of the REACH legislation in the European Union, there is a requirement for property and toxicity data on chemicals produced in or imported into the EU at levels of 1 tonne/year or more. This has meant an increase in the in silico prediction of such data. One of the chief predictive approaches is QSAR (quantitative structure-activity relationships), which is widely used in many fields. A QSAR approach that is increasingly being used is that of artificial neural networks (ANNs), and this chapter discusses its application to the range of physicochemical properties and toxicities required by REACH. ANNs generally outperform the main QSAR approach of multiple linear regression (MLR), although other approaches such as support vector machines sometimes outperform ANNs. Most ANN QSARs reported to date comply with only two of the five OECD Guidelines for the Validation of (Q)SARs.

Virtual design of chemical penetration enhancers for transdermal drug delivery

Traditional drug design is a laborious and expensive process that often challenges the pharmaceutical industries. As a result, researchers have turned to computational methods for computer-assisted molecular design. Recently, genetic and evolutionary algorithms have emerged as efficient methods in solving combinatorial problems associated with computer-aided molecular design. Further, combining genetic algorithms with quantitative structure-property relationship analyses has proved effective in drug design. In this work, we have integrated a new genetic algorithm and nonlinear quantitative structure-property relationship models to develop a reliable virtual screening algorithm for the generation of potential chemical penetration enhancers. The genetic algorithms-quantitative structure-property relationship algorithm has been implemented successfully to identify potential chemical penetration enhancers for transdermal drug delivery of insulin. Validation of the newly identified chemical penetration enhancer molecular structures was conducted through carefully designed experiments, which elucidated the cytotoxicity and permeability of the chemical penetration enhancers.

Perspectives on Non-Animal Alternatives for Assessing Sensitization Potential in Allergic Contact Dermatitis

Skin sensitization remains a major environmental and occupational health hazard. Animal models have been used as the gold standard method of choice for estimating chemical sensitization potential. However, a growing international drive and consensus for minimizing animal usage have prompted the development of in vitro methods to assess chemical sensitivity. In this paper, we examine existing approaches including in silico models, cell and tissue based assays for distinguishing between sensitizers and irritants. The in silico approaches that have been discussed include Quantitative Structure Activity Relationships (QSAR) and QSAR based expert models that correlate chemical molecular structure with biological activity and mechanism based read-across models that incorporate compound electrophilicity. The cell and tissue based assays rely on an assortment of mono and co-culture cell systems in conjunction with 3D skin models. Given the complexity of allergen induced immune responses, and the limited ability of existing systems to capture the entire gamut of cellular and molecular events associated with these responses, we also introduce a microfabricated platform that can capture all the key steps involved in allergic contact sensitivity. Finally, we describe the development of an integrated testing strategy comprised of two or three tier systems for evaluating sensitization potential of chemicals.

An integrated QSAR-PBPK modelling approach for predicting the inhalation toxicokinetics of mixtures of volatile organic chemicals in the rat

The objective of this study was to predict the inhalation toxicokinetics of chemicals in mixtures using an integrated QSAR-PBPK modelling approach. The approach involved: (1) the determination of partition coefficients as well as V(max) and K(m) based solely on chemical structure for 53 volatile organic compounds, according to the group contribution approach; and (2) using the QSAR-driven coefficients as input in interaction-based PBPK models in the rat to predict the pharmacokinetics of chemicals in mixtures of up to 10 components (benzene, toluene, m-xylene, o-xylene, p-xylene, ethylbenzene, dichloromethane, trichloroethylene, tetrachloroethylene, and styrene). QSAR-estimated values of V(max) varied compared with experimental results by a factor of three for 43 out of 53 studied volatile organic compounds (VOCs). K(m) values were within a factor of three compared with experimental values for 43 out of 53 VOCs. Cross-validation performed as a ratio of predicted residual sum of squares and sum of squares of the response value indicates a value of 0.108 for V(max) and 0.208 for K(m). The integration of QSARs for partition coefficients, V(max) and K(m), as well as setting the K(m) equal to K(i) (metabolic inhibition constant) within the mixture PBPK model allowed to generate simulations of the inhalation pharmacokinetics of benzene, toluene, m-xylene, o-xylene, p-xylene, ethylbenzene, dichloromethane, trichloroethylene, tetrachloroethylene and styrene in various mixtures. Overall, the present study indicates the potential usefulness of the QSAR-PBPK modelling approach to provide first-cut evaluations of the kinetics of chemicals in mixtures of increasing complexity, on the basis of chemical structure.

Prediction of skin sensitization with a particle swarm optimized support vector machine

Skin sensitization is the most commonly reported occupational illness, causing much suffering to a wide range of people. Identification and labeling of environmental allergens is urgently required to protect people from skin sensitization. The guinea pig maximization test (GPMT) and murine local lymph node assay (LLNA) are the two most important in vivo models for identification of skin sensitizers. In order to reduce the number of animal tests, quantitative structure-activity relationships (QSARs) are strongly encouraged in the assessment of skin sensitization of chemicals. This paper has investigated the skin sensitization potential of 162 compounds with LLNA results and 92 compounds with GPMT results using a support vector machine. A particle swarm optimization algorithm was implemented for feature selection from a large number of molecular descriptors calculated by Dragon. For the LLNA data set, the classification accuracies are 95.37% and 88.89% for the training and the test sets, respectively. For the GPMT data set, the classification accuracies are 91.80% and 90.32% for the training and the test sets, respectively. The classification performances were greatly improved compared to those reported in the literature, indicating that the support vector machine optimized by particle swarm in this paper is competent for the identification of skin sensitizers.