Prediction of tyrosinase inhibition activity using atom-based bilinear indices

Computational discovery of novel anthelmintic natural compounds from Agave Brittoniana trel. Spp. Brachypus

Helminth infections are a medical problem in the world nowadays. This report used bond-based 2D quadratic indices, a bond-level QuBiLs-MAS molecular descriptor family, and Linear Discriminant Analysis (LDA) to obtain a quantitative linear model that discriminates between anthelmintic and non-anthelmintic drug-like organic-compounds. The model obtained correctly classified 87.46% and 81.82% of the training and external data sets, respectively. The developed model was used in a virtual screening to predict the biological activity of all chemicals (19) previously obtained and chemically characterized by some authors of this report from Agave brittoniana Trel. spp. Brachypus. The model identified several metabolites (12) as possible anthelmintics, and a group of 5 novel natural products was tested in an in vitro assay against Fasciola hepatica (100% effectivity at 500 µg/mL). Finally, the two best hits were evaluated in vivo in bald/c mice and the same helminth parasite using a 25 mg/kg dose. Compound 8 (Karatavinoside A) showed an efficacy of 92.2% in vivo. It is important to remark that this natural compound exhibits similar-to-superior activity as triclabendazole, the best human fasciolicide available in the market against Fasciola hepatica, resulting in a novel lead scaffold with anti-helminthic activity. Keywords: TOMOCOMD-CARDD Software; QuBiLs-MAS, nonstochastic and stochastic bond-based quadratic indices; LDA-based QSAR model; Computational Screening, Anthelmintic Agent; Agave brittoniana Trel. spp. Brachypus, Fasciola hepatica.

The Chemical Structure and Bioactivity of Cycloartane-type Compounds

:

For decades now, compounds in the cycloartane-type series have been shown to have versatile pharmacological activities. However, no extensive review has been written to summarize these health-beneficial activities. Therefore, the purpose of this paper is to systematically highlight the biological activities of these compounds, including their antitumor and anti-osteoporosis effects, their effects on receptors, cytokine release, and chronic renal failure, as well as their tyrosinase inhibitory, anticomplement, anti-parasite, anti-HIV, and antituberculosis activities. In this review, we have summarized the structures of over 200 compounds based on their characteristics and described their structureactivity relationships (SARs), and potential mechanisms of action.

Fragment-based in silico modeling of multi-target inhibitors against breast cancer-related proteins

Breast cancer is the most frequent cancer reported in women, being responsible for hundreds of thousands of deaths. Chemotherapy has proven to be effective against this malignant neoplasm depending on different biological factors such as the histopathology, grade, and stage, among others. However, breast cancer cells have become resistant to current chemotherapeutic regimens, urging the discovery of new anti-breast cancer drugs. Computational approaches have the potential to offer promising alternatives to accelerate the search for potent and versatile anti-breast cancer agents. In the present work, we introduce the first multitasking (mtk) computational model devoted to the in silico fragment-based design of new molecules with high inhibitory activity against 19 different proteins involved in breast cancer. The mtk-computational model was created from a dataset formed by 24,285 cases, and it exhibited accuracy around 93% in both training and prediction (test) sets. Several molecular fragments were extracted from the molecules present in the dataset, and their quantitative contributions to the inhibitory activities against all the proteins under study were calculated. The combined use of the fragment contributions and the physicochemical interpretations of the different molecular descriptors in the mtk-computational model allowed the design of eight new molecular entities not reported in our dataset. These molecules were predicted as potent multi-target inhibitors against all the proteins, and they exhibited a desirable druglikeness according to the Lipinski's rule of five and its variants.

Antitumor Activity of HM781-36B, alone or in Combination with Chemotherapeutic Agents, in Colorectal Cancer Cells

Purpose

HM781-36B is a novel and irreversible pan-human epidermal growth factor receptor (HER) inhibitor with TEC cytoplasmic kinase inhibition. The aim of this study is to evaluate the antitumor activity and mechanism of action for HM781-36B in CRC cell lines.

Materials and Methods

The CRC cell lines were exposed to HM781-36B and/or oxaliplatin (L-OHP), 5-fluorouracil (5-FU), SN-38. The cell viability was examined by Cell Titer-Glo luminescent cell viability assay kit. Change in the cell cycle and protein expression was determined by flow cytometry and immunoblot analysis, respectively. Synergism between 2 drugs was evaluated by the combination index.

Results

