Structure-based design, synthesis, and biological evaluation of imidazo[1,2-b]pyridazine-based p38 MAP kinase inhibitors.

We identified novel potent inhibitors of p38 MAP kinase using structure-based design strategy. X-ray crystallography showed that when p38 MAP kinase is complexed with TAK-715 (1) in a co-crystal structure, Phe169 adopts two conformations, where one interacts with 1 and the other shows no interaction with 1. Our structure-based design strategy shows that these two conformations converge into one via enhanced protein-ligand hydrophobic interactions. According to the strategy, we focused on scaffold transformation to identify imidazo[1,2-b]pyridazine derivatives as potent inhibitors of p38 MAP kinase. Among the herein described and evaluated compounds, N-oxide 16 exhibited potent inhibition of p38 MAP kinase and LPS-induced TNF-α production in human monocytic THP-1 cells, and significant in vivo efficacy in rat collagen-induced arthritis models. In this article, we report the discovery of potent, selective and orally bioavailable imidazo[1,2-b]pyridazine-based p38 MAP kinase inhibitors with pyridine N-oxide group.

Multiscale enhanced sampling of glucokinase: Regulation of the enzymatic reaction via a large scale domain motion.

Enhanced sampling yields a comprehensive structural ensemble or a free energy landscape, which is beyond the capability of a conventional molecular dynamics simulation. Our recently developed multiscale enhanced sampling (MSES) method employs a coarse-grained model coupled with the target physical system for the efficient acceleration of the dynamics. MSES has demonstrated applicability to large protein systems in solution, such as intrinsically disordered proteins and protein-protein and protein-ligand interactions. Here, we applied the MSES simulation to an important drug discovery target, glucokinase (GCK), to elucidate the structural basis of the positive cooperativity of the enzymatic reaction at an atomistic resolution. MSES enabled us to compare two sets of the free energy landscapes of GCK, for the glucose-bound and glucose-unbound forms, and thus demonstrated the drastic change of the free energy surface depending on the glucose concentration. In the glucose-bound form, we found two distinct basins separated by a high energy barrier originating from the domain motion and the folding/unfolding of the α13 helix. By contrast, in the glucose-unbound form, a single flat basin extended to the open and super-open states. These features illustrated the two distinct phases achieving the cooperativity, the fast reaction cycle staying in the closed state at a high glucose concentration and the slow cycle primarily in the open/super-open state at a low concentration. The weighted ensemble simulations revealed the kinetics of the structural changes in GCK with the synergetic use of the MSES results; the rate constant of the transition between the closed state and the open/super-open states, kC/O = 1.1 ms-1, is on the same order as the experimental catalytic rate, kcat = 0.22 ms-1. Finally, we discuss the pharmacological activities of GCK activators (small molecular drugs modulating the GCK activity) in terms of the slight changes in the domain motion, depending on their chemical structures as regulators. The present study demonstrated the capability of the enhanced sampling and the associated kinetic calculations for understanding the atomistic structural dynamics of protein systems in physiological environments.

Novel approach of fragment-based lead discovery applied to renin inhibitors.

A novel approach was conducted for fragment-based lead discovery and applied to renin inhibitors. The biochemical screening of a fragment library against renin provided the hit fragment which showed a characteristic interaction pattern with the target protein. The hit fragment bound only to the S1, S3, and S3SP (S3 subpocket) sites without any interactions with the catalytic aspartate residues (Asp32 and Asp215 (pepsin numbering)). Prior to making chemical modifications to the hit fragment, we first identified its essential binding sites by utilizing the hit fragment's substructures. Second, we created a new and smaller scaffold, which better occupied the identified essential S3 and S3SP sites, by utilizing library synthesis with high-throughput chemistry. We then revisited the S1 site and efficiently explored a good building block attaching to the scaffold with library synthesis. In the library syntheses, the binding modes of each pivotal compound were determined and confirmed by X-ray crystallography and the library was strategically designed by structure-based computational approach not only to obtain a more active compound but also to obtain informative Structure Activity Relationship (SAR). As a result, we obtained a lead compound offering synthetic accessibility as well as the improved in vitro ADMET profiles. The fragments and compounds possessing a characteristic interaction pattern provided new structural insights into renin's active site and the potential to create a new generation of renin inhibitors. In addition, we demonstrated our FBDD strategy integrating highly sensitive biochemical assay, X-ray crystallography, and high-throughput synthesis and in silico library design aimed at fragment morphing at the initial stage was effective to elucidate a pocket profile and a promising lead compound.

