Enhancing Hit Discovery in Virtual Screening through Absolute Protein-Ligand Binding Free-Energy Calculations.

In the hit identification stage of drug discovery, a diverse chemical space needs to be explored to identify initial hits. Contrary to empirical scoring functions, absolute protein-ligand binding free-energy perturbation (ABFEP) provides a theoretically more rigorous and accurate description of protein-ligand binding thermodynamics and could, in principle, greatly improve the hit rates in virtual screening. In this work, we describe an implementation of an accurate and reliable ABFEP method in FEP+. We validated the ABFEP method on eight congeneric compound series binding to eight protein receptors including both neutral and charged ligands. For ligands with net charges, the alchemical ion approach is adopted to avoid artifacts in electrostatic potential energy calculations. The calculated binding free energies correlate with experimental results with a weighted average of R2 = 0.55 for the entire dataset. We also observe an overall root-mean-square error (RMSE) of 1.1 kcal/mol after shifting the zero-point of the simulation data to match the average experimental values. Through ABFEP calculations using apo versus holo protein structures, we demonstrated that the protein conformational and protonation state changes between the apo and holo proteins are the main physical factors contributing to the protein reorganization free energy manifested by the overestimation of raw ABFEP calculated binding free energies using the holo structures of the proteins. Furthermore, we performed ABFEP calculations in three virtual screening applications for hit enrichment. ABFEP greatly improves the hit rates as compared to docking scores or other methods like metadynamics. The good performance of ABFEP in rank ordering compounds demonstrated in this work confirms it as a useful tool to improve the hit rates in virtual screening, thus facilitating hit discovery.

Antimalarial proteasome inhibitor reveals collateral sensitivity from intersubunit interactions and fitness cost of resistance.

We describe noncovalent, reversible asparagine ethylenediamine (AsnEDA) inhibitors of the Plasmodium falciparum proteasome (Pf20S) ß5 subunit that spare all active subunits of human constitutive and immuno-proteasomes. The compounds are active against erythrocytic, sexual, and liver-stage parasites, against parasites resistant to current antimalarials, and against P. falciparum strains from patients in Africa. The ß5 inhibitors synergize with a ß2 inhibitor in vitro and in mice and with artemisinin. P. falciparum selected for resistance to an AsnEDA ß5 inhibitor surprisingly harbored a point mutation in the noncatalytic ß6 subunit. The ß6 mutant was resistant to the species-selective Pf20S ß5 inhibitor but remained sensitive to the species-nonselective ß5 inhibitors bortezomib and carfilzomib. Moreover, resistance to the Pf20S ß5 inhibitor was accompanied by increased sensitivity to a Pf20S ß2 inhibitor. Finally, the ß5 inhibitor-resistant mutant had a fitness cost that was exacerbated by irradiation. Thus, used in combination, multistage-active inhibitors of the Pf20S ß5 and ß2 subunits afford synergistic antimalarial activity with a potential to delay the emergence of resistance to artemisinins and each other.

Development of a Highly Selective Plasmodium falciparum Proteasome Inhibitor with Anti-malaria Activity in Humanized Mice.

Plasmodium falciparum proteasome (Pf20S) inhibitors are active against Plasmodium at multiple stages-erythrocytic, gametocyte, liver, and gamete activation stages-indicating that selective Pf20S inhibitors possess the potential to be therapeutic, prophylactic, and transmission-blocking antimalarials. Starting from a reported compound, we developed a noncovalent, macrocyclic peptide inhibitor of the malarial proteasome with high species selectivity and improved pharmacokinetic properties. The compound demonstrates specific, time-dependent inhibition of the ß5 subunit of the Pf20S, kills artemisinin-sensitive and artemisinin-resistant P. falciparum isolates in vitro and reduces parasitemia in humanized, P. falciparum-infected mice.

The E2.65A mutation disrupts dynamic binding poses of SB269652 at the dopamine D2 and D3 receptors.