The addition of HM781-36B induced potent growth inhibition in both DiFi cells with EGFR overexpression and SNU-175 cells (IC₅₀ = 0.003 and 0.005 μM, respectively). Furthermore, HM781-36B induced G1 arrest of the cell cycle and apoptosis, and reduced the levels of HER family and downstream signaling molecules, pERK and pAKT, as well as nonreceptor/cytoplasmic tyrosine kinase, BMX. The combination of HM781-36B with 5-FU, L-OHP, or SN-38 showed an additive or synergistic effect in most CRC cells.

Conclusion

These findings suggest the potential roles of HM781-36B as the treatment for EGFR-overexpressing colon cancer, singly or in combination with chemotherapeutic agents. The role of BMX expression as a marker of response to HM781-36B should be further explored.

Discovery of highly potent tyrosinase inhibitor, T1, with significant anti-melanogenesis ability by zebrafish in vivo assay and computational molecular modeling

Tyrosinase is involved in melanin biosynthesis and the abnormal accumulation of melanin pigments leading to hyperpigmentation disorders that can be treated with depigmenting agents. A natural product T1, bis(4-hydroxybenzyl)sulfide, isolated from the Chinese herbal plant, Gastrodia elata, is a strong competitive inhibitor against mushroom tyrosinase (IC₅₀ = 0.53 μM, Ki = 58 ± 6 nM), outperforms than kojic acid. The cell viability and melanin quantification assay demonstrate that 50 μM of T1 apparently attenuates 20% melanin content of human normal melanocytes without significant cell toxicity. Moreover, the zebrafish in vivo assay reveals that T1 effectively reduces melanogenesis with no adverse side effects. The acute oral toxicity study evidently confirms that T1 molecule is free of discernable cytotoxicity in mice. Furthermore, the molecular modeling demonstrates that the sulfur atom of T1 coordinating with the copper ions in the active site of tyrosinase is essential for mushroom tyrosinase inhibition and the ability of diminishing the human melanin synthesis. These results evident that T1 isolated from Gastrodia elata is a promising candidate in developing pharmacological and cosmetic agents of great potency in skin-whitening.

On the interpretation of tyrosinase inhibition kinetics

Abstract Tyrosinases (monophenol monooxygenases, EC 1.14.18.1) utilize a reaction mechanism that involves two intertwined catalytic cycles and at least three different ligand-binding enzyme forms. Therefore a variety of different inhibition types may arise in inhibition experiments, depending on the binding mode of the compound studied. Here we discuss a steady-state equation that describes inhibition of the diphenolase cycle of tyrosinase catalysis in a general way. In addition, we employ numerical simulations to explore the kinetic outcome of various binding schemes. As the full equation is far too complex to be applicable for data evaluation by curve fitting, we propose to use the general modifier scheme of Botts-Morales for fitting and demonstrate that especially the value of parameter α of the equation allows conclusions about the binding mode of the inhibitor. The approach is exemplified by selected data describing the inhibition of human tyrosinase by typical inhibitors.

Neonicotinoid insecticides: oxidative stress in planta and metallo-oxidase inhibition

Neonicotinoids not only control insect pests but also sometimes independently alter plant growth and response to stress. We find that imidacloprid, thiacloprid, acetamiprid, thiamethoxam, and clothianidin but not nitenpyram and dinotefuran induce foliar lesions and peroxidative damage in soybean ( Glycine max ) seedlings assayed with the 3,3'-diaminobenzidine stain. The chloropyridinyl-carboxylic acid (COOH) but not the -carboxaldehyde (CHO) metabolites induce peroxidative damage but in a different pattern. Surprisingly, the chlorothiazolyl -CHO and -COOH metabolites induce chlorosis but no clear superimposable peroxidative damage or cell death. Four metallo-oxidases known to modulate reactive oxygen species were not sensitive in vitro to the parent neonicotinoid itself but were to several CHO and COOH metabolites and related compounds, with a sensitivity order of CHO > COOH and tyrosinase > xanthine oxidase and aldehyde oxidase > catalase. Although metallo-oxidase inhibition does not correlate overall with lesion formation, it may play an as yet unknown role in plant response to neonicotinoids.

Non-linear models based on simple topological indices to identify RNase III protein members

