Rational Design of Highly Potent, Selective, and Bioavailable SGK1 Protein Kinase Inhibitors for the Treatment of Osteoarthritis.

The serine/threonine kinase SGK1 is an activator of the ß-catenin pathway and a powerful stimulator of cartilage degradation that is found to be upregulated under genomic control in diseased osteoarthritic cartilage. Today, no oral disease-modifying treatments are available and chronic treatment in this indication sets high requirements for the drug selectivity, pharmacokinetic, and safety profile. We describe the identification of a highly selective druglike 1H-pyrazolo[3,4-d]pyrimidine SGK1 inhibitor 17a that matches both safety and pharmacokinetic requirements for oral dosing. Rational compound design was facilitated by a novel hSGK1 co-crystal structure, and multiple ligand-based computer models were applied to guide the chemical optimization of the compound ADMET and selectivity profiles. Compounds were selected for subchronic proof of mechanism studies in the mouse femoral head cartilage explant model, and compound 17a emerged as a druglike SGK1 inhibitor, with a highly optimized profile suitable for oral dosing as a novel, potentially disease-modifying agent for osteoarthritis.

Estimation of Protein-Ligand Unbinding Kinetics Using Non-Equilibrium Targeted Molecular Dynamics Simulations.

We here report on nonequilibrium targeted molecular dynamics simulations as a tool for the estimation of protein-ligand unbinding kinetics. Correlating simulations with experimental data from SPR kinetics measurements and X-ray crystallography on two small molecule compound libraries bound to the N-terminal domain of the chaperone Hsp90, we show that the mean nonequilibrium work computed in an ensemble of trajectories of enforced ligand unbinding is a promising predictor for ligand unbinding rates. We furthermore investigate the molecular basis determining unbinding rates within the compound libraries. We propose ligand conformational changes and protein-ligand nonbonded interactions to impact on unbinding rates. Ligands may remain longer at the protein if they exhibit strong electrostatic and/or van der Waals interactions with the target. In the case of ligands with a rigid chemical scaffold that exhibit longer residence times, transient electrostatic interactions with the protein appear to facilitate unbinding. Our results imply that understanding the unbinding pathway and the protein-ligand interactions along this path is crucial for the prediction of small molecule ligands with defined unbinding kinetics.

The signalling conformation of the insulin receptor ectodomain.

Understanding the structural biology of the insulin receptor and how it signals is of key importance in the development of insulin analogs to treat diabetes. We report here a cryo-electron microscopy structure of a single insulin bound to a physiologically relevant, high-affinity version of the receptor ectodomain, the latter generated through attachment of C-terminal leucine zipper elements to overcome the conformational flexibility associated with ectodomain truncation. The resolution of the cryo-electron microscopy maps is 3.2 Å in the insulin-binding region and 4.2 Å in the membrane-proximal region. The structure reveals how the membrane proximal domains of the receptor come together to effect signalling and how insulin's negative cooperativity of binding likely arises. Our structure further provides insight into the high affinity of certain super-mitogenic insulins. Together, these findings provide a new platform for insulin analog investigation and design.

Estimation of Drug-Target Residence Times by τ-Random Acceleration Molecular Dynamics Simulations.

Drug-target residence time (τ), one of the main determinants of drug efficacy, remains highly challenging to predict computationally and, therefore, is usually not considered in the early stages of drug design. Here, we present an efficient computational method, τ-random acceleration molecular dynamics (τRAMD), for the ranking of drug candidates by their residence time and obtaining insights into ligand-target dissociation mechanisms. We assessed τRAMD on a data set of 70 diverse drug-like ligands of the N-terminal domain of HSP90α, a pharmaceutically important target with a highly flexible binding site, obtaining computed relative residence times with an accuracy of about 2.3τ for 78% of the compounds and less than 2.0τ within congeneric series. Analysis of dissociation trajectories reveals features that affect ligand unbinding rates, including transient polar interactions and steric hindrance. These results suggest that τRAMD will be widely applicable as a computationally efficient aid to improving drug residence times during lead optimization.

