Structure-based Inhibitor Design Research Articles

Synthesis, In Silico and Kinetics Evaluation of N-(β-d-glucopyranosyl)-2-arylimidazole-4(5)-carboxamides and N-(β-d-glucopyranosyl)-4(5)-arylimidazole-2-carboxamides as Glycogen Phosphorylase Inhibitors.

Recently studied N-(β-d-glucopyranosyl)-3-aryl-1,2,4-triazole-5-carboxamides have proven to be low micromolar inhibitors of glycogen phosphorylase (GP), a validated target for the treatment of type 2 diabetes mellitus. Since in other settings, the bioisosteric replacement of the 1,2,4-triazole moiety with imidazole resulted in significantly more efficient GP inhibitors, in silico calculations using Glide molecular docking along with unbound state DFT calculations were performed on N-(β-d-glucopyranosyl)-arylimidazole-carboxamides, revealing their potential for strong GP inhibition. The syntheses of the target compounds involved the formation of an amide bond between per-O-acetylated β-d-glucopyranosylamine and the corresponding arylimidazole-carboxylic acids. Kinetics experiments on rabbit muscle GPb revealed low micromolar inhibitors, with the best inhibition constants (Kis) of ~3-4 µM obtained for 1- and 2-naphthyl-substituted N-(β-d-glucopyranosyl)-imidazolecarboxamides, 2b-c. The predicted protein-ligand interactions responsible for the observed potencies are discussed and will facilitate the structure-based design of other inhibitors targeting this important therapeutic target. Meanwhile, the importance of the careful consideration of ligand tautomeric states in binding calculations is highlighted, with the usefulness of DFT calculations in this regard proposed.

Structure-based design of pan-coronavirus inhibitors targeting host cathepsin L and calpain-1

Respiratory disease caused by coronavirus infection remains a global health crisis. Although several SARS-CoV-2-specific vaccines and direct-acting antivirals are available, their efficacy on emerging coronaviruses in the future, including SARS-CoV-2 variants, might be compromised. Host-targeting antivirals provide preventive and therapeutic strategies to overcome resistance and manage future outbreak of emerging coronaviruses. Cathepsin L (CTSL) and calpain-1 (CAPN1) are host cysteine proteases which play crucial roles in coronaviral entrance into cells and infection-related immune response. Here, two peptidomimetic α-ketoamide compounds, 14a and 14b, were identified as potent dual target inhibitors against CTSL and CAPN1. The X-ray crystal structures of human CTSL and CAPN1 in complex with 14a and 14b revealed the covalent binding of α-ketoamide groups of 14a and 14b to C25 of CTSL and C115 of CAPN1. Both showed potent and broad-spectrum anticoronaviral activities in vitro, and it is worth noting that they exhibited low nanomolar potency against SARS-CoV-2 and its variants of concern (VOCs) with EC50 values ranging from 0.80 to 161.7 nM in various cells. Preliminary mechanistic exploration indicated that they exhibited anticoronaviral activity through blocking viral entrance. Moreover, 14a and 14b exhibited good oral pharmacokinetic properties in mice, rats and dogs, and favorable safety in mice. In addition, both 14a and 14b treatments demonstrated potent antiviral potency against SARS-CoV-2 XBB 1.16 variant infection in a K18-hACE2 transgenic mouse model. And 14b also showed effective antiviral activity against HCoV-OC43 infection in a mouse model with a final survival rate of 60%. Further evaluation showed that 14a and 14b exhibited excellent anti-inflammatory effects in Raw 264.7 mouse macrophages and in mice with acute pneumonia. Taken together, these results suggested that 14a and 14b are promising drug candidates, providing novel insight into developing pan-coronavirus inhibitors with antiviral and anti-inflammatory properties.

Structure-based design of small molecule inhibitors of the cagT4SS ATPase Cagα of Helicobacter pylori.

