CADRE-SS, an in Silico Tool for Predicting Skin Sensitization Potential Based on Modeling of Molecular Interactions

Tale of Three N-Nitrosamines and the Variables Needed to Assess Their Carcinogenicity In Silico Incorporated into a Single Workflow

N-Nitrosamine impurities in pharmaceuticals present a considerable challenge for regulators and industry alike, where the absence of carcinogenic-potency studies has left a gap that must be adequately filled to protect public health. In the interim, this means balancing risk assessment with the necessity to continue research, development, and supply of pharmaceuticals. In the long term, we need a cost-effective solution that optimizes both. As if beholden to Newton's Third Law, every crisis breeds an opportunity of equal magnitude. Consequently, cross-industry consortia have been racing to find a solution by advancing our current science. Recent spotlight has been on in silico tools, as a fast and increasingly reliable alternative to in vivo and in vitro testing. Because N-nitrosamine bioactivation lends itself uniquely to quantum mechanics (QM) approaches, the integration of electronic-structure considerations has emerged as the dominant in silico approach. This signifies a considerable leap in predictive toxicology, which has, for much of its existence, relied on atomistic (quantitative) structure-activity relationships, i.e., (Q)SARs. Here we present a validation of an integrated docking-QM approach within the CADRE program and demonstrate its utility on three different impurities, N-nitroso-7-monomethylamino-6-deoxytetracycline, N-nitroso-dabigatran etexilate, and 1-methyl-4-nitrosopiperazine. We show that a combined in silico strategy, which considers bioavailability, transport, cytochrome P450 binding, and reactivity, can be leveraged to supplement the overly conservative Carcinogenic Potency Categorization Approach (CPCA) in setting the daily acceptable intake (AI) using defensible, highly mechanistic, and quantitative drivers of N-nitrosamine metabolism. To that end, we argue that while N-nitroso-7-monomethylamino-6-deoxytetracycline and 1-methyl-4-nitrosopiperazine are cohort-of-concern impurities, N-nitroso-dabigatran etexilate is not a potent carcinogen (TD50 > 1.5 mg/kg/day), contrasting the CPCA-derived AI. Lastly, we discuss how the CADRE tool can be integrated with the broader landscape of QM methods and the CPCA into a single harmonized in silico strategy for carcinogenicity assessment of N-nitrosamine impurities.

Capturing Differential Quality of Experimental Evidence in a Predictive Quantum-Mechanical Model for Respiratory Sensitization

Asthma is of concern in occupational toxicology with significant public-health and economic costs. In the absence of benchmark in vivo and in vitro tests, the use of mechanistically sound in silico models is critical to inform hazard and to protect workers from exposure to potentially harmful substances. We recently reported on the computer-aided discovery and REdesign (CADRE) model for respiratory sensitization, which relies on a tiered structure of expert rules, molecular simulations, quantum-mechanics calculations and advanced statistics to accurately identify respiratory sensitizers from first principles. Here, we present an update to this model based on two years of testing in the pharmaceutical space, where we captured the heterogeneity of the underlying experimental evidence in two predictive tiers, thus allowing the practitioner to select an outcome based on their expert assessment of the data reliability and relevance. This user-based tuning of predictive models is critical for end points that lack consensus on what constitutes satisfactory evidence to support a decision in the handling of chemicals for occupational safety.

Molecular similarity in chemical informatics and predictive toxicity modeling: from quantitative read-across (q-RA) to quantitative read-across structure-activity relationship (q-RASAR) with the application of machine learning

This article aims to provide a comprehensive critical, yet readable, review of general interest to the chemistry community on molecular similarity as applied to chemical informatics and predictive modeling with a special focus on read-across (RA) and read-across structure-activity relationships (RASAR). Molecular similarity-based computational tools, such as quantitative structure-activity relationships (QSARs) and RA, are routinely used to fill the data gaps for a wide range of properties including toxicity endpoints for regulatory purposes. This review will explore the background of RA starting from how structural information has been used through to how other similarity contexts such as physicochemical, absorption, distribution, metabolism, and elimination (ADME) properties, and biological aspects are being characterized. More recent developments of RA's integration with QSAR have resulted in the emergence of novel models such as ToxRead, generalized read-across (GenRA), and quantitative RASAR (q-RASAR). Conventional QSAR techniques have been excluded from this review except where necessary for context.

