Computational Approaches and Drug Discovery: Where Are We Going?

Science is a point of view [...].

Natural products as non-covalent and covalent modulators of the KEAP1/NRF2 pathway exerting antioxidant effects.

By controlling several antioxidant and detoxifying genes at the transcriptional level, including NAD(P)H quinone oxidoreductase 1 (NQO1), multidrug resistance-associated proteins (MRPs), UDP-glucuronosyltransferase (UGT), glutamate-cysteine ligase catalytic (GCLC) and modifier (GCLM) subunits, glutathione S-transferase (GST), sulfiredoxin1 (SRXN1), and heme-oxygenase-1 (HMOX1), the KEAP1/NRF2 pathway plays a crucial role in the oxidative stress response. Accordingly, the discovery of modulators of this pathway, activating cellular signaling through NRF2, and targeting the antioxidant response element (ARE) genes is pivotal for the development of effective antioxidant agents. In this context, natural products could represent promising drug candidates for supplementation to provide antioxidant capacity to human cells. In recent decades, by coupling in silico and experimental methods, several natural products have been characterized to exert antioxidant effects by targeting the KEAP1/NRF2 pathway. In this review article, we analyze several natural products that were investigated experimentally and in silico for their ability to modulate KEAP1/NRF2 by non-covalent and covalent mechanisms. These latter represent the two main sections of this article. For each class of inhibitors, we reviewed their antioxidant effects and potential therapeutic applications, and where possible, we analyzed the structure-activity relationship (SAR). Moreover, the main computational techniques used for the most promising identified compounds are detailed in this survey, providing an updated view on the development of natural products as antioxidant agents.

Virtual Screening Strategy and In Vitro Tests to Identify New Inhibitors of the Immunoproteasome.

Immunoproteasome inhibition is a promising strategy for the treatment of hematological malignancies, autoimmune diseases, and inflammatory diseases. The design of non-covalent inhibitors of the immunoproteasome ß1i/ß5i catalytic subunits could be a novel approach to avoid the drawbacks of the known covalent inhibitors, such as toxicity due to off-target binding. In this work, we report the biological evaluation of thirty-four compounds selected from a commercially available collection. These hit compounds are the outcomes of a virtual screening strategy including a dynamic pharmacophore modeling approach onto the ß1i subunit and a pharmacophore/docking approach onto the ß5i subunit. The computational studies were first followed by in vitro enzymatic assays at 100 µM. Only compounds capable of inhibiting the enzymatic activity by more than 50% were characterized in detail using Tian continuous assays, determining the dissociation constant (Ki) of the non-covalent complex where Ki is also the measure of the binding affinity. Seven out of thirty-four hits showed to inhibit ß1i and/or ß5i subunit. Compound 3 is the most active on the ß1i subunit with Ki = 11.84 ± 1.63 µM, and compound 17 showed Ki = 12.50 ± 0.77 µM on the ß5i subunit. Compound 2 showed inhibitory activity on both subunits (Ki = 12.53 ± 0.18 and Ki = 31.95 ± 0.81 on the ß1i subunit and ß5i subunit, respectively). The induced fit docking analysis revealed interactions with Thr1 and Phe31 of ß1i subunit and that represent new key residues as reported in our previous work. Onto ß5i subunit, it interacts with the key residues Thr1, Thr21, and Tyr169. This last hit compound identified represents an interesting starting point for further optimization of ß1i/ß5i dual inhibitors of the immunoproteasome.

Investigating the Inhibition of FTSJ1, a Tryptophan tRNA-Specific 2'-O-Methyltransferase by NV TRIDs, as a Mechanism of Readthrough in Nonsense Mutated CFTR.

Cystic Fibrosis (CF) is an autosomal recessive genetic disease caused by mutations in the CFTR gene, coding for the CFTR chloride channel. About 10% of the CFTR gene mutations are "stop" mutations that generate a premature termination codon (PTC), thus synthesizing a truncated CFTR protein. A way to bypass PTC relies on ribosome readthrough, which is the ribosome's capacity to skip a PTC, thus generating a full-length protein. "TRIDs" are molecules exerting ribosome readthrough; for some, the mechanism of action is still under debate. We investigate a possible mechanism of action (MOA) by which our recently synthesized TRIDs, namely NV848, NV914, and NV930, could exert their readthrough activity by in silico analysis and in vitro studies. Our results suggest a likely inhibition of FTSJ1, a tryptophan tRNA-specific 2'-O-methyltransferase.

