Structural mechanism of a drug-binding process involving a large conformational change of the protein target

Proteins often undergo large conformational changes when binding small molecules, but atomic-level descriptions of such events have been elusive. Here, we report unguided molecular dynamics simulations of Abl kinase binding to the cancer drug imatinib. In the simulations, imatinib first selectively engages Abl kinase in its autoinhibitory conformation. Consistent with inferences drawn from previous experimental studies, imatinib then induces a large conformational change of the protein to reach a bound complex that closely resembles published crystal structures. Moreover, the simulations reveal a surprising local structural instability in the C-terminal lobe of Abl kinase during binding. The unstable region includes a number of residues that, when mutated, confer imatinib resistance by an unknown mechanism. Based on the simulations, NMR spectra, hydrogen-deuterium exchange measurements, and thermostability measurements and estimates, we suggest that these mutations confer imatinib resistance by exacerbating structural instability in the C-terminal lobe, rendering the imatinib-bound state energetically unfavorable.

Molecular mechanisms of fentanyl mediated β-arrestin biased signaling

The development of novel analgesics with improved safety profiles to combat the opioid epidemic represents a central question to G protein coupled receptor structural biology and pharmacology: What chemical features dictate G protein or β-arrestin signaling? Here we use adaptively biased molecular dynamics simulations to determine how fentanyl, a potent β-arrestin biased agonist, binds the μ-opioid receptor (μOR). The resulting fentanyl-bound pose provides rational insight into a wealth of historical structure-activity-relationship on its chemical scaffold. Following an in-silico derived hypothesis we found that fentanyl and the synthetic opioid peptide DAMGO require M153 to induce β-arrestin coupling, while M153 was dispensable for G protein coupling. We propose and validate an activation mechanism where the n-aniline ring of fentanyl mediates μOR β-arrestin through a novel M153 “microswitch” by synthesizing fentanyl-based derivatives that exhibit complete, clinically desirable, G protein biased coupling. Together, these results provide molecular insight into fentanyl mediated β-arrestin biased signaling and a rational framework for further optimization of fentanyl-based analgesics with improved safety profiles.

The global opioid crisis has drawn significant attention to the risks associated with over-use of synthetic opioids. Despite the public attention, and perhaps in-line with the profit-based incentives of the pharmaceutical industry, there is no public structure of mu-opioid receptor bound to fentanyl or fentanyl derivatives. A publicly available structure of the complex would allow open-source development of safer painkillers and synthetic antagonists. Current overdose antidotes, antagonists, require natural products in their synthesis which persists a sizable barrier to market and develop better antidotes. In this work we use advance molecular dynamics techniques to obtain the bound geometry of mu-opioid receptor with fentanyl (and derivatives) and uncovered a novel molecular switch involved in receptor activation. Based on our in-silico structure, we synthesized and tested novel compounds to validate our predicted structure. Herein we report the bound state of several dangerous fentanyl derivatives and introduce new derivatives with signaling profiles that may lead to lower risk of respiratory depression.

Structure of formylpeptide receptor 2-Gi complex reveals insights into ligand recognition and signaling

Formylpeptide receptors (FPRs) as G protein-coupled receptors (GPCRs) can recognize formylpeptides derived from pathogens or host cells to function in host defense and cell clearance. In addition, FPRs, especially FPR2, can also recognize other ligands with a large chemical diversity generated at different stages of inflammation to either promote or resolve inflammation in order to maintain a balanced inflammatory response. The mechanism underlying promiscuous ligand recognition and activation of FPRs is not clear. Here we report a cryo-EM structure of FPR2-Gi signaling complex with a peptide agonist. The structure reveals a widely open extracellular region with an amphiphilic environment for ligand binding. Together with computational docking and simulation, the structure suggests a molecular basis for the recognition of formylpeptides and a potential mechanism of receptor activation, and reveals conserved and divergent features in Gi coupling. Our results provide a basis for understanding the molecular mechanism of the functional promiscuity of FPRs.

Molecular Basis for Hormone Recognition and Activation of Corticotropin-Releasing Factor Receptors

Corticotropin-releasing factor (CRF) and the three related peptides urocortins 1-3 (UCN1-UCN3) are endocrine hormones that control the stress responses by activating CRF1R and CRF2R, two members of class B G-protein-coupled receptors (GPCRs). Here, we present two cryoelectron microscopy (cryo-EM) structures of UCN1-bound CRF1R and CRF2R with the stimulatory G protein. In both structures, UCN1 adopts a single straight helix with its N terminus dipped into the receptor transmembrane bundle. Although the peptide-binding residues in CRF1R and CRF2R are different from other members of class B GPCRs, the residues involved in receptor activation and G protein coupling are conserved. In addition, both structures reveal bound cholesterol molecules to the receptor transmembrane helices. Our structures define the basis of ligand-binding specificity in the CRF receptor-hormone system, establish a common mechanism of class B GPCR activation and G protein coupling, and provide a paradigm for studying membrane protein-lipid interactions for class B GPCRs.