Quantum-Mechanics Calculations Elucidate Skin-Sensitizing Pharmaceutical Compounds

Skin sensitization is a critical end point in occupational toxicology that necessitates the use of fast, accurate, and affordable models to aid in establishing handling guidance for worker protection. While many in silico models have been developed, the scarcity of reliable data for active pharmaceutical ingredients (APIs) and their intermediates (together regarded as pharmaceutical compounds) brings into question the reliability of these tools, which are largely constructed using publicly available nonspecialty chemicals. Here, we present the quantum-mechanical (QM) Computer-Aided Discovery and REdesign (CADRE) model, which was developed with the bioactive and structurally complex chemical space in mind by relying on the fundamentals of chemical interactions in key events (versus structural attributes of training-set data). Validated in this study on 345 APIs and intermediates, CADRE achieved 95% accuracy, sensitivity, and specificity and a combined 79% accuracy in assigning potency categories compared to the mouse local lymph node assay data. We show how historical outcomes from CADRE testing in the pharmaceutical space, generated over the past 10 years on ca. 2500 chemicals, can be used to probe the relationships between sensitization mechanisms (or the underlying chemical classes) and the probability of eliciting a sensitization response in mice of a given potency. We believe this information to be of value to both practitioners, who can use it to quickly screen and triage their data sets, as well as to model developers to fine-tune their structure-based tools. Lastly, we leverage our experimentally validated subset of APIs and intermediates to show the importance of dermal permeability on the sensitization potential and potency. We demonstrate that common physicochemical properties used to assess permeation, such as the octanol-water partition coefficient and molecular weight, are poor proxies for the more accurate energy-pair distributions that can be computed from mixed QM and classical simulations using model representations of the stratum corneum.

Computational screening of pathogenic missense nsSNPs in heme oxygenase 1 (HMOX1) gene and their structural and functional consequences

Heme Oxygenase 1 (HMOX1) is a cytoprotective enzyme, exhibiting the highest activity in the spleen, catalyzing the heme ring breakdown into products of biological significance- biliverdin, CO, and Fe2+. In vascular cells, HMOX1 possesses strong anti-apoptotic, antioxidant, anti-proliferative, anti-inflammatory, and immunomodulatory actions. The majority of these activities are crucial for the prevention of atherogenesis. Single amino acid substitutions in proteins generated by missense non-synonymous single nucleotide polymorphism (nsSNPs) in the protein-encoding regions of genes are potent enough to cause significant medical challenges due to the alteration of protein structure and function. The current study aimed at characterizing and analyzing high-risk nsSNPs associated with the human HMOX1 gene. Preliminary screening of the total available 288 missense SNPs was performed through the lens of deleteriousness and stability prediction tools. Finally, a total of seven nsSNPs (Y58D, A131T, Y134H, F166S, F167S, R183S and M186V) were found to be most deleterious by all tools that are present at highly conserved positions. Molecular dynamics simulations (MDS) analysis explained the mutational effects on the dynamic action of the wild-type and mutant proteins. In a nutshell, R183S (rs749644285) was identified as a highly detrimental mutation that could significantly render the enzymatic activity of HMOX1. The finding of this computational analysis might help subject the experimental confirmatory analysis to characterize the role of nsSNPs in HMOX1.Communicated by Ramaswamy H. Sarma.

Low-cost quantum mechanical descriptors for data efficient skin sensitization QSAR models