We here describe the structure-based design of small molecule inhibitors of the type IV secretion system of Helicobacter pylori. The secretion system is encoded by the cag pathogenicity island, and we chose Cagα, a hexameric ATPase and member of the family of VirB11-like proteins, as target for inhibitor design. We first solved the crystal structure of Cagα in a complex with the previously identified small molecule inhibitor 1G2. The molecule binds at the interface between two Cagα subunits and mutagenesis of the binding site identified Cagα residues F39 and R73 as critical for 1G2 binding. Based on the inhibitor binding site we synthesized 98 small molecule derivates of 1G2 to improve binding of the inhibitor. We used the production of interleukin-8 of gastric cancer cells during H. pylori infection to screen the potency of inhibitors and we identified five molecules (1G2_1313, 1G2_1338, 1G2_2886, 1G2_2889, and 1G2_2902) that have similar or higher potency than 1G2. Differential scanning fluorimetry suggested that these five molecules bind Cagα, and enzyme assays demonstrated that some are more potent ATPase inhibitors than 1G2. Finally, scanning electron microscopy revealed that 1G2 and its derivatives inhibit the assembly of T4SS-determined extracellular pili suggesting a mechanism for their anti-virulence effect.

Ion Transport and Inhibitor Binding by Human NHE1: Insights from Molecular Dynamics Simulations and Free Energy Calculations.

The human Na+/H+ exchanger (NHE1) plays a crucial role in maintaining intracellular pH by regulating the electroneutral exchange of a single intracellular H+ for one extracellular Na+ across the plasma membrane. Understanding the molecular mechanisms governing ion transport and the binding of inhibitors is of importance in the development of anticancer therapeutics targeting NHE1. In this context, we performed molecular dynamics (MD) simulations based on the recent cryo-electron microscopy (cryo-EM) structures of outward- and inward-facing conformations of NHE1. These simulations allowed us to explore the dynamics of the protein, examine the ion-translocation pore, and confirm that Asp267 is the ion-binding residue. Our free energy calculations did not show a significant difference between Na+ and K+ binding at the ion-binding site. Consequently, Na+ over K+ selectivity cannot be solely explained by differences in ion binding. Our MD simulations involving NHE1 inhibitors (cariporide and amiloride analogues) maintained stable interactions with Asp267 and Glu346. Our study highlights the importance of the salt bridge between the positively charged acylguanidine moiety and Asp267, which appears to play a role in the competitive inhibitory mechanism for this class of inhibitors. Our computational study provides a detailed mechanistic interpretation of experimental data and serves the basis of future structure-based inhibitor design.

Uncovering New Conformational States of the Substrate Binding Pocket of LSD1 Potential for Inhibitor Design via Funnel Metadynamics.

Lysine-specific demethylase 1 (LSD1) is a promising therapeutic target for cancer therapy. So far, over 80 crystal structures of LSD1 in different complex states have been deposited in the Protein Data Bank, which are valuable resources for performing structure-based drug design. However, among all of the crystal structures of LSD1, the substrate binding pocket, which is the most efficient druggable site for designing LSD1 inhibitors at present, is very similar no matter whether LSD1 is in the apo or any holo forms, which is inconsistent with its versatile demethylase functions. To investigate whether the substrate binding pocket is rigid or exhibits other representative conformations different from the crystal conformations that are feasible for designing new LSD1 inhibitors, we performed funnel metadynamics simulations to study the conformation dynamics of LSD1 in the binding process of two effective LSD1 inhibitors (CC-90011 and 6X0, CC-90011 undergoing clinical trials). Our results showed that the entrance of the substrate binding pocket is very flexible. Two representative entrance conformations of LSD1 counting against binding with the substrate of histone H3 were detected, which may be used for structure-based LSD1 inhibitor design. Besides, alternative optimal binding modes and prebinding modes for both inhibitors were also detected, which depicted that the key interactions changed along with the binding process. Our results should provide great help for LSD1 inhibitor design.

Structure-Based Design and Evaluation of Reversible KRAS G13D Inhibitors.

Oncogenic KRAS mutations were identified decades ago, yet the selective inhibition of specific KRAS mutant proteins represents an ongoing challenge. Recent progress has been made in targeting certain P-loop mutant proteins, in particular KRAS G12C, for which the covalent inhibition of the GDP state via the Switch II pocket is now a clinically validated strategy. Inhibition of other KRAS mutant proteins such as KRAS G13D, on the other hand, still requires clinical validation. The remoteness of the D13 residue relative to the Switch II pocket in combination with the solvent exposure and conformational flexibility of the D13 side chain, as well as the difficulties of targeting carboxylate residues covalently, renders this specific protein particularly challenging to target selectively. In this report, we describe the design and evaluation of potent and KRAS G13D-selective reversible inhibitors. Subnanomolar binding to the GDP state Switch II pocket and biochemical selectivity over WT KRAS are achieved by leveraging a salt bridge with D13.