Alignment-free classifiers are especially useful in the functional classification of protein classes with variable homology and different domain structures. Thus, the Topological Indices to BioPolymers (TI2BioP) methodology (Agüero-Chapin et al., 2010) inspired in both the TOPS-MODE and the MARCH-INSIDE methodologies allows the calculation of simple topological indices (TIs) as alignment-free classifiers. These indices were derived from the clustering of the amino acids into four classes of hydrophobicity and polarity revealing higher sequence-order information beyond the amino acid composition level. The predictability power of such TIs was evaluated for the first time on the RNase III family, due to the high diversity of its members (primary sequence and domain organization). Three non-linear models were developed for RNase III class prediction: Decision Tree Model (DTM), Artificial Neural Networks (ANN)-model and Hidden Markov Model (HMM). The first two are alignment-free approaches, using TIs as input predictors. Their performances were compared with a non-classical HMM, modified according to our amino acid clustering strategy. The alignment-free models showed similar performances on the training and the test sets reaching values above 90% in the overall classification. The non-classical HMM showed the highest rate in the classification with values above 95% in training and 100% in test. Although the higher accuracy of the HMM, the DTM showed simplicity for the RNase III classification with low computational cost. Such simplicity was evaluated in respect to HMM and ANN models for the functional annotation of a new bacterial RNase III class member, isolated and annotated by our group.

Vanilloid derivatives as tyrosinase inhibitors driven by virtual screening-based QSAR models

A number of vanilloids have been tested as tyrosinase inhibitors using Ligand-Based Virtual Screening (LBVS) driven by QSAR (Quantitative Structure-Activity Relationship) models as the multi-agent classification system. A total of 81 models were used to screen this family. Then, a preliminary cluster analysis of the selected chemicals was carried out based on their bioactivity to detect possible similar substructural features among these compounds and the active database used in the QSAR model construction. The compounds identified were tested in vitro to corroborate the results obtained in silico. Among them, two chemicals, isovanillin (K(M) (app) = 1.08 mM) near to kojic acid (reference drug) in one cluster and isovanillyl alcohol (K(M) (app) = 0.88 mM) at the same distance as hydroquinone (reference drug) in another cluster showed inhibitory activity against tyrosinase. The algorithm proposed here could result in a suitable approach for faster and more effective identification of hit and/or lead compounds with tyrosinase inhibitory activity, helping to shorten the long pipeline in the research of novel depigmenting agents to treat skin disorders.

Bond-based 2D quadratic fingerprints in QSAR studies: virtual and in vitro tyrosinase inhibitory activity elucidation

In this report, we show the results of quantitative structure-activity relationship (QSAR) studies of tyrosinase inhibitory activity, by using the bond-based quadratic indices as molecular descriptors (MDs) and linear discriminant analysis (LDA), to generate discriminant functions to predict the anti-tyrosinase activity. The best two models [Eqs (6) and (12)] out of the total 12 QSAR models developed here show accuracies of 93.51% and 91.21%, as well as high Matthews correlation coefficients (C) of 0.86 and 0.82, respectively, in the training set. The validation external series depicts values of 90.00% and 89.44% for these best two equations (6) and (12), respectively. Afterwards, a second external prediction data are used to perform a virtual screening of compounds reported in the literature as active (tyrosinase inhibitors). In a final step, a series of lignans is analysed using the in silico-developed models, and in vitro corroboration of the activity is carried out. An issue of great importance to remark here is that all compounds present greater inhibition values than Kojic acid (standard tyrosinase inhibitor: IC₅₀ = 16.67 μm). The current obtained results could be used as a framework to increase the speed, in the biosilico discovery of leads for the treatment of skin disorders.

Novel coumarin-based tyrosinase inhibitors discovered by OECD principles-validated QSAR approach from an enlarged, balanced database

The present work is devoted to the development and application of a multi-agent Quantitative Structure-Activity Relationship (QSAR) classification system for tyrosinase inhibitor identification, in which the individual QSAR outputs are the inputs of a fusion approach based on the voting mechanism. The individual models are based on TOMOCOMD-CARDD (TOpological Molecular COMputational Design-Computer Aided Rational Drug Design) atom-based bilinear descriptors and Linear Discriminant Analysis (LDA) on a novel enlarged, balanced database of 1,429 compounds within 701 greatly dissimilar molecules presenting anti-tyrosinase activity. A total of 21 adequate models are obtained taking into account the requirements of the Organization for Economic Cooperation and Development (OECD) principles for QSAR validation and present global accuracies (Q) above 84.50 and 79.27% in the training and test sets, respectively. The resulted fusion system is used for the in silico identification of synthesized coumarin derivatives as novel tyrosinase inhibitors. The 7-hydroxycoumarin (compound C07) shows potent activity for the inhibition of monophenolase activity of mushroom tyrosinase giving a value of inhibition percentage close to 100% in vitro assays, by means of spectrophotometric analysis. The current report could help to shed some clues in the identification of new chemicals that inhibit tyrosinase enzyme, for entering in the pipeline of drug discovery development.

