Navigating the conformational landscape of G protein-coupled receptor kinases during allosteric activation

Fungal alkaloid malbrancheamide reorients the lipid binding domain of GRK5

G protein-coupled receptors (GPCRs) are the largest group of receptors involved in various types of signaling. GPCR signaling is regulated via receptor phosphorylation by G protein-coupled receptor kinases 5 (GRK5). Calmodulin (CaM), a universal Ca2+ sensor, inhibits receptor phosphorylation by binding to GRK5. However, the inhibitor malbrancheamide (MBC), which binds at CaM C-lobe, allows for receptor phosphorylation. To understand the phosphorylation mechanism by GRK5, we carried out a MD simulation of the CaM/GRK5 complex in the presence and absence of the MBC inhibitor. The lipid binding domain (LBD) of GRK5 adopted different positions in the presence and absence of inhibitor. Furthermore, the inhibitor MBC restricted the movement of the N-lobe tether (NLT) loop, probably blocking the autophosphorylation of GRK5.Communicated by Ramaswamy H. Sarma.

Multilevel superposition for deciphering the conformational variability of protein ensembles

The dynamics and variability of protein conformations are directly linked to their functions. Many comparative studies of X-ray protein structures have been conducted to elucidate the relevant conformational changes, dynamics and heterogeneity. The rapid increase in the number of experimentally determined structures has made comparison an effective tool for investigating protein structures. For example, it is now possible to compare structural ensembles formed by enzyme species, variants or the type of ligands bound to them. In this study, the author developed a multilevel model for estimating two covariance matrices that represent inter- and intra-ensemble variability in the Cartesian coordinate space. Principal component analysis using the two estimated covariance matrices identified the inter-/intra-enzyme variabilities, which seemed to be important for the enzyme functions, with the illustrative examples of cytochrome P450 family 2 enzymes and class A $\beta$-lactamases. In P450, in which each enzyme has its own active site of a distinct size, an active-site motion shared universally between the enzymes was captured as the first principal mode of the intra-enzyme covariance matrix. In this case, the method was useful for understanding the conformational variability after adjusting for the differences between enzyme sizes. The developed method is advantageous in small ensemble-size problems and hence promising for use in comparative studies on experimentally determined structures where ensemble sizes are smaller than those generated, for example, by molecular dynamics simulations.

CMV-encoded GPCRs in infection, disease, and pathogenesis

G protein coupled receptors (GPCRs) are seven-transmembrane domain proteins that modulate cellular processes in response to external stimuli. These receptors represent the largest family of membrane proteins, and in mammals, their signaling regulates important physiological functions, such as vision, taste, and olfaction. Many organisms, including yeast, slime molds, and viruses encode GPCRs. Cytomegaloviruses (CMVs) are large, betaherpesviruses, that encode viral GPCRs (vGPCRs). Human CMV (HCMV) encodes four vGPCRs, including UL33, UL78, US27, and US28. Each of these vGPCRs, as well as their rodent and primate orthologues, have been investigated for their contributions to viral infection and disease. Herein, we discuss how the CMV vGPCRs function during lytic and latent infection, as well as our understanding of how they impact viral pathogenesis.

G protein-coupled receptor interactions with arrestins and GPCR kinases: the unresolved issue of signal bias

G protein–coupled receptor (GPCR) kinases (GRKs) and arrestins interact with agonist-bound GPCRs to promote receptor desensitization and downregulation. They also trigger signaling cascades distinct from those of heterotrimeric G proteins. Biased agonists for GPCRs that favor either heterotrimeric G protein or GRK/arrestin signaling are of profound pharmacological interest because they could usher in a new generation of drugs with greatly reduced side effects. One mechanism by which biased agonism might occur is by stabilizing receptor conformations that preferentially bind to GRKs and/or arrestins. In this review, we explore this idea by comparing structures of GPCRs bound to heterotrimeric G proteins with those of the same GPCRs in complex with arrestins and GRKs. The arrestin and GRK complexes all exhibit high conformational heterogeneity, which is likely a consequence of their unusual ability to adapt and bind to hundreds of different GPCRs. This dynamic behavior, along with the experimental tactics required to stabilize GPCR complexes for biophysical analysis, confounds these comparisons, but some possible molecular mechanisms of bias are beginning to emerge. We also examine if and how the recent structures advance our understanding of how arrestins parse the “phosphorylation barcodes” installed in the intracellular loops and tails of GPCRs by GRKs. In the future, structural analyses of arrestins in complex with intact receptors that have well-defined native phosphorylation barcodes, such as those installed by the two nonvisual subfamilies of GRKs, will be particularly illuminating.

