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 13C and 15N 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.

3D-SDAR modeling of hERG potassium channel affinity: A case study in model design and toxicophore identification.

A dataset of 237 human Ether-à-go-go Related Gene (hERG) potassium channel inhibitors (180 of which were used for model building and validation, whereas 57 constituted the "true" external prediction set) collected from 22 literature sources was modeled by 3D-SDAR. To produce reliable and reproducible classification models for hERG blocking, the initial set of 180 chemicals was split into two subsets: a balanced modeling set consisting of 118 compounds and an unbalanced validation set comprised of 62 compounds. A PLS bagging-like algorithm written in Matlab was used to process the data and assign each compound to one of the two (hERG+ or hERG-) activity classes. The best predictive model evaluated on the basis of a fully randomized hold-out test set (comprising 20% of the modeling set) used 4 latent variables and a grid of 6ppm×6ppm×1Å in the C-C region, 6ppm×30ppm×1Å in the C-N region, and 30ppm×30ppm×1Å in the N-N region. An overall accuracy of 0.84 was obtained for both the hold-out test set and the validation set. Further, an external prediction set consisting of 57 drugs and drug derivatives was used to estimate the true predictive power of the reported 3D-SDAR model - a slight reduction of the overall accuracy down to 0.77 was observed. 3D-SDAR map of the most frequently occurring bins and their projection on the standard coordinate space of the chemical structures allowed identification of a three-center toxicophore composed of two aromatic rings and an amino group. A U test along the distance axis of the most frequently occurring 3D-SDAR bins was used to set the distance limits of the toxicophore. This toxicophore was found to be similar to an earlier reported phospholipidosis (PLD) toxicophore.

Rigorous 3-dimensional spectral data activity relationship approach modeling strategy for ToxCast estrogen receptor data classification, validation, and feature extraction.

The estrogenic potential (expressed as a score composite of 18 high throughput screening bioassays) of 1528 compounds from the ToxCast database was modeled by a 3-dimensional spectral data activity relationship approach (3D-SDAR). Due to a lack of 17 O nuclear magnetic resonance (NMR) simulation software, the most informative carbon-carbon 3D-SDAR fingerprints were augmented with indicator variables representing oxygen atoms from carbonyl and carboxamide, ester, sulfonyl, nitro, aliphatic hydroxyl, and phenolic hydroxyl groups. To evaluate the true predictive performance of the authors' model the United States Environmental Protection Agency provided them with a blind test set consisting of 2008 compounds. Of these, 543 had available literature data-their binding affinity served to estimate the external classification accuracy of the developed model: predictive accuracy of 0.62, sensitivity of 0.71, and specificity of 0.53 were obtained. Compared with alternative modeling techniques, the authors' model displayed very little reduction in performance between the modeling and the prediction set. A 3D-SDAR mapping technique allowed identification of structural features essential for estrogenicity: 1) the presence of a phenolic OH group or cyclohexenone, 2) a second aromatic or phenolic ring at a distance of 6 Što 8 Šfrom the oxygen of the first phenol ring, 3) the presence of a methyl group approximately 6 Šaway from the centroid of a phenol ring, and 4) a carbonyl group in close proximity (∼4 Šmeasured to the centroid) to 1 of the phenol rings. Environ Toxicol Chem 2017;36:823-830. Published 2016 Wiley Periodicals Inc. on behalf of SETAC. This article is a US government work and, as such, is in the public domain in the United States of America.

4-Fluoro-3',4',5'-trimethoxychalcone as a new anti-invasive agent. From discovery to initial validation in an in vivo metastasis model.

Invasion and metastasis are responsible for 90% of cancer-related mortality. Herein, we report on our quest for novel, clinically relevant inhibitors of local invasion, based on a broad screen of natural products in a phenotypic assay. Starting from micromolar chalcone hits, a predictive QSAR model for diaryl propenones was developed, and synthetic analogues with a 100-fold increase in potency were obtained. Two nanomolar hits underwent efficacy validation and eADMET profiling; one compound was shown to increase the survival time in an artificial metastasis model in nude mice. Although the molecular mechanism(s) by which these substances mediate efficacy remain(s) unrevealed, we were able to eliminate the major targets commonly associated with antineoplastic chalcones.

Computational identification of a phospholipidosis toxicophore using (13)C and (15)N NMR-distance based fingerprints.

