Computational Approaches and Drug Discovery: Where Are We Going?

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

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.

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.

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 [...].

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.

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.

Micelles of the chiral biocompatible surfactant (1R,2S)-dodecyl(2-hydroxy-1-methyl-2-phenylethyl)dimethylammonium bromide (DMEB): molecular dynamics and fragmentation patterns in the gas phase.

RATIONALE: The study of self-assembly processes of surfactant molecules in the gas phase is of great interest for several theoretical and technological reasons related to their possible exploitation as drug carriers, protein shields and cleaning agents in the gas phase. METHODS: The stability and fragmentation patterns of singly and multiply charged (either positively or negatively) aggregates of the surfactant (1R,2S)-dodecyl(2-hydroxy-1-methyl-2-phenylethyl)dimethyl ammonium bromide (DMEB) in the gas phase have been studied by ion mobility mass spectrometry and tandem mass spectrometry. Molecular dynamics (MD) simulations of positively and negatively singly and multiply charged DMEB aggregates have been performed to obtain structural and energetics information. Finally, in order to ascertain some clues on the DMEB growth mechanism, quantum mechanics calculations were carried out. RESULTS: It has been evidenced that positively and negatively singly charged aggregates at low collision energy decompose preferentially by loss of only one DMEB molecule. Increasing the collision energy, the loss of neutrals becomes increasingly abundant. Multiply charged DMEB aggregates are unstable and decompose forming singly charged monomers or dimers. MD simulations show reverse micelle-like structures with polar heads somewhat segregated into the aggregate interior. Finally, a good correlation between experimental and calculated collisional cross sections (CCS) was found. CONCLUSIONS: The fragmentation pathways of DMEB charged species evidenced for singly charged aggregates exhibit features similar to that of other detergent aggregates, but multiply charged aggregates showed a system-specific behavior. QM calculations on the optimized structures (21+ , 31+ , 11- and 21- ) indicate that the most determinant interactions are due to an OH---Br hydrogen bonding that is also involved in the link between monomeric DMEB units. The MD models gave CCS values in good agreement with experimental ones, evidenced by a less strict reverse micelle-like structure and a reasonably spread bromine anion distribution Copyright © 2017 John Wiley & Sons, Ltd.

Common Hits Approach: Combining Pharmacophore Modeling and Molecular Dynamics Simulations.

We present a new approach that incorporates flexibility based on extensive MD simulations of protein-ligand complexes into structure-based pharmacophore modeling and virtual screening. The approach uses the multiple coordinate sets saved during the MD simulations and generates for each frame a pharmacophore model. Pharmacophore models with the same pharmacophore features are pooled. In this way the high number of pharmacophore models that results from the MD simulation is reduced to only a few hundred representative pharmacophore models. Virtual screening runs are performed with every representative pharmacophore model; the screening results are combined and rescored to generate a single hit-list. The score for a particular molecule is calculated based on the number of representative pharmacophore models which classified it as active. Hence, the method is called common hits approach (CHA). The steps between the MD simulation and the final hit-list are performed automatically and without user interaction. We test the performance of CHA for virtual screening using screening databases with active and inactive compounds for 40 protein-ligand systems. The results of the CHA are compared to the (i) median screening performance of all representative pharmacophore models of protein-ligand systems, as well as to the virtual screening performance of (ii) a random classifier, (iii) the pharmacophore model derived from the experimental structure in the PDB, and (iv) the representative pharmacophore model appearing most frequently during the MD simulation. For the 34 (out of 40) protein-ligand complexes, for which at least one of the approaches was able to perform better than a random classifier, the highest enrichment was achieved using CHA in 68% of the cases, compared to 12% for the PDB pharmacophore model and 20% for the representative pharmacophore model appearing most frequently. The availabilithy of diverse sets of different pharmacophore models is utilized to analyze some additional questions of interest in 3D pharmacophore-based virtual screening.

New benzothieno[3,2-d]-1,2,3-triazines with antiproliferative activity: synthesis, spectroscopic studies, and biological activity.

New benzothieno[3,2-d]-1,2,3-triazines, together with precursors triazenylbenzo[b]thiophenes, were designed, synthesized and screened as anticancer agents. The structural features of these compounds prompted us to investigate their DNA binding capability through UV-vis absorption titrations, circular dichroism, and viscometry, pointing out the occurrence of groove-binding. The derivative 3-(4-methoxy-phenyl)benzothieno[3,2-d]-1,2,3-triazin-4(3H)-one showed the highest antiproliferative effect against HeLa cells and was also tested in cell cycle perturbation experiments. The obtained results assessed for the first time the anticancer activity of benzothieno[3,2-d]-1,2,3-triazine nucleus, and we related it to its DNA-binding properties.

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.

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.

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.

Molecular dynamics studies on Mdm2 complexes: an analysis of the inhibitor influence.

p53 is a powerful anti-tumoral molecule frequently inactivated by mutations or deletions in cancer. However, half of all human tumors expresses wild-type p53, and its activation, by antagonizing its negative regulator Mdm2, might offer a new strategy for therapeutic protocol. In this work, we present a molecular dynamics study on Mdm2 structure bound to two different known inhibitors with the aim to investigate the structural transitions between apo-Mdm2 and Mdm2-inhibitor complexes. We tried to gain information about conformational changes binding a benzodiazepine derivative inhibitor with respect the known nutlin and the apo form. The conformational changes alter the size of the cleft and were mainly in the linker regions, suggesting that the overall dynamic nature of Mdm2 is related to dynamic movements in these regions.

In vitro and in silico studies of polycondensed diazine systems as anti-parasitic agents.

Parasitic diseases caused by protozoarian agents are still relevant today more than ever. Recently, we synthesized several polycondensed diazine derivatives by means 1,3-dipolar cycloaddition reactions. A broad selection of these compounds were submitted to in vitro biological screening against Plasmodium falciparum, Leishmania infantum, Trypanosoma brucei, and Trypanosoma cruzi, resulting active at micromolar level. Induced Fit Docking/MM-GBSA studies were performed giving interesting indications about the probable mechanism of action of the most active compounds.

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.

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.

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.