Discovery of novel benzothiophene derivatives as potent and narrow spectrum inhibitors of DYRK1A and DYRK1B.

The discovery of potent and selective inhibitors for understudied kinases can provide relevant pharmacological tools to illuminate their biological functions. DYRK1A and DYRK1B are protein kinases linked to chronic human diseases. Current DYRK1A/DYRK1B inhibitors also antagonize the function of related protein kinases, such as CDC2-like kinases (CLK1, CLK2, CLK4) and DYRK2. Here, we reveal narrow spectrum dual inhibitors of DYRK1A and DYRK1B based on a benzothiophene scaffold. Compound optimization exploited structural differences in the ATP-binding sites of the DYRK1 kinases and resulted in the discovery of 3n, a potent and cell-permeable DYRK1A/DYRK1B inhibitor. This compound has a different scaffold and a narrower off-target profile compared to current DYRK1A/DYRK1B inhibitors. We expect the benzothiophene derivatives described here to aid establishing DYRK1A/DYRK1B cellular functions and their role in human pathologies.

Novel Dihydropteridinone Derivatives As Potent Inhibitors of the Understudied Human Kinases Vaccinia-Related Kinase 1 and Casein Kinase 1δ/ε.

Vaccinia-related kinase 1 (VRK1) and the δ and ε isoforms of casein kinase 1 (CK1) are linked to various disease-relevant pathways. However, the lack of tool compounds for these kinases has significantly hampered our understanding of their cellular functions and therapeutic potential. Here, we describe the structure-based development of potent inhibitors of VRK1, a kinase highly expressed in various tumor types and crucial for cell proliferation and genome integrity. Kinome-wide profiling revealed that our compounds also inhibit CK1δ and CK1ε. We demonstrate that dihydropteridinones 35 and 36 mimic the cellular outcomes of VRK1 depletion. Complementary studies with existing CK1δ and CK1ε inhibitors suggest that these kinases may play overlapping roles in cell proliferation and genome instability. Together, our findings highlight the potential of VRK1 inhibition in treating p53-deficient tumors and possibly enhancing the efficacy of existing cancer therapies that target DNA stability or cell division.

From partial to full agonism: identification of a novel 2,4,5,6-tetrahydropyrrolo[3,4-c]pyrazole as a full agonist of the human GPR119 receptor.

A novel GPR119 agonist based on the 2,4,5,6-tetrahydropyrrolo[3,4-c]pyrazole scaffold was designed through lead optimization starting from pyrazole-based GPR119 agonist 1. The design is centered on the conformational restriction of the core scaffold, while minimizing the change in spatial relationships of two key pharmacophoric elements (piperidine-carbamate and aryl sulfone).

Deconstruction of activity-dependent covalent modification of heme in human neutrophil myeloperoxidase by multistage mass spectrometry (MS(4)).

Myeloperoxidase (MPO) is known to be inactivated and covalently modified by treatment with hydrogen peroxide and agents similar to 3-(2-ethoxypropyl)-2-thioxo-2,3-dihydro-1H-purin-6(9H)-one (1), a 254.08 Da derivative of 2-thioxanthine. Peptide mapping by liquid chromatography and mass spectrometry detected modification by 1 in a labile peptide-heme-peptide fragment of the enzyme, accompanied by a mass increase of 252.08 Da. The loss of two hydrogen atoms was consistent with mechanism-based oxidative coupling. Multistage mass spectrometry (MS(4)) of the modified fragment in an ion trap/Orbitrap spectrometer demonstrated that 1 was coupled directly to heme. Use of a 10 amu window delivered the full isotopic envelope of each precursor ion to collision-induced dissociation, preserving definitive isotopic profiles for iron-containing fragments through successive steps of multistage mass spectrometry. Iron isotope signatures and accurate mass measurements supported the structural assignments. Crystallographic analysis confirmed linkage between the methyl substituent of the heme pyrrole D ring and the sulfur atom of 1. The final orientation of 1 perpendicular to the plane of the heme ring suggested a mechanism consisting of two consecutive one-electron oxidations of 1 by MPO. Multistage mass spectrometry using stage-specific collision energies permits stepwise deconstruction of modifications of heme enzymes containing covalent links between the heme group and the polypeptide chain.

Identification of novel series of pyrazole and indole-urea based DFG-out PYK2 inhibitors.