The dopamine D2 and D3 receptors (D2R and D3R) are important targets for antipsychotics and for the treatment of drug abuse. SB269652, a bitopic ligand that simultaneously binds both the orthosteric binding site (OBS) and a secondary binding pocket (SBP) in both D2R and D3R, was found to be a negative allosteric modulator. Previous studies identified Glu2.65 in the SBP to be a key determinant of both the affinity of SB269652 and the magnitude of its cooperativity with orthosteric ligands, as the E2.65A mutation decreased both of these parameters. However, the proposed hydrogen bond (H-bond) between Glu2.65 and the indole moiety of SB269652 is not a strong interaction, and a structure activity relationship study of SB269652 indicates that this H-bond may not be the only element that determines its allosteric properties. To understand the structural basis of the observed phenotype of E2.65A, we carried out molecular dynamics simulations with a cumulative length of ~77 µs of D2R and D3R wild-type and their E2.65A mutants bound to SB269652. In combination with Markov state model analysis and by characterizing the equilibria of ligand binding modes in different conditions, we found that in both D2R and D3R, whereas the tetrahydroisoquinoline moiety of SB269652 is stably bound in the OBS, the indole-2-carboxamide moiety is dynamic and only intermittently forms H-bonds with Glu2.65. Our results also indicate that the E2.65A mutation significantly affects the overall shape and size of the SBP, as well as the conformation of the N terminus. Thus, our findings suggest that the key role of Glu2.65 in mediating the allosteric properties of SB269652 extends beyond a direct interaction with SB269652, and provide structural insights for rational design of SB269652 derivatives that may retain its allosteric properties.

Inhibition of 3-phosphoglycerate dehydrogenase (PHGDH) by indole amides abrogates de novo serine synthesis in cancer cells.

Cancer cells reprogram their metabolism to support growth and to mitigate cellular stressors. The serine synthesis pathway has been identified as a metabolic pathway frequently altered in cancers and there has been considerable interest in developing pharmacological agents to target this pathway. Here, we report a series of indole amides that inhibit human 3-phosphoglycerate dehydrogenase (PHGDH), the enzyme that catalyzes the first committed step of the serine synthesis pathway. Using X-ray crystallography, we show that the indole amides bind the NAD+ pocket of PHGDH. Through structure-based optimization we were able to develop compounds with low nanomolar affinities for PHGDH in an enzymatic IC50 assay. In cellular assays, the most potent compounds inhibited de novo serine synthesis with low micromolar to sub-micromolar activities and these compounds successfully abrogated the proliferation of cancer cells in serine free media. The indole amide series reported here represent an important improvement over previously published PHGDH inhibitors as they are markedly more potent and their mechanism of action is better defined.

What can crystal structures of aminergic receptors tell us about designing subtype-selective ligands?

G protein-coupled receptors (GPCRs) are integral membrane proteins that represent an important class of drug targets. In particular, aminergic GPCRs interact with a significant portion of drugs currently on the market. However, most drugs that target these receptors are associated with undesirable side effects, which are due in part to promiscuous interactions with close homologs of the intended target receptors. Here, based on a systematic analysis of all 37 of the currently available high-resolution crystal structures of aminergic GPCRs, we review structural elements that contribute to and can be exploited for designing subtype-selective compounds. We describe the roles of secondary binding pockets (SBPs), as well as differences in ligand entry pathways to the orthosteric binding site, in determining selectivity. In addition, using the available crystal structures, we have identified conformational changes in the SBPs that are associated with receptor activation and explore the implications of these changes for the rational development of selective ligands with tailored efficacy.

Computational approaches to detect allosteric pathways in transmembrane molecular machines.

Many of the functions of transmembrane proteins involved in signal processing and transduction across the cell membrane are determined by allosteric couplings that propagate the functional effects well beyond the original site of activation. Data gathered from breakthroughs in biochemistry, crystallography, and single molecule fluorescence have established a rich basis of information for the study of molecular mechanisms in the allosteric couplings of such transmembrane proteins. The mechanistic details of these couplings, many of which have therapeutic implications, however, have only become accessible in synergy with molecular modeling and simulations. Here, we review some recent computational approaches that analyze allosteric coupling networks (ACNs) in transmembrane proteins, and in particular the recently developed Protein Interaction Analyzer (PIA) designed to study ACNs in the structural ensembles sampled by molecular dynamics simulations. The power of these computational approaches in interrogating the functional mechanisms of transmembrane proteins is illustrated with selected examples of recent experimental and computational studies pursued synergistically in the investigation of secondary active transporters and GPCRs. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

A new mechanism of allostery in a G protein-coupled receptor dimer.

SB269652 is to our knowledge the first drug-like allosteric modulator of the dopamine D2 receptor (D2R), but it contains structural features associated with orthosteric D2R antagonists. Using a functional complementation system to control the identity of individual protomers within a dimeric D2R complex, we converted the pharmacology of the interaction between SB269652 and dopamine from allosteric to competitive by impairing ligand binding to one of the protomers, indicating that the allostery requires D2R dimers. Additional experiments identified a 'bitopic' pose for SB269652 extending from the orthosteric site into a secondary pocket at the extracellular end of the transmembrane (TM) domain, involving TM2 and TM7. Engagement of this secondary pocket was a requirement for the allosteric pharmacology of SB269652. This suggests a new mechanism whereby a bitopic ligand binds in an extended pose on one G protein-coupled receptor protomer to allosterically modulate the binding of a ligand to the orthosteric site of a second protomer.