Abstract B090: Partially open conformation of the G323E mutated HIF-2α PASB domain captured by X-ray crystallography

Abstract Under hypoxic conditions, vital cell processes such as angiogenesis, erythropoiesis, cell proliferation and metabolism are regulated by a group of transcription factors known as hypoxia-inducible factors (HIFs). These proteins belong to the helix-loop-helix family of transcription factors and require heterodimerization to function. Constitutively expressed HIF-1β protein forms a heterodimer with one of the oxygen-labile counterparts: HIF-1α, HIF-2α, or HIF-3α. The HIF-αs are degraded under normoxic conditions via ubiquitination by the von Hippel Lindau (VHL) E3 ubiquitin ligase and become available for dimerization under either hypoxia or VHL dysfunction. Certain cancer types take advantage of this HIF-regulated transcription for survival and growth, as exemplified by VHL-deficient clear cell renal cell carcinoma (ccRCC). Inhibitors aimed at disruption of the HIF-1β:HIF-2α heterodimer have been advanced into clinical development, including belzutifan (approved by the Food and Drug Administration for adult patients with VHL disease who require therapy for RCC among other indications) and AB521 (Arcus Biosciences). This class of small molecules functions by binding to and causing a conformational change in the internal cavity of the PasB domain of HIF-2α. A mutation within this internal cavity, G323E, was reported in a patient undergoing treatment with HIF-2α inhibitor. This mutation restored the ability of HIF-1β to heterodimerize with HIF-2α in presence of the inhibitor. To better understand this potential resistance mechanism and how the G323E mutation affects inhibitor binding, we used X-ray crystallography to examine the structure of the mutant HIF-2α PasB domain in complex with a novel small molecule inhibitor active in vitro against the G323E mutant form. We report a co-crystal structure that reveals a partially disordered open conformation of the HIF-2α PasB domain. This structure strengthens existing hypotheses regarding the route of compound entry into the HIF-2α PasB cavity and lays the foundation for structure-based design of HIF-2α PasB G323E-active inhibitors. Citation Format: Steven Shia, Ian Malgapo, Jeffrey Chen, Amber Pham, Cesar A Meleza, Lunda Shen, Mariya Morar, Kenneth V Lawson, Kelsey E Sivick Gauthier, Hyock J Kwon. Partially open conformation of the G323E mutated HIF-2α PASB domain captured by X-ray crystallography [abstract]. In: Proceedings of the AACR-NCI-EORTC Virtual International Conference on Molecular Targets and Cancer Therapeutics; 2023 Oct 11-15; Boston, MA. Philadelphia (PA): AACR; Mol Cancer Ther 2023;22(12 Suppl):Abstract nr B090.

Structure-based design of potent FABP4 inhibitors with high selectivity against FABP3

Fatty-acid binding protein 4 (FABP4) presents an attractive target for therapeutic intervention in metabolic and inflammatory diseases in recent years. However, highly similar three-dimensional structures and fatty acid binding ability of multiple FABP family members pose a significant challenge in design of FABP4-selective inhibitors. Particularly, inhibition of FABP3 raises safety concerns such as cardiac dysfunction and exercise intolerance. Here, we reported the discovery of new FABP4 inhibitors with high selectivity over FABP3 by exploiting the little structural difference in the ligand binding pockets of FABP4 and FABP3. On the basis of our previously reported FABP4 inhibitors with nanomolar potency, different substituents were further introduced to perfectly occupy two sub-pockets of FABP4 that are distinct from those of FABP3. Remarkably, a single methyl group introduction leads to the discovery of compound C3 that impressively exhibits a 601-fold selectivity over FABP3 when maintained nanomolar binding affinity for FABP4. Moreover, C3 also shows good metabolic stability and potent cellular anti-inflammatory activity, making it a promising inhibitor for further development. Therefore, the present study highlights the utility of the structure-based rational design strategy for seeking highly selective and potent inhibitors of FABP4 and the importance of identifying the appropriate subsite as well as substituent for gaining the desired selectivity.