Previous drug discovery efforts identified classical PYK2 kinase inhibitors such as 2 and 3 that possess selectivity for PYK2 over its intra-family isoform FAK. Efforts to identify more kinome-selective chemical matter that stabilize a DFG-out conformation of the enzyme are described herein. Two sub-series of PYK2 inhibitors, an indole carboxamide-urea and a pyrazole-urea have been identified and found to have different binding interactions with the hinge region of PYK2. These leads proved to be more selective than the original classical inhibitors.

Use of 3D properties to characterize beyond rule-of-5 property space for passive permeation.

The application of conformationally dependent measures of size and polarity to characterize beyond rule-of-5 (Ro5) space for passive permeation was investigated. Specifically, radius of gyration, an alternative to molecular weight, and three-dimensional polar surface area and the generalized Born/surface area dehydration free energy, alternatives to hydrogen-bond donor and acceptor counts, were computed on models of the permeating conformations of over 35 000 molecules. The resulting guidelines for size and polarity, described by the 3D properties, should aid the design of Ro5 violators with passive permeability.

Exploring aromatic chemical space with NEAT: novel and electronically equivalent aromatic template.

In this paper, we describe a lead transformation tool, NEAT (Novel and Electronically equivalent Aromatic Template), which can help identify novel aromatic rings that are estimated to have similar electrostatic potentials, dipoles, and hydrogen bonding capabilities to a query template; hence, they may offer similar bioactivity profiles. In this work, we built a comprehensive heteroaryl database, and precalculated high-level quantum mechanical (QM) properties, including electrostatic potential charges, hydrogen bonding ability, dipole moments, chemical reactivity, and othe properties. NEAT bioisosteric similarities are based on the electrostatic potential surface calculated by Brood, using the precalculated QM ESP charges and other QM properties. Compared with existing commercial lead transformation software, (1) NEAT is the only one that covers the comprehensive heteroaryl chemical space, and (2) NEAT offers a better characterization of novel aryl cores by using high-evel QM properties that are relevant to molecular interactions. NEAT provides unique value to medicinal chemists quickly exploring the largely uncharted aromatic chemical space, and one successful example of its application is discussed herein.

Design and evaluation of a 2-(2,3,6-trifluorophenyl)acetamide derivative as an agonist of the GPR119 receptor.

The design and synthesis of a GPR119 agonist bearing a 2-(2,3,6-trifluorophenyl)acetamide group is described. The design capitalized on the conformational restriction found in N-ß-fluoroethylamide derivatives to help maintain good levels of potency while driving down both lipophilicity and oxidative metabolism in human liver microsomes. The chemical stability and bioactivation potential are discussed.

Understanding the impact of the P-loop conformation on kinase selectivity.

This work addresses the link between selectivity and an unusual, folded conformation for the P-loop observed initially for MAP4K4 and subsequently for other kinases. Statistical and computational analyses of our crystal structure database demonstrate that inhibitors that induce the P-loop folded conformation tend to be more selective, especially if they take advantage of this specific conformation by interacting more favorably with a conserved Tyr or Phe residue from the P-loop.

Addressing limitations with the MM-GB/SA scoring procedure using the WaterMap method and free energy perturbation calculations.

The MM-GB/SA scoring technique has become an important computational approach in drug design. We, and others, have demonstrated that for congeneric molecules the correlation with experimental data obtained with the physics-based scoring is usually superior to scoring functions from typical docking algorithms. Despite showing good accuracy when applied within a series, much work is necessary to improve the MM-GB/SA method in order to gain greater efficiency in drug design. Here, we investigate the poor estimation of protein desolvation provided by the GB/SA solvation model and the large dynamic range observed in the MM-GB/SA scoring compared to that of the experimental data. In the former, replacing the GB/SA protein desolvation in the MM-GB/SA method by the free energy associated with displacing binding site waters upon ligand binding estimated by WaterMap provides the best results when ranking congeneric series of factor Xa and cyclin-dependent kinase 2 (CDK2) inhibitors. However, the improvement is modest over results obtained with the MM-GB/SA and WaterMap methods individually, apparently due to the high correlation between the free energy liberation of the displaced solvent and the protein-ligand van der Waals interactions, which in turn may be interpretable as estimates of the hydrophobic effect and hydrophobic-like interactions, respectively. As for the large dynamic range, comparisons between MM-GB/SA and FEP calculations indicate that for the factor Xa test set this problem has its origin in the lack of shielding effects of protein--ligand electrostatic interactions; that overly favors ligands that engage in hydrogen bonds with the protein.

Thermodynamic analysis of mRNA cap binding by the human initiation factor eIF4E via free energy perturbations.

