Extending (Q)SARs to incorporate proprietary knowledge for regulatory purposes: A case study using aromatic amine mutagenicity

Local QSAR based on quantum chemistry calculations for the stability of nitrenium ions to reduce false positive outcomes from standard QSAR systems for the mutagenicity of primary aromatic amines

BACKGROUND

Primary aromatic amines (PAAs) present significant challenges in the prediction of mutagenicity using current standard quantitative structure activity relationship (QSAR) systems, which are knowledge-based and statistics-based, because of their low positive prediction values (PPVs). Previous studies have suggested that PAAs are metabolized into genotoxic nitrenium ions. Moreover, ddE, a relative-energy based index derived from quantum chemistry calculations that measures the stability nitrenium ions, has been correlated with mutagenicity. This study aims to further examine the ability of the ddE-based approach in improving QSAR mutagenicity predictions for PAAs and to develop a refined method to decrease false positive predictions.

RESULTS

Information on 1,177 PAAs was collected, of which 420 were from public databases and 757 were from in-house databases across 16 laboratories. The total dataset included 465 Ames test-positive and 712 test-negative chemicals. For internal PAAs, detailed Ames test data were scrutinized and final decisions were made using common evaluation criteria. In this study, ddE calculations were performed using a convenient and consistent protocol. An optimal ddE cutoff value of -5 kcal/mol, combined with a molecular weight ≤ 500 and ortho substitution groups yielded well-balanced prediction scores: sensitivity of 72.0%, specificity of 75.9%, PPV of 65.6%, negative predictive value of 80.9% and a balanced accuracy of 74.0%. The PPV of the ddE-based approach was greatly reduced by the presence of two ortho substituent groups of ethyl or larger, as because almost all of them were negative in the Ames test regardless of their ddE values, probably due to steric hindrance affecting interactions between the PAA and metabolic enzymes. The great majority of the PAAs whose molecular weights were greater than 500 were also negative in Ames test, despite ddE predictions indicating positive mutagenicity.

CONCLUSIONS

This study proposes a refined approach to enhance the accuracy of QSAR mutagenicity predictions for PAAs by minimizing false positives. This integrative approach incorporating molecular weight, ortho substitution patterns, and ddE values, substantially can provide a more reliable basis for evaluating the genotoxic potential of PAAs.

Computational Strategies for Assessing Adverse Outcome Pathways: Hepatic Steatosis as a Case Study

The evolving landscape of chemical risk assessment is increasingly focused on developing tiered, mechanistically driven approaches that avoid the use of animal experiments. In this context, adverse outcome pathways have gained importance for evaluating various types of chemical-induced toxicity. Using hepatic steatosis as a case study, this review explores the use of diverse computational techniques, such as structure-activity relationship models, quantitative structure-activity relationship models, read-across methods, omics data analysis, and structure-based approaches to fill data gaps within adverse outcome pathway networks. Emphasizing the regulatory acceptance of each technique, we examine how these methodologies can be integrated to provide a comprehensive understanding of chemical toxicity. This review highlights the transformative impact of in silico techniques in toxicology, proposing guidelines for their application in evidence gathering for developing and filling data gaps in adverse outcome pathway networks. These guidelines can be applied to other cases, advancing the field of toxicological risk assessment.

Mutagenicity assessment of two potential impurities in preparations of 5-amino-2,4,6 triiodoisophthalic acid, a key intermediate in the synthesis of the iodinated contrast agent iopamidol

5-Aminoisophthalic acid and 5-nitroisophthalic acid (5-NIPA) are potential impurities in preparations of 5-amino-2,4,6-triiodoisophthalic acid, which is a key intermediate in the synthesis of the iodinated contrast agent iopamidol. We have studied their mutagenicity in silico (quantitative structure-activity relationships, QSAR) and by the bacterial reverse mutation assay (Ames test). First, the compounds were screened with the tools Derek Nexus™ and Leadscope®. Both compounds were flagged as potentially mutagenic (class 3 under ICH M7). However, contrary to the in silico prediction, neither chemical was mutagenic in the Ames test (plate incorporation method) with or without S9 metabolic activation.