Quantitative Structure Activity Relationship modelling methodologies need to incorporate relevant mechanistic information to have high predictive performance and validity. Electrophilic reactivity is a common mechanistic feature of skin sensitization endpoints which could be concisely characterized with electronic descriptors which is key to enabling the modelling of small datasets in this domain. However, quantum mechanical methodologies have previously featured high computational costs which would exclude the use of large datasets. Consequently, we investigate the use of electronic descriptors calculated using the Hartree Fock with 3 corrections (Hf-3c) method, a low-cost ab initio methodology that has higher chemical accuracy than previous semiempirical methodologies for modelling in vitro skin sensitization assay outcomes. We also model the Ames assay as a surrogate for determining skin sensitization outcomes. The quantum chemical descriptors calculated using the Hf-3c method with conductor-like polarizable continuum model (CPCM) implicit solvation found improved QSAR model performance for the in vitro Ames (n = 6049, 0.770 AUC), KeratinoSens (n = 164, 0.763 AUC), and Direct Peptide Reactivity Assay (n = 122, 0.750 AUC) datasets, with their combination producing high predictive performance for unseen in vivo Local Lymph Node Assay (n = 86, 0.789 AUC) and Human Repeated Insult Patch Test (n = 86, 0.791 AUC) assay toxicant outcomes.

Development of a hierarchical support vector regression-based in silico model for the prediction of the cysteine depletion in DPRA

Topical and transdermal treatments have been dramatically growing recently and it is crucial to consider skin sensitization during the drug discovery and development process for these administration routes. Various tests, including animal and non-animal approaches, have been devised to assess the potential for skin sensitization. Furthermore, numerous in silico models have been created, providing swift and cost-effective alternatives to traditional methods such as in vivo, in vitro, and in chemico methods for categorizing compounds. In this study, a quantitative structure-activity relationship (QSAR) model was developed using the innovative hierarchical support vector regression (HSVR) scheme. The aim was to quantitatively predict the potential for skin sensitization by analyzing the percent of cysteine depletion in Direct Peptide Reactivity Assay (DPRA). The results demonstrated accurate, consistent, and robust predictions in the training set, test set, and outlier set. Consequently, this model can be employed to estimate skin sensitization potential of novel or virtual compounds.

Comprehensive analysis of non-synonymous missense SNPs of human galactose mutarotase (GALM) gene: an integrated computational approach

Missense Non-synonymous single nucleotide polymorphisms (nsSNPs) of Galactose Mutarotase (GALM) are associated with the Novel type of Galactosemia (Galactosemia type 4) together with symptoms such as high blood galactose levels and eye cataracts. The objective of the present study was to identify deleterious nsSNPs of GALM recorded on the dbSNP database through comprehensive insilico analysis. Among the 319 missense nsSNPs reported, various insilco tools predicted R78S, R82G, A163E, P210S, Y281C, E307G and F339C as the most deleterious mutations. Structural analysis, PTM analysis and molecular dynamics simulations (MDS) were carried out to understand the effect of these mutations on the structural and physicochemical properties of the GALM protein. The residues R82G and E307G were found to be part of the binding site that resulted in decreased surface accessibility. Replacing the charged wild-type residue with a neutral mutant type affected its substrate binding. All 7 mutations were found to increase the rigidity of the protein structure, which is unfavorable during ligand binding. The mutation F339E made the protein structure more rigid than all the other mutations. Y281 is a phosphorylated site, and therefore, less significant structural changes were observed when compared to other mutations; however, it may have significant differences in the usual functioning of the protein. In summary, the structural and functional analysis of missense SNPs of GALM is important to reduce the number of potential mutations to be evaluated in vitro to understand the association with some genetic diseases.Communicated by Ramaswamy H. Sarma.

Structural characterization and in silico toxicity prediction of degradation impurities of roxadustat

