Non-Linear QSAR Treatment of Genotoxicity

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

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

Modelling of antiproliferative activity measured in HeLa cervical cancer cells in a series of xanthene derivatives

Cancer remains one of the leading causes of death in humans, and new drug substances are therefore being developed. Thus, the anti-cancer activity of xanthene derivatives has become an important topic in the development of new and potent anti-cancer drug substances. Previously published novel series of xanthen-3-one and xanthen-1,8-dione derivatives have been synthesized in one of our laboratories and showed anti-proliferative activity in HeLa cancer cell lines. This series serves as a good basis to develop quantitative structure-activity relationship (QSAR), to study the relations between anti-proliferative activity and chemical structures. A QSAR model has been derived that relies only on two-dimensional molecular descriptors, providing mechanistic insight into the anti-proliferative activity of xanthene derivatives. The model is validated internally and externally and additionally with the set of inactive compounds of the original data, confirming model applicability for the design and discovery of novel xanthene derivatives. The QSAR model is available at the QsarDB repository (http://dx.doi.10.15152/QDB.237).

Chemical structure and correlation analysis of HIV-1 NNRT and NRT inhibitors and database-curated, published inhibition constants with chemical structure in diverse datasets

Human immunodeficiency virus (HIV-1) reverse transcriptase is a major target for designing anti-HIV drugs. Developed inhibitors are divided into non-nucleoside analog reverse-transcriptase inhibitors (NNRTIs) and nucleoside analog reverse-transcriptase inhibitors (NRTIs) depending on their mechanism. Given that many inhibitors have been studied and for many of them binding affinity constants have been calculated, it is beneficial to analyze the chemical landscape of these families of inhibitors and correlate these inhibition constants with molecular structure descriptors. For this, the HIV-1 RT data was retrieved from the ChEMBL database, carefully curated, and original literature verified, grouped into NRTIs and NNRTIs, analyzed using a hierarchical scaffold classification method and modelled with best multi-linear regression approach. Analysis of the HIV-1 NNRTIs subset results in ten different common structural parent types of oxazepanone, piperazinone, pyrazine, oxazinanone, diazinanone, pyridine, pyrrole, diazepanone, thiazole, and triazine. The same analysis for HIV-1 NRTIs groups structures into four different parent types of uracil, pyrimide, pyrimidione, and imidazole. Each scaffold tree corresponding to the parent types has been carefully analyzed and examined, and changes in chemical structure favorable to potency and stability are highlighted. For both subsets, descriptive and predictive QSAR models are derived, discussed and externally validated, revealing general trends in relationships between molecular structure and binding affinity constants in structurally diverse datasets. Data and QSAR models are available at the QsarDB repository (http://dx.doi.org/10.15152/QDB.202).

Computational methods in drug discovery

The process for drug discovery and development is challenging, time consuming and expensive. Computer-aided drug discovery (CADD) tools can act as a virtual shortcut, assisting in the expedition of this long process and potentially reducing the cost of research and development. Today CADD has become an effective and indispensable tool in therapeutic development. The human genome project has made available a substantial amount of sequence data that can be used in various drug discovery projects. Additionally, increasing knowledge of biological structures, as well as increasing computer power have made it possible to use computational methods effectively in various phases of the drug discovery and development pipeline. The importance of in silico tools is greater than ever before and has advanced pharmaceutical research. Here we present an overview of computational methods used in different facets of drug discovery and highlight some of the recent successes. In this review, both structure-based and ligand-based drug discovery methods are discussed. Advances in virtual high-throughput screening, protein structure prediction methods, protein-ligand docking, pharmacophore modeling and QSAR techniques are reviewed.

Measurement of baseline toxicity and QSAR analysis of 50 non-polar and 58 polar narcotic chemicals for the alga Pseudokirchneriella subcapitata

In this paper a set of homogenous experimental algal toxicity data was measured for 50 non-polar narcotic chemicals using the alga Pseudokirchneriella subcapitata in a closed test with a growth rate endpoint. Most of the tested compounds are high volume industrial chemicals that so far lacked published REACH-compliant algal growth inhibition values. The test protocol fulfilled the criteria set forth in the OECD guideline 201 and had the same sensitivity as the open test which allowed direct comparison of toxicity values. Baseline QSAR model for non-polar narcotic compounds was established and compared with previous analogous models. Multi-linear QSAR model was derived for the non-polar and 58 previously tested polar (anilines and phenols) narcotic compounds modulating hydrophobicity, molecular size, electronic and molecular stability effects coded in the molecular descriptors. Descriptors in the model were analyzed and applicability domain was assessed providing further guidelines for the in silico prediction purposes in decision support while performing risk assessment. QSAR models in the manuscript are available on-line through QsarDB repository for exploring and prediction services (http://hdl.handle.net/10967/106).

Quantitative relationship between rate constants and molecular structure descriptors for the gas phase hydrogen abstraction reactions

The abstraction of hydrogen by general radicals has a wide role in environmental and also in technological processes because it results in reactive free radicals that play a vital role in atmospheric chemistry and also in biochemical processes. In addition to experimental studies, the theoretical modelling of this elementary reaction has been important for understanding and predicting respective rate constants. In this paper, molecular descriptors in the context of a QSAR approach are used to codify the relationship between molecular structure and rate constants. Unique experimental data is collected from the literature for the reaction R(i)• + R(j)H → R(i)H + R(j)•, where R(i)• = H• and R(j)• are diverse radicals. The four-parameter QSAR model (n = 34, r(2) = 0.81, r(2)(CV) = 0.74, r(2)(scr) = 0.12, s(2) = 0.19) is presented for the bimolecular rate constants, accompanied with model diagnostics and analysis of descriptors in the model.

Comparative analysis of local and consensus quantitative structure-activity relationship approaches for the prediction of bioconcentration factor

Quantitative structure-activity relationships (QSARs) are broadly classified as global or local, depending on their molecular constitution. Global models use large and diverse training sets covering a wide range of chemical space. Local models focus on smaller structurally or chemically similar subsets that are conventionally selected by human experts or alternatively using clustering analysis. The current study focuses on the comparative analysis of different clustering algorithms (expectation-maximization, K-means and hierarchical) for seven different descriptor sets as structural characteristics and two rule-based approaches to select subsets for designing local QSAR models. A total of 111 local QSAR models are developed for predicting bioconcentration factor. Predictions from local models were compared with corresponding predictions from the global model. The comparison of coefficients of determination (r(2)) and standard deviations for local models with similar subsets from the global model show improved prediction quality in 97% of cases. The descriptor content of derived QSARs is discussed and analyzed. Local QSAR models were further consolidated within the framework of consensus approach. All different consensus approaches increased performance over the global and local models. The consensus approach reduced the number of strongly deviating predictions by evening out prediction errors, which were produced by some local QSARs.

Biopartitioning micellar chromatography to predict mutagenicity of aromatic amines

Mutagenicity is a toxicity endpoint associated with the chronic exposure to chemicals. Aromatic amines have considerable industrial and environmental importance due to their widespread use in industry and their mutagenic capacity. Biopartitioning micellar chromatography (BMC), a mode of micellar liquid chromatography that uses micellar mobile phases of Brij35 in adequate experimental conditions, has demonstrated to be useful in mimicking the drug partitioning process into biological systems. In this paper, the usefulness of BMC for predicting mutagenicity of aromatic amines is demonstrated. A multiple linear regression (MLR) model based on BMC retention data is proposed and compared with other ones reported in bibliography. The proposed model present better or similar descriptive and predictive capability.

Neural Networks Convergence Using Physicochemical Data

An investigation of the neural network convergence and prediction based on three optimization algorithms, namely, Levenberg-Marquardt, conjugate gradient, and delta rule, is described. Several simulated neural networks built using the above three algorithms indicated that the Levenberg-Marquardt optimizer implemented as a back-propagation neural network converged faster than the other two algorithms and provides in most of the cases better prediction. These conclusions are based on eight physicochemical data sets, each with a significant number of compounds comparable to that usually used in the QSAR/QSPR modeling. The superiority of the Levenberg-Marquardt algorithm is revealed in terms of functional dependence of the change of the neural network weights with respect to the gradient of the error propagation as well as distribution of the weight values. The prediction of the models is assessed by the error of the validation sets not used in the training process.

Reparameterized Austin Model 1 for quantitative structure–property relationships in liquid media

A reparameterization of the quantum-chemical AM1 (Austin Model 1) model has been carried out using a nonlinear optimization based on a modification of the Levenberg-Marquardt technique. The optimum numerical values for the one-electron resonance integral parameters (beta (s) and beta (p)) and core-core repulsion atomic parameters alpha were obtained for the elements H, C, N, O, Cl and Br using the statistical fit of a two-parameter QSPR equation for the boiling points of organic compounds. A substantially improved two-parameter correlation (R2=0.9685, s=13.48 K) was obtained by using the new optimized parameters. The QSPR equation employs two molecular descriptors, a bulk cohesiveness descriptor, [Formula: see text] and the area-weighted surface charge of hydrogen-bonding donor atom(s) in the molecule. The model developed shows remarkably accurate predictions of the normal boiling points for nine additional simple inorganic compounds. The new parameters were tested on the critical temperatures of 165 organic compounds. The new QSPR model obtained for this property was found to be statistically significantly better than the original model. [Figure: see text].

Structure-Mutagenicity Study of 12 Trimethylimidazopyridine Isomers Using Orbital Energies and “Spectrum-like Representation” As Descriptors

The set of 12 trimethylimidazopyridine isomers with mutagenic potency toward two strains of Salmonella was treated in this study. Ten isomers with known mutagenic properties were taken to build the models. Fifteen molecular orbital energies, or a "spectrum-like" representation of 3D structures, were taken as descriptors. As modeling techniques the multiple linear regression and the counter propagation neural network were applied. Models were tested with the recall ability test and the leave-one-out cross-validation tests. For two isomers, which have not been synthesized yet, we report predicted values for both mutagenic potencies obtained with different models. The best models were found when unoccupied molecular orbital energies are among the descriptors.

Predicting the genotoxicity of secondary and aromatic amines using data subsetting to generate a model ensemble

Binary quantitative structure-activity relationship (QSAR) models are developed to classify a data set of 334 aromatic and secondary amine compounds as genotoxic or nongenotoxic based on information calculated solely from chemical structure. Genotoxic endpoints for each compound were determined using the SOS Chromotest in both the presence and absence of an S9 rat liver homogenate. Compounds were considered genotoxic if assay results indicated a positive genotoxicity hit for either the S9 inactivated or S9 activated assay. Each compound in the data set was encoded through the calculation of numerical descriptors that describe various aspects of chemical structure (e.g. topological, geometric, electronic, polar surface area). Furthermore, five additional descriptors that focused on the secondary and aromatic nitrogen atoms in each molecule were calculated specifically for this study. Descriptor subsets were examined using a genetic algorithm search engine interfaced with a k-Nearest Neighbor fitness evaluator to find the most information-rich subsets, which ultimately served as the final predictive models. Models were chosen for their ability to minimize the total number of misclassifications, with special attention given to those models that possessed fewer occurrences of positive toxicity hits being misclassified as nontoxic (false negatives). In addition, a subsetting procedure was used to form an ensemble of models using different combinations of compounds in the training and prediction sets. This was done to ensure that consistent results could be obtained regardless of training set composition. The procedure also allowed for each compound to be externally validated three times by different training set data with the resultant predictions being used in a "majority rules" voting scheme to produce a consensus prediction for each member of the data set. The individual models produced an average training set classification rate of 71.6% and an average prediction set classification rate of 67.7%. However, the model ensemble was able to correctly classify the genotoxicity of 72.2% of all prediction set compounds.

The computational prediction of toxicity

Recent developments in the prediction of toxicity from chemical structure have been reviewed. Attention has been drawn to some of the problems that can be encountered in the area of predictive toxicology, including the need for a multi-disciplinary approach and the need to address mechanisms of action. Progress has been hampered by the sparseness of good quality toxicological data. Perhaps too much effort has been devoted to exploring new statistical methods rather than to the creation of data sets for hitherto uninvestigated toxicological endpoints and/or classes of chemicals.