Prediction of multi-target networks of neuroprotective compounds with entropy indices and synthesis, assay, and theoretical study of new asymmetric 1,2-rasagiline carbamates.

In a multi-target complex network, the links (L(ij)) represent the interactions between the drug (d(i)) and the target (t(j)), characterized by different experimental measures (K(i), K(m), IC50, etc.) obtained in pharmacological assays under diverse boundary conditions (c(j)). In this work, we handle Shannon entropy measures for developing a model encompassing a multi-target network of neuroprotective/neurotoxic compounds reported in the CHEMBL database. The model predicts correctly >8300 experimental outcomes with Accuracy, Specificity, and Sensitivity above 80%-90% on training and external validation series. Indeed, the model can calculate different outcomes for >30 experimental measures in >400 different experimental protocolsin relation with >150 molecular and cellular targets on 11 different organisms (including human). Hereafter, we reported by the first time the synthesis, characterization, and experimental assays of a new series of chiral 1,2-rasagiline carbamate derivatives not reported in previous works. The experimental tests included: (1) assay in absence of neurotoxic agents; (2) in the presence of glutamate; and (3) in the presence of H2O2. Lastly, we used the new Assessing Links with Moving Averages (ALMA)-entropy model to predict possible outcomes for the new compounds in a high number of pharmacological tests not carried out experimentally.

New Markov-Shannon Entropy models to assess connectivity quality in complex networks: from molecular to cellular pathway, Parasite-Host, Neural, Industry, and Legal-Social networks.

Graph and Complex Network theory is expanding its application to different levels of matter organization such as molecular, biological, technological, and social networks. A network is a set of items, usually called nodes, with connections between them, which are called links or edges. There are many different experimental and/or theoretical methods to assign node-node links depending on the type of network we want to construct. Unfortunately, the use of a method for experimental reevaluation of the entire network is very expensive in terms of time and resources; thus the development of cheaper theoretical methods is of major importance. In addition, different methods to link nodes in the same type of network are not totally accurate in such a way that they do not always coincide. In this sense, the development of computational methods useful to evaluate connectivity quality in complex networks (a posteriori of network assemble) is a goal of major interest. In this work, we report for the first time a new method to calculate numerical quality scores S(L(ij)) for network links L(ij) (connectivity) based on the Markov-Shannon Entropy indices of order k-th (θ(k)) for network nodes. The algorithm may be summarized as follows: (i) first, the θ(k)(j) values are calculated for all j-th nodes in a complex network already constructed; (ii) A Linear Discriminant Analysis (LDA) is used to seek a linear equation that discriminates connected or linked (L(ij)=1) pairs of nodes experimentally confirmed from non-linked ones (L(ij)=0); (iii) the new model is validated with external series of pairs of nodes; (iv) the equation obtained is used to re-evaluate the connectivity quality of the network, connecting/disconnecting nodes based on the quality scores calculated with the new connectivity function. This method was used to study different types of large networks. The linear models obtained produced the following results in terms of overall accuracy for network reconstruction: Metabolic networks (72.3%), Parasite-Host networks (93.3%), CoCoMac brain cortex co-activation network (89.6%), NW Spain fasciolosis spreading network (97.2%), Spanish financial law network (89.9%) and World trade network for Intelligent & Active Food Packaging (92.8%). In order to seek these models, we studied an average of 55,388 pairs of nodes in each model and a total of 332,326 pairs of nodes in all models. Finally, this method was used to solve a more complicated problem. A model was developed to score the connectivity quality in the Drug-Target network of US FDA approved drugs. In this last model the θ(k) values were calculated for three types of molecular networks representing different levels of organization: drug molecular graphs (atom-atom bonds), protein residue networks (amino acid interactions), and drug-target network (compound-protein binding). The overall accuracy of this model was 76.3%. This work opens a new door to the computational reevaluation of network connectivity quality (collation) for complex systems in molecular, biomedical, technological, and legal-social sciences as well as in world trade and industry.