Roxadustat is the first drug approved for anemia due to chronic kidney disease. Drug degradation profile is very crucial for assessing the quality and safety of the drug substances and their formulations. Forced degradation studies are conducted for quick prediction of drug degradation products. Forced degradation of roxadustat was carried out as per ICH guidelines, and nine degradation products (DPs) were observed. These DPs (DP-1 to DP-9) were separated using the reverse phase HPLC gradient method with an XBridge column (250 mm × 4.6 mm, 5 µm). The mobile phase consisted of 0.1% formic acid (solvent A) and acetonitrile (solvent B) at a flow rate of 1.0 ml/min. The chemical structures of all the DPs were proposed by using LC-Q-TOF/MS. DP-4 and DP-5, the two major degradation impurities, were isolated, and NMR was used to confirm their chemical structures. Based on our experiments, the roxadustat was found stable to thermal degradation in solid state and oxidative conditions. However, it was unstable in acidic, basic, and photolytic conditions. A very remarkable observation was made about DP-4 impurity. DP-4 was generated as a common degradation impurity in alkaline hydrolysis, neutral hydrolysis as well as photolysis conditions. DP-4 has a similar molecular mass to roxadustat but is structurally different. DP-4 is chemically, (1a-methyl-6-oxo-3-phenoxy-1,1a,6,6a-tetrahydroindeno [1,2-b] aziridine-6a-carbonyl) glycine. In silico toxicity study was conducted using Dereck software to gain the best knowledge of the drug and its degradation products towards carcinogenicity, mutagenicity, teratogenicity, and skin sensitivity. A further study using molecular docking confirmed the potential interaction of DPs with proteins responsible for toxicity. DP-4 shows a toxicity alert due to the presence of aziridine moiety.

Making the Case for Quantum Mechanics in Predictive Toxicology─Nearly 100 Years Too Late?

The use of quantum mechanics (QM) has long been the norm to study covalent-binding phenomena in chemistry and biochemistry. The pharmaceutical industry leverages QM models explicitly in covalent drug discovery and implicitly to characterize short-range interactions in noncovalent binding. Predictive toxicology has resisted widespread adoption of QM, including in the pharmaceutical industry, despite its obvious relevance to the metabolic processes in the upstream of adverse outcome pathways and advances in both QM methods and computational resources, which support fit-for-purpose applications in reasonable timeframes. Here, we make the case for embracing QM as an indispensable part of a toxicologist's toolkit. We argue that QM provides the necessary orthogonality to alert-based expert systems and traditional QSARs, consistent with calls for animal-free integrated testing strategies for safety assessments of commercial chemicals. We outline existing roadblocks to this transition, including the need to train model developers in QM and the shift toward service-based toxicity models that utilize high-performance computing clusters. Lastly, we describe recent examples of successful implementations of QM in hazard assessments and propose how in silico toxicology can be further advanced by integrating QM with artificial intelligence.

Toward Property-Based Regulation

An expanding web of adverse impacts on people and the environment has been steadily linked to anthropogenic chemicals and their proliferation. Central to this web are the regulatory structures intended to protect human and environmental health through the control of new molecules. Through chronically insufficient and inefficient action, the current chemical-by-chemical regulatory approach, which considers regulation at the level of chemical identity, has enabled many adverse impacts to develop and persist. Recognizing the link between fundamental physicochemical properties and hazards, we describe a new paradigm─property-based regulation. By regulating physicochemical properties, we show how governments can delineate and enforce safe chemical spaces, increasing the scalability of chemical assessments, reducing the time and resources to regulate a substance, and providing transparency for chemical designers. We highlight sparse existing property-based approaches and demonstrate their applicability using bioaccumulation as an example. Finally, we present a path to implementation in the United States, prescribing roles and steps for government, nongovernmental organizations, and industry to accelerate this transition, to the benefit of all.

The Good, The Bad, and The Perplexing: Structural Alerts and Read-Across for Predicting Skin Sensitization Using Human Data

In our earlier work (Golden et al., 2021), we showed 70-80% accuracies for several skin sensitization computational tools using human data. Here, we expanded the data set using the NICEATM human skin sensitization database to create a final data set of 1355 discrete chemicals (largely negative, ∼70%). Using this expanded data set, we analyzed model performance and evaluated mispredictions using Toxtree (v 3.1.0), OECD QSAR Toolbox (v 4.5), VEGA's (1.2.0 BETA) CAESAR (v 2.1.7), and a k-nearest-neighbor (kNN) classification approach. We show that the accuracy on this data set was lower than previous estimates, with balanced accuracies being 63% and 65% for Toxtree and OECD QSAR Toolbox, respectively, 46% for VEGA, and 59% for a kNN approach, with the lower accuracy likely due to the higher percentage of nonsensitizing chemicals. Two hundred eighty seven chemicals were mispredicted by both Toxtree and OECD QSAR Toolbox, which was approximately 20% of the entire data set, and 84% of these were false positives. The absence or presence of metabolic simulation in OECD QSAR Toolbox made no overall difference. While Toxtree is known for overpredicting, 60% of the chemicals in the data set had no alert for skin sensitization, and a substantial number of these chemicals were in fact sensitizers, pointing to sensitization mechanisms not recognized by Toxtree. Interestingly, we observed that chemicals with more than one Toxtree alert were more likely to be nonsensitizers. Finally, a kNN approach tended to mispredict different chemicals than either OECD QSAR Toolbox or Toxtree, suggesting that there was additional information to be garnered from a kNN approach. Overall, the results demonstrate that while there is merit in structural alerts as well as QSAR or read-across approaches (perhaps even more so in their combination), additional improvement will require a more nuanced understanding of mechanisms of skin sensitization.

