Prediction of aromatic amines mutagenicity from theoretical molecular descriptors

QSARINS-Chem standalone version: A new platform-independent software to profile chemicals for physico-chemical properties, fate, and toxicity

The new software QSARINS-Chem standalone version is a multiplatform tool, freely downloadable, for the in silico profiling of multiple properties and activities of organic chemicals. This software, which is based on the concept of the QSARINS-chem module embedded in the QSARINS software, has been fully redesigned and redeveloped in the Java™ language. In addition to a selection of models included in the old module, the new software predicts biotransformation rates and aquatic toxicities of pharmaceuticals and personal care products in multiple organisms, and offers a suite of tools for the analysis of predictions. Furthermore, a comprehensive and transparent database of molecular structures is provided. The new QSARINS-Chem standalone version is an informative and solid tool, which is useful to support the assessment of the potential hazard and risks related to organic chemicals and is dedicated to users which are interested in the application of QSARs to generate reliable predictions.

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.

QSAR modeling for predicting mutagenic toxicity of diverse chemicals for regulatory purposes

The safety assessment process of chemicals requires information on their mutagenic potential. The experimental determination of mutagenicity of a large number of chemicals is tedious and time and cost intensive, thus compelling for alternative methods. We have established local and global QSAR models for discriminating low and high mutagenic compounds and predicting their mutagenic activity in a quantitative manner in Salmonella typhimurium (TA) bacterial strains (TA98 and TA100). The decision treeboost (DTB)-based classification QSAR models discriminated among two categories with accuracies of >96% and the regression QSAR models precisely predicted the mutagenic activity of diverse chemicals yielding high correlations (R ²) between the experimental and model-predicted values in the respective training (>0.96) and test (>0.94) sets. The test set root mean squared error (RMSE) and mean absolute error (MAE) values emphasized the usefulness of the developed models for predicting new compounds. Relevant structural features of diverse chemicals that were responsible and influence the mutagenic activity were identified. The applicability domains of the developed models were defined. The developed models can be used as tools for screening new chemicals for their mutagenicity assessment for regulatory purpose.

A knowledge-based expert rule system for predicting mutagenicity (Ames test) of aromatic amines and azo compounds

Cancer is one of the main causes of death in Western countries, and a major issue for human health. Prolonged exposure to a number of chemicals was observed to be one of the primary causes of cancer in occupationally exposed persons. Thus, the development of tools for identifying hazardous chemicals and the increase of mechanistic understanding of their toxicity is a major goal for scientific research. We constructed a new knowledge-based expert system accounting the effect of different substituents for the prediction of mutagenicity (Ames test) of aromatic amines, a class of compounds of major concern because of their widespread application in industry. The herein presented model implements a series of user-defined structural rules extracted from a database of 616 primary aromatic amines, with their Ames test outcomes, aimed at identifying mutagenic and non-mutagenic chemicals. The chemical rationale behind such rules is discussed. Besides assessing the model's ability to correctly classify aromatic amines, its predictivity was further evaluated on a second database of 354 azo dyes, another class of chemicals of major concern, whose toxicity has been predicted on the basis of the toxicity of aromatic amines potentially generated from the metabolic reduction of the azo bond. Good performance in classification on both the amine (MCC, Matthews Correlation Coefficient=0.743) and the azo dye (MCC=0.584) datasets confirmed the predictive power of the model, and its suitability for use on a wide range of chemicals. Finally, the model was compared with a series of well-known mutagenicity predicting software. The good performance of our model compared with other mutagenicity models, especially in predicting azo dyes, confirmed the usefulness of this expert system as a reliable support to in vitro mutagenicity assays for screening and prioritization purposes. The model has been fully implemented as a KNIME workflow and is freely available for downstream users.

QSARINS-chem: Insubria datasets and new QSAR/QSPR models for environmental pollutants in QSARINS

A database of environmentally hazardous chemicals, collected and modeled by QSAR by the Insubria group, is included in the updated version of QSARINS, software recently proposed for the development and validation of QSAR models by the genetic algorithm-ordinary least squares method. In this version, a module, named QSARINS-Chem, includes several datasets of chemical structures and their corresponding endpoints (physicochemical properties and biological activities). The chemicals are accessible in different ways (CAS, SMILES, names and so forth) and their three-dimensional structure can be visualized. Some of the QSAR models, previously published by our group, have been redeveloped using the free online software for molecular descriptor calculation, PaDEL-Descriptor. The new models can be easily applied for future predictions on chemicals without experimental data, also verifying the applicability domain to new chemicals. The QSAR model reporting format (QMRF) of these models is also here downloadable. Additional chemometric analyses can be done by principal component analysis and multicriteria decision making for screening and ranking chemicals to prioritize the most dangerous.

Electron-correlation based externally predictive QSARs for mutagenicity of nitrated-PAHs in Salmonella typhimurium TA100