ANN multiplexing model of drugs effect on macrophages; theoretical and flow cytometry study on the cytotoxicity of the anti-microbial drug G1 in spleen.

Multiplexed biological assays provide multiple measurements of cellular parameters in the same test. In this work, we have trained and tested an Artificial Neural Network (ANN) model for the first time, in order to perform a multiplexing prediction of drugs effect on macrophage populations. In so doing, we have used the TOPS-MODE approach to calculate drug molecular descriptors and the software STATISTICA to seek different ANN models such as: Linear Neural Network (LNN), Radial Basis Function (RBF), Probabilistic Neural Networks (PNN) and Multi-Layer Perceptrons (MLP). The best model found was the LNN, which correctly classified 8258 out of 9000 (Accuracy = 93.0%) multiplexing assay endpoints of 7903 drugs (including both training and test series). Each endpoint corresponds to one out of 1418 assays, 36 molecular or cellular targets, 46 standard type measures, in two possible organisms (human and mouse). Secondly, we have determined experimentally, for the first time, the values of EC(50) = 11.41 µg/mL and Cytotoxicity = 27.1% for the drug G1 over Balb/C mouse spleen macrophages using flow cytometry. In addition, we have used the LNN model to predict the G1 activity in 1265 multiplexing assays not measured experimentally (including 152 cytotoxicity assay endpoints). Both experimental and theoretical results point out a low macrophage cytotoxicity of G1. This work breaks new ground for the 'in silico' multiplexing screening of large libraries of compounds. The results obtained are very significant because they complement the immunotoxicology studies of this important anti-microbial/anti-parasite drug.

MIND-BEST: Web server for drugs and target discovery; design, synthesis, and assay of MAO-B inhibitors and theoretical-experimental study of G3PDH protein from Trichomonas gallinae.

Many drugs with very different affinity to a large number of receptors are described. Thus, in this work, we selected drug-target pairs (DTPs/nDTPs) of drugs with high affinity/nonaffinity for different targets. Quantitative structure-activity relationship (QSAR) models become a very useful tool in this context because they substantially reduce time and resource-consuming experiments. Unfortunately, most QSAR models predict activity against only one protein target and/or they have not been implemented on a public Web server yet, freely available online to the scientific community. To solve this problem, we developed a multitarget QSAR (mt-QSAR) classifier combining the MARCH-INSIDE software for the calculation of the structural parameters of drug and target with the linear discriminant analysis (LDA) method in order to seek the best model. The accuracy of the best LDA model was 94.4% (3,859/4,086 cases) for training and 94.9% (1,909/2,012 cases) for the external validation series. In addition, we implemented the model into the Web portal Bio-AIMS as an online server entitled MARCH-INSIDE Nested Drug-Bank Exploration & Screening Tool (MIND-BEST), located at http://miaja.tic.udc.es/Bio-AIMS/MIND-BEST.php . This online tool is based on PHP/HTML/Python and MARCH-INSIDE routines. Finally, we illustrated two practical uses of this server with two different experiments. In experiment 1, we report for the first time a MIND-BEST prediction, synthesis, characterization, and MAO-A and MAO-B pharmacological assay of eight rasagiline derivatives, promising for anti-Parkinson drug design. In experiment 2, we report sampling, parasite culture, sample preparation, 2-DE, MALDI-TOF and -TOF/TOF MS, MASCOT search, 3D structure modeling with LOMETS, and MIND-BEST prediction for different peptides as new protein of the found in the proteome of the bird parasite Trichomonas gallinae, which is promising for antiparasite drug targets discovery.

NL MIND-BEST: a web server for ligands and proteins discovery--theoretic-experimental study of proteins of Giardia lamblia and new compounds active against Plasmodium falciparum.

