Photomutagenicity of N-nitrosoproline dissolved in non-aqueous solvent, oleic acid.

In the present study, we investigated the genotoxicity of the active products formed from N-nitrosoproline (NPRO) dissolved in oleic acid following ultraviolet A (UVA) irradiation, bypassing the need for metabolic activation. We previously demonstrated the photomutagenicity of NPRO dissolved in a phosphate-buffered solution. It has been suggested that the association of the nitrosamine group with acid ions facilitates rapid photodissociation and photoactivation. We hypothesized that NPRO's inherent carboxyl group may mimic an acid, inducing photodissociation and photomutagenicity, even in a non-aqueous solvent lacking acidic ions. Following UVA irradiation, NPRO dissolved in oleic acid exhibited a dose-dependent mutagenic activity. Similar results were obtained when NPRO was dissolved in linoleic acid and triolein. Nitric oxide formation, which is dependent on NPRO concentration, is accompanied by mutagenic activity. The mutagenicity spectrum obtained in response to NPRO irradiation followed the absorption curve of NPRO dissolved in oleic acid. Irradiated NPRO in oleic acid displayed relative stability, retaining approximately 18, 36, and 63 % of initial mutagenicity after 10 days of storage at 25, 4, and -20 °C, respectively. Thus NPRO stored in a fatty environment undergoes photoactivation upon irradiation, leading to genotoxicity.

Investigation into the genotoxic impurity, 1-methyl-4-nitrosopiperazine, in rifampicin.

This study assessed the presence of the genotoxic impurity 1-methyl-4-nitrosopiperazine (MNP) in 27 batches of rifampicin capsules obtained from 11 manufacturers in China. While they were below the temporary limit of 5 ppm set by the US Food and Drug Administration, the observed levels (0.33-2.36 ppm) exceeded the acceptable threshold of 0.16 ppm. Building upon preliminary findings and degradation experiments, we concluded that MNP is a by-product of the oxidative degradation of rifampicin or is introduced via oxidation or nitrosation during the synthesis process involving 1-methyl-4-aminopiperazine. The pathways of MNP formation were confirmed in this study. Furthermore, we observed that the addition of antioxidants, sealed storage, and selection of dominant crystal forms can aid in controlling MNP levels.

Design and synthesis of large Stokes shift DNA dyes with reduced genotoxicity.

Despite intensive search over the past decades, only a few small-molecule DNA fluorescent dyes were found with large Stokes shifts. These molecules, however, are often too toxic for widespread usage. Here, we designed DNA-specific fluorescent dyes rooted in benzimidazole architectures with a hitherto unexplored molecular framework based on thiazole-benzimidazole scaffolding. We further incorporated a pyrazole ring with an extended sidechain to prevent cell penetration. These novel benzimidazole derivatives were predicted by quantum calculations and subsequently validated to have large Stokes shifts ranging from 135 to 143 nm, with their emission colors changed from capri blue for the Hoechst reference compound to iguana green. These readily-synthesized compounds, which displayed improved DNA staining intensity and detection limits along with a complete loss of capability for cellular membrane permeation and negligible mutagenic effects as designed, offer a safer alternative to the existing high-performance small-molecule DNA fluorescent dyes.

In Vitro Metabolism and p53 Activation of Genotoxic Chemicals: Abiotic CYP Enzyme vs Liver Microsomes.

Chemicals often require metabolic activation to become genotoxic. Established test guidelines recommend the use of the rat liver S9 fraction or microsomes to introduce metabolic competence to in vitro cell-based bioassays, but the use of animal-derived components in cell culture raises ethical concerns and may lead to quality issues and reproducibility problems. The aim of the present study was to compare the metabolic activation of cyclophosphamide (CPA) and benzo[a]pyrene (BaP) by induced rat liver microsomes and an abiotic cytochrome P450 (CYP) enzyme based on a biomimetic porphyrine catalyst. For the detection of genotoxic effects, the chemicals were tested in a reporter gene assay targeting the activation of the cellular tumor protein p53. Both chemicals were metabolized by the abiotic CYP enzyme and the microsomes. CPA showed no activation of p53 and low cytotoxicity without metabolic activation, but strong activation of p53 and increased cytotoxicity upon incubation with liver microsomes or abiotic CYP enzyme. The effect concentration causing a 1.5-fold induction of p53 activation was very similar with both metabolization systems (within a factor of 1.5), indicating that genotoxic metabolites were formed at comparable concentrations. BaP also showed low cytotoxicity and no p53 activation without metabolic activation. The activation of p53 was detected for BaP upon incubation with active and inactive microsomes at similar concentrations, indicating experimental artifacts caused by the microsomes or NADPH. The activation of BaP with the abiotic CYP enzyme increased the cytotoxicity of BaP by a factor of 8, but no activation of p53 was detected. The results indicate that abiotic CYP enzymes may present an alternative to rat liver S9 fraction or microsomes for the metabolic activation of test chemicals, which are completely free of animal-derived components. However, an amendment of existing test guidelines would require testing of more chemicals and genotoxicity end points.