Structures of rhodopsin in complex with G-protein-coupled receptor kinase 1

G-protein-coupled receptor (GPCR) kinases (GRKs) selectively phosphorylate activated GPCRs, thereby priming them for desensitization¹. Although it is unclear how GRKs recognize these receptors^2-4, a conserved region at the GRK N terminus is essential for this process^5-8. Here we report a series of cryo-electron microscopy single-particle reconstructions of light-activated rhodopsin (Rho*) bound to rhodopsin kinase (GRK1), wherein the N terminus of GRK1 forms a helix that docks into the open cytoplasmic cleft of Rho*. The helix also packs against the GRK1 kinase domain and stabilizes it in an active configuration. The complex is further stabilized by electrostatic interactions between basic residues that are conserved in most GPCRs and acidic residues that are conserved in GRKs. We did not observe any density for the regulator of G-protein signalling homology domain of GRK1 or the C terminus of rhodopsin. Crosslinking with mass spectrometry analysis confirmed these results and revealed dynamic behaviour in receptor-bound GRK1 that would allow the phosphorylation of multiple sites in the receptor tail. We have identified GRK1 residues whose mutation augments kinase activity and crosslinking with Rho*, as well as residues that are involved in activation by acidic phospholipids. From these data, we present a general model for how a small family of protein kinases can recognize and be activated by hundreds of different GPCRs.

The Open Question of How GPCRs Interact with GPCR Kinases (GRKs)

G protein-coupled receptors (GPCRs), which regulate a vast number of eukaryotic processes, are desensitized by various mechanisms but, most importantly, by the GPCR kinases (GRKs). Ever since GRKs were first identified, investigators have sought to determine which structural features of GRKs are used to select for the agonist-bound states of GPCRs and how this binding event in turn enhances GRK catalytic activity. Despite a wealth of molecular information from high-resolution crystal structures of GRKs, the mechanisms driving activation have remained elusive, in part because the GRK N-terminus and active site tether region, previously proposed to serve as a receptor docking site and to be key to kinase domain closure, are often disordered or adopt inconsistent conformations. However, two recent studies have implicated other regions of GRKs as being involved in direct interactions with active GPCRs. Atomic resolution structures of GPCR–GRK complexes would help refine these models but are, so far, lacking. Here, we assess three distinct models for how GRKs recognize activated GPCRs, discuss limitations in the approaches used to generate them, and then experimentally test a hypothetical GPCR interaction site in GRK2 suggested by the two newest models.

Targeting G protein-coupled receptor kinases (GRKs) to G protein-coupled receptors

G protein-coupled receptors (GPCRs) interact with three protein families following agonist binding: heterotrimeric G proteins, G protein-coupled receptor kinases (GRKs) and arrestins. GRK-mediated phosphorylation of GPCRs promotes arrestin binding to uncouple the receptor from G protein, a process called desensitization, and for many GPCRs, arrestin binding also promotes receptor endocytosis and intracellular signaling. Thus, GRKs play a central role in modulating GPCR signaling and localization. Here we review recent advances in this field which include additional insight into how GRKs target GPCRs and bias signaling, and the development of specific inhibitors to dissect GRK function in model systems.

Stabilization of the Hinge Region of Human E-selectin Enhances Binding Affinity to Ligands Under Force

Introduction

E-selectin is a member of the selectin family of cell adhesion molecules expressed on the plasma membrane of inflamed endothelium and facilitates initial leukocyte tethering and subsequent cell rolling during the early stages of the inflammatory response via binding to glycoproteins expressing sialyl Lewis^X and sialyl Lewis^A (sLe^X/A). Existing crystal structures of the extracellular lectin/EGF-like domain of E-selectin complexed with sLe^X have revealed that E-selectin can exist in two conformation states, a low affinity (bent) conformation, and a high affinity (extended) conformation. The differentiating characteristic of the two conformations is the interdomain angle between the lectin and the EGF-like domain.

