Discovery and characterization of highly potent and selective allosteric USP7 inhibitors.

Given the importance of ubiquitin-specific protease 7 (USP7) in oncogenic pathways, identification of USP7 inhibitors has attracted considerable interest. Despite substantial efforts, however, the development of validated deubiquitinase (DUB) inhibitors that exhibit drug-like properties and a well-defined mechanism of action has proven particularly challenging. In this article, we describe the identification, optimization and detailed characterization of highly potent (IC50 < 10 nM), selective USP7 inhibitors together with their less active, enantiomeric counterparts. We also disclose, for the first time, co-crystal structures of a human DUB enzyme complexed with small-molecule inhibitors, which reveal a previously undisclosed allosteric binding site. Finally, we report the identification of cancer cell lines hypersensitive to USP7 inhibition (EC50 < 30 nM) and demonstrate equal or superior activity in these cell models compared to clinically relevant MDM2 antagonists. Overall, these findings demonstrate the tractability and druggability of DUBs, and provide important tools for additional target validation studies.

Life-history variation along environmental and harvest clines of a northern freshwater fish: Plasticity and adaptation.

Plasticity, local adaptation and evolutionary trade-offs drive clinal variation in traits associated with lifetime growth. Disentangling the processes and determinants that cause these traits to vary helps to understand species' responses to changing environments. This is particularly urgent for exploited populations, where size-selective harvest can induce life-history evolution. Lake trout (Salvelinus namaycush) are an exploited fish with a life history adapted to low-productivity freshwaters of northern North America, which makes them highly vulnerable to ecosystem changes and overfishing. We characterized life-history variation across a broad and diverse landscape for this iconic northern freshwater fish and evaluated whether clinal variation was consistent with hypotheses for local adaptation or growth plasticity. We estimated growth-associated traits for 90 populations exposed to a diversity of environments using a Bayesian multivariate hierarchical model. We tested for clinal variation in their somatic growth, size at maturity and reproductive allocation along environmental gradients of lake productivity, climate, prey and exploitation clines under competing hypotheses of plasticity and local adaptation. Clinal life-history variation was consistent with growth plasticity and local adaptations but not harvest-induced evolution. Variation in somatic growth was explained by exploitation, climate and prey fish occurrence. Increased exploitation, from pristine to fully exploited conditions, led to increased somatic growth (from 32 to 45 mm/year) and adult life spans, and reduced age at maturity (from 11 to 8 years). Variation in size at maturity was explained by climate and, less certainly, prey fish occurrence, while reproductive allocation was explained by evolutionary trade-offs with mortality and other traits, but not environment. Lake trout life-history variation within this range was as wide as that observed across dozens of other freshwater species. Lake trout life histories resulted from evolutionary trade-offs, growth plasticity and local adaptations along several environmental clines. Presuming a plastic response, we documented ~1.4-fold growth compensation to exploitation-lower growth compensation than observed in many freshwater fishes. These results suggest that harvested species exposed to spatially structured and diverse environments may have substantial clinal variation on different traits, but due to different processes, and this has implications for their resilience and successful management.

Toward an understanding of agonist binding to human Orexin-1 and Orexin-2 receptors with G-protein-coupled receptor modeling and site-directed mutagenesis.

The class A G-protein-coupled receptors (GPCRs) Orexin-1 (OX1) and Orexin-2 (OX2) are located predominantly in the brain and are linked to a range of different physiological functions, including the control of feeding, energy metabolism, modulation of neuro-endocrine function, and regulation of the sleep-wake cycle. The natural agonists for OX1 and OX2 are two neuropeptides, Orexin-A and Orexin-B, which have activity at both receptors. Site-directed mutagenesis (SDM) has been reported on both the receptors and the peptides and has provided important insight into key features responsible for agonist activity. However, the structural interpretation of how these data are linked together is still lacking. In this work, we produced and used SDM data, homology modeling followed by MD simulation, and ensemble-flexible docking to generate binding poses of the Orexin peptides in the OX receptors to rationalize the SDM data. We also developed a protein pairwise similarity comparing method (ProS) and a GPCR-likeness assessment score (GLAS) to explore the structural data generated within a molecular dynamics simulation and to help distinguish between different GPCR substates. The results demonstrate how these newly developed methods of structural assessment for GPCRs can be used to provide a working model of neuropeptide-Orexin receptor interaction.

