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

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

The main protease (M^pro) of SARS-CoV-2 is a key antiviral drug target. While most M^pro inhibitors have a γ-lactam glutamine surrogate at the P1 position, we recently found that several M^pro 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 M^pro in complex with calpain inhibitors II and XII and three analogs of GC-376 The structure of M^pro 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.

Chemical Probes for Blocking of Influenza A M2 Wild-type and S31N Channels

We report on using the synthetic aminoadamantane-CH₂-aryl derivatives 1-6 as sensitive probes for blocking M2 S31N and influenza A virus (IAV) M2 wild-type (WT) channels as well as virus replication in cell culture. The binding kinetics measured using electrophysiology (EP) for M2 S31N channel are very dependent on the length between the adamantane moiety and the first ring of the aryl headgroup realized in 2 and 3 and the girth and length of the adamantane adduct realized in 4 and 5. Study of 1-6 shows that, according to molecular dynamics (MD) simulations and molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) calculations, all bind in the M2 S31N channel with the adamantyl group positioned between V27 and G34 and the aryl group projecting out of the channel with the phenyl (or isoxazole in 6) embedded in the V27 cluster. In this outward binding configuration, an elongation of the ligand by only one methylene in rimantadine 2 or using diamantane or triamantane instead of adamantane in 4 and 5, respectively, causes incomplete entry and facilitates exit, abolishing effective block compared to the amantadine derivatives 1 and 6. In the active M2 S31N blockers 1 and 6, the phenyl and isoxazolyl head groups achieve a deeper binding position and high k_on/low k_off and high k_on/high k_off rate constants, compared to inactive 2-5, which have much lower k_on and higher k_off. Compounds 1-5 block the M2 WT channel by binding in the longer area from V27-H37, in the inward orientation, with high k_on and low k_off rate constants. Infection of cell cultures by influenza virus containing M2 WT or M2 S31N is inhibited by 1-5 or 1-4 and 6, respectively. While 1 and 6 block infection through the M2 block mechanism in the S31N variant, 2-4 may block M2 S31N virus replication in cell culture through the lysosomotropic effect, just as chloroquine is thought to inhibit SARS-CoV-2 infection.

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

The main protease (Mpro) of SARS-CoV-2, the pathogen responsible for the COVID-19 pandemic, is a key antiviral drug target. While most SARS-CoV-2 Mpro inhibitors have a γ-lactam glutamine surrogate at the P1 position, we recently discovered several Mpro inhibitors have hydrophobic moieties at the P1 site, including calpain inhibitors II/XII, which are also active against human cathepsin L, a host-protease that is important for viral entry. To determine the binding mode of these calpain inhibitors and establish a structure-activity relationship, we solved X-ray crystal structures of Mpro in complex with calpain inhibitors II and XII, and three analogues of GC-376, one of the most potent Mpro inhibitors in vitro. The structure of Mpro with calpain inhibitor II confirmed the S1 pocket of Mpro can accommodate a hydrophobic methionine side chain, challenging the idea that a hydrophilic residue is necessary at this position. Interestingly, the structure of calpain inhibitor XII revealed an unexpected, inverted binding pose where the P1' pyridine inserts in the S1 pocket and the P1 norvaline is positioned in the S1' pocket. The overall conformation is semi-helical, wrapping around the catalytic core, in contrast to the extended conformation of other peptidomimetic inhibitors. Additionally, the structures of three GC-376 analogues UAWJ246, UAWJ247, and UAWJ248 provide insight to the sidechain preference of the S1', S2, S3 and S4 pockets, and the superior cell-based activity of the aldehyde warhead compared with the α-ketoamide. Taken together, the biochemical, computational, structural, and cellular data presented herein provide new directions for the development of Mpro inhibitors as SARS-CoV-2 antivirals.

Investigation of the Drug Resistance Mechanism of M2-S31N Channel Blockers through Biomolecular Simulations and Viral Passage Experiments

Recent efforts in drug development against influenza A virus (IAV) M2 proton channel S31N mutant resulted in conjugates of amantadine linked with aryl head heterocycles. To understand the mechanism of drug resistance, we chose a representative M2-S31N inhibitor, compound 3, as a chemical probe to identify resistant mutants. To increase the possibility of identifying novel resistant mutants, serial viral passage experiments were performed with multiple strains of H1N1 and H3N2 viruses in different cell lines. This approach not only identified M2 mutations around the drug-binding site, including the pore-lining residues (V27A, V27F, N31S, and G34E) and an interhelical residue (I32N), but also a new allosteric mutation (R45H), in addition to L46P previously identified, located at the C-terminus of M2 that is more than 10 Å away from the drug-binding site. The effects of each mutation were next investigated using electrophysiology, recombinant viruses, and molecular dynamics (MD) simulations. The reduced sensitivity in channel blockage correlated with increased drug resistance in antiviral assays using recombinant viruses. The MD simulations show that the V27A, V27F, G34E, and R45H mutations increase the diameter and hydration state of the pore in complex with compound 3. The Molecular Mechanics Generalized Born (MM-GBSA) calculations result in more positive binding free energies for the complexes of resistant M2 (V27A, V27F, G34E, R45H) with compound 3 compared to the stable complexes (S31N and I32N). Overall, this is the first systematic study of the drug resistance mechanism of M2-S31N channel blockers using multiple viruses in different cell lines.

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 A₃ receptor (A₃R) 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 A₃R. 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 A₃R. MD simulations and MM-GBSA calculations of the WT A₃R 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 A₃R is defined not only by direct interactions with residues within the orthosteric binding area but also by remote residues playing a significant role. Although V169^5.30 is considered to be a selectivity filter for A₃R 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 L90^3.32 in the low region and the remote L264^7.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 A₃R, which we previously performed for agonists binding to A₃R, 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 A₃ receptor (A₃R) binds adenosine and is a drug target against cancer cell proliferation. Currently, there is no experimental structure of A₃R. Here, we have generated a molecular model of A₃R 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 A₃R 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 A₃R requires not only direct interactions with residues within the orthosteric binding area but also with remote residues. Although V169^5.30 is considered to be a selectivity filter for A₃R 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 A₃R agonists.

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 A_2A and A₁ARs 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 A_2AAR 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 A_2AAR, the identified antagonists exhibited low micromolar or micromolar affinity against A_2A/A₃, ARs, or A₃AR, 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 (A_2A/A₃, A₁/A₃, A₁/A_2A/A₃, and A₃). 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 A₃ receptor.

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