A local QSAR model based on the stability of nitrenium ions to support the ICH M7 expert review on the mutagenicity of primary aromatic amines

BACKGROUND

Aromatic amines, often used as intermediates for pharmaceutical synthesis, may be mutagenic and therefore pose a challenge as metabolites or impurities in drug development. However, predicting the mutagenicity of aromatic amines using commercially available, quantitative structure-activity relationship (QSAR) tools is difficult and often requires expert review. In this study, we developed a shareable QSAR tool based on nitrenium ion stability.

RESULTS

The evaluation using in-house aromatic amine intermediates revealed that our model has prediction accuracy of aromatic amine mutagenicity comparable to that of commercial QSAR tools. The effect of changing the number and position of substituents on the mutagenicity of aromatic amines was successfully explained by the change in the nitrenium ion stability. Furthermore, case studies showed that our QSAR tool can support the expert review with quantitative indicators.

CONCLUSIONS

This local QSAR tool will be useful as a quantitative support tool to explain the substituent effects on the mutagenicity of primary aromatic amines. By further refinement through method sharing and standardization, our tool can support the International Council for Harmonisation of Technical Requirements for Pharmaceuticals for Human Use (ICH) M7 expert review with quantitative indicators.

Implementation of in silico toxicology protocols within a visual and interactive hazard assessment platform

Mechanistically-driven alternative approaches to hazard assessment invariably require a battery of tests, including both in silico models and experimental data. The decision-making process, from selection of the methods to combining the information based on the weight-of-evidence, is ideally described in published guidelines or protocols. This ensures that the application of such approaches is defendable to reviewers within regulatory agencies and across the industry. Examples include the ICH M7 pharmaceutical impurities guideline and the published in silico toxicology protocols. To support an efficient, transparent, consistent and fully documented implementation of these protocols, a new and novel interactive software solution is described to perform such an integrated hazard assessment based on public and proprietary information.

Oxidation of Drugs during Drug Product Development: Problems and Solutions

Oxidation is the second most common degradation pathway for pharmaceuticals, after hydrolysis. However, in contrast to hydrolysis, oxidation is mechanistically more complex and produces a wider range of degradation products; oxidation is thus harder to control. The propensity of a drug towards oxidation is established during forced degradation studies. However, a more realistic insight into degradation in the solid state can be achieved with accelerated studies of mixtures of drugs and excipients, as the excipients are the most common sources of impurities that have the potential to initiate oxidation of a solid drug product. Based on the results of these studies, critical parameters can be identified and appropriate measures can be taken to avoid the problems that oxidation poses to the quality of a drug product. This article reviews the most common types of oxidation mechanisms, possible sources of reactive oxygen species, and how to minimize the oxidation of a solid drug product based on a well-planned accelerated study.

Artificial Intelligence in Drug Safety and Metabolism

The use of artificial intelligence methods in drug safety began in the early 2000s with applications such as predicting bacterial mutagenicity and hERG inhibition. The field has been endlessly expanding ever since and the models have become more complex. These approaches are now integrated into molecule risk assessment processes along with in vitro and in vivo methods. Today, artificial intelligence can be used in every phase of drug discovery and development, from profiling chemical libraries in early discovery, to predicting off-target effects in the mid-discovery phase, to assessing potential mutagenic impurities in development and degradants as part of life cycle management. This chapter provides an overview of artificial intelligence in drug safety and describes its application throughout the entire discovery and development process.

In Silico Approaches In Carcinogenicity Hazard Assessment: Current Status and Future Needs

Historically, identifying carcinogens has relied primarily on tumor studies in rodents, which require enormous resources in both money and time. In silico models have been developed for predicting rodent carcinogens but have not yet found general regulatory acceptance, in part due to the lack of a generally accepted protocol for performing such an assessment as well as limitations in predictive performance and scope. There remains a need for additional, improved in silico carcinogenicity models, especially ones that are more human-relevant, for use in research and regulatory decision-making. As part of an international effort to develop in silico toxicological protocols, a consortium of toxicologists, computational scientists, and regulatory scientists across several industries and governmental agencies evaluated the extent to which in silico models exist for each of the recently defined 10 key characteristics (KCs) of carcinogens. This position paper summarizes the current status of in silico tools for the assessment of each KC and identifies the data gaps that need to be addressed before a comprehensive in silico carcinogenicity protocol can be developed for regulatory use.