There are many protein ligands and/or drugs described with very different affinity to a large number of target proteins or receptors. In this work, we selected Ligands or Drug-target pairs (DTPs/nDTPs) of drugs with high affinity/non-affinity for different targets. Quantitative Structure-Activity Relationships (QSAR) models become a very useful tool in this context to substantially reduce time and resources consuming experiments. Unfortunately most QSAR models predict activity against only one protein target and/or have not been implemented in the form of public web server freely accessible online to the scientific community. To solve this problem, we developed here a multi-target QSAR (mt-QSAR) classifier using the MARCH-INSIDE technique to calculate structural parameters of drug and target plus one Artificial Neuronal Network (ANN) to seek the model. The best ANN model found is a Multi-Layer Perceptron (MLP) with profile MLP 20:20-15-1:1. This MLP classifies correctly 611 out of 678 DTPs (sensitivity=90.12%) and 3083 out of 3408 nDTPs (specificity=90.46%), corresponding to training accuracy=90.41%. The validation of the model was carried out by means of external predicting series. The model classifies correctly 310 out of 338 DTPs (sensitivity=91.72%) and 1527 out of 1674 nDTP (specificity=91.22%) in validation series, corresponding to total accuracy=91.30% for validation series (predictability). This model favorably compares with other ANN models developed in this work and Machine Learning classifiers published before to address the same problem in different aspects. We implemented the present model at web portal Bio-AIMS in the form of an online server called: Non-Linear MARCH-INSIDE Nested Drug-Bank Exploration & Screening Tool (NL MIND-BEST), which is located at URL: http://miaja.tic.udc.es/Bio-AIMS/NL-MIND-BEST.php. This online tool is based on PHP/HTML/Python and MARCH-INSIDE routines. Finally we illustrated two practical uses of this server with two different experiments. In experiment 1, we report by first time Quantum QSAR study, synthesis, characterization, and experimental assay of antiplasmodial and cytotoxic activities of oxoisoaporphine alkaloids derivatives as well as NL MIND-BEST prediction of potential target proteins. In experiment 2, we report sampling, parasite culture, sample preparation, 2-DE, MALDI-TOF, and -TOF/TOF MS, MASCOT search, MM/MD 3D structure modeling, and NL MIND-BEST prediction for different peptides a new protein of the found in the proteome of the human parasite Giardia lamblia, which is promising for anti-parasite drug-targets discovery.

Theoretical study of GSK-3α: neural networks QSAR studies for the design of new inhibitors using 2D descriptors.

Glycogen synthase kinase-3 (GSK-3) targets encompass proteins implicated in AD and neurological disorders. The functions of GSK-3 and its implication in various human diseases have triggered an active search for potent and selective GSK-3 inhibitors. In this sense, QSAR could play an important role in studying these GSK-3 inhibitors. For this reason, we developed QSAR models for GSK-3α, linear discriminant analysis (LDA), and artificial neural networks (ANNs) from nearly 50,000 cases with more than 700 different GSK-3α inhibitors obtained from ChEMBL database server; in total we used more than 20,000 different molecules to develop the QSAR models. The model correctly classified 237 out of 275 active compounds (86.2%) and 14,870 out of 15,970 non-active compounds (93.2%) in the training series. The overall training performance was 93.0%. Validation of the model was carried out using an external predicting series. In these series, the model classified correctly 458 out of 549 (83.4%) compounds and 29,637 out of 31,927 non-active compounds (83.4%). The overall predictability performance was 92.7%. In this study, we propose three types of non-linear ANN as alternative to already existing models, such as LDA. Linear neural network: LNN: 236:236-1-1:1 which had an overall training performance of 96% proved to be the best model. In addition, we did a study of the different fragments of the molecules of the database to see which fragments had more influence in the activity. This can help design new inhibitors of GSK-3α. This study reports the attempts to calculate, within a unified framework probabilities of GSK-3α inhibitors against different molecules found in the literature.

Trypano-PPI: a web server for prediction of unique targets in trypanosome proteome by using electrostatic parameters of protein-protein interactions.

