Literature DB >> 32528175

Structure of a D2 dopamine receptor-G-protein complex in a lipid membrane.

Jie Yin1, Kuang-Yui M Chen2, Mary J Clark3, Mahdi Hijazi2, Punita Kumari1, Xiao-Chen Bai1, Roger K Sunahara4, Patrick Barth5, Daniel M Rosenbaum6.   

Abstract

The D2 dopamine receptor (DRD2) is a therapeutic target for Parkinson's disease1 and antipsychotic drugs2. DRD2 is activated by the endogenous neurotransmitter dopamine and synthetic agonist drugs such as bromocriptine3, leading to stimulation of Gi and inhibition of adenylyl cyclase. Here we used cryo-electron microscopy to elucidate the structure of an agonist-bound activated DRD2-Gi complex reconstituted into a phospholipid membrane. The extracellular ligand-binding site of DRD2 is remodelled in response to agonist binding, with conformational changes in extracellular loop 2, transmembrane domain 5 (TM5), TM6 and TM7, propagating to opening of the intracellular Gi-binding site. The DRD2-Gi structure represents, to our knowledge, the first experimental model of a G-protein-coupled receptor-G-protein complex embedded in a phospholipid bilayer, which serves as a benchmark to validate the interactions seen in previous detergent-bound structures. The structure also reveals interactions that are unique to the membrane-embedded complex, including helix 8 burial in the inner leaflet, ordered lysine and arginine side chains in the membrane interfacial regions, and lipid anchoring of the G protein in the membrane. Our model of the activated DRD2 will help to inform the design of subtype-selective DRD2 ligands for multiple human central nervous system disorders.

Entities:  

Mesh:

Substances:

Year:  2020        PMID: 32528175      PMCID: PMC7415663          DOI: 10.1038/s41586-020-2379-5

Source DB:  PubMed          Journal:  Nature        ISSN: 0028-0836            Impact factor:   49.962


Dopamine signaling through dopamine receptors

Dopamine is an essential neurotransmitter that regulates multiple physiological and cognitive functions including movement, reward, and emotion. DRD2 signaling is impaired in Parkinson’s disease, in which dopaminergic neurons are damaged and a deficit of dopamine leads to a combination of movement and psychiatric pathologies[4]. Parkinson’s disease is currently treated by dopamine replacement therapies, either with the prodrug L-dopa or analogs that can act directly on DRD2. Agonists such as bromocriptine have also been effectively used against pituitary adenomas[5]. Conversely, inhibitors of DRD2 are used in schizophrenia and psychosis, in order to reduce elevated dopaminergic signaling[6]. The dopamine family of G protein-coupled receptors (GPCRs) contains 5 subtypes, D1-D5, which have partially overlapping functions and pharmacology[1]. The receptors D1 and D5 couple to the stimulatory G protein Gs, while D2/D3/D4 couple mainly to Gi/o (inhibitory toward cAMP). At a sequence level, D1 and D5 share high overall homology (58% identity, 79% similarity) and D2/D3/D4 are likewise closely related (40–57% identity, 60–70% similarity)[7]. Despite the division between these receptor families, key residues that bind dopamine are highly conserved in all 5 receptors. For therapeutic ligands that target DRD2 in Parkinson’s and schizophrenia, subtype and/or family selectivity has been pursued to maximize efficacy and minimize off-target effects[8].

Structural biology of the dopamine receptors and G protein complexes

The structure of DRD2 bound to the inverse agonist risperidone was previously determined using an engineered T4L fusion protein[9], elucidating the overall inactive conformation and the molecular basis for inhibitor binding. The structure also demonstrated how risperidone can bind to D2-family receptors (DRD2, DRD3, DRD4) with selectivity over the D1 family (DRD1 and DRD5), based on interactions with an extended binding pocket including the extracellular loops. The structures of inhibitor-bound inactive conformations of other dopamine receptors DRD3[10] and DRD4[11] have also been reported using engineered fusion proteins. To date, no structure of an activated dopamine receptor has been reported. In contrast, activated GPCR-G protein complexes have recently been structurally characterized by cryo-electron microscopy (cryo-EM) for μ opioid[12] and CB1 cannabinoid[13] receptors, among others. Common to all these studies is the use of detergent micelles to solubilize the receptor and allow co-purification with G protein as single particles, with the caveat that the complexes lack the phospholipid bilayer environment important in receptor regulation and G protein recruitment. Indeed, early purification of endogenous DRD2 from bovine pituitary showed that the detergent-solubilized receptor has compromised ligand binding properties[14]. Previous work established the use of reconstituted HDL particles[15,16] (also known as lipid nanodiscs[17]) as effective host membranes for biochemical and biophysical studies of GPCRs and their G protein complexes. These systems have been exploited for cryo-EM structure determination of other membrane proteins such as ion channels[18], as well as for the full-length mGlu5 Family C GPCR[19] and the M2 acetylcholine receptor-arrestin complex[20], however no such structures have been reported for GPCR-G protein complexes. In this work, we report the first structure of an agonist-bound D2 dopamine receptor bound to Gi in an HDL/nanodisc system containing native lipids.

DRD2 construct design and cryo-EM structure determination

To obtain a human DRD2 construct optimized for active-state structure characterization, we computationally designed mutations that would shift the ligand-free receptor towards the active state by lowering the free energy of a modeled active conformation and raising the free energy of the inactive conformation (Methods, Extended Data Figs. 1a–e). The resulting construct bears 5 mutations relative to the wild-type human DRD2, none of which are at the orthosteric pocket. Because the native hDRD2 has a 120-residue intracellular loop 3 (ICL3) that could hinder purification and structural characterization, we de novo designed a 7-residue ICL3 compatible with Gi binding (Methods, Extended Data Fig. 1f). We also incorporated an N-terminal T4L (nT4L) domain, as previously described for β2AR[21], producing the final construct hDRD2EM. Inclusion of the nT4L and truncation/modification of ICL3 did not alter dopamine binding (Extended Data Fig. 2a) or agonist-mediated (dopamine or bromocriptine) stimulation of Gi (Extended Data Fig. 2b). As would be predicted based on its altered allosteric equilibrium, our designed receptor hDRD2EM has elevated constitutive activity, higher affinity for agonists such as dopamine and bromocriptine, and lower affinity for the inverse agonist spiperone (Extended Data Fig. 3a–c). In addition, the designed mutations confer enhanced thermostability compared to the wild-type hDRD2 (Extended Data Fig. 3e, f).
Extended Data Figure 1.

Rational design of receptor variants shifted in the active state. (a) Energy landscape of 3 receptor variants with various energy differences in the resting (ligand-free) state ΔGresting, in the allosteric coupling controlling the inverse agonist efficacy ΔGinverse agonist and the stabilization of the receptor inactive state upon ligand stimulus. (b) Active state occupancy as a function of ΔGresting and ΔGinverse agonist (the energy of inverse agonist binding), as determined by the Boltzmann law. (c) Sensitivity of the receptors to inverse agonist binding. (d-f). Rational design of DRD2 receptor stabilized in the active state. (d) DRD2 WT active state homology model. The 120 residue long intracellular loop 3 (ICL3) is represented as a dotted blue line. (e) DRD2 designed with residue microswitches shown in spheres stabilizing selectively the active state structure. (f) DRD2 further stabilized in the active state through a de novo designed 7-residue stable ICL3.

Extended Data Figure 2.

Effect of T4L, modified ICL3, and scFv16 on hDRD2WT ligand binding properties. (a) Competition binding of dopamine for [3H]spiperone on membranes containing human DRD2 wild type (hDRD2WT, blue) or human DRD2 modified with N-terminal T4L and de novo designed ICL3, nT4L-hDRD2(ΔICL3) (magenta). Log Ki ± S.E.M. = −4.6 ± 0.04 for hDRD2WT. Log Ki ± S.E. = −4.6 ± 0.04 for nT4L-hDRD2(ΔICL3). n = 4 biologically independent experiments on Sf9 cell membranes. (b) Agonist-stimulated [35S]GTPγS binding of nT4L-hDRD2(ΔICL3) (magenta symbols) compared to hDRD2WT (blue symbols). Dose-response using dopamine (closed magenta) or bromocriptine (open magenta) show log EC50 ± S.E.M. = −7.5 ± 0.02 M and −7.8 ± 0.06 M, respectively, for nT4L-hDRD2(ΔICL3) compared to log EC50 ± S.E.M. −7.0 ± 0.02 M for dopamine on hDRD2WT (open blue). n = 3 biologically independent experiments on HEK293T cell membranes. c) Saturation binding of [3H]spiperone to membranes containing DRD2WT with or without scFv16. Kd ± S.E.M. = 0.15 nM ± 0.03 nM for control. Kd ± S.E.M. = 0.14 nM ± 0.02 nM with scFv16. n = 4 biologically independent experiments on Sf9 cell membranes. (d) Competition binding of dopamine for [3H]spiperone on membranes containing DRD2WT with or without scFv16. Log Khigh ± S.E.M. = −6.4 ± 0.3 and Log Klow ± S.E.M. = −4.3 ± 0.08 for control. Log Khigh ± S.E.M. = −6.5 ± 0.2 and Log Klow ± S.E.M. = −4.3 ± 0.07 with scFv16. n = 4 biologically independent experiments on Sf9 cell membranes. In all panels, each data point is displayed as mean with error bars showing ± S.E.M., with individual data points for each repeat represented in small symbols.