Kinetics of human cannabinoid 1 (CB1) receptor antagonists: Structure-kinetics relationships (SKR) and implications for insurmountable antagonism.

While equilibrium binding affinities and in vitro functional antagonism of CB1 receptor antagonists have been studied in detail, little is known on the kinetics of their receptor interaction. In this study, we therefore conducted kinetic assays for nine 1-(4,5-diarylthiophene-2-carbonyl)-4-phenylpiperidine-4-carboxamide derivatives and included the CB1 antagonist rimonabant as a comparison. For this we newly developed a dual-point competition association assay with [3H]CP55940 as the radioligand. This assay yielded Kinetic Rate Index (KRI) values from which structure-kinetics relationships (SKR) of hCB1 receptor antagonists could be established. The fast dissociating antagonist 6 had a similar receptor residence time (RT) as rimonabant, i.e. 19 and 14 min, respectively, while the slowest dissociating antagonist (9) had a very long RT of 2222 min, i.e. pseudo-irreversible dissociation kinetics. In functional assays, 9 displayed insurmountable antagonism, while the effects of the shortest RT antagonist 6 and rimonabant were surmountable. Taken together, this study shows that hCB1 receptor antagonists can have very divergent RTs, which are not correlated to their equilibrium affinities. Furthermore, their RTs appear to define their mode of functional antagonism, i.e. surmountable vs. insurmountable. Finally, based on the recently resolved hCB1 receptor crystal structure, we propose that the differences in RT can be explained by a different binding mode of antagonist 9 from short RT antagonists that is able to displace unfavorable water molecules. Taken together, these findings are of importance for future design and evaluation of potent and safe hCB1 receptor antagonists.

Accounting for conformational variability in protein-ligand docking with NMR-guided rescoring.

A key component to success in structure-based drug design is reliable information on protein-ligand interactions. Recent development in NMR techniques has accelerated this process by overcoming some of the limitations of X-ray crystallography and computational protein-ligand docking. In this work we present a new scoring protocol based on NMR-derived interligand INPHARMA NOEs to guide the selection of computationally generated docking modes. We demonstrate the performance in a range of scenarios, encompassing traditionally difficult cases such as docking to homology models and ligand dependent domain rearrangements. Ambiguities associated with sparse experimental information are lifted by searching a consensus solution based on simultaneously fitting multiple ligand pairs. This study provides a previously unexplored integration between molecular modeling and experimental data, in which interligand NOEs represent the key element in the rescoring algorithm. The presented protocol should be widely applicable for protein-ligand docking also in a different context from drug design and highlights the important role of NMR-based approaches to describe intermolecular ligand-receptor interactions.

Influence of heparin mimetics on assembly of the FGF.FGFR4 signaling complex.

Fibroblast growth factor (FGF) signaling regulates mammalian development and metabolism, and its dysregulation is implicated in many inherited and acquired diseases, including cancer. Heparan sulfate glycosaminoglycans (HSGAGs) are essential for FGF signaling as they promote FGF.FGF receptor (FGFR) binding and dimerization. Using novel organic synthesis protocols to prepare homogeneously sulfated heparin mimetics (HM), including hexasaccharide (HM(6)), octasaccharide (HM(8)), and decasaccharide (HM(10)), we tested the ability of these HM to support FGF1 and FGF2 signaling through FGFR4. Biological assays show that both HM(8) and HM(10) are significantly more potent than HM(6) in promoting FGF2-mediated FGFR4 signaling. In contrast, all three HM have comparable activity in promoting FGF1.FGFR4 signaling. To understand the molecular basis for these differential activities in FGF1/2.FGFR4 signaling, we used NMR spectroscopy, isothermal titration calorimetry, and size-exclusion chromatography to characterize binding interactions of FGF1/2 with the isolated Ig-domain 2 (D2) of FGFR4 in the presence of HM, and binary interactions of FGFs and D2 with HM. Our data confirm the existence of both a secondary FGF1.FGFR4 interaction site and a direct FGFR4.FGFR4 interaction site thus supporting the formation of the symmetric mode of FGF.FGFR dimerization in solution. Moreover, our results show that the observed higher activity of HM(8) relative to HM(6) in stimulating FGF2.FGFR4 signaling correlates with the higher affinity of HM(8) to bind and dimerize FGF2. Notably FGF2.HM(8) exhibits pronounced positive binding cooperativity. Based on our findings we propose a refined symmetric FGF.FGFR dimerization model, which incorporates the differential ability of HM to dimerize FGFs.