Trypanosoma brucei causes African trypanosomiasis in humans (HAT or African sleeping sickness) and Nagana in cattle. The disease threatens over 60 million people and uncounted numbers of cattle in 36 countries of sub-Saharan Africa and has a devastating impact on human health and the economy. On the other hand, Trypanosoma cruzi is responsible in South America for Chagas disease, which can cause acute illness and death, especially in young children. In this context, the discovery of novel drug targets in Trypanosome proteome is a major focus for the scientific community. Recently, many researchers have spent important efforts on the study of protein-protein interactions (PPIs) in pathogen Trypanosome species concluding that the low sequence identities between some parasite proteins and their human host render these PPIs as highly promising drug targets. To the best of our knowledge, there are no general models to predict Unique PPIs in Trypanosome (TPPIs). On the other hand, the 3D structure of an increasing number of Trypanosome proteins is reported in databases. In this regard, the introduction of a new model to predict TPPIs from the 3D structure of proteins involved in PPI is very important. For this purpose, we introduced new protein-protein complex invariants based on the Markov average electrostatic potential xi(k)(R(i)) for amino acids located in different regions (R(i)) of i-th protein and placed at a distance k one from each other. We calculated more than 30 different types of parameters for 7866 pairs of proteins (1023 TPPIs and 6823 non-TPPIs) from more than 20 organisms, including parasites and human or cattle hosts. We found a very simple linear model that predicts above 90% of TPPIs and non-TPPIs both in training and independent test subsets using only two parameters. The parameters were (d)xi(k)(s) = |xi(k)(s(1)) - xi(k)(s(2))|, the absolute difference between the xi(k)(s(i)) values on the surface of the two proteins of the pairs. We also tested nonlinear ANN models for comparison purposes but the linear model gives the best results. We implemented this predictor in the web server named TrypanoPPI freely available to public at http://miaja.tic.udc.es/Bio-AIMS/TrypanoPPI.php. This is the first model that predicts how unique a protein-protein complex in Trypanosome proteome is with respect to other parasites and hosts, opening new opportunities for antitrypanosome drug target discovery.

Unified QSAR & network-based computational chemistry approach to antimicrobials. II. Multiple distance and triadic census analysis of antiparasitic drugs complex networks.

In the previous work, we reported a multitarget Quantitative Structure-Activity Relationship (mt-QSAR) model to predict drug activity against different fungal species. This mt-QSAR allowed us to construct a drug-drug multispecies Complex Network (msCN) to investigate drug-drug similarity (González-Díaz and Prado-Prado, J Comput Chem 2008, 29, 656). However, important methodological points remained unclear, such as follows: (1) the accuracy of the methods when applied to other problems; (2) the effect of the distance type used to construct the msCN; (3) how to perform the inverse procedure to study species-species similarity with multidrug resistance CNs (mdrCN); and (4) the implications and necessary steps to perform a substructural Triadic Census Analysis (TCA) of the msCN. To continue the present series with other important problem, we developed here a mt-QSAR model for more than 700 drugs tested in the literature against different parasites (predicting antiparasitic drugs). The data were processed by Linear Discriminate Analysis (LDA) and the model classifies correctly 93.62% (1160 out of 1239 cases) in training. The model validation was carried out by means of external predicting series; the model classified 573 out of 607, that is, 94.4% of cases. Next, we carried out the first comparative study of the topology of six different drug-drug msCNs based on six different distances such as Euclidean, Chebychev, Manhattan, etc. Furthermore, we compared the selected drug-drug msCN and species-species mdsCN with random networks. We also introduced here the inverse methodology to construct species-species msCN based on a mt-QSAR model. Last, we reported the first substructural analysis of drug-drug msCN using Triadic Census Analysis (TCA) algorithm.

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.

QSAR for RNases and theoretic-experimental study of molecular diversity on peptide mass fingerprints of a new Leishmania infantum protein.