Extended Data Figure 3.

Function and thermostability of hDRD2EM. (a) Gi stimulation by dopamine in Sf9 cells transformed with Gi and hDRD2 baculoviruses. Log EC50 ± S.E.M. = −6.3 ± 0.1 for hDRD2WT (blue). Log EC50 ± S.E.M. = −6.6 ± 0.3 for hDRD2EM (orange). n = 4 biologically independent experiments on Sf9 cell membranes. (b) Inhibition of dopamine-stimulated (10 μM) Gi activation by the inverse agonist spiperone. Log IC50 ± S.E.M = −8.6 ± 0.1 for hDRD2WT (blue). Log IC50 ± S.E.M = −6.3 ± 0.2 for hDRD2EM (orange). n = 4 biologically independent experiments on Sf9 cell membranes. (c) Modulation of basal Gi activation by dopamine (open), bromocriptine (closed) or spiperone (squares) with the hDRD2EM (orange symbols) relative to hDRD2WT(blue symbols). Log EC50 ± S.E.M. for dopamine (−6.6 ± 0.2 M), bromocriptine (−7.6 ± 0.1 M) and spiperone (−7.6 ± 0.2 M) for hDRD2EM. Log EC50 ± S.E.M. for dopamine (−7.0 ± 0.1 M) and bromocriptine (−7.8 ± 0.2 M) for hDRD2WT. n = 3 biologically independent experiments on Sf9 cell membranes. (d) Apparent functional stability of agonist-bound hDRD2 constructs assessed by measuring the fraction of partially-purified receptor activating Gi as a function of incubation time at 37ºC. Half-life ± S.E.M.= 27.6 min ± 1.0 min for hDRD2WT (blue). Half-life ± S.E.M. = 49.3 min ± 1.1 min for hDRD2 with 5 designed mutations and designed truncated ICL3 (orange). n = 3 biologically independent experiments (e) Apparent thermostability of agonist-bound hDRD2 constructs assessed by measuring the fraction of partially-purified receptor binding agonist as a function of temperature. Tm ± S.E.M.= 30.9 °C ± 0.9 °C for hDRD2WT(blue). Tm ± S.E.M. = 39.6 °C ± 0.3 °C for hDRD2 with 5 designed mutations and designed truncated ICL3 (orange). n = 3 biologically independent experiments. In all panels, each data point is displayed as mean with error bars showing ± S.E.M., with individual data points for each repeat represented in small symbols.

After screening commercially available agonists, we found that addition of the high-affinity agonist bromocriptine[3] permitted solubilization and purification of hDRD2EM. Importantly, bromocriptine displays full agonism towards hDRD2EM, and a lower EC50 relative to dopamine consistent with its higher affinity[22] (Extended Data Fig. 3c). We could efficiently co-purify a stable complex between hDRD2EM and Gi (Extended Data Fig. 4a, b), however constructs without the 5 active-state mutations or with intact full-length ICL3 displayed poor biochemical behavior upon detergent solubilization (Extended Data Fig. 4c). The purified hDRD2EM-Gi complex was further stabilized by addition of scFv16[12,13] that binds to the Gα and Gβ subunits away from the receptor interface, which also does not perturb ligand binding to wild-type hDRD2 (Extended Data Fig. 2b, c). This sample was reconstituted from detergent into nanodiscs containing a lipid mixture of palmitoyl-oleoyl-phosphatidylcholine (POPC), palmitoyl-oleoyl-phosphatidylglycerol (POPG), and cholesterol (3:2:1 molar ratio) along with the MSP1D1 scaffold peptide[23]. The final reconstituted complex was frozen in vitreous ice for single-particle cryo-EM analysis. We collected 9760 images on a Titan Krios microscope and obtained a reconstruction of the reconstituted hDRD2EM-Gi complex at 3.7 Å resolution. Representative class averages and data analysis flowchart are included in Extended Data Figure 5 & 6, and associated statistics are included in Extended Data Table 1. The cryo-EM maps show strong density for the TM domains of DRD2 including sidechains (Extended Data Fig. 6d), without large variation in local resolution (Extended Data Fig. 5c).
Extended Data Figure 4.

Purification and reconstitution of the hDRD2EM-Gi complex. (a) Superose 6 gel filtration profile of hDRD2EM-Gi-scFv16-HDL purified by M1 FLAG affinity chromatography. (b) Coomassie stained PAGE of the isolated peak fraction from gel filtration. (c) Overlaid Superose 200 gel filtration profiles (in LMNG detergent) of different hDRD2 constructs purified by M1 FLAG chromatography with saturating bromocriptine present. nT4L-hDRD2–5mut has the 5 thermostabilizing mutations of the EM construct, but has full-length ICL3. nT4L-hDRD2-ΔICL3 has the wild-type TM region (no mutations) but the truncated ICL3. Receptors were purified as in Methods, but in the absence of co-transduced G protein baculovirus.

Extended Data Figure 5.

Cryo-EM data analysis. (a) Representative 2D class averages. (b) Gold-standard FSC curves of the 3D reconstructions. (c) Image processing procedure with the final maps colored based on local resolution.

Extended Data Figure 6.

3D classification with signal subtraction, model validation and representative map density. (a) Focused 3D classification subtracting T4L, AH domain, scFv16, and nanodisc density. (b) Focused 3D classification subtracting all but the G protein. (c) FSC curve showing the model-map correlation (focused map without T4L, AH domain, scFv16, and nanodisc density). (d) representative DRD2 TM and ligand density.

Extended Data Table 1.

Cryo-EM data collection, refinement and validation statistics.

hDRD2EM-DNGαi1Gβ1Gλ2-scFv16-HDL (EMD-21243)hDRD2EM-DNGαi1Gβ1Gλ2 focused (-nT4L, AH domain, scFv16 and HDL) (EMD-21244) (PDB 6vms)DNGαi1Gβ1Gλ2 focused (-Receptor, scFv16 and HDL) (EMDB-21245)
Data collection and processing
Magnification105,000105,000105,000
Voltage (kV)300300300
Electron exposure (e/Å2)646464
Defocus range (μm)1.5–2.51.5–2.51.5–2.5
Pixel size (Å)0.830.830.83
Symmetry imposedC1C1C1
Initial particle images (no.)2,001,0232,001,0232,001,023
Final particle images (no.)783,98443,73358,356
Map resolution (Å)3.83.73.7
 FSC threshold0.1430.1430.143
Refinement
Initial model used (PDB code)6dde, 1gp2
Model Resolution (Å)3.9
 FSC threshold0.5
Map sharpening B factor (Å2)−155
Model composition
 Non-hydrogen atoms8953
 Protein residues1135 (8910 atoms)
 Ligands1 (43 atoms)
B factors (Å2)
 Protein81.57
 Ligand67.82
R.m.s deviations
 Bond lengths (Å)0.005
 Bond angles (°)0.703
Validation
 MolProbity score1.37
 Clashscore5.56
 Poor rotamers (%)0
Ramachandran plot
 Favored (%)97.68
 Allowed (%)2.32
 Disallowed (%)0

Global structure of the DRD2-Gi complex reconstituted into a phospholipid bilayer

