Legitimate Utilization of Large Descriptor Pools for QSPR/QSAR Models

New Antagonists of the Membrane Androgen Receptor OXER1 from the ZINC Natural Product Database

OXER1 (oxoeicosanoid receptor 1) was deorphanized in 1993 and found to be the specific receptor for the arachidonic acid metabolite 5-oxo-ETE. Recently, we have reported that androgen binds to this receptor also, being a membrane androgen receptor, triggering a number of its membrane-mediated actions (cell migration, apoptosis, cell proliferation, Ca²⁺ movements). In addition, our previous work suggested that a number of natural monomeric and oligomeric polyphenols interact with OXER1, acting similar to testosterone. Here, we interrogated the natural product chemical space and identified nine polyphenolic molecules with interesting in silico pharmacological activities as putative OXER1 antagonists. The molecule with the best pharmacokinetic-pharmacodynamic properties (ZINC15959779) was purchased and tested on OXER1, in prostate cancer cell cultures. It showed that it has actions similar to those of testosterone in inhibiting cAMP, while it had no action in intracellular Ca²⁺ mobilization or actin cytoskeleton rearrangement/migration. These results are discussed under the prism of structure-activity relationships and in silico models of the OXER1 binding groove. We suggest that these compounds, together with the previously reported (poly)phenolic compounds, can be lead structures for the exploration of the anti-inflammatory and antiproliferative effects of OXER1 antagonists.

Development and interpretation of a QSAR model for in vitro breast cancer (MCF-7) cytotoxicity of 2-phenylacrylonitriles

The Arylhydrocarbon Receptor (AhR), a member of the Per-ARNT-SIM transcription factor family, has been as a potential new target to treat breast cancer sufferers. A series of 2-phenylacrylonitriles targeting AhR has been developed that have shown promising and selective activity against cancerous cell lines while sparing normal non-cancerous cells. A quantitative structure–activity relationship (QSAR) modeling approach was pursued in order to generate a predictive model for cytotoxicity to support ongoing synthetic activities and provide important structure-activity information for new structure design. Recent work conducted by us has identified a number of compounds that exhibited false positive cytotoxicity values in the standard MTT assay. This work describes a good quality model that not only predicts the activity of compounds in the MCF-7 breast cancer cell line, but was also able to identify structures that subsequently gave false positive values in the MTT assay by identifying compounds with aberrant biological behavior. This work not only allows the design of future breast cancer cytotoxic activity in vitro, but allows the avoidance of the synthesis of those compounds anticipated to result in anomalous cytotoxic behavior, greatly enhancing the design of such compounds.

Principles of QSAR Modeling

At the end of her academic career, the author summarizes the main aspects of QSAR modeling, giving comments and suggestions according to her 23 years' experience in QSAR research on environmental topics. The focus is mainly on Multiple Linear Regression, particularly Ordinary Least Squares, using a Genetic Algorithm for variable selection from various theoretical molecular descriptors, but the comments can be useful also for other QSAR methods. The need for rigorous validation, also external, and for applicability domain check to guarantee predictivity and reliability of QSAR models is particularly highlighted. The commented approach is the “predictive” one, based on chemometrics, and is usefully applied to the prioritization of environmental pollutants. All the discussed points and the author's ideas are implemented in the software QSARINS, as a legacy to the QSAR community.

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.

QSAR modelling of integrin antagonists using enhanced Bayesian regularised genetic neural networks