The toxicity and low success of current treatments for Leishmaniosis determines the search of new peptide drugs and/or molecular targets in Leishmania pathogen species (L. infantum and L. major). For example, Ribonucleases (RNases) are enzymes relevant to several biologic processes; then, theoretical and experimental study of the molecular diversity of Peptide Mass Fingerprints (PMFs) of RNases is useful for drug design. This study introduces a methodology that combines QSAR models, 2D-Electrophoresis (2D-E), MALDI-TOF Mass Spectroscopy (MS), BLAST alignment, and Molecular Dynamics (MD) to explore PMFs of RNases. We illustrate this approach by investigating for the first time the PMFs of a new protein of L. infantum. Here we report and compare new versus old predictive models for RNases based on Topological Indices (TIs) of Markov Pseudo-Folding Lattices. These group of indices called Pseudo-folding Lattice 2D-TIs include: Spectral moments pi ( k )(x,y), Mean Electrostatic potentials xi ( k )(x,y), and Entropy measures theta ( k )(x,y). The accuracy of the models (training/cross-validation) was as follows: xi ( k )(x,y)-model (96.0%/91.7%)>pi ( k )(x,y)-model (84.7/83.3) > theta ( k )(x,y)-model (66.0/66.7). We also carried out a 2D-E analysis of biological samples of L. infantum promastigotes focusing on a 2D-E gel spot of one unknown protein with M<20, 100 and pI <7. MASCOT search identified 20 proteins with Mowse score >30, but not one >52 (threshold value), the higher value of 42 was for a probable DNA-directed RNA polymerase. However, we determined experimentally the sequence of more than 140 peptides. We used QSAR models to predict RNase scores for these peptides and BLAST alignment to confirm some results. We also calculated 3D-folding TIs based on MD experiments and compared 2D versus 3D-TIs on molecular phylogenetic analysis of the molecular diversity of these peptides. This combined strategy may be of interest in drug development or target identification.

Complex network spectral moments for ATCUN motif DNA cleavage: first predictive study on proteins of human pathogen parasites.

The development of methods that can predict the metal-mediated biological activity based only on the 3D structure of metal-unbound proteins has become a goal of major importance. This work is dedicated to the amino terminal Cu(II)- and Ni(II)-binding (ATCUN) motifs that participate in the DNA cleavage and have antitumor activity. We have calculated herein, for the first time, the 3D electrostatic spectral moments for 415 different proteins, including 133 potential ATCUN antitumor proteins. Using these parameters as input for Linear Discriminant Analysis, we have found a model that discriminates between ATCUN-DNA cleavage proteins and nonactive proteins with 91.32% Accuracy (379 out of 415 of proteins including both training and external validation series). Finally, the model has predicted for the first time the DNA cleavage function of proteins from the pathogen parasites. We have predicted possible ATCUN-like proteins with a probability higher than 99% in nine parasite families such as Trypanosoma, Plasmodium, Leishmania, or Toxoplasma. The distribution by biological function of the ATCUN proteins predicted has been the following: oxidoreductases 70.5%, signaling proteins 62.5%, lyases 58.2%, membrane proteins 45.5%, ligases 44.4%, hydrolases 41.3%, transferases 39.2%, cell adhesion proteins 34.5%, metal binders 33.5%, translation proteins 25.0%, transporters 16.7%, structural proteins 9.1%, and isomerases 8.2%. The model is implemented at http://miaja.tic.udc.es/Bio-AIMS/ATCUNPred.php.

Unified QSAR approach to antimicrobials. 4. Multi-target QSAR modeling and comparative multi-distance study of the giant components of antiviral drug-drug complex networks.