Methods

Using molecular dynamics (MD) simulations we observed that in the absence of tensile force E-selectin undergoes spontaneous switching between the two conformational states at equilibrium. A single amino acid substitution at residue 2 (serine to tyrosine) on the lectin domain favors the extended conformation.

Results

Steered molecular dynamics (SMD) simulations of E-selectin and PSGL-1 in conjunction with experimental cell adhesion assays show a longer binding lifetime of E-selectin (S2Y) to PSGL-1 compared to wildtype protein.

Conclusions

The findings in this study advance our understanding into how the structural makeup of E-selectin allosterically influences its adhesive dynamics.

Structure of a GRK5-Calmodulin Complex Reveals Molecular Mechanism of GRK Activation and Substrate Targeting

The phosphorylation of G protein-coupled receptors (GPCRs) by GPCR kinases (GRKs) facilitates arrestin binding and receptor desensitization. Although this process can be regulated by Ca²⁺-binding proteins such as calmodulin (CaM) and recoverin, the molecular mechanisms are poorly understood. Here, we report structural, computational, and biochemical analysis of a CaM complex with GRK5, revealing how CaM shapes GRK5 response to calcium. The CaM N and C domains bind independently to two helical regions at the GRK5 N and C termini to inhibit GPCR phosphorylation, though only the C domain interaction disrupts GRK5 membrane association, thereby facilitating cytoplasmic translocation. The CaM N domain strongly activates GRK5 via ordering of the amphipathic αN-helix of GRK5 and allosteric disruption of kinase-RH domain interaction for phosphorylation of cytoplasmic GRK5 substrates. These results provide a framework for understanding how two functional effects, GRK5 activation and localization, can cooperate under control of CaM for selective substrate targeting by GRK5.

Detecting Functional Dynamics in Proteins with Comparative Perturbed-Ensembles Analysis

Recent advances have made all-atom molecular dynamics (MD) a powerful tool to sample the conformational energy landscape. There are still however three major challenges in the application of MD to biological systems: accuracy of force field, time scale, and the analysis of simulation trajectories. Significant progress in addressing the first two challenges has been made and extensively reviewed previously. This Account focuses on strategies of analyzing simulation data of biomolecules that also covers ways to properly design simulations and validate simulation results. In particular, we examine an approach named comparative perturbed-ensembles analysis, which we developed to efficiently detect dynamics in protein MD simulations that can be linked to biological functions. In our recent studies, we implemented this approach to understand allosteric regulations in several disease-associated human proteins. The central task of a comparative perturbed-ensembles analysis is to compare two or more conformational ensembles of a system generated by MD simulations under distinct perturbation conditions. Perturbations can be different sequence variations, ligand-binding conditions, and other physical/chemical modifications of the system. Each simulation is long enough (e.g., microsecond-long) to ensure sufficient sampling of the local substate. Then, sophisticated bioinformatic and statistical tools are applied to extract function-related information from the simulation data, including principal component analysis, residue-residue contact analysis, difference contact network analysis (dCNA) based on the graph theory, and statistical analysis of side-chain conformations. Computational findings are further validated with experimental data. By comparing distinct conformational ensembles, functional micro- to millisecond dynamics can be inferred. In contrast, such a time scale is difficult to reach in a single simulation; even when reached for a single condition of a system, it is elusive as to what dynamical motions are related to functions without, for example, comparing free and substrate-bound proteins at the minimum. We illustrate our approach with three examples. First, we discuss using the approach to identify allosteric pathways in cyclophilin A (CypA), a member of a ubiquitous class of peptidyl-prolyl cis-trans isomerase enzymes. By comparing side-chain torsion-angle distributions of CypA in wild-type and mutant forms, we identified three pathways: two are consistent with recent nuclear magnetic resonance experiments, whereas the third is a novel pathway. Second, we show how the approach enables a dynamical-evolution analysis of the human cyclophilin family. In the analysis, both conserved and divergent conformational dynamics across three cyclophilin isoforms (CypA, CypD, and CypE) were summarized. The conserved dynamics led to the discovery of allosteric networks resembling those found in CypA. A residue wise determinant underlying the unique dynamics in CypD was also detected and validated with additional mutational MD simulations. In the third example, we applied the approach to elucidate a peptide sequence-dependent allosteric mechanism in human Pin 1, a phosphorylation-dependent peptidyl-prolyl isomerase. We finally present our outlook of future directions. Especially, we envisage how the approach could help open a new avenue in drug discovery.

