"Hit" to lead optimization and chemoinformatic studies for a new series of Autotaxin inhibitors.

Robust experimental evidence has highlighted the role of Autotaxin (ATX)/Lysophosphatidic acid (LPA) axis not only in the pathogenesis of chronic inflammatory conditions and especially in fibroproliferative diseases but also in several types of cancer. As a result, different series of substrate-, lipid-based and small-molecule ATX inhibitors have been identified thus far by both academia and pharma. The "crowning achievement" of these drug discovery campaigns was the development and entry of the first-in-class ATX inhibitor (ziritaxestat, GLPG-1690) in advanced clinical trials against idiopathic pulmonary fibrosis. Herein, the potency optimization efforts of a new series of Autotaxin inhibitors, namely 2-substituted-2,6-dihydro-4H-thieno[3,4-c]pyrazol-1-substituted amide, is described using a previously identified novel chemical scaffold as a "hit". The mode of inhibition of the most promising ATX inhibitors was investigated, while their cellular activity, aqueous solubility and cytotoxicity were evaluated. Our pharmacological results were corroborated by chemoinformatic tools (molecular docking and molecular dynamics simulations) deployed, to provide insight into the binding mechanism of the synthesized inhibitors to ATX.

In Silico Identification and Evaluation of Natural Products as Potential Tumor Necrosis Factor Function Inhibitors Using Advanced Enalos Asclepios KNIME Nodes.

Tumor necrosis factor (TNF) is a regulator of several chronic inflammatory diseases, such as rheumatoid arthritis. Although anti-TNF biologics have been used in clinic, they render several drawbacks, such as patients' progressive immunodeficiency and loss of response, high cost, and intravenous administration. In order to find new potential anti-TNF small molecule inhibitors, we employed an in silico approach, aiming to find natural products, analogs of Ampelopsin H, a compound that blocks the formation of TNF active trimer. Two out of nine commercially available compounds tested, Nepalensinol B and Miyabenol A, efficiently reduced TNF-induced cytotoxicity in L929 cells and production of chemokines in mice joints' synovial fibroblasts, while Nepalensinol B also abolished TNF-TNFR1 binding in non-toxic concentrations. The binding mode of the compounds was further investigated by molecular dynamics and free energy calculation studies, using and advancing the Enalos Asclepios pipeline. Conclusively, we propose that Nepalensinol B, characterized by the lowest free energy of binding and by a higher number of hydrogen bonds with TNF, qualifies as a potential lead compound for TNF inhibitors' drug development. Finally, the upgraded Enalos Asclepios pipeline can be used for improved identification of new therapeutics against TNF-mediated chronic inflammatory diseases, providing state-of-the-art insight on their binding mode.

Ebselen, Disulfiram, Carmofur, PX-12, Tideglusib, and Shikonin Are Nonspecific Promiscuous SARS-CoV-2 Main Protease Inhibitors.

Among the drug targets being investigated for SARS-CoV-2, the viral main protease (Mpro) is one of the most extensively studied. Mpro is a cysteine protease that hydrolyzes the viral polyprotein at more than 11 sites. It is highly conserved and has a unique substrate preference for glutamine in the P1 position. Therefore, Mpro inhibitors are expected to have broad-spectrum antiviral activity and a high selectivity index. Structurally diverse compounds have been reported as Mpro inhibitors. In this study, we investigated the mechanism of action of six previously reported Mpro inhibitors, ebselen, disulfiram, tideglusib, carmofur, shikonin, and PX-12, using a consortium of techniques including FRET-based enzymatic assay, thermal shift assay, native mass spectrometry, cellular antiviral assays, and molecular dynamics simulations. Collectively, the results showed that the inhibition of Mpro by these six compounds is nonspecific and that the inhibition is abolished or greatly reduced with the addition of reducing reagent 1,4-dithiothreitol (DTT). Without DTT, these six compounds inhibit not only Mpro but also a panel of viral cysteine proteases including SARS-CoV-2 papain-like protease and 2Apro and 3Cpro from enterovirus A71 (EV-A71) and EV-D68. However, none of the compounds inhibits the viral replication of EV-A71 or EV-D68, suggesting that the enzymatic inhibition potency IC50 values obtained in the absence of DTT cannot be used to faithfully predict their cellular antiviral activity. Overall, we provide compelling evidence suggesting that these six compounds are nonspecific SARS-CoV-2 Mpro inhibitors and urge the scientific community to be stringent with hit validation.