Synthesis of new antiproliferative 1,3,4-substituted-pyrrolo[3,2-c]quinoline derivatives, biological and in silico insights.

A series of biologically unexplored substituted 1,3,4-subtituted-pyrrolo[3,2-c]quinoline derivatives (PQs) was evaluated against a panel of about 60 tumor cells (NCI). Based on the preliminary antiproliferative data, the optimizations efforts permitted us to design and synthesize a new series of derivatives allowing us to individuate a promising hit (4g). The insertion of a 4-benzo[d] [1,3]dioxol-5-yl moiety on increased and extended the activity towards five panel tumor cell lines such as leukemia, CNS, melanoma, renal and breast cancer, reaching IG50 in the low µM range. Replacement of this latter with a 4-(OH-di-Cl-Ph) group (4i) or introduction a Cl-propyl chain in position 1 (5), selectively addressed the activity against the entire leukemia sub-panel (CCRF-CEM, K-562, MOLT-4, RPMI-8226, SR). Preliminary biological assays on MCF-7 such as cell cycle, clonogenic assay, ROS content test alongside a comparison of viability between MCF-7 and non-tumorigenic MCF-10 were investigated. Among the main anticancer targets involved in breast cancer, HSP90 and ER receptors were selected for in silico studies. Docking analysis revealed a valuable affinity for HSP90 providing structural insights on the binding mode, and useful features for optimization.

Fighting Antibiotic Resistance: New Pyrimidine-Clubbed Benzimidazole Derivatives as Potential DHFR Inhibitors.

The present work describes the design and development of seventeen pyrimidine-clubbed benzimidazole derivatives as potential dihydrofolate reductase (DHFR) inhibitors. These compounds were filtered by using ADMET, drug-likeness characteristics calculations, and molecular docking experiments. Compounds 27, 29, 30, 33, 37, 38, and 41 were chosen for the synthesis based on the results of the in silico screening. Each of the synthesized compounds was tested for its in vitro antibacterial and antifungal activities using a variety of strains. All the compounds showed antibacterial properties against Gram-positive bacteria (Staphylococcus aureus and Staphylococcus pyogenes) as well as Gram-negative bacteria (Escherichia coli and Pseudomonas aeruginosa). Most of the compounds either had a higher potency than chloramphenicol or an equivalent potency to ciprofloxacin. Compounds 29 and 33 were effective against all the bacterial and fungal strains. Finally, the 1,2,3,4-tetrahydropyrimidine-2-thiol derivatives with a 6-chloro-2-(chloromethyl)-1H-benzo[d]imidazole moiety are potent enough to be considered a promising lead for the discovery of an effective antibacterial agent.

Sulfonamide Moiety as "Molecular Chimera" in the Design of New Drugs.

BACKGROUND: The -SO2NH- group is of great significance in modern pharmaceutical use since, in sulfa-drugs, it is possible to introduce easily chemical modifications, and even small changes may lead to an improved version of an already existing drug. OBJECTIVE: This paper aims to describe updated information in the sulfonamide field with a particular focus on new mechanisms of action, especially if discovered by employing computational approaches. METHODS: Research articles that focused on the use of the sulfonamide moiety for the design, synthesis, and in vitro/in vivo tests of various diseases were collected from various search engines like PubMed, Science Direct, Google Scholar, and Scopus, using keywords like sulfonamide moiety, aryl/heteroary lsulfonamides, alkyl sulfonamides, in silico drug design, etc. Conclusion: The more relevant reports highlighting the prominent role of sulfonamide moiety in drug discovery have been critically analyzed. Sulfonamides can be considered as "molecular chimera", which are found to form hydrogen bonds as well as interact with unipolar environments within proteins. Therefore, based on the analysis reported herein, it is strongly foresight that new entities can be developed easily to improve the available machinery helpful in the fight against new and emerging diseases.

In Silico Design, Synthesis, and Biological Evaluation of Anticancer Arylsulfonamide Endowed with Anti-Telomerase Activity.

