Application of linear discriminant analysis in the virtual screening of antichagasic drugs through trypanothione reductase inhibition

Ligand and structure-based discovery of phosphorus-containing compounds as potential metalloproteinase inhibitors

In this study, a methodology is proposed, combining ligand- and structure-based virtual screening tools, for the identification of phosphorus-containing compounds as inhibitors of zinc metalloproteases. First, we use Dragon molecular descriptors to develop a Linear Discriminant Analysis classification model, which is widely validated according to the OECD principles. This model is simple, robust, stable and has good discriminating power. Furthermore, it has a defined applicability domain and it is used for virtual screening of the DrugBank database. Second, docking experiments are carried out on the identified compounds that showed good binding energies to the enzyme thermolysin. Considering the potential toxicity of phosphorus-containing compounds, their toxicological profile is evaluated according to Protox II. Of the five molecules evaluated, two show carcinogenic and mutagenic potential at small LD50, not recommended as drugs, while three of them are classified as non-toxic, and could constitute a starting point for the development of new vasoactive metalloprotease inhibitor drugs. According to molecular dynamics simulation, two of them show stable interactions with the active site maintaining coordination with the metal. A high agreement is evident between QSAR, docking and molecular dynamics results, demonstrating the potentialities of the combination of these tools.

Computational identification of chemical compounds with potential anti-Chagas activity using a classification tree

Chagas disease is endemic to 21 Latin American countries and is a great public health problem in that region. Current chemotherapy remains unsatisfactory; consequently the need to search for new drugs persists. Here we present a new approach to identify novel compounds with potential anti-chagasic action. A large dataset of 584 compounds, obtained from the Drugs for Neglected Diseases initiative, was selected to develop the computational model. Dragon software was used to calculate the molecular descriptors and WEKA software to obtain the classification tree. The best model shows accuracy greater than 93.4% for the training set; the tree was also validated using a 10-fold cross-validation procedure and through a test set, achieving accuracy values over 90.5% and 92.2%, correspondingly. The values of sensitivity and specificity were around 90% in all series; also the false alarm rate values were under 10.5% for all sets. In addition, a simulated ligand-based virtual screening for several compounds recently reported as promising anti-chagasic agents was carried out, yielding good agreement between predictions and experimental results. Finally, the present work constitutes an example of how this rational computer-based method can help reduce the cost and increase the rate in which novel compounds are developed against Chagas disease.

Computational approaches for drug discovery against trypanosomatid-caused diseases

During three decades, only about 20 new drugs have been developed for malaria, tuberculosis and all neglected tropical diseases (NTDs). This critical situation was reached because NTDs represent only 10% of health research investments; however, they comprise about 90% of the global disease burden. Computational simulations applied in virtual screening (VS) strategies are very efficient tools to identify pharmacologically active compounds or new indications for drugs already administered for other diseases. One of the advantages of this approach is the low time-consuming and low-budget first stage, which filters for testing experimentally a group of candidate compounds with high chances of binding to the target and present trypanocidal activity. In this work, we review the most common VS strategies that have been used for the identification of new drugs with special emphasis on those applied to trypanosomiasis and leishmaniasis. Computational simulations based on the selected protein targets or their ligands are explained, including the method selection criteria, examples of successful VS campaigns applied to NTDs, a list of validated molecular targets for drug development and repositioned drugs for trypanosomatid-caused diseases. Thereby, here we present the state-of-the-art of VS and drug repurposing to conclude pointing out the future perspectives in the field.

Probing the opportunities for designing anthelmintic leads by sub-structural topology-based QSAR modelling

A quantitative structure-activity (QSAR) model has been developed for enriched tubulin inhibitors, which were retrieved from sequence similarity searches and applicability domain analysis. Using partial least square (PLS) method and leave-one-out (LOO) validation approach, the model was generated with the correlation statistics of [Formula: see text] and [Formula: see text] of 0.68 and 0.69, respectively. The present study indicates that topological descriptors, viz. BIC, CH_3_C, IC, JX and Kappa_2 correlate well with biological activity. ADME and toxicity (or ADME/T) assessment showed that out of 260 molecules, 255 molecules successfully passed the ADME/T assessment test, wherein the drug-likeness attributes were exhibited. These results showed that topological indices and the colchicine binding domain directly influence the aetiology of helminthic infections. Further, we anticipate that our model can be applied for guiding and designing potential anthelmintic inhibitors.

Bond-based bilinear indices for computational discovery of novel trypanosomicidal drug-like compounds through virtual screening

Two-dimensional bond-based bilinear indices and linear discriminant analysis are used in this report to perform a quantitative structure-activity relationship study to identify new trypanosomicidal compounds. A data set of 440 organic chemicals, 143 with antitrypanosomal activity and 297 having other clinical uses, is used to develop the theoretical models. Two discriminant models, computed using bond-based bilinear indices, are developed and both show accuracies higher than 86% for training and test sets. The stochastic model correctly indentifies nine out of eleven compounds of a set of organic chemicals obtained from our synthetic collaborators. The in vitro antitrypanosomal activity of this set against epimastigote forms of Trypanosoma cruzi is assayed. Both models show a good agreement between theoretical predictions and experimental results. Three compounds showed IC50 values for epimastigote elimination (AE) lower than 50 μM, while for the benznidazole the IC50 = 54.7 μM which was used as reference compound. The value of IC50 for cytotoxicity of these compounds is at least 5 times greater than their value of IC50 for AE. Finally, we can say that, the present algorithm constitutes a step forward in the search for efficient ways of discovering new antitrypanosomal compounds.