Fighting obesity with a sugar-based library: discovery of novel MCH-1R antagonists by a new computational-VAST approach for exploration of GPCR binding sites.

Obesity is an increasingly common disease. While antagonism of the melanin-concentrating hormone-1 receptor (MCH-1R) has been widely reported as a promising therapeutic avenue for obesity treatment, no MCH-1R antagonists have reached the market. Discovery and optimization of new chemical matter targeting MCH-1R is hindered by reduced HTS success rates and a lack of structural information about the MCH-1R binding site. X-ray crystallography and NMR, the major experimental sources of structural information, are very slow processes for membrane proteins and are not currently feasible for every GPCR or GPCR-ligand complex. This situation significantly limits the ability of these methods to impact the drug discovery process for GPCR targets in "real-time", and hence, there is an urgent need for other practical and cost-efficient alternatives. We present here a conceptually pioneering approach that integrates GPCR modeling with design, synthesis, and screening of a diverse library of sugar-based compounds from the VAST technology (versatile assembly on stable templates) to provide structural insights on the MCH-1R binding site. This approach creates a cost-efficient new avenue for structure-based drug discovery (SBDD) against GPCR targets. In our work, a primary VAST hit was used to construct a high-quality MCH-1R model. Following model validation, a structure-based virtual screen yielded a 14% hit rate and 10 novel chemotypes of potent MCH-1R antagonists, including EOAI3367472 (IC50 = 131 nM) and EOAI3367474 (IC50 = 213 nM).

Study of human Orexin-1 and -2 G-protein-coupled receptors with novel and published antagonists by modeling, molecular dynamics simulations, and site-directed mutagenesis.

The class A G-protein-coupled receptors (GPCRs) Orexin-1 (OX1) and Orexin-2 (OX2) are located predominantly in the brain and are linked to a range of different physiological functions, including the control of feeding, energy metabolism, modulation of neuro-endocrine function, and regulation of the sleep-wake cycle. Site-directed mutagenesis (SDM) and domain exchange (chimera) studies have provided important insight into key features of the OX1 and OX2 binding sites. However, the precise determinants of antagonist binding and selectivity are still not fully known. In this work, we used homology modeling of OX receptors to direct further SDM studies. These SDM studies were followed by molecular dynamics (MD) simulations to rationalize the full scope of the SDM data and to explain the role of each mutated residue in the binding and selectivity of a set of OX antagonists: Almorexant (dual OX1 and OX2 antagonist), SB-674042 (OX1 selective antagonist), EMPA (OX2 selective antagonist), and others. Our primary interest was focused on transmembrane helix 3 (TM3), which is identified as being of great importance for the selectivity of OX antagonists. These studies revealed conformational differences between the TM3 helices of OX1 and OX2, resulting from differences in amino acid sequences of the OX receptors that affect key interhelical interactions formed between TM3 and neighboring TM domains. The MD simulation protocol used here, which was followed by flexible docking studies, went beyond the use of static models and allowed for a more detailed exploration of the OX structures. In this work, we have demonstrated how even small differences in the amino acid sequences of GPCRs can lead to significant differences in structure, antagonist binding affinity, and selectivity of these receptors. The MD simulations allowed refinement of the OX receptor models to a degree that was not possible with static homology modeling alone and provided a deeper rationalization of the SDM data obtained. To validate these findings and to demonstrate that they can be usefully applied to the design of novel, very selective OX antagonists, we show here two examples of antagonists designed in house: EP-109-0092 (OX1 selective) and EP-009-0513 (OX2 selective).

Identification and development of a subtype-selective allosteric AKT inhibitor suitable for clinical development.

The serine/threonine protein kinase AKT plays a pivotal role within the PI3K pathway in regulating cellular proliferation and apoptotic cellular functions, and AKT hyper-activation via gene amplification and/or mutation has been implicated in multiple human malignancies. There are 3 AKT isoenzymes (AKT1-3) which mediate critical, non-redundant functions. We present the discovery and development of ALM301, a novel, allosteric, sub-type selective inhibitor of AKT1/2. ALM301 binds in an allosteric pocket created by the combined movement of the PH domain and the catalytic domain, resulting in a DFG out conformation. ALM301 was shown to be highly selective against a panel of over 450 kinases and potently inhibited cellular proliferation. These effects were particularly pronounced in MCF-7 cells containing a PI3KCA mutation. Subsequent cellular downstream pathway analysis in this sensitive cell line revealed potent inhibition of pAKT signalling up to 48 h post dosing. ALM301 treatment was well tolerated in an MCF-7 xenograft model and led to a dose-dependent reduction in tumour growth. Enhanced efficacy was observed in combination with tamoxifen. In summary, ALM301 is a highly specific AKT 1/2 inhibitor with an excellent pharmacological profile suitable for further clinical development.