Two-dimensional replica-exchange method for predicting protein-ligand binding structures.

We have developed a two-dimensional replica-exchange method for the prediction of protein-ligand binding structures. The first dimension is the umbrella sampling along the reaction coordinate, which is the distance between a protein binding pocket and a ligand. The second dimension is the solute tempering, in which the interaction between a ligand and a protein and water is weakened. The second dimension is introduced to make a ligand follow the umbrella potential more easily and enhance the binding events, which should improve the sampling efficiency. As test cases, we applied our method to two protein-ligand complex systems (MDM2 and HSP 90-alpha). Starting from the configuration in which the protein and the ligand are far away from each other in each system, our method predicted the ligand binding structures in excellent agreement with the experimental data from Protein Data Bank much faster with the improved sampling efficiency than the replica-exchange umbrella sampling method that we have previously developed.

Design and synthesis of novel pyrimido[4,5-b]azepine derivatives as HER2/EGFR dual inhibitors.

A novel 7,6 fused bicyclic scaffold, pyrimido[4,5-b]azepine was designed to fit into the ATP binding site of the HER2/EGFR proteins. The synthesis of this scaffold was accomplished by an intramolecular Claisen-type condensation. As the results of optimization lead us to 4-anilino and 6-functional groups, we discovered 6-substituted amide derivative 19b, which has a 1-benzothiophen-4-yloxy group attached to the 4-anilino group. An X-ray co-crystal structure of 19b with EGFR demonstrated that the N-1 and N-3 nitrogens of the pyrimido[4,5-b]azepine scaffold make hydrogen-bonding interactions with the main chain NH of Met793 and the side chain of Thr854 via a water-mediated hydrogen bond network, respectively. In addition, the NH proton at the 9-position makes an additional hydrogen bond with the carbonyl group of Met793, as we expected. Compound 19b revealed potent HER2/EGFR kinase (IC50: 24/36 nM) and BT474 cell growth (GI50: 18 nM) inhibitory activities based on its pseudo-irreversible (PI) profile.

Synthetic studies of centromere-associated protein-E (CENP-E) inhibitors: 1.Exploration of fused bicyclic core scaffolds using electrostatic potential map.

Centromere-associated protein-E (CENP-E), a mitotic kinesin that plays an important role in mitotic progression, is an attractive target for cancer therapeutic drugs. For the purpose of developing novel CENP-E inhibitors as cancer therapeutics, we investigated a fused bicyclic compound identified by high throughput screening, 4-oxo-4,5-dihydrothieno[3,4-c]pyridine-6-carboxamide 1a. Based on this scaffold, we designed inhibitors for efficient binding at the L5 site in CENP-E utilizing homology modeling as well as electrostatic potential map (EPM) analysis to enhance CENP-E inhibitory activity. This resulted in a new lead, 5-bromoimidazo[1,2-a]pyridine 7, which showed potent CENP-E enzyme inhibition (IC50: 50nM) and cellular activity with accumulation of phosphorylated histone H3 in HeLa cells. Our homology model and EPM analysis proved to be useful tools for the rational design of CENP-E inhibitors.

Structural analysis of the mechanism of inhibition and allosteric activation of the kinase domain of HER2 protein.

Aberrant signaling of ErbB family members human epidermal growth factor 2 (HER2) and epidermal growth factor receptor (EGFR) is implicated in many human cancers, and HER2 expression is predictive of human disease recurrence and prognosis. Small molecule kinase inhibitors of EGFR and of both HER2 and EGFR have received approval for the treatment of cancer. We present the first high resolution crystal structure of the kinase domain of HER2 in complex with a selective inhibitor to understand protein activation, inhibition, and function at the molecular level. HER2 kinase domain crystallizes as a dimer and suggests evidence for an allosteric mechanism of activation comparable with previously reported activation mechanisms for EGFR and HER4. A unique Gly-rich region in HER2 following the α-helix C is responsible for increased conformational flexibility within the active site and could explain the low intrinsic catalytic activity previously reported for HER2. In addition, we solved the crystal structure of the kinase domain of EGFR in complex with a HER2/EGFR dual inhibitor (TAK-285). Comparison with previously reported inactive and active EGFR kinase domain structures gave insight into the mechanism of HER2 and EGFR inhibition and may help guide the design and development of new cancer drugs with improved potency and selectivity.