Using topological indices to predict anti-Alzheimer and anti-parasitic GSK-3 inhibitors by multi-target QSAR in silico screening

Plasmodium falciparum, Leishmania, Trypanosomes, are the causers of diseases such as malaria, leishmaniasis and African trypanosomiasis that nowadays are the most serious parasitic health problems worldwide. The great number of deaths and the few drugs available against these parasites, make necessary the search for new drugs. Some of these antiparasitic drugs also are GSK-3 inhibitors. GSKI-3 are candidates to develop drugs for the treatment of Alzheimer's disease. In this work topological descriptors for a large series of 3,370 active/non-active compounds were initially calculated with the ModesLab software. Linear Discriminant Analysis was used to fit the classification function and it predicts heterogeneous series of compounds like paullones, indirubins, meridians, etc. This study thus provided a general evaluation of these types of molecules.

Computational discovery of novel trypanosomicidal drug-like chemicals by using bond-based non-stochastic and stochastic quadratic maps and linear discriminant analysis

Herein we present results of a quantitative structure-activity relationship (QSAR) studies to classify and design, in a rational way, new antitrypanosomal compounds by using non-stochastic and stochastic bond-based quadratic indices. A data set of 440 organic chemicals, 143 with antitrypanosomal activity and 297 having other clinical uses, is used to develop QSAR models based on linear discriminant analysis (LDA). Non-stochastic model correctly classifies more than 93% and 95% of chemicals in both training and external prediction groups, respectively. On the other hand, the stochastic model shows an accuracy of about the 87% for both series. As an experiment of virtual lead generation, the present approach is finally satisfactorily applied to the virtual evaluation of 9 already synthesized in house compounds. The in vitro antitrypanosomal activity of this series against epimastigote forms of Trypanosoma cruzi is assayed. The model is able to predict correctly the behaviour for the majority of these compounds. Four compounds (FER16, FER32, FER33 and FER 132) showed more than 70% of epimastigote inhibition at a concentration of 100 microg/mL (86.74%, 78.12%, 88.85% and 72.10%, respectively) and two of these chemicals, FER16 (78.22% of AE) and FER33 (81.31% of AE), also showed good activity at a concentration of 10 microg/mL. At the same concentration, compound FER16 showed lower value of cytotoxicity (15.44%), and compound FER33 showed very low value of 1.37%. Taking into account all these results, we can say that these three compounds can be optimized in forthcoming works, but we consider that compound FER33 is the best candidate. Even though none of them resulted more active than Nifurtimox, the current results constitute a step forward in the search for efficient ways to discover new lead antitrypanosomals.

Variations in IC(50) values with purity of mushroom tyrosinase

The effects of various inhibitors on crude, commercial and partially purified commercial mushroom tyrosinase were examined by comparing IC(50) values. Kojic acid, salicylhydroxamic acid, tropolone, methimazole, and ammonium tetrathiomolybdate had relatively similar IC(50) values for the crude, commercial and partially purified enzyme. 4-Hexylresorcinol seemed to have a somewhat higher IC(50) value using crude extracts, compared to commercial or purified tyrosinase. Some inhibitors (NaCl, esculetin, biphenol, phloridzin) showed variations in IC(50) values between the enzyme samples. In contrast, hydroquinone, lysozyme, Zn(2+), and anisaldehyde showed little or no inhibition in concentration ranges reported to be effective inhibitors. Organic solvents (DMSO and ethanol) had IC(50) values that were similar for some of the tyrosinase samples. Depending of the source of tyrosinase and choice of inhibitor, variations in IC(50) values were observed.

Nucleotide's bilinear indices: novel bio-macromolecular descriptors for bioinformatics studies of nucleic acids. I. Prediction of paromomycin's affinity constant with HIV-1 Psi-RNA packaging region