Structure-Based Design of Selective, Covalent G Protein-Coupled Receptor Kinase 5 Inhibitors

The ability of G protein-coupled receptor (GPCR) kinases (GRKs) to regulate desensitization of GPCRs has made GRK2 and GRK5 attractive targets for treating heart failure and other diseases such as cancer. Although advances have been made toward developing inhibitors that are selective for GRK2, there have been far fewer reports of GRK5 selective compounds. Herein, we describe the development of GRK5 subfamily selective inhibitors, 5 and 16d that covalently interact with a nonconserved cysteine (Cys474) unique to this subfamily. Compounds 5 and 16d feature a highly amenable pyrrolopyrimidine scaffold that affords high nanomolar to low micromolar activity that can be easily modified with Michael acceptors with various reactivities and geometries. Our work thereby establishes a new pathway toward further development of subfamily selective GRK inhibitors and establishes Cys474 as a new and useful covalent handle in GRK5 drug discovery.

Characterization of a hyperphosphorylated variant of G protein-coupled receptor kinase 5 expressed in E. coli

G protein-coupled receptors (GPCRs) are the largest family of cell-surface receptors in humans and regulate numerous physiological processes through the activation of heterotrimeric G proteins. GPCR kinases (GRKs) selectively phosphorylate active GPCRs, which promotes arrestin binding, receptor internalization, and initiation of alternative signaling pathways. GRK5 is a representative member of one of three GRK subfamilies that does not need post-translational lipidation or other binding partners to exhibit full activity against GPCRs, rendering it a useful tool for biophysical studies directed at characterizing GRK function. However, recombinant expression of GRK5 has thus far been limited to insect and mammalian systems. Here, we describe the expression of functional GRK5 in E. coli and its purification and biochemical characterization. Bacterially expressed GRK5 is hyperphosphorylated, primarily in regions known to be flexible from prior crystal structures, which slightly decreases its catalytic activity toward receptor substrates. Mutation of a single phosphorylation site, Thr10, restores kinetic parameters to those of GRK5 purified from insect cells. Consequently, bacterial expression will allow for production of GRK5 at a reduced cost and faster pace and would facilitate production of isotopically labeled kinase for NMR studies or for the incorporation of unnatural amino acids.

Computational study of paroxetine-like inhibitors reveals new molecular insight to inhibit GRK2 with selectivity over ROCK1

The G-protein coupled receptor kinase 2 (GRK2) regulates the desensitization of beta-adrenergic receptors (β-AR), and its overexpression has been implicated in heart failure. Hence, the inhibition of GRK2 is considered to be an important drug target for the treatment of heart failure. Due to the high sequence similarity of GRK2 with the A, G, and C family (AGC family) of kinases, the inhibition of GRK2 also leads to the inhibition of AGC kinases such as Rho-associated coiled-coil kinase 1 (ROCK1). Therefore, unraveling the mechanisms to selectively inhibit GRK2 poses an important challenge. We have performed molecular docking, three dimensional quantitative structure activity relationship (3D-QSAR), molecular dynamics (MD) simulation, and free energy calculations techniques on a series of 53 paroxetine-like compounds to understand the structural properties desirable for enhancing the inhibitory activity for GRK2 with selectivity over ROCK1. The formation of stable hydrogen bond interactions with the residues Phe202 and Lys220 of GRK2 seems to be important for selective inhibition of GRK2. Electropositive substituents at the piperidine ring and electronegative substituents near the amide linker between the benzene ring and pyrazole ring showed a higher inhibitory preference for GRK2 over ROCK1. This study may be used in designing more potent and selective GRK2 inhibitors for therapeutic intervention of heart failure.

