Avoidance of the Ames test liability for aryl-amines via computation

Bacterial mutagenicity test data: collection by the task force of the Japan pharmaceutical manufacturers association

Background

Ames test is used worldwide for detecting the bacterial mutagenicity of chemicals. In silico analyses of bacterial mutagenicity have recently gained acceptance by regulatory agencies; however, current in silico models for prediction remain to be improved. The Japan Pharmaceutical Manufacturers Association (JPMA) organized a task force in 2017 in which eight Japanese pharmaceutical companies had participated. The purpose of this task force was to disclose a piece of pharmaceutical companies’ proprietary Ames test data.

Results

Ames test data for 99 chemicals of various chemical classes were collected for disclosure in this study. These chemicals are related to the manufacturing process of pharmaceutical drugs, including reagents, synthetic intermediates, and drug substances. The structure-activity (mutagenicity) relationships are discussed in relation to structural alerts for each chemical class. In addition, in silico analyses of these chemicals were conducted using a knowledge-based model of Derek Nexus (Derek) and a statistics-based model (GT1_BMUT module) of CASE Ultra. To calculate the effectiveness of these models, 89 chemicals for Derek and 54 chemicals for CASE Ultra were selected; major exclusions were the salt form of four chemicals that were tested both in the salt and free forms for both models, and 35 chemicals called “known” positives or negatives for CASE Ultra. For Derek, the sensitivity, specificity, and accuracy were 65% (15/23), 71% (47/66), and 70% (62/89), respectively. The sensitivity, specificity, and accuracy were 50% (6/12), 60% (25/42), and 57% (31/54) for CASE Ultra, respectively. The ratio of overall disagreement between the CASE Ultra “known” positives/negatives and the actual test results was 11% (4/35). In this study, 19 out of 28 mutagens (68%) were detected with TA100 and/or TA98, and 9 out of 28 mutagens (32%) were detected with either TA1535, TA1537, WP2uvrA, or their combination.

Conclusion

The Ames test data presented here will help avoid duplicated Ames testing in some cases, support duplicate testing in other cases, improve in silico models, and enhance our understanding of the mechanisms of mutagenesis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s41021-021-00206-1.

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.

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.

Discovery of a Potent, Selective T-type Calcium Channel Blocker as a Drug Candidate for the Treatment of Generalized Epilepsies

We report here the discovery and pharmacological characterization of N-(1-benzyl-1H-pyrazol-3-yl)-2-phenylacetamide derivatives as potent, selective, brain-penetrating T-type calcium channel blockers. Optimization focused mainly on solubility, brain penetration, and the search for an aminopyrazole metabolite that would be negative in an Ames test. This resulted in the preparation and complete characterization of compound 66b (ACT-709478), which has been selected as a clinical candidate.

The Necessary Nitrogen Atom: A Versatile High-Impact Design Element for Multiparameter Optimization

There is a continued desire in biomedical research to reduce the number and duration of design cycles required to optimize lead compounds into high-quality chemical probes or safe and efficacious drug candidates. The insightful application of impactful molecular design elements is one approach toward achieving this goal. The replacement of a CH group with a N atom in aromatic and heteroaromatic ring systems can have many important effects on molecular and physicochemical properties and intra- and intermolecular interactions that can translate to improved pharmacological profiles. In this Perspective, the "necessary nitrogen atom" is shown to be a versatile high-impact design element for multiparameter optimization, wherein ≥10-, 100-, or 1000-fold improvement in a variety of key pharmacological parameters can be realized.

Mutagenicity of heteroaromatic amines: Computational study on the influence of methyl substituents

Quantum mechanical calculations were performed to elucidate the factors determining the variations in mutagenic activity within groups of isomeric heteroaromatic amines that differ in the position of methyl substituents. Formation energies for noncovalent complexes and covalent DNA adducts were evaluated by means of high level quantum chemical methods. According to the computational results in this work, covalent adduct stability is proposed to influence the relative mutagenicities of structurally related heterocyclic amines. The stability of covalent C8-dG DNA adducts was found to be mainly determined by π-stacking interactions between the fused ring system of the heteroaromatic amines and the flanking nucleobases. Relative mutagenicity of amines of very related structure is proposed to be regulated by both nitrenium ion and covalent adduct stabilities.