Quantum-Mechanical Approach to Predicting the Carcinogenic Potency of N-Nitroso Impurities in Pharmaceuticals

N-Nitroso contaminants in medicinal products are of concern due to their high carcinogenic potency; however, not all these compounds are created equal, and some are relatively benign chemicals. Understanding the structure-activity relationships (SARs) that drive hazards in one molecule versus another is key to both protecting human health and alleviating costly and sometimes inaccurate animal testing. Here, we report on an extension of the CADRE (computer-aided discovery and REdesign) platform, which is used broadly by the pharmaceutical and personal care industries to assess environmental and human health endpoints, to predict the carcinogenic potency of N-nitroso compounds. The model distinguishes compounds in three potency categories with 77% accuracy in external testing, which surpasses the reproducibility of rodent cancer bioassays and constraints imposed by limited (high-quality) data. The robustness of predictions for more complex pharmaceuticals is maximized by capturing key SARs using quantum mechanics, that is, by hinging the model on the underlying chemistry versus chemicals in the training set. To this end, the present approach can be leveraged in a quantitative hazard assessment and to offer qualitative guidance using electronic structure comparisons between well-studied analogues and unknown contaminants.

Computer-Aided Discovery and Redesign for Respiratory Sensitization: A Tiered Mechanistic Model to Deliver Robust Performance Across a Diverse Chemical Space

Asthma is among the most common occupational diseases with considerable public health and economic costs. Chemicals that induce hypersensitivity in the airways can cause respiratory distress and comorbidities with respiratory infections such as COVID. Robust predictive models for this end point are still elusive due to the lack of an experimental benchmark and the over-reliance of existing in silico tools on structural alerts and structural (vs chemical) similarities. The Computer-Aided Discovery and REdesign (CADRE) platform is a proven strategy for providing robust computational predictions for hazard end points using a tiered hybrid system of expert rules, molecular simulations, and quantum mechanics calculations. The recently developed CADRE model for respiratory sensitization is based on a highly curated data set of structurally diverse chemicals with high-fidelity biological data. The model evaluates absorption kinetics in lung mucosa using Monte Carlo simulations, assigns reactive centers in a molecule and possible biotransformations via expert rules, and determines subsequent reactivity with cell proteins via quantum-mechanics calculations using a multi-tiered regression. The model affords an accuracy above 0.90, with a series of external validations based on literature data in the range of 0.88-0.95. The model is applicable to all low-molecular-weight organics and can inform not only chemical substitution but also chemical redesign to advance development of safer alternatives.

Framework for sensitization assessment of extractables and leachables in pharmaceuticals

During the toxicological assessment of extractables and leachables in drug products, localized hazards such as irritation or sensitization may be identified. Typically, because of the low concentration at which leachables occur in pharmaceuticals, irritation is of minimal concern; therefore, this manuscript focuses on sensitization potential. The primary objective of performing a leachable sensitization assessment is protection against Type IV induction of sensitization, rather than prevention of an elicitation response, as it is not possible to account for the immunological state of every individual. Sensitizers have a wide range of potencies and those which induce sensitization upon exposure at a low concentration (i.e. strong, or extreme sensitizers) pose the highest risk to patients and should be the focus of the risk assessment. The Extractables and Leachables Safety Information Exchange (ELSIE) consortium has reviewed the status of dermal, respiratory, and systemic risk assessment in cosmetic and pharmaceutical industries, and proposes a framework to evaluate the safety of known or potential dermal sensitizers in pharmaceuticals. Due to the lack of specific regulatory guidance on this topic, the science-driven risk-based approach proposed by ELSIE encourages consistency in the toxicological assessment of extractables and leachables to maintain high product quality and ensure patient safety.