GaMF1.39's antibiotic efficacy and its enhanced antitubercular activity in combination with clofazimine, Telacebec, ND-011992, or TBAJ-876.

New drugs are needed to combat multidrug-resistant tuberculosis. The electron transport chain (ETC) maintains the electrochemical potential across the cytoplasmic membrane and allows the production of ATP, the energy currency of any living cell. The mycobacterial engine F-ATP synthase catalyzes the formation of ATP and has come into focus as an attractive and rich drug target. Recent deep insights into these mycobacterial F1FO-ATP synthase elements opened the door for a renaissance of structure-based target identification and inhibitor design. In this study, we present the GaMF1.39 antimycobacterial compound, targeting the rotary subunit γ of the biological engine. The compound is bactericidal, inhibits infection ex vivo, and displays enhanced anti-tuberculosis activity in combination with ETC inhibitors, which promises new strategies to shorten tuberculosis chemotherapy.

Structure-Based Design and Synthesis of Potent and Selective KRAS G12D Inhibitors.

KRAS G12D mutation has been found in approximately 45% of pancreatic ductal adenocarcinoma (PDAC) cases, making it an attractive therapeutic target. Through structure-based drug design, a series of potent and selective KRAS G12D inhibitors were designed. The lead compound, ERAS-5024, inhibited ERK1/2 phosphorylation and cell proliferation in three-dimensional Cell-Titer Glo assays in AsPC-1 PDAC cells with single-digit nanomolar potency and caused tumor regression in the in vivo efficacy studies. We describe here the details of the design and synthesis program that led to the discovery of ERAS-5024.

Molecular basis for bacterial N-glycosylation by a soluble HMW1C-like N-glycosyltransferase

Soluble HMW1C-like N-glycosyltransferases (NGTs) catalyze the glycosylation of Asn residues in proteins, a process fundamental for bacterial autoaggregation, adhesion and pathogenicity. However, our understanding of their molecular mechanisms is hindered by the lack of structures of enzymatic complexes. Here, we report structures of binary and ternary NGT complexes of Aggregatibacter aphrophilus NGT (AaNGT), revealing an essential dyad of basic/acidic residues located in the N-terminal all α-domain (AAD) that intimately recognizes the Thr residue within the conserved motif Asn0-X+1-Ser/Thr+2. Poor substrates and inhibitors such as UDP-galactose and UDP-glucose mimetics adopt non-productive conformations, decreasing or impeding catalysis. QM/MM simulations rationalize these results, showing that AaNGT follows a SN2 reaction mechanism in which the acceptor asparagine uses its imidic form for catalysis and the UDP-glucose phosphate group acts as a general base. These findings provide key insights into the mechanism of NGTs and will facilitate the design of structure-based inhibitors to treat diseases caused by non-typeable H. influenzae or other Gram-negative bacteria.

Structural basis of lipid-droplet localization of 17-beta-hydroxysteroid dehydrogenase 13

Hydroxysteroid 17-beta-dehydrogenase 13 (HSD17B13) is a hepatic lipid droplet-associated enzyme that is upregulated in patients with non-alcoholic fatty liver disease. Recently, there have been several reports that predicted loss of function variants in HSD17B13 protect against the progression of steatosis to non-alcoholic steatohepatitis with fibrosis and hepatocellular carcinoma. Here we report crystal structures of full length HSD17B13 in complex with its NAD+ cofactor, and with lipid/detergent molecules and small molecule inhibitors from two distinct series in the ligand binding pocket. These structures provide insights into a mechanism for lipid droplet-associated proteins anchoring to membranes as well as a basis for HSD17B13 variants disrupting function. Two series of inhibitors interact with the active site residues and the bound cofactor similarly, yet they occupy different paths leading to the active site. These structures provide ideas for structure-based design of inhibitors that may be used in the treatment of liver disease.

Molecular cloning, overexpression, biophysical studies and structure-based inhibitor design for the Raetz pathway peripheral membrane protein UDP-diacylglucosamine pyrophosphohydrolase (LpxH) from Salmonella typhi

Molecular cloning, overexpression, biophysical studies and structure-based inhibitor design for the Raetz pathway peripheral membrane protein UDP-diacylglucosamine pyrophosphohydrolase (LpxH) from <i>Salmonella typhi</i>

