Bond-based 2D TOMOCOMD-CARDD approach for drug discovery: aiding decision-making in 'in silico' selection of new lead tyrosinase inhibitors

Computational studies of tyrosinase inhibitors

Computational studies have significantly advanced the understanding of tyrosinase (TYR) function, mechanism, and inhibition, accelerating the development of more effective and selective inhibitors. This chapter provides an overview of in silico studies on TYR inhibitors, emphasizing key inhibitory chemotypes and the main residues involved in ligand-target interactions. The chapter discusses tools applied in the context of TYR inhibitor development, e.g., structure-based virtual screening, molecular docking, artificial intelligence, and machine learning algorithms.

Serendipitous discovery of short peptides from natural products as tyrosinase inhibitors

Tyrosinase, which is the crucial copper-containing enzyme involved in melanin synthesis, is strongly associated with hyperpigmentation disorders, cancer, and neurodegenerative disease; thus, it has attracted considerable interest in the fields of medicine and cosmetics. The known tyrosinase inhibitors show numerous adverse side effects, and there is a lack of safety regulations governing their use. As a result, there is a need to develop novel inhibitors with no toxicity and long-term stability. In this study, we use molecular docking and pharmacophore modeling to construct a reasonable and reliable pharmacophore model, called Hypo 1, that could be used for identifying potent natural products with crucial complementary functional groups for mushroom tyrosinase inhibition. It was observed that, out of 47,263 natural compounds, A5 structurally resembles a dipeptide (WY) and natural compound B16 is the equivalent of a tripeptide (KFY), revealing that the C-terminus tyrosine residues play a key role in tyrosinase inhibition. Tripeptides RCY and CRY, which show high tyrosinase inhibitory potency, revealed a positional and functional preference for the cysteine residue at the N-terminus of the tripeptides, essentially determining the capacity of tyrosinase inhibition. CRY and RCY used the thiol group of cysteine residues to coordinate with the Cu ions in the active site of tyrosinase and showed reduced tyrosinase activity. We discovered the novel tripeptide CRY that shows the most striking inhibitory potency against mushroom tyrosinase (IC50 = 6.16 μM); this tripeptide is more potent than the known oligopeptides and comparable with kojic acid-tripeptides. Our study provides an insight into the structural and functional roles of key amino acids of tripeptides derived from the natural compound B16, and the results are expected to be useful for the development of tyrosinase inhibitors.

Rational design of new agrochemical fungicides using substructural descriptors

BACKGROUND

The increasing resistance of several phytopathogenic fungal species to existing agrochemical fungicides has alarmed the worldwide scientific community. In an attempt to overcome this problem, a discriminant model based on substructural descriptors was developed from a heterogeneous database of compounds for the design of, search for and prediction of agrochemical fungicides.

RESULTS

The discriminant model classifies correctly 81.95% of the fungicides and 81.54% of the inactive compounds in the training series, with an accuracy of 81.72%. In the prediction series, the percentage of correct classification was 80.59 and 85.56% for fungicides and inactive compounds respectively, with an accuracy of 83.44%. Some fragments were extracted and their contributions were calculated. From the fragments that were determined to make positive contributions to the fungicidal activity, new molecules such as pyrrole derivatives were designed and the probabilities of their being fungicides were calculated. These molecules were correctly classified as potential fungicides.

CONCLUSION

The discriminant model based on substructural descriptors provides a promising methodology for the development of molecular patterns to be used in the design of, search for and prediction of agrochemical fungicides of wide spectrum. This constitutes an alternative for the discovery of compounds that are able to decrease crop losses caused by phytopathogenic fungal species.

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.

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.

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.

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).