Eukaryotic mRNAs are appended at the 5' end, with the 7-methylguanosine cap linked by a 5'-5'-triphosphate bridge to the first transcribed nucleoside (m7GpppX). Initiation of cap-dependent translation of mRNA requires direct interaction between the cap structure and the eukaryotic translation initiation factor eIF4E. Biophysical studies of the association between eIF4E and various cap analogs have demonstrated that m(7)GTP binds to the protein ca. -5.0 kcal/mol more favorably than unmethylated GTP. In this work, a thermodynamic analysis of the binding process between eIF4E and several cap analogs has been conducted using Monte Carlo (MC) simulations in conjunction with free energy perturbation (FEP) calculations. To address the role of the 7-methyl group in the eIF4E/m7GpppX cap interaction, binding free energies have been computed for m(7)GTP, GTP, protonated GTP at N(7), the 7-methyldeazaguanosine 5'-triphosphate (m(7)DTP), and 7-deazaguanosine 5'-triphosphate (DTP) cap analogs. The MC/FEP simulations for the GTP-->m(7)DTP transformation demonstrate that half of the binding free energy gain of m(7)GTP with respect to GTP can be attributed to favorable van der Waals interactions with Trp166 and reduced desolvation penalty due to the N(7) methyl group. The methyl group both eliminates the desolvation penalty of the N(7) atom upon binding and creates a larger cavity within the solvent that further facilitates the desolvation step. Analysis of the pair m(7)GTP-m(7)DTP suggests that the remaining gain in affinity is related to the positive charge created on the guanine moiety due to the N(7) methylation. The charge provides favorable cation-pi interactions with Trp56 and Trp102 and decreases the negative molecular charge, which helps the transfer from the solvent, a more polar environment, to the protein.

Development of Pyridine-based Inhibitors for the Human Vaccinia-related Kinases 1 and 2.

Vaccinia-related kinases 1 and 2 (VRK1 and VRK2) are human Ser/Thr protein kinases associated with increased cell division and neurological disorders. Nevertheless, the cellular functions of these proteins are not fully understood. Despite their therapeutic potential, there are no potent and specific inhibitors available for VRK1 or VRK2. We report here the discovery and elaboration of an aminopyridine scaffold as a basis for VRK1 and VRK2 inhibitors. The most potent compound for VRK1 (26) displayed an IC50 value of 150 nM and was fairly selective in a panel of 48 human kinases (selectivity score S(50%) of 0.04). Differences in compound binding mode and substituent preferences between the two VRKs were identified by the structure-activity relationship combined with the crystallographic analysis of key compounds. We expect our results to serve as a starting point for the design of more specific and potent inhibitors against each of the two VRKs.

2-Aminopyridine-Based Mitogen-Activated Protein Kinase Kinase Kinase Kinase 4 (MAP4K4) Inhibitors: Assessment of Mechanism-Based Safety.

Studies have linked the serine-threonine kinase MAP4K4 to the regulation of a number of biological processes and/or diseases, including diabetes, cancer, inflammation, and angiogenesis. With a majority of the members of our lead series (e.g., 1) suffering from time-dependent inhibition (TDI) of CYP3A4, we sought design avenues that would eliminate this risk. One such approach arose from the observation that carboxylic acid-based intermediates employed in our discovery efforts retained high MAP4K4 inhibitory potency and were devoid of the TDI risk. The medicinal chemistry effort that led to the discovery of this central nervous system-impaired inhibitor together with its preclinical safety profile is described.

From docking false-positive to active anti-HIV agent.

Virtual screening of the Maybridge library of ca. 70 000 compounds was performed using a similarity filter, docking, and molecular mechanics-generalized Born/surface area postprocessing to seek potential non-nucleoside inhibitors of human immunodeficiency virus-1 (HIV-1) reverse transcriptase (NNRTIs). Although known NNRTIs were retrieved well, purchase and assaying of representative, top-scoring compounds from the library failed to yield any active anti-HIV agents. However, the highest-ranked library compound, oxadiazole 1, was pursued as a potential "near-miss" with the BOMB program to seek constructive modifications. Subsequent synthesis and assaying of several polychloro-analogs did yield anti-HIV agents with EC50 values as low as 310 nM. The study demonstrates that it is possible to learn from a formally unsuccessful virtual-screening exercise and, with the aid of computational analyses, to efficiently evolve a false positive into a true active.

Structural characterization of human Vaccinia-Related Kinases (VRK) bound to small-molecule inhibitors identifies different P-loop conformations.

