Experimental and Theoretical Evaluation of the Ethynyl Moiety as a Halogen Bioisostere

First fragment-based screening identifies new chemotypes inhibiting ERAP1-metalloprotease

Inhibition of endoplasmic reticulum aminopeptidase 1 (ERAP1) by small-molecules is being eagerly investigated for the treatment of various autoimmune diseases and in the field of immuno-oncology after its active involvement in antigen presentation and processing. Currently, ERAP1 inhibitors are at different stages of clinical development, which highlights its significance as a promising drug target. In the present work, we describe the first-ever successful identification of several ERAP1 inhibitors derived from a fragment-based screening approach. We applied an enzymatic activity assay to a large library of ∼3000 fragment entries in order to retrieve 32 hits. After a multi-faceted selection process, we prioritized 3 chemotypes for SAR optimization and strategic modifications provided 2 series (2-thienylacetic acid and rhodanine scaffolds) with improved analogues at the low micromolar range of ERAP1 inhibition. We report also evidence of selectivity against homologous aminopeptidase IRAP, combined with complementary in silico docking studies to predict the binding mode and site of inhibition. Our compounds can be the starting point for future fragment growing and rational drug development, incorporating new chemical modalities.

Computational Design, Synthesis, and Biological Evaluation of Diimidazole Analogues Endowed with Dual PCSK9/HMG-CoAR-Inhibiting Activity

Proprotein convertase subtilisin/kexin 9 (PCSK9) is responsible for the degradation of the hepatic low-density lipoprotein receptor (LDLR), which regulates circulating cholesterol levels. Consequently, the PCSK9 inhibition is a valuable therapeutic approach for the treatment of hypercholesterolemia and cardiovascular diseases. In our studies, we discovered Rim13, a polyimidazole derivative reducing the protein-protein interaction between PCSK9 and LDLR with an IC₅₀ of 1.6 μM. The computational design led to the optimization of the shape of the PCSK9/ligand complementarity, enabling the discovery of potent diimidazole derivatives. In fact, carrying out biological assays to fully characterize the cholesterol-lowering activity of the new analogues and using both biochemical and cellular techniques, compound Dim16 displayed improved PCSK9 inhibitory activity (IC₅₀ 0.9 nM). Interestingly, similar to other lupin-derived peptides and their synthetic analogues, some compounds in this series showed dual hypocholesterolemic activity since some of them complementarily inhibited the 3-hydroxy-3-methylglutaryl coenzyme A reductase.

Exploring the Antitubercular Activity of Anthranilic Acid Derivatives: From MabA (FabG1) Inhibition to Intrabacterial Acidification

Mycobacterium tuberculosis, the pathogen that causes tuberculosis, is responsible for the death of 1.5 million people each year and the number of bacteria resistant to the standard regimen is constantly increasing. This highlights the need to discover molecules that act on new M. tuberculosis targets. Mycolic acids, which are very long-chain fatty acids essential for M. tuberculosis viability, are synthesized by two types of fatty acid synthase (FAS) systems. MabA (FabG1) is an essential enzyme belonging to the FAS-II cycle. We have recently reported the discovery of anthranilic acids as MabA inhibitors. Here, the structure–activity relationships around the anthranilic acid core, the binding of a fluorinated analog to MabA by NMR experiments, the physico-chemical properties and the antimycobacterial activity of these inhibitors were explored. Further investigation of the mechanism of action in bacterio showed that these compounds affect other targets than MabA in mycobacterial cells and that their antituberculous activity is due to the carboxylic acid moiety which induces intrabacterial acidification.

Synthesis of Amino Acids Bearing Halodifluoromethyl Moieties and Their Application to p53-Derived Peptides Binding to Mdm2/Mdm4

Introduction

Therapeutic peptides are a significant class of drugs in the treatment of a wide range of diseases. To enhance their properties, such as stability or binding affinity, they are usually chemically modified. This includes, among other techniques, cyclization of the peptide chain by bridging, modifications to the backbone, and incorporation of unnatural amino acids. One approach previously established, is the use of halogenated aromatic amino acids. In principle, they are thereby enabled to form halogen bonds (XB). In this study, we focus on the -R-CF₂X moiety (R = O, NHCO; X = Cl, Br) as an uncommon halogen bond donor. These groups enable more spatial variability in protein-protein interactions. The chosen approach via Fmoc-protected building blocks allows for the incorporation of these modified amino acids in peptides using solid-phase peptide synthesis.