An Evaluation of the Occupational Health Hazards of Peptide Couplers

Peptide couplers (also known as amide bond-forming reagents or coupling reagents) are broadly used in organic chemical syntheses, especially in the pharmaceutical industry. Yet, occupational health hazards associated with this chemical class are largely unexplored, which is disconcerting given the intrinsic reactivity of these compounds. Several case studies involving occupational exposures reported adverse respiratory and dermal health effects, providing initial evidence of chemical sensitization. To address the paucity of toxicological data, a pharmaceutical cross-industry task force was formed to evaluate and assess the potential of these compounds to cause eye and dermal irritation as well as corrosivity and dermal sensitization. The goal of our work was to inform health and safety professionals as well as pharmaceutical and organic chemists of the occupational health hazards associated with this chemical class. To that end, 25 of the most commonly used peptide couplers and five hydrolysis products were selected for in vivo, in vitro, and in silico testing. Our findings confirmed that dermal sensitization is a concern for this chemical class with 21/25 peptide couplers testing positive for dermal sensitization and 15 of these being strong/extreme sensitizers. We also found that dermal corrosion and irritation (8/25) as well as eye irritation (9/25) were health hazards associated with peptide couplers and their hydrolysis products (4/5 were dermal irritants or corrosive and 4/5 were eye irritants). Resulting outcomes were synthesized to inform decision making in peptide coupler selection and enable data-driven hazard communication to workers. The latter includes harmonized hazard classifications, appropriate handling recommendations, and accurate safety data sheets, which support the industrial hygiene hierarchy of control strategies and risk assessment. Our study demonstrates the merits of an integrated, in vivo -in silico analysis, applied here to the skin sensitization endpoint using the Computer-Aided Discovery and REdesign (CADRE) and Derek Nexus programs. We show that experimental data can improve predictive models by filling existing data gaps while, concurrently, providing computational insights into key initiating events and elucidating the chemical structural features contributing to adverse health effects. This interactive, interdisciplinary approach is consistent with Green Chemistry principles that seek to improve the selection and design of less hazardous reagents in industrial processes and applications.

Structure-Energetics-Property Relationships Support Computational Design of Photodegradable Pesticides

Development of high-performing pesticides with tunable degradation properties is vital to increasing the safety and effectiveness of tomorrow's analogs. Chromophoric dissolved organic matter in the excited triple state (³CDOM*) is known to play a key role in the removal of pesticides via indirect photodegradation. However, the potential of these transformations to guide the design of safer chemicals has not yet been fully realized. Here, we report a two-tier computational framework developed to probe and predict both kinetics and thermodynamics of ³CDOM*-pesticide interactions. In the first tier, robust in silico models were constructed by fitting free energies obtained from density functional theory (DFT) calculations to cell potentials and second-order rate constants for the ³CDOM*-pesticide electron transfer. In the second tier, Gibbs free energies and corresponding free energy barriers, determined in solution using the Marcus theory, were applied to develop a quick yet accurate screening approach based on the frontier molecular orbital (FMO) Theory. Being highly mechanistic and spanning ca. 1500 unique ³CDOM*-pesticide interactions, our approach is both robust and broadly applicable. To that end, the outcomes of our computational models were integrated into an easy-to-use decision framework that can guide structure-based design of less persistent pesticide analogs.

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.

Going All In: A Strategic Investment in In Silico Toxicology