Structure and inhibition of the SARS-CoV-2 main protease reveal strategy for developing dual inhibitors against Mpro and cathepsin L.

The main protease (Mpro) of SARS-CoV-2 is a key antiviral drug target. While most Mpro inhibitors have a γ-lactam glutamine surrogate at the P1 position, we recently found that several Mpro inhibitors have hydrophobic moieties at the P1 site, including calpain inhibitors II and XII, which are also active against human cathepsin L, a host protease that is important for viral entry. In this study, we solved x-ray crystal structures of Mpro in complex with calpain inhibitors II and XII and three analogs of GC-376 The structure of Mpro with calpain inhibitor II confirmed that the S1 pocket can accommodate a hydrophobic methionine side chain, challenging the idea that a hydrophilic residue is necessary at this position. The structure of calpain inhibitor XII revealed an unexpected, inverted binding pose. Together, the biochemical, computational, structural, and cellular data presented herein provide new directions for the development of dual inhibitors as SARS-CoV-2 antivirals.

Pharmacological characterisation of novel adenosine A3 receptor antagonists.

The adenosine A3 receptor (A3R) belongs to a family of four adenosine receptor (AR) subtypes which all play distinct roles throughout the body. A3R antagonists have been described as potential treatments for numerous diseases including asthma. Given the similarity between (adenosine receptors) orthosteric binding sites, obtaining highly selective antagonists is a challenging but critical task. Here we screen 39 potential A3R, antagonists using agonist-induced inhibition of cAMP. Positive hits were assessed for AR subtype selectivity through cAMP accumulation assays. The antagonist affinity was determined using Schild analysis (pA2 values) and fluorescent ligand binding. Structure-activity relationship investigations revealed that loss of the 3-(dichlorophenyl)-isoxazolyl moiety or the aromatic nitrogen heterocycle with nitrogen at α-position to the carbon of carboximidamide group significantly attenuated K18 antagonistic potency. Mutagenic studies supported by molecular dynamic simulations combined with Molecular Mechanics-Poisson Boltzmann Surface Area calculations identified the residues important for binding in the A3R orthosteric site. We demonstrate that K18, which contains a 3-(dichlorophenyl)-isoxazole group connected through carbonyloxycarboximidamide fragment with a 1,3-thiazole ring, is a specific A3R (< 1 µM) competitive antagonist. Finally, we introduce a model that enables estimates of the equilibrium binding affinity for rapidly disassociating compounds from real-time fluorescent ligand-binding studies. These results demonstrate the pharmacological characterisation of a selective competitive A3R antagonist and the description of its orthosteric binding mode. Our findings may provide new insights for drug discovery.

Ebselen, disulfiram, carmofur, PX-12, tideglusib, and shikonin are non-specific promiscuous SARS-CoV-2 main protease inhibitors.

