An inverse agonist of orphan receptor GPR61 acts by a G protein-competitive allosteric mechanism.

GPR61 is an orphan GPCR related to biogenic amine receptors. Its association with phenotypes relating to appetite makes it of interest as a druggable target to treat disorders of metabolism and body weight, such as obesity and cachexia. To date, the lack of structural information or a known biological ligand or tool compound has hindered comprehensive efforts to study GPR61 structure and function. Here, we report a structural characterization of GPR61, in both its active-like complex with heterotrimeric G protein and in its inactive state. Moreover, we report the discovery of a potent and selective small-molecule inverse agonist against GPR61 and structural elucidation of its allosteric binding site and mode of action. These findings offer mechanistic insights into an orphan GPCR while providing both a structural framework and tool compound to support further studies of GPR61 function and modulation.

Small molecule branched-chain ketoacid dehydrogenase kinase (BDK) inhibitors with opposing effects on BDK protein levels.

Branched chain amino acid (BCAA) catabolic impairments have been implicated in several diseases. Branched chain ketoacid dehydrogenase (BCKDH) controls the rate limiting step in BCAA degradation, the activity of which is inhibited by BCKDH kinase (BDK)-mediated phosphorylation. Screening efforts to discover BDK inhibitors led to identification of thiophene PF-07208254, which improved cardiometabolic endpoints in mice. Structure-activity relationship studies led to identification of a thiazole series of BDK inhibitors; however, these inhibitors did not improve metabolism in mice upon chronic administration. While the thiophenes demonstrated sustained branched chain ketoacid (BCKA) lowering and reduced BDK protein levels, the thiazoles increased BCKAs and BDK protein levels. Thiazoles increased BDK proximity to BCKDH-E2, whereas thiophenes reduced BDK proximity to BCKDH-E2, which may promote BDK degradation. Thus, we describe two BDK inhibitor series that possess differing attributes regarding BDK degradation or stabilization and provide a mechanistic understanding of the desirable features of an effective BDK inhibitor.

Structural studies identify angiotensin II receptor blocker-like compounds as branched-chain ketoacid dehydrogenase kinase inhibitors.

The mammalian mitochondrial branched-chain ketoacid dehydrogenase (BCKD) complex is a multienzyme complex involved in the catabolism of branched-chain amino acids. BCKD is regulated by the BCKD kinase, or BCKDK, which binds to the E2 subunit of BCKD, phosphorylates its E1 subunit, and inhibits enzymatic activity. Inhibition of the BCKD complex results in increased levels of branched-chain amino acids and branched-chain ketoacids, and this buildup has been associated with heart failure, type 2 diabetes mellitus, and nonalcoholic fatty liver disease. To find BCKDK inhibitors for potential treatment of these diseases, we performed both NMR and virtual fragment screening and identified tetrazole-bearing fragments that bind BCKDK at multiple sites. Through structure-based virtual screening expanding from these fragments, the angiotensin receptor blocker class antihypertension drugs and angiotensin receptor blocker-like compounds were discovered to be potent BCKDK inhibitors, suggesting potential new avenues for heart failure treatment combining BCKDK inhibition and antihypertension.

NMR spectroscopy: the swiss army knife of drug discovery.

Nuclear magnetic resonance (NMR) spectroscopy has evolved into a powerful tool within drug discovery over the last two decades. While traditionally being used by medicinal chemists for small molecule structure elucidation, it can also be a valuable tool for the identification of small molecules that bind to drug targets, for the characterization of target-ligand interactions and for hit-to-lead optimization. Here, we describe how NMR spectroscopy is integrated into the Pfizer drug discovery pipeline and how we utilize this approach to identify and validate initial hits and generate leads.

Author Correction: Dopamine D3 receptor antagonist reveals a cryptic pocket in aminergic GPCRs.

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has not been fixed in the paper.

Late-Stage Microsomal Oxidation Reduces Drug-Drug Interaction and Identifies Phosphodiesterase 2A Inhibitor PF-06815189.

Late-stage oxidation using liver microsomes was applied to phosphodiesterase 2 inhibitor 1 to reduce its clearance by cytochrome P450 enzymes, introduce renal clearance, and minimize the risk for victim drug-drug interactions. This approach yielded PF-06815189 (2) with improved physicochemical properties and a mixed metabolic profile. This example highlights the importance of C-H diversification methods to drug discovery.

Identification of a Potent, Highly Selective, and Brain Penetrant Phosphodiesterase 2A Inhibitor Clinical Candidate.