The overall reconstruction of DRD2-Gi in phospholipid bilayers is shown in Figure 1a, with the resulting model for receptor, ligand and Gi heterotrimer in Figure 1b. The structure is similar in global arrangement to previous GPCR-Gi complexes, with a notable difference in the asymmetric position of the receptor in the rHDL/nanodisc particle (Fig. 1a). From the receptor on one side to the G protein-bilayer contacts on the other, the protein complex fits tightly into the rHDL/nanodisc. Previous biochemical studies demonstrated that monomeric β2 adrenergic receptor (β2AR) embedded in rHDL fully retains its pharmacological properties, couples to Gs and stimulates nucleotide exchange[15]. Our reconstruction provides a structural basis for such activity, as the G protein footprint is accommodated within the ~80 Å inner diameter of the rHDL particle. Compared to the CB1-Gi complex in micelles[13], DRD2 in our reconstruction is tilted 10 degrees relative to the membrane-proximal surface of Gi (i.e. Gα’s N-terminal helix αN) which positions the receptor’s helix 8 in the membrane interfacial region (a larger translation of receptor relative to Gi is present in comparison with μOR-Gi[12]) (Fig. 1c). These differences reflect the fact that the phospholipid membrane places a physical constraint on the elements of the receptor and G protein that can come into contact. The structure of hDRD2EM in the complex is in close agreement with our computationally designed active-state conformation (Fig. 1d, Cα rmsd 2.1 Å overall and 1 Å for the TM helices).
Figure 1.

Overall structure of the hDRD2EM-Gi complex

(a) Cryo-EM reconstruction of the complex (density of T4L, AH domain, and scFv16 not shown). Receptor density is in purple, bromocriptine is yellow, Gα is orange, Gβ is cyan, Gγ is pink.

(b) Model of complex as a cartoon, with the bromocriptine agonist as spheres with yellow carbons. scFv16, T4L, and AH domain are not shown.

(c) Superposition of the hDRD2EM-Gi complex (purple) and the CB1-Gi complex (light orange).

(d) Comparison of the designed active state conformation of DRD2 (green cartoon) with the experimental structure of hDRD2EM (receptor as purple cartoon). The root mean squared deviation (rmsd) for the Cα positions in the superposition is 2.1 Å.

Agonist binding to DRD2

Bromocriptine is an ergot alkaloid that acts as a full agonist at DRD2. In studies using subtype-selective radioligand antagonists as probes, bromocriptine displayed 2–5 nM affinity towards mammalian DRD2 and roughly 100-fold selectivity over DRD1, DRD4, and DRD5 (with little discrimination against DRD3)[22]. The residues interacting directly with bromocriptine in the active-state structure are shown in Figure 2a, and the cryo-EM density surrounding the ligand is displayed in Figure 2b. The ergoline core of bromocriptine is sandwiched between TM3 and TM6, and contains atoms analogous to the dopamine skeleton (Fig. 2c). As in other catecholamine receptors[24], the binding pocket is flanked on one side by an aspartate on TM3 (Asp1143.32) forming a salt bridge with the agonist amino group, and on the other side by a cluster of serines on TM5 (Ser1935.42, Ser1945.43, Ser1975.47), which are important for dopamine binding and activation[25,26]. The dopamine-like core of bromocriptine packs against Phe3896.51 and Phe3906.52 on TM6, in good agreement with biochemical precedent[27]. The bicyclic tripeptide group of the ligand extends toward the extracellular surface, primarily contacting residues on TM7 and extracellular loop 2 (ECL2) including Ile184ECL2 previously characterized as lining the ligand binding pocket[28]. The ECL2 sequences of DRD1 and DRD5 are completely divergent from DRD2, which helps to explain the agonist’s selectivity against those receptors. Bromocriptine and related ergopeptines are unique among DRD2 agonists in displaying single-site high affinity binding that is independent of G protein coupling (i.e. with or without the presence of GTP)[22,29,30]. The bivalent nature of bromocriptine’s interaction with DRD2 (Fig. 2a) prevents closure of the extracellular lid of the receptor, which was shown to be critical for high-affinity agonist binding in β2AR[31]. The extensive interface of bromocriptine with DRD2 may allow for slow dissociation and resulting low-nM affinity even when the receptor is in an inactive conformation. Bromocriptine’s binding overlaps with the interaction of risperidone (Fig. 2c), however the inverse agonist extends into a deeper subpocket in the inactive conformation[9], which contracts in the active conformation (Extended Data Fig. 7).
Figure 2.

Binding of bromocriptine to DRD2

(a) Contact residues (sticks with purple carbons) within 4.1 Å of bromocriptine (yellow carbons). Val1153.33 is not shown to better illustrate the polar interaction between Asp1143.32 and bromocriptine’s amine (dotted line, 2.8 Å).

(b) Cryo-EM density surrounding the bromocriptine ligand.

(c) Overlay of bromocriptine and risperidone (transparent green sticks) when superimposing the receptor in the agonist-bound active conformation (this work) and the inverse agonist-bound inactive conformation (PDB 6cm4). Box denotes atoms of bromocriptine analogous to dopamine (shown above as 2D structure).

Extended Data Figure 7.

Surfaces of the orthosteric pocket in the active and inactive conformations of DRD2. (a) Solvent accessible surface for the active conformation (purple) bound to bromocriptine (yellow sticks). (b) Solvent accessible surface for the inactive conformation (orange) bound to risperidone (green sticks), showing the deeper subpocket (PDB 6cm4).

Propagation of conformational changes from dopamine to Gi

Comparison of the activated receptor in the Gi complex with the inverse agonist-bound structure shows how changes at the dopamine binding pocket are translated through the receptor to facilitate nucleotide exchange on the bound Gi. Bromocriptine binding at the orthosteric pocket occurs with rotation of the upper half of TM5 toward TM4 (Fig. 3a), concomitant with rigid body shifts of the extracellular halves of TM6 and TM7 inward toward the ligand and movement of ECL2 over the pocket. Residues Phe1895.38 and Ile184ECL2 play an important role in this conformational change, since the swivel of the Phe1895.38 phenyl group makes room for a downward shift of the Ile184ECL2 sidechain. Recent studies have implicated these two residues’ movement in DRD2 activation, as well as in variable interaction with arrestin-biased agonists[32,33]. Repacking of sidechains in these upper TM segments (Fig. 3a) results in a downward shift in Trp3866.48, formation of an active PIF transmission switch (Fig. 3b), and outward movement of the bottom half of TM6 by 11.9 Å (measured between E3686.30 Cα’s), all hallmarks of GPCR activation[34]. The designed V3786.40Y mutation at the apex of the Gi binding site helps to stabilize the TM6-out active conformation by packing against Tyr4267.53 and Arg1323.50 of the DRY motif (Extended Data Fig. 8), interactions that are not present in the inactive conformation.
Figure 3.

Propagated changes in DRD2 activation

(a) Rotation of the upper half of TM5, comparing the active conformation (purple sticks, transparent gray cartoon) and the inactive conformation (light orange sticks and transparent cartoon, PDB 6cm4). Bromocriptine is shown as yellow sticks and risperidone is in green sticks.

(b) PIF transmission switch in the active conformation of DRD2 (upper panel, purple sticks and transparent spheres) and the active conformation of β2AR (lower panel, PDB 3sn6, green sticks and transparent spheres).

(c) Surface representation of the receptor-G protein interface in the hDRD2EM-Gi complex. Receptor is purple, Gα is orange, Gβ is cyan, and Gγ is pink.

(d) Comparison of the Gα subunit in the complex (purple cartoon) and Gα in the GDP-bound free Gi (yellow cartoon). Receptor in the complex is shown as a purple transparent surface.

Extended Data Figure 8.

Interaction stabilizing the designed hDRD2EM active state conformation. (a) Zoomed view of the TM6-TM7 interface on the cytoplasmic side of the receptor. Atomic interactions involving the designed Y6.40 were accurately predicted in the design model. The conformation of the R3.50 is displaced in the experimental structure due to the binding of Giα which was absent in the model. Experimental structure is in purple and predicted DRD2 conformation is in green. (b, c) Conformational energy landscape of the WT (b) and designed (c) DRD2 receptors constructed from multiple all-atom MD simulation trajectories of the receptor active state homology models (orange dot) in absence of bound agonist. The conformational space is reported along 2 canonical active state structure metrics (TM3–6 distance and NPxxY rmsd to inactive state). The DRD2 WT relaxes primarily to an inactive state conformation while the designed variant remains in partially active state conformation displaying large TM3–6 distance. The antagonist-bound DRD2 crystal structure (PDB 6CM4) is used as a reference for the inactive state (green square).