There is an urgent need for vaccines and antiviral drugs to combat the COVID-19 pandemic. Encouraging progress has been made in developing antivirals targeting SARS-CoV-2, the etiological agent of COVID-19. Among the drug targets being investigated, the viral main protease (M pro ) is one of the most extensively studied drug targets. M pro is a cysteine protease that hydrolyzes the viral polyprotein at more than 11 sites and it is highly conserved among coronaviruses. In addition, M pro has a unique substrate preference for glutamine in the P1 position. Taken together, it appears that M pro inhibitors can achieve both broad-spectrum antiviral activity and a high selectivity index. Structurally diverse compounds have been reported as M pro inhibitors, with several of which also showed antiviral activity in cell culture. In this study, we investigated the mechanism of action of six previously reported M pro inhibitors, ebselen, disulfiram, tideglusib, carmofur, shikonin, and PX-12 using a consortium of techniques including FRET-based enzymatic assay, thermal shift assay, native mass spectrometry, cellular antiviral assays, and molecular dynamics simulations. Collectively, the results showed that the inhibition of M pro by these six compounds is non-specific and the inhibition is abolished or greatly reduced with the addition of reducing reagent DTT. In the absence of DTT, these six compounds not only inhibit M pro , but also a panel of viral cysteine proteases including SARS-CoV-2 papain-like protease, the 2A pro and 3C pro from enterovirus A71 (EV-A71) and EV-D68. However, none of the compounds inhibits the viral replication of EV-A71 or EV-D68, suggesting that the enzymatic inhibition potency IC 50 values obtained in the absence of DTT cannot be used to faithfully predict their cellular antiviral activity. Overall, we provide compelling evidence suggesting that ebselen, disulfiram, tideglusib, carmofur, shikonin, and PX-12 are non-specific SARS-CoV-2 M pro inhibitors, and urge the scientific community to be stringent with hit validation.

Insights to the Binding of a Selective Adenosine A3 Receptor Antagonist Using Molecular Dynamic Simulations, MM-PBSA and MM-GBSA Free Energy Calculations, and Mutagenesis.

Adenosine A3 receptor (A3R) is a promising drug target cancer and for a number of other conditions like inflammatory diseases, including asthma and rheumatoid arthritis, glaucoma, chronic obstructive pulmonary disease, and ischemic injury. Currently, there is no experimentally determined structure of A3R. We explored the binding profile of O4-{[3-(2,6-dichlorophenyl)-5-methylisoxazol-4-yl]carbonyl}-2-methyl-1,3-thiazole-4-carbohydroximamide (K18), which is a new specific and competitive antagonist at the orthosteric binding site of A3R. MD simulations and MM-GBSA calculations of the WT A3R in complex with K18 combined with in vitro mutagenic studies show that the most plausible binding conformation for the dichlorophenyl group of K18 is oriented toward trans-membrane helices (TM) 5, 6 and reveal important residues for binding. Further, MM-GBSA calculations distinguish mutations that reduce or maintain or increase antagonistic activity. Our studies show that selectivity of K18 toward A3R is defined not only by direct interactions with residues within the orthosteric binding area but also by remote residues playing a significant role. Although V1695.30 is considered to be a selectivity filter for A3R binders, when it was mutated to glutamic acid, K18 maintained antagonistic potency, in agreement with our previous results obtained for agonists binding profile investigation. Mutation of the direct interacting residue L903.32 in the low region and the remote L2647.35 in the middle/upper region to alanine increases antagonistic potency, suggesting an empty space in the orthosteric area available for increasing antagonist potency. These results approve the computational model for the description of K18 binding at A3R, which we previously performed for agonists binding to A3R, and the design of more effective antagonists based on K18.

Structural Characterization of Agonist Binding to an A3 Adenosine Receptor through Biomolecular Simulations and Mutagenesis Experiments.