Design and synthesis of pyrrolo[3,2-d]pyrimidine HER2/EGFR dual inhibitors: improvement of the physicochemical and pharmacokinetic profiles for potent in vivo anti-tumor efficacy.

During the course of our studies on a novel HER2/EGFR dual inhibitor (TAK-285), we found an alternative potent pyrrolo[3,2-d]pyrimidine compound (1a). To enhance the pharmacokinetic (PK) profile of this compound, we conducted chemical modifications into its N-5 side chain and conversion of the chemically modified compounds into their salts. Among them, 2cb, the tosylate salt of compound 2c, showed potent HER2/EGFR kinase inhibitory activity (IC(50): 11/11 nM) and cellular growth inhibitory activity (BT-474 cell GI(50): 56 nM) with a good drug metabolism and PK (DMPK) profile. Furthermore, 2cb exhibited significant in vivo antitumor efficacy in both mouse and rat xenograft models with transplanted 4-1ST gastric cancer cell lines (mouse, T/C=0%, 2cb po bid at 100 mg/kg; rat, T/C: -1%, 2cb po bid at 25 mg/kg).

A new class of non-thiazolidinedione, non-carboxylic-acid-based highly selective peroxisome proliferator-activated receptor (PPAR) γ agonists: design and synthesis of benzylpyrazole acylsulfonamides.

Herein, we describe the design, synthesis, and structure-activity relationships of novel benzylpyrazole acylsulfonamides as non-thiazolidinedione (TZD), non-carboxylic-acid-based peroxisome proliferator-activated receptor (PPAR) γ agonists. Docking model analysis of in-house weak agonist 2 bound to the reported PPARγ ligand binding domain suggested that modification of the carboxylic acid of 2 would help strengthen the interaction of 2 with the TZD pocket and afford non-carboxylic-acid-based agonists. In this study, we used an acylsulfonamide group as the ring-opening analog of TZD as an isosteric replacement of carboxylic acid moiety of 2; further, preliminary modification of the terminal alkyl chain on the sulfonyl group gave the lead compound 3c. Subsequent optimization of the resulting compound gave the potent agonists 25c, 30b, and 30c with high metabolic stability and significant antidiabetic activity. Further, we have described the difference in binding mode of the carboxylic-acid-based agonist 1 and acylsulfonamide 3d.

Structure-activity relationships and key structural feature of pyridyloxybenzene-acylsulfonamides as new, potent, and selective peroxisome proliferator-activated receptor (PPAR) γ Agonists.

In our search for a novel class of non-TZD, non-carboxylic acid peroxisome proliferator-activated receptor (PPAR) γ agonists, we explored alternative lipophilic templates to replace benzylpyrazole core of the previously reported agonist 1. Introduction of a pentylsulfonamide group into arylpropionic acids derived from previous in-house PPARγ ligands succeeded in the identification of 2-pyridyloxybenzene-acylsulfonamide 2 as a lead compound. Docking studies of compound 2 suggested that a substituent para to the central benzene ring should be incorporated to effectively fill the Y-shaped cavity of the PPARγ ligand-binding domain (LBD). This strategy led to significant improvement of PPARγ activity. Further optimization to balance in vitro activity and metabolic stability allowed the discovery of the potent, selective and orally efficacious PPARγ agonist 8f. Structure-activity relationship study as well as detailed analysis of the binding mode of 8f to the PPARγ-LBD revealed the essential structural features of this series of ligands.

Ab initio prediction of protein-ligand binding structures by replica-exchange umbrella sampling simulations.

We have developed a prediction method for the binding structures of ligands with proteins. Our method consists of three steps. First, replica-exchange umbrella sampling simulations are performed along the distance between a putative binding site of a protein and a ligand as the reaction coordinate. Second, we obtain the potential of mean force (PMF) of the unbiased system using the weighted histogram analysis method and determine the distance that corresponds to the global minimum of PMF. Third, structures that have this global-minimum distance and energy values around the average potential energy are collected and analyzed using the principal component analysis. We predict the binding structure as the global-minimum free energy state on the free energy landscapes along the two major principal component axes. As test cases, we applied our method to five protein-ligand complex systems. Starting from the configuration in which the protein and the ligand are far away from each other in each system, our method predicted the ligand binding structures in excellent agreement with the experimental data from Protein Data Bank.

17,20-lyase inhibitors. Part 4: design, synthesis and structure-activity relationships of naphthylmethylimidazole derivatives as novel 17,20-lyase inhibitors.