The interaction between DRD2 and Gi in rHDL/nanodiscs buries 2,856 Å2 (Fig. 3c). Like the CB1-Gi[13], μOR-Gi[12], or β2AR-Gs[35], the most extensive contacts occur between TM5, TM6, TM7, and ICL2 of the receptor and the α5 helix of the G protein. As in the μOR-Gi complex (but not in CB1-Gi), DRD2’s ICL2 makes multiple contacts with the Gαi αN helix that extends along the membrane surface. Similar to the β2AR-Gs crystal structure in monoolein bilayers, DRD2 also contacts the Gβ subunit through residues on ICL1 and ICL2 (Fig. 3c). Engagement by DRD2 causes a 5 Å translation and 60º rotation of the α5 helix relative to the GDP-bound free G[36], similar to other complexes, resulting in movement of the TCAT motif on the β6-α5 loop away from the nucleotide and movement of the P-loop into space sterically overlapping with GDP binding (Fig. 3d). During 3D classification of our cryo-EM dataset (Extended Data Fig. 5), we also observed weaker density corresponding to the alpha-helical (AH) domain of Gαi, and after focused refinement on the G protein (Extended Data Fig. 6b), we could dock the AH domain into density as a rigid body (Fig. 3d, Extended Data Fig. 9a). The AH domain in DRD2-Gi is positioned similarly to its site in the rhodopsin-Gi[37] (Extended Data Fig. 9b–d), implying that these two structures may represent a low-energy conformation within nucleotide-free GPCR-Gi complexes despite the AH domain’s well-characterized flexibility[21,35].
Extended Data Figure 9.

AH domain density and comparison with rhodopsin-Gi. (a) AH domain density in the current structure. (b) AH domain density in the rhodopsin-Gi complex (PDB 6cmo). (c) Alignment of Gαi subunit in hDRD2EM-Gi (purple) and rhodopsin-Gi complex (blue). (d) Alignment of AH domain in hDRD2EM-Gi (purple) and rhodopsin-Gi complex (blue).

Interactions specific to the DRD2-Gi complex in the rHDL/nanodisc

G protein subunits are specifically acylated, which is essential for their trafficking to the plasma membrane and the formation of functional heterotrimers. Gαi is first myristoylated at the N-terminal Gly2 residue[38] and then reversibly palmitoylated at the following Cys3 position[39], whereas Gγ is geranylgeranylated at the C-terminal CAAX motif which is important for membrane association of Gβγ and Gα recruitment[40]. Disruption of the lipid modifications on Gαi by mutagenesis impairs agonist-dependent inhibition of adenylyl cyclase[41]. In our reconstruction of DRD2-Gi in rHDL/nanodiscs, the αN helix runs along the outer surface of the bilayer (Fig. 4), where positively charged residues Lys10, Lys17, and Arg24 project up to interact with the polar membrane headgroups. Weak density could be observed at the ends of the Gαi and Gγ subunits by cryo-EM, however we could not accurately model the acylated residues or lipid modifications themselves (presumably due to disorder). What has not been apparent in structures of GPCR-G protein complexes reported to date is the lipid association of Gβ, other than its tight association with the acylated Gγ. In the reconstituted DRD2-Gi complex, we observe density for basic residues Arg42, Arg46 and Arg52 of Gβ, which interact with the polar lipid headgroups (Fig. 4). Likewise, Arg62 and Lys64 on the C-terminal portion of Gγ are also well-ordered and interact with the lipid headgroups of the phospholipid bilayer. These collective polar interactions between the Gi heterotrimer and the lipid headgroups, together with Gαi and Gγ acylation, highlight the critical role of the phospholipid bilayer in receptor-G protein coupling.
Figure 4.

Interactions of the hDRD2EM-Gi complex with the rHDL/nanodisc membrane.

Cryo-EM density for the hDRD2EM-Gi complex reconstituted in rHDL/nanodisc membrane highlighting densities corresponding to the receptor (purple), Gαi1 (orange), Gβ1 (cyan), Gγ2 (pink) and rHDL/nanodisc (gray). (top right inset, red) illustrates the position of helix 8 of DRD2 (tan) embedded in the phospholipid bilayer. Also rendered is the putative palmitoyl-cysteine based on the palmitoylated helix 8 of the β2AR (PDB 2RH1). (bottom right inset, green) The Gαi1-rHDL/nanodisc interface is highlighted by the strong EM density on the extreme amino-terminus of Gαi1 (orange) and Gγ2 (magenta). Also illustrated is the Gαi1 C-terminal (α5-helix) interacting with the DRD2 core. (left inset, light blue) Interaction of the N-terminal region of Gβ1 with the phospholipid bilayer. Highlighted are residues Arg46 (R46) and Arg52 (R52).

The amphipathic helix 8 in Class A GPCRs often contains a cysteine residue that can be reversibly palmitoylated, represented by Cys443 in DRD2. Lipidation of the receptor is important for partitioning of helix 8 into the membrane[42], and in cases such as β2AR depalmitoylation is regulated by agonist stimulation[43]. We were unable to model the palmitate moiety in DRD2, most likely because it is disordered within the membrane (although it is possible that some fraction of the palmitate was hydrolyzed during purification). Strikingly, helix 8 in the activated DRD2-Gi complex is entirely embedded within the phospholipid membrane interfacial region (Fig. 4), in agreement with model studies of rhodopsin helix 8 peptides in membranes using fluorescence quenching by brominated lipids[44]. The attitude of helix 8 within the membrane positions residues Phe4297.56, Asn430, and Glu432 at the TM7-helix 8 junction to participate in the interface with the C-terminus of Gαi’s α5 helix. The degree to which helix 8 is embedded in the bilayer again highlights the importance of the lipid membrane in receptor structure and function. DRD2, like many GPCRs, has an abundance of positively charged basic residues on the cytoplasmic side of the membrane, in accordance with the ‘positive-inside rule’ for α-helical membrane proteins[45]. We observe strong density for a number of these arginine and lysine sidechains in the inner leaflet interfacial region (e.g. Arg611.59, Arg1514.41, Arg2175.66, Lys3706.32, Lys439), which may reflect energetically favorable interactions with lipid headgroups and conformational restriction previously characterized in MD simulations of basic residues in model TM helices[46]. These basic residues create a positively charged ‘ring’ at the membrane surface that is complementary to electrostatically negative surfaces on the Gi heterotrimer (Extended Data Fig. 10).
Extended Data Figure 10.

Electrostatic surfaces in the hDRD2EM-Gi complex. (a) Surface potential of the DRD2 and Gi components shown separately. (b) Surface potential of the complex, with scale in eV. Calculation was done using the APBS plugin in Pymol.

The distinctive features of this complex may be specific to the DRD2-Gi interaction, but more likely reflect general features of GPCR-G protein complexes in membranes. As we elucidate additional structures of receptors in a more native-like phospholipid bilayer environment (rHDL/nanodiscs), more light will be shed on how the membrane contributes to the molecular mechanism of G protein coupling and specificity.

Methods

Computational modeling of DRD2 active and inactive states

The inactive and active state structures of DRD2 WT were modeled from the close homolog inactive state structure of DRD3 (3PBL) and the more distant homolog active state structure of β2AR, respectively, using the software iPHoLD[47]. Low energy representative models of the DRD2 WT were selected for design calculations, as previously described.

Computational design of residue microswitches

The TM region of DRD2 was subjected to complete multi-state design using a fitness function selecting residue microswitches stabilizing the active state structure while decreasing the allosteric coupling between the extracellular ligand binding site and the inactive state conformation of the intracellular G protein binding site. Conformational stability was assessed using an all-atom energy function developed for membrane protein modeling and design. Allosteric couplings were approximated by the correlated dynamic fluctuations between allosteric residues using normal mode calculations[48].

Computational de novo design of minimal intracellular loop 3

The kinematic loop design protocol of Rosetta was used to de novo design short loop structures and sequences connecting the cytoplasmic ends of TM5 and TM6. Loop structures and sequences were selected to maximize stability and optimize conformation for Gαi binding. Loop lengths of three to nine amino acids were tested and the predicted models were clustered, selected by energy and structurally analyzed for distortions of the TM5/6 helical tips. Several designed loops of five to seven residues adopted stable helical structures enabling proper Gαi binding. Selected sequences from these loops with further refined with RosettaMembrane on active state DRD2 models onto which a peptide corresponding to the Gαi C-terminal α5 helix was docked. The final selected 7 residue long designed ICL3 was predicted to retain native receptor-G protein TM contacts while stabilizing the DRD2 active conformation.

Mutagenesis and expression of receptors in HEK293 cells for stability assays