Assessing the impact of expert knowledge on ICH M7 (Q)SAR predictions. Is expert review still needed?

The ICH M7 (R1) guideline recommends the use of complementary (Q)SAR models to assess the mutagenic potential of drug impurities as a state-of-the-art, high-throughput alternative to empirical testing. Additionally, it includes a provision for the application of expert knowledge to increase prediction confidence and resolve conflicting calls. Expert knowledge, which considers structural analogs and mechanisms of activity, has been valuable when models return an indeterminate (equivocal) result or no prediction (out-of-domain). A retrospective analysis of 1002 impurities evaluated in drug regulatory applications between April 2017 and March 2019 assessed the impact of expert review on (Q)SAR predictions. Expert knowledge overturned the default predictions for 26% of the impurities and resolved 91% of equivocal predictions and 75% of out-of-domain calls. Of the 261 overturned default predictions, 15% were upgraded to equivocal or positive and 79% were downgraded to equivocal or negative. Chemical classes with the most overturns were primary aromatic amines (46%), aldehydes (45%), Michael-reactive acceptors (37%), and non-primary alkyl halides (33%). Additionally, low confidence predictions were the most often overturned. Collectively, the results suggest that expert knowledge continues to play an important role in an ICH M7 (Q)SAR prediction workflow and triaging predictions based on chemical class and probability can improve (Q)SAR review efficiency.

Mechanism-Based Insights into Removing the Mutagenicity of Aromatic Amines by Small Structural Alterations

Aromatic and heteroaromatic amines (ArNH₂) are activated by cytochrome P450 monooxygenases, primarily CYP1A2, into reactive N-arylhydroxylamines that can lead to covalent adducts with DNA nucleobases. Hereby, we give hands-on mechanism-based guidelines to design mutagenicity-free ArNH₂. The mechanism of N-hydroxylation of ArNH₂ by CYP1A2 is investigated by density functional theory (DFT) calculations. Two putative pathways are considered, the radicaloid route that goes via the classical ferryl-oxo oxidant and an alternative anionic pathway through Fenton-like oxidation by ferriheme-bound H₂O₂. Results suggest that bioactivation of ArNH₂ follows the anionic pathway. We demonstrate that H-bonding and/or geometric fit of ArNH₂ to CYP1A2 as well as feasibility of both proton abstraction by the ferriheme-peroxo base and heterolytic cleavage of arylhydroxylamines render molecules mutagenic. Mutagenicity of ArNH₂ can be removed by structural alterations that disrupt geometric and/or electrostatic fit to CYP1A2, decrease the acidity of the NH₂ group, destabilize arylnitrenium ions, or disrupt their pre-covalent transition states with guanine.

Current opinion on risk assessment of cosmetics

Risk assessment of cosmetic ingredients is a useful scientific method to characterize potential adverse effects resulting from using cosmetics. The process of risk assessment consists of four steps: hazard identification, dose-response assessment, exposure assessment, and risk characterization. Hazard identification of chemicals refers to the initial stage of risk assessment and generally utilizes animal studies to evaluate toxicity. Since 2013, however, toxicity studies of cosmetic ingredients using animals have not been permitted in the EU and alternative toxicity test methods for animal studies have momentum to be developed for cosmetic ingredients. In this paper, we briefly review the alternative test methods that are available for cosmetic ingredients including read-across, in silico, in chemico, and invitro methods. In addition, new technologies such as omics and artificial intelligence (AI) have been discussed to expand or improve the knowledge and hazard identification of cosmetic ingredients. Aggregate exposure of cosmetic ingredients is another safety issue and methods for its improvement were reviewed. There have been concerns over the safety of nano-cosmetics for a long time, but the risk of nano-cosmetics remains unclear. Therefore, current issues of cosmetic risk assessment are discussed and expert opinion will be provided for the safety of cosmetics.

A cross-industry collaboration to assess if acute oral toxicity (Q)SAR models are fit-for-purpose for GHS classification and labelling