One limitation of almost all antiviral Quantitative Structure-Activity Relationships (QSAR) models is that they predict the biological activity of drugs against only one species of virus. Consequently, the development of multi-tasking QSAR models (mt-QSAR) to predict drugs activity against different species of virus is of the major vitally important. These mt-QSARs offer also a good opportunity to construct drug-drug Complex Networks (CNs) that can be used to explore large and complex drug-viral species databases. It is known that in very large CNs we can use the Giant Component (GC) as a representative sub-set of nodes (drugs) and but the drug-drug similarity function selected may strongly determines the final network obtained. In the three previous works of the present series we reported mt-QSAR models to predict the antimicrobial activity against different fungi [Gonzalez-Diaz, H.; Prado-Prado, F. J.; Santana, L.; Uriarte, E. Bioorg.Med.Chem.2006, 14, 5973], bacteria [Prado-Prado, F. J.; Gonzalez-Diaz, H.; Santana, L.; Uriarte E. Bioorg.Med.Chem.2007, 15, 897] or parasite species [Prado-Prado, F.J.; González-Díaz, H.; Martinez de la Vega, O.; Ubeira, F.M.; Chou K.C. Bioorg.Med.Chem.2008, 16, 5871]. However, including these works, we do not found any report of mt-QSAR models for antivirals drug, or a comparative study of the different GC extracted from drug-drug CNs based on different similarity functions. In this work, we used Linear Discriminant Analysis (LDA) to fit a mt-QSAR model that classify 600 drugs as active or non-active against the 41 different tested species of virus. The model correctly classifies 143 of 169 active compounds (specificity=84.62%) and 119 of 139 non-active compounds (sensitivity=85.61%) and presents overall training accuracy of 85.1% (262 of 308 cases). Validation of the model was carried out by means of external predicting series, classifying the model 466 of 514, 90.7% of compounds. In order to illustrate the performance of the model in practice, we develop a virtual screening recognizing the model as active 92.7%, 102 of 110 antivirus compounds. These compounds were never use in training or predicting series. Next, we obtained and compared the topology of the CNs and their respective GCs based on Euclidean, Manhattan, Chebychey, Pearson and other similarity measures. The GC of the Manhattan network showed the more interesting features for drug-drug similarity search. We also give the procedure for the construction of Back-Projection Maps for the contribution of each drug sub-structure to the antiviral activity against different species.

Study of peptide fingerprints of parasite proteins and drug-DNA interactions with Markov-Mean-Energy invariants of biopolymer molecular-dynamic lattice networks.

Since the advent of Molecular Dynamics (MD) in biopolymers science with the study by Karplus et al. on protein dynamics, MD has become the by foremost well established, computational technique to investigate structure and function of biomolecules and their respective complexes and interactions. The analysis of the MD trajectories (MDTs) remains, however, the greatest challenge and requires a great deal of insight, experience, and effort. Here, we introduce a new class of invariants for MDTs based on the spatial distribution of Mean-Energy values ξk (L) on a 2D Euclidean space representation of the MDTs. The procedure forces one MD trajectory to fold into a 2D Cartesian coordinates system using a step-by-step procedure driven by simple rules. The ξk (L) values are invariants of a Markov matrix (1 Π), which describes the probabilities of transition between two states in the new 2D space; which is associated to a graph representation of MDTs similar to the lattice networks (LNs) of DNA and protein sequences. We also introduce a new algorithm to perform phylogenetic analysis of peptides based on MDTs instead of the sequence of the polypeptide. In a first experiment, we illustrate this algorithm for 35 peptides present on the Peptide Mass Fingerprint (PMF) of a new protein of Leishmania infantum studied in this work. We report, by the first time, 2D Electrophoresis isolation, MALDI TOF Mass Spectroscopy characterization, and MASCOT search results for this PMF. In a second experiment, we construct the LNs for 422 MDTs obtained in DNA-Drug Docking simulations of the interaction of 57 anticancer furocoumarins with a DNA oligonucleotide. We calculated the respective ξk (L) values for all these LNs and used them as inputs to train a new classifier with Accuracy = 85.44% and 84.91% in training and validation respectively. The new model can be used as scoring function to guide DNA-Drug Docking studies in drug design of new coumarins for PUVA therapy. The new phylogenetics analysis algorithms encode information different from sequence similarity and may be used to analyze MDTs obtained in Docking or modeling experiments for any classes of biopolymers. The work opens new perspective on the analysis and applications of MD in polymer sciences.

Unified QSAR and network-based computational chemistry approach to antimicrobials, part 1: multispecies activity models for antifungals.