Results and Discussion

Using a competitive fluorescence polarization assay to monitor binding to Mdm4, we demonstrate that a p53-derived peptide with Lys24Nle(εNHCOCF₂X) exhibits an improved inhibition constant K_i compared to the unmodified peptide. Decreasing K_i values observed with the increasing XB capacity of the halogen atoms (F ≪ Cl < Br) indicates the formation of a halogen bond. By reducing the side chain length of Nle(εNHCOCF₂X) to Abu(γNHCOCF₂X) as control experiments and through quantum mechanical calculations, we suggest that the observed affinity enhancement is related to halogen bond-induced intramolecular stabilization of the α-helical binding mode of the peptide or a direct interaction with His54 in human Mdm4.

Revisiting a challenging p53 binding site: a diversity-optimized HEFLib reveals diverse binding modes in T-p53C-Y220C

The cellular tumor antigen p53 is a key component in cell cycle control. The mutation Y220C heavily destabilizes the protein thermally but yields a druggable crevice. We have screened the diversity-optimized halogen-enriched fragment library against T-p53C-Y220C with STD-NMR and DSF to identify hits, which we validated by 1H,15N-HSQC NMR. We could identify four hits binding in the Y220C cleft, one hit binding covalently and four hits binding to an uncharacterized binding site. Compound 1151 could be crystallized showing a flip of C220 and thus opening subsite 3. Additionally, 4482 was identified to alkylate cysteines. Data shows that the diversity-optimized HEFLib leads to multiple diverse hits. The identified scaffolds can be used to further optimize interactions with T-p53C-Y220C and increase thermal stability.

Predicting Blood–Brain Barrier Permeation of Erlotinib and JCN037 by Molecular Simulation

Glioblastoma (GBM) is a highly malignant primary brain tumor, and epidermal growth factor receptor (EGFR) is a well characterized biomaker on GBM. Treatment of GBM with EGFR inhibitors achieved limited efficacy due to low blood-brain barrier (BBB) permeability, and BBB-penetrant drugs are required. In this study, the BBB penetration of erlotinib and JN037 were studied using molecular dynamics method with explicit membrane model. The free energy profiles indicate that JCN037 has a lower central energy barrier than erlotinib, and it has a local minimum at lipid-water interface while erlotinib has not. Unconstrained MD simulations found that erlotinib prefers staying in water while JCN037 tends to interact with lipid molecules. Further analysis reveals that the Br atom of JCN037 plays an important role in its interaction with lipid molecules, and the adjacent F atom enhances the interaction of Br. The two flexible methoxyethoxy chains of erlotinib are responsible for its poor penetration. Our computational results agree well with the experimental results, providing useful information in the design and improvement of drugs with good BBB permeation.

Probing Factor Xa Protein-Ligand Interactions: Accurate Free Energy Calculations and Experimental Validations of Two Series of High-Affinity Ligands

The accurate prediction of protein-ligand binding affinity belongs to one of the central goals in computer-based drug design. Molecular dynamics (MD)-based free energy calculations have become increasingly popular in this respect due to their accuracy and solid theoretical basis. Here, we present a combined study which encompasses experimental and computational studies on two series of factor Xa ligands, which enclose a broad chemical space including large modifications of the central scaffold. Using this integrated approach, we identified several new ligands with different heterocyclic scaffolds different from the previously identified indole-2-carboxamides that show superior or similar affinity. Furthermore, the so far underexplored terminal alkyne moiety proved to be a suitable non-classical bioisosteric replacement for the higher halogen-π aryl interactions. With this challenging example, we demonstrated the ability of the MD-based non-equilibrium free energy calculation approach for guiding crucial modifications in the lead optimization process, such as scaffold replacement and single-site modifications at molecular interaction hot spots.

Rationalizing the binding and α subtype selectivity of synthesized imidazodiazepines and benzodiazepines at GABAA receptors by using molecular docking studies