Bayesian regularised genetic neural network (BRGNN) has been used for modelling the inhibition activity of 141 biphenylalanine derivatives as integrin antagonists. Three local pattern search (PS) methods, simulated annealing and threshold acceptance were combined with BRGNN in the form of a hybrid genetic algorithm (HGA). The results obtained revealed that PS is a suitable method for improving the ability of BRGNN to break out from the local minima. The proposed HGA technique is able to retrieve important variables from complex systems and nonlinear search spaces for optimisation. Two models with 8-3-1 artificial neural network (ANN) architectures were developed for describing α₄β₇ and α₄β₁ modulatory activities of integrin antagonists. Monte Carlo cross-validation was performed to validate the models and Q₂ values of 0.75 and 0.74 were obtained for α₄β₇ and α₄β₁ inhibitory activities, respectively. The scrambling technique was used for sensitivity analysis of descriptors appearing in ANN models. Frequencies of repetition and sensitivity analysis of molecular descriptors revealed that 3D-Morse descriptors are influential factors for describing α₄β₇ inhibitory activity, while in the case of α₄β₁ inhibitory activity, the Randic shape index, the lowest eigenvalue of the Burden matrix and the number of rotatable bonds are important parameters.

A Quantitative Structure-Property Relationship (QSPR) Study of aliphatic alcohols by the method of dividing the molecular structure into substructure

A quantitative structure–property relationship (QSPR) analysis of aliphatic alcohols is presented. Four physicochemical properties were studied: boiling point (BP), n-octanol–water partition coefficient (lg P_OW), water solubility (lg W) and the chromatographic retention indices (RI) on different polar stationary phases. In order to investigate the quantitative structure–property relationship of aliphatic alcohols, the molecular structure ROH is divided into two parts, R and OH to generate structural parameter. It was proposed that the property is affected by three main factors for aliphatic alcohols, alkyl group R, substituted group OH, and interaction between R and OH. On the basis of the polarizability effect index (PEI), previously developed by Cao, the novel molecular polarizability effect index (MPEI) combined with odd-even index (OEI), the sum eigenvalues of bond-connecting matrix (SX_1CH) previously developed in our team, were used to predict the property of aliphatic alcohols. The sets of molecular descriptors were derived directly from the structure of the compounds based on graph theory. QSPR models were generated using only calculated descriptors and multiple linear regression techniques. These QSPR models showed high values of multiple correlation coefficient (R > 0.99) and Fisher-ratio statistics. The leave-one-out cross-validation demonstrated the final models to be statistically significant and reliable.

QSAR model for the prediction of bio-concentration factor using aqueous solubility and descriptors considering various electronic effects

The in silico modelling of bio-concentration factor (BCF) is of considerable interest in environmental sciences, because it is an accepted indicator for the accumulation potential of chemicals in organisms. Numerous QSAR models have been developed for the BCF, and the majority utilize the octanol/water partition coefficient (log P) to account for the penetration characteristics of the chemicals. The present work used descriptors from a variety of software packages for the development of a multi-linear regression model to estimate BCF. The modelled data set of 473 diverse compounds covers a wide range of log BCF values. In the proposed QSAR model, most of the variation is described by the calculated solubility in water. Other contributing descriptors describe, for instance, hydrophobic surface area, hydrogen bonding and other electronic effects. The model was validated internally by using a variety of statistical approaches. Two external validations were also performed. For the former validation, a subset from the same data source was used. The 2nd external validation was based on an independent data set collected from different resources. All validations showed the consistency of the model. The applicability domain of the model was discussed and described and a thorough outlier analysis was performed.

Estimating the toxicities of organic chemicals in activated sludge process

The experimental logEC50 toxicity values of 104 compounds causing bioluminescent repression of the bacterium strain Pseudomonas isolated from an industrial wastewater were studied. Using the Best Multilinear Regression method implemented in CODESSA PRO, models with up to 8 theoretical descriptors were obtained. Utilizing a rigorous descriptor selection and validation procedure a reliable QSAR model with four parameters was selected as best. The proposed model emphasizes the importance of the halogen atoms presented in each compound, the possibility of H-bond formation and the flexibility and degree of branching of the molecules. As pointed out by many researchers, the contribution of the octanol-water partition coefficient to the explanation of the toxicity effect was also found to be significant. In addition, the model currently proposed was compared to those reported earlier and its advantages were discussed in detail.