Computational modeling was used to direct the synthesis of analogs of previously reported phosphodiesterase 2A (PDE2A) inhibitor 1 with an imidazotriazine core to yield compounds of significantly enhanced potency. The analog PF-05180999 (30) was subsequently identified as a preclinical candidate targeting cognitive impairment associated with schizophrenia. Compound 30 demonstrated potent binding to PDE2A in brain tissue, dose responsive mouse brain cGMP increases, and reversal of N-methyl-d-aspartate (NMDA) antagonist-induced (MK-801, ketamine) effects in electrophysiology and working memory models in rats. Preclinical pharmacokinetics revealed unbound brain/unbound plasma levels approaching unity and good oral bioavailability resulting in an average concentration at steady state (Cav,ss) predicted human dose of 30 mg once daily (q.d.). Modeling of a modified release formulation suggested that 25 mg twice daily (b.i.d.) could maintain plasma levels of 30 at or above targeted efficacious plasma levels for 24 h, which became part of the human clinical plan.

Dopamine D3 receptor antagonist reveals a cryptic pocket in aminergic GPCRs.

The recent increase in the number of X-ray crystal structures of G-protein coupled receptors (GPCRs) has been enabling for structure-based drug design (SBDD) efforts. These structures have revealed that GPCRs are highly dynamic macromolecules whose function is dependent on their intrinsic flexibility. Unfortunately, the use of static structures to understand ligand binding can potentially be misleading, especially in systems with an inherently high degree of conformational flexibility. Here, we show that docking a set of dopamine D3 receptor compounds into the existing eticlopride-bound dopamine D3 receptor (D3R) X-ray crystal structure resulted in poses that were not consistent with results obtained from site-directed mutagenesis experiments. We overcame the limitations of static docking by using large-scale high-throughput molecular dynamics (MD) simulations and Markov state models (MSMs) to determine an alternative pose consistent with the mutation data. The new pose maintains critical interactions observed in the D3R/eticlopride X-ray crystal structure and suggests that a cryptic pocket forms due to the shift of a highly conserved residue, F6.52. Our study highlights the importance of GPCR dynamics to understand ligand binding and provides new opportunities for drug discovery.

Structure-inspired design of ß-arrestin-biased ligands for aminergic GPCRs.

Development of biased ligands targeting G protein-coupled receptors (GPCRs) is a promising approach for current drug discovery. Although structure-based drug design of biased agonists remains challenging even with an abundance of GPCR crystal structures, we present an approach for translating GPCR structural data into ß-arrestin-biased ligands for aminergic GPCRs. We identified specific amino acid-ligand contacts at transmembrane helix 5 (TM5) and extracellular loop 2 (EL2) responsible for Gi/o and ß-arrestin signaling, respectively, and targeted those residues to develop biased ligands. For these ligands, we found that bias is conserved at other aminergic GPCRs that retain similar residues at TM5 and EL2. Our approach provides a template for generating arrestin-biased ligands by modifying predicted ligand interactions that block TM5 interactions and promote EL2 interactions. This strategy may facilitate the structure-guided design of arrestin-biased ligands at other GPCRs, including polypharmacological biased ligands.

Application of Structure-Based Design and Parallel Chemistry to Identify a Potent, Selective, and Brain Penetrant Phosphodiesterase 2A Inhibitor.

Phosphodiesterase 2A (PDE2A) inhibitors have been reported to demonstrate in vivo activity in preclinical models of cognition. To more fully explore the biology of PDE2A inhibition, we sought to identify potent PDE2A inhibitors with improved brain penetration as compared to current literature compounds. Applying estimated human dose calculations while simultaneously leveraging synthetically enabled chemistry and structure-based drug design has resulted in a highly potent, selective, brain penetrant compound 71 (PF-05085727) that effects in vivo biochemical changes commensurate with PDE2A inhibition along with behavioral and electrophysiological reversal of the effects of NMDA antagonists in rodents. This data supports the ability of PDE2A inhibitors to potentiate NMDA signaling and their further development for clinical cognition indications.

Dopamine D3/D2 Receptor Antagonist PF-4363467 Attenuates Opioid Drug-Seeking Behavior without Concomitant D2 Side Effects.

Dopamine receptor antagonism is a compelling molecular target for the treatment of a range of psychiatric disorders, including substance use disorders. From our corporate compound file, we identified a structurally unique D3 receptor (D3R) antagonist scaffold, 1. Through a hybrid approach, we merged key pharmacophore elements from 1 and D3 agonist 2 to yield the novel D3R/D2R antagonist PF-4363467 (3). Compound 3 was designed to possess CNS drug-like properties as defined by its CNS MPO desirability score (≥4/6). In addition to good physicochemical properties, 3 exhibited low nanomolar affinity for the D3R (D3 Ki = 3.1 nM), good subtype selectivity over D2R (D2 Ki = 692 nM), and high selectivity for D3R versus other biogenic amine receptors. In vivo, 3 dose-dependently attenuated opioid self-administration and opioid drug-seeking behavior in a rat operant reinstatement model using animals trained to self-administer fentanyl. Further, traditional extrapyramidal symptoms (EPS), adverse side effects arising from D2R antagonism, were not observed despite high D2 receptor occupancy (RO) in rodents, suggesting that compound 3 has a unique in vivo profile. Collectively, our data support further investigation of dual D3R and D2R antagonists for the treatment of drug addiction.