The pharmacological actions exerted by benzodiazepines are dependent on the discrete α protein subunits of the γ-aminobutyric acid type A receptor (GABA_A R). Recent developments via a cryo-EM structure of the α1β3γ2L GABA_A R ion channel provide crucial insights into ligand efficacy and binding affinity at this subtype. We investigated the molecular interactions of diazepam and alprazolam bound GABA_A R structures (6HUP and 6HUO) to determine key binding interaction domains. A halogen bond between the chlorine atoms of diazepam and alprazolam with the group on the backbone of the α1 histidine amino acid 102 is important to the positive allosteric modulatory actions of diazepam and alprazolam in the α1β3γ2L GABA_A R ion channel. In order to gain insight into α subtype selectivity we designed and synthesized close structural analogs of diazepam and alprazolam. These compounds were then docked into the recently publish cryo-EM structures of GABA_A Rs (6HUP and 6HUO). This modeling along with radio-ligand binding data resulted in the conclusion that the non-classical bioisosteric replacement of the chlorine atom at C7 with an ethinyl group (compound 5) resulted in an 11-fold gain in α5 binding selectivity over the α1 subtype. Moreover, the potency of compound 5 resulted in a ligand with less sedation than diazepam, while still maintaining the same anxiolytic potency. These modeling data extend our understanding of the structural requirements for α-subtype-selective compounds that can be utilized to achieve improved medical treatments. It is clear that the ethinyl group in place of a halogen atom decreases the affinity and efficacy of benzodiazepines and imidazodiazepines at α1 subtypes, which results in less sedation and ataxia.

Identification of MRTX1133, a Noncovalent, Potent, and Selective KRASG12D Inhibitor

KRAS^G12D, the most common oncogenic KRAS mutation, is a promising target for the treatment of solid tumors. However, when compared to KRAS^G12C, selective inhibition of KRAS^G12D presents a significant challenge due to the requirement of inhibitors to bind KRAS^G12D with high enough affinity to obviate the need for covalent interactions with the mutant KRAS protein. Here, we report the discovery and characterization of the first noncovalent, potent, and selective KRAS^G12D inhibitor, MRTX1133, which was discovered through an extensive structure-based activity improvement and shown to be efficacious in a KRAS^G12D mutant xenograft mouse tumor model.

Development of a Potent Brain-Penetrant EGFR Tyrosine Kinase Inhibitor against Malignant Brain Tumors

The epidermal growth factor receptor (EGFR) is genetically altered in nearly 60% of glioblastoma tumors; however, tyrosine kinase inhibitors (TKIs) against EGFR have failed to show efficacy for patients with these lethal brain tumors. This failure is attributed to the inability of clinically tested EGFR TKIs to cross the blood-brain barrier (BBB) and achieve adequate pharmacological levels to inhibit various oncogenic forms of EGFR that drive glioblastoma. Through SAR analysis, we developed compound 5 (JCN037) from an anilinoquinazoline scaffold by ring fusion of the 6,7-dialkoxy groups to reduce the number of rotatable bonds and polar surface area and by introduction of an ortho-fluorine and meta-bromine on the aniline ring for improved potency and BBB penetration. Relative to the conventional EGFR TKIs erlotinib and lapatinib, JCN037 displayed potent activity against EGFR amplified/mutant patient-derived cell cultures, significant BBB penetration (2:1 brain-to-plasma ratio), and superior efficacy in an EGFR-driven orthotopic glioblastoma xenograft model.

A structure-guided molecular chaperone approach for restoring the transcriptional activity of the p53 cancer mutant Y220C

Aim:

The p53 cancer mutation Y220C creates a conformationally unstable protein with a unique elongated surface crevice that can be targeted by molecular chaperones. We report the structure-guided optimization of the carbazole-based stabilizer PK083.

Materials & methods:

Biophysical, cellular and x-ray crystallographic techniques have been employed to elucidate the mode of action of the carbazole scaffolds.

Results:

Targeting an unoccupied subsite of the surface crevice with heterocycle-substituted PK083 analogs resulted in a 70-fold affinity increase to single-digit micromolar levels, increased thermal stability and decreased rate of aggregation of the mutant protein. PK9318, one of the most potent binders, restored p53 signaling in the liver cancer cell line HUH-7 with homozygous Y220C mutation.

Conclusion:

The p53-Y220C mutant is an excellent paradigm for the development of mutant p53 rescue drugs via protein stabilization. Similar rescue strategies may be applicable to other cavity-creating p53 cancer mutations.

Roles of computational modelling in understanding p53 structure, biology, and its therapeutic targeting

