Common activation mechanism of class A GPCRs
Abstract
Class A G-protein-coupled receptors (GPCRs) influence virtually every aspect of human physiology. Understanding receptor activation mechanism is critical for discovering novel therapeutics since about one-third of all marketed drugs target members of this family. GPCR activation is an allosteric process that couples agonist binding to G-protein recruitment, with the hallmark outward movement of transmembrane helix 6 (TM6). However, what leads to TM6 movement and the key residue level changes of this movement remain less well understood. Here, we report a framework to quantify conformational changes. By analyzing the conformational changes in 234 structures from 45 class A GPCRs, we discovered a common GPCR activation pathway comprising of 34 residue pairs and 35 residues. The pathway unifies previous findings into a common activation mechanism and strings together the scattered key motifs such as CWxP, DRY, Na+ pocket, NPxxY and PIF, thereby directly linking the bottom of ligand-binding pocket with G-protein coupling region. Site-directed mutagenesis experiments support this proposition and reveal that rational mutations of residues in this pathway can be used to obtain receptors that are constitutively active or inactive. The common activation pathway provides the mechanistic interpretation of constitutively activating, inactivating and disease mutations. As a module responsible for activation, the common pathway allows for decoupling of the evolution of the ligand binding site and G-protein-binding region. Such an architecture might have facilitated GPCRs to emerge as a highly successful family of proteins for signal transduction in nature.
Introduction
As the largest and most diverse group of membrane receptors in eukaryotes, GPCRs mediate a wide variety of physiological functions (Lagerström and Schiöth, 2008; Rosenbaum et al., 2009; Katritch et al., 2012; Venkatakrishnan et al., 2013; Katritch et al., 2013), including vision, olfaction, taste, neurotransmission, endocrine and immune responses via more than 800 family members, and are involved in many diseases (Rana et al., 2001; Smit et al., 2007; Vassart and Costagliola, 2011; Thompson et al., 2014; Hauser et al., 2018). Therefore, GPCRs are important drug targets. There are 475 marketed drugs (~34% of all FDA-approved therapeutic agent agents) targeting 108 members of the GPCR superfamily (Hauser et al., 2018; Hauser et al., 2017; Allen and Roth, 2011). Class A is the largest and most diverse GPCR subfamily in humans (Kolakowski, 1994; Bockaert and Pin, 1999; Fredriksson et al., 2003; Isberg et al., 2016), including 388 olfactory (Krautwurst et al., 1998; Spehr and Munger, 2009) and 286 non-olfactory receptors (Pándy-Szekeres et al., 2018; Munk et al., 2019) (Figure 1a). They share a seven-transmembrane (7TM) helices domain, with ligand binding pocket and G-protein-binding region located in the extracellular and intracellular ends of the helix bundle. Responding to a wide variety of extracellular signals ranged from small molecules to peptides even proteins, the extracellular facing ligand-binding pockets have evolved to be highly diverse in both shape and sequences (Venkatakrishnan et al., 2013; Ngo et al., 2017; Vass et al., 2018). Similarly, the G-protein-binding regions are also quite diverse in sequences, modulating the activity of different signalling pathways by recruiting dozens of heterotrimeric G proteins (Rasmussen et al., 2011a; Du et al., 2019; Liu et al., 2019a), arrestins (Gainetdinov et al., 2004; Zhou et al., 2017; Yang et al., 2018; Latorraca et al., 2018), GPCR kinases (Komolov et al., 2017) in a ligand-specific manner. Residues that connect the ligand-binding pocket to the G-protein-coupling region are significantly more conserved (Ballesteros and Weinstein, 1995; Isberg et al., 2015), with evolutionarily conserved sequence motifs (CWxP [Eddy et al., 2018; Filipek, 2019; Wescott et al., 2016; Tehan et al., 2014; Holst et al., 2010; Nygaard et al., 2009; Hofmann et al., 2009; Trzaskowski et al., 2012], PIF [Ballesteros and Weinstein, 1995; Ishchenko et al., 2017; Schönegge et al., 2017; Kato et al., 2019] Na+ pocket [Eddy et al., 2018; Filipek, 2019; Liu et al., 2012; Yuan et al., 2013; Fenalti et al., 2014; Katritch et al., 2014; Vickery et al., 2018; White et al., 2018; Ye et al., 2018; Chen et al., 2019], NPxxY [Rasmussen et al., 2011a; Filipek, 2019; Wescott et al., 2016; Nygaard et al., 2009; Hofmann et al., 2009; Trzaskowski et al., 2012; Schönegge et al., 2017; Chen et al., 2019; Venkatakrishnan et al., 2016] and DRY [Schönegge et al., 2017; Alhadeff et al., 2018; Jacobson et al., 2014; Feng et al., 2017; Roth et al., 2017; Shihoya et al., 2017; Yuan et al., 2014]) scattered in the intracellular half of the 7TM domain (Figure 1—figure supplement 1).
GPCR activation is agonist binding induced G-protein recruitment (Eddy et al., 2018; Christopoulos and Kenakin, 2002; Manglik et al., 2015; Staus et al., 2016; DeVree et al., 2016). It is an allosteric process (Gilchrist, 2007; May et al., 2007; Christopoulos, 2014; Latorraca et al., 2017), transducing various external stimuli into cellular responses. Understanding the activation mechanism of GPCR is of paramount importance in pharmacology research and drug discovery. Tremendous previous efforts, involving sequence analysis, structural, biophysical, biochemical and computational approaches such as X-ray (Popov et al., 2018; Roth et al., 2008; Warne et al., 2009; Carpenter et al., 2016; Nehmé et al., 2017; Tsai et al., 2018; Rasmussen et al., 2011b; Pardon et al., 2014; Rosenbaum et al., 2007; Kobilka and Schertler, 2008; Cherezov et al., 2007; Ghosh et al., 2015; Cherezov et al., 2004; Caffrey and Cherezov, 2009; Liu et al., 2013; Weierstall et al., 2014; Stauch and Cherezov, 2018; Zhang et al., 2015), NMR (Ye et al., 2018; Chen et al., 2019; Nygaard et al., 2013; Lamichhane et al., 2015; Isogai et al., 2016; Sounier et al., 2015; Ye et al., 2016; Shimada et al., 2019), Cryo-EM (Liang et al., 2017; Zhang et al., 2017; Cheng, 2018; Renaud et al., 2018), labeling biosensors (Irannejad et al., 2013; Tian et al., 2017), FRET (Gregorio et al., 2017; Halls and Canals, 2018; Sandhu et al., 2019), BRET (Lan et al., 2012; Lee et al., 2016; Okashah et al., 2019), DEER (Wingler et al., 2019; Van Eps et al., 2018; Dror et al., 2015), molecular dynamic simulations (Yuan et al., 2013; Dror et al., 2015; Dror et al., 2011a; Dror et al., 2011b; Dror et al., 2013; Miao et al., 2013; Bhattacharya and Vaidehi, 2014; Kohlhoff et al., 2014; Ciancetta et al., 2015; Alhadeff et al., 2018), evolutionary tracing (Madabushi et al., 2004; Schöneberg et al., 2007; Rodriguez et al., 2010), molecular docking (Kufareva et al., 2011; Jacobson et al., 2014; Kooistra et al., 2016; Miao et al., 2016; Feng et al., 2017; Roth et al., 2017; Lyu et al., 2019; Cooke et al., 2015) and mutagenesis (Schönegge et al., 2017; Sung et al., 2016; Massink et al., 2015; Ragnarsson et al., 2019; Hulme, 2013) have been made to study the allosteric nature of GPCRs including but not limited to receptor activation (Wescott et al., 2016; Tehan et al., 2014; Trzaskowski et al., 2012; Venkatakrishnan et al., 2016; Isom and Dohlman, 2015; Okada et al., 2001; Hunyady et al., 2003; Dalton et al., 2015; Lans et al., 2015), G protein activation (Trzaskowski et al., 2012; Flock et al., 2017; Flock et al., 2015; Furness et al., 2016; Glukhova et al., 2018; Ilyaskina et al., 2018; Inoue et al., 2019; Weis and Kobilka, 2018), biased agonism (McCorvy et al., 2018; Onaran et al., 2014; Schmid et al., 2017; Smith et al., 2018; Tan et al., 2018; Whalen et al., 2011; Wootten et al., 2018; Wootten et al., 2016; Masureel et al., 2018), ligand efficiency (Gregorio et al., 2017; Furness et al., 2016; Livingston et al., 2018; Solt et al., 2017; Yao et al., 2009), allosteric modulators (Thal et al., 2018; Kruse et al., 2013; Leach et al., 2007; Liu et al., 2019b; Lu and Zhang, 2019; Robertson et al., 2018; Shao et al., 2019; Wu et al., 2019; Zheng et al., 2016; Ortiz Zacarías et al., 2018; Jaeger et al., 2019; Hollingsworth et al., 2019; Liu et al., 2017; Chaturvedi et al., 2018; Oswald et al., 2016), and inverse agonism (Chaturvedi et al., 2018; Oswald et al., 2016; Hori et al., 2018; Nagiri et al., 2019; Peng et al., 2018; Shihoya et al., 2017). Starting from the first structure of a GPCR-G protein complex (β2AR-Gs) (Rasmussen et al., 2011a), the rapidly growing structures of receptor-G-protein complex have provided excellent opportunity to better understand receptor conformation changes upon activation. Meanwhile, mutagenesis studies on different receptors also identified functional roles of key residues in receptor activation, one good example is CXCR4 (Wescott et al., 2016), where a large-scale mutagenesis study covering all 352 residues of the receptor identified 41 amino acids that are required for signalling induced by agonist CXCL12. Notably, family-wide analysis on GPCR activation with the concept of residue contacts (Venkatakrishnan et al., 2013; Flock et al., 2017; Kayikci et al., 2018) have revealed the converged activation pathway near the G-protein-coupling region (Venkatakrishnan et al., 2016) and selectivity determinants of GPCR–G-protein binding (Flock et al., 2017). While these studies have provided key insights into GPCR activation mechanism for individual receptors or specific motifs, a family-wide common activation mechanism that directly connect ligand-binding pocket and G-protein-coupling region has yet to be discovered. Although it is well established that outward movement of transmembrane helix 6 (TM6) upon ligand binding is a common feature of receptor activation, the residue level changes that trigger the movement of TM6 remain less well understood (Figure 1b).
Receptor activation requires global reorganization of residue contacts as well as water-mediated interactions (Yuan et al., 2014; Yuan et al., 2016; Venkatakrishnan et al., 2019). Since prior studies primarily investigated conformational changes through visual inspection (Tehan et al., 2014; Trzaskowski et al., 2012) or through the presence or absence of non-covalent contacts between residues (Venkatakrishnan et al., 2016; Flock et al., 2017), we reasoned that one could gain comprehensive knowledge about mechanism of receptor activation by developing approaches that can capture not just the presence or absence of a contact but also subtle, and potentially important alterations in conformations upon receptor activation.
Results
A residue-residue contact score-based framework to characterize GPCR conformational changes
To address this, we developed an approach to rigorously quantify residue contacts in proteins structures and infer statistically significant conformational changes. We first defined a residue-residue contact score (RRCS) which is an atomic distance-based calculation that quantifies the strength of contact between residue pairs (Ngo et al., 2017) by summing up all possible inter-residue heavy atom pairs (Figure 2a and Figure 2—figure supplement 1a). We then defined ∆RRCS, which is the difference in RRCS of a residue pair between any two conformational states of a receptor that quantitatively describes the rearrangements of residue contacts (Figure 2b and Figure 2—figure supplement 1b). While RRCS can be 0 (no contact) or higher (stronger contact), ∆RRCS can be negative (loss in strength of residue contact), positive (gain in strength of residue contact) or 0 (no change in strength of residue contact). To capture the entirety of conformational changes in receptor structure upon activation, we computed the ∆RRCS between the active and inactive states of a receptor and defined two types of conformational changes (Figure 2c): (i) switching contacts: these are contacts that are present in the inactive state but lost in the active state (or vicw versa) such as loss of intra-helical contacts between D/E3×49 (GPCRdb numbering [Isberg et al., 2016]) and R3×50, and gain of inter-helical hydrophobic contacts between residues at 3×40 and 6×48 upon receptor activation; and (ii) repacking contacts: these are contacts that result in an increase or decrease in residue packing such as the decreased packing of intra-helical side-chain contacts between W6×48 and F6×44, and the increase in inter-helical residue packing due to the translocation of N7×49 toward D2×50 upon receptor activation. In this manner, we quantified the global, local, major and subtle conformational changes in a systematic way (i.e. inter-helical and intra-helical, switching and repacking contacts).
We then analyzed 234 structures of 45 class A GPCRs that were grouped into three categories (Figure 1a): (i) antagonist- or inverse agonist-bound (inactive; 142 structures from 38 receptors); (ii) both agonist- and G protein/G protein mimetic-bound (fully active; 27 structures from eight receptors); and (iii) agonist-bound (intermediate; 65 structures from 15 receptors). Among them, six receptors [rhodopsin (bRho) (Li et al., 2004; Choe et al., 2011), β2-adrenergic receptor (β2AR) (Rasmussen et al., 2011a; Cherezov et al., 2007), M2 muscarinic receptor (M2R) (Kruse et al., 2013; Haga et al., 2012), μ-opioid receptor (μOR) (Manglik et al., 2012; Huang et al., 2015), adenosine A2A receptor (A2AR) (Carpenter et al., 2016; Jaakola et al., 2008) and κ-opioid receptor (κ-OR) (Wu et al., 2012; Che et al., 2018)] have both inactive- and active-state crystal structures available. Given that ∆RRCS can capture major and subtle conformational changes, we computed RRCS for all structures and ∆RRCS for the six pairs of receptors and investigated the existence of a common activation pathway (i.e. a common set of residue contact changes) across class A GPCRs. Two criteria (Figure 2d; further details in Materials and methods) were applied to identify conserved rearrangements of residue contacts: (i) equivalent residue pairs show a similar and substantial change in RRCS between the active and inactive state structures of each of the six receptors (i.e. the same sign of ∆RRCS and |∆RRCS| > cut-off for all receptors) and (ii) family-wide comparison of the RRCS for the 142 inactive and 27 active state structures shows a statistically significant difference (p<0.001; two sample t-test). This allowed us to reliably capture both the major rearrangements as well as subtle but conserved conformational changes at the level of individual residues in diverse GPCRs in a statistically robust and significant manner. Consistent with this, a comparison with earlier studies revealed that the RRCS based approach is able to capture a larger number of conserved large-scale and subtle changes in residues contacts (Figure 2d) that would have been missed by visual inspection or residue contact presence/absence criteria alone (see Materials and methods for conceptual advance of this approach and detailed comparison).
Discovery of a common and conserved receptor activation pathway
Remarkably, for the first time, our analysis of the structures allowed the discovery of a common and conserved activation pathway that directly links ligand-binding pocket and G protein-coupling regions in class A GPCRs (Figure 3). The pathway is comprised of 34 residue pairs (formed by 35 residues) with conserved rearrangement of residue contacts upon activation (Figure 2d), connecting several well-known but structurally and spatially disconnected motifs (CWxP, PIF, Na+ pocket, NPxxY and DRY) all the way from the extracellular side (where the ligand binds) to the intracellular side (where the G-protein binds). Inspection of the rewired contacts as a ∆RRCS network reveals that the conserved receptor activation pathway is of modular nature and involves conformational changes in four layers. In layer 1, there is a conserved signal initiation step involving changes in residue contacts at the bottom of the ligand-binding pocket and Na+ pocket. In layer 2, critical hydrophobic contacts are broken (i.e., opening of the hydrophobic lock). In layer 3, microswitch residues (6×37, Y7×53) are rewired and in layer 4, the residue R3×50 and G protein contacting positions are rewired, making them competent to bind to G protein on the cytosolic side (Figure 3). Strikingly, recently released cryo-EM structures of four receptors (5-HT1B, rhodopsin, A1R and µOR) in complex with Gi/o(Glukhova et al., 2018; García-Nafría et al., 2018; Kang et al., 2018; Koehl et al., 2018; Draper-Joyce et al., 2018; Tsai et al., 2018) also support the conservation of contacts involving these 34 residue pairs (Figure 4, Figure 4—figure supplements 1 and 2). These observations highlight the conserved and common nature of a previously undescribed activation pathway linking ligand binding to G-protein coupling, regardless of the subtypes of intracellular effectors (i.e., Gs (Rasmussen et al., 2011a; Carpenter et al., 2016), Gi/o(Glukhova et al., 2018; García-Nafría et al., 2018; Kang et al., 2018; Koehl et al., 2018; Draper-Joyce et al., 2018; Tsai et al., 2018), arrestin (Zhou et al., 2017; Kang et al., 2015) or G-protein mimetic nanobody/peptide (Rasmussen et al., 2011b; Kruse et al., 2013; Choe et al., 2011; Huang et al., 2015; Che et al., 2018), Figure 4a).
Collectively, these findings illustrate how a combination of intra-helical and inter-helical switching contacts as well as repacking contacts underlies the common activation mechanism of GPCRs.
Molecular insights into key steps of the common activation pathway
Receptor activation is triggered by ligand binding and is characterised by movements of different transmembrane helices. How does ligand-induced receptor activation connect the different and highly conserved motifs, rewire residue contacts and result in the observed changes in transmembrane helices? To this end, we analyzed the common activation pathway in detail and mapped, where possible, how they influence helix packing, rotation and movement (Figure 3). A qualitative analysis suggests the presence of four layers of residues in the pathway linking the ligand binding residues to the G-protein-binding region.
Layer 1: We did not see a single ligand-residue contact that exhibits conserved rearrangement, which accurately reflects the diverse repertoire of ligands that bind GPCRs (Katritch et al., 2012; Venkatakrishnan et al., 2013; Ngo et al., 2017) (Figure 3—figure supplement 1). Instead, as a first common step, extracellular binding of diverse agonists converges to trigger an identical alteration of the transmission switch (3×40, 5×51, 6×44 and 6×48) and Na+ pocket (2×50, 3×39, 7×45 and 7×49). Specifically, the repacking of an intra-helical contact between residues at 6×48 and 6×44, together with the switching contacts of residue at 3×40 toward 6×48 and residue at 5×51 toward 6×44, contract the TM3-5-6 interface in this layer. This reorganization initializes the rotation of the cytoplasmic end of TM6. The collapse of Na+ pocket leads to a denser repacking of the four residues (2×50, 3×39, 7×45 and 7×49), initiating the movement of TM7 toward TM3. Remarkably, a recent NMR study on A2AR (Eddy et al., 2018) demonstrated the strong coupling between allosteric switch D2×50 and toggle switch W6×48, which is consistent with the present observation.
Layer 2: In parallel with these movements, two residues (6×40 and 6×41) switch their contacts with residue at 3×43, and form new contacts instead with residues at 5×58 and 5×55, respectively. Residues at 3×43, 6×40 and 6×41 are mainly composed of hydrophobic amino acids and referred as hydrophobic lock (Wescott et al., 2016; Tehan et al., 2014; Han et al., 2012; Martí-Solano et al., 2014). Its opening loosens the packing of TM3-TM6 and facilitates the outward movement of the cytoplasmic end of TM6, which is necessary for receptor activation. Additionally, N7×49 develops contacts with residue at 3×43 from nothing, facilitating the movement of TM7 toward TM3.
Layer 3: Upon receptor activation, Y7×53 loses its inter-helical contacts (Venkatakrishnan et al., 2016) with residues at 1×53 and 8×50, and forms new contacts with residues at 3×43, 3×46 and R3×50, which were closely packed with residues in TM6. Thus, the switching of contacts by Y7×53 strengthens the packing of TM3-TM7, while the packing of TM3-TM6 is further loosened with the outward movement of TM6.
Layer 4: Finally, the restrains on R3×50, including more conserved, local intra-helical contacts with D(E)3×49 and less conserved ionic lock with D(E)6×30, are eliminated and R3×50 is released. Notably, the switching contacts between R3×50 and residue at 6×37 are essential for the release of R3×50, which breaks the remaining contacts between TM3 and TM6 in the cytoplasmic end and drives the outward movement of TM6. The rewired contacts of R3×50 and other G-protein contacting positions (3×53, 3×54, 5×61 and 6×33) make the receptor competent to bind to G protein on the cytosolic side.
Together, these findings demonstrate that the intra-helical/inter-helical and switching/repacking contacts between residues is not only critical to reveal the continuous and modular nature of the activation pathway, but also to link residue-level changes to transmembrane helix-level changes in the receptor.
Common activation pathway induced changes in TM helix packing
To capture the patterns in the global movements of transmembrane helices, all inter-helical residues pairs in the common activation pathway were used to describe the inter-helical contacts between the cytoplasmic end of TM3 and TM6 as well as TM3 and TM7 (Figure 5a). Analysis of the RRCSTM3-TM7 (X-axis) and RRCSTM3-TM6 (Y-axis) for each of the 234 class A GPCR structures revealed distinct compact clusters of inactive and active states. Surprisingly, the inactive state has zero or close to zero RRCSTM3-TM7 regardless of the wide distribution of RRCSTM3-TM6. In contrast, the active state has a high RRCSTM3-TM7 and strictly zero RRCSTM3-TM6. Thus, receptor activation from inactive to active state occurs as a harmonious process of inter-helical contact changes: elimination of TM3-TM6 contacts, formation of TM3-TM7 contacts and repacking of TM5-TM6 (Figure 5b and Figure 5—figure supplement 1). In terms of global conformational changes, the binding of diverse agonists converges to trigger outward movement of the cytoplasmic end of TM6 and inward movement of TM7 toward TM3 (Rasmussen et al., 2011a; Nygaard et al., 2009; Venkatakrishnan et al., 2016), thereby creating an intracellular crevice for G protein coupling (Figure 5b).
It is noteworthy that the common activation pathway we discovered in this study is not the only pathway that connecting extracellular ligand-binding and intracellular effector coupling for class A GPCRs ‒ it is likely to be a shared portion of various activation pathways of GPCR members belonging to this class ‒ each receptor still has its unique receptor-, ligand- and effector-specific activation pathways. In fact, research on this subject has boosted the discovery of selective and biased ligands (McCorvy et al., 2018; Onaran et al., 2014; Schmid et al., 2017; Smith et al., 2018; Tan et al., 2018; Whalen et al., 2011; Wootten et al., 2018; Wootten et al., 2016).
As shown in Figure 5—figure supplement 2, collapse of the Na+ pocket leads to a denser repacking of six residues (five residue pairs), reflected by higher RRCSsodium_pocket scores in active state than that in inactive state structures. Recently, the crystal structures (5X33 [Hori et al., 2018], 6BQG [Peng et al., 2018] and 6K1Q [Nagiri et al., 2019]) whose ligand (inverse agonist) diffuses deep in ligand-binding pocket or even occupies the sodium binding pocket (below D2×50) were reported. These inverse agonists disrupt the collapse of Na+ pocket by blocking the rotation of W6×48 and/or taking the space of Na+, and stabilize the receptors in an inactive state. Indeed, these inactive state structures showed zero RRCSTM3-TM7 but high RRCSTM3-TM6 scores. The inverse agonism are not only consistent with both our activation model and mutagenesis experiments, but also supported by the NMR study of A2AR (Eddy et al., 2018). This study demonstrated the role of D522×50 as an allosteric link between the orthosteric ligand-binding site and the intracellular signaling surface, revealing strong interactions with the toggle switch W2466×48.
Experimental validation of the modular nature of the common activation pathway
Based on the knowledge of the common activation pathway, one would expect that mutations of residues in the pathway are likely to severely affect receptor activation. The two extreme consequences are constitutive activation (without agonist binding) or inactivation (abolished signalling). To experimentally test this hypothesis, we systematically designed site-directed mutagenesis for residues in the pathway on a prototypical receptor A2AR, aiming to create constitutively activating/inactivating mutations (CAM/CIM), by promoting/blocking residue and helix level conformational changes revealed in the pathway. 6/15 designed CAMs and 15/20 designed CIMs were validated by functional cAMP accumulation assays, and none of them were reported before for A2AR (Figure 6, Figure 6—figure supplement 1 and Figure 6—source data 1). The design of functional active/inactive mutants has been very challenging. However, the knowledge of common activation pathway of GPCRs presented here greatly improves the success rate. The mechanistic interpretation of 21 successful predicted mutants is explained below. We also discussed the 14 unsuccessful predictions in Figure 6—source data 2. Besides, we extended mutagenesis studies to Gs-coupled 5-HT7 and Gi-coupled 5-HT1B receptors by designing CAM/CIMs in residues at 3×40, 3×43, 6×40, 6×44, and 7×49 (Figure 6—figure supplement 2).
In layer 1, the mutation I923×40N likely stabilizes the active state by forming amide-π interactions with W2466×48 and hydrogen bond with the backbone of C1855×461, which rewires the packing at the transmission switch and initiates the outward movement of the cytoplasmic end of TM6; this mutation elevated the basal cAMP level by sevenfold. Conversely, I923×40A would reduce the favorable contacts with W6×48 upon activation, which retards the initiation of the outward movement of TM6; this mutation resulted in a decrease in both basal cAMP level [71% of wild-type (WT)] and agonist potency (eightfold). Another example is the residue at 6×44, the mutation F2426×44R would stabilize the inactive state by forming salt bridge with D522×50, which blocks the rotation of TM6 and thus abolishes Gs coupling; indeed this mutation greatly reduced basal cAMP level (to 63% WT) and agonist potency (by 374-fold). In contrast, F2426×44A would reduce contacts with W2466×48, loosen TM3-TM6 contacts, diminish the energy barrier of TM6 release and make outward movement of TM6 easier; consistently this mutation elevated the basal cAMP level (by twofold) and increased the agonist potency (by eightfold). Mutations of residues forming the Na+ pocket, such as D522×50A and N2807×45R, would destroy the hydrogen bond network at the Na+ pocket and retard the initiation of the inward movement of TM7. These mutations completely abolished agonist potency and greatly reduced the basal cAMP level (to 80% and 78% of WT, respectively).
In layer 2, the mutations L953×43A/R and I2386×40Y would loosen the hydrophobic lock, weaken TM3-TM6 contacts, promote the outward movement of cytoplasmic end of TM6 and eventually make receptor constitutively active; this is reflected by remarkably high basal cAMP production (28-, 2- and 11-fold increase, respectively). Notably, mutations at/near the Na+ pocket, L482×46R and N2847×49K, could lock the Na+ pocket at inactive packing mode by introducing salt bridge with D522×50, thus blocking the inward movement of TM7 toward TM3. As expected, these mutations completely abolished agonist potency. The CIMs at/near the Na+ pocket (from both layers 1 and 2) reflect that the subtle inward movement of TM7 towards TM3 is essential for receptor activation, which is often underappreciated and overshadowed by the movement of TM6. In line with this, two mutations on TM7, N2847×49A and Y2887×53A, attenuate the TM3-TM7 contacts upon activation and completely abolished or greatly reduced (by 16-fold) agonist potency, respectively.
In layer 3, I983×46A likely reduces contacts with Y2887×53, weakens the packing between TM3-TM7, and retards the movement of TM7 toward TM3; similarly, L2356×37A would reduce contacts with F2015×62, weaken the packing between TM5-TM6, and makes the TM6 movement toward TM5 more difficult. In line with the interpretation, these mutations resulted in reduced basal cAMP level (72% and 71% WT, respectively) and decreased agonist potency (23- and 4-fold, respectively). These results are consistent with previous findings on vasopressin type-2 receptor (V2R) (Venkatakrishnan et al., 2016).
In layer 4, D1013×49N likely diminishes its intra-helical interaction with R1023×50 and thus makes the release of the latter easier, which in turn promotes G-protein recruitment. Consistent with this possibility, this mutation led to a greatly elevated basal cAMP level (eightfold).
Despite these A2AR mutants greatly affect receptor activation, our radioligand binding assay shows that they generally retain the agonist binding ability, with the exception of two CIMs: W2466×48A and N2847×45K (Figure 6b,c and Figure 6—source data 1). This suggests that the common activation pathway is of modular nature and that such an organization allows for a significant number of residues involved in agonist binding to be uncoupled from receptor activation/inactivation and G-protein binding.
As shown in Figure 6—figure supplements 2 and 5-HT7 receptor mutations F3366×44R and N3807×49K completely abolished agonist potency and greatly reduced the basal cAMP level, which is remarkably consistent with the observation on A2AR, highlighting the crucial role of the highly conserved residues F6×44 and N7×49. Beyond Gs-coupled A2AR and 5-HT7 receptor, we also validated this mutation design in Gi-coupled 5-HT1B receptor. Indeed, two CIMs, I1373×40N and F3236×44H greatly reduced receptor-mediated Gi activity compared to WT, whereas three CAMs, L1733×43A in Gs-coupled 5-HT7 receptor, F3236×44A and I1373×40A in Gi-coupled 5-HT1B receptor, were verified to promote their basal activities, consistent with the observation on CAMs (L953×43A, F2426×44A and I923×40A) designed for A2AR.
The common pathway allows mechanistic interpretation of mutations
Four hundred thirty five disease-associated mutations were collected, among which 28% can be mapped to the common activation pathway, much higher than that to the ligand-binding and G-protein-binding regions (20% and 7%, respectively) (Figure 7a,b). Furthermore, 272 CAMs/CIMs from 41 receptors (Figure 7c) were mined from the literature for the 14 hub residues (i.e. residues that have more than one edges in the pathway).
The average number of disease-associated mutations in the common activation pathway is much higher than that of ligand-binding pocket, G-protein-binding site, and residues in other regions (2.5-, 3.5- and 3.5-fold, respectively), reflecting the enrichment of disease-associated mutations on the pathway (Figure 7a). Within this pathway, the enrichment of disease mutations and CAMs/CIMs in layers 1 and 2 is noteworthy, which highlights the importance of signal initiation and hydrophobic lock opening, and further supports the modular and hierarchical nature of GPCR activation (Figures 3 and 5b). Notably, for certain residues, such as D2×50 and Y7×53, only loss-of-function disease mutations or CIMs were observed (Figure 7), implying they are indispensable for receptor activation and the essential role of TM7 movement (Figures 3 and 5).
The functional consequence of these single point mutations can be rationalized by analysing if they are stabilizing/destabilizing the contacts in the common activation pathway or promoting/retarding the required helix movement upon activation (Figure 7b and Figure 7—figure supplement 1). For example, I1303×43N/F (layer 2) in V2R was reported as a gain-/loss-of-function mutation that causes nephrogenic syndrome of inappropriate antidiuresis (Erdélyi et al., 2015) or nephrogenic diabetes insipidus (Pasel, 2000), respectively. I1303×43N/F likely loosens/stabilizes the hydrophobic lock, weakens/strengthens the TM3-TM6 packing and leads to constitutively active/inactive receptors. Another example is T581×53M in rhodopsin, which was reported as a loss-of-function mutation that causes retinitis pigmentosa 4 (Napier et al., 2015). T581×53M likely increases hydrophobic contacts with Y3067×53 and P3037×50, which retards the inward movement of TM7 towards TM3 and eventually decreases G-protein recruitment. As in the case of disease-associated mutations, CAMs/CIMs that have been previously reported in the literature can also be interpreted by the framework of common activation pathway (Figure 7—figure supplement 1b). For example, F2486×44Y in CXCR4 (Wescott et al., 2016) was reported as a CIM. This residue likely forms hydrogen bond with S1233×39, which blocks the rotation of the cytoplasmic end of TM6, and decreases G-protein engagement.
Not surprisingly, the 35 residues constituting the pathway are highly conserved across class A GPCRs, dominated by physiochemically similar amino acids (Figure 7—figure supplement 2). The average sequence similarity of these positions across 286 non-olfactory class A receptors is 66.2%, significantly higher than that of ligand-binding pockets (31.9%) or signaling protein-coupling regions (35.1%). Together, these observations suggest that the modular and hierarchical nature of the activation pathway allows decoupling of the ligand-binding pocket, G-protein-binding site and the residues contributing to the common activation mechanism. Such an organization of the receptor might facilitate the uneven sequence conservation between different regions of GPCRs, confers their functional diversity in ligand recognition and G-protein binding while still retaining a common activation mechanism.
Discussion
Using a novel, quantitative residue contact descriptor, RRCS, and a family-wide comparison across 234 structures from 45 class A GPCRs, we reveal a common activation pathway that directly links ligand-binding pocket and G-protein-coupling region. Key residues that connect the different modules allows for the decoupling of a large number of residues in the ligand-binding site, G-protein contacting region and residues involved in the activation pathway. Such an organization may have facilitated the rapid expansion of GPCRs through duplication and divergence, allowing them to evolve independently and bind to diverse ligands due to removal of the constraint (i.e., between a large number of ligand-binding residues and those involved in receptor activation). This model unifies many previous important motifs and observations on GPCR activation in the literature (CWxP, PIF, Na+ pocket, NPxxY, DRY and hydrophobic lock) and is consistent with numerous experimental findings.
We focused on the common activation pathway (i.e. the common part of activation mechanism shared by all class A GPCRs and various intracellular effectors) in this study. Obviously, individual class A receptor naturally has its intrinsic activation mechanism(s), as a result of the diversified sequences, ligands and physiological functions. Indeed, receptor-specific activation pathways (including mechanisms of orthosteric, positive or negative allosteric modulators, biased signaling/selectivity of downstream effectors) have been revealed by both experimental studies including biophysical (such as X-ray, cryo-EM, NMR, FRET/BRET and DEER), biochemical and computational approaches (such as evolutionary trace analysis and molecular dynamics simulations), especially for the prototypical receptors such as rhodopsin, β2-adrenergic and A2A receptors. These studies demonstrated the complexity and plasticity of signal transduction of GPCRs. The computational framework we have developed may assist us in better understanding the mechanism of allosteric modulation, G-protein selectivity and diverse activation processes via intermediate states as more GPCR structures become available. While we interpret the changes as a linear set of events, future studies aiming at understanding dynamics could provide further insights into how the common activation mechanism operates in individual receptors.
Given the common nature of this pathway, we envision that the knowledge obtained from this study can not only be used to mechanistically interpret the effect of mutations in biological and pathophysiological context but also to rationally introduce mutations in other receptors by promoting/blocking residue and helix level movements that are essential for activation. Such protein engineering approaches may enable us to make receptors in specific conformational states to accelerate structure determination studies using X-ray crystallography or electron microscopy and functional investigation in the future. The method developed here could also be readily adapted to map allosteric pathways and reveal mechanisms of action for other key biological systems such as kinases, ion channels and transcription factors.
Materials and methods
Glossary
Request a detailed protocolTransmembrane domains (TMD): the core domain exists in all GPCRs, and consists of seven-transmembrane helices (TM1–7) that are linked by three extracellular loops (ECL1-3) and three intracellular loops (ICL1-3).
GPCRdb numbering scheme: a structure-based numbering system for GPCRs (Isberg et al., 2016; Isberg et al., 2015), an improved version of sequence-based Ballesteros–Weinstein numbering (Ballesteros and Weinstein, 1995) that considers structural distortions such as helical bulges or constrictions. The most conserved residue in a helix n is designated n×50, while other residues on the helix are numbered relative to this position.
Node: a point in a network at which lines intersect, branch or terminate. In this case, nodes represent amino acid residues.
Edge: a connection between the nodes in a network. In this case, an edge represents a residue-residue contact.
Hub: a node with two or more edges in a network.
Constitutively activating mutation (CAM): a mutant that could increase the inherent basal activity of the receptor by activating the G-protein-signaling cascade in the absence of agonist.
Constitutively inactivating mutation (CIM): a mutant completely abolishes receptor signalling.
GPCR structure data set
Request a detailed protocolAs of October 1, 2018, there are 234 released structures of 45 class A GPCRs with resolution better than 3.8 Å (Figure 1—source data 1), which covers 71% (203 out of 286 receptors, including 158 receptors that have no structures but share >50% sequence similarity in the TMD with the 45 structure-determined receptors) of class A GPCRs (Figure 1a). Based on the type of bound ligand and effector, these structures could be classified into three states: inactive state (antagonist or inverse agonist-bound, 142 structures from 38 receptors), active state (both agonist- and G protein/G protein mimetic-bound, 27 structures from eight receptors) and intermediate state (only agonist-bound, 65 structures from 15 receptors). In this study, we primarily focused on conformational comparison between inactive- and active- state structures, while also investigating the intermediate state structures. In the structure data set, seven receptors have both inactive and active structures: rhodopsin (bRho), β2-adrenergic receptor (β2AR), M2 muscarinic receptor (M2R), μ-opioid receptor (μOR), adenosine A2A receptor (A2AR), κ-opioid receptor (κ-OR) and adenosine A1A receptor (A1AR), the active state structure of which was recently determined by cryo-EM. In addition, 32 receptors have either inactive or active structures (Figure 1—source data 1).
Calculation of residue-residue contact score (RRCS)
View detailed protocolWe developed a much finer distance-based method (than coarse-grained Boolean descriptors such as contact map and residues contact [Venkatakrishnan et al., 2016; Kayikci et al., 2018; Eldridge et al., 1997; Verdonk et al., 2003; Wang et al., 2017; Adhikari and Cheng, 2016]), namely residue-residue contact score (RRCS). For a pair of residues, RRCS is calculated by summing up a plateau-linear-plateau form atomic contact score adopted from GPCR–CoINPocket (Ngo et al., 2017; Kufareva et al., 2011; Kufareva and Abagyan, 2012; Kufareva et al., 2014a; Kufareva et al., 2013; Kufareva et al., 2014b; Marsden and Abagyan, 2004) for each possible inter-residue heavy atom pairs (Figure 2—figure supplement 1a). GPCR–CoINPocket is a modified version of the hydrophobic term of ChemScore (Eldridge et al., 1997; Verdonk et al., 2003) that has been successfully used to describe hydrophobic contribution to binding free energy between ligand and protein. RRCS can describe the strength of residue-residue contact quantitatively in a much more accurate manner than Boolean descriptors (Venkatakrishnan et al., 2016; Flock et al., 2017). For example, Boolean descriptors do not capture side chain repacking if the backbone atoms of the two residues are close to each other (e.g. translocation of Y7×53 away from residue at 2×43 upon GPCR activation) and local contacts involving adjacent residues (residues within four/six amino acids in protein sequence) (e.g., disengagement between D/E3×49 and R3×50), while both cases can be well reflected by the change of RRCS (Figure 2c and Figure 2—figure supplement 1b).
All RRCS data can be found in Figure 2—source data 1. The computational details are described as below:
For the residue pairs between adjacent residues that are within four amino acids in protein sequence, only side chain heavy atom pairs were considered, atom pairs involving in backbone atoms (Cα, C, O, N) were excluded, since the latter seldom change during GPCR activation. For other residue pairs, all possible heavy atom pairs (including backbone atoms) were included when calculating RRCS.
Atomic contact scores are solely based on interatomic distance, and they are treated equally without weighting factors such as atom type or contact orientation. In principle, weighting of atomic contact by atom type and/or orientation would improve residue-residue contact score. However, parameterization of atom type or contact orientation is relatively arbitrary, subjective and complicated, especially considering the lipid bilayer environment surrounding GPCRs. Our preliminary study for 12 structures from six receptors (bRho, β2AR, M2R, µOR, A2AR and κ-OR) revealed that amino acids with hydrophobic side chains (one-letter code: A, V, I, L, M, P, F, Y, W) contribute to the majority (~88%) of residue pairs. Meanwhile, ionic lock opening of well-known motif DRY upon receptor activation can be adequately reflected by RRCS change between D/E3×49 and R3×50. These results suggest that interatomic distance-dependent residue pair contact score may represent an acceptable approximation of actual (either hydrophobic or charge-charge) interaction energies (Ngo et al., 2017) and is accurate enough for identifying conserved rearrangements of residue contacts upon receptor activation.
The quality of structures is extremely important for RRCS calculation. We adopted two criteria to exclude unreliable structures and residues: (a) crystal structures whose resolution is ≥3.8 Å. Structures in this category are: 5DGY (7.70 Å), 2I37 (4.20 Å), 2I36 (4.10 Å), 5TE5 (4.00 Å), 4GBR (4.00 Å), 5NJ6 (4.00 Å), 5V54 (3.90 Å), 2I35 (3.80 Å), 5D5B (3.80 Å), 4XT3 (3.80 Å); (b) residues whose residue-based real-space R-value (RSR [Jones et al., 1991]) is greater than 0.35. RSR is measure of how well ‘observed’ and calculated electron densities agree for a residue. RSR ranges from 0 (perfect match) to 1 (no match); RSR greater than 0.4 indicates a poor fit (Smart et al., 2018). Here we adopted a stricter cut-off, 0.35. Among the 234 class A GPCR structures, 156 have available RSR information (Kleywegt et al., 2004) (http://eds.bmc.uu.se), with 8.8% residues have RSR >0.35 and they are omitted in our analysis. For the 35 residues that constitute the common activation pathway, 255 out of 5460 RSR data points (~4.7%, lower than 8.8% for all residues) were omitted for having RSR values > 0.35.
For structures with multiple chains, RRCS were the average over all chains. For residues with multiple alternative conformations, RRCS was the sum of individual values multiplied by the weighting factor: occupancy value extracted from PDB files. Small molecule/peptide ligand, or intracellular binding partner (G protein or its mimetic) was treated as a single residue.
For the family-wide comparison of conformational changes upon activation, structurally equivalent residues are numbered by GPCRdb numbering scheme (Isberg et al., 2016; Isberg et al., 2015). Of the 35 residues in the common activation pathway, their GPCRdb numbering in all structures is almost identical to the Ballesteros–Weinstein numbering (Ballesteros and Weinstein, 1995), the exceptions are residues at 6×37, 6×41 and 6×44 for five receptors: FFAR1, P2Y1, P2Y12, F2R and PAR2, which are all from the delta branch of class A family.
Identification of conserved rearrangements of residue contacts upon activation
Request a detailed protocolUsing RRCS, structural information of TMD and helix eight in each structure can be decomposed into 400 ~ 500 residue pairs with positive RRCS. ∆RRCS, defined as RRCSactive − RRCSinactive, reflects the change of RRCS for a residue pair from inactive- to active- state (Figure 2b–d and Figure 2—figure supplement 1b). To identify residue pairs with conserved conformational rearrangements upon activation across class A GPCRs, two rounds of selections (Figure 2d and Figure 2—source data 1) were performed: (i) identification of conserved rearrangements of residue contacts upon activation for six receptors (bRho, β2AR, M2R, μOR, A2AR and κ-OR), that is equivalent residue pairs show a similar and substantial change in RRCS between the active and inactive state structure of each of the six receptors (the same sign of ∆RRCS and |∆RRCS| > 0.2 for all receptors) and (ii) family-wide RRCS comparison between the 142 inactive and 27 active state structures to identify residues pairs of statistically significant different (p<0.001; two sample t-test) RRCS upon activation.
Round 1. Identification of conserved rearrangements of residue contacts. Six receptors with available inactive- and active- state structures were analyzed using ∆RRCS to identify residue pairs that share similar conformational changes. Twelve representative crystal structures (high-resolution, no mutation or one mutation in TMD without affecting receptor signalling) were chosen in this stage: six inactive state structures (PDB codes 1GZM for bRho, 2RH1 for β2AR, 3UON for M2R, 4DKL for μOR, 3EML for A2AR and 4DJH for κ-OR) and six active state structures (3PQR for bRho, 3SN6 for β2AR, 4MQS for M2R, 5C1M for μOR, 5G53 for A2AR and 6B73 for κ-OR) (Figure 2d, Figure 2—figure supplement 1c and Figure 2—source data 1). Each receptor has approximately 600 residues pairs that have positive RRCS. Roughly one quarter are newly formed during receptor activation (RRCSinactive = 0 and RRCSactive >0); another quarter lose their contacts upon receptor activation (RRCSinactive >0 and RRCSactive = 0); and the remaining appear in both the inactive- or active- state structures (RRCSinactive >0 and RRCSactive >0), the contact rearrangement of which can only be reflected by ∆RRCS, but not Boolean descriptors.
To identify residue pairs that share conserved rearrangements of residue contacts upon activation, two steps are performed to qualify residue pairs for the next round. Firstly, residue pairs with same sign of ∆RRCS and |∆RRCS| > 0.2 for all six receptors were identified. There are 32 intra-receptor residues pairs (1×49:7×50, 1×53:7×53, 1×53:7×54, 2×37:2×40, 2×42:4×45, 2×43:7×53, 2×45:4×50, 2×46:2×50, 2×50:3×39, 2×57:7×42, 3×40:6×48, 3×43:6×40, 3×43:6×41, 3×43:7×49, 3×43:7×53, 3×46:6×37, 3×46:7×53, 3×49:3×50, 3×50:3×53, 3×50:6×37, 350:7×53, 3×51:5×57, 5×51:6×44, 5×58:6×40, 5×62:6×37, 6×40:7×49, 6×44:6×48, 7×50:7×55, 7×52:7×53, 7×53:8×50, 7×54:8×50 and 7×54:8×51) and five receptor-G protein/its mimetic residue pairs (3×50:G protein, 3×53:G protein, 3×54:G protein, 5×61:G protein and 6×33:G protein) that meet this criterion. Secondly, we also investigated residue pairs with ∆RRCS that are conserved in five receptors (i.e., with one receptor as exception). Considering there is no Na+ pocket for rhodopsin, three residue pairs (2×50:7×49, 6×44:7×45 and 6×48:7×45) around Na+ pocket were analyzed for five receptors but not bRho. Additionally, three residue pairs have 0 (3×46:3×50, 5×55:6×41) or negative (7×45:7×49) ∆RRCS for κ-OR but positive ∆RRCS for the other five receptors. As for 3×46:3×50, nanobody-stabilized active structures (β2AR: 3P0G, 4LDO, 4LDL, 4LDE, 4QKX; and μOR: 5C1M) generally have lower contact scores (<0.4) compared with G-protein-bound active-state structures (2.17 for 3SN6 of β2AR, 2.57 for 5G53 of A2AR and 6.93 for 3PQR of bRho). For these residue pairs, we added newly determined Gi-bound active A1AR and 5-HT1B receptor and found that they have positive ∆RRCS, like other five receptors (Figure 4—figure supplements 1 and 2). Thus, these three residue pairs (3×46:3×50, 5×55:6×41 and 7×45:7×49) were retained. Totally, six residue pairs with conserved ∆RRCS in five receptors were rescued. Taken together, 38 intra-receptor residue pairs and five receptor-G protein/its mimetic residue pairs were identified to have conserved rearrangements of residue contacts upon activation.
Round 2. Family-wide conservation analysis of residue contact pattern. To investigate the conservation of residue contact pattern for the 38 intra-receptor residue pairs across these functionally diverse receptors, two-tailed unpaired t-test between inactive state (142 inactive structures from 38 receptors) and active state (27 active structures from eight receptors) groups were performed (Figure 2d and Figure 2—source data 2). Thirty one residue pairs have significantly different RRCS between inactive- and active-state (p<10−5). As rhodopsin lacks the Na+ pocket, all rhodopsin structures were neglected in the analysis of 3 residue pairs around the pocket (2×50:7×49, 6×44:7×45 and 6×48:7×45), which have good p value (<10−3) for these non-rhodopsin class A GPCRs. Four residue pairs were filtered out in this round due to their poor p value, that is there are no statistically significant difference in RRCS between inactive and active states (p=0.01 for 2×37:2×40, 0.96 for 2×42:4×45, 0.02 for 2×45:4×50 and 0.014 for 2×57:7×42).
Finally, 34 intra-receptor residue pairs (Figure 2d, Figure 4—figure supplements 1 and 2) and five receptor-G-protein residue pairs were identified with conserved rearrangements of residue contacts upon activation, including all six residues pairs identified by the previous RC approaches (Venkatakrishnan et al., 2016).
Sequence analysis of class A GPCRs
Request a detailed protocolThe alignment of 286 non-olfactory, class A human GPCRs were obtained from the GPCRdb (Isberg et al., 2016; Isberg et al., 2015). The distribution of sequence similarity/identity across class A GPCRs were extracted from the sequence similarity/identity matrix for different structural regions using ‘Similarity Matrix’ tool in GPCRdb. The sequence conservation score (Figure 1—figure supplement 1) for all residue positions across 286 non-olfactory class A GPCRs were evaluated by the Protein Residue Conservation Prediction (Capra and Singh, 2007) tool with scoring method 'Property Entropy' (Mirny and Shakhnovich, 1999). Sequence conservation analysis (Figure 7—figure supplement 2) were visualized by WebLogo3 (Crooks et al., 2004) with sequence alignment files from GPCRdb as the input.
CAM/CIM in class A GPCRs
Request a detailed protocolFor the 14 hub residues in the common activation pathway, we collected the functional mutation data from the literature and GPCRdb (Isberg et al., 2016; Isberg et al., 2015). Mutations with ‘more than two fold-increase in basal activity/constitutively active’ or ‘abolished effect’ compared to the wild-type receptor were selected. Together, 272 mutations from 41 class A GPCRs on the 14 hub residues were collected, including the mutations we designed and validated in this study (Figure 7—source data 1).
Disease-associated mutations in class A GPCRs
Request a detailed protocolTo reveal the relationship between disease-associated mutations and related phenotypes of different transmembrane regions (Vassart and Costagliola, 2011; Thompson et al., 2014; Tao, 2006; Tao, 2008), we collected disease-associated mutation information for all 286 non-olfactory class A GPCRs by database integration and literature investigation. Four commonly used databases (UniProt [The UniProt Consortium, 2017], OMIM [Amberger et al., 2011], Ensembl [Zerbino et al., 2018] and GPCRdb [Isberg et al., 2016; Isberg et al., 2015]) were first filtered by disease mutations and then merged. Totally 435 disease mutations from 61 class A GPCRs were collected (Figure 1—source data 2).
Pathway-guided CAM/CIM design in A2AR
Request a detailed protocolWe designed mutations for a prototypical receptor A2AR, guided by the common activation pathway, aiming to get constitutively active/inactive receptor. Mutations that can either stabilize active or inactive state structures of A2AR or promote/block conformational changes upon activation were designed (Figure 6c and Figure 6—figure supplement 1) and tested by a functional cAMP accumulation assay. The inactive state structure 3EML and active state structure 5G53 were used. In silico mutagenesis was performed by Residue Scanning module in BioLuminate (Beard et al., 2013). Side-chain prediction with backbone sampling and a cut-off value of 6 Å were applied during the scanning. ΔStability is the change of receptor stability when introducing a mutation. We filtered the mutations by one of the following criteria: (i) ΔStability in active and inactive structures have opposite signs; or (ii) ΔStability in active and inactive structures have the same sign, but favorable interactions such as hydrogen bonds, salt bridge or pi-pi stacking exist in only one structure that can promote/block conformational changes upon activation. Totally, 15 and 20 mutations were predicted to be CAMs and CIMs, respectively. (Figure 6c and Figure 6—figure supplement 1).
cAMP accumulation assays
Request a detailed protocol(i)A2AR. The desired mutations were introduced into amino-terminally Flag tag-labeled human A2AR in the pcDNA3.1 vector (Invitrogen, Carlsbad, CA). This construct displayed equivalent pharmacological features to that of untagged human receptor based on radioligand binding and cAMP assays (Massink et al., 2015). The mutants were constructed by PCR-based site-directed mutagenesis (Muta-directTM kit, Beijing SBS Genetech Co., Ltd., China). Sequences of receptor clones were confirmed by DNA sequencing. HEK-293 cells (obtained from ATCC and confirmed as negative for mycoplasma contamination) were seeded onto 6-well cell culture plates. After overnight culture, the cells were transiently transfected with WT or mutant DNA using Lipofectamine 2000 transfection reagent (Invitrogen). After 24 hr, the transfected cells were seeded onto 384-well plates (3,000 cells per well). cAMP accumulation was measured using the LANCE cAMP kit (PerkinElmer, Boston, MA) according to the manufacturer’s instructions. Briefly, transfected cells were incubated for 40 min in assay buffer (DMEM, 1 mM 3-isobutyl-1-methylxanthine) with different concentrations of agonist [CGS21680 (179 pM to 50 μM)]. The reactions were stopped by addition of lysis buffer containing LANCE reagents. Plates were then incubated for 60 min at room temperature and time-resolved FRET signals were measured at 625 nm and 665 nm by an EnVision multilabel plate reader (PerkinElmer). The cAMP response is depicted relative to the maximal response of CGS21680 (100%) at the WT A2AR. (ii) 5-HT1B receptor. cAMP accumulation was measured using LANCE cAMP kit (PerkinElmer). Briefly, HEK293T (obtained from and certified by the Cell Bank at the Chinese Academy of Science and confirmed as negative for mycoplasma contamination) cells were transfected with plasmids bearing WT or mutant 5-HT1B receptor. Cells were collected 24 hr post-transfection and used to seed white poly-D-lysine coated 384-well plates at a density of 2,000 cells per well. Cells were incubated for a further 24 hr at 37°C. Cells were then incubated for 30 min in assay buffer (HBSS, 5 mM HEPES, 0.1% BSA, 0.5 mM 3-isobutyl-1-methylxanthine) with constant Forsklin (800 nM) and different concentrations of dihydroergotamine (DHE, 0.64 pM to 50 nM) at 37°C. The reactions were stopped by addition of lysis buffer containing LANCE reagents. Plates were then incubated for 60 min at room temperature, and time-resolved FRET signals were measured after excitation at 620 nm and 650 nm by EnVision (PerkinElmer).
CGS21680 binding assay
Request a detailed protocolCGS21680 (a specific adenosine A2A subtype receptor agonist) binding was analyzed using plasma membranes prepared from HEK-293 cells transiently expressing WT and mutant A2ARs. Approximately 1.2×108 transfected HEK-293 cells were harvested, suspended in 10 ml ice-cold membrane buffer (50 mM Tris-HCl, pH 7.4) and centrifuged for 5 min at 700 g. The resulting pellet was resuspended in ice-cold membrane buffer, homogenized by Dounce Homogenizer (Wheaton, Millville, NJ) and centrifuged for 20 min at 50,000 g. The pellet was resuspended, homogenized, centrifuged again and the precipitate containing the plasma membranes was then suspended in the membrane buffer containing protease inhibitor (Sigma-Aldrich, St. Louis, MO) and stored at −80°C. Protein concentration was determined using a protein BCA assay kit (Pierce Biotechnology, Pittsburgh, PA). For homogeneous binding, cell membrane homogenates (10 µg protein per well) were incubated in membrane binding buffer (50 mM Tris-HCl, 10 mM NaCl, 0.1 mM EDTA, pH 7.4) with constant concentration of [3H]-CGS21680 (1 nM, PerkinElmer) and serial dilutions of unlabeled CGS21680 (0.26 nM to 100 µM) at room temperature for 3 hr. Nonspecific binding was determined in the presence of 100 µM CGS21680. Following incubation, the samples were filtered rapidly in vacuum through glass fiber filter plates (PerkinElmer). After soaking and rinsing four times with ice-cold PBS, the filters were dried and counted for radioactivity in a MicroBeta2 scintillation counter (PerkinElmer).
Surface expression of A2ARs
Request a detailed protocolHEK293 cells were seeded into six-well plate and incubated overnight. After transient transfection with WT or mutant plasmids for 24 hr, the cells were collected and blocked with 5% BSA in PBS at room temperature for 15 min and incubated with primary anti-Flag antibody (1:100, Sigma-Aldrich) at room temperature for 1 hr. The cells were then washed three times with PBS containing 1% BSA followed by 1 hr incubation with anti-rabbit Alexa-488-conjugated secondary antibody (1:1000, Cell Signaling Technology, Danvers, MA) at 4°C in the dark. After three washes, the cells were resuspended in 200 µl of PBS containing 1% BSA for detection in a NovoCyte flow cytometer (ACEA Biosciences, San Diego, CA) utilizing laser excitation and emission wavelengths of 488 nm and 519 nm, respectively. For each assay point, approximately 15,000 cellular events were collected, and the total fluorescence intensity of positive expression cell population was calculated.
Plasmid constructs of 5-HT7 receptor
Request a detailed protocolA plasmid encoding the 5-HT7 receptor was obtained from PRESTO-Tango Kit produced by Addgene (Watertown, MA). 5-HT7 coding sequence was amplified and ligated into the pRluc8-N1 vector to produce WT 5-HT7-Rluc8. Mutant 5-HT7-Rluc8 receptors were generated from this plasmid using the Quikchange mutagenesis kit (Agilent Technologies, Santa Clara, CA). A plasmid encoding the Nluc-EPAC-VV cAMP sensor was kindly provided by Kirill Martemyanov (The Scripps Research Institute, Jupiter, FL) and has been described previously. (Masuho et al., 2015) All plasmid constructs were verified by DNA sequencing.
Cell transfection
Request a detailed protocolHEK293 cells cultured in 6-well plates were transiently transfected with the above plasmids (3.0 μg DNA) in growth medium using linear polyethyleneimine MAX (Polysciences, Warrington, PA) at an N/P ratio of 20 and used for experimentation 12–48 hr thereafter.
BRET cAMP and trafficking assays
Request a detailed protocolHEK 293 cells transiently transfected with WT or mutant 5-HT7 receptor plasmids plus either the Nluc-EPAC-VV cAMP sensor (at a 15:1 ratio) or Venus-kras (at a 1:8 ratio) (Tian et al., 2017) were incubated for 24 hr. After washing twice with PBS, they were transferred to opaque black 96-well plates. Steady-state BRET measurements were made using a Mithras LB940 photon-counting plate reader (Berthold Technologies GmbH, Bad Wildbad, Germany). Furimazine (NanoGlo; 1:1000, Promega) for cAMP measurement or coelenterazine h (5 µM; Nanolight, Pinetop, AZ) for trafficking assay was added followed by BRET signal detection and calculation at an emission intensity of 520–545 nm divided by that of 475–495 nm.
Data and materials availability
Request a detailed protocolThe open source code is available at GitHub (Zhou, 2019; copy archived at https://github.com/elifesciences-publications/RRCS). For availability of codes that were developed in-house, please contacts the corresponding authors. All data are available in the main text or the source data.
Data availability
All data generated or analysed during this study are included in the manuscript and supporting files. Source data files have been provided for Figures 1, 2, 6 and 7.
References
-
Protein residue contacts and prediction methodsMethods in Molecular Biology 1415:463–476.https://doi.org/10.1007/978-1-4939-3572-7_24
-
Strategies to discover unexpected targets for drugs active at G protein-coupled receptorsAnnual Review of Pharmacology and Toxicology 51:117–144.https://doi.org/10.1146/annurev-pharmtox-010510-100553
-
Molecular tinkering of G protein-coupled receptors: an evolutionary successThe EMBO Journal 18:1723–1729.https://doi.org/10.1093/emboj/18.7.1723
-
Crystallizing membrane proteins using lipidic mesophasesNature Protocols 4:706–731.https://doi.org/10.1038/nprot.2009.31
-
Emerging paradigm of intracellular targeting of G Protein-Coupled receptorsTrends in Biochemical Sciences 43:533–546.https://doi.org/10.1016/j.tibs.2018.04.003
-
A robotic system for crystallizing membrane and soluble proteins in lipidic mesophasesActa Crystallographica Section D Biological Crystallography 60:1795–1807.https://doi.org/10.1107/S0907444904019109
-
Advances in G protein-coupled receptor allostery: from function to structureMolecular Pharmacology 86:463–478.https://doi.org/10.1124/mol.114.094342
-
G protein-coupled receptor allosterism and complexingPharmacological Reviews 54:323–374.https://doi.org/10.1124/pr.54.2.323
-
Advances in computational techniques to study GPCR-Ligand recognitionTrends in Pharmacological Sciences 36:878–890.https://doi.org/10.1016/j.tips.2015.08.006
-
Structures of G protein-coupled receptors reveal new opportunities for drug discoveryDrug Discovery Today 20:1355–1364.https://doi.org/10.1016/j.drudis.2015.08.003
-
Empirical scoring functions: I. the development of a fast empirical scoring function to estimate the binding affinity of ligands in receptor complexesJournal of Computer-Aided Molecular Design 11:425–445.https://doi.org/10.1023/a:1007996124545
-
Computational design of ligand-binding membrane receptors with high selectivityNature Chemical Biology 13:715–723.https://doi.org/10.1038/nchembio.2371
-
Molecular switches in GPCRsCurrent Opinion in Structural Biology 55:114–120.https://doi.org/10.1016/j.sbi.2019.03.017
-
Desensitization of G protein-coupled receptors and neuronal functionsAnnual Review of Neuroscience 27:107–144.https://doi.org/10.1146/annurev.neuro.27.070203.144206
-
Methodological advances: the unsung heroes of the GPCR structural revolutionNature Reviews Molecular Cell Biology 16:69–81.https://doi.org/10.1038/nrm3933
-
Modulating G-protein-coupled receptors: from traditional pharmacology to allostericsTrends in Pharmacological Sciences 28:431–437.https://doi.org/10.1016/j.tips.2007.06.012
-
Rules of engagement: gpcrs and G proteinsACS Pharmacology & Translational Science 1:73–83.https://doi.org/10.1021/acsptsci.8b00026
-
Genetically encoded FRET biosensors to illuminate compartmentalised GPCR signallingTrends in Pharmacological Sciences 39:148–157.https://doi.org/10.1016/j.tips.2017.09.005
-
Trends in GPCR drug discovery: new agents, targets and indicationsNature Reviews Drug Discovery 16:829–842.https://doi.org/10.1038/nrd.2017.178
-
A G protein-coupled receptor at work: the rhodopsin modelTrends in Biochemical Sciences 34:540–552.https://doi.org/10.1016/j.tibs.2009.07.005
-
A conserved aromatic lock for the tryptophan rotameric switch in TM-VI of seven-transmembrane receptorsJournal of Biological Chemistry 285:3973–3985.https://doi.org/10.1074/jbc.M109.064725
-
Na+-mimicking ligands stabilize the inactive state of leukotriene B4 receptor BLT1Nature Chemical Biology 14:262–269.https://doi.org/10.1038/nchembio.2547
-
GPCR activation: a mutagenic spotlight on crystal structuresTrends in Pharmacological Sciences 34:67–84.https://doi.org/10.1016/j.tips.2012.11.002
-
Agonist induction and conformational selection during activation of a G-protein-coupled receptorTrends in Pharmacological Sciences 24:81–86.https://doi.org/10.1016/S0165-6147(02)00050-0
-
Generic GPCR residue numbers - aligning topology maps while minding the gapsTrends in Pharmacological Sciences 36:22–31.https://doi.org/10.1016/j.tips.2014.11.001
-
GPCRdb: an information system for G protein-coupled receptorsNucleic Acids Research 44:D356–D364.https://doi.org/10.1093/nar/gkv1178
-
Computational studies to predict or explain G protein coupled receptor polypharmacologyTrends in Pharmacological Sciences 35:658–663.https://doi.org/10.1016/j.tips.2014.10.009
-
Improved methods for building protein models in electron density maps and the location of errors in these modelsActa Crystallographica Section a Foundations of Crystallography 47:110–119.https://doi.org/10.1107/S0108767390010224
-
Diversity and modularity of G protein-coupled receptor structuresTrends in Pharmacological Sciences 33:17–27.https://doi.org/10.1016/j.tips.2011.09.003
-
Structure-function of the G protein-coupled receptor superfamilyAnnual Review of Pharmacology and Toxicology 53:531–556.https://doi.org/10.1146/annurev-pharmtox-032112-135923
-
Allosteric sodium in class A GPCR signalingTrends in Biochemical Sciences 39:233–244.https://doi.org/10.1016/j.tibs.2014.03.002
-
Visualization and analysis of non-covalent contacts using the protein contacts atlasNature Structural & Molecular Biology 25:185–194.https://doi.org/10.1038/s41594-017-0019-z
-
The Uppsala Electron-Density serverActa Crystallographica Section D Biological Crystallography 60:2240–2249.https://doi.org/10.1107/S0907444904013253
-
New G-protein-coupled receptor crystal structures: insights and limitationsTrends in Pharmacological Sciences 29:79–83.https://doi.org/10.1016/j.tips.2007.11.009
-
Homology Modeling: Methods and Protocols231–257, Methods of Protein Structure Comparison, Homology Modeling: Methods and Protocols, Humana Press.
-
Structural diversity of G protein-coupled receptors and significance for drug discoveryNature Reviews Drug Discovery 7:339–357.https://doi.org/10.1038/nrd2518
-
GPCR dynamics: structures in motionChemical Reviews 117:139–155.https://doi.org/10.1021/acs.chemrev.6b00177
-
Allosteric GPCR modulators: taking advantage of permissive receptor pharmacologyTrends in Pharmacological Sciences 28:382–389.https://doi.org/10.1016/j.tips.2007.06.004
-
Structure of bovine rhodopsin in a trigonal crystal formJournal of Molecular Biology 343:1409–1438.https://doi.org/10.1016/j.jmb.2004.08.090
-
Small molecule allosteric modulators of G-Protein-Coupled receptors: drug-target interactionsJournal of Medicinal Chemistry 62:24–45.https://doi.org/10.1021/acs.jmedchem.7b01844
-
Evolutionary trace of G protein-coupled receptors reveals clusters of residues that determine global and class-specific functionsJournal of Biological Chemistry 279:8126–8132.https://doi.org/10.1074/jbc.M312671200
-
Sodium ion binding pocket mutations and adenosine A2A receptor functionMolecular Pharmacology 87:305–313.https://doi.org/10.1124/mol.114.095737
-
Structural insights into binding specificity, efficacy and Bias of a β2AR partial agonistNature Chemical Biology 14:1059–1066.https://doi.org/10.1038/s41589-018-0145-x
-
Allosteric modulation of G protein-coupled receptorsAnnual Review of Pharmacology and Toxicology 47:1–51.https://doi.org/10.1146/annurev.pharmtox.47.120505.105159
-
Structure-inspired design of β-arrestin-biased ligands for aminergic GPCRsNature Chemical Biology 14:126–134.https://doi.org/10.1038/nchembio.2527
-
Mutational analysis of the Rhodopsin gene in sector retinitis pigmentosaOphthalmic Genetics 36:239–243.https://doi.org/10.3109/13816810.2014.958862
-
Orphan receptor ligand discovery by pickpocketing pharmacological neighborsNature Chemical Biology 13:235–242.https://doi.org/10.1038/nchembio.2266
-
Ligand binding and micro-switches in 7tm receptor structuresTrends in Pharmacological Sciences 30:249–259.https://doi.org/10.1016/j.tips.2009.02.006
-
Activation of rhodopsin: new insights from structural and biochemical studiesTrends in Biochemical Sciences 26:318–324.https://doi.org/10.1016/S0968-0004(01)01799-6
-
What is biased efficacy? defining the relationship between intrinsic efficacy and free energy couplingTrends in Pharmacological Sciences 35:639–647.https://doi.org/10.1016/j.tips.2014.09.010
-
Intracellular receptor modulation: novel approach to target GPCRsTrends in Pharmacological Sciences 39:547–559.https://doi.org/10.1016/j.tips.2018.03.002
-
GPCRdb in 2018: adding GPCR structure models and ligandsNucleic Acids Research 46:D440–D446.https://doi.org/10.1093/nar/gkx1109
-
Functional characterization of the molecular defects causing nephrogenic diabetes insipidus in eight familiesJournal of Clinical Endocrinology & Metabolism 85:1703–1710.https://doi.org/10.1210/jc.85.4.1703
-
Genetic variations and polymorphisms of G protein-coupled receptors: functional and therapeutic implicationsAnnual Review of Pharmacology and Toxicology 41:593–624.https://doi.org/10.1146/annurev.pharmtox.41.1.593
-
Cryo-EM in drug discovery: achievements, limitations and prospectsNature Reviews Drug Discovery 17:471–492.https://doi.org/10.1038/nrd.2018.77
-
Discovery of new GPCR ligands to illuminate new biologyNature Chemical Biology 13:1143–1151.https://doi.org/10.1038/nchembio.2490
-
Learning from the past: evolution of GPCR functionsTrends in Pharmacological Sciences 28:117–121.https://doi.org/10.1016/j.tips.2007.01.001
-
Structure of an allosteric modulator bound to the CB1 cannabinoid receptorNature Chemical Biology 15:1199–1205.https://doi.org/10.1038/s41589-019-0387-2
-
X-ray structures of endothelin ETB receptor bound to clinical antagonist bosentan and its analogNature Structural & Molecular Biology 24:758–764.https://doi.org/10.1038/nsmb.3450
-
GPCR drug discovery: integrating solution NMR data with crystal and cryo-EM structuresNature Reviews Drug Discovery 18:59–82.https://doi.org/10.1038/nrd.2018.180
-
Validation of ligands in macromolecular structures determined by X-ray crystallographyActa Crystallographica Section D Structural Biology 74:228–236.https://doi.org/10.1107/S2059798318002541
-
Pharmacogenomic and structural analysis of constitutive g protein-coupled receptor activityAnnual Review of Pharmacology and Toxicology 47:53–87.https://doi.org/10.1146/annurev.pharmtox.47.120505.105126
-
Biased signalling: from simple switches to allosteric microprocessorsNature Reviews Drug Discovery 17:243–260.https://doi.org/10.1038/nrd.2017.229
-
Olfactory receptors: g protein-coupled receptors and beyondJournal of Neurochemistry 109:1570–1583.https://doi.org/10.1111/j.1471-4159.2009.06085.x
-
Serial femtosecond crystallography of G Protein-Coupled receptorsAnnual Review of Biophysics 47:377–397.https://doi.org/10.1146/annurev-biophys-070317-033239
-
Unifying family A GPCR theories of activationPharmacology & Therapeutics 143:51–60.https://doi.org/10.1016/j.pharmthera.2014.02.004
-
UniProt: the universal protein knowledgebaseNucleic Acids Research 45:D158–D169.https://doi.org/10.1093/nar/gkw1099
-
G protein-coupled receptor mutations and human genetic diseaseMethods in Molecular Biology 1175:153–187.https://doi.org/10.1007/978-1-4939-0956-8_8
-
Action of molecular switches in GPCRs--theoretical and experimental studiesCurrent Medicinal Chemistry 19:1090–1109.https://doi.org/10.2174/092986712799320556
-
Chemical Diversity in the G Protein-Coupled Receptor SuperfamilyTrends in Pharmacological Sciences 39:494–512.https://doi.org/10.1016/j.tips.2018.02.004
-
G protein-coupled receptors: mutations and endocrine diseasesNature Reviews Endocrinology 7:362–372.https://doi.org/10.1038/nrendo.2011.20
-
Improved protein-ligand docking using GOLDProteins: Structure, Function, and Bioinformatics 52:609–623.https://doi.org/10.1002/prot.10465
-
Accurate de novo prediction of protein contact map by ultra-deep learning modelPLOS Computational Biology 13:e1005324.https://doi.org/10.1371/journal.pcbi.1005324
-
Development and crystallization of a minimal thermostabilised G protein-coupled receptorProtein Expression and Purification 65:204–213.https://doi.org/10.1016/j.pep.2009.01.014
-
The molecular basis of G Protein-Coupled receptor activationAnnual Review of Biochemistry 87:897–919.https://doi.org/10.1146/annurev-biochem-060614-033910
-
Therapeutic potential of β-arrestin- and G protein-biased agonistsTrends in Molecular Medicine 17:126–139.https://doi.org/10.1016/j.molmed.2010.11.004
-
Mechanisms of signalling and biased agonism in G protein-coupled receptorsNature Reviews Molecular Cell Biology 19:638–653.https://doi.org/10.1038/s41580-018-0049-3
-
Protein Allostery in Drug Discovery225–251, GPCR Allosteric Modulator Discovery, Protein Allostery in Drug Discovery, Springer.
-
Allosteric mechanisms underlie GPCR signaling to SH3-domain proteins through arrestinNature Chemical Biology 14:876–886.https://doi.org/10.1038/s41589-018-0115-3
-
The role of water and sodium ions in the activation of the μ-opioid receptorAngewandte Chemie International Edition 52:10112–10115.https://doi.org/10.1002/anie.201302244
Article and author information
Author details
Funding
Medical Research Council (MC_U105185859)
- M Madan Babu
Novo Nordisk-CAS Research (NNCAS-2017-1-CC)
- Dehua Yang
Young Talent Program of Shanghai
- Suwen Zhao
Shanghai Science and Technology Development Fund (16ZR1448500)
- Suwen Zhao
Shanghai Science and Technology Development Fund (16ZR1407100)
- Antao Dai
National Natural Science Foundation of China (21704064)
- Qingtong Zhou
National Natural Science Foundation of China (81573479)
- Dehua Yang
National Natural Science Foundation of China (81773792)
- Dehua Yang
National Natural Science Foundation of China (31971178)
- Suwen Zhao
National Mega R&D Program for Drug Discovery (2018ZX09735-001)
- Ming-Wei Wang
National Key R&D Program of China (2016YFC0905900)
- Suwen Zhao
National Mega R&D Program for Drug Discovery (2018ZX09711002-002-005)
- Dehua Yang
National Key R&D Program of China (2018YFA0507000)
- Ming-Wei Wang
- Suwen Zhao
National Natural Science Foundation of China (81872915)
- Ming-Wei Wang
National Institute of General Medical Sciences (GM130142)
- Nevin A Lambert
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
This work was partially supported by National Natural Science Foundation of China grants 31971178 (SZ), 81872915 (M-WW), 21704064 (QZ), 81573479 (DHY) and 81773792 (DHY), National Key R and D Program of China grants 2016YFC0905900 (SZ) and 2018YFA0507000 (SZ and M-WW), National Mega R and D Program for Drug Discovery grants 2018ZX09711002–002–005 (DHY) and 2018ZX09735–001 (M-WW), Shanghai Science and Technology Development Fund grant 16ZR1407100 (ATD), Novo Nordisk-CAS Research Fund grant NNCAS-2017–1-CC (DHY), the Medical Research Council MC_U105185859 (MMB), National Institutes of Health Grants GM130142 (NAL) and annual overhead support from ShanghaiTech University and Chinese Academy of Sciences. We thank A Sali, MA Hanson and AJ Venkatakrishnan for valuable discussions, YM Xu for technical assistance, and AP IJzerman for providing the WT A2AR plasmid.
Copyright
© 2019, Zhou et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 25,522
- views
-
- 3,809
- downloads
-
- 396
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Computational and Systems Biology
- Neuroscience
Diffusional kurtosis imaging (DKI) is a methodology for measuring the extent of non-Gaussian diffusion in biological tissue, which has shown great promise in clinical diagnosis, treatment planning, and monitoring of many neurological diseases and disorders. However, robust, fast, and accurate estimation of kurtosis from clinically feasible data acquisitions remains a challenge. In this study, we first outline a new accurate approach of estimating mean kurtosis via the sub-diffusion mathematical framework. Crucially, this extension of the conventional DKI overcomes the limitation on the maximum b-value of the latter. Kurtosis and diffusivity can now be simply computed as functions of the sub-diffusion model parameters. Second, we propose a new fast and robust fitting procedure to estimate the sub-diffusion model parameters using two diffusion times without increasing acquisition time as for the conventional DKI. Third, our sub-diffusion-based kurtosis mapping method is evaluated using both simulations and the Connectome 1.0 human brain data. Exquisite tissue contrast is achieved even when the diffusion encoded data is collected in only minutes. In summary, our findings suggest robust, fast, and accurate estimation of mean kurtosis can be realised within a clinically feasible diffusion-weighted magnetic resonance imaging data acquisition time.
-
- Computational and Systems Biology
- Developmental Biology
The initially homogeneous epithelium of the early Drosophila embryo differentiates into regional subpopulations with different behaviours and physical properties that are needed for morphogenesis. The factors at top of the genetic hierarchy that control these behaviours are known, but many of their targets are not. To understand how proteins work together to mediate differential cellular activities, we studied in an unbiased manner the proteomes and phosphoproteomes of the three main cell populations along the dorso-ventral axis during gastrulation using mutant embryos that represent the different populations. We detected 6111 protein groups and 6259 phosphosites of which 3398 and 3433 were differentially regulated, respectively. The changes in phosphosite abundance did not correlate with changes in host protein abundance, showing phosphorylation to be a regulatory step during gastrulation. Hierarchical clustering of protein groups and phosphosites identified clusters that contain known fate determinants such as Doc1, Sog, Snail, and Twist. The recovery of the appropriate known marker proteins in each of the different mutants we used validated the approach, but also revealed that two mutations that both interfere with the dorsal fate pathway, Toll10B and serpin27aex do this in very different manners. Diffused network analyses within each cluster point to microtubule components as one of the main groups of regulated proteins. Functional studies on the role of microtubules provide the proof of principle that microtubules have different functions in different domains along the DV axis of the embryo.