Effect-directed analysis supporting monitoring of aquatic environments--An in-depth overview

Aquatic environments are often contaminated with complex mixtures of chemicals that may pose a risk to ecosystems and human health. This contamination cannot be addressed with target analysis alone but tools are required to reduce this complexity and identify those chemicals that might cause adverse effects. Effect-directed analysis (EDA) is designed to meet this challenge and faces increasing interest in water and sediment quality monitoring. Thus, the present paper summarizes current experience with the EDA approach and the tools required, and provides practical advice on their application. The paper highlights the need for proper problem formulation and gives general advice for study design. As the EDA approach is directed by toxicity, basic principles for the selection of bioassays are given as well as a comprehensive compilation of appropriate assays, including their strengths and weaknesses. A specific focus is given to strategies for sampling, extraction and bioassay dosing since they strongly impact prioritization of toxicants in EDA. Reduction of sample complexity mainly relies on fractionation procedures, which are discussed in this paper, including quality assurance and quality control. Automated combinations of fractionation, biotesting and chemical analysis using so-called hyphenated tools can enhance the throughput and might reduce the risk of artifacts in laboratory work. The key to determining the chemical structures causing effects is analytical toxicant identification. The latest approaches, tools, software and databases for target-, suspect and non-target screening as well as unknown identification are discussed together with analytical and toxicological confirmation approaches. A better understanding of optimal use and combination of EDA tools will help to design efficient and successful toxicant identification studies in the context of quality monitoring in multiply stressed environments.

Multicriteria approach to select polyaromatic river mutagen candidates

The identification of unknown compounds remains one of the most challenging tasks to link observed toxic effects in complex environmental mixtures to responsible toxicants in effect-directed analysis (EDA). Here, a workflow is presented based on nontarget liquid chromatography-high resolution mass spectrometry (LC-HRMS) starting with molecular formulas determined in a previous study. A compound database search (ChemSpider) was performed to retrieve candidates for each formula. Subsequently, the number of candidates was reduced by applying MS-, physical-chemical, and chromatography-based selection criteria including HRMS/MS fragmentation and plausibility, ionization efficiency with different ion sources and detection modes, acid/base behavior, octanol/water partitioning, retention time prediction and finally toxic effects (mutagenicity caused by aromatic amines). The workflow strongly decreased the number of possible candidates and resulted in the tentative identification of possible mutagens and the positive identification of the nonmutagen benzyl(diphenyl) phosphine oxide in a mutagenic fraction. The positive identification of mutagens was hampered by a lack of commercially available standards. The workflow is an innovative and promising approach and forms an excellent basis for possible further advancements.

Aza follow-ups to BI 207524, a thumb pocket 1 HCV NS5B polymerase inhibitor. Part 1: Mitigating the genotoxic liability of an aniline metabolite

A series of heterocyclic aza-analogs of BI 207524 (2), a potent HCV NS5B polymerase thumb pocket 1 inhibitor, was investigated with the goal to reduce the liability associated with the release of a genotoxic aniline metabolite in vivo. Analog 4, containing a 2-aminopyridine aniline isostere that is negative in the Ames test was identified, and was found to provide comparable GT1a/1b potency to 2. Although the cross-species PK profile, poor predicted human liver distribution of analog 4 and allometry principles projected high doses to achieve a strong antiviral response in patients, this work has provided a path forward toward the design of novel thumb pocket 1 NS5B polymerase inhibitors with improved safety profiles.

Comparability of mixed IC₅₀ data - a statistical analysis