The human genome encodes two active Vaccinia-related protein kinases (VRK), VRK1 and VRK2. These proteins have been implicated in a number of cellular processes and linked to a variety of tumors. However, understanding the cellular role of VRKs and establishing their potential use as targets for therapeutic intervention has been limited by the lack of tool compounds that can specifically modulate the activity of these kinases in cells. Here we identified BI-D1870, a dihydropteridine inhibitor of RSK kinases, as a promising starting point for the development of chemical probes targeting the active VRKs. We solved co-crystal structures of both VRK1 and VRK2 bound to BI-D1870 and of VRK1 bound to two broad-spectrum inhibitors. These structures revealed that both VRKs can adopt a P-loop folded conformation, which is stabilized by different mechanisms on each protein. Based on these structures, we suggest modifications to the dihydropteridine scaffold that can be explored to produce potent and specific inhibitors towards VRK1 and VRK2.

Discovery of a potent, selective, and efficacious class of reversible alpha-ketoheterocycle inhibitors of fatty acid amide hydrolase effective as analgesics.

Fatty acid amide hydrolase (FAAH) degrades neuromodulating fatty acid amides including anandamide (endogenous cannabinoid agonist) and oleamide (sleep-inducing lipid) at their sites of action and is intimately involved in their regulation. Herein we report the discovery of a potent, selective, and efficacious class of reversible FAAH inhibitors that produce analgesia in animal models validating a new therapeutic target for pain intervention. Key to the useful inhibitor discovery was the routine implementation of a proteomics-wide selectivity screen against the serine hydrolase superfamily ensuring selectivity for FAAH coupled with systematic in vivo examinations of candidate inhibitors.

Discovery of an in Vivo Tool to Establish Proof-of-Concept for MAP4K4-Based Antidiabetic Treatment.

Recent studies in adipose tissue, pancreas, muscle, and macrophages suggest that MAP4K4, a serine/threonine protein kinase may be a viable target for antidiabetic drugs. As part of the evaluation of MAP4K4 as a novel antidiabetic target, a tool compound, 16 (PF-6260933) and a lead 17 possessing excellent kinome selectivity and suitable properties were delivered to establish proof of concept in vivo. The medicinal chemistry effort that led to the discovery of these lead compounds is described herein together with in vivo pharmacokinetic properties and activity in a model of insulin resistance.

Enantioselective hydroarylation of bridged [3.2.1] heterocycles: an efficient entry into the homoepibatidine skeleton.

Achiral [3.2.1] bridged heterocycles containing a bridging amide can undergo enantioselective hydroarylation reactions under rhodium(I) catalysis. These reactions proceed in high yield and enantioselectivity in most cases, under mild reaction conditions and using commercially available Josiphos ligands. The phosphine ligand structure and the protecting group on the nitrogen both have significant effects on the selectivity and yield of the reactions.

Design and synthesis of diazatricyclodecane agonists of the G-protein-coupled receptor 119.

A series of GPR119 agonists based on a 2,6-diazatricyclo[3.3.1.1∼3,7∼]decane ring system is described. Also provided is a detailed account of the development of a multigram scale synthesis of the diazatricyclic ring system, which was achieved using a Hofmann-Löffler-Freytag reaction as the key step. The basis for the use of this complex framework lies in an attempt to constrain one end of the molecule in the "agonist conformation" as was previously described for 3-oxa-7-aza-bicyclo[3.3.1]nonanes. Optimization of carbamate analogues of the diazatricylic compounds led to the identification of 32i as a potent agonist of the GPR119 receptor with low unbound human liver microsomal clearance. The use of an agonist response weighted ligand lipophilic efficiency (LLE) termed AgLLE is discussed along with the issues of applying efficiency measures to agonist programs. Ultimately, solubility limited absorption and poor exposure reduced further interest in these molecules.

MM-GB/SA rescoring of docking poses.

The critical issues in docking include the prediction of the correct binding pose and the accurate estimation of the corresponding binding affinity. Different docking methodologies have all been successful in reproducing the crystallographic binding modes, but struggle when predicting the corresponding binding affinities. The rescoring of docking poses using the MM-GB/SA technique has emerged as an important computational approach in structure-based lead optimization as it provides for congeneric molecules, clearly superior correlations with experimental data to those obtained with typical docking scoring functions. Although the technique has been collectively referred as MM-GB/SA, there are in fact many flavors in the literature. Here we describe the details of our MM-GB/SA scoring protocol, highlighting not only its strengths but also the limitations.