The adenosine A3 receptor (A3R) binds adenosine and is a drug target against cancer cell proliferation. Currently, there is no experimental structure of A3R. Here, we have generated a molecular model of A3R in complex with two agonists, the nonselective 1-(6-amino-9H-purin-9-yl)-1-deoxy-N-ethyl-ß-d-ribofuranuronamide (NECA) and the selective 1-deoxy-1-[6-[[(3-iodophenyl)methyl]amino]-9H-purin-9-yl]-N-methyl-ß-d-ribofuranuronamide (IB-MECA). Molecular dynamics simulations of the wild-type A3R in complex with both agonists, combined with in vitro mutagenic studies revealed important residues for binding. Further, molecular mechanics-generalized Born surface area calculations were able to distinguish mutations that reduce or negate agonistic activity from those that maintained or increased the activity. Our studies reveal that selectivity of IB-MECA toward A3R requires not only direct interactions with residues within the orthosteric binding area but also with remote residues. Although V1695.30 is considered to be a selectivity filter for A3R binders, when it was mutated to glutamic acid or alanine, the activity of IB-MECA increased by making new van der Waals contacts with TM5. This result may have implications in the design of new A3R agonists.

Molecular Dynamics Simulations on the Bioactive Molecule of hIAPP22-29 (NFGAILSS) and Rational Drug Design.

This chapter includes information about the structure in equilibrium of the bioactive molecule hIAPP22-29 (NFGAILSS). The experimental structure was derived using X-ray and its 2D NOESY NMR experiments in d 6-DMSO and d-HFIP solvents. This molecule contains eight of the ten amino acids of the 20-29 region of the human islet amyloid polypeptide (hIAPP) often referred as the "amyloidogenic core." Amyloid deposits are well-known to cause as many as 20 pathological neurodegenerative disorders such as Alzheimer, Parkinson, Huntington, and Creutzfeldt-Jakob. The experimental structure was relaxed using molecular dynamics (MD) in simulation boxes consisting in DMSO and HFIP; the latter not provided by the applied software. The calculations were performed in GPUs and supercomputers, and some basic scripting is described for reference. The simulations confirmed the inter- and intramolecular forces that led to an "amyloidogenic core" observed from NOE experiments. The results showed that in DMSO and HFIP environment, Phe is not in spatial proximity with Leu or Ile, and this is consistent with an amyloidogenic core. However, in an amphipathic environment such as the model lipid bilayers, this communication is possible and may influence peptide amyloidogenic properties. The knowledge gained through this study may contribute to the rational drug design of novel peptides or organic molecules acting by modifying preventing amyloidogenic properties of the hIAPP peptide.

Discovery of Novel Adenosine Receptor Antagonists through a Combined Structure- and Ligand-Based Approach Followed by Molecular Dynamics Investigation of Ligand Binding Mode.

An intense effort is made by pharmaceutical and academic research laboratories to identify and develop selective antagonists for each adenosine receptor (AR) subtype as potential clinical candidates for "soft" treatment of various diseases. Crystal structures of subtypes A2A and A1ARs offer exciting opportunities for structure-based drug design. In the first part of the present work, Maybridge HitFinder library of 14400 compounds was utilized to apply a combination of structure-based against the crystal structure of A2AAR and ligand-based methodologies. The docking poses were rescored by CHARMM energy minimization and calculation of the desolvation energy using Poisson-Boltzmann equation electrostatics. Out of the eight selected and tested compounds, five were found positive hits (63% success). Although the project was initially focused on targeting A2AAR, the identified antagonists exhibited low micromolar or micromolar affinity against A2A/A3, ARs, or A3AR, respectively. Based on these results, 19 compounds characterized by novel chemotypes were purchased and tested. Sixteen of them were identified as AR antagonists with affinity toward combinations of the AR family isoforms (A2A/A3, A1/A3, A1/A2A/A3, and A3). The second part of this work involves the performance of hundreds of molecular dynamics (MD) simulations of complexes between the ARs and a total of 27 ligands to resolve the binding interactions of the active compounds, which were not achieved by docking calculations alone. This computational work allowed the prediction of stable and unstable complexes which agree with the experimental results of potent and inactive compounds, respectively. Of particular interest is that the 2-amino-thiophene-3-carboxamides, 3-acylamino-5-aryl-thiophene-2-carboxamides, and carbonyloxycarboximidamide derivatives were found to be selective and possess a micromolar to low micromolar affinity for the A3 receptor.