Telomerase, a reverse transcriptase enzyme involved in DNA synthesis, has a tangible role in tumor progression. Several studies have evidenced telomerase as a promising target for developing cancer therapeutics. The main reason is due to the overexpression of telomerase in cancer cells (85-90%) compared with normal cells where it is almost unexpressed. In this paper, we used a structure-based approach to design potential inhibitors of the telomerase active site. The MYSHAPE (Molecular dYnamics SHared PharmacophorE) approach and docking were used to screen an in-house library of 126 arylsulfonamide derivatives. Promising compounds were synthesized using classical and green methods. Compound 2C revealed an interesting IC50 (33 ± 4 µM) against the K-562 cell line compared with the known telomerase inhibitor BIBR1532 IC50 (208 ± 11 µM) with an SI ~10 compared to the BALB/3-T3 cell line. A 100 ns MD simulation of 2C in the telomerase active site evidenced Phe494 as the key residue as well as in BIBR1532. Each moiety of compound 2C was involved in key interactions with some residues of the active site: Arg557, Ile550, and Gly553. Compound 2C, as an arylsulfonamide derivative, is an interesting hit compound that deserves further investigation in terms of optimization of its structure to obtain more active telomerase inhibitors.

Computational Approaches: Drug Discovery and Design in Medicinal Chemistry and Bioinformatics.

To date, computational approaches have been recognized as a key component in drug design and discovery workflows [...].

Evaluation of the IKKß Binding of Indicaxanthin by Induced-Fit Docking, Binding Pose Metadynamics, and Molecular Dynamics.

Background: Indicaxanthin, a betaxanthin belonging to the betalain class of compounds, has been recently demonstrated to exert significant antiproliferative effects inducing apoptosis of human melanoma cells through the inhibition of NF-κB as the predominant pathway. Specifically, Indicaxanthin inhibited IκBα degradation in A375 cells. In resting cells, NF-κB is arrested in the cytoplasm by binding to its inhibitor protein IκBα. Upon stimulation, IκBα is phosphorylated by the IKK complex, and degraded by the proteasome, liberating free NF-κB into the nucleus to initiate target gene transcription. Inhibition of the IKK complex leads to the arrest of the NF-κB pathway. Methods: To acquire details at the molecular level of Indicaxanthin's inhibitory activity against hIKKß, molecular modeling and simulation techniques including induced-fit docking (IFD), binding pose metadynamics (BPMD), molecular dynamics simulations, and MM-GBSA (molecular mechanics-generalized Born surface area continuum solvation) have been performed. Results: The computational calculations performed on the active and inactive form, and the allosteric binding site of hIKKß, revealed that Indicaxanthin inhibits prevalently the active form of the hIKKß. MM-GBSA computations provide further evidence of Indicaxanthin's stability inside the active binding pocket with a binding free energy of -22.2 ± 4.3 kcal/mol with respect to the inactive binding pocket with a binding free energy of -20.7 ± 4.7 kcal/mol. BPMD and MD simulation revealed that Indicaxanthin is likely not an allosteric inhibitor of hIKKß. Conclusion: As a whole, these in silico pieces of evidence show that Indicaxanthin can inhibit the active form of the hIKKß adding novel mechanistic insights on its recently discovered ability to impair NF-κB signaling in melanoma A375 cells. Moreover, our results suggest the phytochemical as a new lead compound for novel, more potent IKKß inhibitors to be employed in the treatment of cancer and inflammation-related conditions.

Immunoproteasome and Non-Covalent Inhibition: Exploration by Advanced Molecular Dynamics and Docking Methods.