Modified 3D-SDAR fingerprints combining (13)C and (15)N NMR chemical shifts augmented with inter-atomic distances were used to model the potential of chemicals to induce phospholipidosis (PLD). A curated dataset of 328 compounds (some of which were cationic amphiphilic drugs) was used to generate 3D-QSDAR models based on tessellations of the 3D-SDAR space with grids of different density. Composite PLS models averaging the aggregated predictions from 100 fully randomized individual models were generated. On each of the 100 runs, the activities of an external blind test set comprised of 294 proprietary chemicals were predicted and averaged to provide composite estimates of their PLD-inducing potentials (PLD+ if PLD is observed, otherwise PLD-). The best performing 3D-QSDAR model utilized a grid with a density of 8ppm×8ppm in the C-C region, 8ppm×20ppm in the C-N region and 20ppm×20ppm in the N-N region. The classification predictive performance parameters of this model evaluated on the basis of the external test set were as follows: accuracy=0.70, sensitivity=0.73 and specificity=0.66. A projection of the most frequently occurring bins on the standard coordinate space suggested a toxicophore composed of an aromatic ring with a centroid 3.5-7.5Å distant from an amino-group. The presence of a second aromatic ring separated by a 4-5Å spacer from the first ring and at a distance of between 5.5Å and 7Å from the amino-group was also associated with a PLD+ effect. These models provide comparable predictive performance to previously reported models for PLD with the added benefit of being based entirely on non-confidential, publicly available training data and with good predictive performance when tested in a rigorous, external validation exercise.

Partial least square and k-nearest neighbor algorithms for improved 3D quantitative spectral data-activity relationship consensus modeling of acute toxicity.

A diverse set of 154 chemicals that included US Food and Drug Administration-regulated compounds tested for their aquatic toxicity in Daphnia magna were modeled by a 3-dimensional quantitative spectral data-activity relationship (3D-QSDAR). Two distinct algorithms, partial least squares (PLS) and Tanimoto similarity-based k-nearest neighbors (KNN), were used to process bin occupancy descriptor matrices obtained after tessellation of the 3D-QSDAR space into regularly sized bins. The performance of models utilizing bins ranging in size from 2 ppm × 2 ppm × 0.5 Å to 20 ppm × 20 ppm × 2.5 Å was explored. Rigorous quality-control criteria were imposed: 1) 100 randomized 20% hold-out test sets were generated and the average R(2) test of the respective models was used as a measure of their performance, and 2) a Y-scrambling procedure was used to identify chance correlations. A consensus between the best-performing composite PLS model using 0.5 Å × 14 ppm × 14 ppm bins and 10 latent variables (average R(2) test = 0.770) and the best composite KNN model using 0.5 Å × 8 ppm × 8 ppm and 2 neighbors (average R(2) test = 0.801) offered an improvement of about 7.5% (R(2) test consensus = 0.845). Projection of the most frequently occurring bins on the standard coordinate space indicated that the presence of a primary or secondary amino group-substituted aromatic systems-would result in an increased toxic effect in Daphnia. The presence of a second aromatic ring with highly electronegative substituents 5 Å to 7 Å apart from the first ring would lead to a further increase in toxicity.

Complementary PLS and KNN algorithms for improved 3D-QSDAR consensus modeling of AhR binding.

Multiple validation techniques (Y-scrambling, complete training/test set randomization, determination of the dependence of R2test on the number of randomization cycles, etc.) aimed to improve the reliability of the modeling process were utilized and their effect on the statistical parameters of the models was evaluated. A consensus partial least squares (PLS)-similarity based k-nearest neighbors (KNN) model utilizing 3D-SDAR (three dimensional spectral data-activity relationship) fingerprint descriptors for prediction of the log(1/EC50) values of a dataset of 94 aryl hydrocarbon receptor binders was developed. This consensus model was constructed from a PLS model utilizing 10 ppm x 10 ppm x 0.5 Å bins and 7 latent variables (R2test of 0.617), and a KNN model using 2 ppm x 2 ppm x 0.5 Å bins and 6 neighbors (R2test of 0.622). Compared to individual models, improvement in predictive performance of approximately 10.5% (R2test of 0.685) was observed. Further experiments indicated that this improvement is likely an outcome of the complementarity of the information contained in 3D-SDAR matrices of different granularity. For similarly sized data sets of Aryl hydrocarbon (AhR) binders the consensus KNN and PLS models compare favorably to earlier reports. The ability of 3D-QSDAR (three dimensional quantitative spectral data-activity relationship) to provide structural interpretation was illustrated by a projection of the most frequently occurring bins on the standard coordinate space, thus allowing identification of structural features related to toxicity.

¹³C NMR-distance matrix descriptors: optimal abstract 3D space granularity for predicting estrogen binding.

An improved three-dimensional quantitative spectral data-activity relationship (3D-QSDAR) methodology was used to build and validate models relating the activity of 130 estrogen receptor binders to specific structural features. In 3D-QSDAR, each compound is represented by a unique fingerprint constructed from (13)C chemical shift pairs and associated interatomic distances. Grids of different granularity can be used to partition the abstract fingerprint space into congruent "bins" for which the optimal size was previously unexplored. For this purpose, the endocrine disruptor knowledge base data were used to generate 50 3D-QSDAR models with bins ranging in size from 2 ppm × 2 ppm × 0.5 Å to 20 ppm × 20 ppm × 2.5 Å, each of which was validated using 100 training/test set partitions. Best average predictivity in terms of R(2)test was achieved at 10 ppm ×10 ppm × Z Å (Z = 0.5, ..., 2.5 Å). It was hypothesized that this optimum depends on the chemical shifts' estimation error (±4.13 ppm) and the precision of the calculated interatomic distances. The highest ranked bins from partial least-squares weights were found to be associated with structural features known to be essential for binding to the estrogen receptor.