Identification of Phosphorylation Consensus Sequences and Endogenous Neuronal Substrates of the Psychiatric Risk Kinase TNIK.

Traf2- and Nck-interacting kinase (TNIK) is a serine/threonine kinase highly expressed in the brain and enriched in the postsynaptic density of glutamatergic synapses in the mammalian brain. Accumulating genetic evidence and functional data have implicated TNIK as a risk factor for psychiatric disorders. However, the endogenous substrates of TNIK in neurons are unknown. Here, we describe a novel selective small molecule inhibitor of the TNIK kinase family. Using this inhibitor, we report the identification of endogenous neuronal TNIK substrates by immunoprecipitation with a phosphomotif antibody followed by mass spectrometry. Phosphorylation consensus sequences were defined by phosphopeptide sequence analysis. Among the identified substrates were members of the delta-catenin family including p120-catenin, δ-catenin, and armadillo repeat gene deleted in velo-cardio-facial syndrome (ARVCF), each of which is linked to psychiatric or neurologic disorders. Using p120-catenin as a representative substrate, we show TNIK-induced p120-catenin phosphorylation in cells requires intact kinase activity and phosphorylation of TNIK at T181 and T187 in the activation loop. Addition of the small molecule TNIK inhibitor or knocking down TNIK by two shRNAs reduced endogenous p120-catenin phosphorylation in cells. Together, using a TNIK inhibitor and phosphomotif antibody, we identify endogenous substrates of TNIK in neurons, define consensus sequences for TNIK, and suggest signaling pathways by which TNIK influences synaptic development and function linked to psychiatric and neurologic disorders.

Discovery and preclinical profiling of 3-[4-(morpholin-4-yl)-7H-pyrrolo[2,3-d]pyrimidin-5-yl]benzonitrile (PF-06447475), a highly potent, selective, brain penetrant, and in vivo active LRRK2 kinase inhibitor.

Leucine rich repeat kinase 2 (LRRK2) has been genetically linked to Parkinson's disease (PD) by genome-wide association studies (GWAS). The most common LRRK2 mutation, G2019S, which is relatively rare in the total population, gives rise to increased kinase activity. As such, LRRK2 kinase inhibitors are potentially useful in the treatment of PD. We herein disclose the discovery and optimization of a novel series of potent LRRK2 inhibitors, focusing on improving kinome selectivity using a surrogate crystallography approach. This resulted in the identification of 14 (PF-06447475), a highly potent, brain penetrant and selective LRRK2 inhibitor which has been further profiled in in vivo safety and pharmacodynamic studies.

Kinase domain inhibition of leucine rich repeat kinase 2 (LRRK2) using a [1,2,4]triazolo[4,3-b]pyridazine scaffold.

Leucine rich repeat kinase 2 (LRRK2) has been genetically linked to Parkinson's disease (PD). The most common mutant, G2019S, increases kinase activity, thus LRRK2 kinase inhibitors are potentially useful in the treatment of PD. We herein disclose the structure, potential ligand-protein binding interactions, and pharmacological profiling of potent and highly selective kinase inhibitors based on a triazolopyridazine chemical scaffold.

Improving the odds of success in drug discovery: choosing the best compounds for in vivo toxicology studies.

A set of molecules that advanced into exploratory animal toxicology studies (two species) was examined to determine what properties contributed to success in these safety studies. Compounds were rigorously evaluated across numerous safety end points and classified as "pass" if a suitable in vivo therapeutic index (TI) was achieved for advancement into regulatory toxicology studies. The most predictive end point contributing to compound survival was a predicted human efficacious concentration (Ceff) of ≤250 nM (total drug) and ≤40 nM (free drug). This trend held across a wide range of CNS modes of action, encompassing targets such as enzymes, G-protein-coupled receptors, ion channels, and transporters.

Application of the bicyclo[1.1.1]pentane motif as a nonclassical phenyl ring bioisostere in the design of a potent and orally active γ-secretase inhibitor.