Next-generation neuropeptide Y receptor small-molecule agonists inhibit mosquito-biting behavior.

BACKGROUND: Female Aedes aegypti mosquitoes can spread disease-causing pathogens when they bite humans to obtain blood nutrients required for egg production. Following a complete blood meal, host-seeking is suppressed until eggs are laid. Neuropeptide Y-like receptor 7 (NPYLR7) plays a role in endogenous host-seeking suppression and previous work identified small-molecule NPYLR7 agonists that inhibit host-seeking and blood-feeding when fed to mosquitoes at high micromolar doses. METHODS: Using structure-activity relationship analysis and structure-guided design we synthesized 128 compounds with similarity to known NPYLR7 agonists. RESULTS: Although in vitro potency (EC50) was not strictly predictive of in vivo effect, we identified three compounds that reduced blood-feeding from a live host when fed to mosquitoes at a dose of 1 µM-a 100-fold improvement over the original reference compound. CONCLUSIONS: Exogenous activation of NPYLR7 represents an innovative vector control strategy to block mosquito biting behavior and prevent mosquito-human host interactions that lead to pathogen transmission.

Next Generation Neuropeptide Y Receptor Small Molecule Agonists Inhibit Mosquito Biting Behavior.

Female Aedes aegypti mosquitoes can spread disease-causing pathogens when they bite humans to obtain blood nutrients required for egg production. Following a complete blood meal, host-seeking is suppressed until eggs are laid. Neuropeptide Y-like Receptor 7 (NPYLR7) plays a role in endogenous host-seeking suppression and previous work identified small molecule NPYLR7 agonists that suppress host-seeking and blood feeding when fed to mosquitoes at high micromolar doses. Using structure activity relationship analysis and structure-guided design we synthesized 128 compounds with similarity to known NPYLR7 agonists. Although in vitro potency (EC50) was not strictly predictive of in vivo effect, we identified 3 compounds that suppressed blood feeding from a live host when fed to mosquitoes at a 1 µM dose, a 100-fold improvement over the original reference compound. Exogenous activation of NPYLR7 represents an innovative vector control strategy to block mosquito biting behavior and prevent mosquito/human host interactions that lead to pathogen transmission.

Lead Optimization of Small Molecule ENL YEATS Inhibitors to Enable In Vivo Studies: Discovery of TDI-11055.

Eleven-nineteen leukemia (ENL) is an epigenetic reader protein that drives oncogenic transcriptional programs in acute myeloid leukemia (AML). AML is one of the deadliest hematopoietic malignancies, with an overall 5-year survival rate of 27%. The epigenetic reader activity of ENL is mediated by its YEATS domain that binds to acetyl and crotonyl marks on histone tails and colocalizes with promoters of actively transcribed genes that are essential for leukemia. Prior to the discovery of TDI-11055, existing inhibitors of ENL YEATS showed in vitro potency, but had not shown efficacy in in vivo animal models. During the course of the medicinal chemistry campaign described here, we identified ENL YEATS inhibitor TDI-11055 that has an improved pharmacokinetic profile and is appropriate for in vivo evaluation of the ENL YEATS inhibition mechanism in AML.

A single glycine in extracellular loop 1 is the critical determinant for pharmacological specificity of dopamine D2 and D3 receptors.

Subtype-selective agents for the dopamine D3 receptor (D3R) have been considered as potential medications for drug addiction and other neuropsychiatric disorders. Medicinal chemistry efforts have led to the discovery of 4-phenylpiperazine derivatives that are >100-fold selective for the dopamine D3 receptor over dopamine D2 receptor (D2R), despite high sequence identity (78% in the transmembrane domain). Based on the recent crystal structure of D3R, we demonstrated that the 4-phenylpiperazine moiety in this class of D3R-selective compounds binds to the conserved orthosteric binding site, whereas the extended aryl amide moiety is oriented toward a divergent secondary binding pocket (SBP). In an effort to further characterize molecular determinants of the selectivity of these compounds, we modeled their binding modes in D3R and D2R by comparative ligand docking and molecular dynamics simulations. We found that the aryl amide moiety in the SBP differentially induces conformational changes in transmembrane segment 2 and extracellular loop 1 (EL1), which amplify the divergence of the SBP in D3R and D2R. Receptor chimera and site-directed mutagenesis studies were used to validate these binding modes and to identify a divergent glycine in EL1 as critical to D3R over D2R subtype selectivity. A better understanding of drug-dependent receptor conformations such as these is key to the rational design of compounds targeting a specific receptor among closely related homologs, and may also lead to discovery of novel chemotypes that exploit subtle differences in protein conformations.