The transcription factor p53 plays pivotal roles in numerous biological processes, including the suppression of tumours. The rich availability of biophysical data aimed at understanding its structure–function relationships since the 1990s has enabled the application of a variety of computational modelling techniques towards the establishment of mechanistic models. Together they have provided deep insights into the structure, mechanics, energetics, and dynamics of p53. In parallel, the observation that mutations in p53 or changes in its associated pathways characterize several human cancers has resulted in a race to develop therapeutic modulators of p53, some of which have entered clinical trials. This review describes how computational modelling has played key roles in understanding structural-dynamic aspects of p53, formulating hypotheses about domains that are beyond current experimental investigations, and the development of therapeutic molecules that target the p53 pathway.

Recent Synthetic Approaches towards Small Molecule Reactivators of p53

The tumor suppressor protein p53 is often called "the genome guardian" and controls the cell cycle and the integrity of DNA, as well as other important cellular functions. Its main function is to trigger the process of apoptosis in tumor cells, and approximately 50% of all cancers are related to the inactivation of the p53 protein through mutations in the TP53 gene. Due to the association of mutant p53 with cancer therapy resistance, different forms of restoration of p53 have been subject of intense research in recent years. In this sense, this review focus on the main currently adopted approaches for activation and reactivation of p53 tumor suppressor function, focusing on the synthetic approaches that are involved in the development and preparation of such small molecules.

Acetylene Group, Friend or Foe in Medicinal Chemistry

The use of an acetylene (ethynyl) group in medicinal chemistry coincides with the launch of the Journal of Medicinal Chemistry in 1959. Since then, the acetylene group has been broadly exploited in drug discovery and development. As a result, it has become recognized as a privileged structural feature for targeting a wide range of therapeutic target proteins, including MAO, tyrosine kinases, BACE1, steroid receptors, mGlu5 receptors, FFA1/GPR40, and HIV-1 RT. Furthermore, a terminal alkyne functionality is frequently introduced in chemical biology probes as a click handle to identify molecular targets and to assess target engagement. This Perspective is divided into three parts encompassing: (1) the physicochemical properties of the ethynyl group, (2) the advantages and disadvantages of the ethynyl group in medicinal chemistry, and (3) the impact of the ethynyl group on chemical biology approaches.

Discovery of Small-Molecule Cyclic GMP-AMP Synthase Inhibitors

Cyclic guanosine monophosphate-adenosine monophosphate (GMP-AMP) (cGAS), a cytosolic DNA sensor, plays an important role in the type I interferon response. DNA from either invading microbes or self-origin triggers the enzymatic activity of cGAS. Aberrant activation of cGAS is associated with various autoimmune disorders. Only one selective probe exists for inhibiting cGAS in cells, while others are limited by their poor cellular activity or specificity, which underscores the urgency for discovering new cGAS inhibitors. Here, we describe the development of new small-molecule human cGAS (hcGAS) inhibitors (80 compounds synthesized) with high binding affinity in vitro and cellular activity. Our studies show CU-32 and CU-76 selectively inhibit the DNA pathway in human cells but have no effect on the RIG-I-MAVS or Toll-like receptor pathways. CU-32 and CU-76 represent a new class of hcGAS inhibitors with activity in cells and provide a new chemical scaffold for designing probes to study cGAS function and development of autoimmune therapeutics.

Halogen bonding in halocarbon-protein complexes and computational tools for rational drug design

Introduction: Halogens have a prominent role in drug design. Often used as a mean to improve ADME properties, they are also becoming a tool in protein-ligand recognition given their ability to form a non-covalent interaction, termed halogen bond, where halogens act as electrophilic species interacting with electron-rich partners. Rational drug design of halogen-bonding lead molecules requires an accurate description of halocarbon-protein complexes by computational tools though not all methods are able to tackle this non-covalent interaction. Areas covered: The authors present a review of computational methodologies that can be used to properly describe halogen bonds in the context of protein-ligand complexes, providing also insights on how these methods can be used in the context of computer-aided drug design. Expert opinion: Although in the last few years many computational tools, ranging from fast screening methods to the more expensive QM calculations, have been developed to tackle the halogen bonding phenomenon, they are not yet standard in the literature. This will eventually change as official software distributions are including support for halogen bonding in their methods. Tackling desolvation of halogenated species seems to be a good strategy to improve the accuracy of computational methods, that will be more commonly used prior to laboratory work in the future.

Troxerutin subdues hepatic tumorigenesis via disrupting the MDM2-p53 interaction