Prediction of the maximum absorption wavelength of azobenzene dyes by QSPR tools

The maximum absorption wavelength (λ(max)) of a large data set of 191 azobenzene dyes was predicted by quantitative structure-property relationship (QSPR) tools. The λ(max) was correlated with the 4 molecular descriptors calculated from the structure of the dyes alone. The multiple linear regression method (MLR) and the non-linear radial basis function neural network (RBFNN) method were applied to develop the models. The statistical parameters provided by the MLR model were R(2)=0.893, R(adj)(2)=0.893, q(LOO)(2)=0.884, F=1214.871, RMS=11.6430 for the training set; and R(2)=0.849, R(adj)(2)=0.845, q(ext)(2)=0.846, F=207.812, RMS=14.0919 for the external test set. The RBFNN model gave even improved statistical results: R(2)=0.920, R(adj)(2)=0.919, q(LOO)(2)=0.898, F=1664.074, RMS=9.9215 for the training set, and R(2)=0.895, R(adj)(2)=0.892, q(ext)(2)=0.895, F=314.256, RMS=11.6427 for the external test set. This theoretical method provides a simple, precise and an alternative method to obtain λ(max) of azobenzene dyes.

Dihydroquinazolines as a novel class of Trypanosoma brucei trypanothione reductase inhibitors: discovery, synthesis, and characterization of their binding mode by protein crystallography

Trypanothione reductase (TryR) is a genetically validated drug target in the parasite Trypanosoma brucei , the causative agent of human African trypanosomiasis. Here we report the discovery, synthesis, and development of a novel series of TryR inhibitors based on a 3,4-dihydroquinazoline scaffold. In addition, a high resolution crystal structure of TryR, alone and in complex with substrates and inhibitors from this series, is presented. This represents the first report of a high resolution complex between a noncovalent ligand and this enzyme. Structural studies revealed that upon ligand binding the enzyme undergoes a conformational change to create a new subpocket which is occupied by an aryl group on the ligand. Therefore, the inhibitor, in effect, creates its own small binding pocket within the otherwise large, solvent exposed active site. The TryR-ligand structure was subsequently used to guide the synthesis of inhibitors, including analogues that challenged the induced subpocket. This resulted in the development of inhibitors with improved potency against both TryR and T. brucei parasites in a whole cell assay.

Prediction of enzyme classes from 3D structure: a general model and examples of experimental-theoretic scoring of peptide mass fingerprints of Leishmania proteins

The number of protein and peptide structures included in Protein Data Bank (PDB) and Gen Bank without functional annotation has increased. Consequently, there is a high demand for theoretical models to predict these functions. Here, we trained and validated, with an external set, a Markov Chain Model (MCM) that classifies proteins by their possible mechanism of action according to Enzyme Classification (EC) number. The methodology proposed is essentially new, and enables prediction of all EC classes with a single equation without the need for an equation for each class or nonlinear models with multiple outputs. In addition, the model may be used to predict whether one peptide presents a positive or negative contribution of the activity of the same EC class. The model predicts the first EC number for 106 out of 151 (70.2%) oxidoreductases, 178/178 (100%) transferases, 223/223 (100%) hydrolases, 64/85 (75.3%) lyases, 74/74 (100%) isomerases, and 100/100 (100%) ligases, as well as 745/811 (91.9%) nonenzymes. It is important to underline that this method may help us predict new enzyme proteins or select peptide candidates that improve enzyme activity, which may be of interest for the prediction of new drugs or drug targets. To illustrate the model's application, we report the 2D-Electrophoresis (2DE) isolation from Leishmania infantum as well as MADLI TOF Mass Spectra characterization and theoretical study of the Peptide Mass Fingerprints (PMFs) of a new protein sequence. The theoretical study focused on MASCOT, BLAST alignment, and alignment-free QSAR prediction of the contribution of 29 peptides found in the PMF of the new protein to specific enzyme action. This combined strategy may be used to identify and predict peptides of prokaryote and eukaryote parasites and their hosts as well as other superior organisms, which may be of interest in drug development or target identification.

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.

Discovery of anticonvulsant activity of abietic acid through application of linear discriminant analysis

Linear discriminant analysis was performed to derive discriminant functions based on 2D descriptors and capable of classifying anticonvulsant from non-anticonvulsant compounds. Through application in virtual screening of the discriminant function which performed best in the validation steps, abietic acid was identified as a potential new anticonvulsant agent. The anticonvulsant activity of abietic acid at 30 and 100mg/kg was confirmed in the Maximal Electroshock Test, both orally and intraperitoneally.