There are many pathogen microbial species with very different antimicrobial drugs susceptibility. In this work, we selected pairs of antifungal drugs with similar/dissimilar species predicted-activity profile and represented it as a large network, which may be used to identify drugs with similar mechanism of action. Computational chemistry prediction of the biological activity based on quantitative structure-activity relationships (QSAR) susbtantially increases the potentialities of this kind of networks, avoiding time and resource-consuming experiments. Unfortunately, most QSAR models are unspecific or predict activity against only one species. To solve this problem we developed a multispecies QSAR classification model, in which the outputs were the inputs of the aforementioned network. Overall model classification accuracy was 87.0% (161/185 compounds) in training, 83.4% (50/61) in validation, and 83.7% for 288 additional antifungal compounds used to extend model validation for network construction. The network predicted has 59 nodes (compounds), 648 edges (pairs of compounds with similar activity), low coverage density d = 37.8%, and distribution more close to normal than to exponential. These results are more characteristic of a not-overestimated random network, clustering different drug mechanisms of actions, than of a less useful power law network with few mechanisms (network hubs).

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.

Review of Theoretical Models to Study Natural Products with Antiprotozoal Activity.

In nature, pathogenic parasite species with different susceptibility patterns of antiparasitic drugs abound. In this sense, natural products derived from plants are a potency for drugs with potential antiparasitic activity. Unfortunately, there are many metabolites and studying all of them would be costly in terms of money and resources. To this end, theoretical studies such as QSAR models could be useful. These, for the most part, predict the biological activity of the drugs against a single species of parasite. Consequently, foretell the probability with which a drug is active against many different species with a single QSAR model is an important achievement. This review consists of three parts: the first part is a review of metabolites found in nature that have antiparasitic activity, in particular the antiprotozoal (Leishmania and Trypanosoma); the second part includes a review of theoretical studies looking for a model that predicts the antiprotozoal activity of natural products; the third and final part concerns the study of theoretical models focused on the interaction between drug and receptor, analyzing new metabolites with antiprotozoal activity.

Review of Theoretical Prediction Models for Organic Extract Metabolites, Effect of Drying Temperature on Smooth Muscle Relaxing Activity Induced by Organic Extracts Specially Cecropia Obtusifolia Portal and Web Server Predictors of Drug-Protein Interaction.

Cecropia obtusifolia bertol is medicinal specie used in the treatment of diabetes mellitus and hypertension and it has scientific studies that support the traditional use. However, it is required to understand the influence of drying temperature on the yield and pharmacological activity. Drying rate, extraction efficiency, changes in the UV-Vis spectrum and estimating chlorophylls were stimulated with the increasing temperature. Finally, relaxant activity of vascular smooth muscle is increased by 70ºC and reducing ability by the method of CARF increases with temperature. Analytical studies are required to identify changes in the metabolic content and those that ensure the safety and efficacy for human consumption. In this sense, bioinformatic studies may be helpful. Studies such as QSAR can help us to study these metabolites derived from natural products. MIND-BETS model and NL MIND-BETS model to predict DPIs was introduced using MARCH-INSIDE (MI) software to calculate structural parameters for drugs and enzymes respectively. We firstly revised the state-of-art on the design with review of previous works with hypertension activity based on theoretical studies. A study, evaluating the effect of drying temperature of leaves of C. obtusifolia on the relaxing of vascular smooth muscle, antioxidant activity and the presence of chlorophylls, with a focus on Cecropia metabolites. Last, we carried out QSAR studies using MIND-BEST and NL MIND-BEST web servers in order to understand the essential metabolites structural requirement for binding with receptors for FDA proteins.

Galvez-Markov network transferability indices: review of classic theory and new model for perturbations in metabolic reactions.

Topological Indices (TIs) are numerical parameters useful to carry out Quantitative Structure-Property Relationships (QSPR) analysis and predict the effect of perturbations in many types of Complex Networks. This work, focuses on a very powerful class of TIs called Galvez charge transfer indices. First, we review the classic concept and some applications of these indices. Next, we review the Galvez-Markov TIs of order k (GMk), a recent generalization to these TIs introduced by us. We also reviewed some previous examples of calculation of GMk values for different classes of networks, including metabolic networks. Here, we also demonstrated that Galvez- Markov TIs are useful to predict perturbations and the transferability of biochemical patterns forms metabolic networks of species to others. We report a linear QSPR-Perturbation theory model that predicts more than 300,000 perturbations in metabolic networks with 85 - 99% of good classification in training and validation series.