The synthesis and evaluation of indolylureas as PKCα inhibitors.

PKCα and PKA have crucial but opposing roles in the regulation of calcium handling within myocytes. Identification of compounds that inhibit PKCα, but not PKA, are potential therapeutic targets for the treatment of heart disease. The synthesis of indolylureas are described, and a compound displaying nanomolar inhibition towards PKCα with significant selectivity over PKA has been identified.

DNA gyrase (GyrB)/topoisomerase IV (ParE) inhibitors: synthesis and antibacterial activity.

The synthesis and antibacterial activities of three chemotypes of DNA supercoiling inhibitors based on imidazolo[1,2-a]pyridine and [1,2,4]triazolo[1,5-a]pyridine scaffolds that target the ATPase subunits of DNA gyrase and topoisomerase IV (GyrB/ParE) is reported. The most potent scaffold was selected for optimization leading to a series with potent Gram-positive antibacterial activity and a low resistance frequency.

The multiple roles of computational chemistry in fragment-based drug design.

Fragment-based drug discovery (FBDD) represents a change in strategy from the screening of molecules with higher molecular weights and physical properties more akin to fully drug-like compounds, to the screening of smaller, less complex molecules. This is because it has been recognised that fragment hit molecules can be efficiently grown and optimised into leads, particularly after the binding mode to the target protein has been first determined by 3D structural elucidation, e.g. by NMR or X-ray crystallography. Several studies have shown that medicinal chemistry optimisation of an already drug-like hit or lead compound can result in a final compound with too high molecular weight and lipophilicity. The evolution of a lower molecular weight fragment hit therefore represents an attractive alternative approach to optimisation as it allows better control of compound properties. Computational chemistry can play an important role both prior to a fragment screen, in producing a target focussed fragment library, and post-screening in the evolution of a drug-like molecule from a fragment hit, both with and without the available fragment-target co-complex structure. We will review many of the current developments in the area and illustrate with some recent examples from successful FBDD discovery projects that we have conducted.

Identification and Structure-Guided Development of Pyrimidinone Based USP7 Inhibitors.

Ubiquitin specific protease 7 (USP7, HAUSP) has become an attractive target in drug discovery due to the role it plays in modulating Mdm2 levels and consequently p53. Increasing interest in USP7 is emerging due to its potential involvement in oncogenic pathways as well as possible roles in both metabolic and immune disorders in addition to viral infections. Potent, novel, and selective inhibitors of USP7 have been developed using both rational and structure-guided design enabled by high-resolution cocrystallography. Initial hits were identified via fragment-based screening, scaffold-hopping, and hybridization exercises. Two distinct subseries are described along with associated structure-activity relationship trends, as are initial efforts aimed at developing compounds suitable for in vivo experiments. Overall, these discoveries will enable further research into the wider biological role of USP7.

Binding mode and structure-activity relationships around direct inhibitors of the Nrf2-Keap1 complex.

An X-ray crystal structure of Kelch-like ECH-associated protein (Keap1) co-crystallised with (1S,2R)-2-[(1S)-1-[(1,3-dioxo-2,3-dihydro-1H-isoindol-2-yl)methyl]-1,2,3,4-tetrahydroisoquinolin-2-carbonyl]cyclohexane-1-carboxylic acid (compound (S,R,S)-1 a) was obtained. This X-ray crystal structure provides breakthrough experimental evidence for the true binding mode of the hit compound (S,R,S)-1 a, as the ligand orientation was found to differ from that of the initial docking model, which was available at the start of the project. Crystallographic elucidation of this binding mode helped to focus and drive the drug design process more effectively and efficiently.

Fragment-based identification of Hsp90 inhibitors.

Heat shock protein 90 (Hsp90) plays a key role in stress response and protection of the cell against the effects of mutation. Herein we report the identification of an Hsp90 inhibitor identified by fragment screening using a high-concentration biochemical assay, as well as its optimisation by in silico searching coupled with a structure-based drug design (SBDD) approach.