A novel series of naphthylmethylimidazole derivatives and related compounds have been investigated as selective 17,20-lyase inhibitors. Optimization of the substituent at the 6-position on the naphthalene ring was performed to yield a methylcarbamoyl derivative, which exhibited potent inhibitory activity against human 17,20-lyase and promising selectivity (>200-fold) for 17,20-lyase over CYP3A4. Further modifications of the methylcarbamoyl derivative led to the discovery of the corresponding tricyclic compound, which showed highly potent activity against human 17,20-lyase (IC(50) 19 nM) and good selectivity (>1000-fold) for inhibition of 17,20-lyase over CYP3A4. Additional biological evaluation revealed that the tricyclic compound had potent in vivo efficacy in monkeys and favorable pharmacokinetic profiles when administered in rats. Asymmetric synthesis of the selective tricyclic inhibitor was also achieved using a chiral α-hydroxy ketone.

17,20-Lyase inhibitors. Part 3: Design, synthesis, and structure-activity relationships of biphenylylmethylimidazole derivatives as novel 17,20-lyase inhibitors.

A novel series of biphenylylmethylimidazole derivatives and related compounds were synthesized as inhibitors of 17,20-lyase, a key enzyme in the production of steroid hormones, and their biological activities were evaluated. In an attempt to identify potent and selective inhibitors of 17,20-lyase over the related CYP3A4 enzyme, a homology model for human 17,20-lyase was developed using the X-ray crystallographic structure of the mammalian CYP2C5 enzyme. With the aid of molecular modeling, optimization of the biphenyl moiety was performed to give an acetamide derivative, which was resolved by HPLC to give the active (-)-enantiomer. The obtained active enantiomer showed not only potent inhibition of both rat and human 17,20-lyase,with IC(50) values of 14 and 26 nM, respectively, but also excellent selectivity (>300-fold) for inhibition of 17,20-lyase over CYP3A4. Moreover, the active enantiomer significantly reduced both serum testosterone and DHEA concentrations in a monkey model after single oral administration. Asymmetric synthesis of the active enantiomer was also developed via a chiral intermediate using a diastereoselective Grignard reaction.

Discovery of orteronel (TAK-700), a naphthylmethylimidazole derivative, as a highly selective 17,20-lyase inhibitor with potential utility in the treatment of prostate cancer.

A novel naphthylmethylimidazole derivative 1 and its related compounds were identified as 17,20-lyase inhibitors. Based on the structure-activity relationship around the naphthalene scaffold and the results of a docking study of 1a in the homology model of 17,20-lyase, the 6,7-dihydro-5H-pyrrolo[1,2-c]imidazole derivative (+)-3c was synthesized and identified as a potent and highly selective 17,20-lyase inhibitor. Biological evaluation of (+)-3c at a dose of 1mg/kg in a male monkey model revealed marked reductions in both serum testosterone and dehydroepiandrosterone concentrations. Therefore, (+)-3c (termed orteronel [TAK-700]) was selected as a candidate for clinical evaluation and is currently in phase III clinical trials for the treatment of castration-resistant prostate cancer.

Dependency of ligand free energy landscapes on charge parameters and solvent models.

We performed replica-exchange molecular dynamics (REMD) simulations of six ligands to examine the dependency of their free energy landscapes on charge parameters and solvent models. Six different charge parameter sets for each ligand were first generated by RESP and AM1-BCC methods using three different conformations independently. RESP charges showed some conformational dependency. On the other hand, AM1-BCC charges did not show conformational dependency and well reproduced the overall trend of RESP charges. The free energy landscapes obtained from the REMD simulations of ligands in vacuum, Generalized-Born (GB), and TIP3P solutions were then analyzed. We found that even small charge differences can produce qualitatively different landscapes in vacuum condition, but the differences tend to be much smaller under GB and TIP3P conditions. The simulations in the GB model well reproduced the landscapes in the TIP3P model using only a fraction of the computational cost. The protein-bound ligand conformations were rarely the global minimum states, but similar conformations were found to exist in aqueous solution without proteins in regions close to the global minimum, local minimum or intermediate states.

Design, synthesis and structure-activity relationships of 1,3,4-oxadiazole derivatives as novel inhibitors of glycogen synthase kinase-3beta.

