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.

Fullerite crystal thermodynamic characteristics and the law of corresponding states

The existence of single-wall carbon nanotubes in organic solvents in the form of clusters is discussed. A theory is developed based on a bundlet model, which enables describing the cluster-size distribution function. Comparison of calculated solubilities with experiments would permit obtaining energetic parameters, characterizing the interaction of a nanotube with its surrounding. Fullerenes and nanotubes are objects whose behaviour in many physical situations is characterized by peculiarities, which show up in that these systems represent the only soluble forms of carbon, what is related to their molecular structures. The fullerene molecule is a virtually uniform closed spherical-spheroidal surface and a nanotube is a smooth cylindrical unit. Both structures give rise to weak interactions between the neighbouring units in a crystal and promote their interaction with solvent molecules. The phenomena have a unified explanation in the bundlet model, in which the free energy of a nanotube in a cluster is combined from two components: a volume one proportional to the number of molecules n in a cluster and a surface one proportional to n1/2. Growth mechanisms of fractal clusters in fullerene solutions are analyzed along with similarity laws, determining the thermodynamic characteristics of fullerite crystals. In accordance with the similarity laws, the dimensionless Debye temperatures theta0 for all crystals belonging to the considered class should be close. Temperatures theta0 are determined by a similarity relation from experimental and estimated data. Fullerite theta0 is twice that for inert-gas crystals because, near the Debye point, the fullerite crystal is orientationally ordered so that its structure is dissimilar to face-centred cubic. A fullerene molecule whose thermal rotation is frozen cannot be considered as a spherically symmetric particle. The fulfilment of the similarity laws, which are valuable for particles with spherically symmetric interaction potential, would hardly be expected.

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.

Atom- and bond-based 2D TOMOCOMD-CARDD approach and ligand-based virtual screening for the drug discovery of new tyrosinase inhibitors

Two-dimensional atom- and bond-based TOMOCOMD-CARDD descriptors and linear discriminant analysis (LDA) are used in this report to perform a quantitative structure-activity relationship (QSAR) study of tyrosinase-inhibitory activity. A database of inhibitors of the enzyme is collected for this study, within 246 highly dissimilar molecules presenting antityrosinase activity. In total, 7 discriminant functions are obtained by using the whole set of atom- and bond-based 2D indices. All the LDA-based QSAR models show accuracies above 90% in the training set and values of the Matthews correlation coefficient (C) varying from 0.85 to 0.90. The external validation set shows globally good classifications between 89% and 91% and C values ranging from 0.75 to 0.81. Finally, QSAR models are used in the selection/identification of the 20 new dicoumarins subset to search for tyrosinase inhibitory activity. Theoretical and experimental results show good correspondence between one another. It is important to remark that most compounds in this series exhibit a more potent inhibitory activity against the mushroom tyrosinase enzyme than the reference compound, Kojic acid (IC(50) = 16.67 muM), resulting in a novel nucleus base (lead) with antityrosinase activity, and this could serve as a starting point for the drug discovery of novel tyrosinase inhibitor lead compounds. ( Journal of Biomolecular Screening 2008:1014-1024).

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 QSAR's 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.

New antitrichomonal drug-like chemicals selected by bond (edge)-based TOMOCOMD-CARDD descriptors

Bond-based quadratic indices, new TOMOCOMD-CARDD molecular descriptors, and linear discriminant analysis (LDA) were used to discover novel lead trichomonacidals. The obtained LDA-based quantitative structure-activity relationships (QSAR) models, using nonstochastic and stochastic indices, were able to classify correctly 87.91% (87.50%) and 89.01% (84.38%) of the chemicals in training (test) sets, respectively. They showed large Matthews correlation coefficients of 0.75 (0.71) and 0.78 (0.65) for the training (test) sets, correspondingly. Later, both models were applied to the virtual screening of 21 chemicals to find new lead antitrichomonal agents. Predictions agreed with experimental results to a great extent because a correct classification for both models of 95.24% (20 of 21) of the chemicals was obtained. Of the 21 compounds that were screened and synthesized, 2 molecules (chemicals G-1, UC-245) showed high to moderate cytocidal activity at the concentration of 10 microg/ml, another 2 compounds (G-0 and CRIS-148) showed high cytocidal activity only at the concentration of 100 microg/ml, and the remaining chemicals (from CRIS-105 to CRIS-153, except CRIS-148) were inactive at these assayed concentrations. Finally, the best candidate, G-1 (cytocidal activity of 100% at 10 microg/ml) was in vivo assayed in ovariectomized Wistar rats achieving promising results as a trichomonacidal drug-like compound.