Shape-Based Virtual Screening of a Billion-Compound Library Identifies Mycobacterial Lipoamide Dehydrogenase Inhibitors.

Lpd (lipoamide dehydrogenase) in Mycobacterium tuberculosis (Mtb) is required for virulence and is a genetically validated tuberculosis (TB) target. Numerous screens have been performed over the last decade, yet only two inhibitor series have been identified. Recent advances in large-scale virtual screening methods combined with make-on-demand compound libraries have shown the potential for finding novel hits. In this study, the Enamine REAL library consisting of ∼1.12 billion compounds was efficiently screened using the GPU Shape screen method against Mtb Lpd to find additional chemical matter that would expand on the known sulfonamide inhibitor series. We identified six new inhibitors with IC50 in the range of 5-100 µM. While these compounds remained chemically close to the already known sulfonamide series inhibitors, some diversity was found in the cores of the hits. The two most potent hits were further validated by one-step potency optimization to submicromolar levels. The co-crystal structure of optimized analogue TDI-13537 provided new insights into the potency determinants of the series.

Structure-Activity Relationship Studies of Antimalarial Plasmodium Proteasome Inhibitors─Part II.

With increasing reports of resistance to artemisinins and artemisinin-combination therapies, targeting the Plasmodium proteasome is a promising strategy for antimalarial development. We recently reported a highly selective Plasmodium falciparum proteasome inhibitor with anti-malarial activity in the humanized mouse model. To balance the permeability of the series of macrocycles with other drug-like properties, we conducted further structure-activity relationship studies on a biphenyl ether-tethered macrocyclic scaffold. Extensive SAR studies around the P1, P3, and P5 groups and peptide backbone identified compound TDI-8414. TDI-8414 showed nanomolar antiparasitic activity, no toxicity to HepG2 cells, high selectivity against the Plasmodium proteasome over the human constitutive proteasome and immunoproteasome, improved solubility and PAMPA permeability, and enhanced metabolic stability in microsomes and plasma of both humans and mice.

Characterization of Pyrrolidinyl-hexahydro-pyranopiperazines as a Novel Kappa Opioid Receptor Agonist Scaffold.

The kappa agonist structure-activity relationship around the novel, pyrrolidinyl substituted pyranopiperazine scaffold was developed. More specifically, the dichloroPhenylAcetamide-Pyrrolidinyl-PyranoPiperazine (PAPPP) core A was the focus of our work. The modulation of kappa receptor potency/G-protein activation and arrestin recruitment with respect to changes of the piperazine R group in A was demonstrated. Reduced ß2-arrestin recruitment and differential G-protein bias were observed for select analogues. To better understand the subtlety in receptor signaling, analogues were profiled as the resolved enantiomers. To determine in vivo target engagement, a subset of compounds was tested in mice for stimulation of serum prolactin, a neuroendocrine biomarker of KOR-agonist effects. Additional in vivo characterization included measurement of potential unwanted effects of kappa receptor activation such as sedation. These studies demonstrate a novel kappa receptor agonist scaffold with potential for G-protein signaling bias to probe in vivo pharmacology.

Design, Synthesis, and Optimization of Macrocyclic Peptides as Species-Selective Antimalaria Proteasome Inhibitors.

With over 200 million cases and close to half a million deaths each year, malaria is a threat to global health, particularly in developing countries. Plasmodium falciparum, the parasite that causes the most severe form of the disease, has developed resistance to all antimalarial drugs. Resistance to the first-line antimalarial artemisinin and to artemisinin combination therapies is widespread in Southeast Asia and is emerging in sub-Saharan Africa. The P. falciparum proteasome is an attractive antimalarial target because its inhibition kills the parasite at multiple stages of its life cycle and restores artemisinin sensitivity in parasites that have become resistant through mutation in Kelch K13. Here, we detail our efforts to develop noncovalent, macrocyclic peptide malaria proteasome inhibitors, guided by structural analysis and pharmacokinetic properties, leading to a potent, species-selective, metabolically stable inhibitor.

A Chemical Strategy toward Novel Brain-Penetrant EZH2 Inhibitors.