Glycogen synthase kinase-3beta (GSK-3beta) is implicated in abnormal hyperphosphorylation of tau protein and its inhibitors are expected to be a promising therapeutic agents for the treatment of Alzheimer's disease. Here we report design, synthesis and structure-activity relationships of a novel series of oxadiazole derivatives as GSK-3beta inhibitors. Among these inhibitors, compound 20x showed highly selective and potent GSK-3beta inhibitory activity in vitro and its binding mode was determined by obtaining the X-ray co-crystal structure of 20x and GSK-3beta.

Cesium carbonate-catalyzed alpha-phenylchalcogenation of carbonyl compounds with diphenyl dichalcogenide.

It was found that cesium carbonate has a unique catalytic ability on the reaction of carbonyl compounds with diphenyl diselenide to give the corresponding alpha-phenylseleno carbonyl compounds in moderate to good yields. Similarly, the alpha-phenylthiolation of carbonyl compounds with diphenyl disulfide was promoted by the cesium carbonate catalyst.

Structure-Based Design, Synthesis, and Biological Evaluation of Imidazo[4,5-b]pyridin-2-one-Based p38 MAP Kinase Inhibitors: Part 1.

We identified a lead series of p38 mitogen-activated protein kinase inhibitors using a structure-based design strategy from high-throughput screening of hit compound 1. X-ray crystallography of 1 with the kinase showed an infrequent flip of the peptide bond between Met109 and Gly110, which was considered to lead to high kinase selectivity. Our structure-based design strategy was to conduct scaffold transformation of 1 with maintenance of hydrogen bond interactions with the flipped hinge backbone of the enzyme. In accordance with this strategy, we focused on scaffold transformation to identify imidazo[4,5-b]pyridin-2-one derivatives as potent inhibitors of the p38 MAP kinase. Of the compounds evaluated, 21 was found to be a potent inhibitor of the p38 MAP kinase, lipopolysaccharide-induced tumor necrosis factor-α (TNF-α) production in human monocytic leukemia cells, and TNF-α-induced production of interleukin-8 in human whole blood cells. Herein we describe the discovery of potent and orally bioavailable imidazo[4,5-b]pyridin-2-one-based p38 MAP kinase inhibitors that suppressed cytokine production in a human whole blood cell-based assay.

Structure-Based Design, Synthesis, and Biological Evaluation of Imidazo[4,5-b]Pyridin-2-one-Based p38 MAP Kinase Inhibitors: Part†2.

We identified novel potent inhibitors of p38 mitogen-activated protein (MAP) kinase using a structure-based design strategy, beginning with lead compound, 3-(butan-2-yl)-6-(2,4-difluoroanilino)-1,3-dihydro-2H-imidazo[4,5-b]pyridin-2-one (1). To enhance the inhibitory activity of 1 against production of tumor necrosis factor-α (TNF-α) in human whole blood (hWB) cell assays, we designed and synthesized hybrid compounds in which the imidazo[4,5-b]pyridin-2-one core was successfully linked with the p-methylbenzamide fragment. Among the compounds evaluated, 3-(3-tert-butyl-2-oxo-2,3-dihydro-1H-imidazo[4,5-b]pyridin-6-yl)-4-methyl-N-(1-methyl-1H-pyrazol-3-yl)benzamide (25) exhibited potent p38 inhibition, superior suppression of TNF-α production in hWB cells, and also significant in vivo efficacy in a rat model of collagen-induced arthritis (CIA). In this paper, we report the discovery of potent, selective, and orally bioavailable imidazo[4,5-b]pyridin-2-one-based p38 MAP kinase inhibitors.

Discovery of imidazo[1,5-c]imidazol-3-ones: weakly basic, orally active factor Xa inhibitors.

The coagulation enzyme factor Xa (FXa) has been recognized as a promising target for the development of new antithrombotic agents. We previously found compound 1 to be an orally bioavailable FXa inhibitor in fasted monkeys; however, 1 showed poor bioavailability in rats and fed monkeys. To work out the pharmacokinetic problems, we focused our synthetic efforts on the chemical conversion of the 4-(imidazo[1,2- a]pyridin-5-yl)piperazine moiety of 1 to imidazolylpiperidine derivatives (fused and nonfused), which resulted in the discovery of the weakly basic imidazo[1,5- c]imidazol-3-one 3q as a potent and selective FXa inhibitor. Compound 3q showed favorable oral bioavailability in rats and monkeys under both fasted and fed conditions and antithrombotic efficacy in a rat model of venous thrombosis after oral administration, without a significant increase in bleeding time (unlike warfarin). On the basis of these promising properties, compound 3q was selected for further evaluation.