Structure-Based Design of Transport-Specific Multitargeted One-Carbon Metabolism Inhibitors in Cytosol and Mitochondria.

Multitargeted agents provide tumor selectivity with reduced drug resistance and dose-limiting toxicities. We previously described the multitargeted 6-substituted pyrrolo[3,2-d]pyrimidine antifolate 1 with activity against early- and late-stage pancreatic tumors with limited tumor selectivity. Structure-based design with our human serine hydroxymethyl transferase (SHMT) 2 and glycinamide ribonucleotide formyltransferase (GARFTase) structures, and published X-ray crystal structures of 5-aminoimidazole-4-carboxamide ribonucleotide formyltransferase/inosine monophosphate cyclohydrolase (ATIC), SHMT1, and folate receptor (FR) α and β afforded 11 analogues. Multitargeted inhibition and selective tumor transport were designed by providing promiscuous conformational flexibility in the molecules. Metabolite rescue identified mitochondrial C1 metabolism along with de novo purine biosynthesis as the targeted pathways. We identified analogues with tumor-selective transport via FRs and increased SHMT2, SHMT1, and GARFTase inhibition (28-, 21-, and 11-fold, respectively) compared to 1. These multitargeted agents represent an exciting new structural motif for targeted cancer therapy with substantial advantages of selectivity and potency over clinically used antifolates.

Structural basis for the development of potential inhibitors targeting FadD23 from Mycobacterium tuberculosis.

Sulfolipid-1 (SL-1) is a lipid that is abundantly found in the cell wall of Mycobacterium tuberculosis (Mtb). MtbFadD23 is crucial in the SL-1 synthesis pathway. Previously, 5'-O-[N-(11-phenoxyundecanoyl)sulfamoyl]adenosine (PhU-AMS) has been shown to be a general inhibitor of fatty-acid-adenylating enzymes (FadDs) in Mtb. However, the fatty acyl-AMP ligase (FAAL) class of FadDs, which includes MtbFadD23, appears to be functionally nonredundant in the production of multiple fatty acids. In this study, the ability of PhU-AMS to bind to MtbFadD23 was examined under in vitro conditions. The crystal structure of the MtbFadD23-PhU-AMS complex was determined at a resolution of 2.64 Å. Novel features were identified by structural analysis and comparison. Although PhU-AMS could bind to MtbFadD23, it did not inhibit the FAAL adenylation activity of MtbFadD23. However, PhU-AMS improved the main Tm value in a differential scanning fluorimetry assay, and a structural comparison of MtbFadD23-PhU-AMS with FadD32 and PA1221 suggested that PhU-AMS blocks the loading of the acyl chain onto Pks2. This study sheds light on the structure-based design of specific inhibitors of MtbFadD23 and general inhibitors of FAALs.

Structure-Based Design of Novel MAO-B Inhibitors: A Review.

With the significant growth of patients suffering from neurodegenerative diseases (NDs), novel classes of compounds targeting monoamine oxidase type B (MAO-B) are promptly emerging as distinguished structures for the treatment of the latter. As a promising function of computer-aided drug design (CADD), structure-based virtual screening (SBVS) is being heavily applied in processes of drug discovery and development. The utilization of molecular docking, as a helping tool for SBVS, is providing essential data about the poses and the occurring interactions between ligands and target molecules. The current work presents a brief discussion of the role of MAOs in the treatment of NDs, insight into the advantages and drawbacks of docking simulations and docking software, and a look into the active sites of MAO-A and MAO-B and their main characteristics. Thereafter, we report new chemical classes of MAO-B inhibitors and the essential fragments required for stable interactions focusing mainly on papers published in the last five years. The reviewed cases are separated into several chemically distinct groups. Moreover, a modest table for rapid revision of the revised works including the structures of the reported inhibitors together with the utilized docking software and the PDB codes of the crystal targets applied in each study is provided. Our work could be beneficial for further investigations in the search for novel, effective, and selective MAO-B inhibitors.