This study assesses whether currently available acute oral toxicity (AOT) in silico models, provided by the widely employed Leadscope software, are fit-for-purpose for categorization and labelling of chemicals. As part of this study, a large data set of proprietary and marketed compounds from multiple companies (pharmaceutical, plant protection products, and other chemical industries) was assembled to assess the models' performance. The absolute percentage of correct or more conservative predictions, based on a comparison of experimental and predicted GHS categories, was approximately 95%, after excluding a small percentage of inconclusive (indeterminate or out of domain) predictions. Since the frequency distribution across the experimental categories is skewed towards low toxicity chemicals, a balanced assessment was also performed. Across all compounds which could be assigned to a well-defined experimental category, the average percentage of correct or more conservative predictions was around 80%. These results indicate the potential for reliable and broad application of these models across different industrial sectors. This manuscript describes the evaluation of these models, highlights the importance of an expert review, and provides guidance on the use of AOT models to fulfill testing requirements, GHS classification/labelling, and transportation needs.

Screening for Ames mutagenicity of food flavor chemicals by (quantitative) structure-activity relationship

Background

(Quantitative) Structure-Activity Relationship ((Q)SAR) is a promising approach to predict the potential adverse effects of chemicals based on their structure without performing toxicological studies. We evaluate the mutagenicity of food flavor chemicals by (Q) SAR tools, identify potentially mutagenic chemicals, and verify their mutagenicity by actual Ames test.

Results

The Ames mutagenicity of 3942 food flavor chemicals was predicted using two (Q)SAR) tools, DEREK Nexus and CASE Ultra. Three thousand five hundred seventy-five chemicals (91%) were judged to be negative in both (Q) SAR tools, and 75 chemicals (2%) were predicted to be positive in both (Q) SAR tools. When the Ames test was conducted on ten of these positive chemicals, nine showed positive results.

Conclusion

The (Q) SAR method can be used for screening the mutagenicity of food flavors.

Supplementary Information

The online version contains supplementary material available at 10.1186/s41021-020-00171-1.

Management of pharmaceutical ICH M7 (Q)SAR predictions - The impact of model updates

Pharmaceutical applicants conduct (Q)SAR assessments on identified and theoretical impurities to predict their mutagenic potential. Two complementary models—one rule-based and one statistical-based—are used, followed by expert review. (Q)SAR models are continuously updated to improve predictions, with new versions typically released on a yearly basis. Numerous releases of (Q)SAR models will occur during the typical 6–7 years of drug development until new drug registration. Therefore, it is important to understand the impact of model updates on impurity mutagenicity predictions over time. Compounds representative of pharmaceutical impurities were analyzed with three rule- and three statistical-based models covering a 4–8 year period, with the individual time frame being dependent on when the individual models were initially made available. The largest changes in the combined outcome of two complementary models were from positive or equivocal to negative and from negative to equivocal. Importantly, the cumulative change of negative to positive predictions was small in all models (<5%) and was further reduced when complementary models were combined in a consensus fashion. We conclude that model updates of the type evaluated in this manuscript would not necessarily require re-running a (Q)SAR prediction unless there is a specific need. However, original (Q)SAR predictions should be evaluated when finalizing the commercial route of synthesis for marketing authorization.

Exploring the SAR of the β-Ketoacyl-ACP Synthase Inhibitor GSK3011724A and Optimization around a Genotoxic Metabolite

In the course of optimizing a novel indazole sulfonamide series that inhibits β-ketoacyl-ACP synthase (KasA) of Mycobacterium tuberculosis, a mutagenic aniline metabolite was identified. Further lead optimization efforts were therefore dedicated to eliminating this critical liability by removing the embedded aniline moiety or modifying its steric or electronic environment. While the narrow SAR space against the target ultimately rendered this goal unsuccessful, key structural knowledge around the binding site of this underexplored target for TB was generated to inform future discovery efforts.

Development of improved QSAR models for predicting the outcome of the in vivo micronucleus genetic toxicity assay