A complex structure of arrestin-2 bound to a G protein-coupled receptor

Arrestins comprise a family of signal regulators of G-protein-coupled receptors (GPCRs), which include arrestins 1 to 4. While arrestins 1 and 4 are visual arrestins dedicated to rhodopsin, arrestins 2 and 3 (Arr2 and Arr3) are β-arrestins known to regulate many nonvisual GPCRs. The dynamic and promiscuous coupling of Arr2 to nonvisual GPCRs has posed technical challenges to tackle the basis of arrestin binding to GPCRs. Here we report the structure of Arr2 in complex with neurotensin receptor 1 (NTSR1), which reveals an overall assembly that is strikingly different from the visual arrestin-rhodopsin complex by a 90° rotation of Arr2 relative to the receptor. In this new configuration, intracellular loop 3 (ICL3) and transmembrane helix 6 (TM6) of the receptor are oriented toward the N-terminal domain of the arrestin, making it possible for GPCRs that lack the C-terminal tail to couple Arr2 through their ICL3. Molecular dynamics simulation and crosslinking data further support the assembly of the Arr2‒NTSR1 complex. Sequence analysis and homology modeling suggest that the Arr2‒NTSR1 complex structure may provide an alternative template for modeling arrestin-GPCR interactions.

Structure and dynamics of the active human parathyroid hormone receptor-1

The parathyroid hormone receptor-1 (PTH1R) is a class B G protein-coupled receptor central to calcium homeostasis and a therapeutic target for osteoporosis and hypoparathyroidism. Here we report the cryo-electron microscopy structure of human PTH1R bound to a long-acting PTH analog and the stimulatory G protein. The bound peptide adopts an extended helix with its amino terminus inserted deeply into the receptor transmembrane domain (TMD), which leads to partial unwinding of the carboxyl terminus of transmembrane helix 6 and induces a sharp kink at the middle of this helix to allow the receptor to couple with G protein. In contrast to a single TMD structure state, the extracellular domain adopts multiple conformations. These results provide insights into the structural basis and dynamics of PTH binding and receptor activation.

Conformational heterogeneity of the allosteric drug and metabolite (ADaM) site in AMP-activated protein kinase (AMPK)

AMP-activated protein kinase (AMPK) is a master regulator of energy homeostasis and a promising drug target for managing metabolic diseases such as type 2 diabetes. Many pharmacological AMPK activators, and possibly unidentified physiological metabolites, bind to the allosteric drug and metabolite (ADaM) site at the interface between the kinase domain (KD) in the α-subunit and the carbohydrate-binding module (CBM) in the β-subunit. Here, using double electron-electron resonance (DEER) spectroscopy, we demonstrate that the CBM-KD interaction is partially dissociated and the interface highly disordered in the absence of pharmacological ADaM site activators as inferred from a low depth of modulation and broad DEER distance distributions. ADaM site ligands such as 991, and to a lesser degree phosphorylation, stabilize the KD-CBM association and strikingly reduce conformational heterogeneity in the ADaM site. Our findings that the ADaM site, formed by the KD-CBM interaction, can be modulated by diverse ligands and by phosphorylation suggest that it may function as a hub for integrating regulatory signals.

Crystal structure of the human 5-HT1B serotonin receptor bound to an inverse agonist

5-hydroxytryptamine (5-HT, also known as serotonin) regulates many physiological processes through the 5-HT receptor family. Here we report the crystal structure of 5-HT1B subtype receptor (5-HT1BR) bound to the psychotropic serotonin receptor inverse agonist methiothepin (MT). Crystallization was facilitated by replacing ICL3 with a novel optimized variant of BRIL (OB1) that enhances the formation of intermolecular polar interactions, making OB1 a potential useful tool for structural studies of membrane proteins. Unlike the agonist ergotamine (ERG), MT occupies only the conserved orthosteric binding pocket, explaining the wide spectrum effect of MT on serotonin receptors. Compared with ERG, MT shifts toward TM6 and sterically pushes residues W3276.48, F3306.50 and F3316.51 from inside the orthosteric binding pocket, leading to an outward movement of the extracellular end and a corresponding inward shift of the intracellular end of TM6, a feature shared by other reported inactive G protein-coupled receptor (GPCR) structures. Together with the previous agonist-bound serotonin receptor structures, the inverse agonist-bound 5-HT1BR structure identifies a basis for the ligand-mediated switch of 5-HT1BR activity and provides a structural understanding of the inactivation mechanism of 5-HT1BR and some other class A GPCRs, characterized by ligand-induced outward movement of the extracellular end of TM6 that is coupled with inward movement of the cytoplasmic end of this helix.