Structure-Based Drug Design of RdRp Inhibitors against SARS-CoV-2.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has caused a worldwide pandemic since 2019, spreading rapidly and posing a significant threat to human health and life. With over 6 billion confirmed cases of the virus, the need for effective therapeutic drugs has become more urgent than ever before. RNA-dependent RNA polymerase (RdRp) is crucial in viral replication and transcription, catalysing viral RNA synthesis and serving as a promising therapeutic target for developing antiviral drugs. In this article, we explore the inhibition of RdRp as a potential treatment for viral diseases, analysing the structural information of RdRp in virus proliferation and summarizing the reported inhibitors' pharmacophore features and structure-activity relationship profiles. We hope that the information provided by this review will aid in structure-based drug design and aid in the global fight against SARS-CoV-2 infection.

Structure-Based Design of Inhibitors of the m6A-RNA Writer Enzyme METTL3.

Methyltransferase-like 3 (METTL3) and METTL14 form a heterodimeric complex that catalyzes the most abundant internal mRNA modification, N6-methyladenosine (m6A). METTL3 is the catalytic subunit that binds the co-substrate S-adenosyl methionine (SAM), while METTL14 is involved in mRNA binding. The m6A modification provides post-transcriptional level control over gene expression as it affects almost all stages of the mRNA life cycle, including splicing, nuclear export, translation, and decay. There is increasing evidence for an oncogenic role of METTL3 in acute myeloid leukemia. Here, we use structural and dynamic details of the catalytic subunit METTL3 for developing small-molecule inhibitors that compete with SAM. Starting from a hit identified by high-throughput docking, protein crystallography and molecular dynamics simulations were employed to guide the optimization of inhibitory activity. The potency was successfully improved by 8000-fold as measured by a homogeneous time-resolved fluorescence assay. The optimized compound is selective against the off-targets RNA methyltransferases METTL1 and METTL16.

Discovery and Structure-Based Design of Inhibitors of the WD Repeat-Containing Protein 5 (WDR5)-MYC Interaction.

WDR5 is a critical chromatin cofactor of MYC. WDR5 interacts with MYC through the WBM pocket and is hypothesized to anchor MYC to chromatin through its WIN site. Blocking the interaction of WDR5 and MYC impairs the recruitment of MYC to its target genes and disrupts the oncogenic function of MYC in cancer development, thus providing a promising strategy for the treatment of MYC-dysregulated cancers. Here, we describe the discovery of novel WDR5 WBM pocket antagonists containing a 1-phenyl dihydropyridazinone 3-carboxamide core that was identified from high-throughput screening and subsequent structure-based design. The leading compounds showed sub-micromolar inhibition in the biochemical assay. Among them, compound 12 can disrupt WDR5-MYC interaction in cells and reduce MYC target gene expression. Our work provides useful probes to study WDR5-MYC interaction and its function in cancers, which can also be used as the starting point for further optimization toward drug-like small molecules.

Computational analysis of hypothetical proteins from Mycobacterium orygis identifies proteins with therapeutic and diagnostic potentials

The genome of Mycobacterium orygis encodes 3767 proteins that includes 442 functionally uncharacterized proteins, known as hypothetical proteins (HPs). These functionally unknown proteins may play significant biological role in bacterium's survivability, virulence etc., hence functional annotation of these HPs with bioinformatics techniques might shed light on these aspects. So, the present study aimed to assign functions to those HPs in M. orygis genome. With the help of an extensive in silico workflow involving several bioinformatic tools, physicochemical characteristics, protein-protein interaction, sub-cellular localization, functional classification, and antigenicity of the HPs were analyzed. Out of 4039 genes and 3767 encoded proteins, 422 were annotated as “hypothetical protein, unknown function” in NCBI database. Out of 422 HPs, the FASTA sequences of 400 proteins, available in the NCBI database, were extracted. Preliminarily 400 proteins were screened for presence of conserved domain by different bioinformatic tools and 32 shortlisted proteins were further analyzed. The analysis revealed that these proteins may be involved in enzymatic, toxin-antitoxin system, cellular and other processes. Virulence and antigenicity of the complete set of HPs were also predicted. Based on antigenicity scores and sub-cellular localization, four proteins were selected for secondary and tertiary structure modelling. Molecular docking and MD-simulation of four proteins with TLR1 showed stable binding. B-cell and T-cell epitopes were identified for two proteins (QOY47331.1 and QOY49361.1) with highest antigenicity scores. The result of the present study may facilitate better understanding of virulence, pathogenesis and therapeutic target. The information about the structures of the HPs with highest antigenicity score could aid in structure-based inhibitor design and devising novel diagnostic assays.