All drugs entering clinical trials are expected to undergo a series of in vitro and in vivo genotoxicity tests as outlined in the International Council on Harmonization (ICH) S2 (R1) guidance. Among the standard battery of genotoxicity tests used for pharmaceuticals, the in vivo micronucleus assay, which measures the frequency of micronucleated cells mostly from blood or bone marrow, is recommended for detecting clastogens and aneugens. (Quantitative) structure-activity relationship [(Q)SAR] models may be used as early screening tools by pharmaceutical companies to assess genetic toxicity risk during drug candidate selection. Models can also provide decision support information during regulatory review as part of the weight-of-evidence when experimental data are insufficient. In the present study, two commercial (Q)SAR platforms were used to construct in vivo micronucleus models from a recently enhanced in-house database of non-proprietary study findings in mice. Cross-validated performance statistics for the new models showed sensitivity of up to 74% and negative predictivity of up to 86%. In addition, the models demonstrated cross-validated specificity of up to 77% and coverage of up to 94%. These new models will provide more reliable predictions and offer an investigational approach for drug safety assessment with regards to identifying potentially genotoxic compounds.

Transitioning to composite bacterial mutagenicity models in ICH M7 (Q)SAR analyses

The International Council on Harmonisation (ICH) M7(R1) guideline describes the use of complementary (quantitative) structure-activity relationship ((Q)SAR) models to assess the mutagenic potential of drug impurities in new and generic drugs. Historically, the CASE Ultra and Leadscope software platforms used two different statistical-based models to predict mutations at G-C (guanine-cytosine) and A-T (adenine-thymine) sites, to comprehensively assess bacterial mutagenesis. In the present study, composite bacterial mutagenicity models covering multiple mutation types were developed. These new models contain more than double the number of chemicals (n = 9,254 and n = 13,514) than the corresponding non-composite models and show better toxicophore coverage. Additionally, the use of a single composite bacterial mutagenicity model simplifies impurity analysis in an ICH M7 (Q)SAR workflow by reducing the number of model outputs requiring review. An external validation set of 388 drug impurities representing proprietary pharmaceutical chemical space showed performance statistics ranging from of 66-82% in sensitivity, 91-95% in negative predictivity and 96% in coverage. This effort represents a major enhancement to these (Q)SAR models and their use under ICH M7(R1), leading to improved patient safety through greater predictive accuracy, applicability, and efficiency when assessing the bacterial mutagenic potential of drug impurities.

Genetic toxicology in silico protocol

In silico toxicology (IST) approaches to rapidly assess chemical hazard, and usage of such methods is increasing in all applications but especially for regulatory submissions, such as for assessing chemicals under REACH as well as the ICH M7 guideline for drug impurities. There are a number of obstacles to performing an IST assessment, including uncertainty in how such an assessment and associated expert review should be performed or what is fit for purpose, as well as a lack of confidence that the results will be accepted by colleagues, collaborators and regulatory authorities. To address this, a project to develop a series of IST protocols for different hazard endpoints has been initiated and this paper describes the genetic toxicity in silico (GIST) protocol. The protocol outlines a hazard assessment framework including key effects/mechanisms and their relationships to endpoints such as gene mutation and clastogenicity. IST models and data are reviewed that support the assessment of these effects/mechanisms along with defined approaches for combining the information and evaluating the confidence in the assessment. This protocol has been developed through a consortium of toxicologists, computational scientists, and regulatory scientists across several industries to support the implementation and acceptance of in silico approaches.

Improvement of quantitative structure-activity relationship (QSAR) tools for predicting Ames mutagenicity: outcomes of the Ames/QSAR International Challenge Project