Hepatocellular carcinoma (HCC) is the leading cause of cancer death worldwide that lacks proper medical prognosis and treatment. In the present study, the anti-tumoral potential of troxerutin (TX), an ethnomedicine, was examined in relation to its effects on the promoter 2-acetylaminofluorene (2-AAF) in N-nitrosodiethylamine (NDEA) initiated HCC, as compared to its effects on HCC induced by NDEA alone. Liver samples from each experimental group were collected and evaluated for histological, biochemical and cellular characterization. The protein expressions of apoptotic and cell proliferation markers were determined via immunohistochemistry and western blotting. Molecular docking was also performed to delineate the inhibitory mechanism of TX on HCC. The results show that only higher doses of TX showed a significant reduction in the incidence of hepatic nodule formation, and they also counteracted NDEA plus 2-AAF induced alterations in the enzymic status. The frequencies of glutathione-S-transferase and proliferating cell nuclear antigen, markers of S phase progression, were markedly reduced during TX treatment. TX also modulated the imbalance in the MDM2-p53 interaction. The molecular docking results confirmed the interaction of TX with the upstream kinases that regulate apoptosis. This study provides evidence that a copious dose of TX is required to counteract the differential mitoinhibitory effect of 2-AAF in NDEA initiated hepatomas, and TX exhibits an anti-tumoral effect via suppressing oxidative stress, regulating liver function enzymes, inhibiting inflammatory responses and modulating MDM2-p53 interactions, thus inducing apoptosis, and thereby suggesting that TX may provide promising therapeutic effects for the chemoprevention of HCC.

Synthesis and pre-clinical evaluation of a new class of high-affinity 18F-labeled PSMA ligands for detection of prostate cancer by PET imaging

Purpose

Current clinical imaging of PSMA-positive prostate cancer by positron emission tomography (PET) mainly features ⁶⁸Ga-labeled tracers, notably [⁶⁸Ga]Ga-PSMA-HBED-CC. The longer half-life of fluorine-18 offers significant advantages over Ga-68, clinically and logistically. We aimed to develop high-affinity PSMA inhibitors labeled with fluorine-18 as alternative tracers for prostate cancer.

Methods

Six triazolylphenyl ureas and their alkyne precursors were synthesized from the Glu-urea-Lys PSMA binding moiety. PSMA affinity was determined in a competitive binding assay using LNCaP cells. The [¹⁸F]triazoles were isolated following a Cu(I)-catalyzed click reaction between the alkynes and [¹⁸F]fluoroethylazide. The ¹⁸F-labeled compounds were evaluated in nude mice bearing LNCaP tumors and compared to [⁶⁸Ga]Ga-PSMA-HBED-CC and [¹⁸F]DCFPyL. Biodistribution studies of the two tracers with the highest imaged-derived tumor uptake and highest PSMA affinity were undertaken at 1 h, 2 h and 4 h post-injection (p.i.), and co-administration of PMPA was used to determine whether uptake was PSMA-specific.

Results

F-18-labeled triazolylphenyl ureas were prepared with a decay-corrected RCY of 20–40 %, >98 % radiochemical and chemical purity, and specific activity of up to 391 GBq/μmol. PSMA binding (IC₅₀) ranged from 3–36 nM. The position of the triazole influenced tumor uptake (3 > 4 > 2), and direct conjugation of the triazole with the phenylurea moiety was preferred to insertion of a spacer group. Image-derived tumor uptake ranged from 6–14 %ID/g at 2 h p.i., the time of maximum tumor uptake; uptake of [⁶⁸Ga]Ga-PSMA-HBED-CC and [¹⁸F]DCFPyL was 5–6 %ID/g at 1–3 h p.i., the time of maximum tumor uptake. Biodistribution studies of the two most promising compounds gave maximum tumor uptakes of 10.9 ± 1.0 % and 14.3 ± 2.5 %ID/g, respectively, as compared to 6.27 ± 1.44 %ID/g for [⁶⁸Ga]Ga-PSMA-HBED-CC.

Conclusions

Six [¹⁸F]triazolylphenyl ureas were prepared in good radiochemical yield. Compounds showed PSMA-specific uptake in LNCaP tumors as high as 14 % ID/g, more than a 2-fold increase over [⁶⁸Ga]Ga-PSMA-HBED-CC. The facile and high-yielding radiosynthesis of these ¹⁸F-labeled triazoles as well as their promising in vitro and in vivo characteristics make them worthy of clinical development for PET imaging of prostate cancer.

Electronic supplementary material

The online version of this article (doi:10.1007/s00259-016-3556-5) contains supplementary material, which is available to authorized users.