Preparation and optimization of dummy molecularly imprinted polymer-based solid-phase extraction system for selective enrichment of p-toluene sulfonate esters genotoxic impurities.

Sulfonate esters, one class of genotoxic impurities (GTIs), have gained significant attention in recent years due to their potential to cause genetic mutations and cancer. In the current study, we employed the dummy template molecular imprinting technology with a dummy template molecule replacing the target molecule to establish a pretreatment method for samples containing p-toluene sulfonate esters. Through computer simulation and ultraviolet-visible spectroscopy analysis, the optimal functional monomer acrylamide and polymerization solvent chloroform were selected. Subsequently, a dummy template molecularly imprinted polymer (DMIP) was prepared by the precipitation polymerization method, and the polymer was characterized in morphology, particle size, and composition. The results of the adsorption and enrichment study demonstrated that the DMIP has high adsorption capability (Q = 7.88 mg/g) and favorable imprinting effects (IF = 1.37); Further, it could simultaneously adsorb three p-toluene sulfonate esters. The optimal adsorption conditions were obtained by conditional optimization of solid-phase extraction (SPE). A pH 7 solution was selected as the loading condition, the methanol/1 % phosphoric acid solution (20:80, v/v) was selected as the washing solution, and acetonitrile containing 10 % acetic acid in 6 mL was selected as the elution solvent. Finally, we determined methyl p-toluene sulfonate alkyl esters, ethyl p-toluene sulfonate alkyl esters, and isopropyl p-toluene sulfonate alkyl esters in tosufloxacin toluene sulfonate and capecitabine at the 10 ppm level (relative to 1 mg/mL active pharmaceutical ingredient (API) samples) by using DMIP-based SPE coupled with HPLC. This approach facilitated the selective enrichment of p-toluene sulfonate esters GTIs from complex API samples.

Metabolism-based category formation for the prioritisation of genotoxicity hazard assessment for plant protection product residues (Part 4): α-Chloroacetamides.

In dietary risk assessment of plant protection products, residues of active ingredients and their metabolites need to be evaluated for their genotoxic potential. The European Food Safety Authority recommend a tiered approach focussing assessment and testing on classes of similar chemicals. To characterise similarity, in terms of metabolism, a metabolic similarity profiling scheme has been developed from an analysis of 69 α-chloroacetamide herbicides for which either Ames, chromosomal aberration or micronucleus test results are publicly available. A set of structural space alerts were defined, each linked to a key metabolic transformation present in the α-chloroacetamide metabolic space. The structural space alerts were combined with covalent chemistry profiling to develop categories suitable for chemical prioritisation via read-across. The method is a robust and reproducible approach to such read-across predictions, with the potential to reduce unnecessary testing. The key challenge in the approach was identified as being the need for metabolism data individual groups of plant protection products as the basis for the development of the structural space alerts.

Determination of genotoxic impurities of aromatic aldehydes in pharmaceutical preparations by high performance liquid chromatography after derivatization with N-Cyclohexyl-4-hydrazino-1,8-naphthalenediimide.