As vast numbers of new chemicals are introduced to market annually, we are faced with the grand challenge of protecting humans and the environment while minimizing economically and ethically costly animal testing. In silico models promise to be the solution we seek, but we find ourselves at crossroads of future development efforts that would ensure standalone applicability and reliability of these tools. A conscientious effort that prioritizes experimental testing to support the needs of in silico models (versus regulatory needs) is called for to achieve this goal. Using economic analogy in the title of this work, we argue that a prudent investment is to go all-in to support in silico model development, rather than gamble our future by keeping the status quo of a "balanced portfolio" of testing approaches. We discuss two paths to future in silico toxicology-one based on big-data statistics ("broadsword"), and the other based on direct modeling of molecular interactions ("scalpel")-and offer rationale that the latter approach is more transparent, is better aligned with our quest for fundamental knowledge, and has a greater potential to succeed if we are willing to transform our toxicity-testing paradigm.

In silico prediction of toxicity and its applications for chemicals at work

OBJECTIVE AND METHODS

This study reviewed the concept of in silico prediction of chemical toxicity for prevention of occupational cancer and future prospects in workers' health. In this review, a new approach to determine the credibility of in silico predictions with raw data is explored, and the method of determining the confidence level of evaluation based on the credibility of data is discussed. I searched various papers and books related to the in silico prediction of chemical toxicity and carcinogenicity. The intention was to utilize the most recent reports after 2015 regarding in silico prediction.

RESULTS AND CONCLUSION

The application of in silico methods is increasing with the prediction of toxic risks to human and the environment. The various toxic effects of industrial chemicals have triggered the recognition of the importance of using a combination of in silico models in the risk assessments. In silico occupational exposure models, industrial accidents, and occupational cancers are effectively managed and chemicals evaluated. It is important to identify and manage hazardous substances proactively through the rigorous evaluation of chemicals.

Toward Less Hazardous Industrial Compounds: Coupling Quantum Mechanical Computations, Biomarker Responses, and Behavioral Profiles To Identify Bioactivity of SN2 Electrophiles in Alternative Vertebrate Models

Sustainable molecular design of less hazardous chemicals promises to reduce risks to public health and the environment. Computational chemistry modeling coupled with alternative toxicology models (e.g., larval fish) present unique high-throughput opportunities to understand structural characteristics eliciting adverse outcomes. Numerous environmental contaminants with reactive properties can elicit oxidative stress, an important toxicological response associated with diverse adverse outcomes (i.e., cancer, diabetes, neurodegenerative disorders, etc.). We examined a common chemical mechanism (bimolecular nucleophilic substitution (S_N2)) associated with oxidative stress using property-based computational modeling coupled with acute (mortality) and sublethal (glutathione, photomotor behavior) responses in the zebrafish (Danio rerio) and the fathead minnow (Pimephales promelas) models to identify whether relationships exist among biological responses and molecular attributes of industrial chemicals. Following standardized methods, embryonic zebrafish and larval fathead minnows were exposed separately to eight different S_N2 compounds for 96 h. Acute and sublethal responses were compared to computationally derived in silico chemical descriptors. Specifically, frontier molecular orbital energies were significantly related to acute LC₅₀ values and photomotor response (PMR) no observed effect concentrations (NOECs) in both fathead minnow and zebrafish. This reactivity index, LC₅₀ values, and PMR NOECs were also significantly related to whole body glutathione (GSH) levels, suggesting that acute and chronic toxicity results from protein adduct formation for S_N2 electrophiles. Shared refractory locomotor response patterns among study compounds and two alternative vertebrate models appear informative of electrophilic properties associated with oxidative stress for S_N2 chemicals. Electrophilic parameters derived from frontier molecular orbitals were predictive of experimental in vivo acute and sublethal toxicity. These observations provide important implications for identifying and designing less hazardous industrial chemicals with reduced potential to elicit oxidative stress through bimolecular nucleophilic substitution.

Towards Sustainable Environmental Quality: Priority Research Questions for the Australasian Region of Oceania