Model for high-throughput screening of drug immunotoxicity--study of the anti-microbial G1 over peritoneal macrophages using flow cytometry.

Quantitative Structure-Activity (mt-QSAR) techniques may become an important tool for prediction of cytotoxicity and High-throughput Screening (HTS) of drugs to rationalize drug discovery process. In this work, we train and validate by the first time mt-QSAR model using TOPS-MODE approach to calculate drug molecular descriptors and Linear Discriminant Analysis (LDA) function. This model correctly classifies 8258 out of 9000 (Accuracy = 91.76%) multiplexing assay endpoints of 7903 drugs (including both train and validation series). Each endpoint correspond to one out of 1418 assays, 36 molecular and cellular targets, 46 standard type measures, in two possible organisms (human and mouse). After that, we determined experimentally, by the first time, the values of EC50 = 21.58 µg/mL and Cytotoxicity = 23.6% for the anti-microbial/anti-parasite drug G1 over Balb/C mouse peritoneal macrophages using flow cytometry. In addition, the model predicts for G1 only 7 positive endpoints out 1251 cytotoxicity assays (0.56% of probability of cytotoxicity in multiple assays). The results obtained complement the toxicological studies of this important drug. This work adds a new tool to the existing pool of few methods useful for multi-target HTS of ChEMBL and other libraries of compounds towards drug discovery.

Entropy model for multiplex drug-target interaction endpoints of drug immunotoxicity.

Entropy measures are universal parameters useful to codify biologically-relevant information in many systems. In our previous work, (Gonzalez-Diaz, H., et al. Chem. Res. Toxicol. 2003, 16, 1318-1327), we introduced the molecular structure information indices called 3D-Markovian electronic delocalization entropies (3D-MEDNEs) to study the quantitative structure-toxicity relationships (QSTR) of drugs. In a second part, (Cruz-Monteagudo, M. et al. Chem. Res. Toxicol., 2008, 21 (3), 619-632), we extended 3D-MEDNEs to numerically encode toxicologically-relevant information present in Mass Spectra of the serum proteome. These works demonstrated that the idea behind classic drug QSTR models can be extended to solve more general problems in toxicological chemical research. For instance, there are not many reports of multi-target QSTR (mt-QSTR) models useful to predict multiplexed endpoints of drugs in a high number of cytotoxicity assays. In this work, we train and validate for the first time a QSTR model that correctly classifies 8,806 out of 9,001 (Accuracy = 91.1%) multiplexing assay endpoints of 7903 drugs (including both training and validation series). Each endpoint corresponds to one out of 1443 assays, 32 molecular and cellular targets, 46 standard type measures, in two possible organisms (human and mouse). We have also determined experimentally, for the first time, the values of EC50 = 8.21 µg/mL and Cytotoxicity = 26.25 % for the antimicrobial / antiparasitic drug G1 on Balb/C mouse thymic macrophages using flow cytometry. In addition, we have used the new model to predict G1 endpoints in 1,283 assays finding a low average probability of p(1) = 0.50% in 152 cytotoxicity assays. Last, we have used the model to predict average probability of the interaction of G1 with different proteins in macrophages. Interestingly, the Macrophage colony-stimulating factor receptor, the Macrophage colony-stimulating factor 1 receptor, the Macrophage migration inhibitory factor, Macrophage scavenger receptor types I and II, and the Macrophage-stimulating protein receptor, have also very low average predicted probabilities of p(1) = 0.092, 0.038, 0.077, 0.026, 0.2, 0.106, respectively. Both experimental and theoretical results show a moderate thymic macrophage cytotoxicity of G1. The obtained results are significant because they complement the immunotoxicology studies of this important drug.