The detection of aromatic aldehydes, considered potential genotoxic impurities, holds significant importance during drug development and production. Current analytical methods necessitate complex pre-treatment processes and exhibit insufficient specificity and sensitivity. This study presents the utilization of naphthalenediimide as a pre-column derivatisation reagent to detect aromatic aldehyde impurities in pharmaceuticals via high-performance liquid chromatography (HPLC). We screened a series of derivatisation reagents through density functional theory (DFT) and investigated the phenomenon of photoinduced electron transfer (PET) for both the derivatisation reagents and the resulting products. Optimal experimental conditions for derivatisation were achieved at 40 °C for 60 min. This approach has been successfully applied to detect residual aromatic aldehyde genotoxic impurities in various pharmaceutical preparations, including 4-Nitrobenzaldehyde, 2-Nitrobenzaldehyde, 1,4-Benzodioxane-6-aldehyde, and 5-Hydroxymethylfurfural. The pre-column derivatisation method significantly enhanced detection sensitivity and reduced the limit of detection (LOD), which ranged from 0.002 to 0.008 µg/ml for the analytes, with relative standard deviations < 3 %. The correlation coefficient (R2) >0.998 demonstrated high quality. In chloramphenicol eye drops, the concentration of 4-Nitrobenzaldehyde was measured to be 8.6 µg/mL below the specified concentration, with recoveries ranging from 90.0 % to 119.2 %. In comparison to existing methods, our work simplifies the pretreatment process, enhances the sensitivity and specificity of the analysis, and offers comprehensive insights into impurity detection in pharmaceutical preparations.

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.

Prediction of molecular-specific mutagenic alerts and related mechanisms of chemicals by a convolutional neural network (CNN) model based on SMILES split.

Structural alerts (SAs) are essential to identify chemicals for toxicity evaluation and health risk assessment. We constructed a novel SMILES split-based deep learning model (SSDL) that was trained and verified with 5850 chemicals from the ISSSTY database and 384 external test chemicals from published papers. The training accuracy was above 0.90 and the evaluation metrics (precision, recall and F1-score) all reached 0.78 or above on both internal and external test chemicals. In this model, the molecular-specific fragment importance of chemicals was first quantified independently. Then, the SA identification method based on the importance of these fragments was statistically analyzed and verified with the ISSSTY test and external test chemicals containing one of 28 typical SAs, and most of the performances were better than that of expert rules. Furthermore, a mutagenicity mechanism prediction method was developed using 237 chemicals with four known mutagenic mechanisms based on molecular similarity calibrated by the SSDL method and fragment importance, which significantly improved accuracy in three mechanisms and had comparable accuracy in the other one compared to traditional methods. Overall, the SSDL model quantifying fragment toxicity within molecules would be a novel potentially powerful tool in the determination and visualization of molecular-specific SAs and the prediction of mutagenicity mechanisms for environmental or industrial compounds and drugs.

Ames test study designs for nitrosamine mutagenicity testing: qualitative and quantitative analysis of key assay parameters.

The robust control of genotoxic N-nitrosamine (NA) impurities is an important safety consideration for the pharmaceutical industry, especially considering recent drug product withdrawals. NAs belong to the 'cohort of concern' list of genotoxic impurities (ICH M7) because of the mutagenic and carcinogenic potency of this chemical class. In addition, regulatory concerns exist regarding the capacity of the Ames test to predict the carcinogenic potential of NAs because of historically discordant results. The reasons postulated to explain these discordant data generally point to aspects of Ames test study design. These include vehicle solvent choice, liver S9 species, bacterial strain, compound concentration, and use of pre-incubation versus plate incorporation methods. Many of these concerns have their roots in historical data generated prior to the harmonization of Ames test guidelines. Therefore, we investigated various Ames test assay parameters and used qualitative analysis and quantitative benchmark dose modelling to identify which combinations provided the most sensitive conditions in terms of mutagenic potency. Two alkyl-nitrosamines, N-nitrosodimethylamine (NDMA) and N-nitrosodiethylamine (NDEA) were studied. NDMA and NDEA mutagenicity was readily detected in the Ames test and key assay parameters were identified that contributed to assay sensitivity rankings. The pre-incubation method (30-min incubation), appropriate vehicle (water or methanol), and hamster-induced liver S9, alongside Salmonella typhimurium strains TA100 and TA1535 and Escherichia coli strain WP2uvrA(pKM101) provide the most sensitive combination of assay parameters in terms of NDMA and NDEA mutagenic potency in the Ames test. Using these parameters and further quantitative benchmark dose modelling, we show that N-nitrosomethylethylamine (NMEA) is positive in Ames test and therefore should no longer be considered a historically discordant NA. The results presented herein define a sensitive Ames test design that can be deployed for the assessment of NAs to support robust impurity qualifications.