Identification of Phosphorylation Codes for Arrestin Recruitment by G Protein-Coupled Receptors

G protein-coupled receptors (GPCRs) mediate diverse signaling in part through interaction with arrestins, whose binding promotes receptor internalization and signaling through G protein-independent pathways. High-affinity arrestin binding requires receptor phosphorylation, often at the receptor's C-terminal tail. Here, we report an X-ray free electron laser (XFEL) crystal structure of the rhodopsin-arrestin complex, in which the phosphorylated C terminus of rhodopsin forms an extended intermolecular β sheet with the N-terminal β strands of arrestin. Phosphorylation was detected at rhodopsin C-terminal tail residues T336 and S338. These two phospho-residues, together with E341, form an extensive network of electrostatic interactions with three positively charged pockets in arrestin in a mode that resembles binding of the phosphorylated vasopressin-2 receptor tail to β-arrestin-1. Based on these observations, we derived and validated a set of phosphorylation codes that serve as a common mechanism for phosphorylation-dependent recruitment of arrestins by GPCRs.

Deconvoluting AMP-activated protein kinase (AMPK) adenine nucleotide binding and sensing

AMP-activated protein kinase (AMPK) is a central cellular energy sensor that adapts metabolism and growth to the energy state of the cell. AMPK senses the ratio of adenine nucleotides (adenylate energy charge) by competitive binding of AMP, ADP, and ATP to three sites (CBS1, CBS3, and CBS4) in its γ-subunit. Because these three binding sites are functionally interconnected, it remains unclear how nucleotides bind to individual sites, which nucleotides occupy each site under physiological conditions, and how binding to one site affects binding to the other sites. Here, we comprehensively analyze nucleotide binding to wild-type and mutant AMPK protein complexes by quantitative competition assays and by hydrogen-deuterium exchange MS. We also demonstrate that NADPH, in addition to the known AMPK ligand NADH, directly and competitively binds AMPK at the AMP-sensing CBS3 site. Our findings reveal how AMP binding to one site affects the conformation and adenine nucleotide binding at the other two sites and establish CBS3, and not CBS1, as the high affinity exchangeable AMP/ADP/ATP-binding site. We further show that AMP binding at CBS4 increases AMP binding at CBS3 by 2 orders of magnitude and reverses the AMP/ATP preference of CBS3. Together, these results illustrate how the three CBS sites collaborate to enable highly sensitive detection of cellular energy states to maintain the tight ATP homeostastis required for cellular metabolism.

A fast, open source implementation of adaptive biasing potentials uncovers a ligand design strategy for the chromatin regulator BRD4

In this communication we introduce an efficient implementation of adaptive biasing that greatly improves the speed of free energy computation in molecular dynamics simulations. We investigated the use of accelerated simulations to inform on compound design using a recently reported and clinically relevant inhibitor of the chromatin regulator BRD4 (bromodomain-containing protein 4). Benchmarking on our local compute cluster, our implementation achieves up to 2.5 times more force calls per day than plumed2. Results of five 1 μs-long simulations are presented, which reveal a conformational switch in the BRD4 inhibitor between a binding competent and incompetent state. Stabilization of the switch led to a -3 kcal/mol improvement of absolute binding free energy. These studies suggest an unexplored ligand design principle and offer new actionable hypotheses for medicinal chemistry efforts against this druggable epigenetic target class.

Rearrangement of a polar core provides a conserved mechanism for constitutive activation of class B G protein-coupled receptors

The glucagon receptor (GCGR) belongs to the secretin-like (class B) family of G protein-coupled receptors (GPCRs) and is activated by the peptide hormone glucagon. The structures of an activated class B GPCR have remained unsolved, preventing a mechanistic understanding of how these receptors are activated. Using a combination of structural modeling and mutagenesis studies, we present here two modes of ligand-independent activation of GCGR. First, we identified a GCGR-specific hydrophobic lock comprising Met-338 and Phe-345 within the IC3 loop and transmembrane helix 6 (TM6) and found that this lock stabilizes the TM6 helix in the inactive conformation. Disruption of this hydrophobic lock led to constitutive G protein and arrestin signaling. Second, we discovered a polar core comprising conserved residues in TM2, TM3, TM6, and TM7, and mutations that disrupt this polar core led to constitutive GCGR activity. On the basis of these results, we propose a mechanistic model of GCGR activation in which TM6 is held in an inactive conformation by the conserved polar core and the hydrophobic lock. Mutations that disrupt these inhibitory elements allow TM6 to swing outward to adopt an active TM6 conformation similar to that of the canonical β₂-adrenergic receptor complexed with G protein and to that of rhodopsin complexed with arrestin. Importantly, mutations in the corresponding polar core of several other members of class B GPCRs, including PTH1R, PAC1R, VIP1R, and CRFR1, also induce constitutive G protein signaling, suggesting that the rearrangement of the polar core is a conserved mechanism for class B GPCR activation.