Prediction of tyrosinase inhibition activity using atom-based bilinear indices

A set of novel atom-based molecular fingerprints is proposed based on a bilinear map similar to that defined in linear algebra. These molecular descriptors (MDs) are proposed as a new means of molecular parametrization easily calculated from 2D molecular information. The nonstochastic and stochastic molecular indices match molecular structure provided by molecular topology by using the kth nonstochastic and stochastic graph-theoretical electronic-density matrices, M(k) and S(k), respectively. Thus, the kth nonstochastic and stochastic bilinear indices are calculated using M(k) and S(k) as matrix operators of bilinear transformations. Chemical information is coded by using different pair combinations of atomic weightings (mass, polarizability, vdW volume, and electronegativity). The results of QSAR studies of tyrosinase inhibitors using the new MDs and linear discriminant analysis (LDA) demonstrate the ability of the bilinear indices in testing biological properties. A database of 246 structurally diverse tyrosinase inhibitors was assembled. An inactive set of 412 drugs with other clinical uses was used; both active and inactive sets were processed by hierarchical and partitional cluster analyses to design training and predicting sets. Twelve LDA-based QSAR models were obtained, the first six using the nonstochastic total and local bilinear indices and the last six with the stochastic MDs. The discriminant models were applied; globally good classifications of 99.58 and 89.96 % were observed for the best nonstochastic and stochastic bilinear indices models in the training set along with high Matthews correlation coefficients (C) of 0.99 and 0.79, respectively, in the learning set. External prediction sets used to validate the models obtained were correctly classified, with accuracies of 100 and 87.78 %, respectively, yielding C values of 1.00 and 0.73. This subset contains 180 active and inactive compounds not considered to fit the models. A simulated virtual screen demonstrated this approach in searching tyrosinase inhibitors from compounds never considered in either training or predicting series. These fitted models permitted the selection of new cycloartane compounds isolated from herbal plants as new tyrosinase inhibitors. A good correspondence between theoretical and experimental inhibitory effects on tyrosinase was observed; compound CA6 (IC(50)=1.32 microM) showed higher activity than the reference compounds kojic acid (IC(50)=16.67 microM) and L-mimosine (IC(50)=3.68 microM).

Estimation of ADME Properties in Drug Discovery: Predicting Caco-2 Cell Permeability Using Atom-Based Stochastic and Non-stochastic Linear Indices

The in vitro determination of the permeability through cultured Caco-2 cells is the most often-used in vitro model for drug absorption. In this report, we use the largest data set of measured P(Caco-2), consisting of 157 structurally diverse compounds. Linear discriminant analysis (LDA) was used to obtain quantitative models that discriminate higher absorption compounds from those with moderate-poorer absorption. The best LDA model has an accuracy of 90.58% and 84.21% for training and test set. The percentage of good correlation, in the virtual screening of 241 drugs with the reported values of the percentage of human intestinal absorption (HIA), was greater than 81%. In addition, multiple linear regression models were developed to predict Caco-2 permeability with determination coefficients of 0.71 and 0.72. Our method compares favorably with other approaches implemented in the Dragon software, as well as other methods from the international literature. These results suggest that the proposed method is a good tool for studying the oral absorption of drug candidates.

Improvement of charge-transfer indices for multifunctional amino acids: application to lysozyme

Valence topological Charge-Transfer (CT) indices are applied to the calculation of pH at the isoelectric point (pI). The model is generalized for molecules with heteroatoms. The ability of the indices for the description of molecular charge distribution is established by comparing them with the pI of 21 amino acids. Linear correlation models are obtained. The CT indices improve multivariable regression equations for pI. The variance decreases by 95%. No superimposition of the corresponding G(k)-J(k) and G(k)(V)-J(k)(V) pairs is observed in most fits, which diminishes the risk of collinearity. The inclusion of heteroatoms in pi-electron system is beneficial for the description of pI, the because of either the role of the additional p orbitals provided by heteroatom or the role of steric factors in pi-electron conjugation. The use of CT/valence CT indices {G(k), J(k), G(k)(V), J(k)(V)} gives limited results for modelling amino acid pI. The inclusion of the number of acidic/basic groups improves the models. The effect is specially noticeable for amino acids with more than two functional groups. The fitting line for 21 amino acids is used to estimate lysozyme pI replacing (1 + Delta n/n(T)) with (M + Delta n)/n(T). The results for lysozyme fragments can estimate the pI of the whole protein with 1-13% error.