Replacement of the central, para-substituted fluorophenyl ring in the γ-secretase inhibitor 1 (BMS-708,163) with the bicyclo[1.1.1]pentane motif led to the discovery of compound 3, an equipotent enzyme inhibitor with significant improvements in passive permeability and aqueous solubility. The modified biopharmaceutical properties of 3 translated into excellent oral absorption characteristics (~4-fold ↑ C(max) and AUC values relative to 1) in a mouse model of γ-secretase inhibition. In addition, SAR studies into other fluorophenyl replacements indicate the intrinsic advantages of the bicyclo[1.1.1]pentane moiety over conventional phenyl ring replacements with respect to achieving an optimal balance of properties (e.g., γ-secretase inhibition, aqueous solubility/permeability, in vitro metabolic stability). Overall, this work enhances the scope of the [1.1.1]-bicycle beyond that of a mere "spacer" unit and presents a compelling case for its broader application as a phenyl group replacement in scenarios where the aromatic ring count impacts physicochemical parameters and overall drug-likeness.

U1A protein-stem loop 2 RNA recognition: prediction of structural differences from protein mutations.

Molecular dynamics (MD) simulations were carried out to compare the free and bound structures of wild type U1A protein with several Phe56 mutant U1A proteins that bind the target stem loop 2 (SL2) RNA with a range of affinities. The simulations indicate the free U1A protein is more flexible than the U1A-RNA complex for both wild type and Phe56 mutant systems. A complete analysis of the hydrogen-bonding (HB) and non-bonded (VDW) interactions over the course of the MD simulations suggested that changes in the interactions in the free U1A protein caused by the Phe56Ala and Phe56Leu mutations may stabilize the helical character in loop 3, and contribute to the weak binding of these proteins to SL2 RNA. Compared with wild type, changes in HB and VDW interactions in Phe56 mutants of the free U1A protein are global, and include differences in ß-sheet, loop 1 and loop 3 interactions. In the U1A-RNA complex, the Phe56Ala mutation leads to a series of differences in interactions that resonate through the complex, while the Phe56Leu and Phe56Trp mutations cause local differences around the site of mutation. The long-range networks of interactions identified in the simulations suggest that direct interactions and dynamic processes in both the free and bound forms contribute to complex stability.

Characterization of the dynamics of an essential helix in the U1A protein by time-resolved fluorescence measurements.

The RNA recognition motif (RRM), one of the most common RNA-binding domains, recognizes single-stranded RNA. A C-terminal helix that undergoes conformational changes upon binding is often an important contributor to RNA recognition. The N-terminal RRM of the U1A protein contains a C-terminal helix (helix C) that interacts with the RNA-binding surface of a beta-sheet in the free protein (closed conformation), but is directed away from this beta-sheet in the complex with RNA (open conformation). The dynamics of helix C in the free protein have been proposed to contribute to binding affinity and specificity. We report here a direct investigation of the dynamics of helix C in the free U1A protein on the nanosecond time scale using time-resolved fluorescence anisotropy. The results indicate that helix C is dynamic on a 2-3 ns time scale within a 20 degrees range of motion. Steady-state fluorescence experiments and molecular dynamics simulations suggest that the dynamical motion of helix C occurs within the closed conformation. Mutation of a residue on the beta-sheet that contacts helix C in the closed conformation dramatically destabilizes the complex (Phe56Ala) and alters the steady-state fluorescence, but not the time-resolved fluorescence anisotropy, of a Trp in helix C. Mutation of Asp90 in the hinge region between helix C and the remainder of the protein to Ala or Gly subtly alters the dynamics of the U1A protein and destabilizes the complex. Together these results show that helix C maintains a dynamic closed conformation that is stable to these targeted protein modifications and does not equilibrate with the open conformation on the nanosecond time scale.

Affinity and specificity of protein U1A-RNA complex formation based on an additive component free energy model.

An MM-GBSA computational protocol was used to investigate wild-type U1A-RNA and F56 U1A mutant experimental binding free energies. The trend in mutant binding free energies compared to wild-type is well-reproduced. Following application of a linear-response-like equation to scale the various energy components, the binding free energies agree quantitatively with observed experimental values. Conformational adaptation contributes to the binding free energy for both the protein and the RNA in these systems. Small differences in DeltaGs are the result of different and sometimes quite large relative contributions from various energetic components. Residual free energy decomposition indicates differences not only at the site of mutation, but throughout the entire protein. MM-GBSA and ab initio calculations performed on model systems suggest that stacking interactions may nearly, but not completely, account for observed differences in mutant binding affinities. This study indicates that there may be different underlying causes of ostensibly similar experimentally observed binding affinities of different mutants, and thus recommends caution in the interpretation of binding affinities and specificities purely by inspection.