The selective inhibition of immunoproteasome is a valuable strategy to treat autoimmune, inflammatory diseases, and hematologic malignancies. Recently, a new series of amide derivatives as non-covalent inhibitors of the ß1i subunit with Ki values in the low/submicromolar ranges have been identified. Here, we investigated the binding mechanism of the most potent and selective inhibitor, N-benzyl-2-(2-oxopyridin-1(2H)-yl)propanamide (1), to elucidate the steps from the ligand entrance into the binding pocket to the ligand-induced conformational changes. We carried out a total of 400 ns of MD-binding analyses, followed by 200 ns of plain MD. The trajectories clustering allowed identifying three representative poses evidencing new key interactions with Phe31 and Lys33 together in a flipped orientation of a representative pose. Further, Binding Pose MetaDynamics (BPMD) studies were performed to evaluate the binding stability, comparing 1 with four other inhibitors of the ß1i subunit: N-benzyl-2-(2-oxopyridin-1(2H)-yl)acetamide (2), N-cyclohexyl-3-(2-oxopyridin-1(2H)-yl)propenamide (3), N-butyl-3-(2-oxopyridin-1(2H)-yl)propanamide (4), and (S)-2-(2-oxopyridin-1(2H)-yl)-N,4-diphenylbutanamide (5). The obtained results in terms of free binding energy were consistent with the experimental values of inhibition, confirming 1 as a lead compound of this series. The adopted methods provided a full dynamic description of the binding events, and the information obtained could be exploited for the rational design of new and more active inhibitors.