The International Conference on Harmonization (ICH) M7 guideline allows the use of in silico approaches for predicting Ames mutagenicity for the initial assessment of impurities in pharmaceuticals. This is the first international guideline that addresses the use of quantitative structure–activity relationship (QSAR) models in lieu of actual toxicological studies for human health assessment. Therefore, QSAR models for Ames mutagenicity now require higher predictive power for identifying mutagenic chemicals. To increase the predictive power of QSAR models, larger experimental datasets from reliable sources are required. The Division of Genetics and Mutagenesis, National Institute of Health Sciences (DGM/NIHS) of Japan recently established a unique proprietary Ames mutagenicity database containing 12140 new chemicals that have not been previously used for developing QSAR models. The DGM/NIHS provided this Ames database to QSAR vendors to validate and improve their QSAR tools. The Ames/QSAR International Challenge Project was initiated in 2014 with 12 QSAR vendors testing 17 QSAR tools against these compounds in three phases. We now present the final results. All tools were considerably improved by participation in this project. Most tools achieved >50% sensitivity (positive prediction among all Ames positives) and predictive power (accuracy) was as high as 80%, almost equivalent to the inter-laboratory reproducibility of Ames tests. To further increase the predictive power of QSAR tools, accumulation of additional Ames test data is required as well as re-evaluation of some previous Ames test results. Indeed, some Ames-positive or Ames-negative chemicals may have previously been incorrectly classified because of methodological weakness, resulting in false-positive or false-negative predictions by QSAR tools. These incorrect data hamper prediction and are a source of noise in the development of QSAR models. It is thus essential to establish a large benchmark database consisting only of well-validated Ames test results to build more accurate QSAR models.

A practice of expert review by read-across using QSAR Toolbox

The International Council for Harmonisation of Technical Requirement for Pharmaceuticals for Human Use (ICH) M7 guideline on 'Assessment and Control of DNA Reactive (Mutagenic) Impurities in Pharmaceuticals to Limit Potential Carcinogenic Risk' provides the application of two types of quantitative structure-activity relationship (QSAR) systems (rule- and statistics-based) as an alternative to the Ames test for evaluating the mutagenicity of impurities in pharmaceuticals. M7 guideline also states that the expert reviews can be applied when the outcomes of the two QSAR analyses show any conflicting or inconclusive prediction. However, the guideline does not provide any information of how to conduct expert reviews. Therefore, a conservative approach was chosen in this study, which is based on the intention to capture any mutagenic chemical substances. The 36 chemical substances, which are the model chemical substances in which positive mutagenicity was not observed according to the two types of QSAR analyses (i.e. the results are either conflicting or both negative), were selected from the list of chemical substances with strong mutagenicity known as the reported chemicals under the Industrial Safety and Health Act in Japan. The QSAR Toolbox was used in this study to rationally determine the positive mutagenicity of the 36 model chemical substances by applying a read-across method, a technique to evaluate the endpoint of the model chemical substances using the endpoint information of chemicals that are structurally similar to the model chemical substances. Resulting from the expert review by the read-across method, the 23 model chemical substances (63.8%) were rationally concluded as positive. In addition, 9 out of 11 model chemical substances that were assessed as negative for mutagenicity by both of the QSAR systems had positive analogues, supporting their mutagenicity. These results suggested that the read-across is a useful method, when conducting a conservative approach intended to capture any mutagenic chemical substances.

Extending (Q)SARs to incorporate proprietary knowledge for regulatory purposes: is aromatic N-oxide a structural alert for predicting DNA-reactive mutagenicity?

(Quantitative) structure-activity relationship or (Q)SAR predictions of DNA-reactive mutagenicity are important to support both the design of new chemicals and the assessment of impurities, degradants, metabolites, extractables and leachables, as well as existing chemicals. Aromatic N-oxides represent a class of compounds that are often considered alerting for mutagenicity yet the scientific rationale of this structural alert is not clear and has been questioned. Because aromatic N-oxide-containing compounds may be encountered as impurities, degradants and metabolites, it is important to accurately predict mutagenicity of this chemical class. This article analysed a series of publicly available aromatic N-oxide data in search of supporting information. The article also used a previously developed structure-activity relationship (SAR) fingerprint methodology where a series of aromatic N-oxide substructures was generated and matched against public and proprietary databases, including pharmaceutical data. An assessment of the number of mutagenic and non-mutagenic compounds matching each substructure across all sources was used to understand whether the general class or any specific subclasses appear to lead to mutagenicity. This analysis resulted in a downgrade of the general aromatic N-oxide alert. However, it was determined there were enough public and proprietary data to assign the quindioxin and related chemicals as well as benzo[c][1,2,5]oxadiazole 1-oxide subclasses as alerts. The overall results of this analysis were incorporated into Leadscope's expert-rule-based model to enhance its predictive accuracy.