The biochemical half maximal inhibitory concentration (IC₅₀) is the most commonly used metric for on-target activity in lead optimization. It is used to guide lead optimization, build large-scale chemogenomics analysis, off-target activity and toxicity models based on public data. However, the use of public biochemical IC₅₀ data is problematic, because they are assay specific and comparable only under certain conditions. For large scale analysis it is not feasible to check each data entry manually and it is very tempting to mix all available IC₅₀ values from public database even if assay information is not reported. As previously reported for K_i database analysis, we first analyzed the types of errors, the redundancy and the variability that can be found in ChEMBL IC₅₀ database. For assessing the variability of IC₅₀ data independently measured in two different labs at least ten IC₅₀ data for identical protein-ligand systems against the same target were searched in ChEMBL. As a not sufficient number of cases of this type are available, the variability of IC₅₀ data was assessed by comparing all pairs of independent IC₅₀ measurements on identical protein-ligand systems. The standard deviation of IC₅₀ data is only 25% larger than the standard deviation of K_i data, suggesting that mixing IC₅₀ data from different assays, even not knowing assay conditions details, only adds a moderate amount of noise to the overall data. The standard deviation of public ChEMBL IC₅₀ data, as expected, resulted greater than the standard deviation of in-house intra-laboratory/inter-day IC₅₀ data. Augmenting mixed public IC₅₀ data by public K_i data does not deteriorate the quality of the mixed IC₅₀ data, if the K_i is corrected by an offset. For a broad dataset such as ChEMBL database a K_i- IC₅₀ conversion factor of 2 was found to be the most reasonable.

Flexible use of QSAR models in predictive toxicology: a case study on aromatic amines

Over the last years, predictive toxicology approaches based on Structure-Activity Relationships have emerged as fundamental tools in the regulatory assessments of chemicals, especially in those programs where regulatory constraints and assessment schemes limit the amount of data available from experimental test methods. Both the qualitative (e.g., Structural Alerts) and the quantitative (Quantitative Structure-Activity Relationships, QSAR) approach can play important roles. However, the two approaches are not familiar to the same extent to the regulators that most often use only the qualitative approach, so that the potentiality of the more sophisticated QSAR approach is neglected. In fact, QSAR is a very flexible tool that allows the user to modulate its response according to different goals and requirements. Here, we present a non-naïve approach to the use of a QSAR relative to a dichotomous biological activity (such as mutagen/nonmutagen), and we show how the user can maximize alternatively the reliability of the prediction of negative compounds (i.e., safe chemicals) or that of positive chemicals (i.e., chemicals that pose high hazard). Because of the environmental and industrial importance of the class of aromatic amines, we apply the approach to a previously published QSAR on the Salmonella typhimurium mutagenicity of these chemicals.

An investigation into pharmaceutically relevant mutagenicity data and the influence on Ames predictive potential

Background

In drug discovery, a positive Ames test for bacterial mutation presents a significant hurdle to advancing a drug to clinical trials. In a previous paper, we discussed success in predicting the genotoxicity of reagent-sized aryl-amines (ArNH₂), a structure frequently found in marketed drugs and in drug discovery, using quantum mechanics calculations of the energy required to generate the DNA-reactive nitrenium intermediate (ArNH:+). In this paper we approach the question of what molecular descriptors could improve these predictions and whether external data sets are appropriate for further training.

Results

In trying to extend and improve this model beyond this quantum mechanical reaction energy, we faced considerable difficulty, which was surprising considering the long history and success of QSAR model development for this test. Other quantum mechanics descriptors were compared to this reaction energy including AM1 semi-empirical orbital energies, nitrenium formation with alternative leaving groups, nitrenium charge, and aryl-amine anion formation energy. Nitrenium formation energy, regardless of the starting species, was found to be the most useful single descriptor. External sets used in other QSAR investigations did not present the same difficulty using the same methods and descriptors. When considering all substructures rather than just aryl-amines, we also noted a significantly lower performance for the Novartis set. The performance gap between Novartis and external sets persists across different descriptors and learning methods. The profiles of the Novartis and external data are significantly different both in aryl-amines and considering all substructures. The Novartis and external data sets are easily separated in an unsupervised clustering using chemical fingerprints. The chemical differences are discussed and visualized using Kohonen Self-Organizing Maps trained on chemical fingerprints, mutagenic substructure prevalence, and molecular weight.

Conclusions

Despite extensive work in the area of predicting this particular toxicity, work in designing and publishing more relevant test sets for compounds relevant to drug discovery is still necessary. This work also shows that great care must be taken in using QSAR models to replace experimental evidence. When considering all substructures, a random forest model, which can inherently cover distinct neighborhoods, built on Novartis data and previously reported external data provided a suitable model.