Juan A. Castillo-Garit

A novel approach to predict aquatic toxicity from molecular structure

The main aim of the study was to develop quantitative structure-activity relationship (QSAR) models for the prediction of aquatic toxicity using atom-based non-stochastic and stochastic linear indices. The used dataset consist of 392 benzene derivatives, separated into training and test sets, for which toxicity data to the ciliate Tetrahymena pyriformis were available. Using multiple linear regression, two statistically significant QSAR models were obtained with non-stochastic (R2=0.791 and s=0.344) and stochastic (R2=0.799 and s=0.343) linear indices. A leave-one-out (LOO) cross-validation procedure was carried out achieving values of q2=0.781 (scv=0.348) and q2=0.786 (scv=0.350), respectively. In addition, a validation through an external test set was performed, which yields significant values of Rpred2 of 0.762 and 0.797. A brief study of the influence of the statistical outliers in QSARs model development was also carried out. Finally, our method was compared with other approaches implemented in the Dragon software achieving better results. The non-stochastic and stochastic linear indices appear to provide an interesting alternative to costly and time-consuming experiments for determining toxicity.

Bond-Based 3D-Chiral Linear Indices: Theory and QSAR Applications to Central Chirality Codification

The recently introduced non‐stochastic and stochastic bond‐based linear indices are been generalized to codify chemical structure information for chiral drugs, making use of a trigonometric 3D‐chirality correction factor. These improved modified descriptors are applied to several well‐known data sets to validate each one of them. Particularly, Cramers steroid data set has become a benchmark for the assessment of novel quantitative structure activity relationship methods. This data set has been used by several researchers using 3D‐QSAR approaches such as Comparative Molecular Field Analysis, Molecular Quantum Similarity Measures, Comparative Molecular Moment Analysis, E‐state, Mapping Property Distributions of Molecular Surfaces, and so on. For that reason, it is selected by us for the sake of comparability. In addition, to evaluate the effectiveness of this novel approach in drug design we model the angiotensin‐converting enzyme inhibitory activity of perindoprilates σ‐stereoisomers combinatorial library, as well as codify information related to a pharmacological property highly dependent on the molecular symmetry of a set of seven pairs of chiral N‐alkylated 3‐(3‐hydroxyphenyl)‐piperidines that bind σ‐receptors. The validation of this method is achieved by comparison with earlier publications applied to the same data sets. The non‐stochastic and stochastic bond‐based 3D‐chiral linear indices appear to provide a very interesting alternative to other more common 3D‐QSAR descriptors.

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.

Tyrosinase Enzyme: 1. An Overview on a Pharmacological Target

The tyrosinase enzyme (EC 1.14.18.1) is an oxidoreductase inside the general enzyme classification and is involved in the oxidation and reduction process in the epidermis. These chemical reactions that the enzyme catalyzes are of principal importance in the melanogenesis process. This process of melanogenesis is related to the melanin formation, a heteropolymer of indolic nature that provides the different tonalities in the skin and helps to the protection from the ultraviolet radiation. However, a pigment overproduction, come up by the action of the tyrosinase, can cause different disorders in the skin related to the hyperpigmentation. Several studies mainly focused on the characteristics of the enzyme have been reported. In this work, an approximation to general aspects related to this enzyme is made. Besides, it is treated the researches that have been published in the part of the biochemical anatomy dealing with diseases associated with this protein (melanogenesis), its active place and its physiological states, the molecular mechanism, the methods carried out to detect the inhibitory activity, and the used substrates.

A review of QSAR studies to discover new drug-like compounds actives against leishmaniasis and trypanosomiasis.

The neglected tropical diseases (NTDs) affect more than one billion people (one-sixth of the worlds population) and occur primarily in undeveloped countries in sub-Saharan Africa, Asia, and Latin America. Available drugs for these diseases are decades old and present an important number of limitations, especially high toxicity and, more recently, the emergence of drug resistance. In the last decade several Quantitative Structure-Activity Relationship (QSAR) studies have been developed in order to identify new organic compounds with activity against the parasites responsible for these diseases, which are reviewed in this paper. The topics summarized in this work are: 1) QSAR studies to identify new organic compounds actives against Chagas disease; 2) Development of QSAR studies to discover new antileishmanial drusg; 3) Computational studies to identify new drug-like compounds against human African trypanosomiasis. Each topic include the general characteristics, epidemiology and chemotherapy of the disease as well as the main QSAR approaches to discovery/identification of new actives compounds for the corresponding neglected disease. The last section is devoted to a new approach know as multi-target QSAR models developed for antiparasitic drugs specifically those actives against trypanosomatid parasites. At present, as a result of these QSAR studies several promising compounds, active against these parasites, are been indentify. However, more efforts will be required in the future to develop more selective (specific) useful drugs.

Bond-based linear indices of the non-stochastic and stochastic edge-adjacency matrix. 1. Theory and modeling of ChemPhys properties of organic molecules