Antiproliferative Properties and G-Quadruplex-Binding of Symmetrical Naphtho[1,2-b:8,7-b']dithiophene Derivatives.

Background: G-quadruplex (G4) forming sequences are recurrent in telomeres and promoter regions of several protooncogenes. In normal cells, the transient arrangements of DNA in G-tetrads may regulate replication, transcription, and translation processes. Tumors are characterized by uncontrolled cell growth and tissue invasiveness and some of them are possibly mediated by gene expression involving G-quadruplexes. The stabilization of G-quadruplex sequences with small molecules is considered a promising strategy in anticancer targeted therapy. Methods: Molecular virtual screening allowed us identifying novel symmetric bifunctionalized naphtho[1,2-b:8,7-b']dithiophene ligands as interesting candidates targeting h-Telo and c-MYC G-quadruplexes. A set of unexplored naphtho-dithiophene derivatives has been synthesized and biologically tested through in vitro antiproliferative assays and spectroscopic experiments in solution. Results: The analysis of biological and spectroscopic data highlighted noteworthy cytotoxic effects on HeLa cancer cell line (GI50 in the low µM range), but weak interactions with G-quadruplex c-MYC promoter. Conclusions: The new series of naphtho[1,2-b:8,7-b']dithiophene derivatives, bearing the pharmacophoric assumptions necessary to stabilize G-quadruplexes, have been designed and successfully synthesized. The interesting antiproliferative results supported by computer aided rational approaches suggest that these studies are a significant starting point for a lead optimization process and the isolation of a more efficacious set of G-quadruplexes stabilizers.

Pharmacophore-Based Design of New Chemical Scaffolds as Translational Readthrough-Inducing Drugs (TRIDs).

Translational readthrough-inducing drugs (TRIDs) rescue the functional full-length protein expression in genetic diseases, such as cystic fibrosis, caused by premature termination codons (PTCs). Small molecules have been developed as TRIDs to trick the ribosomal machinery during recognition of the PTC. Herein we report a computational study to identify new TRID scaffolds. A pharmacophore approach was carried out on compounds that showed readthrough activity. The pharmacophore model applied to screen different libraries containing more than 87000 compounds identified four hit-compounds presenting scaffolds with diversity from the oxadiazole lead. These compounds have been synthesized and tested using the Fluc reporter harboring the UGA PTC. Moreover, the cytotoxic effect and the expression of the CFTR protein were evaluated. These compounds, a benzimidazole derivative (NV2899), a benzoxazole derivative (NV2913), a thiazole derivative (NV2909), and a benzene-1,3-disulfonate derivative (NV2907), were shown to be potential new lead compounds as TRIDs, boosting further efforts to address the optimization of the chemical scaffolds.

The In Silico Fischer Lock-and-Key Model: The Combined Use of Molecular Descriptors and Docking Poses for the Repurposing of Old Drugs.

Not always lead compound and/or derivatives are suitable for the specific biological target for which they are designed but, in some cases, discarded compounds proved to be good binders for other biological targets; therefore, drug repurposing constitute a valid alternative to avoid waste of human and financial resources. Our virtual lock-and-key methods, VLKA and Conf-VLKA, furnish a strong support to predict the efficacy of a designed drug a priori its biological evaluation, or the correct biological target for a set of the selected compounds, allowing thus the repurposing of known and unknown, active and inactive compounds.

Indicaxanthin, a multi-target natural compound from Opuntia ficus-indica fruit: From its poly-pharmacological effects to biochemical mechanisms and molecular modelling studies.

Over the latest years phytochemical consumption has been associated to a decreased risk of both the onset and the development of a number of pathological conditions. In this context indicaxanthin, a betalain pigment from Opuntia ficus-indica fruit, has been the object of sound research. Explored, at first, for its mere antioxidant potential, Indicaxanthin is now regarded as a redox-active compound able to exert significant poly-pharmacological effects against several targets in a number of experimental conditions both in vivo and in vitro. This paper aims to provide an overview on the therapeutical effects of indicaxanthin, ranging from the anti-inflammatory to the neuro-modulatory and anti-tumoral ones and favored by its high bioavailability. Moreover, biochemical and molecular modelling investigations are aimed to identify the pharmacological targets the compound is able to interact with and to address the challenging development in the future research.

Deciphering the Nonsense Readthrough Mechanism of Action of Ataluren: An in Silico Compared Study.

Ataluren was reported to suppress nonsense mutations by promoting the readthrough of premature stop codons, although its mechanism of action (MOA) is still debated. The likely interaction of Ataluren with CFTR-mRNA has been previously studied by molecular dynamics. In this work we extended the modeling of Ataluren's MOA by complementary computational approaches such as induced fit docking (IFD), quantum polarized ligand docking (QPLD), MM-GBSA free-energy calculations, and computational mutagenesis. In addition to CFTR-mRNA, this study considered other model targets implicated in the translation process, such as eukaryotic rRNA 18S, prokaryotic rRNA 16S, and eukaryotic Release Factor 1 (eRF1), and we performed a comparison with a new promising Ataluren analogue (NV2445) and with a series of aminoglycosides, known to suppress the normal proofreading function of the ribosome. Results confirmed mRNA as the most likely candidate target for Ataluren and its analogue, and binding energies calculated after computational mutagenesis highlighted how Ataluren's interaction with the premature stop codon could be affected by ancillary nucleotides in the genetic context.

Predicting Skin Permeability by Means of Computational Approaches: Reliability and Caveats in Pharmaceutical Studies.

The skin is the main barrier between the internal body environment and the external one. The characteristics of this barrier and its properties are able to modify and affect drug delivery and chemical toxicity parameters. Therefore, it is not surprising that permeability of many different compounds has been measured through several in vitro and in vivo techniques. Moreover, many different in silico approaches have been used to identify the correlation between the structure of the permeants and their permeability, to reproduce the skin behavior, and to predict the ability of specific chemicals to permeate this barrier. A significant number of issues, like interlaboratory variability, experimental conditions, data set building rationales, and skin site of origin and hydration, still prevent us from obtaining a definitive predictive skin permeability model. This review wants to show the main advances and the principal approaches in computational methods used to predict this property, to enlighten the main issues that have arisen, and to address the challenges to develop in future research.

Reverse screening on indicaxanthin from Opuntia ficus-indica as natural chemoactive and chemopreventive agent.

Indicaxanthin is a bioactive and bioavailable betalain pigment extracted from Opuntia ficus indica fruits. Indicaxanthin has pharmacokinetic proprieties, rarely found in other phytochemicals, and it has been demonstrated that it provides a broad-spectrum of pharmaceutical activity, exerting anti-proliferative, anti-inflammatory, and neuromodulator effects. The discovery of the Indicaxanthin physiological targets plays an important role in understanding the biochemical mechanism. In this study, combined reverse pharmacophore mapping, reverse docking, and text-based database search identified Inositol Trisphosphate 3-Kinase (ITP3K-A), Glutamate carboxypeptidase II (GCPII), Leukotriene-A4 hydrolase (LTA4H), Phosphoserine phosphatase (HPSP), Phosphodiesterase 4D (PDE4D), AMPA receptor (GluA3 and GluA2 subunits) and Kainate receptor (GluK1 isoform) as potential targets for Indicaxanthin. These targets are implicated in neuromodulation, and inflammatory regulation, normally expressed mostly in the CNS, and expressed (or overexpressed) in cancer tissues (i.e. breast, thyroid, and prostate cancer cells). Moreover, this study provides qualitative and quantitative information about dynamic interactions of Indicaxanthin at the binding site of target proteins, through molecular dynamics simulations and MM-GBSA.