Quickchange PCR mutagenesis (Stratagene) was performed on the HA-tagged human DRD2 gene in the pcDNA3.1(+) vector. Sequenced mutant plasmids were transiently transfected using lipofectamine (Invitrogen) into HEK293T cells for DRD2 as previously described[48]. Briefly, 5 × 106 cells were plated on 10 cm tissue culture plates and grown overnight. This was followed by transfection with Lipofectamine 2000 (Invitrogen) and 2 μg of DNA per plate. After 24 hours, the cells were washed with PBS and grown in standard growth medium (DMEM supplemented with L-glutamine (2 mM), penicillin (100 mg/mL), streptomycin (100 mg/mL), and fetal bovine serum (10%)) for an additional 24 hours prior to harvesting. The HEK293T cell line used in this study was obtained from ATCC (catalog # CRL-3216). The cell line was not further authenticated, and was not tested for mycoplasma contamination.

Cloning and expression of hDRD2EM and Gi heterotrimer

The designed hDRD2EM construct contains 5 mutations: T2055.54I, M3746.36L, V3786.40Y, V3816.43L and V4217.48I. Residues 1–29 were replaced with T4 lysozyme followed by a 2-residue (AlaAla) linker, and residues 222–360 in ICL3 were replaced with residues LVNTN, as designed in the optimized 7-residue ICL3 described above (with one wild-type amino acid flanking each side of the insertion). The resulting construct was cloned into a modified pFastBac baculovirus expression vector with a HA signal sequence followed by a FLAG tag at the N terminus. A dominant-negative rat Gαi1 (DNGαi1) was made by incorporating 3 mutations S47N[49], G203A[50] and E245A[51] which can stabilize the GPCR-G protein complex as previously reported[52]. Rat Gαi1 only differs from human Gαi1 by a single residue at position 98 in the AH domain where it is Ala in rat and Ser in human. A second pFastBac baculovirus was made with rat DNGαi1. An additional pFastBac-Dual baculovirus was made with human Gβ1 and human Gγ2 genes. An 8xHis tag was placed at the N-terminus of Gγ2. hDRD2EM, DNGαi1, Gβ1, and His8-tagged Gγ2 were co-expressed in Spodoptera frugiperda (Sf9) with the addition of all three baculoviruses (ratio of 1 hDRD2EM:2 DNGαi1:1 Gβ1Gγ2) to Sf9 cells at a density of 3 × 106 ml−1 with 1 μM bromocriptine mesylate (Sigma) added to the media during growth. After 48 hours, cells were harvested and stored at −80 °C for future use. The Sf9 cell line used in this study was obtained from Expression Systems Inc. (catalog # 94–001S). The cell line was not further authenticated, and was not tested for mycoplasma contamination.

Cloning, expression and purification of scFv16

The single chain antibody scFv16 was described in previous GPCR-Gi heterotrimer complexes[12,13]. A synthesized scFv16 DNA fragment was cloned into a modified pFastBac baculovirus expression vector (Invitrogen) with a melittin signal sequence at the N terminus and a decahistidine tag at the C terminus. The resulting construct was expressed in Sf9 cells with the Bac-to-Bac system (Invitrogen). scFv16 was purified as previously described[12]. In brief, the secreted scFv16 in the cell culture media was separated from Sf9 cells by centrifugation and 10 mM Tris buffer pH 8.0 was added for balancing pH. Then 1 mM Ni2SO4 and 5 mM CaCl2 were added to quench chelating agents. The solution containing scFv16 was loaded onto Ni-NTA affinity resin by a siphon. The column was washed with 20 column volumes of buffer consisting of 50 mM Hepes, pH 7.4, 150 mM NaCl and 10 mM imidazole. Protein was eluted from Ni-NTA resin with 10 column volumes of buffer consisting of 50 mM Hepes, pH 7.4, 150 mM NaCl and 250 mM imidazole. The eluate was concentrated and run on a Superdex 200 gel filtration column. The peak corresponding to monomeric scFv16 was collected, concentrated and frozen until further use.

Purification of the complex of hDRD2EM and Gi heterotrimer

Sf9 cell pellets were resuspended and lysed in a hypotonic buffer consisting of 10 mM Hepes pH 7.4, 160 μg/ml benzamidine, 100 μg/ml leupeptin,1 mM MgCl2, 0.1 mM TCEP, Apyrase (0.5 mU/ml, NEB) and 1 μM bromocriptine mesylate (Sigma) for 30 min at 4 °C. Lysed cells were spun down and were homogenized using a Dounce tissue grinder (Wheaton) in a solubilization buffer consisting of 50 mM Hepes pH 7.4, 150 mM NaCl, 160 μg/ml benzamidine, 100 μg/ml leupeptin, 1 mM MgCl2, 0.1 mM MnCl2, 1% (w/v) lauryl maltose neopentyl glycol (LMNG, Anatrace), 0.1% Na Cholate, 0.1% cholesteryl hemi-succinate (CHS, Steraloids), 10% glycerol, 0.5 mM TCEP, apyrase (15 mU/ml, NEB), λ phosphatase (10 U/ml, NEB) and 10 μM bromocriptine mesylate (Sigma). Solubilization proceeded for 1 hour at 4 °C, followed by centrifugation at 35,000 rpm for 30 min at 4 °C. The supernatant was incubated with M1 anti-Flag affinity beads (Sigma) in batch-binding mode overnight at 4 °C in the presence of 2 mM CaCl2. After binding, the M1 beads were spun down and transferred into a glass column. Then beads were washed with 10 volumes of buffer consisting of 50 mM Hepes, pH 7.4, 150 mM NaCl, 0.05% LMNG, 0.005% sodium cholate, 0.005% CHS, 5% glycerol, 1 mM MgCl2, 2 mM CaCl2, 0.5 mM TCEP and 5 μM bromocriptine mesylate. The complex was eluted from M1 beads with buffer consisting of 50 mM Hepes pH 7.4, 150 mM NaCl, 0.05% LMNG, 0.005% Na Cholate, 0.005% CHS, 5% glycerol, 0.5 mM TCEP, 5 μM bromocriptine mesylate, 200 μg/ml FLAG peptide and 5 mM EDTA. Finally, the eluate was concentrated and applied onto Superdex 200 size exclusion column (GE Healthcare). The peak corresponding to the hDRD2EM-Gi heterotrimer complex was collected and mixed with scFv16 at a molar ratio 1:1.2 and incubated for 30 min on the ice. The isolated complex of hDRD2EM, Gi heterotrimer and scFv16 was concentrated and used in further experiments.

Reconstitution into rHDL/nanodiscs

The complex of hDRD2EM, Gi heterotrimer and scFv16 was inserted into reconstituted high density lipoprotein particles (rHDLs), commonly known as nanodiscs, as previously described[15]. In brief, the complex was mixed with MSP1D1 and a lipid mixture (POPC: POPG: Cholesterol = 3:2:1) at molar ratio 1:5:400 and the mixture was incubated for 1 hr at 4 °C. Bio-beads SM2 (50 mg/ml, Bio-Rad) were added into the mixture and incubated for 1 hr at 4 °C. The sample was transferred to fresh Bio-beads two additional times with the last incubation proceeding for 8 hr. The reconstitution mixture was spun down and the supernatant was loaded onto an M1 anti-Flag affinity column to remove empty nanodiscs. The eluate from the M1 beads was concentrated and applied onto a Superose 6 gel filtration column (GE Healthcare) (Extended Data Fig. 4a). The peak corresponding to the complex reconstituted in rHDL/nanodisc was collected and analyzed by SDS-PAGE gel electrophoresis (Extended Data Fig. 4b).

Cryo-EM data acquisition

The hDRD2EM-Gi-scFv16 rHDL/nanodisc complex was concentrated to 6 mg/ml and supplemented with 2 mM fluorinated fos-choline-8 (Anatrace) immediately prior to freezing. Cryo-EM grids were prepared by applying 3.5 μl of this sample to glow-discharged Quantifoil R1.2/1.3 300-mesh gold holey carbon grids (Quantifoil Micro Tools GmbH, Germany), blotted for 4.5 sec under 90% humidity at 4 °C and plunge frozen in liquid ethane cooled by liquid nitrogen using a Mark IV Vitrobot. Cryo-EM grids were first imaged on a Talos Arctica microscope (FEI) operated at 200kV with a K3 direct electron detector (Gatan) using a Volta phase plate following a protocol as previously described[53]. SerialEM software was used for automated data collection. Images were recorded at a nominal magnification of ×45,000, corresponding to a pixel size of 0.4385 and defocus set at 0.75 μm. Each stack image was dose-fractionated to 50 frames for a total of 2.5 s under a dose rate of 22 e-/Å2/s, resulting in a total dose of 55 e-/Å2. For high-resolution data collection, images were recorded on a Titan Krios microscope (FEI) operated at 300kV with a K3 direct electron detector (Gatan) in super-resolution mode using a slit width of 20 eV on a GIF-Quantum energy filter. SerialEM[54] was used for automated data collection. Images were recorded at a nominal magnification of ×105,000, corresponding to a pixel size of 0.415 Å and the defocus set in a range from 1.5 to 2.5 μm. Each stack image was dose-fractionated to 40 frames for a total of 2 s under a dose rate of 32 e-/Å2/s, resulting in a total dose of 64 e-/ Å2.

Image processing

All datasets were processed in Relion 3.0[55] following the same general protocol as previously described[56]. 697 dose-fractioned images collected on a Talos Arctica microscope were gain normalized, 2 × 2 Fourier binned, motion corrected, dose-weighted, and summed using MotionCorr2[57]. Contrast transfer function parameters were estimated using GCTF[58]. Approximately 1,000 particles were picked manually and subjected to 2D classification. Representative projections of the complex were selected as templates for automated particle picking from all images. The 233,645 extracted particles were 4 × 4 binned and subjected to 2D classification. 2D classes displaying clear domain boundaries in complexes were picked and used to generate an ab initio model in Relion. A total of 42,956 particles were finally selected for 3D classification using the ab initio model as the reference. After 3D classification, the classes showing good secondary structure features were selected, combined and re-extracted at the original pixel size of 0.877 Å. After 3D refinement, the resulting 3D reconstitution from 11,290 particles yielded a map at 19 Å resolution (gold-standard FSC criterion[59]). The same procedure was applied to process data collected on Titan Krios. 9,760 dose-fractioned images collected on Titan Krios were gain normalized, 2 × 2 Fourier binned, motion corrected, dose-weighted, and summed using MotionCor2[57]. Contrast transfer function parameters were estimated using GCTF[58]. Approximately 1,000 particles were picked manually and subjected to 2D classification. Representative projections of the complex were selected as templates for automated particle picking from all images. The 2,001,023 extracted particles were 4 × 4 binned and subjected to 2D classification. A total of 1,341,844 particles were finally selected for 3D classification using the map reconstituted from the data collected on Talos Arctica as the reference. After 2 rounds of 3D classification, the classes showing good secondary structure features were selected, combined and re-extracted at the original pixel size of 0.83 Å. After 3D refinement and postprocessing, the resulting 3D reconstitution from 783,984 particles yielded a map at 3.8 Å resolution (gold-standard FSC criterion). To identify a subset of these particles with better density in the transmembrane domain, we carried out a focused 3D classification[60] subtracting the density of T4 lysozyme, AH domain in the Gαi1 subunit, scFv16 and the belt-like density from the MSP1D1 nanodisc. The subsequent 3D classification separated out a class with 43,733 particles showing better density in the transmembrane domain (Extended Data Fig 6a). The selected particles from this class were subjected to 3D refinement and postprocessing yielding a map with 3.7 Å resolution. To improve the density of the AH domain in the Gαi1 subunit, we carried out another focused 3D classification subtracting the density of T4 lysozyme, DRD2, scFv16 and the belt-like density from the MSP1D1 nanodisc. The subsequent 3D classification separated out a class with 58,356 particles showing partial secondary structures within the AH domain (Extended Data Fig. 6b), which allowed us to dock a crystal structure of the AH domain (PDB: 1GP2).

Model building

The Cryo-EM structure of the μOR–Gi complex with scFv164 (PDB: 6ddf) was used as an initial model for model rebuilding and refinement. A polyalanine model was made from the μOR structure, and together with structures of Gi and sFv16 (PDB: 6ddf) was docked into the EM electron density map in Chimera[61]. The resulting model was subjected to autobuilding in Buccaneer[62], iterative building in Coot[63] and refinement in Phenix.real_space_refine[64]. Initial coordinates and refinement parameters for the ligand bromocriptine were prepared with Grade web server (http://grade.globalphasing.org). The structure of the AH domain from the GDP-bound Gi crystal structure[36] (PDB: 1GP2) was docked into the focused refined EM map. Manual editing in Coot and refinement in Phenix.real_space_refine were carried out to reduce the clashes of the docked AH domain with the Gαi Ras-like domain and Gβ subunit. The density of T4 lysozyme was poor and not modelled. FSC curves were calculated based on the model and maps in Phenix.validatiion_cryoem (Extended Data Fig. 6c). Representative cryo-EM density for the receptor is displayed in Extended Data Fig. 6d. MolProbity[65] was used to evaluate the final structures. In the Ramachandran plot, 97.7% and 2.3% of residues were in favored and allowed regions, respectively. The statistics for data collection and refinement are included in Extended Data Table 1. The consensus cryo-EM density map and two focused maps have been deposited in the Electron Microscopy Data Bank under accession codes EMD-21243, EMD-21244, and EMD-21245. Atomic coordinates have been deposited in the PDB under accession code 6vms.

Radioligand binding assays

Membranes were prepared from Sf9 cells co-infected with Gαi1β1γ2 and wild-type human DRD2long (DRD2WT) as described[31]. Membranes (7–10 μg) were pretreated with or without antibody fragment scFv16 (1 μg/10 μg membranes) for 30 min at room temp and then incubated for 2 hr at 30oC in binding buffer (20 mM HEPES, pH 7.4, 100 mM NaCl, 1 mM ascorbic acid, 0.001% BSA, 100 nM GDP) and 0.03 – 1.6 nM [3H]spiperone (for saturations) or with 0.13 nM [3H]spiperone (for competitions) along with varying concentrations of dopamine (or with 20 μM (+)butaclamol to determine non-specific binding). Samples were subject to rapid filtration through GF/B Unifilter™ plates (Packard) and rinsed with ice cold binding buffer to remove free [3H]spiperone. Filter plates were dried before adding Microscint 0™ (Perkin Elmer) and counting bound [3H]spiperone using a TopCount™ (Perkin Elmer). All data were analyzed using Prism (Graphpad Software Inc, San Diego CA).

[35S]GTPγS binding assay

Membranes were prepared from Sf9 cells expressing Gαi1β1γ2 and hDRD2 wild-type or hDRD2EM. Typically, membranes (~30 μg) were pretreated with GDP (final assay concentration of 1 μM) in G assay buffer (20 mM HEPES, pH 7.4, 100 mM NaCl, 10 mM MgCl2, and 1 mM ascorbic acid) and different concentrations of ligand(s) for 10 min at room temperature before adding [35S]GTPγS (for a final concentration of 0.1 nM). The assay was incubated at 30oC for 1 hr before stopping by rapid filtration through GF/B Unifilter plates and washing with ice-cold assay buffer. Filter plates were dried before adding Microscint 0™ and counting bound [35S]GTPγS using a TopCount™ (Perkin Elmer).

Receptor stability

For stability experiments, partially purified DRD2 receptor (or mutant) was prepared from membranes from transfected cells expressing either hDRD2WT or hDRD2EM as previously described[48]. Briefly, transfected HEK293 cells were gently harvested in PBS solution, pelleted by low speed centrifugation and resuspended in cold hypotonic buffer (1 mM Tris-HCL, pH 6.8, 10 mM EDTA, protease inhibitor cocktail). Cells were disrupted through a 26-gauge needle three times and the lysate was layered onto a 38% sucrose solution in buffer A (150 mM NaCl, 1 mM MgCl, 10 mM EDTA, 20 mM Tris-HCl, pH 6.8, protease inhibitors) and spun at 15,000 rpm for 20 minutes at 4 °C using a swinging bucket rotor (SW28). The fluffy interface band containing enriched membranes were diluted in Buffer A and pelleted by ultracentrifugation. Membrane pellets were resuspended in Buffer A, flash frozen in liquid N2 and stored at −80 °C until further use. Membranes were rapidly thawed, resuspended in 50 mM Tris-HCl, pH 7.4, 150 mM NaCl (plus protease inhibitors) and membrane proteins were solubilized in 1% n-dodecylmaltoside (DDM) for one hour at 4 °C. Insoluble material was removed by ultracentrifugation. Solubilized DRD2 (or mutant) were enriched using anti-HA agarose beads (Pierce/Thermo Scientific) at 4°C. The beads were washed with TBS with 0.1% DDM wash buffer three times and HA-tagged receptor variants were eluted with HA peptide (1 mg/mL in TBS with 0.1% DDM). To measure apparent receptor half-life, 10 nM partially purified DRD2WT or DRD2 with 5 designed mutations and designed truncated ICL3 were pre-incubated at 37 °C for 0–75 minutes prior to the addition and short preincubation with reaction mixture containing 20 μM [35S]GTPγS and 20 μM dopamine in a buffer containing 50 mM Tris-HCl, pH 7.2, 100 mM NaCl, 4 mM MgCl2 and 1 mM dithiothreitol. [35S]GTPγS binding was initiated with the addition of 4 μM Gi2 (purified as described previously[48]) for 30 min at room temperature. Reaction was terminated by rapid filtration on nitrocellulose filterplates. Microscint 0™ was added to dried filterplates and counted on a Topcount™ (Perkin Elmer). To assess the apparent melting temperatures, partially pure receptor preparations (10 nM) were pre-incubated at varying temperatures (0 – 50oC) for 30 min. Functional receptor was determined by assessing radioligand binding using [3H]dopamine in a buffer containing 50 mM Tris-HCl, pH 7.2, 100 mM NaCl and 4 mM MgCl2. Reaction was terminated by rapid filtration on nitrocellulose filterplates. Microscint 0™ was added to dried filterplates and counted on a Topcount™ (Perkin Elmer). Apparent melting temperatures and active state half-life curves were fitted using Graphpad Prism version 8 (Graphpad Software Inc., San Diego CA).

Data Availability

Structural data have been deposited in the Protein Data Bank (PDB) with coordinate accession number 6vms, and maps have been deposited in the Electron Microscopy Data Bank (EMDB) with accession numbers EMD-21243, EMD-21244, and EMD-21245. All other data generated or analyzed during this study are included in this published article or are available from the corresponding authors on reasonable request. Rational design of receptor variants shifted in the active state. (a) Energy landscape of 3 receptor variants with various energy differences in the resting (ligand-free) state ΔGresting, in the allosteric coupling controlling the inverse agonist efficacy ΔGinverse agonist and the stabilization of the receptor inactive state upon ligand stimulus. (b) Active state occupancy as a function of ΔGresting and ΔGinverse agonist (the energy of inverse agonist binding), as determined by the Boltzmann law. (c) Sensitivity of the receptors to inverse agonist binding. (d-f). Rational design of DRD2 receptor stabilized in the active state. (d) DRD2 WT active state homology model. The 120 residue long intracellular loop 3 (ICL3) is represented as a dotted blue line. (e) DRD2 designed with residue microswitches shown in spheres stabilizing selectively the active state structure. (f) DRD2 further stabilized in the active state through a de novo designed 7-residue stable ICL3. Effect of T4L, modified ICL3, and scFv16 on hDRD2WT ligand binding properties. (a) Competition binding of dopamine for [3H]spiperone on membranes containing human DRD2 wild type (hDRD2WT, blue) or human DRD2 modified with N-terminal T4L and de novo designed ICL3, nT4L-hDRD2(ΔICL3) (magenta). Log Ki ± S.E.M. = −4.6 ± 0.04 for hDRD2WT. Log Ki ± S.E. = −4.6 ± 0.04 for nT4L-hDRD2(ΔICL3). n = 4 biologically independent experiments on Sf9 cell membranes. (b) Agonist-stimulated [35S]GTPγS binding of nT4L-hDRD2(ΔICL3) (magenta symbols) compared to hDRD2WT (blue symbols). Dose-response using dopamine (closed magenta) or bromocriptine (open magenta) show log EC50 ± S.E.M. = −7.5 ± 0.02 M and −7.8 ± 0.06 M, respectively, for nT4L-hDRD2(ΔICL3) compared to log EC50 ± S.E.M. −7.0 ± 0.02 M for dopamine on hDRD2WT (open blue). n = 3 biologically independent experiments on HEK293T cell membranes. c) Saturation binding of [3H]spiperone to membranes containing DRD2WT with or without scFv16. Kd ± S.E.M. = 0.15 nM ± 0.03 nM for control. Kd ± S.E.M. = 0.14 nM ± 0.02 nM with scFv16. n = 4 biologically independent experiments on Sf9 cell membranes. (d) Competition binding of dopamine for [3H]spiperone on membranes containing DRD2WT with or without scFv16. Log Khigh ± S.E.M. = −6.4 ± 0.3 and Log Klow ± S.E.M. = −4.3 ± 0.08 for control. Log Khigh ± S.E.M. = −6.5 ± 0.2 and Log Klow ± S.E.M. = −4.3 ± 0.07 with scFv16. n = 4 biologically independent experiments on Sf9 cell membranes. In all panels, each data point is displayed as mean with error bars showing ± S.E.M., with individual data points for each repeat represented in small symbols. Function and thermostability of hDRD2EM. (a) Gi stimulation by dopamine in Sf9 cells transformed with Gi and hDRD2 baculoviruses. Log EC50 ± S.E.M. = −6.3 ± 0.1 for hDRD2WT (blue). Log EC50 ± S.E.M. = −6.6 ± 0.3 for hDRD2EM (orange). n = 4 biologically independent experiments on Sf9 cell membranes. (b) Inhibition of dopamine-stimulated (10 μM) Gi activation by the inverse agonist spiperone. Log IC50 ± S.E.M = −8.6 ± 0.1 for hDRD2WT (blue). Log IC50 ± S.E.M = −6.3 ± 0.2 for hDRD2EM (orange). n = 4 biologically independent experiments on Sf9 cell membranes. (c) Modulation of basal Gi activation by dopamine (open), bromocriptine (closed) or spiperone (squares) with the hDRD2EM (orange symbols) relative to hDRD2WT(blue symbols). Log EC50 ± S.E.M. for dopamine (−6.6 ± 0.2 M), bromocriptine (−7.6 ± 0.1 M) and spiperone (−7.6 ± 0.2 M) for hDRD2EM. Log EC50 ± S.E.M. for dopamine (−7.0 ± 0.1 M) and bromocriptine (−7.8 ± 0.2 M) for hDRD2WT. n = 3 biologically independent experiments on Sf9 cell membranes. (d) Apparent functional stability of agonist-bound hDRD2 constructs assessed by measuring the fraction of partially-purified receptor activating Gi as a function of incubation time at 37ºC. Half-life ± S.E.M.= 27.6 min ± 1.0 min for hDRD2WT (blue). Half-life ± S.E.M. = 49.3 min ± 1.1 min for hDRD2 with 5 designed mutations and designed truncated ICL3 (orange). n = 3 biologically independent experiments (e) Apparent thermostability of agonist-bound hDRD2 constructs assessed by measuring the fraction of partially-purified receptor binding agonist as a function of temperature. Tm ± S.E.M.= 30.9 °C ± 0.9 °C for hDRD2WT(blue). Tm ± S.E.M. = 39.6 °C ± 0.3 °C for hDRD2 with 5 designed mutations and designed truncated ICL3 (orange). n = 3 biologically independent experiments. In all panels, each data point is displayed as mean with error bars showing ± S.E.M., with individual data points for each repeat represented in small symbols. Purification and reconstitution of the hDRD2EM-Gi complex. (a) Superose 6 gel filtration profile of hDRD2EM-Gi-scFv16-HDL purified by M1 FLAG affinity chromatography. (b) Coomassie stained PAGE of the isolated peak fraction from gel filtration. (c) Overlaid Superose 200 gel filtration profiles (in LMNG detergent) of different hDRD2 constructs purified by M1 FLAG chromatography with saturating bromocriptine present. nT4L-hDRD2–5mut has the 5 thermostabilizing mutations of the EM construct, but has full-length ICL3. nT4L-hDRD2-ΔICL3 has the wild-type TM region (no mutations) but the truncated ICL3. Receptors were purified as in Methods, but in the absence of co-transduced G protein baculovirus. Cryo-EM data analysis. (a) Representative 2D class averages. (b) Gold-standard FSC curves of the 3D reconstructions. (c) Image processing procedure with the final maps colored based on local resolution. 3D classification with signal subtraction, model validation and representative map density. (a) Focused 3D classification subtracting T4L, AH domain, scFv16, and nanodisc density. (b) Focused 3D classification subtracting all but the G protein. (c) FSC curve showing the model-map correlation (focused map without T4L, AH domain, scFv16, and nanodisc density). (d) representative DRD2 TM and ligand density. Surfaces of the orthosteric pocket in the active and inactive conformations of DRD2. (a) Solvent accessible surface for the active conformation (purple) bound to bromocriptine (yellow sticks). (b) Solvent accessible surface for the inactive conformation (orange) bound to risperidone (green sticks), showing the deeper subpocket (PDB 6cm4). Interaction stabilizing the designed hDRD2EM active state conformation. (a) Zoomed view of the TM6-TM7 interface on the cytoplasmic side of the receptor. Atomic interactions involving the designed Y6.40 were accurately predicted in the design model. The conformation of the R3.50 is displaced in the experimental structure due to the binding of Giα which was absent in the model. Experimental structure is in purple and predicted DRD2 conformation is in green. (b, c) Conformational energy landscape of the WT (b) and designed (c) DRD2 receptors constructed from multiple all-atom MD simulation trajectories of the receptor active state homology models (orange dot) in absence of bound agonist. The conformational space is reported along 2 canonical active state structure metrics (TM3–6 distance and NPxxY rmsd to inactive state). The DRD2 WT relaxes primarily to an inactive state conformation while the designed variant remains in partially active state conformation displaying large TM3–6 distance. The antagonist-bound DRD2 crystal structure (PDB 6CM4) is used as a reference for the inactive state (green square). AH domain density and comparison with rhodopsin-Gi. (a) AH domain density in the current structure. (b) AH domain density in the rhodopsin-Gi complex (PDB 6cmo). (c) Alignment of Gαi subunit in hDRD2EM-Gi (purple) and rhodopsin-Gi complex (blue). (d) Alignment of AH domain in hDRD2EM-Gi (purple) and rhodopsin-Gi complex (blue). Electrostatic surfaces in the hDRD2EM-Gi complex. (a) Surface potential of the DRD2 and Gi components shown separately. (b) Surface potential of the complex, with scale in eV. Calculation was done using the APBS plugin in Pymol. Cryo-EM data collection, refinement and validation statistics.
  64 in total

Review 1.  Molecular biology of dopamine receptors.

Authors:  D R Sibley; F J Monsma
Journal:  Trends Pharmacol Sci       Date:  1992-02       Impact factor: 14.819

Review 2.  Advances and challenges in the search for D2 and D3 dopamine receptor-selective compounds.

Authors:  Amy E Moritz; R Benjamin Free; David R Sibley
Journal:  Cell Signal       Date:  2017-07-14       Impact factor: 4.315

3.  Structure of the human dopamine D3 receptor in complex with a D2/D3 selective antagonist.

Authors:  Ellen Y T Chien; Wei Liu; Qiang Zhao; Vsevolod Katritch; Gye Won Han; Michael A Hanson; Lei Shi; Amy Hauck Newman; Jonathan A Javitch; Vadim Cherezov; Raymond C Stevens
Journal:  Science       Date:  2010-11-19       Impact factor: 47.728

4.  History of the discovery of the antipsychotic dopamine D2 receptor: a basis for the dopamine hypothesis of schizophrenia.

Authors:  Bertha K Madras
Journal:  J Hist Neurosci       Date:  2013       Impact factor: 0.529

Review 5.  Molecular diversity of the dopamine receptors.

Authors:  O Civelli; J R Bunzow; D K Grandy
Journal:  Annu Rev Pharmacol Toxicol       Date:  1993       Impact factor: 13.820

Review 6.  Targeting the dopamine D2 receptor in schizophrenia.

Authors:  Philip Seeman
Journal:  Expert Opin Ther Targets       Date:  2006-08       Impact factor: 6.902

Review 7.  ANNIVERSARY REVIEW: 50 years since the discovery of bromocriptine.

Authors:  G Michael Besser; Ronald F Pfeiffer; Michael O Thorner
Journal:  Eur J Endocrinol       Date:  2018-05-11       Impact factor: 6.664

Review 8.  Neuropathological assessment of Parkinson's disease: refining the diagnostic criteria.

Authors:  Dennis W Dickson; Heiko Braak; John E Duda; Charles Duyckaerts; Thomas Gasser; Glenda M Halliday; John Hardy; James B Leverenz; Kelly Del Tredici; Zbigniew K Wszolek; Irene Litvan
Journal:  Lancet Neurol       Date:  2009-12       Impact factor: 44.182

9.  Structure of the D2 dopamine receptor bound to the atypical antipsychotic drug risperidone.

Authors:  Sheng Wang; Tao Che; Anat Levit; Brian K Shoichet; Daniel Wacker; Bryan L Roth
Journal:  Nature       Date:  2018-01-24       Impact factor: 49.962

10.  Repurposing of Bromocriptine for Cancer Therapy.

Authors:  Ean-Jeong Seo; Yoshikazu Sugimoto; Henry Johannes Greten; Thomas Efferth
Journal:  Front Pharmacol       Date:  2018-10-08       Impact factor: 5.810

View more
  41 in total

1.  Structural insights into the human D1 and D2 dopamine receptor signaling complexes.

Authors:  Youwen Zhuang; Peiyu Xu; Chunyou Mao; Lei Wang; Brian Krumm; X Edward Zhou; Sijie Huang; Heng Liu; Xi Cheng; Xi-Ping Huang; Dan-Dan Shen; Tinghai Xu; Yong-Feng Liu; Yue Wang; Jia Guo; Yi Jiang; Hualiang Jiang; Karsten Melcher; Bryan L Roth; Yan Zhang; Cheng Zhang; H Eric Xu
Journal:  Cell       Date:  2021-02-10       Impact factor: 41.582

2.  Structure-based design of a novel third-generation antipsychotic drug lead with potential antidepressant properties.

Authors:  Zhangcheng Chen; Luyu Fan; Huan Wang; Jing Yu; Dengyu Lu; Jianzhong Qi; Fen Nie; Zhipu Luo; Zhen Liu; Jianjun Cheng; Sheng Wang
Journal:  Nat Neurosci       Date:  2021-12-09       Impact factor: 24.884

Review 3.  Characteristic structural difference between inactive and active states of orexin 2 receptor determined using molecular dynamics simulations.

Authors:  Shun Yokoi; Ayori Mitsutake
Journal:  Biophys Rev       Date:  2021-11-03

4.  Cryo-EM structure of the human MT1-Gi signaling complex.

Authors:  Hiroyuki H Okamoto; Hirotake Miyauchi; Asuka Inoue; Francesco Raimondi; Hirokazu Tsujimoto; Tsukasa Kusakizako; Wataru Shihoya; Keitaro Yamashita; Ryoji Suno; Norimichi Nomura; Takuya Kobayashi; So Iwata; Tomohiro Nishizawa; Osamu Nureki
Journal:  Nat Struct Mol Biol       Date:  2021-08-05       Impact factor: 15.369

5.  Ligand-directed bias of G protein signaling at the dopamine D2 receptor.

Authors:  Ee Von Moo; Kasper Harpsøe; Alexander S Hauser; Ikuo Masuho; Hans Bräuner-Osborne; David E Gloriam; Kirill A Martemyanov
Journal:  Cell Chem Biol       Date:  2021-07-23       Impact factor: 8.116

6.  A target-agnostic screen identifies approved drugs to stabilize the endoplasmic reticulum-resident proteome.

Authors:  Mark J Henderson; Kathleen A Trychta; Shyh-Ming Yang; Susanne Bäck; Adam Yasgar; Emily S Wires; Carina Danchik; Xiaokang Yan; Hideaki Yano; Lei Shi; Kuo-Jen Wu; Amy Q Wang; Dingyin Tao; Gergely Zahoránszky-Kőhalmi; Xin Hu; Xin Xu; David Maloney; Alexey V Zakharov; Ganesha Rai; Fumihiko Urano; Mikko Airavaara; Oksana Gavrilova; Ajit Jadhav; Yun Wang; Anton Simeonov; Brandon K Harvey
Journal:  Cell Rep       Date:  2021-04-27       Impact factor: 9.423

7.  Crystal structure of dopamine D1 receptor in complex with G protein and a non-catechol agonist.

Authors:  Bingfa Sun; Dan Feng; Matthew Ling-Hon Chu; Inbar Fish; Silvia Lovera; Zara A Sands; Sebastian Kelm; Anne Valade; Martyn Wood; Tom Ceska; Tong Sun Kobilka; Florence Lebon; Brian K Kobilka
Journal:  Nat Commun       Date:  2021-06-03       Impact factor: 14.919

8.  Adaptive Membrane Fluidity Modulation: A Feedback Regulated Homeostatic System and Target for Pharmacological Intervention.

Authors:  Elzbieta Izbicka; Robert T Streeper; Christopher Louden
Journal:  In Vivo       Date:  2021 Nov-Dec       Impact factor: 2.155

9.  Evaluation of Substituted N-Phenylpiperazine Analogs as D3 vs. D2 Dopamine Receptor Subtype Selective Ligands.

Authors:  Boeun Lee; Michelle Taylor; Suzy A Griffin; Tamara McInnis; Nathalie Sumien; Robert H Mach; Robert R Luedtke
Journal:  Molecules       Date:  2021-05-26       Impact factor: 4.927

10.  Allosteric modulation of ghrelin receptor signaling by lipids.

Authors:  Marjorie Damian; Maxime Louet; Antoniel Augusto Severo Gomes; Céline M'Kadmi; Séverine Denoyelle; Sonia Cantel; Sophie Mary; Paulo M Bisch; Jean-Alain Fehrentz; Laurent J Catoire; Nicolas Floquet; Jean-Louis Banères
Journal:  Nat Commun       Date:  2021-06-24       Impact factor: 14.919

View more

北京卡尤迪生物科技股份有限公司 © 2022-2023.