Multi-target spectral moment QSAR versus ANN for antiparasitic drugs against different parasite species

There are many of pathogen parasite species with different susceptibility profile to antiparasitic drugs. Unfortunately, almost QSAR models predict the biological activity of drugs against only one parasite species. Consequently, predicting the probability with which a drug is active against different species with a single unify model is a goal of the major importance. In so doing, we use Markov Chains theory to calculate new multi-target spectral moments to fit a QSAR model that predict by the first time a mt-QSAR model for 500 drugs tested in the literature against 16 parasite species and other 207 drugs no tested in the literature using spectral moments. The data was processed by linear discriminant analysis (LDA) classifying drugs as active or non-active against the different tested parasite species. The model correctly classifies 311 out of 358 active compounds (86.9%) and 2328 out of 2577 non-active compounds (90.3%) in training series. Overall training performance was 89.9%. Validation of the model was carried out by means of external predicting series. In these series the model classified correctly 157 out 190, 82.6% of antiparasitic compounds and 1151 out of 1277 non-active compounds (90.1%). Overall predictability performance was 89.2%. In addition we developed four types of non Linear Artificial neural networks (ANN) and we compared with the mt-QSAR model. The improved ANN model had an overall training performance was 87%. The present work report the first attempts to calculate within a unify framework probabilities of antiparasitic action of drugs against different parasite species based on spectral moment analysis.     

Stochastic entropy QSAR for the in silico discovery of anticancer compounds: prediction, synthesis, and in vitro assay of new purine carbanucleosides

A Markov model based QSAR is introduced for the rational selection of anticancer compounds. The model discriminates 90.3% of 226 structurally heterogeneous anticancer/non-anticancer compounds in training series. External validation series were used to validate the model; the 91.8% containing 85 compounds, not considered to fit the model, were correctly classified. The model developed is afterwards used in a simulation of a virtual search for anticancer compounds never considered either in training or in predicting series. The 87.7% of the 213 anticancer compounds used in this simulated search were correctly classified. The model also shows high values for specificity (0.89), sensitivity (0.91), and Mathews correlation coefficient (0.79). In addition, the present model compares better-to-similar with respect to other four models elsewhere reported if one takes into consideration 26 comparison parameters. Finally, we exemplify the use of the model in practice with the design of a new series of carbanucleosides. The compounds evaluated with the model were synthesized and experimentally assayed for their antitumor effects on the proliferation of murine leukemia cells (L1210/0) and human T-lymphocyte cells (CEM/0 and Molt4/C8). The more interesting activity was detected for the compound 5a with a predicted probability of 80.2% and IC(50) = 27.0, 27.2, and 29.4 microM, respectively, against the above-mentioned cellular lines. These values are comparable to those for the control compound Ara-A.     

Unified QSAR approach to antimicrobials. Part 2: predicting activity against more than 90 different species in order to halt antibacterial resistance

There are many different kinds of pathogenic bacteria species with very different susceptibility profiles to different antibacterial drugs. One limitation of QSAR models is that they consider the biological activity of drugs against only one species of bacteria. In a previous paper, we developed a unified Markov model to describe the biological activity of different drugs tested in the literature against some antimicrobial species. Consequently, predicting the probability with which a drug is active against different species of bacteria with a single unified model is a goal of major importance. The work described here develops the unified Markov model to describe the biological activity of more than 70 drugs from the literature tested against 96 species of bacteria. We applied linear discriminant analysis (LDA) to classify drugs as active or inactive against the different tested bacterial species. The model correctly classified 199 out of 237 active compounds (83.9%) and 168 out of 200 inactive compounds (84%). Overall training predictability was 84% (367 out of 437 cases). Validation of the model was carried out using an external predicting series, with the model classifying 202 out of 243 (i.e., 83.13%) of the compounds. In order to show how the model functions in practice, a virtual screening was carried out and the model recognized as active 84.5% (480 out of 568) antibacterial compounds not used in the training or predicting series. The current study is an attempt to calculate within a unified framework the probabilities of antibacterial action of drugs against many different species.     

Unified QSAR approach to antimicrobials. Part 3: first multi-tasking QSAR model for input-coded prediction, structural back-projection, and complex networks clustering of antiprotozoal compounds

