Toward structure-multiple activity relationships (SMARts) using computational approaches: A polypharmacological perspective.

In the current era of biological big data, which are rapidly populating the biological chemical space, in silico polypharmacology drug design approaches help to decode structure-multiple activity relationships (SMARts). Current computational methods can predict or categorize multiple properties simultaneously, which aids the generation, identification, curation, prioritization, optimization, and repurposing of molecules. Computational methods have generated opportunities and challenges in medicinal chemistry, pharmacology, food chemistry, toxicology, bioinformatics, and chemoinformatics. It is anticipated that computer-guided SMARts could contribute to the full automatization of drug design and drug repurposing campaigns, facilitating the prediction of new biological targets, side and off-target effects, and drug-drug interactions.

Consensus Pharmacophore Strategy For Identifying Novel SARS-Cov-2 Mpro Inhibitors from Large Chemical Libraries.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) main Protease (Mpro) is an enzyme that cleaves viral polyproteins translated from the viral genome and is critical for viral replication. Mpro is a target for anti-SARS-CoV-2 drug development, and multiple Mpro crystals complexed with competitive inhibitors have been reported. In this study, we aimed to develop an Mpro consensus pharmacophore as a tool to expand the search for inhibitors. We generated a consensus model by aligning and summarizing pharmacophoric points from 152 bioactive conformers of SARS-CoV-2 Mpro inhibitors. Validation against a library of conformers from a subset of ligands showed that our model retrieved poses that reproduced the crystal-binding mode in 77% of the cases. Using models derived from a consensus pharmacophore, we screened >340 million compounds. Pharmacophore-matching and chemoinformatics analyses identified new potential Mpro inhibitors. The candidate compounds were chemically dissimilar to the reference set, and among them, demonstrating the relevance of our model. We evaluated the effect of 16 candidates on Mpro enzymatic activity finding that seven have inhibitory activity. Three compounds (1, 4, and 5) had IC50 values in the midmicromolar range. The Mpro consensus pharmacophore reported herein can be used to identify compounds with improved activity and novel chemical scaffolds against Mpro. The method developed for its generation is provided as an open-access code (https://github.com/AngelRuizMoreno/ConcensusPharmacophore) and can be applied to other pharmacological targets.

Identification of SARS-CoV-2 Main Protease Inhibitors Using Chemical Similarity Analysis Combined with Machine Learning.

SARS-CoV-2 Main Protease (Mpro) is an enzyme that cleaves viral polyproteins translated from the viral genome, which is critical for viral replication. Mpro is a target for anti-SARS-CoV-2 drug development. Herein, we performed a large-scale virtual screening by comparing multiple structural descriptors of reference molecules with reported anti-coronavirus activity against a library with >17 million compounds. Further filtering, performed by applying two machine learning algorithms, identified eighteen computational hits as anti-SARS-CoV-2 compounds with high structural diversity and drug-like properties. The activities of twelve compounds on Mpro's enzymatic activity were evaluated by fluorescence resonance energy transfer (FRET) assays. Compound 13 (ZINC13878776) significantly inhibited SARS-CoV-2 Mpro activity and was employed as a reference for an experimentally hit expansion. The structural analogues 13a (ZINC4248385), 13b (ZNC13523222), and 13c (ZINC4248365) were tested as Mpro inhibitors, reducing the enzymatic activity of recombinant Mpro with potency as follows: 13c > 13 > 13b > 13a. Then, their anti-SARS-CoV-2 activities were evaluated in plaque reduction assays using Vero CCL81 cells. Subtoxic concentrations of compounds 13a, 13c, and 13b displayed in vitro antiviral activity with IC50 in the mid micromolar range. Compounds 13a-c could become lead compounds for the development of new Mpro inhibitors with improved activity against anti-SARS-CoV-2.

Toward α-1,3/4 fucosyltransferases targeted drug discovery: In silico uncovering of promising natural inhibitors of fucosyltransferase 6.

Sialyl Lewis X (sLex ) antigen is a fucosylated cell-surface glycan that is normally involved in cell-cell interactions. The enhanced expression of sLex on cell surface glycans, which is attributed to the upregulation of fucosyltransferase 6 (FUT6), has been implicated in facilitating metastasis in human colorectal, lung, prostate, and oral cancers. The role that the upregulated FUT6 plays in the progression of tumor to malignancy, with reduced survival rates, makes it a potential target for anticancer drugs. Unfortunately, the lack of experimental structures for FUT6 has hampered the design and development of its inhibitors. In this study, we used in silico techniques to identify potential FUT6 inhibitors. We first modeled the three-dimensional structure of human FUT6 using AlphaFold. Then, we screened the natural compound libraries from the COCONUT database to sort out potential natural products (NPs) with best affinity toward the FUT6 model. As a result of these simulations, we identified three NPs for which we predicted binding affinities and interaction patterns quite similar to those we calculated for two experimentally tested FUT6 inhibitors, that is, fucose mimetic-1 and a GDP-triazole derived compound. We also performed molecular dynamics (MD) simulations for the FUT6 complexes with identified NPs, to investigate their stability. Analysis of the MD simulations showed that the identified NPs establish stable contacts with FUT6 under dynamics conditions. On these grounds, the three screened compounds appear as promising natural alternatives to experimentally tested FUT6 synthetic inhibitors, with expected comparable binding affinity. This envisages good prospects for future experimental validation toward FUT6 inhibition.

Towards Decoding Hepatotoxicity of Approved Drugs through Navigation of Multiverse and Consensus Chemical Spaces.

Drug-induced liver injury (DILI) is the principal reason for failure in developing drug candidates. It is the most common reason to withdraw from the market after a drug has been approved for clinical use. In this context, data from animal models, liver function tests, and chemical properties could complement each other to understand DILI events better and prevent them. Since the chemical space concept improves decision-making drug design related to the prediction of structure-property relationships, side effects, and polypharmacology drug activity (uniquely mentioning the most recent advances), it is an attractive approach to combining different phenomena influencing DILI events (e.g., individual "chemical spaces") and exploring all events simultaneously in an integrated analysis of the DILI-relevant chemical space. However, currently, no systematic methods allow the fusion of a collection of different chemical spaces to collect different types of data on a unique chemical space representation, namely "consensus chemical space." This study is the first report that implements data fusion to consider different criteria simultaneously to facilitate the analysis of DILI-related events. In particular, the study highlights the importance of analyzing together in vitro and chemical data (e.g., topology, bond order, atom types, presence of rings, ring sizes, and aromaticity of compounds encoded on RDKit fingerprints). These properties could be aimed at improving the understanding of DILI events.

Consensus Virtual Screening Protocol Towards the Identification of Small Molecules Interacting with the Colchicine Binding Site of the Tubulin-microtubule System.

Modification of the tubulin-microtubule (Tub-Mts) system has generated effective strategies for developing different treatments for cancer. A huge amount of clinical data about inhibitors of the tubulin-microtubule system have supported and validated the studies on this pharmacological target. However, many tubulin-microtubule inhibitors have been developed from representative and common scaffolds that cover a small region of the chemical space with limited structural innovation. The main goal of this study is to develop the first consensus virtual screening protocol for natural products (ligand- and structure-based drug design methods) tuned for the identification of new potential inhibitors of the Tub-Mts system. A combined strategy that involves molecular similarity, molecular docking, pharmacophore modeling, and in silico ADMET prediction has been employed to prioritize the selections of potential inhibitors of the Tub-Mts system. Five compounds were selected and further studied using molecular dynamics and binding energy predictions to characterize their possible binding mechanisms. Their structures correspond to 5-[2-(4-hydroxy-3-methoxyphenyl) ethyl]-2,3-dimethoxyphenol (1), 9,10-dihydro-3,4-dimethoxy-2,7-phenanthrenediol (2), 2-(3,4-dimethoxyphenyl)-5,7-dihydroxy-6-methoxy-4H-1-benzopyran-4-one (3), 13,14-epoxyparvifoline-4',5',6'-trimethoxybenzoate (4), and phenylmethyl 6-hydroxy-2,3-dimethoxybenzoate (5). Compounds 1-3 have been associated with literature reports that confirm their activity against several cancer cell lines, thus supporting the utility of this protocol.

Potential Anticancer Activity of Pomegranate (Punica granatum L.) Fruits of Different Color: In Vitro and In Silico Evidence.

Pomegranate (PMG; Punica granatum L.) fruits possess a well-balanced nutrient/phytochemical composition, with proven adjuvant benefits in experimental cancer chemotherapy; however, such bioactivity could be affected by PMG's phenogenotype (varietal). Here, the chemical and phytochemical (UPLC-DAD-MS2) composition, antioxidant capacity and anticancer potential [in vitro (MTT assay) and in silico (foodinformatics)] of three PMG fruits of different aryl color [red (cv. Wonderful), pink (cv. Molar de Elche), and white (cv. Indian)] were evaluated. The macro/micronutrient (ascorbic acid, tocols, carotenoids), organic acid (citric/malic), and polyphenol content were changed by PMG's varietal and total antioxidant activity (ABTS, alcoholic > hexane extract) in the order of red > pink > white. However, their in vitro cytotoxicity was the same (IC50 > 200 µg.mL-1) against normal (retinal) and cancer (breast, lung, colorectal) cell lines. Sixteen major phytochemicals were tentatively identified, four of them with a high GI absorption/bioavailability score [Ellagic (pink), vanillic (red), gallic (white) acids, D-(+)-catechin (white)] and three of them with multiple molecular targets [Ellagic (52) > vanillic (32) > gallic (23)] associated with anticancer (at initiation and promotion stages) activity. The anticancer potential of the PMG fruit is phenogenotype-specific, although it could be more effective in nutraceutical formulations (concentrates).

7-Aminoalkoxy-Quinazolines from Epigenetic Focused Libraries Are Potent and Selective Inhibitors of DNA Methyltransferase 1.

Inhibitors of epigenetic writers such as DNA methyltransferases (DNMTs) are attractive compounds for epigenetic drug and probe discovery. To advance epigenetic probes and drug discovery, chemical companies are developing focused libraries for epigenetic targets. Based on a knowledge-based approach, herein we report the identification of two quinazoline-based derivatives identified in focused libraries with sub-micromolar inhibition of DNMT1 (30 and 81 nM), more potent than S-adenosylhomocysteine. Also, both compounds had a low micromolar affinity of DNMT3A and did not inhibit DNMT3B. The enzymatic inhibitory activity of DNMT1 and DNMT3A was rationalized with molecular modeling. The quinazolines reported in this work are known to have low cell toxicity and be potent inhibitors of the epigenetic target G9a. Therefore, the quinazoline-based compounds presented are attractive not only as novel potent inhibitors of DNMTs but also as dual and selective epigenetic agents targeting two families of epigenetic writers.

Methyl benzoate and cinnamate analogs as modulators of DNA methylation in hepatocellular carcinoma.

Phenolic acids represent a large collection of phytochemical molecules present in the plant kingdom; they have an important role as epigenetic regulators, particularly as inhibitors of DNA methylation. In the present study, 14 methyl benzoate and cinnamate analogs were synthesized (11-24). Their cytotoxic activity on hepatocellular carcinoma cells (Hep3B) and immortalized human hepatocyte cells was then evaluated. In addition, its effect on the inhibition of global DNA methylation in Hep3B was also determined. Our results showed that the cinnamic derivatives 11-14 and 20-22 were more potent than the free caffeic acid (IC50 109.7-364.2 µM), being methyl 3,4-dihydroxycinammate (12) the most active with an IC50  = 109.7 ± 0.8 µM. Furthermore, 11-14, 20-23 compounds decreased overall DNA methylation levels by 63% to 97%. The analogs methyl 4-hydroxycinnamate (11), methyl 3,4,5-trimethoxycinnamate (14), methyl 4-methoxycinnamate (21), and methyl 3,4-dimethoxycinnamate (22) showed relevant activities of both cytotoxicity and global DNA methylation inhibition. The molecular docking of 21 and 14 suggested that they partly bind to the SAH-binding pocket of DNA methyltransferase 1. These results emphasize the importance of natural products and their analogs as potential sources of DNA methylation modulating agents.

Tubulin Inhibitors: A Chemoinformatic Analysis Using Cell-Based Data.

Inhibiting the tubulin-microtubules (Tub-Mts) system is a classic and rational approach for treating different types of cancers. A large amount of data on inhibitors in the clinic supports Tub-Mts as a validated target. However, most of the inhibitors reported thus far have been developed around common chemical scaffolds covering a narrow region of the chemical space with limited innovation. This manuscript aims to discuss the first activity landscape and scaffold content analysis of an assembled and curated cell-based database of 851 Tub-Mts inhibitors with reported activity against five cancer cell lines and the Tub-Mts system. The structure-bioactivity relationships of the Tub-Mts system inhibitors were further explored using constellations plots. This recently developed methodology enables the rapid but quantitative assessment of analog series enriched with active compounds. The constellations plots identified promising analog series with high average biological activity that could be the starting points of new and more potent Tub-Mts inhibitors.

Epigenetic Target Fishing with Accurate Machine Learning Models.

Epigenetic targets are of significant importance in drug discovery research, as demonstrated by the eight approved epigenetic drugs for treatment of cancer and the increasing availability of chemogenomic data related to epigenetics. This data represents many structure-activity relationships that have not been exploited thus far to develop predictive models to support medicinal chemistry efforts. Herein, we report the first large-scale study of 26 318 compounds with a quantitative measure of biological activity for 55 protein targets with epigenetic activity. We built predictive models with high accuracy for small molecules' epigenetic target profiling through a systematic comparison of the machine learning models trained on different molecular fingerprints. The models were thoroughly validated, showing mean precisions of up to 0.952 for the epigenetic target prediction task. Our results indicate that the models reported herein have considerable potential to identify small molecules with epigenetic activity. Therefore, our results were implemented as a freely accessible web application.

An in silico pipeline for the discovery of multitarget ligands: A case study for epi-polypharmacology based on DNMT1/HDAC2 inhibition.

The search for novel therapeutic compounds remains an overwhelming task owing to the time-consuming and expensive nature of the drug development process and low success rates. Traditional methodologies that rely on the one drug-one target paradigm have proven insufficient for the treatment of multifactorial diseases, leading to a shift to multitarget approaches. In this emerging paradigm, molecules with off-target and promiscuous interactions may result in preferred therapies. In this study, we developed a general pipeline combining machine learning algorithms and a deep generator network to train a dual inhibitor classifier capable of identifying putative pharmacophoric traits. As a case study, we focused on dual inhibitors targeting DNA methyltransferase 1 (DNMT) and histone deacetylase 2 (HDAC2), two enzymes that play a central role in epigenetic regulation. We used this approach to identify dual inhibitors from a novel large natural product database in the public domain. We used docking and atomistic simulations as complementary approaches to establish the ligand-interaction profiles between the best hits and DNMT1/HDAC2. By using the combined ligand- and structure-based approaches, we discovered two promising novel scaffolds that can be used to simultaneously target both DNMT1 and HDAC2. We conclude that the flexibility and adaptability of the proposed pipeline has predictive capabilities of similar or derivative methods and is readily applicable to the discovery of small molecules targeting many other therapeutically relevant proteins.

Expanding the Structural Diversity of DNA Methyltransferase Inhibitors.

Inhibitors of DNA methyltransferases (DNMTs) are attractive compounds for epigenetic drug discovery. They are also chemical tools to understand the biochemistry of epigenetic processes. Herein, we report five distinct inhibitors of DNMT1 characterized in enzymatic inhibition assays that did not show activity with DNMT3B. It was concluded that the dietary component theaflavin is an inhibitor of DNMT1. Two additional novel inhibitors of DNMT1 are the approved drugs glyburide and panobinostat. The DNMT1 enzymatic inhibitory activity of panobinostat, a known pan inhibitor of histone deacetylases, agrees with experimental reports of its ability to reduce DNMT1 activity in liver cancer cell lines. Molecular docking of the active compounds with DNMT1, and re-scoring with the recently developed extended connectivity interaction features approach, led to an excellent agreement between the experimental IC50 values and docking scores.

Recent progress on cheminformatics approaches to epigenetic drug discovery.

The ability of epigenetic markers to affect genome function has enabled transformative changes in drug discovery, especially in cancer and other emerging therapeutic areas. Concordant with the introduction of the term 'epi-informatics', the size of the epigenetically relevant chemical space has grown substantially and so did the number of applications of cheminformatic methods to epigenetics. Recent progress in epi-informatics has improved our understanding of the structure-epigenetic activity relationships and boosted the development of models predicting novel epigenetic agents. Herein, we review the advances in computational approaches to drug discovery of small molecules with epigenetic modulation profiles, summarize the current chemogenomics data available for epigenetic targets, and provide a perspective on the greater utility of biomedical knowledge mining as a means to advance the epigenetic drug discovery.

Current advances on the development of BET inhibitors: insights from computational methods.

Epigenetics was coined almost 70 years ago for the description of heritable phenotype without altering DNA sequences. Research on the field has uncovered significant roles of such mechanisms, that account for the biogenesis of several diseases. Further studies have led the way for drug development which targets epi-enzymes, mainly for cancer treatment. Of the numerous epi-targets involved with histone acetylation, bromodomains have captured the spotlight of drug discovery focused on novel therapies. However, due to high sequence identity, the development of potent and selective inhibitors poses a significant challenge. Herein, we discuss recent computational developments on BET inhibitors and other methods that may be applied for drug discovery in general. As a proof-of-concept, we discuss a virtual screening to identify novel BET inhibitors based on coumarin derivatives. From public data, we identified putative structure-activity relationships of coumarin scaffold and propose R-group modifications for BET selectivity. Results showed that the optimization and design of novel coumarins could be further explored.

In silico tools to study molecular targets of neglected diseases: inhibition of TcSir2rp3, an epigenetic enzyme of Trypanosoma cruzi.

There is a growing interest to study and address neglected tropical diseases (NTD). To this end, in silico methods can serve as the bridge that connects academy and industry, encouraging the development of future treatments against these diseases. This chapter discusses current challenges in the development of new therapies, available computational methods and successful cases in computer-aided design with particular focus on human trypanosomiasis. Novel targets are also discussed. As a case study, we identify amentoflavone as a potential inhibitor of TcSir2rp3 (sirtuine) from Trypanosoma cruzi (20.03 µM) with a workflow that integrates chemoinformatic approaches, molecular modeling, and theoretical affinity calculations, as well as in vitro assays.

Dimeric phenalenones from Talaromyces sp. (IQ-313) inhibit hPTP1B1-400: Insights into mechanistic kinetics from in vitro and in silico studies.

A critical biological event that contributes to the appearance and progress of cancer and diabetes is the reversible phosphorylation of proteins, a process controlled by protein tyrosine-kinases (PTKs) and protein tyrosine-phosphatases (PTPs). Within the PTPs, PTP1B has gained significant interest since it is a validated target in drug discovery. Indeed, several PTP1B inhibitors have been developed, from both, synthesis and natural products. However, none have been approved by the FDA, due to their poor selectivity and/or pharmacokinetic properties. One of the most significant challenges to the discovery of PTP1B inhibitors (in vitro or in silico) is the use of truncated structures (PTP1B1-300), missing valuable information about the mechanisms of inhibition, and selectivity of ligands. The present study describes the biochemical characterization of a full-length PTP1B (hPTP1B1-400), as well as the description of phenalenones 1-4 and ursolic acid (5) as allosteric modulators. Compounds 1-5 showed inhibitory potential on hPTP1B1-400, with IC50 values ranging from 12.7 to 82.1 µM. Kinetic studies showed that 1 and 5 behave as mixed and non-competitive inhibitors, respectively. Circular dichroism experiments confirmed that 1 and 5 induced conformational changes to hPTP1B1-400. Further insights into the structure of hPTP1B1-400 were obtained from a homology model, which pointed out that the C-terminus (residues 301-400) is highly disordered. Molecular docking with the homologated model suggested that compounds 1 and 3-5 bind to the C-terminal domain, likely inducing conformational changes on the protein. Docking positions of compounds 1, 4, and 5 were refined with molecular dynamics simulations. Importantly, these simulations confirmed the high flexibility of the C-terminus of hPTP1B1-400, as well as the changes to its rigidity when bound to 1, 4, and 5.

D-Peptide Builder: A Web Service to Enumerate, Analyze, and Visualize the Chemical Space of Combinatorial Peptide Libraries.

Peptide-based drug discovery is re-gaining attention in drug discovery. Similarly, combinatorial chemistry continues to be a useful technique for the rapid exploration of chemical space. A current challenge, however, is the enumeration of combinatorial peptide libraries using freely accessible tools. To facilitate the swift enumeration of combinatorial peptide libraries, we introduce herein D-Peptide Builder. In the current version, the user can build up to pentapeptides, linear or cyclic, using the natural pool of 20 amino acids. The user can use non- and/or N-methylated amino acids. The server also enables the rapid visualization of the chemical space of the newly enumerated peptides in comparison with other libraries relevant to drug discovery and preloaded in the server. D-Peptide Builder is freely accessible at http://dpeptidebuilder. quimica.unam.mx:4000/. It is also accessible through the open D-Tools platform (DIFACQUIM Tools for Chemoinformatics https://www.difacquim.com/d-tools/).

Consistent Cell-selective Analog Series as Constellation Luminaries in Chemical Space.

High-throughput screening data of compounds consistently tested against the same panel of cell lines is a rich source of information for interrogating cell-selectivity of a compound. Nevertheless, there is a high risk of false positives for these rapid-testing strategies. Then, a single cell-inactive compound can be mistakenly labeled as highly cell-selective if a false positive occurs in any of the cell assays. More interesting would be the case of a series of analogs, which are structurally related compounds, that have a trend to be active only against a small number of cells. To this end, it is herein proposed a proof-of-concept of a method for finding consistent cell-selective analog series of chemical compounds through analysis of high-throughput cell-compound screening data systematically obtained. Furthermore, statistics for quantifying cell-promiscuity and consistency of an analog series are presented.

The impact of chemoinformatics on drug discovery in the pharmaceutical industry.

Introduction: Even though there have been substantial advances in our understanding of biological systems, research in drug discovery is only just now beginning to utilize this type of information. The single-target paradigm, which exemplifies the reductionist approach, remains a mainstay of drug research today. A deeper view of the complexity involved in drug discovery is necessary to advance on this field.Areas covered: This perspective provides a summary of research areas where cheminformatics has played a key role in drug discovery, including of the available resources as well as a personal perspective of the challenges still faced in the field.Expert opinion: Although great strides have been made in the handling and analysis of biological and pharmacological data, more must be done to link the data to biological pathways. This is crucial if one is to understand how drugs modify disease phenotypes, although this will involve a shift from the single drug/single target paradigm that remains a mainstay of drug research. Moreover, such a shift would require an increased awareness of the role of physiology in the mechanism of drug action, which will require the introduction of new mathematical, computer, and biological methods for chemoinformaticians to be trained in.