Evaluation of QSAR models for predicting mutagenicity: outcome of the Second Ames/QSAR international challenge project.

Quantitative structure-activity relationship (QSAR) models are powerful in silico tools for predicting the mutagenicity of unstable compounds, impurities and metabolites that are difficult to examine using the Ames test. Ideally, Ames/QSAR models for regulatory use should demonstrate high sensitivity, low false-negative rate and wide coverage of chemical space. To promote superior model development, the Division of Genetics and Mutagenesis, National Institute of Health Sciences, Japan (DGM/NIHS), conducted the Second Ames/QSAR International Challenge Project (2020-2022) as a successor to the First Project (2014-2017), with 21 teams from 11 countries participating. The DGM/NIHS provided a curated training dataset of approximately 12,000 chemicals and a trial dataset of approximately 1,600 chemicals, and each participating team predicted the Ames mutagenicity of each trial chemical using various Ames/QSAR models. The DGM/NIHS then provided the Ames test results for trial chemicals to assist in model improvement. Although overall model performance on the Second Project was not superior to that on the First, models from the eight teams participating in both projects achieved higher sensitivity than models from teams participating in only the Second Project. Thus, these evaluations have facilitated the development of QSAR models.

Inoculum cell count influences separation efficiency and variance in Ames plate incorporation and Ames RAMOS test.

The Ames test is one of the most applied tools in mutagenicity testing of chemicals ever since its introduction by Ames et al. in the 1970s. Its principle is based on histidine auxotrophic bacteria that regain prototrophy through reverse mutations. In the presence of a mutagen, more reverse mutations occur that become visible as increased bacterial growth on medium without histidine. Many miniaturized formats of the Ames test have emerged to enable the testing of environmental water samples, increase experimental throughput, and lower the required amounts of test substances. However, most of these formats still rely on endpoint determinations. In contrast, the recently introduced Ames RAMOS test determines mutagenicity through online monitoring of the oxygen transfer rate. In this study, the oxygen transfer rate of Salmonella typhimurium TA100 during the Ames plate incorporation test was monitored and compared to the Ames RAMOS test to prove its validity further. Furthermore, the Ames RAMOS test in 96-well scale is newly introduced. For both the Ames plate incorporation and the Ames RAMOS test, the influence of the inoculum cell count on the negative control was highlighted: A lower inoculum cell count led to a higher coefficient of variation. However, a lower inoculum cell count also led to a higher separation efficiency in the Ames RAMOS test and, thus, to better detection of a mutagenic substance at lower concentrations.

Mechanistic Task Groupings Enhance Multitask Deep Learning of Strain-Specific Ames Mutagenicity.

The Ames test is a gold standard mutagenicity assay that utilizes various Salmonella typhimurium strains with and without S9 fraction to provide insights into the mechanisms by which a chemical can mutate DNA. Multitask deep learning is an ideal framework for developing QSAR models with multiple end points, such as the Ames test, as the joint training of multiple predictive tasks may synergistically improve the prediction accuracy of each task. This work investigated how toxicology domain knowledge can be used to handcraft task groupings that better guide the training of multitask neural networks compared to a naïve ungrouped multitask neural network developed on a complete set of tasks. Sixteen S. typhimurium ± S9 strain tasks were used to generate groupings based on mutagenic and metabolic mechanisms that were reflected in correlation data analyses. Both grouped and ungrouped multitask neural networks predicted the 16 strain tasks with a higher balanced accuracy compared with single task controls, with grouped multitask neural networks consistently featuring incremental increases in predictivity over the ungrouped approach. We conclude that the main variable driving these performance improvements is the general multitask effect with mechanistic task groupings acting as an enhancement step to further concentrate synergistic training signals united by a common biological mechanism. This approach enables incorporation of toxicology domain knowledge into multitask QSAR model development allowing for more transparent and accurate Ames mutagenicity prediction.

Maximizing use of existing carcinogenicity data to support acceptable intake levels for mutagenic impurities in pharmaceuticals: Learnings from N-nitrosamine case studies.