A new set of nucleotide-based bio-macromolecular descriptors are presented. This novel approach to bio-macromolecular design from a linear algebra point of view is relevant to nucleic acids quantitative structure-activity relationship (QSAR) studies. These bio-macromolecular indices are based on the calculus of bilinear maps on Re(n)[b(mk)(x (m),y (m)):Re(n) x Re(n)-->Re] in canonical basis. Nucleic acid's bilinear indices are calculated from kth power of non-stochastic and stochastic nucleotide's graph-theoretic electronic-contact matrices, M(m)(k) and (s)M(m)(k), respectively. That is to say, the kth non-stochastic and stochastic nucleic acid's bilinear indices are calculated using M(m)(k) and (s)M(m)(k) as matrix operators of bilinear transformations. Moreover, biochemical information is codified by using different pair combinations of nucleotide-base properties as weightings (experimental molar absorption coefficient epsilon(260) at 260 nm and pH=7.0, first (Delta E(1)) and second (Delta E(2)) single excitation energies in eV, and first (f(1)) and second (f(2)) oscillator strength values (of the first singlet excitation energies) of the nucleotide DNA-RNA bases. As example of this approach, an interaction study of the antibiotic paromomycin with the packaging region of the HIV-1 Psi-RNA have been performed and it have been obtained several linear models in order to predict the interaction strength. The best linear model obtained by using non-stochastic bilinear indices explains about 91% of the variance of the experimental Log K (R=0.95 and s=0.08 x 10(-4)M(-1)) as long as the best stochastic bilinear indices-based equation account for 93% of the Log K variance (R=0.97 and s=0.07 x 10(-4)M(-1)). The leave-one-out (LOO) press statistics, evidenced high predictive ability of both models (q(2)=0.86 and s(cv)=0.09 x 10(-4)M(-1) for non-stochastic and q(2)=0.91 and s(cv)=0.08 x 10(-4)M(-1) for stochastic bilinear indices). The nucleic acid's bilinear indices-based models compared favorably with other nucleic acid's indices-based approaches reported nowadays. These models also permit the interpretation of the driving forces of the interaction process. In this sense, developed equations involve short-reaching (k<or=3), middle-reaching (4<k<9), and far-reaching (k=10 or greater) nucleotide's bilinear indices. This situation points to electronic and topologic nucleotide's backbone interactions control of the stability profile of paromomycin-RNA complexes. Consequently, the present approach represents a novel and rather promising way to theoretical-biology studies.

Comparative study of topological indices of macro/supramolecular RNA complex networks

RNA function annotation is often based on alignment to a previously studied template. In contrast to the study of proteins, there are not many alignment-free methods to predict RNA functions if alignment fails. The use of topological indices (TIs) of RNA complex networks (CNs) to find quantitative structure-activity relationships (QSAR) may be an alternative to incorporate secondary structure or sequence-to-sequence similarity. Here, we introduce new QSAR-like techniques using RNA macromolecular CNs (mmCNs), where nodes are nucleotides, or RNA supramolecular CNs (smCNs), where nodes are RNA sequences. We studied a data set of 198 sequences including 18S-rRNAs (important phylogenetic molecular biomarkers). We constructed three types of RNA mmCNs: sequence-linear (SL), Cartesian-lattice (CL), and sequence-folding CNs (SF-CNs) and two smCNs: sequence-sequence disagreement CN (SSD) and sequence-sequence similarity (SSS-smCN). We reported the first comparative QSAR study with all these CIs and CNs, which includes: (i) spectral moments ( ( i )micro d ( w)) of SL-mmCNs (accuracy = 75.3%), (ii) electrostatic CIs (xi d ) of CL-mmCNs (>90%), (iii) thermodynamic parameters (Delta G, Delta H, Delta S, and T m) of SF-mmCNs (64.7%), (iv) disagreement-distribution moments ( M k ) of the SSD-smCN (79.3%), and (v) node centralities of the SSD-smCN (78.0%). Furthermore, we reported the experimental isolation of a new RNA sequence from Psidum guajava leaf tissue and its QSAR and BLAST prediction to illustrate the practical use of these methods. We also investigated the use of these CNs to explore rRNA diversity on bacteria, plants, and parasites from the Dactylogyrus genus. The HPL-mmCNs model was the best of all found. All the CNs and TIs, except SF-mmCNs, were introduced here by the first time for the QSAR study of RNA, which allowed a comparative study for RNA classification.

HP-Lattice QSAR for dynein proteins: experimental proteomics (2D-electrophoresis, mass spectrometry) and theoretic study of a Leishmania infantum sequence