Mechanistic Reactivity Descriptors for the Prediction of Ames Mutagenicity of Primary Aromatic Amines

Pharmaceutical products are often synthesized by the use of reactive starting materials and intermediates. These can, either as impurities or through metabolic activation, bind to the DNA. Primary aromatic amines belong to the critical classes that are considered potentially mutagenic in the Ames test, so there is a great need for good prediction models for risk assessment. How primary aromatic amines exert their mutagenic potential can be rationalized by the widely accepted nitrenium ion hypothesis of covalent binding to the DNA of reactive electrophiles formed out of the aromatic amines. Since the reactive chemical species is different in chemical structure from the actual compound, it is difficult to achieve good predictions via classical descriptor or fingerprint-based machine learning. In this approach, we use a combination of different molecular and atomic descriptors that is able to describe different mechanistic aspects of the metabolic transformation leading from the primary aromatic amine to the reactive metabolite that binds to the DNA. Applied to a test set, the combination shows significantly better performance than models that only use one of these descriptors and complemented the general internal Ames mutagenicity prediction model at Bayer.

Construction and application of (Q)SAR models to predict chemical-induced in vitro chromosome aberrations

In drug development, genetic toxicology studies are conducted using in vitro and in vivo assays to identify potential mutagenic and clastogenic effects, as outlined in the International Council for Harmonisation (ICH) S2 regulatory guideline. (Quantitative) structure-activity relationship ((Q)SAR) models that predict assay outcomes can be used as an early screen to prioritize pharmaceutical candidates, or later during product development to evaluate safety when experimental data are unavailable or inconclusive. In the current study, two commercial QSAR platforms were used to build models for in vitro chromosomal aberrations in Chinese hamster lung (CHL) and Chinese hamster ovary (CHO) cells. Cross-validated CHL model predictive performance showed sensitivity of 80 and 82%, and negative predictivity of 75 and 76% based on 875 training set compounds. For CHO, sensitivity of 61 and 67% and negative predictivity of 68 and 74% was achieved based on 817 training set compounds. The predictive performance of structural alerts in a commercial expert rule-based SAR software was also investigated and showed positive predictivity of 48-100% for selected alerts. Case studies examining incorrectly-predicted compounds, non-DNA-reactive clastogens, and recently-approved pharmaceuticals are presented, exploring how an investigational approach using similarity searching and expert knowledge can improve upon individual (Q)SAR predictions of the clastogenicity of drugs.

Meeting report, ICH M7 relevant workshop: use of (Q)SAR systems and expert judgment

Use of Quantitative Structure-Activity Relationships ((Q)SAR) prediction tools has been increasing since the International Council for Harmonization of Technical Requirements for Pharmaceuticals for Human Use (ICH) M7 guideline was issued in June 2014. The Japanese Environmental Mutagen Society and the Bacterial Mutagenicity Study Group took the initiative of the workshop on (Q)SAR in 2016 to discuss using (Q)SAR to predict mutagenicity. The aim of the workshop was to form a common understanding on the current use of (Q)SAR tools in industry and for regulatory purposes and on the process of expert judgment. This report summarizes the general session that reviewed the use of (Q)SAR tools and the case study session that discussed expert judgment.

Predicting Aromatic Amine Mutagenicity with Confidence: A Case Study Using Conformal Prediction

The occurrence of mutagenicity in primary aromatic amines has been investigated using conformal prediction. The results of the investigation show that it is possible to develop mathematically proven valid models using conformal prediction and that the existence of uncertain classes of prediction, such as both (both classes assigned to a compound) and empty (no class assigned to a compound), provides the user with additional information on how to use, further develop, and possibly improve future models. The study also indicates that the use of different sets of fingerprints results in models, for which the ability to discriminate varies with respect to the set level of acceptable errors.

Semi-correlations combined with the index of ideality of correlation: a tool to build up model of mutagenic potential