The unexpected finding of N-nitrosamine (NA) impurities in many pharmaceutical products raised significant challenges for industry and regulators. In addition to well-studied small molecular weight NAs, many of which are potent rodent carcinogens, novel NAs associated with active pharmaceutical ingredients have been found, many of which have limited or no safety data. A tiered approach to establishing Acceptable Intake (AI) limits for NA impurities has been established using chemical-specific data, read-across, or a class-specific TTC limit. There are ∼140 NAs with some rodent carcinogenicity data, but much of it is older and does not meet current guidelines for what constitutes a 'robust' bioassay. Nevertheless, these data are an important source of information to ensure the best science is used for assessing NA impurities and assuring consumer safety while minimizing impact that can lead to drug shortages. We present several strategies to maximize the use of imperfect data including using a lower confidence limit on a rodent TD50, and leveraging data from multiple NAs. Information on the chemical structure known to impact potency can also support development of an AI or potentially conclude that a particular NA does not fall in the cohort of concern for potent carcinogenicity.

Multiple Instance Learning Improves Ames Mutagenicity Prediction for Problematic Molecular Species.

The prediction of Ames mutagenicity continues to be a concern in both regulatory and pharmacological toxicology. Traditional quantitative structure-activity relationship (QSAR) models of mutagenicity make predictions based on molecular descriptors calculated on a chemical data set used in their training. However, it is known that molecules such as aromatic amines can be non-mutagenic themselves but metabolically activated by S9 rodent liver enzyme in Ames tests forming molecules such as iminoquinones or amine substituents that better stabilize mutagenic nitrenium ions in known pathways of mutagenicity. Modern in silico modeling methods can implicitly model these metabolites through consideration of the structural elements relevant to their formation but do not include explicit modeling of these metabolites' potential activity. These metabolites do not have a known individual mutagenicity label and, in their current state, cannot be fitted into a traditional QSAR model. Multiple instance learning (MIL) however can be applied to a group of metabolites and their parent under a single mutagenicity label. Here we trained MIL models on Ames data, first with an aromatic amines data set (n = 457), a class known to require metabolic activation, and subsequently on a larger data set (n = 6505) incorporating multiple molecular species. MIL was shown to be able to predict Ames mutagenicity with performance in line with previously established models (balanced accuracy = 0.778), suggesting its potential utility in Ames prediction applications. Furthermore, the MIL model predicted well on identified hard-to-predict molecule groups relative to the models in which these molecule groups were identified. These results are presumably due to the increased consideration of the metabolic contribution to the mutagenic outcome. Further exploration of MIL as a supplement to existing models could aid in the prediction of chemicals where implicit modeling of metabolites cannot fully grasp their characteristics. This paper demonstrates the potential of an MIL approach to modeling Ames tests with S9 and is particularly relevant to metabolically activated xenobiotic mutagens.

Augmenting Expert Knowledge-Based Toxicity Alerts by Statistically Mined Molecular Fragments.

Structural alerts are molecular substructures assumed to be associated with molecular initiating events in various toxic effects and an integral part of in silico toxicology. However, alerts derived using the knowledge of human experts often suffer from a lack of predictivity, specificity, and satisfactory coverage. In this work, we present a method to build hybrid QSAR models by combining expert knowledge-based alerts and statistically mined molecular fragments. Our objective was to find out if the combination is better than the individual systems. Lasso regularization-based variable selection was applied on combined sets of knowledge-based alerts and molecular fragments, but the variable elimination was only allowed to happen on the molecular fragments. We tested the concept on three toxicity end points, i.e., skin sensitization, acute Daphnia toxicity, and Ames mutagenicity, which covered both classification and regression problems. Results showed the predictive performance of such hybrid models is, indeed, better than the models based solely on expert alerts or statistically mined fragments alone. The method also enables the discovery of activating and mitigating/deactivating features for toxicity alerts and the identification of new alerts, thereby reducing false positive and false negative outcomes commonly associated with generic alerts and alerts with poor coverage, respectively.

Mutagenic impurities in pharmaceuticals: A critical assessment of the cohort of concern with a focus on N-nitrosamines.