Structural Basis for G Protein-Coupled Receptor Signaling

G protein-coupled receptors (GPCRs), which mediate processes as diverse as olfaction and maintenance of metabolic homeostasis, have become the single most effective class of therapeutic drug targets. As a result, understanding the molecular basis for their activity is of paramount importance. Recent technological advances have made GPCR structural biology increasingly tractable, offering views of these receptors in unprecedented atomic detail. Structural and biophysical data have shown that GPCRs function as complex allosteric machines, communicating ligand-binding events through conformational change. Changes in receptor conformation lead to activation of effector proteins, such as G proteins and arrestins, which are themselves conformational switches. Here, we review how structural biology has illuminated the agonist-induced cascade of conformational changes that culminate in a cellular response to GPCR activation.

The dopamine D2 receptor can directly recruit and activate GRK2 without G protein activation

The dopamine D2 receptor (D2R) is a G protein-coupled receptor (GPCR) that is critical for many central nervous system functions. The D2R carries out these functions by signaling through two transducers: G proteins and β-arrestins (βarrs). Selectively engaging either the G protein or βarr pathway may be a way to improve drugs targeting GPCRs. The current model of GPCR signal transduction posits a chain of events where G protein activation ultimately leads to βarr recruitment. GPCR kinases (GRKs), which are regulated by G proteins and whose kinase action facilitates βarr recruitment, bridge these pathways. Therefore, βarr recruitment appears to be intimately tied to G protein activation via GRKs. Here we sought to understand how GRK2 action at the D2R would be disrupted when G protein activation is eliminated and the effect of this on βarr recruitment. We used two recently developed biased D2R mutants that can preferentially interact either with G proteins or βarrs as well as a βarr-biased D2R ligand, UNC9994. With these functionally selective tools, we investigated the mechanism whereby the βarr-preferring D2R achieves βarr pathway activation in the complete absence of G protein activation. We describe how direct, G protein-independent recruitment of GRK2 drives interactions at the βarr-preferring D2R and also contributes to βarr recruitment at the WT D2R. Additionally, we found an additive interaction between the βarr-preferring D2R mutant and UNC9994. These results reveal that the D2R can directly recruit GRK2 without G protein activation and that this mechanism may have relevance to achieving βarr-biased signaling.

Structural Determinants Influencing the Potency and Selectivity of Indazole-Paroxetine Hybrid G Protein-Coupled Receptor Kinase 2 Inhibitors

G protein-coupled receptor kinases (GRKs) phosphorylate activated receptors to promote arrestin binding, decoupling from heterotrimeric G proteins, and internalization. GRK2 and GRK5 are overexpressed in the failing heart and thus have become therapeutic targets. Previously, we discovered two classes of GRK2-selective inhibitors, one stemming from GSK180736A, a Rho-associated coiled-coil containing kinase 1 (ROCK1) inhibitor, the other from paroxetine, a selective serotonin-reuptake inhibitor. These two classes of compounds bind to the GRK2 active site in a similar configuration but contain different hinge-binding "warheads": indazole and benzodioxole, respectively. We surmised from our prior studies that an indazole would be the stronger hinge binder and would impart increased potency when substituted for benzodioxole in paroxetine derivatives. To test this hypothesis, we synthesized a series of hybrid compounds that allowed us to compare the effects of inhibitors that differ only in the identity of the warhead. The indazole-paroxetine analogs were indeed more potent than their respective benzodioxole derivatives but lost selectivity. To investigate how these two warheads dictate selectivity, we determined the crystal structures of three of the indazole hybrid compounds (CCG224061, CCG257284, and CCG258748) in complex with GRK2-Gβγ Comparison of these structures with those of analogous benzodioxole-containing complexes confirmed that the indazole-paroxetine hybrids form stronger interactions with the hinge of the kinase but also stabilize a distinct conformation of the kinase domain of GRK2 compared with previous complexes with paroxetine analogs. This conformation is analogous to one that can be assumed by GRK5, at least partially explaining the loss in selectivity.