In quantitative modeling, there are two major aspects that decide reliability and real external predictivity of a structure-activity relationship (SAR) based on quantum chemical descriptors. First, the information encoded in employed molecular descriptors, computed through a quantum-mechanical method, should be precisely estimated. The accuracy of the quantum-mechanical method, however, is dependent upon the amount of electron-correlation it incorporates. Second, the real external predictivity of a developed quantitative SAR (QSAR) should be validated employing an external prediction set. In this work, to analyze the role of electron-correlation, QSAR models are developed for a set of 51 ubiquitous pollutants, namely, nitrated monocyclic and polycyclic aromatic hydrocarbons (nitrated-AHs and PAHs) having mutagenic activity in TA100 strain of Salmonella typhimurium. The quality of the models, through state-of-the-art external validation procedures employing an external prediction set, is compared to the best models known in the literature for mutagenicity. The molecular descriptors whose electron-correlation contribution is analyzed include total energy, energy of HOMO and LUMO, and commonly employed electron-density based descriptors such as chemical hardness, chemical softness, absolute electronegativity and electrophilicity index. The electron-correlation based QSARs are also compared with those developed using quantum-mechanical descriptors computed with advanced semi-empirical (SE) methods such as PM6, PM7, RM1, and ab initio methods, namely, the Hartree-Fock (HF) and the density functional theory (DFT). The models, developed using electron-correlation contribution of the quantum-mechanical descriptors, are found to be not only reliable but also satisfactorily predictive when compared to the existing robust models. The robustness of the models based on descriptors computed through advanced SE methods, is also observed to be comparable to those developed with the electron-correlation based descriptors. The work emphasizes that the correlation-energy can serve as a reliable descriptor to explore the origin of biological activities at the level of electron-dynamics.

Comparative QSTR studies for predicting mutagenicity of nitro compounds

Mutagenicity and carcinogenicity are toxicological endpoints which pose a great concern being the major determinants of cancers and tumours. Nitroarenes possess genotoxic properties as they can form various electrophilic intermediates and adducts with biological systems. Different QSTR techniques were employed to develop models for the prediction of mutagenicity of nitroarenes using a diverse set of 197 nitro aromatic and hetero aromatic molecules. The 2D and 3D QSTR methods used for model development gave statistically significant results. The alignment for 3D methods was obtained by maximum common substructures (MCS) approach, by taking the most mutagenic molecule of the dataset as the template. All the QSTR models were developed with the same set of training and test set molecules. The 3D contours and 2D contribution maps along with molecular fingerprints provide useful information about the mutagenic potentials of the molecules. The GFA based model shows thermodynamic and topological descriptors play an important role in characterizing mutagenicity of nitroarenes. Atomic-level thermodynamic descriptor namely AlogP throws light on hydrophobic features and helps to understand the bilinear model. Topological aspects of these classes of compounds were depicted by the fragment fingerprints and Balaban indices obtained from HQSAR and GFA models, respectively. The predictive abilities of 2D and 3D QSTR models may be useful as a vibrant predictive tool to screen out mutagenic nitroarenes and design safer non-mutagenic nitro compounds.

Approaches for externally validated QSAR modelling of Nitrated Polycyclic Aromatic Hydrocarbon mutagenicity

Nitrated Polycyclic Aromatic Hydrocarbons (nitro-PAHs), ubiquitous environmental pollutants, are recognized mutagens and carcinogens. A set of mutagenicity data (TA100) for 48 nitro-PAHs was modeled by the Quantitative Structure-Activity Relationships (QSAR) regression method, and OECD principles for QSAR model validation were applied. The proposed Multiple Linear Regression (MLR) models are based on two topological molecular descriptors. The models were validated for predictivity by both internal and external validation. For the external validation, three different splitting approaches, D-optimal Experimental Design, Self Organizing Maps (SOM) and Random Selection by activity sampling, were applied to the original data set in order to compare these methodologies and to select the best descriptors able to model each prediction set chemicals independently of the splitting method applied. The applicability domain was verified by the leverage approach.

Random Forest Prediction of Mutagenicity from Empirical Physicochemical Descriptors

Fast-to-calculate empirical physicochemical descriptors were investigated for their ability to predict mutagenicity (positive or negative Ames test) from the molecular structure. Fast methods are highly desired for the screening of large libraries of compounds. Global molecular descriptors and MOLMAP descriptors of bond properties were used to train random forests. Error percentages as low as 15% and 16% were achieved for an external test set with 472 compounds and for the training set with 4083 structures, respectively. High sensitivity and specificity were observed. Random forests were able to associate meaningful probabilities to the predictions and to explain the predictions in terms of similarities between query structures and compounds in the training set.

Predicting genotoxicity of aromatic and heteroaromatic amines using electrotopological state indices

A quantitative structure-activity relationship (QSAR) model relating electrotopological state (E-state) indices and mutagenic potency was previously described by Cash [Mutat. Res. 491 (2001) 31-37] using a data set of 95 aromatic amines published by Debnath et al. [Environ. Mol. Mutagen. 19 (1992) 37-52]. Mutagenic potency was expressed as the number of Salmonella typhimurium TA98 revertants per nmol (LogR). Earlier work on the development of QSARs for the prediction of genotoxicity indicated that numerous methods could be effectively employed to model the same aromatic amines data set, namely, Debnath et al.; Maran et al. [Quant. Struct.-Act. Relat. 18 (1999) 3-10]; Basak et al. [J. Chem. Inf. Comput. Sci. 41 (2001) 671-678]; Gramatica et al. [SAR QSAR Environ. Res. 14 (2003) 237-250]. However, results obtained from external validations of those models revealed that the effective predictivity of the QSARs was well below the potential indicated by internal validation statistics (Debnath et al., Gramatica et al.). The purpose of the current research is to externally validate the model published by Cash using a data set of 29 aromatic amines reported by Glende et al. [Mutat. Res. 498 (2001) 19-37; Mutat. Res. 515 (2002) 15-38] and to further explore the potential utility of using E-state sums for the prediction of mutagenic potency of aromatic amines.