Comparison of multilabel and single-label classification applied to the prediction of the isoform specificity of cytochrome p450 substrates.

Each drug can potentially be metabolized by different CYP450 isoforms. In the development of new drugs, the prediction of the metabolic fate is important to prevent drug-drug interactions. In the present study, a collection of 580 CYP450 substrates is deeply analyzed by applying multi- and single-label classification strategies, after the computation and selection of suitable molecular descriptors. Cross-training with support vector machine, multilabel k-nearest-neighbor and counterpropagation neural network modeling methods were used in the multilabel approach, which allows one to classify the compounds simultaneously in multiple classes. In the single-label models, automatic variable selection was combined with various cross-validation experiments and modeling techniques. Moreover, the reliability of both multi- and single-label models was assessed by the prediction of an external test set. Finally, the predicted results of the best models were compared to show that, even if the models present similar performances, the multilabel approach more coherently reflects the real metabolism information.

Exploring potency and selectivity receptor antagonist profiles using a multilabel classification approach: the human adenosine receptors as a key study.

Nowadays, in medicinal chemistry adenosine receptors represent some of the most studied targets, and there is growing interest on the different adenosine receptor (AR) subtypes. The AR subtypes selectivity is highly desired in the development of potent ligands to achieve the therapeutic success. So far, very few ligand-based strategies have been investigated to predict the receptor subtypes selectivity. In the present study, we have carried out a novel application of the multilabel classification approach by combining our recently reported autocorrelated molecular descriptors encoding for the molecular electrostatic potential (autoMEP) with support vector machines (SVMs). Three valuable models, based on decreasing thresholds of potency, have been generated as in series quantitative sieves for the simultaneous prediction of the hA(1)R, hA(2A)R, hA(2B)R, and hA(3)R subtypes potency profile and selectivity of a large collection, more than 500, of known inverse agonists such as xanthine, pyrazolo-triazolo-pyrimidine, and triazolo-pyrimidine analogues. The robustness and reliability of our multilabel classification models were assessed by predicting an internal test set. Finally, we have applied our strategy to 13 newly synthesized pyrazolo-triazolo-pyrimidine derivatives inferring their full adenosine receptor potency spectrum and hAR subtypes selectivity profile.

Combining selectivity and affinity predictions using an integrated Support Vector Machine (SVM) approach: An alternative tool to discriminate between the human adenosine A(2A) and A(3) receptor pyrazolo-triazolo-pyrimidine antagonists binding sites.

G Protein-coupled receptors (GPCRs) selectivity is an important aspect of drug discovery process, and distinguishing between related receptor subtypes is often the key to therapeutic success. Nowadays, very few valuable computational tools are available for the prediction of receptor subtypes selectivity. In the present study, we present an alternative application of the Support Vector Machine (SVM) and Support Vector Regression (SVR) methodologies to simultaneously describe both A(2A)R versus A(3)R subtypes selectivity profile and the corresponding receptor binding affinities. We have implemented an integrated application of SVM-SVR approach, based on the use of our recently reported autocorrelated molecular descriptors encoding for the Molecular Electrostatic Potential (autoMEP), to simultaneously discriminate A(2A)R versus A(3)R antagonists and to predict their binding affinity to the corresponding receptor subtype of a large dataset of known pyrazolo-triazolo-pyrimidine analogs. To validate our approach, we have synthetized 51 new pyrazolo-triazolo-pyrimidine derivatives anticipating both A(2A)R/A(3)R subtypes selectivity and receptor binding affinity profiles.

Prediction of the aqueous solvation free energy of organic compounds by using autocorrelation of molecular electrostatic potential surface properties combined with response surface analysis.

Several quantitative structure-property relationship (QSPR) approaches have been explored for the prediction of aqueous solubility or aqueous solvation free energies, DeltaG(sol), as crucial parameter affecting the pharmacokinetic profile and toxicity of chemical compounds. It is mostly accepted that aqueous solvation free energies can be expressed quantitatively in terms of properties of the molecular surface electrostatic potentials of the solutes. In the present study we have introduced autocorrelation molecular electrostatic potential (autoMEP) vectors in combination with nonlinear response surface analysis (RSA) as alternative 3D-QSPR strategy to evaluate the aqueous solvation free energy of organic compounds. A robust QSPR model (r(cv)=0.93) has been obtained by using a collection of 248 organic chemicals. An external test set based on 23 molecules confirmed the good predictivity of the autoMEP/RSA model suggesting its further applicability in the in silico prediction of water solubility of large organic compound libraries.

Support Vector Machine (SVM) as Alternative Tool to Assign Acute Aquatic Toxicity Warning Labels to Chemicals.

Quantitative structure-activity relationship (QSAR) analysis has been frequently utilized as a computational tool for the prediction of several eco-toxicological parameters including the acute aquatic toxicity. In the present study, we describe a novel integrated strategy to describe the acute aquatic toxicity through the combination of both toxicokinetic and toxicodynamic behaviors of chemicals. In particular, a robust classification model (TOXclass) has been derived by combining Support Vector Machine (SVM) analysis with three classes of toxicokinetic-like molecular descriptors: the autocorrelation molecular electrostatic potential (autoMEP) vectors, Sterimol topological descriptors and logP(o/w) property values. TOXclass model is able to assign chemicals to different levels of acute aquatic toxicity, providing an appropriate answer to the new regulatory requirements. Moreover, we have extended the above mentioned toxicokinetic-like descriptor set with a more toxicodynamic-like descriptors, as for example HOMO and LUMO energies, to generate a valuable SVM classifier (MOAclass) for the prediction of the mode of action (MOA) of toxic chemicals. As preliminary validation of our approach, the toxicokinetic (TOXclass) and the toxicodynamic (MOAclass) models have been applied in series to inspect both aquatic toxicity hazard and mode of action of 296 chemical substances with unknown or uncertain toxicodynamic information to assess the potential ecological risk and the toxic mechanism.

Linear and nonlinear 3D-QSAR approaches in tandem with ligand-based homology modeling as a computational strategy to depict the pyrazolo-triazolo-pyrimidine antagonists binding site of the human adenosine A2A receptor.

The integration of ligand- and structure-based strategies might sensitively increase the success of drug discovery process. We have recently described the application of Molecular Electrostatic Potential autocorrelated vectors (autoMEPs) in generating both linear (Partial Least-Square, PLS) and nonlinear (Response Surface Analysis, RSA) 3D-QSAR models to quantitatively predict the binding affinity of human adenosine A3 receptor antagonists. Moreover, we have also reported a novel GPCR modeling approach, called Ligand-Based Homology Modeling (LBHM), as a tool to simulate the conformational changes of the receptor induced by ligand binding. In the present study, the application of both linear and nonlinear 3D-QSAR methods and LBHM computational techniques has been used to depict the hypothetical antagonist binding site of the human adenosine A2A receptor. In particular, a collection of 127 known human A2A antagonists has been utilized to derive two 3D-QSAR models (autoMEPs/PLS&RSA). In parallel, using a rhodopsin-driven homology modeling approach, we have built a model of the human adenosine A2A receptor. Finally, 3D-QSAR and LBHM strategies have been utilized to predict the binding affinity of five new human A2A pyrazolo-triazolo-pyrimidine antagonists finding a good agreement between the theoretical and the experimental predictions.