Several pathogen parasite species show different susceptibilities to different antiparasite drugs. Unfortunately, almost all structure-based methods are one-task or one-target Quantitative Structure-Activity Relationships (ot-QSAR) that predict the biological activity of drugs against only one parasite species. Consequently, multi-tasking learning to predict drugs activity against different species by a single model (mt-QSAR) is vitally important. In the two previous works of the present series we reported two single mt-QSAR models in order to predict the antimicrobial activity against different fungal (Bioorg. Med. Chem.2006, 14, 5973-5980) or bacterial species (Bioorg. Med. Chem.2007, 15, 897-902). These mt-QSARs offer a good opportunity (unpractical with ot-QSAR) to construct drug-drug similarity Complex Networks and to map the contribution of sub-structures to function for multiple species. These possibilities were unattended in our previous works. In the present work, we continue this series toward other important direction of chemotherapy (antiparasite drugs) with the development of an mt-QSAR for more than 500 drugs tested in the literature against different parasites. The data were processed by Linear Discriminant Analysis (LDA) classifying drugs as active or non-active against the different tested parasite species. The model correctly classifies 212 out of 244 (87.0%) cases in training series and 207 out of 243 compounds (85.4%) in external validation series. In order to illustrate the performance of the QSAR for the selection of active drugs we carried out an additional virtual screening of antiparasite compounds not used in training or predicting series; the model recognized 97 out of 114 (85.1%) of them. We also give the procedures to construct back-projection maps and to calculate sub-structures contribution to the biological activity. Finally, we used the outputs of the QSAR to construct, by the first time, a multi-species Complex Networks of antiparasite drugs. The network predicted has 380 nodes (compounds), 634 edges (pairs of compounds with similar activity). This network allows us to cluster different compounds and identify on average three known compounds similar to a new query compound according to their profile of biological activity. This is the first attempt to calculate probabilities of antiparasitic action of drugs against different parasites.     

TOPS-MODE based QSARs derived from heterogeneous series of compounds. Applications to the design of new anti-inflammatory compounds

A new application of TOPological Sub-structural MOlecular DEsign (TOPS-MODE) was carried out in anti-inflammatory compounds using computer-aided molecular design. Two series of compounds, one containing anti-inflammatory and the other containing nonanti-inflammatory compounds were processed by a k-means cluster analysis in order to design the training and prediction sets. A linear classification function to discriminate the anti-inflammatory from the inactive compounds was developed. The model correctly and clearly classified 88% of active and 91% of inactive compounds in the training set. More specifically, the model showed a good global classification of 90%, that is, (399 cases out of 441). While in the prediction set, they showed an overall predictability of 88% and 84% for active and inactive compounds, being the global percentage of good classification of 85%. Furthermore this paper describes a fragment analysis in order to determine the contribution of several fragments towards anti-inflammatory property, also the present of halogens in the selected fragments were analyzed. It seems that the present TOPS-MODE based QSAR is the first alternate general 'in silico' technique to experimentation in anti-inflammatory discovery.     

3D QSAR Markov model for drug-induced eosinophilia--theoretical prediction and preliminary experimental assay of the antimicrobial drug G1

The application of 3D-MEDNEs as a novel alternative technique to reduce the use of animal experimentation in toxicology in the early stages of medicinal chemistry research has been extended from agranulocytosis to chemically induced eosinophilia. Firstly, a heterogeneous series of organic compounds, which are classified either as eosinophilia inductors or noninductors, was collected. A linear discriminant analysis was subsequently used to obtain a QSTR that gave rise to a very good classification of 91.82% (110 chemicals within training series). Eosinophilia inductors (88.89%) composed the first group while the other one contained only harmless compounds (97.37%). The total predictability (88.1%) was tested by means of an external validation series (42 compounds). The model correctly classifies 88.89% of harmless compounds and 87.5% of toxic ones. Finally, comparison of predicted versus experimental results for G1 [2-bromo-5-(2-bromo-2-nitroethenyl)furan, which is a promising antibacterial-antifungal compound] illustrates the practical application of the method. A dose-dependent study of G1 (9.8-185.6 mg/Kg) at 48, 72 and 96 h after oral administration in rats is reported here for the first time. The study has shown that G1 does not affect the murine eosinophils count under these conditions--a situation in total agreement with the model prediction.     

QSAR-based solubility model for drug-like compounds