Mutagenicity is the ability of a substance to induce mutations. This hazardous ability of a substance is decisive from point of view of ecotoxicology. The number of substances, which are used for practical needs, grows every year. Consequently, methods for at least preliminary estimation of mutagenic potential of new substances are necessary. Semi-correlations are a special case of traditional correlations. These correlations can be named as "correlations along two parallel lines." This kind of correlation has been tested as a tool to predict selected endpoints, which are represented by only two values: "inactive/active" (0/1). Here this approach is used to build up predictive models for mutagenicity of large dataset (n = 3979). The so-called index of ideality of correlation (IIC) has been tested as a statistical criterion to estimate the semi-correlation. Three random splits of experimental data into the training, invisible-training, calibration, and validation sets were analyzed. Two models were built up for each split: the first model based on optimization without the IIC and the second model based on optimization where IIC is involved in the Monte Carlo optimization. The statistical characteristics of the best model (calculated with taking into account the IIC) n = 969; sensitivity = 0.8050; specificity = 0.9069; accuracy = 0.8648; Matthews's correlation coefficient = 0.7196 (using IIC). Thus, the use of IIC improves the statistical quality of the binary classification models of mutagenic potentials (Ames test) of organic compounds.

Computational identification of structural factors affecting the mutagenic potential of aromatic amines: study design and experimental validation

A grid-based, alignment-independent 3D-SDAR (three-dimensional spectral data-activity relationship) approach based on simulated ¹³C and ¹⁵N NMR chemical shifts augmented with through-space interatomic distances was used to model the mutagenicity of 554 primary and 419 secondary aromatic amines. A robust modeling strategy supported by extensive validation including randomized training/hold-out test set pairs, validation sets, "blind" external test sets as well as experimental validation was applied to avoid over-parameterization and build Organization for Economic Cooperation and Development (OECD 2004) compliant models. Based on an experimental validation set of 23 chemicals tested in a two-strain Salmonella typhimurium Ames assay, 3D-SDAR was able to achieve performance comparable to 5-strain (Ames) predictions by Lhasa Limited's Derek and Sarah Nexus for the same set. Furthermore, mapping of the most frequently occurring bins on the primary and secondary aromatic amine structures allowed the identification of molecular features that were associated either positively or negatively with mutagenicity. Prominent structural features found to enhance the mutagenic potential included: nitrobenzene moieties, conjugated π-systems, nitrothiophene groups, and aromatic hydroxylamine moieties. 3D-SDAR was also able to capture "true" negative contributions that are particularly difficult to detect through alternative methods. These include sulphonamide, acetamide, and other functional groups, which not only lack contributions to the overall mutagenic potential, but are known to actively lower it, if present in the chemical structures of what otherwise would be potential mutagens.

Computational Toxicology and Drug Discovery

The use of computational toxicology methods within drug discovery began in the early 2000s with applications such as predicting bacterial mutagenicity and hERG inhibition. The field has been continuously expanding ever since and the tasks at hand have become more complex. These approaches are now strategically integrated into the risk assessment process, as a complement to in vitro and in vivo methods. Today, computational toxicology can be used in every phase of drug discovery and development, from profiling large libraries early on, to predicting off-target effects in the mid-discovery phase, to assessing potential mutagenic impurities in development and degradants as part of life-cycle management. This chapter provides an overview of the field and describes the application of computational toxicology throughout the entire discovery and development process.

An investigation of the mutagenic activity of salamide - a major impurity of hydrochlorothiazide

Hydrochlorothiazide is a widely used antihypertensive agent and one of its major impurities, salamide (4-amino-6-chlorobenzene-1,3-disulphonamide), has a chemical structure containing a primary amino group, a functional group that has previously been reported to be associated with carcinogenic activity. It is known that hydrochlorothiazide purity is a challenging problem for the pharmaceutical industry. As there were no prior mutagenicity data for the impurity salamide, the aim was to investigate its mutagenicity in this study. Salamide was tested for mutagenic potential in Salmonella typhimurium TA98, TA100, TA 1535, TA 1537, and E. coli WP2 uvrA + E. coli WP2 [pKM101] strains at six different concentrations, the highest concentration being the 5000 μg/plate. In both the presence and absence of the metabolic activation system, no mutagenic activity was observed. Results indicated that salamide should be classified as an ordinary impurity and controlled according to Q3A(R2) and Q3B(R2) guidelines.