Environmental challenges persist across the world, including the Australasian region of Oceania, where biodiversity hotspots and unique ecosystems such as the Great Barrier Reef are common. These systems are routinely affected by multiple stressors from anthropogenic activities, and increasingly influenced by global megatrends (e.g., the food-energy-water nexus, demographic transitions to cities) and climate change. Here we report priority research questions from the Global Horizon Scanning Project, which aimed to identify, prioritize, and advance environmental quality research needs from an Australasian perspective, within a global context. We employed a transparent and inclusive process of soliciting key questions from Australasian members of the Society of Environmental Toxicology and Chemistry. Following submission of 78 questions, 20 priority research questions were identified during an expert workshop in Nelson, New Zealand. These research questions covered a range of issues of global relevance, including research needed to more closely integrate ecotoxicology and ecology for the protection of ecosystems, increase flexibility for prioritizing chemical substances currently in commerce, understand the impacts of complex mixtures and multiple stressors, and define environmental quality and ecosystem integrity of temporary waters. Some questions have specific relevance to Australasia, particularly the uncertainties associated with using toxicity data from exotic species to protect unique indigenous species. Several related priority questions deal with the theme of how widely international ecotoxicological data and databases can be applied to regional ecosystems. Other timely questions, which focus on improving predictive chemistry and toxicology tools and techniques, will be important to answer several of the priority questions identified here. Another important question raised was how to protect local cultural and social values and maintain indigenous engagement during problem formulation and identification of ecosystem protection goals. Addressing these questions will be challenging, but doing so promises to advance environmental sustainability in Oceania and globally.

A Fast Ab Initio Predictor Tool for Covalent Reactivity Estimation of Acrylamides

Thanks to their unique mode of action, covalent drugs represent an exceptional opportunity for drug design. After binding to a biologically relevant target system, covalent compounds form a reversible or irreversible covalent bond with a nucleophilic amino acid. Due to the inherently large binding energy of a covalent bond, covalent binders exhibit higher potencies and thus allow potentially lower drug dosages. However, a proper balancing of compound reactivity is key for the design of covalent binders, to achieve high levels of target inhibition while minimizing promiscuous covalent binding to nontarget proteins. In this work, we demonstrated the possibility to apply the electrophilicity index concept to estimate covalent compound reactivity. We tested this approach on acrylamides, one of the most prominent classes of covalent warheads. Our study clearly demonstrated that, for compounds with molecular weight (MW) below 250 Da, the electrophilicity index can be directly used to estimate compound reactivity. On the other hand, for leadlike molecules (MW > 250 Da) we implemented a new truncation algorithm that has to be applied before reactivity calculations. This algorithm can ensure the localization of HOMO/LUMO orbitals on the compound warhead and thus a correct estimation of its reactivity. Our results also indicate that caution should be used when employing the electrophilicity index to estimate the reactivity of nonterminal acrylamides. The nonparametric nature of this method and its reasonable computational cost make it a suitable tool to support covalent drug design.

Kinetics of Glutathione Depletion and Antioxidant Gene Expression as Indicators of Chemical Modes of Action Assessed in Vitro in Mouse Hepatocytes with Enhanced Glutathione Synthesis

Here we report a vertically integrated in vitro - in silico study that aims to elucidate the molecular initiating events involved in the induction of oxidative stress (OS) by seven diverse chemicals (cumene hydroperoxide, t-butyl hydroperoxide, hydroquinone, t-butyl hydroquinone, bisphenol A, Dinoseb, and perfluorooctanoic acid). To that end, we probe the relationship between chemical properties, cell viability, glutathione (GSH) depletion, and antioxidant gene expression. Concentration-dependent effects on cell viability were assessed by MTT assay in two Hepa-1 derived mouse liver cell lines: a control plasmid vector transfected cell line (Hepa-V), and a cell line with increased glutamate-cysteine ligase (GCL) activity and GSH content (CR17). Changes to intracellular GSH content and mRNA expression levels for the Nrf2-driven antioxidant genes Gclc, Gclm, heme oxygenase-1 ( Hmox1), and NADPH quinone oxidoreductase-1 ( Nqo1) were monitored after sublethal exposure to the chemicals. In silico models of covalent and redox reactivity were used to rationalize differences in activity of quinones and peroxides. Our findings show CR17 cells were generally more resistant to chemical toxicity and showed markedly attenuated induction of OS biomarkers; however, differences in viability effects between the two cell lines were not the same for all chemicals. The results highlight the vital role of GSH in protecting against oxidative stress-inducing chemicals as well as the importance of probing molecular initiating events in order to identify chemicals with lower potential to cause oxidative stress.

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.