Solubility plays a very important role in the selection of compounds for drug screening. In this context, a QSAR model was developed for predicting water solubility of drug-like compounds. First, a set of relevant parameters for establishing a drug-like chemical space was defined. The comparison of chemical structures from the FDAMDD and PHYSPROP databases allowed the selection of properties that were more efficient in discriminating drug-like compounds from other chemicals. These filters were later on applied to the PHYSPROP database and 1174 chemicals fulfilling these criteria and with experimental solubility information available at 25 °C were retained. Several QSAR solubility models were developed from this set of compounds, and the best one was selected based on the accuracy of correct classifications obtained for randomly chosen training and validation subsets. Further validation of the model was performed with a set of 102 drugs for which experimental solubility data have been recently reported. A good agreement between the predictions and the experimental values confirmed the reliability of the QSAR model.     

Markov entropy backbone electrostatic descriptors for predicting proteins biological activity

The spherical truncation of electrostatic interactions between aminoacids makes it possible to break down long-range spatial electrostatic interactions, resulting in short-range interactions. As a result, a Markov Chain model may be used to calculate the probabilities with which the effect of a given interaction reaches aminoacids at different distances within the backbone. The entropies of a Markov Chain model of this type may then be used to codify information about the spatial distribution of charges in the protein used in this study exploring the structure-activity relationship. In this paper, a linear discriminant analysis is reported, which correctly classified 92.3% of 26 under investigation in training and leave-one-out cross validation, purely for illustrative purposes. Classification was carried out for three possible activities: lysozymes, dihydrofolate reductases, and alcohol dehydrogenases. The discriminant analysis equations were contracted into two canonical roots. These simple canonical roots have high regression coefficients (R(c1)=0.903 and R(c2)=0.70). Root1 explains the biological activity of alcohol dehydrogenases while Root2 discriminates between lysozymes and dihydrofolate reductases. It was possible to profile the effect of core, middle, and surface aminoacids on biological activity. In contrast, a model considering classic physicochemical parameters such as: polarizability, refractivity, and partition coefficient classify correctly only the 80.8% of the proteins.     

TOPS-MODE model of multiplexing neuroprotective effects of drugs and experimental-theoretic study of new 1,3-rasagiline derivatives potentially useful in neurodegenerative diseases

The interest on computational techniques for the discovery of neuroprotective drugs has increased due to recent fail of important clinical trials. In fact, there is a huge amount of data accumulated in public databases like CHEMBL with respect to structurally heterogeneous series of drugs, multiple assays, drug targets, and model organisms. However, there are no reports of multi-target or multiplexing Quantitative Structure–Property Relationships (mt-QSAR/mx-QSAR) models of these multiplexing assay outcomes reported in CHEMBL for neurotoxicity/neuroprotective effects of drugs. Accordingly, in this paper we develop the first mx-QSAR model for multiplexing assays of neurotoxicity/neuroprotective effects of drugs. We used the method TOPS-MODE to calculate the structural parameters of drugs. The best model found correctly classified 4393 out of 4915 total cases in both training and validation. This is representative of overall train and validation Accuracy, Sensitivity, and Specificity values near to 90%, 98%, and 80%, respectively. This dataset includes multiplexing assay endpoints of 2217 compounds. Every one compound was assayed in at least one out of 338 assays, which involved 148 molecular or cellular targets and 35 standard type measures in 11 model organisms (including human). The second aim of this work is the exemplification of the use of the new mx-QSAR model with a practical case of study. To this end, we obtained again by organic synthesis and reported, by the first time, experimental assays of the new 1,3-rasagiline derivatives 3 different tests: assay (1) in absence of neurotoxic agents, (2) in the presence of glutamate, and (3) in the presence of H₂O₂. The higher neuroprotective effects found for each one of these assays were for the stereoisomers of compound 7: compound 7b with protection = 23.4% in assay (1) and protection = 15.2% in assay (2); and for compound 7a with protection = 46.2% in assay (3). Interestingly, almost all compounds show protection values >10% in assay (3) but not in the other 2 assays. After that, we used the mx-QSAR model to predict the more probable response of the new compounds in 559 unique pharmacological tests not carried out experimentally. The results obtained are very significant because they complement the pharmacological studies of these promising rasagiline derivatives. This work paves the way for further developments in the multi-target/multiplexing screening of large libraries of compounds potentially useful in the treatment of neurodegenerative diseases.