Tumor Targeting with Novel 6-Substituted Pyrrolo [2,3-d] Pyrimidine Antifolates with Heteroatom Bridge Substitutions via Cellular Uptake by Folate Receptor α and the Proton-Coupled Folate Transporter and Inhibition of de Novo Purine Nucleotide Biosynthesis

Targeted antifolates with heteroatom replacements of the carbon vicinal to the phenyl ring in 1 by N (4), O (8), or S (9), or with N-substituted formyl (5), acetyl (6), or trifluoroacetyl (7) moieties, were synthesized and tested for selective cellular uptake by folate receptor (FR) α and β or the proton-coupled folate transporter. Results show increased in vitro antiproliferative activity toward engineered Chinese hamster ovary cells expressing FRs by 4-9 over the CH2 analogue 1. Compounds 4-9 inhibited de novo purine biosynthesis and glycinamide ribonucleotide formyltransferase (GARFTase). X-ray crystal structures for 4 with FRα and GARFTase showed that the bound conformations of 4 required flexibility for attachment to both FRα and GARFTase. In mice bearing IGROV1 ovarian tumor xenografts, 4 was highly efficacious. Our results establish that heteroatom substitutions in the 3-atom bridge region of 6-substituted pyrrolo[2,3-d]pyrimidines related to 1 provide targeted antifolates that warrant further evaluation as anticancer agents.

X-ray laser diffraction for structure determination of the rhodopsin-arrestin complex

Serial femtosecond X-ray crystallography (SFX) using an X-ray free electron laser (XFEL) is a recent advancement in structural biology for solving crystal structures of challenging membrane proteins, including G-protein coupled receptors (GPCRs), which often only produce microcrystals. An XFEL delivers highly intense X-ray pulses of femtosecond duration short enough to enable the collection of single diffraction images before significant radiation damage to crystals sets in. Here we report the deposition of the XFEL data and provide further details on crystallization, XFEL data collection and analysis, structure determination, and the validation of the structural model. The rhodopsin-arrestin crystal structure solved with SFX represents the first near-atomic resolution structure of a GPCR-arrestin complex, provides structural insights into understanding of arrestin-mediated GPCR signaling, and demonstrates the great potential of this SFX-XFEL technology for accelerating crystal structure determination of challenging proteins and protein complexes.

Structural basis of AMPK regulation by adenine nucleotides and glycogen

AMP-activated protein kinase (AMPK) is a central cellular energy sensor and regulator of energy homeostasis, and a promising drug target for the treatment of diabetes, obesity, and cancer. Here we present low-resolution crystal structures of the human α1β2γ1 holo-AMPK complex bound to its allosteric modulators AMP and the glycogen-mimic cyclodextrin, both in the phosphorylated (4.05 Å) and non-phosphorylated (4.60 Å) state. In addition, we have solved a 2.95 Å structure of the human kinase domain (KD) bound to the adjacent autoinhibitory domain (AID) and have performed extensive biochemical and mutational studies. Together, these studies illustrate an underlying mechanism of allosteric AMPK modulation by AMP and glycogen, whose binding changes the equilibria between alternate AID (AMP) and carbohydrate-binding module (glycogen) interactions.

Molecular dynamics of CYP2D6 polymorphisms in the absence and presence of a mechanism-based inactivator reveals changes in local flexibility and dominant substrate access channels

Cytochrome P450 enzymes (CYPs) represent an important enzyme superfamily involved in metabolism of many endogenous and exogenous small molecules. CYP2D6 is responsible for ∼ 15% of CYP-mediated drug metabolism and exhibits large phenotypic diversity within CYPs with over 100 different allelic variants. Many of these variants lead to functional changes in enzyme activity and substrate selectivity. Herein, a molecular dynamics comparative analysis of four different variants of CYP2D6 was performed. The comparative analysis included simulations with and without SCH 66712, a ligand that is also a mechanism-based inactivator, in order to investigate the possible structural basis of CYP2D6 inactivation. Analysis of protein stability highlighted significantly altered flexibility in both proximal and distal residues from the variant residues. In the absence of SCH 66712, *34, *17-2, and *17-3 displayed more flexibility than *1, and *53 displayed more rigidity. SCH 66712 binding reversed flexibility in *17-2 and *17-3, through *53 remained largely rigid. Throughout simulations with docked SCH 66712, ligand orientation within the heme-binding pocket was consistent with previously identified sites of metabolism and measured binding energies. Subsequent tunnel analysis of substrate access, egress, and solvent channels displayed varied bottle-neck radii. Taken together, our results indicate that SCH 66712 should inactivate these allelic variants, although varied flexibility and substrate binding-pocket accessibility may alter its interaction abilities.