Aberrant gene-silencing through dysregulation of polycomb protein activity has emerged as an important oncogenic mechanism in cancer, implicating polycomb proteins as important therapeutic targets. Recently, an inhibitor targeting EZH2, the methyltransferase component of PRC2, received U.S. Food and Drug Administration approval following promising clinical responses in cancer patients. However, the current array of EZH2 inhibitors have poor brain penetrance, limiting their use in patients with central nervous system malignancies, a number of which have been shown to be sensitive to EZH2 inhibition. To address this need, we have identified a chemical strategy, based on computational modeling of pyridone-containing EZH2 inhibitor scaffolds, to minimize P-glycoprotein activity, and here we report the first brain-penetrant EZH2 inhibitor, TDI-6118 (compound 5). Additionally, in the course of our attempts to optimize this compound, we discovered TDI-11904 (compound 21), a novel, highly potent, and peripherally active EZH2 inhibitor based on a 7 member ring structure.

Small-Molecule Inhibition of the Acyl-Lysine Reader ENL as a Strategy against Acute Myeloid Leukemia.

The chromatin reader eleven-nineteen leukemia (ENL) has been identified as a critical dependency in acute myeloid leukemia (AML), but its therapeutic potential remains unclear. We describe a potent and orally bioavailable small-molecule inhibitor of ENL, TDI-11055, which displaces ENL from chromatin by blocking its YEATS domain interaction with acylated histones. Cell lines and primary patient samples carrying MLL rearrangements or NPM1 mutations are responsive to TDI-11055. A CRISPR-Cas9-mediated mutagenesis screen uncovers an ENL mutation that confers resistance to TDI-11055, validating the compound's on-target activity. TDI-11055 treatment rapidly decreases chromatin occupancy of ENL-associated complexes and impairs transcription elongation, leading to suppression of key oncogenic gene expression programs and induction of differentiation. In vivo treatment with TDI-11055 blocks disease progression in cell line- and patient-derived xenograft models of MLL-rearranged and NPM1-mutated AML. Our results establish ENL displacement from chromatin as a promising epigenetic therapy for molecularly defined AML subsets and support the clinical translation of this approach. SIGNIFICANCE: AML is a poor-prognosis disease for which new therapeutic approaches are desperately needed. We developed an orally bioavailable inhibitor of ENL, demonstrated its potent efficacy in MLL-rearranged and NPM1-mutated AML, and determined its mechanisms of action. These biological and chemical insights will facilitate both basic research and clinical translation. This article is highlighted in the In This Issue feature, p. 2483.

Deglycase-activity oriented screening to identify DJ-1 inhibitors.

The oncoprotein and Parkinson's disease-associated enzyme DJ-1/PARK7 has emerged as a promiscuous deglycase that can remove methylglyoxal-induced glycation adducts from both proteins and nucleotides. However, dissecting its structural and enzymatic functions remains a challenge due to the lack of potent, specific, and pharmacokinetically stable inhibitors targeting its catalytic site (including Cys106). To evaluate potential drug-like leads against DJ-1, we leveraged its deglycase activity in an enzyme-coupled, fluorescence lactate-detection assay based on the recent understanding of its deglycation mechanism. In addition, we developed assays to directly evaluate DJ-1's esterase activity using both colorimetric and fluorescent substrates. The resulting optimized assay was used to evaluate a library of potential reversible and irreversible DJ-1 inhibitors. The deglycase activity-oriented screening strategy described herein establishes a new platform for the discovery of potential anti-cancer drugs.

Whole Cell Active Inhibitors of Mycobacterial Lipoamide Dehydrogenase Afford Selectivity over the Human Enzyme through Tight Binding Interactions.

Tuberculosis remains a leading cause of death from a single bacterial infection worldwide. Efforts to develop new treatment options call for expansion into an unexplored target space to expand the drug pipeline and bypass resistance to current antibiotics. Lipoamide dehydrogenase is a metabolic and antioxidant enzyme critical for mycobacterial growth and survival in mice. Sulfonamide analogs were previously identified as potent and selective inhibitors of mycobacterial lipoamide dehydrogenase in vitro but lacked activity against whole mycobacteria. Here we present the development of analogs with improved permeability, potency, and selectivity, which inhibit the growth of Mycobacterium tuberculosis in axenic culture on carbohydrates and within mouse primary macrophages. They increase intrabacterial pyruvate levels, supporting their on-target activity within mycobacteria. Distinct modalities of binding between the mycobacterial and human enzymes contribute to improved potency and hence selectivity through induced-fit tight binding interactions within the mycobacterial but not human enzyme, as indicated by kinetic analysis and crystallography.