The toxicity and inefficacy of actual organic drugs against Leishmaniosis justify research projects to find new molecular targets in Leishmania species including Leishmania infantum (L. infantum) and Leishmaniamajor (L. major), both important pathogens. In this sense, quantitative structure-activity relationship (QSAR) methods, which are very useful in Bioorganic and Medicinal Chemistry to discover small-sized drugs, may help to identify not only new drugs but also new drug targets, if we apply them to proteins. Dyneins are important proteins of these parasites governing fundamental processes such as cilia and flagella motion, nuclear migration, organization of the mitotic splinde, and chromosome separation during mitosis. However, despite the interest for them as potential drug targets, so far there has been no report whatsoever on dyneins with QSAR techniques. To the best of our knowledge, we report here the first QSAR for dynein proteins. We used as input the Spectral Moments of a Markov matrix associated to the HP-Lattice Network of the protein sequence. The data contain 411 protein sequences of different species selected by ClustalX to develop a QSAR that correctly discriminates on average between 92.75% and 92.51% of dyneins and other proteins in four different train and cross-validation datasets. We also report a combined experimental and theoretic study of a new dynein sequence in order to illustrate the utility of the model to search for potential drug targets with a practical example. First, we carried out a 2D-electrophoresis analysis of L. infantum biological samples. Next, we excised from 2D-E gels one spot of interest belonging to an unknown protein or protein fragment in the region M<20,200 and pI<4. We used MASCOT search engine to find proteins in the L. major data base with the highest similarity score to the MS of the protein isolated from L. infantum. We used the QSAR model to predict the new sequence as dynein with probability of 99.99% without relying upon alignment. In order to confirm the previous function annotation we predicted the sequences as dynein with BLAST and the omniBLAST tools (96% alignment similarity to dyneins of other species). Using this combined strategy, we have successfully identified L. infantum protein containing dynein heavy chain, and illustrated the potential use of the QSAR model as a complement to alignment tools.

Enzyme, protein, carbohydrate, and phenolic contaminants in commercial tyrosinase preparations: potential problems affecting tyrosinase activity and inhibition studies

Commercial mushroom tyrosinase contains other proteins, enzymes, carbohydrates, and phenolic material besides tyrosinase. Carbohydrate and phenolic material comprise a large percentage of the powder resuspensions derived from Agaricus bisporus. Enzyme assays identified the presence of tyrosinase, laccase, beta-glucosidase, beta-galactosidase, beta-xylosidase, cellulase, chitinase, xylanase, and mannanase in the commercial tyrosinase. Protein sequencing indicated the presence of tyrosinase, a lectin, and a putative mannanase as well as 10 unidentified protein/peptides in the commercial tyrosinase preparations. Characteristics of tyrosinase isoforms were similar in two different commercial tyrosinase sources. Inhibition studies indicated that I 50 values for some tyrosinase inhibitors were different when the crude powder was compared to a partially purified tyrosinase. The presence of these contaminants has the potential to affect studies using commercial tyrosinase.

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.

Dragon method for finding novel tyrosinase inhibitors: Biosilico identification and experimental in vitro assays

QSAR (quantitative structure-activity relationship) studies of tyrosinase inhibitors employing Dragon descriptors and linear discriminant analysis (LDA) are presented here. A data set of 653 compounds, 245 with tyrosinase inhibitory activity and 408 having other clinical uses were used. The active data set was processed by k-means cluster analysis in order to design training and prediction series. Seven LDA-based QSAR models were obtained. The discriminant functions applied showed a globally good classification of 99.79% for the best model Class=-96.067+1.988 x 10(2)X0Av +9 1.907 BIC3 + 6.853 CIC1 in the training set. External validation processes to assess the robustness and predictive power of the obtained model were carried out. This external prediction set had an accuracy of 99.44%. After that, the developed models were used in ligand-based virtual screening of tyrosinase inhibitors from the literature and never considered in either training or predicting series. In this case, all screened chemicals were correctly classified by the LDA-based QSAR models. As a final point, these fitted models were used in the screening of new bipiperidine series as new tyrosinase inhibitors. These methods are an adequate alternative to the process of selection/identification of new bioactive compounds. The biosilico assays and in vitro results of inhibitory activity on mushroom tyrosinase showed good correspondence. It is important to stand out that compound BP4 (IC(50)=1.72 microM) showed higher activity in the inhibition against the enzyme than reference compound kojic acid (IC(50)=16.67 microM) and l-mimosine (IC(50)=3.68 microM). These results support the role of biosilico algorithm for the identification of new tyrosinase inhibitor compounds.