A computational study of the binding of 3-(arylidene) anabaseines to two major brain nicotinic acetylcholine receptors and to the acetylcholine binding protein.

Nicotinic acetylcholine receptors (nAChRs) have become targets for drug development in recent years. 3-(2,4-dimethoxybenzylidene)-anabaseine (DMXBA), which selectively stimulates the alpha7 nAChR, has been shown to alleviate some cognitive deficits associated with schizophrenia. In this paper we report an analysis of the interactions between 47 arylidene-anabaseines (including 45 benzylidene-anabaseines) and rat brain alpha7 and alpha4beta2 nicotinic acetylcholine receptors, using three different modeling techniques, namely 2D-QSAR, 3D-QSAR and molecular docking to the Aplysia californica acetylcholine binding protein (AChBP), a water soluble, homomeric nAChR surrogate receptor with a known crystal structure. Our investigation indicates the importance of: (1) the nitrogen atom of the tetrahydropyridyl (THP) ring for hydrogen bond formation; (2) pi-pi interactions between the aromatic rings of the ligands and the nAChBP binding site; (3) molecular surface recognition expressed in terms of steric complimentarity. On the basis of the 3D-QSAR results, bulky substituents at positions 2 (and due to the rotational freedom also at position 6) and 4 of the benzylidene moiety, with highly electronegative atoms projecting approximately 3-3.5A away from the benzylidene ring at position 4 seem optimal for enhancing binding affinity to the alpha7 nAChR.

Correlation of the photolysis half-lives of polychlorinated dibenzo-p-dioxins and dibenzofurans with molecular structure.

The photolysis half-lives of 70 polychlorinated dibenzo-p-dioxins and dibenzofurans are correlated with their molecular structures by a QSPR model (R(2) = 0.72) comprising three bond-energy-related descriptors. The photodegradation depends on the stability of the aromatic system and the C-O and C-C bond strengths. Model validation utilized leave-one-out (R(2) = 0.69), leave-many-out (R(2) = 0.72), and scrambling (R(2) = 0.19) procedures. Our results allow estimation of the photolysis half-lives of the remaining possible 140 PCDDs and PCDFs congeners.

Quantitative structure-activity relationship (QSAR) modeling of EC50 of aquatic toxicities for Daphnia magna.

The experimental EC(50) toxicities toward Daphnia magna for a series of 130 benzoic acids, benzaldehydes, phenylsulfonyl acetates, cycloalkane-carboxylates, benzanilides, and other esters were studied using the Best multilinear regression algorithm (BMLR) implemented in CODESSA. A modified quantitative structure-activity relationships (QSAR) procedure was applied guaranteeing the stability and reproducibility of the results. Separating the initial data set into training and test subsets generated three independent models with an average R(2) of .735. A five-descriptor general model including all 130 compounds, constructed using the descriptors found effective for the independent subsets, was characterized by the following statistical parameters: R(2) = .712; R(2)(cv) = .676; F = 61.331; s(2) = 0.6. The removal of two extreme outliers improved significantly the statistical parameters: R(2) = .759; R(2)(cv) = .728; F = 77.032; s(2) = 0.499. The sensitivity of the general model to chance correlations was estimated by applying a scrambling procedure involving 20 randomizations of the original property values. The resulting R(2) = .192 demonstrated the high robustness of the model proposed. The descriptors appearing in the obtained models are related to the biochemical nature of the adverse effects. An additional study of the EC(50)/LC(50) relationship for a series of 28 compounds (part of our general data set) revealed that these endpoints correlated with R(2) = .98.

QSAR modeling of the antifungal activity against Candida albicans for a diverse set of organic compounds.

The molecular structures of 83 diverse organic compounds are correlated by a quantitative structure-activity relationship (QSAR) to their minimum inhibitor concentrations (MIC expressed as log(1/MIC)), involving 6 descriptors with R(2)=0.788, F=47.140, s(2)=0.130. A novel QSAR development technique is utilized combining advantages of the two frequently applied methods. The topological, electronic, geometrical, and hybrid type descriptors for the compounds were calculated by CODESSA PRO software.

QSPR modeling of UV absorption intensities.

Literature UV absorption intensities at 260 nm and 25 degrees C in water of a diverse set of 805 organic compounds when analyzed by CODESSA Pro software using an initial pool of 800 + descriptors provide a significant QSPR correlation (R (2) = 0.692). Concurrently, a neural networks approach was used to develop a corresponding nonlinear model. The descriptors appearing in these models are discussed with respect to the physical nature of the UV absorption phenomenon.