The TTC (Threshold of Toxicological Concern; set at 1.5 µg/day for pharmaceuticals) defines an acceptable patient intake for any unstudied chemical posing a negligible risk of carcinogenicity or other toxic effects. A group of high potency mutagenic carcinogens, defined solely by the presence of particular structural alerts, are referred to as the "cohort of concern" (CoC); aflatoxin-like-, N-nitroso-, and alkyl-azoxy compounds are considered to pose a significant carcinogenic risk at intakes below the TTC. Kroes et al. (2004) derived values for the TTC and CoC in the context of food components, employing a non-transparent dataset never placed in the public domain. Using a reconstructed all-carcinogen dataset from relevant publications, it is now clear that there are exceptions for all three CoC structural classes. N-Nitrosamines represent 62% of the N-nitroso class in the reconstructed dataset. Employing a contemporary dataset, 20% are negative in rodent carcinogenicity bioassays with less than 50% of all N-nitrosamines estimated to fall into the highest risk category. It is recommended that CoC nitrosamines are identified by compound-specific data rather than structural alerts. Thus, it should be possible to distinguish CoC from non-CoC N-nitrosamines in the context of mutagenic impurities described in ICH M7 (R1).

Prioritization of mycotoxins based on mutagenicity and carcinogenicity evaluation using combined in silico QSAR methods.

Mycotoxins and their metabolites are a family of compounds that contains a great diversity of both structure and biological properties. Information on their toxicity is spread within several databases and in scientific literature. Due to the number of molecules and their structure diversity, the cost and time required for hazard evaluation of each compound is unrealistic. In that purpose, new approach methodologies (NAMs) can be applied to evaluate such large set of molecules. Among them, quantitative structure-activity relationship (QSAR) in silico models could be useful to predict the mutagenic and carcinogenic properties of mycotoxins. First, a complete list of 904 mycotoxins and metabolites was built. Then, some known mycotoxins were used to determine the best QSAR tools for mutagenicity and carcinogenicity predictions. The best tool was further applied to the whole list of 904 mycotoxins. At the end, 95 mycotoxins were identified as both mutagen and carcinogen and should be prioritized for further evaluation.

From Structural Alerts to Signature Fragment Alerts: A Case Study on Pyrrolizidine Alkaloids.

Even though modeling is considered a valid alternative to mutagenicity testing for substances with known structures, it can be applied for mixtures only if all of the single chemical structures are identified. Within the present work, we investigate a new avenue to exploit computational toxicology for mixtures, such as plant-based food ingredients. Indeed, considering that in the absence of toxicological information, an important early consideration is whether any substance may be genotoxic through the mutagenic mechanism of action, we tried to establish a correspondence between genotoxic structural alerts (SAs) and so-called signature fragment alerts (SFAs). Once this correspondence is established, chromatograms could be screened for chemical features associated with genotoxic alerts. Pyrrolizidine alkaloids (PAs), a large group of natural toxins (several of them known as genotoxic) were used as a case study because their early identification would bring significant benefits. The method was built using 56 PA pure standards, resulting in the characterization of signature fragment alerts. Finally, the approach was verified in real plant-based samples such as herbal tea and alfalfa, where the screening of signature fragment alerts allowed highlighting quickly the presence of genotoxic PAs in plant-based mixtures. Therefore, the SFA analysis can be used for risk prioritization of newly identified PAs and for their identification in mixtures, contributing to the unnecessary use of animal experimentation for genotoxicity testing.

Automated detection of toxicophores and prediction of mutagenicity using PMCSFG algorithm.

Maximum common substructures (MCS) have received a lot of attention in the chemoinformatics community. They are typically used as a similarity measure between molecules, showing high predictive performance when used in classification tasks, while being easily explainable substructures. In the present work, we applied the Pairwise Maximum Common Subgraph Feature Generation (PMCSFG) algorithm to automatically detect toxicophores (structural alerts) and to compute fingerprints based on MCS. We present a comparison between our MCS-based fingerprints and 12 well-known chemical fingerprints when used as features in machine learning models. We provide an experimental evaluation and discuss the usefulness of the different methods on mutagenicity data. The features generated by the MCS method have a state-of-the-art performance when predicting mutagenicity, while they are more interpretable than the traditional chemical fingerprints.