Novel bond-level molecular descriptors are proposed, based on linear maps similar to the ones defined in algebra theory. The kth edge-adjacency matrix (Ek) denotes the matrix of bond linear indices (non-stochastic) with regard to canonical basis set. The kth stochastic edge-adjacency matrix, ESk, is here proposed as a new molecular representation easily calculated from Ek. Then, the kth stochastic bond linear indices are calculated using ESk as operators of linear transformations. In both cases, the bond-type formalism is developed. The kth non-stochastic and stochastic total linear indices are calculated by adding the kth non-stochastic and stochastic bond linear indices, respectively, of all bonds in molecule. First, the new bond-based molecular descriptors (MDs) are tested for suitability, for the QSPRs, by analyzing regressions of novel indices for selected physicochemical properties of octane isomers (first round). General performance of the new descriptors in this QSPR studies is evaluated with regard to the well-known sets of 2D/3D MDs. From the analysis, we can conclude that the non-stochastic and stochastic bond-based linear indices have an overall good modeling capability proving their usefulness in QSPR studies. Later, the novel bond-level MDs are also used for the description and prediction of the boiling point of 28 alkyl-alcohols (second round), and to the modeling of the specific rate constant (log k), partition coefficient (log P), as well as the antibacterial activity of 34 derivatives of 2-furylethylenes (third round). The comparison with other approaches (edge- and vertices-based connectivity indices, total and local spectral moments, and quantum chemical descriptors as well as E-state/biomolecular encounter parameters) exposes a good behavior of our method in this QSPR studies. Finally, the approach described in this study appears to be a very promising structural invariant, useful not only for QSPR studies but also for similarity/diversity analysis and drug discovery protocols.

Comparative study to predict toxic modes of action of phenols from molecular structures.

Quantitative structure–activity relationship models for the prediction of mode of toxic action (MOA) of 221 phenols to the ciliated protozoan Tetrahymena pyriformis using atom-based quadratic indices are reported. The phenols represent a variety of MOAs including polar narcotics, weak acid respiratory uncouplers, pro-electrophiles and soft electrophiles. Linear discriminant analysis (LDA), and four machine learning techniques (ML), namely k-nearest neighbours (k-NN), support vector machine (SVM), classification trees (CTs) and artificial neural networks (ANNs), have been used to develop several models with higher accuracies and predictive capabilities for distinguishing between four MOAs. Most of them showed global accuracy of over 90%, and false alarm rate values were below 2.9% for the training set. Cross-validation, complementary subsets and external test set were performed, with good behaviour in all cases. Our models compare favourably with other previously published models, and in general the models obtained with ML techniques show better results than those developed with linear techniques. We developed unsupervised and supervised consensus, and these results were better than our ML models, the results of rule-based approach and other ensemble models previously published. This investigation highlights the merits of ML-based techniques as an alternative to other more traditional methods for modelling MOA.

Identification In Silico and In Vitro of Novel Trypanosomicidal Drug‐Like Compounds

Atom‐based bilinear indices and linear discriminant analysis are used to discover novel trypanosomicidal compounds. The obtained linear discriminant analysis‐based quantitative structure–activity relationship models, using non‐stochastic and stochastic indices, provide accuracies of 89.02% (85.11%) and 89.60% (88.30%) of the chemicals in the training (test) sets, respectively. Later, both models were applied to the virtual screening of 18 in‐house synthesized compounds to find new pro‐lead antitrypanosomal agents. The in vitro antitrypanosomal activity of this set against epimastigote forms of Trypanosoma cruzi is assayed. Predictions agree with experimental results to a great extent (16/18) of the chemicals. Sixteen compounds show more than 70% of epimastigote inhibition at a concentration 100 μg/mL. In addition, three compounds (CRIS 112, CRIS 140 and CRIS 147) present more than 70% of epimastigote inhibition at a concentration of 10 μg/mL (79.95%, 73.97% and 78.13%, respectively) with low values of cytotoxicity (19.7%, 7.44% and 20.63%, correspondingly).Taking into account all these results, we could say that these three compounds could be optimized in forthcoming works. 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.

Prediction of acute toxicity of phenol derivatives using multiple linear regression approach for Tetrahymena pyriformis contaminant identification in a median-size database

In this article, the modeling of inhibitory grown activity against Tetrahymena pyriformis is described. The 0-2D Dragon descriptors based on structural aspects to gain some knowledge of factors influencing aquatic toxicity are mainly used. Besides, it is done by some enlarged data of phenol derivatives described for the first time and composed of 358 chemicals. It overcomes the previous datasets with about one hundred compounds. Moreover, the results of the model evaluation by the parameters in the training, prediction and validation give adequate results comparable with those of the previous works. The more influential descriptors included in the model are: X3A, MWC02, MWC10 and piPC03 with positive contributions to the dependent variable; and MWC09, piPC02 and TPC with negative contributions. In a next step, a median-size database of nearly 8000 phenolic compounds extracted from ChEMBL was evaluated with the quantitative-structure toxicity relationship (QSTR) model developed providing some clues (SARs) for identification of ecotoxicological compounds. The outcome of this report is very useful to screen chemical databases for finding the compounds responsible of aquatic contamination in the biomarker used in the current work.

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.