Rational design of potent GSK3beta inhibitors with selectivity for Cdk1 and Cdk2.

From an HTS hit, a series of potent and selective inhibitors of GSK3beta have been designed based on a Cdk2-homology model and with the help of several crystal structures of the compounds within Cdk2.

Thermodynamic characterization of allosteric glycogen phosphorylase inhibitors.

Glycogen phosphorylase (GP) is a validated target for the treatment of type 2 diabetes. Here we describe highly potent GP inhibitors, AVE5688, AVE2865, and AVE9423. The first two compounds are optimized members of the acyl urea series. The latter represents a novel quinolone class of GP inhibitors, which is introduced in this study. In the enzyme assay, both inhibitor types compete with the physiological activator AMP and act synergistically with glucose. Isothermal titration calorimetry (ITC) shows that the compounds strongly bind to nonphosphorylated, inactive GP (GPb). Binding to phosphorylated, active GP (GPa) is substantially weaker, and the thermodynamic profile reflects a coupled transition to the inactive (tense) conformation. Crystal structures confirm that the three inhibitors bind to the AMP site of tense state GP. These data provide the first direct evidence that acyl urea and quinolone compounds are allosteric inhibitors that selectively bind to and stabilize the inactive conformation of the enzyme. Furthermore, ITC reveals markedly different thermodynamic contributions to inhibitor potency that can be related to the binding modes observed in the cocrystal structures. For AVE5688, which occupies only the lower part of the bifurcated AMP site, binding to GPb (Kd = 170 nM) is exclusively enthalpic (Delta H = -9.0 kcal/mol, TDelta S = 0.3 kcal/mol). The inhibitors AVE2865 (Kd = 9 nM, Delta H = -6.8 kcal/mol, TDelta S = 4.2 kcal/mol) and AVE9423 (Kd = 24 nM, Delta H = -5.9 kcal/mol, TDelta S = 4.6 kcal/mol) fully exploit the volume of the binding pocket. Their pronounced binding entropy can be attributed to the extensive displacement of solvent molecules as well as to ionic interactions with the phosphate recognition site.

APRV - a program for automated data processing, refinement and visualization.

APRV (Automatic Processing, Refinement and Visualization) is a new program that enables high-throughput batch processing of crystallographic data. The program combines processing of raw diffraction images, initial structure refinement and visual inspection of resulting electron density into a seamless one-step procedure, during which all relevant parameters are refined automatically. It is controlled by a user-friendly graphical interface, facilitating operation by non-experts.

Crystal structure of MJ1247 protein from M. jannaschii at 2.0 A resolution infers a molecular function of 3-hexulose-6-phosphate isomerase.

The crystal structure of the hypothetical protein MJ1247 from Methanococccus jannaschii at 2 A resolution, a detailed sequence analysis, and biochemical assays infer its molecular function to be 3-hexulose-6-phosphate isomerase (PHI). In the dissimilatory ribulose monophosphate (RuMP) cycle, ribulose-5-phosphate is coupled to formaldehyde by the 3-hexulose-6-phosphate synthase (HPS), yielding hexulose-6-phosphate, which is then isomerized to fructose-6-phosphate by the enzyme 3-hexulose-6-phosphate isomerase. MJ1247 is an alpha/beta structure consisting of a five-stranded parallel beta sheet flanked on both sides by alpha helices, forming a three-layered alpha-beta-alpha sandwich. The fold represents the nucleotide binding motif of a flavodoxin type. MJ1247 is a tetramer in the crystal and in solution and each monomer has a folding similar to the isomerase domain of glucosamine-6-phosphate synthase (GlmS).