Pharmacological hallmarks of allostery at the M4 muscarinic receptor elucidated through structure and dynamics

  1. Ziva Vuckovic
  2. Jinan Wang
  3. Vi Pham
  4. Jesse I Mobbs
  5. Matthew J Belousoff
  6. Apurba Bhattarai
  7. Wessel AC Burger
  8. Geoff Thompson
  9. Mahmuda Yeasmin
  10. Vindhya Nawaratne
  11. Katie Leach
  12. Emma T van der Westhuizen
  13. Elham Khajehali
  14. Yi-Lynn Liang
  15. Alisa Glukhova
  16. Denise Wootten
  17. Craig W Lindsley
  18. Andrew Tobin
  19. Patrick Sexton
  20. Radostin Danev
  21. Celine Valant  Is a corresponding author
  22. Yinglong Miao  Is a corresponding author
  23. Arthur Christopoulos  Is a corresponding author
  24. David M Thal  Is a corresponding author
  1. Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Australia
  2. Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, United States
  3. ARC Centre for Cryo-electron Microscopy of Membrane Proteins, Monash Institute of Pharmaceutical Sciences, Monash University, Australia
  4. Department of Pharmacology, Warren Center for Neuroscience Drug Discovery and Department of Chemistry, Warren Center for Neuroscience Drug Discovery, Vanderbilt University, United States
  5. The Centre for Translational Pharmacology, Advanced Research Centre (ARC), College of Medical, Veterinary and Life Sciences, University of Glasgow, United Kingdom
  6. Graduate School of Medicine, University of Tokyo, Japan
  7. Neuromedicines Discovery Centre, Monash University, Australia

Abstract

Allosteric modulation of G protein-coupled receptors (GPCRs) is a major paradigm in drug discovery. Despite decades of research, a molecular-level understanding of the general principles that govern the myriad pharmacological effects exerted by GPCR allosteric modulators remains limited. The M4 muscarinic acetylcholine receptor (M4 mAChR) is a validated and clinically relevant allosteric drug target for several major psychiatric and cognitive disorders. In this study, we rigorously quantified the affinity, efficacy, and magnitude of modulation of two different positive allosteric modulators, LY2033298 (LY298) and VU0467154 (VU154), combined with the endogenous agonist acetylcholine (ACh) or the high-affinity agonist iperoxo (Ipx), at the human M4 mAChR. By determining the cryo-electron microscopy structures of the M4 mAChR, bound to a cognate Gi1 protein and in complex with ACh, Ipx, LY298-Ipx, and VU154-Ipx, and applying molecular dynamics simulations, we determine key molecular mechanisms underlying allosteric pharmacology. In addition to delineating the contribution of spatially distinct binding sites on observed pharmacology, our findings also revealed a vital role for orthosteric and allosteric ligand–receptor–transducer complex stability, mediated by conformational dynamics between these sites, in the ultimate determination of affinity, efficacy, cooperativity, probe dependence, and species variability. There results provide a holistic framework for further GPCR mechanistic studies and can aid in the discovery and design of future allosteric drugs.

Editor's evaluation

This important work advances our understanding of the structural basis of allosteric modulation of the M4 muscarinic receptor but has broad implications for GPCRs. The evidence supporting the conclusions is exceptional, with multiple cryo-EM structures that are complemented by excellent pharmacological and dynamics studies.

https://doi.org/10.7554/eLife.83477.sa0

Introduction

Over the past 40 y, there have been major advances to the analytical methods that allow for the quantitative determination of the pharmacological parameters that characterize G protein-coupled receptor (GPCR) signaling and allosteric modulation (Figure 1A and B). These analytical methods are based on the operational model of agonism (Black and Leff, 1983) and have been extended or modified to account for allosteric modulation (Leach et al., 2007), biased agonism (Kenakin, 2012), and even biased allosteric modulation (Slosky et al., 2021). Collectively, these models and subsequent key parameters (Figure 1B) are used to guide allosteric drug screening, selectivity, efficacy, and ultimately, clinical utility, and provide the foundation for modern GPCR drug discovery (Wootten et al., 2013). Yet, a systematic understanding of how these pharmacological parameters relate to the molecular structure and dynamics of GPCRs remains elusive.

Figure 1 with 2 supplements see all
Pharmacological characterization of the positive allosteric modulators (PAMs), LY298 and VU154, with acetylcholine (ACh) and iperoxo (Ipx) at the human M4 muscarinic acetylcholine receptor (mAChR).

(A) Schematic of the pharmacological parameters that define effects of orthosteric and allosteric ligands on a G protein-coupled receptor (GPCR). (B) A simplified schematic diagram of the Black–Leff operational model to quantify agonism, allosteric modulation, and agonist bias with pharmacological parameters defined (Black and Leff, 1983). (C) 2D chemical structures of the orthosteric and allosteric ligands used in this study. (D–G) Key pharmacological parameters for interactions between orthosteric and allosteric ligands in [3H]-N-methylscopolamine ([3H]-NMS) binding assays. (D) Equilibrium binding affinities (pKi and pKB) and (E) the degree of binding modulation (α) between the agonists and PAMs resulting in the modified binding affinities (F) α/KA and (G) α/KB. (H–K) Key pharmacological parameters relating to Gαi1 activation for interactions between orthosteric and allosteric ligands measured with the TruPath assay (Figure 1—figure supplement 1). (H) The signaling efficacy (τA and τB) and (I) transduction coupling coefficients (log (τ/K)) of each ligand. (J) The functional cooperativity (αβ) between ligands and (K) the efficacy modulation (β) between ligands. All data are mean ± SEM of three or more independent experiments performed in duplicate or triplicate with the pharmacological parameters determined using a global fit of the data. The error in (F, G, K) was propagated using the square root of the sum of the squares. See Table 1. Concentration–response curves are shown in Figure 1—figure supplement 1.

The muscarinic acetylcholine receptors (mAChRs) are an important family of five Class A GPCRs that have long served as model systems for understanding GPCR allostery (Conn et al., 2009). The mAChRs have been notoriously difficult to exploit therapeutically and selectively due to high-sequence conservation within their orthosteric binding domains (Burger et al., 2018). However, the discovery of highly selective positive allosteric modulators (PAMs) for some mAChR subtypes has paved the way for novel approaches to exploit these high-value drug targets (Chan et al., 2008; Gentry et al., 2014; Marlo et al., 2009). X-ray crystallography and cryo-electron microscopy (cryo-EM) have been used to determine inactive state structures for all five mAChR subtypes (Haga et al., 2012; Kruse et al., 2012; Thal et al., 2016; Vuckovic et al., 2019) and active state structures of the M1 and M2 mAChRs (Maeda et al., 2019). For the M2 mAChR, this includes structures co-bound with the high-affinity agonist iperoxo (Ipx) and the PAM LY2119620 in complex with a G protein mimetic nanobody (Kruse et al., 2013) and the transducers Go (Maeda et al., 2019) and β-arrestin1 (Staus et al., 2020). These M2 mAChR structures were foundational to validating the canonical mAChR allosteric site but are limited to only one agonist (iperoxo) and one PAM (LY2119620) and do not account for the vast pharmacological properties of ligands targeting mAChRs. A recent nuclear magnetic resonance (NMR) study of the M2 mAChR revealed differences in the conformational landscape of the M2 mAChR when bound to different agonists, but no clear link was established between the properties of the ligands and the conformational states of the receptor (Xu et al., 2019). The M4 mAChR subtype is of major therapeutic interest due to its expression in regions of the brain that are rich in dopamine and dopamine receptors, where it regulates dopaminergic neurons involved in cognition, psychosis, and addiction (Bymaster et al., 2003; Dencker et al., 2011; Foster et al., 2016; Tzavara et al., 2004). Importantly, these findings have been supported by studies utilizing novel PAMs that are highly selective for the M4 mAChR (Bubser et al., 2014; Chan et al., 2008; Leach et al., 2010; Suratman et al., 2011). Among these, LY2033298 (LY298) was the first reported highly selective PAM of the M4 mAChR and displayed antipsychotic efficacy in a preclinical animal model of schizophrenia (Chan et al., 2008). Despite LY298 being one of the best characterized M4 mAChR PAMs, its therapeutic potential has been limited by numerous factors, including its chemical scaffold, which has been difficult to optimize with respect to its molecular allosteric parameters (Figure 1C) and variability of response between species (Suratman et al., 2011; Wood et al., 2017b). In the search for better chemical scaffolds, the PAM, VU0467154 (VU154), was subsequently discovered. VU154 showed robust efficacy in preclinical rodent models; however, it also exhibited species selectivity that prevented its clinical translation (Bubser et al., 2014). Collectively, LY298 and VU154 are exemplar tool molecules that highlight the promises and the challenges in understanding and optimizing allosteric GPCR drug activity for translational and clinical applications.

Herein, by examining the pharmacology of the PAMs LY298 and VU154 with the agonists ACh and Ipx across radioligand binding assays and two different signaling assays and analyzing these results with modern analytical methods, we determined the key parameters that describe signaling and allostery for these ligands. To investigate a structural basis for these pharmacological parameters, we used cryo-EM to determine high-resolution structures of the M4 mAChR in complex with a cognate Gi1 heterotrimer and ACh and Ipx. We also determined structures of receptor complexes with Ipx co-bound with the PAMs LY298 or VU154. Moreover, because protein allostery is a dynamic process (Changeux and Christopoulos, 2016), we performed all-atom simulations using the Gaussian accelerated molecular dynamics (GaMD) enhanced sampling method (Draper-Joyce et al., 2021; Miao et al., 2015; Wang et al., 2021a) on the M4 mAChR using the cryo-EM structures. The structures and GaMD simulations, in combination with detailed molecular pharmacology and receptor mutagenesis experiments, provide fundamental insights into the molecular mechanisms underpinning the hallmarks of GPCR allostery. To further validate these findings, we investigated the differences in the selectivity of VU154 between the human and mouse receptors and established a structural basis for species selectivity. Collectively, these results will enable future GPCR drug discovery research and potentially lead to the development of next generation M4 mAChR PAMs.

Results

Pharmacological characterization of M4 mAChR PAMs with ACh and Ipx

The pharmacology of LY298 or VU154 interacting with ACh has been well characterized in binding and functional assays at the M4 mAChR (Bubser et al., 2014; Chan et al., 2008; Gould et al., 2016; Leach et al., 2010; Suratman et al., 2011; Thal et al., 2016). However, their pharmacology with Ipx has not been reported. Therefore, we characterized both PAMs with ACh and Ipx in binding and in two different functional assays to provide a thorough foundational comparative characterization of the pharmacological parameters of these ligands from the same study.

We first used radioligand binding assays (Figure 1—figure supplement 1A) to determine the binding affinities (i.e., equilibrium dissociation constants) of ACh and Ipx (KA) for the orthosteric site and of LY298 and VU154 (KB) for the allosteric site of the unoccupied human M4 mAChR (Figure 1D), along with the degree of binding cooperativity (α) between the agonists and PAMs when the two are co-bound (Figure 1E). Analysis of these experiments revealed that LY298 and VU154 have very similar binding affinities for the allosteric site with values (expressed as negative logarithms; pKB) of 5.65 ± 0.07 and 5.83 ± 0.12, respectively (Table 1), in accordance with previous studies (Bubser et al., 2014; Leach et al., 2011). Both PAMs potentiated the binding affinity of ACh and Ipx (Figure 1E), with the effect being greatest between LY298 and ACh (~400-fold increase in binding affinity). Comparatively, the positive cooperativity between VU154 and ACh was only 40-fold. When Ipx was used as the agonist, the binding affinity modulation mediated by both PAMs was more modest, characterized by an approximately 72-fold potentiation for the combination of Ipx and LY298, and 10-fold potentiation for the combination of Ipx and VU154. These results indicate probe-dependent effects (Valant et al., 2012) with respect to the ability of either PAM to modulate the affinity of each agonist (Figure 1F and G). A probe-dependent effect was also observed with the radioligand, [3H]-NMS, evidenced by a reduction in specific radioligand binding due to negative cooperativity between the antagonist probe and LY298, which has been previously reported (Chan et al., 2008; Leach et al., 2010; Suratman et al., 2011; Thal et al., 2016). It is important to note that binding affinity modulation is thermodynamically reciprocal at equilibrium, and the affinities of LY298 and VU154 were thus also increased in the agonist bound state (Figure 1—figure supplement 1A). This results in LY298 having a fivefold higher binding affinity than VU154 when agonists are bound (Table 1).

Table 1
Pharmacological parameters from radioligand binding and functional experiments.
[3H]-NMS saturation binding on stable M4 mAChR CHO cells
ConstructsSites per cell*pKD
Human WT M4 mAChR598,111 ± 43,067 (7)9.76 ± 0.05 (7)
Mouse WT M4 mAChR21,027 ± 2188 (3)9.76 ± 0.05 (3)
Human D432E M4 mAChR126,377 ± 10,066 (3)9.60 ± 0.07 (3)
Human T433R M4 mAChR157,442 ± 36,658 (6)9.64 ± 0.09 (6)
Human V91L, D432E, T433R M4 mAChR205,771 ± 20,975 (4)9.58 ± 0.08 (4)
[3H]-NMS interaction binding assays between ACh or Ipx and LY298 or VU154 on stable M4 mAChR constructs in Flp-In CHO cells
ConstructsPAMpKi ACh pKi Ipx pKB PAM log αACh §log αIpx §
Human WT M4 mAChRLY2984.50 ± 0.06 (4)8.30 ± 0.06 (4)5.65 ± 0.07 (8) 2.59 ± 0.10 (4)1.86 ± 0.10 (4)
VU1544.40 ± 0.09 (4)8.19 ± 0.06 (8)5.83 ± 0.11 (12) 1.61 ± 0.13 (4)1.03 ± 0.10 (8)
Mouse WT M4 mAChRLY2984.52 ± 0.07 (4)8.55 ± 0.06 (4)5.74 ± 0.07 (8) 1.78 ± 0.10 (4)1.30 ± 0.11 (4)*
VU1544.59 ± 0.06 (4)8.57 ± 0.06 (3)6.07 ± 0.09 (7) 2.43 ± 0.10 (4)1.75 ± 0.12 (3)*
Human D432E M4 mAChRLY298N.T.8.28 ± 0.04 (5)5.86 ± 0.07 (5)N.T.1.59 ± 0.06 (5)
VU154N.T.8.27 ± 0.06 (6)6.21 ± 0.12 (6)N.T.1.04 ± 0.09 (6)
Human T433R M4 mAChRLY298N.T.8.05 ± 0.08 (5)5.04 ± 0.04 (5)*N.T.1.91 ± 0.11 (5)
VU154N.T.7.88 ± 0.04 (5)5.50 ± 0.08 (5)N.T.1.67 ± 0.07 (5)*
Human V91L, D432E, T433R M4 mAChRLY298N.T.7.95 ± 0.10 (4)5.29 ± 0.26 (4)N.T.1.80 ± 0.22 (4)
VU154N.T.7.89 ± 0.12 (4)6.34 ± 0.16 (4)*N.T.1.35 ± 0.16 (4)
i1 activation (TruPath) interaction assays between ACh or Ipx and LY298 or VU154 on transiently expressed M4 mAChR constructs in HEK293A cells
ConstructsPAMlog τ ACh**log τ Ipx**pKB PAM log τ PAM**log αβACh††log αβIpx††
Human WT M4 mAChRLY2982.71 ± 0.14 (4)1.49 ± 0.12 (4)= 5.651.02 ± 0.03 (8) 2.01 ± 0.14 (4)1.96 ± 0.16 (4)
VU154= 5.83–0.55 ± 0.08 (8) 1.22 ± 0.13 (4)0.20 ± 0.13 (4)
pERK1/2 interaction assays between ACh or Ipx and LY298 or VU154 on stable M4 mAChR constructs in Flp-In CHO cells
ConstructsPAMlog τ ACh**log τ Ipx**pKB PAM log τC PAM ‡ ‡log αβACh††log αβIpx††
Human WT M4 mAChRLY2983.27 ± 0.06 (8) 1.74 ± 0.03 (16) = 5.651.19 ± 0.05 (12)**2.29 ± 0.22 (4)1.08 ± 0.28 (8)
VU154= 5.830.11 ± 0.05 (12)**0.88 ± 0.23 (4)0.66 ± 0.15 (8)
Mouse WT M4 mAChRLY298N.T.N.D.= 5.741.32 ± 0.07 (5)N.T.1.24 ± 0.12 (4)
VU154N.T.N.D.= 6.071.47 ± 0.08 (5) § §N.T.2.08 ± 0.15 (5) § §
Human D432E M4 mAChRLY298N.T.N.D.= 5.861.34 ± 0.08 (5)N.T.1.37 ± 0.28 (5)
VU154N.T.N.D.= 6.210.78 ± 0.08 (5) § §N.T.1.02 ± 0.15 (5)
Human T433R M4 mAChRLY298N.T.N.D.= 5.041.73 ± 0.13 (5) § §N.T.1.85 ± 0.28 (5)
VU154N.T.N.D.= 5.500.95 ± 0.12 (5) § §N.T.1.18 ± 0.14 (5)
Human V91L, D432E, T433R M4 mAChRLY298N.T.N.D.= 5.291.62 ± 0.09 (5) § §N.T.1.64 ± 0.30 (5)
VU154N.T.N.D.= 6.340.68 ± 0.06 (5) § §N.T.1.34 ± 0.11 (5) § §
  1. Values represent the mean ± SEM with the number of independent experiments shown in parenthesis.

  2. N.T.: not tested; N.D.: not determined; Ach, acetylcholine; Ipx: iperoxo; PAM: positive allosteric modulator.

  3. *

    Number of [3H]-NMS binding sites per cell.

  4. Negative logarithm of the radioligand equilibrium dissociation constant.

  5. Negative logarithm of the orthosteric (pKi) or allosteric (pKB) equilibrium dissociation constant.

  6. §

    Logarithm of the binding cooperativity factor between the agonist (ACh or Ipx) and the PAM (LY298 or VU154).

  7. Parameter was determined in a shared global analysis between agonists.

  8. **

    Logarithm of the operational efficacy parameter determined using the Operational Model of Agonism.

  9. ††

    Logarithm of the functional cooperativity factor between the agonist (ACh or Ipx) and the PAM (LY298 or VU154).

  10. ‡ ‡

    logτC = logarithm of the operational efficacy parameter corrected for receptor expression (methods in Appendix 1).

  11. § §

    Values from pKB PAM, log αIpx, log τC PAM, and log αβIpx that are significantly different from human WT M4 mAChR (p<0.05) calculated by a one-way ANOVA with a Dunnett’s post-hoc test.

We subsequently used the BRET-based TruPath assay (Olsen et al., 2020) as a proximal measure of G protein activation with Gαi1 (Figure 1—figure supplement 1B). We also used a more amplified downstream signaling assay, extracellular signal-regulated kinases 1/2 phosphorylation (pERK1/2), that is also dependent on Gi activation (Figure 1—figure supplement 2A), to measure the cell-based activity of each PAM with each agonist. These signaling assays allowed us to determine the efficacy of the agonists (τA) and the PAMs (τB) (Figure 1H, Figure 1—figure supplement 2B). Importantly, efficacy (τ), as defined from the Black–Leff operational model of agonism (Black and Leff, 1983), is determined by the ability of an agonist to promote an active receptor conformation, the receptor density (Bmax), and the subsequent ability of a cellular system to generate a response (Figure 1B). Notably, in both signaling assays, the rank order of efficacy was ACh > Ipx > LY298 > VU154. We subsequently calculated the transducer coupling coefficient (τ/K) (Figure 1I, Figure 1—figure supplement 2C), a parameter often used as a starting point to quantify biased agonism (Kenakin et al., 2012) and that is specific to the intact cellular environment in which a given response occurs. Thus the dissociation constant (K) in the transduction coefficient subsumes the affinity for the ground state (non-bound) receptor, in addition to any isomerization states of the receptor that ultimately yield cellular responses (Kenakin and Christopoulos, 2013). Consequently, in both assays, the rank order of transducer coupling was Ipx >> ACh ~ LY298 > VU154 due to Ipx having a higher binding affinity for the receptor. Overall, these results indicate that although ACh is a more efficacious agonist than Ipx, it has lower transducer coupling coefficient. In contrast, LY298 has both better efficacy and transducer coupling coefficient than VU154 (Table 1).

The signaling assays and use of an operational model of allosterism also allowed for the determination of the functional cooperativity (αβ) exerted by the PAMs (Figure 1J, Figure 1—figure supplement 2D), which is a composite parameter accounting for both binding (α) and efficacy (β) modulation. Notably, VU154 displayed lower positive functional cooperativity with ACh than LY298. Strikingly, VU154 had negligible functional modulation with Ipx, in contrast to the cooperativity observed with ACh in the TruPath assay. The tenfold difference in αβ values for VU154 between ACh and Ipx highlights the dependence of the orthosteric probe used in the assay (i.e. probe dependence); on this basis, VU154 would be classified as a ‘neutral’ allosteric ligand (not a PAM) with Ipx in the TruPath assay, that is, VU154 still binds to the allosteric site, but displays neutral cooperativity (αβ = 1) with Ipx (Table 1).

The degree of efficacy modulation (β) that the PAMs have on the agonists can be calculated by subtracting the binding modulation (α) from the functional modulation (αβ) (Figure 1K, Figure 1—figure supplement 2E). A caveat of this analysis is that errors for β are higher due to the error being propagated between calculations. Ideally, the degree of efficacy modulation would be determined in an experimental system where the maximal efficacy of system is not reached by the agonists alone (Berizzi et al., 2016). Nevertheless, our analysis shows the PAMs LY298 and VU154 appear to have a slight negative to neutral effect on agonist efficacy in the Gi1 TruPath and pERK1/2 assays (Table 1), suggesting that the predominant allosteric effect exerted by these PAMs is mediated through modulation of binding affinity.

Collectively, our extensive analysis on the pharmacology of LY298 and VU154 with ACh and Ipx offers detailed insight into the key differences between these ligands across a range of pharmacological properties: ligand binding, probe dependence, efficacy, agonist–receptor–transducer interactions, and allosteric modulation (Figure 1, Table 1). We hypothesized that structures of the human M4 mAChR in complex with different agonists and PAMs combined with molecular dynamic simulations could provide high-resolution molecular insights into the different pharmacological profiles of these ligands.

Determination of M4R-Gi1 complex structures

Similar to the approach used in prior determination of active-state structures of the M1 and M2 mAChRs (Maeda et al., 2019), we used a human M4 mAChR construct that lacked residues 242–387 of the third intracellular loop to improve receptor expression and purification, and made complexes of the receptor with Gi1 protein and either the endogenous agonist, ACh, or Ipx. Due to the higher affinity of Ipx compared to ACh (Schrage et al., 2013), we utilized Ipx to form additional M4R-Gi1 complexes with or without the co-addition of either LY298 or VU154. In all instances, complex formation was initiated by combining purified M4 mAChR immobilized on anti-FLAG resin with detergent solubilized Gi1 membranes, a single-chain variable fragment (scFv16) that binds Gαi and Gβ, and the addition of apyrase to remove guanosine 5′-diphosphate (Maeda et al., 2018). For this study, we used a Gi1 heterotrimer composed of a dominant negative form of human Gαi1, and human Gβ1 and Gγ2. (Liang et al., 2018b). Vitrified samples of each complex were imaged using conventional cryo-TEM on a Titan Krios microscope (Danev et al., 2021).

The structures of ACh-, Ipx-, LY298-Ipx-, and VU154-Ipx-bound M4R-Gi1 complexes were determined to resolutions of 2.8, 2.8, 2.4, and 2.5 Å, respectively (Figure 2A, Figure 2—figure supplement 1, Table 2). For the ACh-bound M4R-Gi1 complex, an additional focus refinement yielded an improved map of the receptor and binding site (2.75 Å) for modeling (Figure 2—figure supplements 2 and 3). The cryo-EM density maps for all complexes were sufficient for confident placement of backbone and sidechains for most of the receptor, Gi1, and scFv16, and the bound ligands with exception of the alkyne bond of Ipx, which was consistent with prior cryo-EM studies (Maeda et al., 2019; Figure 2B, Figure 2—figure supplement 3).

Figure 2 with 5 supplements see all
Cryo-electron microscopy (cryo-EM) structures of the M4R-Gi1-scFv16 complexes.

(A) Cryo-EM maps of Ipx-bound M4R-Gi1-scFv16 complex with views from the membrane and the extracellular surface. Cryo-EM maps of the other ligand-bound structures are shown in Figure 2—figure supplement 1. (B) Representative EM density around the ligands in this study. EM-maps of Ipx-, LY298-Ipx-, and VU154-Ipx were set to a contour level of 0.011 and the receptor-focused map of ACh- was set to 0.32. (C–E) Comparison of the receptor models with bound ligands and views from the (C) membrane, (D) extracellular surface, and (E) intracellular surface.

Table 2
Cryo-electron microscopy (cryo-EM) data collection, refinement, and validation statistics.
M4R-Gi1-IpxM4R-Gi1-Ipx-LY298M4R-Gi1-Ipx-VU154M4R-Gi1-ACh
Data collection & refinement
EMD code26,09926,10026,10126,102
Micrographs5056512160215913
Electron dose (e-/A2)666659.553.6
Voltage (kV)300300300300
Pixel size (Å)0.830.830.830.83
Spot size
Exposure time4435
Movie frames76767571
K3 CDS modeNoNoNoYes
Defocus range (µm)0.5–1.50.5–1.50.5–1.50.5–1.5
Symmetry imposedC1C1C1C1
Particles (final map)415,743617,793677,392315,595
Resolution @0.143 FSC (Å)2.82.42.52.8
Refinement
CCmap–model0.870.870.880.82
Map sharpening B factor (Å2)–80.9–60.8–46.6–85.1
Model quality
PDB code7TRK7TRP7TRQ7TRS
R.M.S. deviations
 Bond length (Å)0.0040.0040.0050.006
 Bond angles (o)0.8490.8110.8260.773
Ramachandran
 Favored (%)98.3899.1498.0298.10
 Outliers (%)0000
Rotamer outliers (%)0.110.2100
C-beta deviations (%)0000
Clashscore2.692.622.264.08
MolProbity score1.061.051.001.19
  1. mAChR: muscarinic acetylcholine receptor; ACh: acetylcholine; Ipx: iperoxo; FSC: Fourier shell correlation.

In all four structures, EM density beyond the top of transmembrane helix 1 (TM1) and the third intracellular loop (ICL3) of the receptor was poorly observed and not modeled. Similarly, the EM density of the α-helical domain of Gαi1 was poor and not modeled. These regions are highly dynamic and typically not modeled in many class A GPCR-G protein complex structures. Apart from these regions, most amino acid side chains were well resolved in the final EM density maps (Figure 2—figure supplement 3).

Structure and dynamics of agonist binding

Recently, cryo-EM structures of M4R-Gi1 complexes bound to Ipx, Ipx, and the PAM, LY2119620, and a putative novel allosteric agonist, c110, were determined (Wang et al., 2022). Surprisingly, comparison of the M4R-Gi1 complex structures revealed larger differences in the position of key orthosteric and allosteric site residues than the M1R-G11 and M2R-GoA complex structures (Figure 2—figure supplement 4). Unfortunately, the quality of density in the EM maps around the orthosteric and allosteric sites of these M4R-Gi1 structures (Wang et al., 2022) was poor, resulting in several key residues being mismodeled in each site (Figure 2—figure supplement 5). Therefore, differences between the M4R-Gi1 structures described herein and those by Wang et al., 2022 are highly likely to not be due to genuine differences and, as such, we compared the prior M1R-G11 and M2R-GoA complex structures (Maeda et al., 2019) in this study.

Overall, our M4R-Gi1 complex structures are similar in architecture to that of other activated class A GPCRs, including the M1R-G11 and M2R-GoA complexes (Figure 2—figure supplement 4). Superposition of the M4R-Gi1 complexes revealed nearly identical structures with root mean square deviations (RMSD) of 0.4–0.5 Å for the full complexes and 0.3–0.4 Å for the receptors alone (Figure 2C). The largest differences occur around the extracellular surface of the receptors (Figure 2D) along with slight displacements in the position of the αN helix of Gαi1 and Gβ1, Gγ2, and scFv16 with respect to the receptor (Figure 2—figure supplement 1D). The EM density of side chains surrounding the ACh and Ipx binding sites (Figure 3A and B) was well resolved providing the opportunity to understand structural determinants of orthosteric agonist binding. The orthosteric site of the M4 mAChR, in common with the other mAChR subtypes, is buried within the TM bundle in an aromatic cage that is composed of four tyrosine residues, two tryptophan residues, one phenylalanine residue, and seven other polar and nonpolar residues (Figure 3C). Notably, all 14 of these residues are absolutely conserved across all five mAChR subtypes, underscoring the difficulty in developing highly subtype-selective orthosteric agonists (Burger et al., 2018). Both ACh and Ipx have a positively charged trimethyl ammonium ion that makes cation-π interactions with Y1133.33, Y4166.51, Y4397.39, and Y4437.43 (Figure 3C; superscript refers to the Ballesteros and Weinstein scheme for conserved class A GPCR residues; Ballesteros and Weinstein, 1995). Likewise, both ACh and Ipx have a polar oxygen atom that can form a hydrogen bond to the indole nitrogen of W1644.57 with the oxygen of Ipx also being in position to interact with the backbone of N1173.37 (Figure 3D). Mutation of any of these contact residues reduces the affinity of ACh, validating their importance for agonist binding (Leach et al., 2011; Thal et al., 2016). The largest chemical difference between ACh and Ipx is the bulkier heterocyclic isoazoline group of Ipx that makes a π-π interaction with the conserved residue W4136.48 (Figure 3D). The residue W4136.48 is part of the CWxP motif, also known as the rotamer toggle switch, a residue that typically undergoes a change in rotamer between the inactive and active states of class A GPCRs (Shi et al., 2002).

Figure 3 with 1 supplement see all
Interactions of acetylcholine (ACh) and iperoxo (Ipx) with the receptor.

(A, B) Cryo-electron microscopy (cryo-EM) density of the (A) ACh- and (B) Ipx-bound structures. (C, D) Interactions at the orthosteric binding site comparing the active state ACh- and Ipx-bound structures with the inactive state tiotropium-bound structure (PDB: 5DSG). Arrows denote relative movement of residues between the inactive and active states. (D) Detailed interactions of ACh and Ipx. Hydrogen bonds are shown as black dashed lines. (E, F) Time courses from Gaussian accelerated molecular dynamics (GaMD) simulations of the ACh- and Ipx- bound M4R-Gi1 cryo-EM structures, each performed with three separate replicates. Individual replicate simulations are illustrated with different colors. The heading of each plot refers to the specific model used in the simulations. Root mean square deviations (RMSDs) of (E) ACh and (F) Ipx from simulations of the cryo-EM structures. (G, H) Cross-sections through the ACh- and Ipx-bound structures denoting the relative size of the binding pockets outlined in black.

To investigate the structural dynamics of the M4 mAChR, we performed three independent 500 ns GaMD simulations on the ACh- and Ipx-bound M4R-Gi1 cryo-EM structures (Table 3). GaMD simulations revealed that ACh undergoes higher fluctuations in the orthosteric site than Ipx (Figure 3E and F, Videos 1 and 2). Similarly, the interactions of N1173.37, W1644.57, and W4136.48 with Ipx were more stable than those with ACh (Figure 3—figure supplement 1). In the ACh-bound structure, W4136.48 was in a conformation matching the inactive-state tiotropium-bound structure (Figure 3C and D). GaMD simulations also showed that W4136.48 sampled a larger conformational space in the ACh-bound structure than in the Ipx-bound structure (Figure 3—figure supplement 1C and G). The predominate χ2 angle of W4136.48 was approximately 60 and 105 in the ACh-bound and Ipx-bound simulations, respectively, corresponding to the cryo-EM conformations.

Table 3
Gaussian accelerated molecular dynamics (GaMD) simulations of the M4 muscarinic acetylcholine receptor (mAChR).
SystemMethod
M4-Gi1-Ipx (cryo-EM structure)GaMD (3 × 500 ns)
M4-Gi1-Ipx-VU154 (cryo-EM structure)GaMD (3 × 500 ns)
M4-Gi1-Ipx-LY298 (cryo-EM structure)GaMD (3 × 500 ns)
M4-Gi1-ACh (cryo-EM structure)GaMD (3 × 500 ns)
M4-D432E-Gi1-Ipx-VU154GaMD (3 × 500 ns)
M4-T433R-Gi1-Ipx-VU154GaMD (3 × 500 ns)
M4-Gi1-ACh-VU154GaMD (3 × 500 ns)
M4-Gi1-ACh-LY298GaMD (3 × 500 ns)
M4-Gi1-VU154GaMD (3 × 500 ns)
M4-Gi1-LY298GaMD (3 × 500 ns)
M4-VU154GaMD (3 × 1000 ns)
M4-LY298GaMD (3 ×1000 ns)
Video 1
Movie from one Ipx-M4R-Gi1 Gaussian accelerated molecular dynamics (GaMD) simulation.
Video 2
Movie from one ACh-M4R-Gi1 Gaussian accelerated molecular dynamics (GaMD) simulation.

Located above ACh and Ipx is a tyrosine lid formed by three residues (Y1133.33, Y4166.51, and Y4397.39) that separate the orthosteric binding site from an extracellular vestibule (ECV) at the top of the receptor and the bulk solvent (Figure 3C). In the inactive conformation, the tyrosine lid is partially open due to Y4166.51 rotating away from the binding pocket to accommodate the binding of bulkier inverse agonists such as tiotropium. In contrast, mAChR agonists are typically smaller in size than antagonists and inverse agonists, and this is reflected in a contraction of the size of the orthosteric binding pocket from 115 Å3 when bound to tiotropium to 77 and 63 Å3 when bound to ACh and Ipx, respectively (Figure 3G and H; Tian et al., 2018). Together, the smaller binding pocket of Ipx and more stable binding interactions with nearby residues that include W4136.48 likely explain why Ipx has greater than 1000-fold higher binding affinity than ACh.

Structure and dynamics of PAM binding and allosteric modulation of agonist affinity

The M4R-Gi1 structures of LY298 and VU154 co-bound with Ipx are very similar to the Ipx- and ACh-bound structures, as well as to prior structures of the M2 mAChR bound to Ipx and the PAM, LY2119620 (Figure 2—figure supplement 4; Kruse et al., 2013; Maeda et al., 2019). Both LY298 and VU154 bind directly above the orthosteric site in the ECV that is composed of a floor delineated by the tyrosine lid, and ‘walls’ formed by residues from TM2, TM6, TM7, ECL2, and ECL3 (Figure 4A and B). The EM density surrounding the PAM binding site and the ECV of the M4 mAChR were clearly resolved with one exception; in the VU154-bound structure, the EM density begins to weaken around the trifluoromethylsulfonyl moiety (Figure 2B, Figure 4B). This was likely due to the moiety’s ability to freely rotate and a lack of strong interactions with the receptor.

Figure 4 with 3 supplements see all
Binding and dynamics of LY298 and VU154.

(A, B) Cryo-electron microscopy (cryo-EM) density of the (A) VU154- and (B) LY298-binding sites. (C) The root mean square deviations (RMSDs) between receptor models of the respective cryo-EM structures that were refined into the first and last frames of the EM maps from each principal component (PC1-PC3) of the 3D variability analysis. Values shown are mean ± SEM. (D, E) Top representative binding conformations of (D) VU154 and (E) LY298 obtained from structural clustering with frame populations ≥1% and time courses of the RMSDs of each positive allosteric modulator (PAM) relative to the cryo-EM structures. (F, G) Binding interactions of VU154 and LY298 with views from the (F) membrane and (G) extracellular surface. (H) Position and χ2 angle of W4357.35 in the tiotropium-, ACh-, Ipx-, VU154-Ipx-, and LY298-Ipx bound structures. (I–K) Time courses of the W4357.35 χ2 angle obtained from Gaussian accelerated molecular dynamics (GaMD) simulations on the (I) Ipx-, (J) VU154-Ipx-, and (K) LY298-Ipx-bound cryo-EM structures. See Table 3.

Given the overall similarities revealed by our four cryo-EM structures, we examined whether there were further differences in the dynamics between the PAM-bound structures by performing a 3D multivariance analysis (3DVA) of the principal components of motion within the Ipx-, LY298-Ipx, VU154-Ipx, and ACh-bound M4R-Gi1 cryo-EM data sets using Cryosparc (Punjani and Fleet, 2021); a similar analysis performed previously on cryo-EM structures of class A and class B GPCRs provided important insights into the allosteric motions of extracellular domains and receptor interactions with G proteins (Josephs et al., 2021; Liang et al., 2020; Mobbs et al., 2021; Zhang et al., 2020).

In the 3DVA of the Ipx-bound complex, the M4 mAChR appeared less flexible than the receptor in the ACh-bound complex (Videos 3 and 4) consistent with Ipx having a higher binding affinity and more stable pose during the GaMD simulations (Figure 3E and F). The LY298-Ipx-bound complex appeared similar to the Ipx-bound complex with LY298 being bound in the ECV (Video 5). In contrast, the 3DVA of the VU154 structure had more dynamic movements in the allosteric pocket that could reflect partial binding of VU154 (Video 6). This observation was in line with our findings that VU154 had lower binding modulation (Figure 1E) and functional modulation with agonists than LY298 (Figure 1J, Figure 1—figure supplement 2D, Table 1). To quantify the differences from the 3DVA, we rigid body fitted and refined the respective M4R-Gi1 models into the first and last frames of the EM maps from each principal component of the 3DVA and then calculated the RMSD between the receptor models (Figure 4C). In agreement with our prior observations, the VU154-Ipx-bound and ACh-bound complexes had greater RMSDs with values of 0.06 and 0.09 Å, respectively. Comparatively, the Ipx-bound and LY298-Ipx-bound complexes had lower RMSD values of 0.02 and 0.001 Å, respectively. The results of the 3DVA do not represent bona fide measures of receptor dynamics, rather they are suggestive of differences between the collected data sets that led to the structures. To support these findings, we compared the GaMD simulations of all four cryo-EM structures (Table 3). Notably, VU154 underwent considerably higher fluctuations than LY298 with RMSDs ranging from 1.5 to 15 Å for VU154 (Video 7) and 0.8–2.1 Å for LY298 (Video 8) relative to the cryo-EM structures (Figure 4D and E). Therefore, the GaMD simulations corroborate our 3DVA results and suggest that complexes bound to agonists with high affinity or co-bound with agonists and PAMs with high positive cooperativity will exhibit lower dynamic fluctuations.

Video 3
3D variability analysis of the Ipx-M4R-Gi1 cryo-electron microscopy (cryo-EM) structure.
Video 4
3D variability analysis of the ACh-M4R- Gi1 cryo-electron microscopy (cryo-EM) structure.
Video 5
3D variability analysis of the LY298-Ipx-M4R- Gi1 cryo-electron microscopy (cryo-EM) structure.
Video 6
3D variability analysis of the VU154-Ipx-M4R- Gi1 cryo-electron microscopy (cryo-EM) structure.
Video 7
Movie from one VU154-Ipx-M4R-Gi1 Gaussian accelerated molecular dynamics (GaMD) simulation.
Video 8
Movie from one LY298-Ipx-M4R-Gi1 Gaussian accelerated molecular dynamics (GaMD) simulation.

To investigate why the binding of LY298 was more stable than VU154, we examined the ligand interactions with the receptor. There are three key binding interactions that are shared between both PAMs and the M4 mAChR: (1) a three-way π-stacking interaction between F18645.51 (ECL2 residues have been numbered 45.X denoting their position between TM4 and TM5 with X.50 being a conserved cysteine residue), the aromatic core of the PAMs, and W4357.35; (2) a hydrogen bond between Y4397.39 of the tyrosine lid and the primary amine of the PAMs; and (3) a hydrogen bond between Y892.61 and the carbonyl oxygen of the PAMs (Figure 4F and G). While these interactions are conserved for both PAMs in the consensus cryo-EM maps, during GaMD simulations these interactions were more stable with LY298 than VU154 (Figure 4H–K, Figure 4—figure supplement 1). The importance of these interactions was validated pharmacologically (Figure 4—figure supplement 2, Table 4), whereby mutation of any of these residues completely abolished the binding affinity modulation mediated by LY298 and VU154 at the M4 mAChR with both Ipx and ACh as agonists.

Table 4
Pharmacological parameters of LY298 and VU154 at key M4 muscarinic acetylcholine receptor (mAChR) mutants.
[3H]-NMS saturation binding on stable M4 mAChR Flp-In CHO cells
ConstructsSites per cell*pKD
Human WT M4 mAChR (from Table 1)598,111 ± 43,067 (7)9.76 ± 0.05 (7)
Y89A2.6132,674 ± 4174 (4)9.88 ± 0.06 (4)
Q184A45.4988,728 ± 3056 (3)9.99 ± 0.06 (3)
F186A45.5136,907 ± 4170 (4)9.75 ± 0.16 (4)
W435A7.3534,861 ± 3510 (3)9.81 ± 0.22 (3)
Y439A7.3942,690 ± 4547 (3)8.31 ± 0.14 (3)
[3H]-NMS interaction binding assays between ACh or Ipx and LY298 or VU154 on stable M4 mAChR constructs in Flp-In CHO cells
ConstructsPAMpKi AChpKi IpxpKB PAMlog αACh§log αIpx§log αNMS
Human WT M4LY2985.09 ± 0.07 (7)8.54 ± 0.04 (11)= 5.651.57 ± 0.111.71 ± 0.09= 0
VU1545.06 ± 0.05 (7)8.54 ± 0.03 (11)= 5.831.44 ± 0.071.11 ± 0.06= 0
Y89A2.61LY2985.25 ± 0.05 (6)8.48 ± 0.05 (6)N.D.N.D.N.D.N.D.
VU1545.27 ± 0.05 (6)8.47 ± 0.05 (6)N.D.N.D.N.D.N.D.
Q184A45.49LY2985.24 ± 0.06 (6)8.74 ± 0.04 (10)6.23 ± 0.061.28 ± 0.131.27 ± 0.11–1.10 ± 0.07
VU1545.28 ± 0.05 (6)8.69 ± 0.04 (10)5.87 ± 0.171.07 ± 0.090.81 ± 0.07= 0
F186A45.51LY2984.91 ± 0.05 (6)8.12 ± 0.05 (8)N.D.N.D.N.D.N.D.
VU1544.91 ± 0.05 (6)8.12 ± 0.05 (8)N.D.N.D.N.D.N.D.
W4357.35LY2983.79 ± 0.07 (7)6.88 ± 0.07 (7)N.D.N.D.N.D.N.D.
VU1543.79 ± 0.07 (7)6.88 ± 0.07 (7)N.D.N.D.N.D.N.D.
Y439A7.39LY2983.23 ± 0.22 (8)5.36 ± 0.25 (8)N.D.N.D.N.D.N.D.
VU1543.23 ± 0.22 (8)5.36 ± 0.25 (8)N.D.N.D.N.D.N.D.
  1. Values represent the mean ± SEM with the number of independent experiments shown in parenthesis.

  2. N.D.: not determined; ACh: acetylcholine; Ipx: iperoxo; PAM: positive allosteric modulator.

  3. *

    Number of [3H]-NMS binding sites per cell.

  4. Negative logarithm of the radioligand equilibrium dissociation constant.

  5. Negative logarithm of the orthosteric (pKi) or allosteric (pKB) equilibrium dissociation constant. pKi values for ACh and Ipx are shared at each M4 mAChR construct. pKB values for the PAMs at Q184A are shared across the agonist data sets.

  6. §

    Logarithm of the binding cooperativity factor between the agonist (ACh or Ipx) and the PAM (LY298 or VU154).

  7. Logarithm of the binding cooperativity factor between the [3H]-NMS and the PAM (LY298 or VU154).

A potential fourth interaction was observed with residue Q18445.49 and the amide nitrogen of the PAMs; however, the GaMD simulations suggest that this interaction is relatively weak (Figure 4—figure supplement 1D and H), consistent with the fact that mutation of Q18445.49 to alanine had no effect on the binding affinity modulation of LY298 or VU154 (Figure 4—figure supplement 2, Table 4). In addition, each PAM had at least one potential unique binding interaction with the receptor (Figure 4F and G). For LY298, this is an interaction between the fluorine atom and N4236.58 that appeared to be stable during simulation and, when mutated to alanine reduced the binding modulation of LY298 (Figure 4—figure supplement 1K, Figure 4—figure supplement 2, Table 4; Thal et al., 2016). For VU154, there were two additional possible hydrogen bonding interactions with residues Y922.64 and T4337.33 (Figure 4G); however, these interactions were highly fluctuating during GaMD simulations, suggesting they were – at best – transient interactions (Figure 4—figure supplement 1I and J). Finally, W4357.35 is a key residue in the ECV that changes from a planar rotamer in the agonist-bound structures to a vertical rotamer that π stacks against the PAMs (Figure 4H). In GaMD simulations of the Ipx-bound structure, W4357.35 is predominantly in a planar conformation that corresponds to its conformation in the cryo-EM structure (Figure 4I). In contrast, the binding of LY298 stabilizes W4357.35 into a vertical position (Figure 4K). However, in the VU154-bound receptor, W4357.35 appears to alternate between the planar and vertical positions, consistent with VU154 having a less stable binding pose (Figure 4J). These results indicate that the binding of LY298 is more stable than VU154 due to LY298 being able to form stable binding interactions with key residues in the ECV. This provides a likely explanation for why LY298 was able to exert greater positive binding cooperativity on orthosteric agonists than VU154.

A molecular mechanism of probe dependence

As highlighted above, PAMs, LY298 and VU154, displayed stronger allosteric binding affinity modulation with ACh than Ipx, an example of probe dependence (Figure 1E, Table 1). These findings are in accord with previous studies where we identified probe dependence in the actions of LY298 when tested against other orthosteric agonists (Chan et al., 2008; Suratman et al., 2011). To investigate a mechanism for probe dependence at the M4 mAChR, we performed GaMD simulations with LY298 and VU154 co-bound with ACh by replacing Ipx with ACh in the corresponding cryo-EM structures (Table 3, Figure 4—figure supplement 3). In the absence of PAM, ACh was more dynamic than Ipx with root-mean-square fluctuations (RMSF) of 2.13 Å versus 0.88 Å, reflective of the fact Ipx binds with higher affinity than ACh (Figure 4—figure supplement 3F). In the presence of LY298 or VU154, the dynamics of ACh binding was decreased, with RMSFs reduced to 1.23 Å and 1.82 Å, respectively, and with LY298 having the greatest effect (Figure 4—figure supplement 3F). This is in line with LY298 having more cooperativity with ACh than VU154 (Figure 1E). In comparison to ACh, there was a modest increase in the dynamics of Ipx with the addition of LY298 or VU154, likely reflecting the fact Ipx binding to the receptor was already stable (Figure 3—figure supplement 1I and J, Figure 4—figure supplement 3F). These results provide a plausible mechanism for probe dependence, at least with regard to differences in the magnitude of the allosteric effect depending on the ligand bound. Namely, PAMs manifest higher cooperativity when interacting with agonists, such as ACh, that are inherently less stable on their own when bound to the receptor, in contrast to more stable ligands such as Ipx.

Structural and dynamic insights into orthosteric and allosteric agonism

In addition to the ability to allosterically modulate the function of orthosteric ligands, it has become increasingly appreciated that allosteric ligands may display variable degrees of direct agonism in their own right, over and above any allosteric modulatory effects (Changeux and Christopoulos, 2016). Prior studies have established that the activation process of GPCRs involves conformational changes that extend from the extracellular domains through to the intracellular surface (Nygaard et al., 2009). Comparison of the active state ACh-, Ipx-, LY298-Ipx-, and VU154-Ipx-bound M4R-Gi1 structures to the inactive state tiotropium-bound M4 mAChR structure (Protein Data Bank accession 5DSG) (Thal et al., 2016) thus affords an opportunity to gain new insights into the activation process mediated by multiple orthosteric agonists in the presence and absence of two different PAMs that display high (LY298) and low (VU154) degrees of direct allosteric agonism (Figures 1H and 5A–C).

Structural and dynamic insights into orthosteric and allosteric agonism.

(A) Cartoon of the receptor models indicating regions of interest for panels (B, C) shown within the red boxes. (B) View of the tiotropium-bound, agonist-bound, and positive allosteric modulator (PAM)-agonist-bound conformations from the extracellular surface. (C) Membrane view of residues and activation motifs involved in signaling. Residues colored red in (B, C) indicate residues of investigated in Gaussian accelerated molecular dynamics (GaMD) simulations. (D–G) Time course of the root mean square deviations (RMSDs) of the PAMs (D, E) from GaMD simulations of the M4R bound to G protein and no orthosteric agonist, (F, G) and in the absence of both G protein and agonist. (H–K) Similar to (D–G) the time courses of (H–K) the W4357.35 χ2 angle, (L–O) the W4136.48 χ2 angle, and (P–S) the TM3-TM6 distance measured by distance between R1303.50 and T3996.34. See Table 3.

As discussed previously, agonist binding decreases the size of the orthosteric binding site (Figure 3G and H). The primary driver of this decrease was the tyrosine lid residue Y4166.51, which underwent a large rotation toward Y1133.33 creating a hydrogen bond that seals off the tyrosine lid (Figure 3C). The closure of the tyrosine lid was further reinforced by a change in the rotamer of W4357.35 to a planar position that sits parallel to the tyrosine lid allowing for a π-π interaction with Y4166.51 and a positioning of the indole nitrogen of W4357.35 to potentially form a hydrogen bond with the hydroxyl of Y892.61 (Figure 5B). The contraction of the orthosteric pocket by the inward movement of Y4166.51 also led to a contraction of the ECV with a 5 Å inward movement of the top of TM6 and ECL3. As a consequence, the top of TM5 was displaced outward by 4 Å forming a new interface between TM5 and TM6 that was stabilized by a hydrogen bond between T4246.59 and the backbone nitrogen of P1935.36 along with aromatic interactions between F1975.40 and F4256.60 (Figure 5B). These interactions were specific to the active state structures and appear to be conserved as they were also present in the M1 and M2 mAChR active state structures (Maeda et al., 2019). In addition to the movements of TM5 and TM6, there was a smaller 1 Å inward movement of ECL2 (Figure 5B). The binding of LY298 and VU154 had a minimal impact on the conformation of most ECL residues, implying that the reorganization of residues in the ECV by orthosteric agonists contributes to the increased affinity of the PAMs (Figure 1G). There was a slight further inward shift of ECL2 toward the PAMs to facilitate the 3-way π-stacking interaction with F18645.51 and W4357.35. In addition, in the PAM-bound structures, Y892.61 rotated away from its position in the ACh- and Ipx-bound structures either due to a loss of an interaction with W4357.35 or to form a better hydrogen bond with the carbonyl oxygen of the PAMs (Figure 5B).

Below the orthosteric binding site are several signaling motifs that are important for the activation of class A GPCRs, including the PIF motif (Rasmussen et al., 2011; Wacker et al., 2013), the Na+ binding site (Liu et al., 2012a; White et al., 2018), the NPxxY motif (Fritze et al., 2003), and the DRY motif (Figure 5C; Ballesteros et al., 2001). The conformations of these activation motifs were very similar across all four active-state M4 mAChR structures and were consistent with the position of these motifs across other active-state class A GPCR structures (Zhou et al., 2019). Collectively, all of the described activation motifs facilitate an 11 Å outward movement of TM6 that typifies GPCR activation and creation of the G protein binding site. In comparison to the ECV residues (Figure 5B), beyond the rotamer toggle switch residue W4136.48, there are no discernible differences between the agonist and PAM-agonist-bound structures, suggesting a shared activation mechanism for residues below W4136.48 (Figure 5C).

As indicated above, LY298 also displays robust allosteric agonism in comparison to VU154 (Figure 1H, Figure 1—figure supplement 2B). To probe whether the allosteric agonism of LY298 could be related to its ability to better stabilize the M4 mAChR in an active conformation in comparison to VU154, we performed additional GaMD simulations on the LY298-Ipx- and VU154-Ipx-bound M4R-Gi1 structures with the agonist Ipx removed (3 × 500 ns) and with both Ipx and the G protein removed (3 × 1000 ns) (Figure 5D–S, Table 3). In GaMD simulations, LY298 underwent lower RMSD fluctuations than VU154 before dissociating from the receptor (Figure 5D–G). Similarly, the conformations of W4357.35 and W4136.48 were better stabilized in the LY298-Ipx-bound systems, indicating that LY298 more strongly promotes an active receptor conformation (Figure 5H–K). In the presence of the G protein, both PAMs stabilized an active conformation of the receptor based on the distances between TM3 and TM6 (Figure 5P and Q). Upon removal of the G protein, the VU154-bound M4 mAChR quickly transitioned toward the inactive conformation, while the LY298-bound M4 mAChR was more resistant to deactivation in the GaMD simulations (Figure 5R and S). This observation supports LY298 having greater efficacy than VU154 (Table 1) as it better stabilizes the active conformation of the M4 mAChR. Overall, the GaMD simulations show that in the absence of agonist alone, or agonist and G protein, LY298 better stabilizes activation motifs from the top of the receptor (W4357.35) all the way down to the intracellular G protein binding pocket (DRY-TM6), providing mechanistic insights into the function of LY298 as a stronger PAM-agonist than VU154.

Structural insights into allosteric modulation of agonist signaling

In a previous study, we characterized over 40 distinct mutations of M4 mAChR residues that span from the orthosteric site up to the extracellular surface (Table 5; Leach et al., 2011; Nawaratne et al., 2010; Thal et al., 2016). As expected, these studies revealed that mutation of residues around the orthosteric and allosteric sites often resulted in a reduction in the binding affinity of either ACh or LY298 at their respective binding sites, though the allosteric site was typically less affected (Figure 6A and B, Table 5). In contrast, the binding affinity modulation between ACh and LY298 was largely affected by mutation of aromatic residues that link the orthosteric and allosteric sites (Figure 6C), implying a network of residues that were responsible for transmitting binding cooperativity between these two sites (Thal et al., 2016). Analyzing unpublished data from prior studies allowed an examination of the signaling efficacy of ACh (τA) and LY298 (τB), but also the functional cooperativity (αβ) in the context of active state structures of the co-complexes (Figure 6D–F, Figure 6—figure supplement 1, Table 5).

Figure 6 with 2 supplements see all
Residues involved in binding, agonism, and modulation of acetylcholine (ACh) and LY298.

(A–F) M4 muscarinic acetylcholine receptor (mAChR) alanine point mutations that increase (green colored sticks) or decrease (pink colored sticks) (A) ACh binding, (B) LY298 binding, (C) binding modulation between ACh and LY298, (D) ACh efficacy, (E) LY298 efficacy, (F) and functional modulation by values more than tenfold. Efficacy values are corrected for receptor expression (Gregory et al., 2010) using receptor expression data from Thal et al., 2016. Quantitative data used to identify key residues are from both the current study and previous studies as summarized in Table 5 (Leach et al., 2011; Nawaratne et al., 2010; Thal et al., 2016). (G–I) pERK1/2 concentration response curves for interaction of ACh and LY298 at (G) WT and (H) W413A6.48 M4 mAChR with (I) values of efficacy and functional modulation. *Indicates statistical significance (p<0.05) relative to WT as determined by a one-way ANOVA with a Dunnett’s post-hoc test that includes the other M4 mAChR mutants. Data shown are mean ± SEM from three or more experiments performed in duplicate with the pharmacological parameters determined from a global fit of the data.

Table 5
Pharmacological parameters of M4 muscarinic acetylcholine receptor (mAChR) mutants.
pERK1/2 interaction assays*[3H]-QNB interaction binding assaysStudy
Constructslog τC ACh log τC LY298 log αβ §pKi ACh pKB LY298 log α**
WT M4 mAChR2.96 ± 0.14 (4)1.10 ± 0.092.43 ± 0.144.51 ± 0.154.89 ± 0.121.97 ± 0.11Current/Thal††
S85A2.573.15 ± 0.11(4)0.91 ± 0.071.75 ± 0.094.09 ± 0.115.44 ± 0.141.43 ± 0.06Current/Thal
Y89A2.612.53 ± 0.17 (5)= –3–0.43 ± 0.27*5.07 ± 0.455.36 ± 0.03–0.13 ± 0.08 ‡ ‡Current/Thal
Y92A2.642.26 ± 0.15 (4)–0.06 ± 0.16*2.25 ± 0.114.15 ± 0.254.53 ± 0.151.2 ± 0.19 ‡ ‡Current/Thal
I93T, I94V, K95I2.57 ± 0.112.27 ± 0.19 ‡ ‡N.T.4.69 ± 0.114.82 ± 0.362.14 ± 0.17 ‡ ‡Nawaratne § §
I93T2.652.34 ± 0.092.38 ± 0.22 ‡ ‡N.T.4.97 ± 0.045.36 ± 0.092.42 ± 0.16 ‡ ‡Nawaratne
I94V2.662.34 ± 0.091.24 ± 0.09N.T.4.71 ± 0.065.17 ± 0.081.74 ± 0.07Nawaratne
K95I2.672.00 ± 0.070.61 ± 0.09 ‡ ‡N.T.4.86 ± 0.055.20 ± 0.141.24 ± 0.04 ‡ ‡Nawaratne
Y97A23.492.94 ± 0.11 (4)2.19 ± 0.07*3.43 ± 0.12*4.69 ± 0.174.25 ± 0.10 ‡ ‡2.33 ± 0.12Current/Thal
W98A23.501.45 ± 0.15* (5)= –32.29 ± 0.103.65 ± 0.11 ‡ ‡4.39 ± 0.040.73 ± 0.07 ‡ ‡Current/Thal
G101A23.532.58 ± 0.09 (4)0.43 ± 0.08*2.10 ± 0.074.37 ± 0.195.03 ± 0.171.47 ± 0.01Current/Thal
D106A3.261.24 ± 0.11= –3N.T.3.95 ± 0.09 ‡ ‡5.29 ± 0.111.51 ± 0.15Leach ¶ ¶
W108A3.281.49 ± 0.17= –3N.T.4.01 ± 0.06 ‡ ‡4.24 ± 0.07 ‡ ‡1.23 ± 0.01 ‡ ‡Leach
L109A3.291.17 ± 0.14= –3N.T.3.11 ± 0.09 ‡ ‡4.28 ± 0.14 ‡ ‡2.54 ± 0.10 ‡ ‡Leach
D112E3.32–0.80 ± 0.16 ‡ ‡= –3N.T.<25.56 ± 0.130.39 ± 0.11 ‡ ‡Leach
D112N3.32N.D.N.D.N.T.3.19 ± 0.02 ‡ ‡5.79 ± 0.2 ‡ ‡0.74 ± 0.08Leach
Y113A3.33N.T.N.T.N.T.2.98 ± 0.12 ‡ ‡4.97 ± 0.150.80 ± 0.10 ‡ ‡Leach
S116A3.360.82 ± 0.17 ‡ ‡–0.35 ± 0.45 ‡ ‡N.T.3.61 ± 0.10 ‡ ‡5.12 ± 0.081.54 ± 0.05Leach
N117A3.370.80 ± 0.27 ‡ ‡–0.27 ± 0.16 ‡ ‡N.T.3.64 ± 0.04 ‡ ‡5.30 ± 0.151.57 ± 0.13Leach
V120A3.401.47 ± 0.111.20 ± 0.19N.T.5.63 ± 0.05 ‡ ‡5.41 ± 0.101.83 ± 0.11Leach
D129E3.491.45 ± 0.240.78 ± 0.16N.T.5.04 ± 0.075.59 ± 0.121.61 ± 0.16Leach
D129N3.492.56 ± 0.391.86 ± 0.12N.T.5.54 ± 0.10 ‡ ‡5.37 ± 0.201.86 ± 0.22Leach
W164A4.57N.D. (3)= –32.17 ± 0.64***3.95 ± 0.245.15 ± 0.28NDCurrent/Thal
F170A4.633.13 ± 0.17 (5)2.66 ± 0.12*3.58 ± 0.17*4.77 ± 0.24.53 ± 0.062.23 ± 0.13Current/Thal
W171A4.642.59 ± 0.17 (5)1.31 ± 0.113.11 ± 0.133.91 ± 0.214.56 ± 0.152.00 ± 0.09Current/Thal
Q172A4.653.05 ± 0.33 (5)1.18 ± 0.312.71 ± 0.164.02 ± 0.094.99 ± 0.031.54 ± 0.08Current/Thal
F173A4.663.39 ± 0.11(4)2.03 ± 0.10*3.31 ± 0.23*4.09 ± 0.014.78 ± 0.191.90 ± 0.14Current/Thal
Q184A45.493.01 ± 0.10 (4)1.05 ± 0.082.08 ± 0.124.25 ± 0.125.36 ± 0.041.70 ± 0.05Current/Thal
F186A45.511.99 ± 0.11N.D.N.T.4.85 ± 0.06NRNRNawaratne
I187A45.522.57 ± 0.100.62 ± 0.092.07 ± 0.153.71 ± 0.125.46 ± 0.291.07 ± 0.29 ‡ ‡Current/Thal
Q188A45.532.48 ± 0.150.99 ± 0.112.35 ± 0.164.6 ± 0.224.94 ± 0.081.49 ± 0.04Current/Thal
F189A45.542.28 ± 0.111.25 ± 0.092.67 ± 0.114.65 ± 0.025.09 ± 0.131.99 ± 0.08Current/Thal
L190A45.552.50 ± 0.141.32 ± 0.122.81 ± 0.124.20 ± 0.064.92 ± 0.12.06 ± 0.09Current/Thal
W413A6.480.66 ± 0.12***,* (4)0.61 ± 0.13***3.54 ± 0.09***,*3.47 ± 0.06 ‡ ‡4.51 ± 0.372.45 ± 0.36Current/Thal
Y416A6.51N.T.N.T.N.T.2.85 ± 0.10 ‡ ‡NRNRThal
N423A6.583.44 ± 0.15 (3)0.82 ± 0.101.43 ± 0.19*4.41 ± 0.155.02 ± 0.061.18 ± 0.08 ‡ ‡Current/Thal
Q427A6.623.15 ± 0.14 (3)0.99 ± 0.121.64 ± 0.124.46 ± 0.035.43 ± 0.061.36 ± 0.04Current/Thal
S428P6.631.99 ± 0.091.40 ± 0.19N.T.5.14 ± 0.03 ‡ ‡5.17 ± 0.151.81 ± 0.11Nawaratne
D432N7.322.26 ± 0.121.25 ± 0.18N.T.5.19 ± 0.04 ‡ ‡5.21 ± 0.21.37 ± 0.04Nawaratne
W435A7.352.58 ± 0.17 (4)= –3N.R3.37 ± 0.08 ‡ ‡NRNRCurrent/Thal
Y439A7.390.60 ± 0.18 ‡ ‡N.D.N.T.3.33 ± 0.10 ‡ ‡5.84 ± 0.120.49 ± 0.03 ‡ ‡Nawaratne
W440A7.403.69 ± 0.17†††(4)0.84 ± 0.111.52 ± 0.20†††4.29 ± 0.244.94 ± 0.060.96 ± 0.04 ‡ ‡Current/Thal
C442A7.421.49 ± 0.16 ‡ ‡0.82 ± 0.31N.T.4.04 ± 0.07 ‡ ‡5.35 ± 0.061.81 ± 0.03Nawaratne
Y443A7.430.50 ± 0.16 ‡ ‡N.D.N.T.3.36 ± 0.01 ‡ ‡6.22 ± 0.05 ‡ ‡1.16 ± 0.01 ‡ ‡Nawaratne
  1. Values represent the mean ± SEM from three or more independent experiments with the number of individual experimental replicates from the current study shown in parenthesis.

  2. N.T.: not tested; N.D.: not determined; N.R.: no response; ACh, acetylcholine.

  3. *

    Data and analysis from pERK1/2 assays were generated in the current study, Nawaratne et al, J. Bio. Chem. 2010, and Leach et al, Mol. Pharm. 2011. logτ ACh were calculated from the operational model of agonism. log τ LY298 and log αβ were calculated using a simplified operational model of allosterism.

  4. Data and analysis from [3H]-QNB interaction binding assays were generated in Nawaratne et al, J. Bio. Chem. 2010, Leach et al, Mol. Pharm., and Thal et al, Nature 2016.

  5. logτC = logarithm of the operational efficacy parameter corrected for receptor expression using the maximum number of receptor binding sites as previously determined from Nawaratne et al, J. Bio. Chem. 2010, Leach et al, Mol. Pharm., and Thal et al, Nature 2016.

  6. §

    Logarithm of the functional cooperativity factor between ACh and LY298.

  7. Negative logarithm of the orthosteric (pKi) or allosteric (pKB) equilibrium dissociation constant.

  8. **

    Logarithm of the binding cooperativity factor between ACh and LY298.

  9. ††

    Values of logτC ACh, log τC LY298, and log αβ that were calculated in this study. Other parameters are from Thal et al., 2016.

  10. ‡ ‡

    Values are significantly different from WT M4 mAChR as determined in previous studies.

  11. § §

    All values are from Nawaratne et al, J. Bio. Chem. 2010 with logτ corrected for receptor expression.

  12. ¶ ¶

    All values are from Leach et al, Mol. Pharm. 2011.

  13. ***

    Parameters determined from the full Operational Model of Allosterism.

  14. †††

    Values are significantly different from WT M4 mAChR (p<0.05) calculated by a one-way ANOVA with a Dunnett’s post-hoc test.

Mutation of residues that directly surround ACh primarily decreased the efficacy of ACh (Figure 6D, Table 5). One exception was W9823.50 (an ECL1 residue numbered 23.X denoting its position between TM2 and TM3 with X.50 denoting the most conserved residue), a residue that was recently identified in a deep scanning mutagenesis study as a conserved class A residue that is intolerant to mutation (Jones et al., 2020) and stabilizes the conserved disulfide bridge between ECL1 and TM3 that is important for the stability of the active state of many GPCRs including mAChRs (Hulme, 2013). Interestingly, residues that affect the efficacy of LY298 include nearly all of the residues that also affect ACh efficacy, along with residues that link to the allosteric site and surround the LY298 binding site (Figure 6E, Table 5). This suggests that the direct signaling of LY298 via the allosteric site is nonetheless linked through a similar network of residues and requires a functional orthosteric site for the transduction of signaling, and that mechanism involves equivalent closure of the orthosteric binding site, consistent with the thermodynamic reciprocity of cooperativity (Canals et al., 2011).

Residues Y892.61, N4326.58, W4357.35, and W4407.40 were identified as residues that, when mutated to alanine, significantly decreased the functional modulation between ACh and LY298 (Figure 6F, Table 5). In prior work, all four residues were also shown to contribute to LY298 binding or affinity modulation (Thal et al., 2016). Surprisingly, three mutations resulted in increased functional modulation by LY298. Of particular interest was, again, the rotamer toggle switch residue W4136.48. Mutation of W4136.48 to alanine significantly impaired the efficacy of ACh but only reduced the efficacy of LY298 by twofold, such that ACh and LY298 had similar efficacy for this mutant (Figure 6G–I, Table 5). Interestingly, the functional modulation (αβ) between ACh and LY298 increased to over 3600 (a 20-fold increase vs. WT) at W413A6.48. Similar results were observed in the TruPath assay with ACh, Ipx, and LY298 (Figure 6—figure supplement 2 [mutant], Figure 1—figure supplement 1B [WT]). However, with VU154, the functional modulation was considerably reduced with ACh and non-existent with Ipx, in line with our TruPath experiments at the WT M4 mAChR. These results show that, at the M4 mAChR, the rotamer toggle switch residue is important for the signaling efficacy of orthosteric agonists and PAM-agonists but does not impair the process of functional allosteric modulation. Thus, suggesting that the stability of LY298 co-binding with agonists can restore impaired function, while the less stable binding of VU154 does not. Together with the observation that most of the structural differences between the active-state M4 mAChR structures occur at or above W4136.48, we propose that this residue has a strong role in maintaining the conformational dynamics of the receptor and is a key trigger for robust signal transduction.

A molecular basis of species selectivity

One of the main advantages of allosteric modulators is the ability to selectivity target highly conserved proteins. The mAChRs are the prime example where allosteric modulators have been designed to selectively target specific subtypes. To date, the only PAM-bound mAChR structures are ones with LY2119620, a PAM that has activity at both the M2 and M4 mAChRs. Similarly, LY298 has activity at the M2 mAChR. However, the allosteric properties of VU154 are differentially affected by the species of the receptor (Wood et al., 2017b; Wood et al., 2017a). At the human M4 mAChR, LY298 displays robust binding affinity modulation, functional modulation, and allosteric agonism, while VU154 has comparatively weaker allosteric properties (Figure 1, Table 1). Conversely, at the mouse M4 mAChR, VU154 has a high degree of positive binding modulation, functional modulation, and allosteric agonism that is comparable to LY298 at the human M4 mAChR (Figure 7—figure supplements 1 and 2, Table 1). Therefore, we aimed to determine whether our prior findings could be used to explain the selectivity of VU154 between the human and mouse receptors.

The amino acid sequences of the human and mouse M4 mAChRs are highly conserved, with most of the differences occurring between the long third intracellular loop and the N- and C- termini. As shown in Figure 7A, only three residues differ between the human and mouse M4 mAChR with respect to the transmembrane domain. Specifically, residue V91 (L in mouse) at the top of TM2 points into the lipid bilayer, and D432 and T433 (E and R in mouse), which are located at the top of TM7 and form part of the allosteric binding site near VU154.

Figure 7 with 3 supplements see all
A molecular mechanism for the species selectivity for VU154.

(A) Comparison of the cryo-electron microscopy (cryo-EM) structure of the human M4 muscarinic acetylcholine receptor (mAChR) bound to Ipx-VU154 with the AlphaFold model of the mouse M4 mAChR (Jumper et al., 2021; Varadi et al., 2022). The three residues that differ between species and within the core 7TM bundle from the human receptor (V91, D432, and T433) are shown as sticks along with the corresponding residues from the mouse receptor. (B) The binding affinity of VU154 for the Ipx-bound conformation (pKB-Ipx = pKB + α) determined from [3H]-NMS binding experiments. Values calculated with data from Figure 7—figure supplement 1 with propagated error. (C) Efficacy of VU154 (τB – corrected for receptor expression) of pERK1/2 signaling from data in Figure 7—figure supplement 2. (D–K) Time courses of obtained from Gaussian accelerated molecular dynamics (GaMD) simulations of the (D–G) D432E and (H–K) T433R mutant M4R-Ipx-Gi1-VU154 systems with (D, H) Ipx RMSDs, (E, I), VU154 root mean square deviations (RMSDs), (F, J) W4357.35 χ2 angle, and (G, K) W4136.48 χ2 angle. Data shown are mean ± SEM from three or more experiments performed in duplicate with the pharmacological parameters determined from a global fit of the data. *Indicates statistical significance (p<0.05) relative to WT as determined by a one-way ANOVA with a Dunnett’s post-hoc test.

Previous work suggested that residues D432 and T433 were important for differences in the species selectivity of LY298 (Chan et al., 2008). As such, we examined two single D432E and T433R mutants and a V91L/D432E/T433R triple mutant of the human receptor, along with the mouse M4 mAChR in radioligand binding and pERK1/2 experiments using Ipx and both PAMs (Figure 7—figure supplements 1 and 2, Table 1). For LY298, there were no statistically significant differences in binding or function between species and across the mutants that were more than threefold in effect. In contrast, VU154 had a tenfold higher binding affinity for the Ipx-bound mouse M4 mAChR (compare Figure 1G with Figure 7B). The affinity of VU154 increased by 2.5-fold at the D432E and T433R mutants and the triple mutant matched the affinity of the mouse receptor (Figure 7B). In functional assays, similar results were observed for VU154 with Ipx at the mouse M4 mAChR, with significant increases in the efficacy (τB – corrected for receptor expression), transduction coefficients (τB/KB), and functional modulation (αβ) (Figure 7B, Figure 7—figure supplements 1 and 2, Table 1). Relative to the WT M4 mAChR, the efficacy (Figure 7C), transduction coefficients, and functional modulation of VU154 increased for all of the mutants (Figure 7—figure supplements 1 and 2, Table 1); however, none of the values fully matched the mouse receptor. Nevertheless, these results indicate that V91L, D432E, and T433R play a key role in mediating the species selectivity of VU154.

Our prior findings suggest the robust allosteric activity of LY298 at the human M4 mAChR was due to stable interactions with the receptor. As a proof-of-principle, we questioned whether GaMD simulations would produce a stable binding mode for VU154 with D432E and T433R mutations to the VU154-Ipx-bound M4R-Gi1 cryo-EM structure that was similar to our previously observed stable binding pose of LY298 (Figure 4). Excitingly, both the D432E and T433R mutants resulted in a dynamic profile of VU154 that matched our GaMD simulations of LY298 from the LY298-Ipx-bound M4R-Gi1 cryo-EM structure, including stabilized VU154 binding, constrained χ2 rotamer conformations of W4357.35 and W4136.48, and stable binding interactions with Y892.61, Y4397.39, Q18445.49, and F18645.51 (Figure 7D–K, Figure 7—figure supplement 3, Videos 9 and 10). The GaMD simulations also suggest that a potential interaction between the mutant residue T433R and the sulfoxide group of VU154 was more stable (5.2 ± 1.5 Å; Figure 7—figure supplement 3I) versus the WT residue T433 (6.56 ± 2.1 Å, Figure 4—figure supplement 1J), albeit the distance of this interaction was far apart and would be better validated by structure determination of VU154 with the mouse M4 mAChR.

Video 9
Movie from one VU154-Ipx-M4R(D432E)-Gi1 Gaussian accelerated molecular dynamics (GaMD) simulation.
Video 10
Movie from one VU154-Ipx-M4R(T433R)-Gi1 Gaussian accelerated molecular dynamics (GaMD) simulation.

Collectively, these findings reiterate the importance of receptor dynamics in the determination of allosteric modulator selectivity as even subtle differences in amino acid residues between species may result in profound changes in overall stability of the same PAM-agonist-receptor complex.

Discussion

Major advances have been made in recent years in the appreciation of the role of GPCR allostery and its relevance to modern drug discovery (Changeux and Christopoulos, 2016; Wootten et al., 2013). Despite an increase in the number of reported high-resolution GPCR structures bound to allosteric ligands (Thal et al., 2018), there remains a paucity of molecular-level details about the interplay between the complex chemical and pharmacological parameters that define allostery at GPCRs. By combining detailed pharmacology studies, multiple high-resolution cryo-EM structures of the M4 mAChR bound to two pharmacologically different agonists and PAMs, and GaMD simulations, we have now provided exquisite in-depth insights into the relationship between both structure and dynamics that govern multiple facets of GPCR allostery (Figure 8A).

Conformational dynamics of the allostery at M4 muscarinic acetylcholine receptor (mAChR) signaling complexes.

(A) A schematic cartoon illustrating the conformational states of the ligands and the M4 mAChR when bound to different types of ligands and transducer, along with the resulting dynamic profiles. Pharmacological parameters related to each conformational change are shown. Stable ligand–receptor interactions are denoted by a straight line and less-stable (more dynamic) interactions are denoted by a wavy line. (B) Iperoxo (Ipx) bound the M4 mAChR with a higher affinity and more stability than ACh but had lower efficacy. ACh being more loosely bound and coupled to G protein may facilitate more G protein turnover accounting for its higher efficacy. (C) LY298 and VU154 bound to the M4 mAChR with similar affinity for the receptor, but LY298 was found to bind more stably. LY298 had a higher efficacy than VU154, suggesting that allosteric agonism at the M4 mAChR is mediated by stabilization of the extracellular vestibule (ECV). (D) The positive allosteric modulators (PAMs) LY298 and VU154 display robust binding modulation at the M4 mAChR with LY298 having a stronger allosteric effect. Both PAMs displayed stronger binding modulation with the agonist ACh versus Ipx, an example of probe dependence. Both PAMs also displayed a slight negative to neutral effect on the efficacy of the agonists, suggesting that their mechanism of action is largely through binding.

Comparison of the ACh- and Ipx-bound M4 mAChR structures revealed that Ipx bound in a smaller binding pocket (Figure 3G and H), and GaMD simulations showed that Ipx formed more stable interactions with the receptor (Figure 3—figure supplement 1). These observations likely explained why Ipx exhibited greater than 1000-fold higher binding affinity than ACh (Figure 1D), being consistent with studies of other agonists at the β1-adrenoceptor and the M1 mAChR (Brown et al., 2021; Warne et al., 2019; Figure 8B). The observation that ACh was a more efficacious agonist than Ipx (Table 1) yet bound with lower affinity and less stable interactions than Ipx was paradoxical. Kenakin and Onaran, 2002 previously opined on the paradox between ligand binding affinity and efficacy and showed via simulations that, in general, there was a negative correlation between binding affinity and efficacy. One interpretation of these results was that the ACh-bound M4 mAChR more readily sampled receptor conformations that engaged with the transducers (Manglik et al., 2015). Similarly, the ACh-bound M4 mAChR may also have faster G protein turnover than Ipx due to Ipx-M4R-Gi1 forming a more stable ternary complex (Furness et al., 2016; Figure 8B).

It is worth noting that structures of GPCRs bound to agonists with different pharmacological properties (full, partial, and biased agonists) have now been reported for some GPCRs (Liang et al., 2018a; Masureel et al., 2018; McCorvy et al., 2018; Ring et al., 2013; Wacker et al., 2013; Warne et al., 2012; Wingler et al., 2019). However, insights gained from such cryo-EM and X-ray crystallography structures may be limited due to the role that the bound transducer plays on the observed final receptor conformation, and not necessarily due solely to the properties of the ligand. The ultimate underlying conformational differences, therefore, are likely to be subtle and dynamic (Seyedabadi et al., 2022), requiring application of additional techniques such as NMR spectroscopy, single-molecule FRET and MD simulations for furthering our understanding (Cao et al., 2021; Cong et al., 2021; Gregorio et al., 2017; Huang et al., 2021; Katayama et al., 2021; Liu et al., 2012b; Solt et al., 2017; Sušac et al., 2018; Xu et al., 2023; Ye et al., 2016).

Indeed, if considering this issue from the perspective of allosteric modulators of GPCRs, our study highlights that two PAMs with distinctly different pharmacological profiles (Figure 1) may bind to and stabilize receptor conformations that were very similar when viewed as static structures (Figure 4). Yet, in contrast, the 3DVA analysis from our cryo-EM structures suggested differences in the dynamics of the cryo-EM structures that were explored further in GaMD simulations (Figure 4C) and revealed that LY298 had a more stable binding pose and interactions with the receptor than VU154 in the PAM-agonist–receptor–transducer-bound conformation. These observations were consistent with LY298 having greater positive binding cooperativity than VU154 (Figure 1E) and suggest that GaMD simulations of GPCRs bound to allosteric ligands could be an extremely valuable tool for drug discovery and optimization (Bhattarai and Miao, 2018).

Pharmacological analysis revealed that LY298 is a better PAM-agonist than VU154 with respect to efficacy (Figure 1H) in the Gi1 TruPath and pERK1/2 signaling assays (Figure 1—figure supplement 2B). GaMD simulations of the PAM–receptor–transducer and PAM–receptor bound complexes, again showed that LY298 more stably interacted with the receptor (Figure 4) and in the absence of G protein better stabilized the duration of the active conformation of the receptor (Figure 5). These findings were not contradictory to our above findings that ACh was more efficacious than Ipx despite having weaker interactions with the receptor because when the affinity of the ligands was accounted for in the transduction coupling coefficients, the rank order was Ipx >> ACh ~ LY298 > VU154 (Figure 1I). Furthermore, these results were in accordance with the observations of Kenakin and Onaran that ligands with the same binding affinity can also have differing efficacies (and vice versa). In addition, the mechanism of agonism for allosteric ligands that bind to the ECV may differ (Xu et al., 2021). Prior work by DeVree et al., 2016 established that allosteric coupling of G proteins to the unliganded active receptor conformation promoted closure of the ECV region. This allosteric coupling is reciprocal and stabilizing the ECV region by PAMs likely leads to increased efficacy (Figure 8).

The PAMs, LY298 and VU154, also displayed stronger allosteric effects with ACh than with Ipx, an observation known as probe dependence (Figure 1E–G). Probe dependence can have substantial implications on how allosteric ligands are detected, validated, and their potential therapeutic utility (Kenakin, 2005). Examples of probe dependence are not limited to studies on mAChRs and have been observed across multiple receptor families (Christopoulos, 2014; Gentry et al., 2015; Pani et al., 2021; Slosky et al., 2020; Wang et al., 2021b). GaMD simulations comparing the PAMs co-bound with either Ipx or ACh showed that the PAMs had a stabilizing effect on ACh, whereas the stability of Ipx was slightly reduced by the PAMs likely because the binding of Ipx was already stable. This is a sensible explanation from thermodynamic principles. Another explanation invokes the two-state receptor model (Canals et al., 2011), which stipulates that the degree of positive modulation for PAMs increases with an increase in the efficacy of the agonists. The pharmacology data support this model as ACh was more efficacious than Ipx and was better modulated by both PAMs (Figure 8D). These observations are also consistent with recent studies that suggest that conformational dynamics between agonist and receptor are important for functional signaling (Bumbak et al., 2020; Cary et al., 2022; Deganutti et al., 2022; O’Connor et al., 2015).

The findings presented here provide new insights into the allosteric signaling and allosteric modulation of GPCRs by combining the analytical analysis of multiple pharmacology assays with cryo-EM structures and GaMD simulations. Overall, these results provide a framework for future mechanistic studies and, ultimately, can aid in the discovery, design, and optimization of allosteric drugs as novel therapeutic candidates for clinical progression.

Limitations of the study

The complexities of GPCR signaling cannot be fully explained by any single receptor or set of experiments. This study was limited to the investigation of two agonists and two PAMs at the human M4 mAChR. Future studies will be required to determine how these results extrapolate to other classes of ligand, mAChR subtypes, and GPCRs. For instance, this study determined the structures of the M4 mAChR bound with the ligands ACh, Ipx, Ipx-LY298, and Ipx-VU154. It is possible that structures of the M4 mAChR bound with ACh-LY298 and ACh-VU154 could reveal different receptor conformations (although GaMD simulations already performed on their docked complexes and the conformational differences between the Ipx-bound cryo-EM structures suggest otherwise). Similarly, structures of the M4 mAChR bound in complex with either PAM alone may provide better insights into direct allosteric agonism. However, we note that our attempt at determining an LY298-bound complex did not have sufficient stability for the determination of a high-resolution structure, as also supported by our GaMD simulations. Additionally, our cryo-EM structures and MD-simulations utilized an M4 mAChR sequence with a large portion of the third intracellular loop removed and were complexed with a dominant negative mutant of Gαi1 and stabilized with the antibody scFv16. This contrasts with our pharmacological characterization of the ligands that were performed on the WT M4 mAChR. Further investigation into the molecular determinants of species selectivity is also warranted, as is the need for future experiments that incorporate the combined interplay between dynamics/kinetics of ligands, receptor, transducer recruitment and activation.

Materials and methods

Bacterial strains

Request a detailed protocol

DH5α (New England Biolabs) and DH10bac (Thermo Fisher Scientific) Escherichia coli cells were grown in LB at 37°C.

Cell culture

Request a detailed protocol

Tni and Sf9 cells (Expression Systems) were maintained in ESF-921 media (Expression Systems) at 27°C. Flp-In Chinese hamster ovary (CHO) (Thermo Fisher Scientific) cells stably expressing human M4 mAChR or mutant constructs were maintained in Dulbecco’s modified Eagle’s medium (DMEM, Invitrogen) containing 5% fetal bovine serum (FBS; ThermoTrace) and 0.6 μg/ml of Hygromycin (Roche) in a humidified incubator (37°C, 5% CO2, 95% O2). HEK293A cells were grown in DMEM supplemented with 5% FBS at 37°C in 5% CO2. Cell lines were authenticated by vendor and confirmed negative for mycoplasma contamination using the Lonza MycoAlert Mycoplasma Detection Kit (#LT07-318).

Radioligand binding assays

Request a detailed protocol

Flp-In CHO cells stably expressing M4 mAChR constructs were seeded at 10,000 cells/well in 96-well white clear bottom isoplates (Greiner Bio-one) and allowed to adhere overnight at 37°C, 5% CO2, and 95% O2. Saturation binding assay was performed to quantify the receptor expression and equilibrium dissociation constant of the radioligand [3H]-NMS (PerkinElmer, specific activity 80 Ci/mmol). Briefly, plates were washed once with phosphate-buffered saline (PBS) and incubated overnight at room temperature (RT) with 0.01–10 nM [3H]-NMS in Hanks’s balanced salt solution (HBSS)/10 mM HEPES (pH 7.4) in a final volume of 100 μl. For binding interaction assays, cells were incubated overnight at RT with a specific concentration of [3H]-NMS (pKD determined at each receptor in saturation binding) and various concentrations of ACh or Ipx in the absence or presence of increasing concentrations of each allosteric modulator. In all cases, nonspecific binding was determined by the coaddition of 10 μM atropine (Sigma). The following day, the assays were terminated by washing the plates twice with ice-cold 0.9% NaCl to remove the unbound radioligand. Cells were solubilized in 100 μl per well of Ultima Gold (PerkinElmer), and radioactivity was measured with a MicroBeta plate reader (PerkinElmer).

G protein activation assay

Request a detailed protocol

Upon 60–80% confluence, HEK293A cells were transfected transiently using polyethylenimine (PEI, Polysciences) and 10 ng per well of each of pcDNA3.1-hM4 mAChR (WT or mutant), pcDNA5/FRT/TO-Gαi1-RLuc8, pcDNA3.1-β3, and pcDNA3.1-Gγ9-GFP2 at a ratio of 1:1:1:1 ratio with 40 ng of total DNA per well. Cells were plated at 30,000 cells per well into 96-well Greiner CELLSTAR white-walled plates (Sigma-Aldrich). 48 hr later, cells were washed with 200 μl phosphate buffer saline (PBS) and replaced with 70 μL of 1× HBSS with 10 mM HEPES. Cells were incubated for 30 min at 37°C before addition of 10 μl of 1.3 μM Prolume Purple coelenterazine (Nanolight Technology). Cells were further incubated for 10 min at 37C° before BRET measurements were performed on a PHERAstar plate reader (BMG Labtech) using 410/80 nm and 515/30 nm filters. Baseline measurements were taken for 8 min before addition of drugs or vehicle to give a final assay volume of 100 μl and further reading for 30 min. BRET signal was calculated as the ratio of 515/30 nm emission over 410/80 nm emission. The ratio was vehicle corrected using the initial 8 min of baseline measurements and then baseline corrected again using the vehicle-treated wells. Data were normalized using the maximum agonist response to allow for grouping of results using an area under the curve analysis in Prism. Data were analyzed at timepoints of 4, 10, and 30 min yielding similar results.

Phospho-ERK1/2 assay

Request a detailed protocol

The level of phosphorylated extracellular signal-regulated protein kinase 1/2 (pERK1/2) was detected using the AlphaScreen SureFire Kit (PerkinElmer Life and Analytical Sciences). Briefly, FlpIn CHO cells stably expressing the receptor were seeded into transparent 96-well plates at a density of 20,000 cells/well and grown overnight at 37°C, 5% CO2. Cells were washed with PBS and incubated in serum-free DMEM at 37°C for 4 hr to allow FBS-stimulated pERK1/2 levels to subside. Cells were stimulated with increasing concentrations of ACh or Ipx in the absence or presence of increasing concentrations of the allosteric modulator at 37°C for 5 min (the time required to maximally promote ERK phosphorylation for each ligand at each M4 mAChR construct in the initial time-course study; data not shown). For all experiments, stimulation with 10% (v/v) FBS for 5 min was used as a positive control. The reaction was terminated by the removal of media and lysis of cells with 50 μl of the SureFire lysis buffer (TGR Biosciences). Plates were then agitated for 5 min and 5 μl of the cell lysate was transferred to a white 384-well ProxiPlate (Greiner Bio-one) followed by the addition of 5 μl of the detection buffer (a mixture of activation buffer:reaction buffer:acceptor beads:donor beads at a ratio of 50:200:1:1). Plates were incubated in the dark for 1 hr at 37°C followed by measurement of fluorescence using an Envision plate reader (PerkinElmer) with standard AlphaScreen settings. Data were normalized to the maximal response mediated by 10 μM ACh, Ipx, or 10% FBS.

Purification of scFv16

Request a detailed protocol

Tni insect cells were infected with scFv16 baculovirus at a density of 4 million cells per ml and harvested at 60 hr post infection by centrifugation for 10 min at 10,000 × g. The supernatant was pH balanced to pH 7.5 by the addition of Tris pH 7.5, and 5 mM CaCl2 was added to quench any chelating agents, then left to stir for 1.5 hr at RT. The supernatant was then centrifuged at 30,000 × g for 15 min to remove any precipitates. 5 ml of EDTA-resistant Ni resin (Cytivia) was added and incubated for 2 hr at 4oC while stirring. Resin was collected in a glass column and washed with 20 column volumes (CVs) of high salt buffer (20 mM HEPES pH 7.5, 500 mM NaCl, 20 mM imidazole) followed by 20 CVs of low salt buffer (20 mM HEPES pH 7.5, 100 mM NaCl, 20 mM imidazole). Protein was then eluted using 8 CV of elution buffer (20 mM HEPES pH 7.5, 100 mM NaCl, 250 mM imidazole) until no more protein was detected using Bradford reagent (Bio-Rad Laboratories). Protein was concentrated using a 10 kDa Amicon filter device (Millipore) and aliquoted into 1 mg aliquots for further use.

Expression and purification of M4R-Gi1-scFv16 complexes

Request a detailed protocol

The human M4 mAChR with residues 242–387 of the third intracellular loop removed and the N-terminal glycosylation sites (N3, N9, N13) mutated to D was expressed in Sf9 insect cells, and human DNGαi1 and His6-tagged human Gβ1γ2 were co-expressed in Tni insect cells. Cell cultures were grown to a density of 4 million cell per ml for Sf9 cells and 3.6 million per ml for Tni cells and then infected with either M4 mAChR baculovirus or both Gαi1 and Gβ1γ2 baculovirus, at a ratio of 1:1. M4 mAChR expression was supplemented with 10 mM atropine. Cultures were grown at 27°C and harvested by centrifugation 60–72 hr (48 hr for Hi5 cells) post infection. Cells were frozen and stored at –80°C for later use. 1–2 l of the frozen cells were used for each purification.

Cells expressing M4 mAChR were thawed at RT and then dounced in the solubilization buffer containing 20 mM HEPES pH 7.5, 10% glycerol, 750 mM NaCl, 5 mM MgCl2, 5 mM CaCl2, 0.5% LMNG, 0.02% CHS, 10 µM atropine, and cOmplete Protease Inhibitor Cocktail (Roche) until homogeneous. The receptor was solubilized for 2 hr at 4°C while stirring. The insoluble material was removed by centrifugation at 30,000 × g for 30 min followed by filtering the supernatant and batch-binding immobilization to M1 anti-flag affinity resin, previously equilibrated with high salt buffer, for 1 hr at RT. The resin with immobilized receptor was then washed using a peristaltic pump for 30 min at 2 ml/min with high salt buffer: 20 mM HEPES pH 7.5, 750 mM NaCl, 5 mM MgCl2, 5 mM CaCl2, 0.5% lauryl maltose neopentyl glycol (LMNG, Anatrace), 0.02% cholesterol hemisuccinate (CHS, Anatrace) followed by low salt buffer: 20 mM HEPES pH 7.5, 100 mM NaCl, 5 mM MgCl2, 5 mM CaCl2, 0.5% LMNG, 0.02% CHS, and an agonist (5 µM Ipx, 1 µM Ipx with 10 µM VU154, or 100 µM ACh). While the receptor was immobilized on anti-FLAG resin, the DNGαi1 cell pellet was thawed, dounced, and solubilized in the solubilization buffer containing 20 mM HEPES pH 7.5, 100 mM NaCl, 5 mM MgCl2, 5 mM CaCl2, 0.5% LMNG, 0.02% CHS, apyrase (five units), and cOmplete Protease Inhibitor Cocktail. DNGαi1 was solubilized for 2 hr at 4°C followed by the centrifugation at 30,000 × g for 30 min to remove the insoluble material. Supernatant was filtered through a glass fiber filter (Millipore) and then added to the receptor bound to anti-Flag resin. Apyrase (five units), scFv16, and agonist (either 1 µM Ipx, 1 µM Ipx with 10 µM VU154, or 100 µM ACh) were added and incubated for 1 hr at RT with gentle mixing. The anti-FLAG resin was then loaded onto a glass column and washed with approximately 20 CVs of washing buffer: 20 mM HEPES pH 7.4, 100 mM NaCl, 5 mM MgCl2, 5 mM CaCl2, 0.01% LMNG, 0.001% CHS, agonist (1 µM Ipx, 1 µM Ipx with 10 µM VU154, or 100 µM ACh). Complex was eluted with size-exclusion chromatography (SEC) buffer: 20 mM HEPES pH 7.5, 100 mM NaCl, 5 mM MgCl2, 0.01% LMNG, 0.001% CHS and agonist (1 µM Ipx, or 1 µM Ipx with 10 µM VU154, or 100 µM ACh) with the addition of 10 mM EGTA and 0.1 mg/mL FLAG peptide. After the elution, an additional 1–2 mg of scFv16 was added and shortly incubated on ice before concentrating using a 100 kDa Amicon filter to a final volume of 500 µl. The sample was filtered using a 0.22 µm filter followed by SEC using a Superdex 200 increase 10/300 column (Cytivia) using SEC buffer. For the ACh- and VU154-Ipx-bound samples, the fractions containing protein were concentrated again and re-run over SEC using a buffer with half the amount of detergent in order to remove empty micelles. Samples were concentrated and flash frozen using liquid nitrogen. In case of the LY298-Ipx-bound sample, the sample was purified with 1 µM Ipx only. After SEC, the sample was then split in half, where one half was incubated with approximately 1.6 µM LY298 at 4°C overnight, and then concentrated and flash frozen in liquid nitrogen.

EM sample preparation and data acquisition

Request a detailed protocol

Samples (3 µl) were applied to glow-discharged Quantifoil R1.2/1.3 Cu/Rh 200 mesh grids (Quantifoil) (M4R-Gi1-Ipx and M4R-Gi1-Ipx-LY298) or UltrAuFoil R1.2/1.3 Au 300 mesh grids (Quantifoil) (M4R-Gi1-Ipx-VU154 and M4R-Gi1-Ach) and were vitrified on a Vitrobot Mark IV (Thermo Fisher Scientific) set to 4°C and 100% humidity and 10 s blot time. Data were collected on a Titan Krios G3i 300 kV electron microscope (Thermo Fisher Scientific) equipped with GIF Quantum energy filter and K3 detector (Gatan). Data acquisition was performed in EFTEM NanoProbe mode with a 50 µM C2 aperture at an indicated magnification of ×105,000 with zero-loss slit width of 25 eV. The data were collected automatically with homemade scripts for SerialEM performing a nine-hole beam-image shift acquisition scheme with one exposure in the center of each hole. Experimental parameters specific to each collected data set is listed in Table 2.

Image processing

Request a detailed protocol

Specific details for the processing of each cryo-EM data set are shown in Figure 2—figure supplement 2. Image frames for each movie were motion corrected using MotionCor2 (Zheng et al., 2017) and contrast transfer function (CTF)-estimated using GCTF (Zhang, 2016). Particles were picked from corrected micrographs using crYOLO (Wagner et al., 2019) or RELION-3.1 software Zivanov et al., 2018 followed by reference-free 2D and 3D classifications. Particles within bad classes were removed and remaining particles subjected to further analysis. Resulting particles were subjected to Bayesian polishing, CTF refinement, 3D auto-refinement in RELION, followed by another round of 3D classification and 3D refinement that yielded the final maps (Zivanov et al., 2018). Local resolution was determined from RELION using half-reconstructions as input maps. Due to the high degree of conformational flexibility between the receptor and G protein, a further local refinement was performed in cryoSPARC for the ACh-bound M4R-complex. A receptor-focused map was generated (2.75 Å), which was used to generate a PDB model of the ACh-bound M4R.

Model building and refinement

Request a detailed protocol

An initial M4R template model was generated from our prior modeling studies of the M4 mAChR that was based on an active state M2 mAChR structure (PBD: 4MQT) (Kruse et al., 2013). An initial model for dominant negative Gαi112 was from a structure in complex with Smoothend (PDB: 6OT0) (Qi et al., 2019) and scFv16 from the X-ray crystal structure in complex with heterotrimeric G protein (PDB: 6CRK) (Maeda et al., 2018). Models were fit into EM maps using UCSF Chimera (Pettersen et al., 2004), and then rigid-body-fit using PHENIX (Liebschner et al., 2019), followed by iterative rounds of model rebuilding in Coot (Casañal et al., 2020) and ISOLDE (Croll, 2018), and real-space refinement in PHENIX. Restrains for all ligands were generated from the GRADE server (https://grade.globalphasing.org). Model validation was performed with MolProbity (Williams et al., 2018) and the wwPDB validation server (Berman et al., 2003). Figures were generated using UCSF Chimera (Pettersen et al., 2004), Chimera X (Pettersen et al., 2021), and PyMOL (Schrödinger).

Cryo-EM 3D variability analysis

Request a detailed protocol

3D variability analysis (3DVAR) was performed to access and visualize the dynamics within the cryo-EM datasets of the M4 mAChR complexes, as previously described using cryoSPARC (Punjani and Fleet, 2021). The polished particle stacks were imported into cryoSPARC, followed by 2D classification and 3D refinement using the respective low-pass-filtered RELION consensus maps as an initial model. 3DVA was analyzed in three components with 20 volume frames of data per component of motion. Output files were visualized using UCSF Chimera (Pettersen et al., 2004).

Gaussian accelerated molecular dynamics (GaMD)

Request a detailed protocol

GaMD enhances the conformational sampling of biomolecules by adding a harmonic boost potential to reduce the system energy barriers (Miao et al., 2015). When the system potential Vr is lower than a reference energy E, the modified potential Vr of the system is calculated as

Vr=Vr+Vr
(1) ΔV(r)={12k(EV(r))2,V(r)<E0,V(r)E,

where k is the harmonic force constant. The two adjustable parameters E and k are automatically determined on three enhanced sampling principles. First, for any two arbitrary potential values v1r and v2r found on the original energy surface, if V1(r)<V2(r) , V should be a monotonic function that does not change the relative order of the biased potential values; that is, V1(r)<V2(r) . Second, if V1(r)<V2(r) , the potential difference observed on the smoothened energy surface should be smaller than that of the original; i.e., V2(r)V1(r)<V2(r)V1(r) . By combining the first two criteria and plugging in the formula of Vr and V, we obtain

(2) VmaxEVmin+1k

where Vmin and Vmax are the system minimum and maximum potential energies. To ensure that Equation 2 is valid, k has to satisfy k1/(VmaxVmin) . Let us define k=k01/Vmax-Vmin , then 0k01. Third, the standard deviation (SD) of V needs to be small enough (i.e. narrow distribution) to ensure accurate reweighting using cumulant expansion to the second order: σΔV=k(EVavg)σVσ0 , where Vavg and σV are the average and SD of V with σ0 as a user-specified upper limit (e.g. 10kBT) for accurate reweighting. When E is set to the lower bound E=Vmax according to Equation 2, k0 can be calculated as

(3) k0=min1.0,k0`=min1.0,σ0σVVmax-VminVmax-Vavg

Alternatively, when the threshold energy E is set to its upper bound E=Vmin+1/k, k0 is set to

(4) k0=k0``1-σ0σVVmax-VminVavg-Vmin

If k0`` is calculated between 0 and 1. Otherwise, k0 is calculated using Equation 3.

Energetic reweighting of GaMD simulations

Request a detailed protocol

For energetic reweighting of GaMD simulations to calculate potential of mean force (PMF), the probability distribution along a reaction coordinate is written as pA . Given the boost potential Vr of each frame, pA can be reweighted to recover the canonical ensemble distribution pA , as

(5) pAj=pAjeβVrji=1MpAieβVri,j=1,,M

where M is the number of bins, β=kBT, and eβVrj is the ensemble-averaged Boltzmann factor of Vr for simulation frames found in the jth bin. The ensemble-averaged reweighting factor can be approximated using cumulant expansion:

(6) eβVr=expk=1βkk!Ck

where the first two cumulants are given by

(7) C1=V,C2=V2-V2=σv2.

The boost potential obtained from GaMD simulations usually follows near-Gaussian distribution (Miao and McCammon, 2017). Cumulant expansion to the second order thus provides a good approximation for computing the reweighting factor (Miao et al., 2015; Miao et al., 2014). The reweighted free energy FA=-kBTlnpA is calculated as

(8) FA=FA-k=12βkk!Ck+Fc

where FA=-kBTlnpA is the modified free energy obtained from GaMD simulation and Fc is a constant.

System setup

Request a detailed protocol

The M4R-ACh-Gi1, M4R-Ipx-Gi1, M4R-Ipx-Gi1-VU154, and M4R-Ipx-Gi1-LY298 cryo-EM structures were used for setting up simulation systems. The scFv16 in the cryo-EM structures was omitted in all simulations. The initial structures of single mutant D432E and T433R mutant of M4R-Ipx-Gi1-VU154 were obtained by mutating the corresponding residues in the M4R-Ipx-Gi1-VU154 cryo-EM structure. The initial structures of M4R-ACh-Gi1-VU154 and M4R-ACh-Gi1-LY298 were obtained from M4R-Ipx-Gi1-VU154 and M4R-Ipx-Gi1-LY298 cryo-EM structures by replacing Ipx with ACh through alignment of receptors to the M4R-ACh-Gi1 cryo-EM structure. The initial structures of M4R-Gi1-VU154 and M4R-Gi1-LY298 were obtained by removing the corresponding Ipx agonist from the M4R-Ipx-Gi1-VU154 and M4R-Ipx-Gi1-LY298 cryo-EM structures. The initial structures of M4R-VU154 and M4R-LY298 were obtained by removing the corresponding Ipx agonist and Gi1 protein from the M4R-Ipx-Gi1-VU154 and M4R-Ipx-Gi1-LY298 cryo-EM structures. According to previous findings, intracellular loop (ICL) 3 is highly flexible and removal of ICL3 does not appear to affect GPCR function (Dror et al., 2015; Dror et al., 2011). The ICL3 was thus omitted as in the current GaMD simulations. Similar to a previous study, helical domains of the Gi1 protein missing in the cryo-EM structures were not included in the simulation models. This was based on earlier simulation of the β2AR-Gs complex, which showed that the helical domain fluctuated substantially (Dror et al., 2015). All chain termini were capped with neutral groups (acetyl and methylamide). All the disulfide bonds in the complexes (i.e. Cys1083.25-Cys18545x50 and Cys426ECL3-Cys429ECL3 in the M4R) that were resolved in the cryo-EM structures were maintained in the simulations. Using the psfgen plugin in VMD (Humphrey et al., 1996), missing atoms in protein residues were added and all protein residues were set to the standard CHARMM protonation states at neutral pH. For each of the complex systems, the receptor was inserted into a palmitoyl-oleoyl-phosphatidyl-choline (POPC) bilayer with all overlapping lipid molecules removed using the membrane plugin in VMD. The system charges were then neutralized at 0.15 M NaCl using the solvate plugin in VMD (Humphrey et al., 1996). The simulation systems were summarized in Table 3.

Simulation protocol

Request a detailed protocol

The CHARMM36M parameter set (Huang et al., 2017; Klauda et al., 2010; Vanommeslaeghe and MacKerell, 2015) was used for the M4 mAChRs, Gi1 proteins, and POPC lipids. Force field parameters of agonists ACh and Ipx, PAMs LY298 and VU154 were obtained from the CHARMM ParamChem web server (Vanommeslaeghe et al., 2012b; Vanommeslaeghe and MacKerell, 2012a). Force field parameters with high penalty were optimized with FFParm (Kumar et al., 2020). GaMD simulations of these systems followed a similar protocol used in previous studies of GPCRs (Draper-Joyce et al., 2021; Miao and McCammon, 2018; Miao and McCammon, 2016). For each of the complex systems, initial energy minimization, thermalization, and 20 ns cMD equilibration were performed using NAMD2.12 (Phillips et al., 2005). A cutoff distance of 12 Å was used for the van der Waals and short-range electrostatic interactions and the long-range electrostatic interactions were computed with the particle-mesh Ewald summation method (Darden et al., 1993). A 2-fs integration time step was used for all MD simulations, and a multiple-time-stepping algorithm was used with bonded and short-range non-bonded interactions computed every time step and long-range electrostatic interactions every two-time steps. The SHAKE algorithm (Ryckaert et al., 1977) was applied to all hydrogen-containing bonds. The NAMD simulation started with equilibration of the lipid tails. With all other atoms fixed, the lipid tails were energy minimized for 1000 steps using the conjugate gradient algorithm and melted with a constant number, volume, and temperature (NVT) run for 0.5 ns at 310 K. The 12 systems were further equilibrated using a constant number, pressure, and temperature (NPT) run at 1 atm and 310 K for 10 ns with 5 kcal/(mol. Å2) harmonic position restraints applied to the protein and ligand atoms. Final equilibration of each system was performed using a NPT run at 1 atm pressure and 310 K for 0.5 ns with all atoms unrestrained. After energy minimization and system equilibration, conventional MD simulations were performed on each system for 20 ns at 1 atm pressure and 310 K with a constant ratio constraint applied on the lipid bilayer in the X-Y plane.

With the NAMD output structure, along with the system topology and CHARMM36M force field files, the ParmEd tool in the AMBER package was used to convert the simulation files into the AMBER format. The GaMD module implemented in the GPU version of AMBER20 (Case et al. 2020) was then applied to perform the GaMD simulation. GaMD simulations of systems with Gi1 protein (M4R-ACh-Gi1, M4R-Ipx-Gi1, M4R-Ipx-Gi1-VU154, M4R-Ipx-Gi1-LY298, M4R-ACh-Gi1-VU154, M4R-ACh-Gi1-LY298, single mutant D432E and T433R mutants of M4R-Ipx-Gi1-VU154) included an 8-ns short cMD simulation used to collect the potential statistics for calculating GaMD acceleration parameters, a 48-ns equilibration after adding the boost potential, and finally three independent 500-ns GaMD production simulations with randomized initial atomic velocities. The average and SD of the system potential energies were calculated every 800,000 steps (1.6 ns). GaMD simulations of M4R-VU154 and M4R-LY298 included a 2.4-ns short cMD simulation used to collect the potential statistics for calculating GaMD acceleration parameters, a 48-ns equilibration after adding the boost potential, and finally three independent 1000-ns GaMD production simulations with randomized initial atomic velocities. The average and SD of the system potential energies were calculated every 240,000 steps (0.48 ns). All GaMD simulations were run at the ‘dual-boost’ level by setting the reference energy to the lower bound. One boost potential is applied to the dihedral energetic term and the other to the total potential energetic term. The upper limit of the boost potential SD, σ0 was set to 6.0 kcal/mol for both the dihedral and the total potential energetic terms. Similar temperature and pressure parameters were used as in the NAMD simulations.

Simulation analysis

Request a detailed protocol

CPPTRAJ (Roe and Cheatham, 2013) and VMD (Humphrey et al., 1996) were used to analyze the GaMD simulations. The RMSDs of the agonist ACh and Ipx, PAM VU154 and LY298 relative to the simulation starting structures, the interactions between receptor and agonists/PAMs, distances between the receptor TM3 and TM6 intracellular ends were selected as reaction coordinates. Particularly, distances were calculated between the Cα atoms of residues Arg3.50 and Thr6.30, N atom of residue N1173.37 and carbon atom (C5) in the acetyl group of ACh or oxygen atom (O09) in the ether bond of Ipx, NE1 atom of residue W1644.67 and carbon atom (C5) in the acetyl group of ACh or oxygen atom (O09) in the ether bond of Ipx, indole ring of residue W4136.48 and acetyl group of ACh or heterocyclic isoazoline group of Ipx, OH atom of residue Y892.61 and oxygen atom in the amide group of VU154/LY298, benzene ring of residue F18645.51 and aromatic core of the PAMs VU154/LY298, OH atom of residue Y4397.39 and nitrogen atoms in the amine group of the PAMs VU154/LY298, CD atom of residue Q18445.49 and nitrogen atom in the amide group of VU154/LY298, CG atom of residue N4236.58 and chlorine atom in PAM LY298, OH atom of residue Y922.64 and nitrogen atom in the amide group of VU154, OG1 atom of residue T4337.33 and sulfur atom in the trifluoromethylsulfonyl group of VU154. In addition, the χ2 angle of residue W4136.48 and W4357.35 were calculated. Time courses of these reaction coordinates obtained from the GaMD simulation were plotted in the respective figures. The PyReweighting (Miao et al., 2014) toolkit was applied to reweight GaMD simulations to recover the original free energy or PMF profiles of the simulation systems. PMF profiles were computed using the combined trajectories from all the three independent 500 ns GaMD simulations for each system. A bin size of 1.0 Å was used for RMSD. The cutoff was set to 500 frames for 2D PMF calculations. The 2D PMF profiles were obtained for wildtype M4R-Ipx-Gi1-LY298, M4R-Ipx-Gi1-VU154, and the D432E and T433R single mutants of the M4R-Ipx-Gi1-VU154 system regarding the RMSDs of the agonist Ipx and the RMSDs of the PAMs relative to the cryo-EM conformation.

Data analysis

Request a detailed protocol

All pharmacological data was fit using GraphPad Prism 9.2.0. Saturation binding experiments to determine Bmax and pKd values were determined as previously described (Leach et al., 2011; Nawaratne et al., 2010; Thal et al., 2016). Detailed equations and analysis details can be found in Appendix 1. Interaction inhibition binding curves between [3H]-NMS, agonists (ACh or Ipx), and PAMs (LY298 or VU154) were analyzed using the allosteric ternary complex model to calculate binding affinity values for each ligand (pKA – for ACh/Ipx and pKB for LY298/VU154) and the degree of binding modulation between agonist and PAM (log α) (Christopoulos and Kenakin, 2002). The pKB values for LY298 and VU154 were determined from global fits of the ACh and Ipx curves to generate one pKB value per ligand (Ehlert, 1988; Leach et al., 2011; Nawaratne et al., 2010; Thal et al., 2016). All pERK1/2 and TruPath assays were analyzed using the operational model allosterism and agonism to determine values of orthosteric (τA) or allosteric efficacy (τB) and the functional modulation (log αβ) between the agonists and PAMs (Leach et al., 2011; Nawaratne et al., 2010). Binding affinities of the agonists and the PAMs were fixed to values determined from equilibrium binding assays. The τB values for LY298 and VU154 were determined from global fits of the ACh and Ipx curves (when possible) to generate one value per ligand. For comparison between WT human M4 mAChR and other M4 mAChR constructs, the log τ values were corrected (denoted log τC) by normalizing to Bmax values from saturation binding experiments (Leach et al., 2011; Nawaratne et al., 2010; Thal et al., 2016). All affinity, potency, and cooperativity values were estimated as logarithms, and statistical analysis between WT and mutant M4 mAChR was determined by one-way ANOVA using a Dunnett’s post-hoc test with a value of p<0.05 considered as significant in this study.

Appendix 1

Data analysis

Pharmacological parameters related to ligand binding

Interaction radioligand binding data were analyzed according to the following adapted form of an allosteric ternary complex model that accounts for the interaction of two orthosteric ligands and one allosteric ligand on a receptor (Christopoulos and Kenakin, 2002; Leach et al., 2007; Leach et al., 2011):

Y=Bmax[A][A]+[(KAKB)/(αA[B]+KB)]{1+([I]/KI)+([B]/KB)+[(αI[I][B]/KIKB])}
  • [A], [B], and [I] represent the concentrations of the radioligand ([3H]-NMS), allosteric ligand, and orthosteric inhibitor, respectively.

  • KA, KB, and KI represent their respective equilibrium dissociation constants. The value KA was fixed to the value determined from saturation binding experiments.

  • Bmax is the total number of receptors.

  • αA and αI represent the affinity cooperativity values between the allosteric ligand and the radioligand or orthosteric inhibitor, respectively. Values greater than 1 indicate positive cooperativity, values <1 (but >0) indicate negative cooperativity, and values of unity indicate neutral cooperativity.

  • All potency, affinity, and cooperativity parameters were estimated as logarithms.

Agonist binding affinity with a PAM bound:

Y=pKI+log αI
ΔY=((ΔpKI)2+(Δlog αI)2)1/2

PAM binding affinity for the agonist-bound state:

Y=pKB+log αI
ΔY=((ΔpKB)2+(Δlog αI)2)1/2

Global analysis for sharing the pKB value between interaction experiments using two different agonists:

Define:

Ag1=agonist 1 (e.g. ACh)

Ag2=agonist 2 (e.g. Ipx)

alpha = αA

betta = αI

Data for agonist 1 in columns <A:F>

Data for agonist 2 in columns <G:L>

Prism equation:

KA=10LogKA
KB=10LogKB
I=10X
A=10LogHotA
KIAg1=10LogKIAg1 (Ki for agonist 1)
KIAg2=10LogKIAg2 (Ki for agonist 2)
alphaAg1=10LogalphaAg1
bettaAg1=10LogbettaAg1
alphaAg2=10LogalphaAg2
bettaAg2=10LogbettaAg2
Part1Ag1=(KAKB)/(alphaAg1B+KB)
Part2Ag1=1+I/KIAg1+B/KB+(bettaAg1IB)/(KIAg1KB)
KAppAg1=Part1Ag1Part2Ag1
Part1Ag2=(KAKB)/(alphaAg2B+KB)
Part2Ag2=1+I/KIAg2+B/KB+(bettaAg2IB)/(KIAg2KB)
KAppAg2=Part1Ag2Part2Ag2
<A:F>Y=(BmaxAg1A)/(A+KAppAg1)
<G:L>Y=(BmaxAg1A)/(A+KAppAg2)

Pharmacological parameters related to function

To determine efficacy values (τA) of ACh and Ipx in TruPath and pERK1/2 assays, data were directly fit to the operational model of agonism (Black and Leff, 1983):

Y=Basal+Em-Basal1+KA+AτA×A

where Em is the maximal response of the system, Basal is the basal level of response in the absence of agonist, [A] is the concentration of agonist, KA denotes the equilibrium dissociation constant of the agonist (A), and τA is the index of the coupling efficiency/efficacy of the agonist. Values of KA were constrained to the corresponding KI values from interaction radioligand binding experiments.

The determination of τB and αβ values for LY298 and VU154 in TruPath and pERK1/2 assays data was directly fit to the following operational model of allosterism and agonism (Leach et al., 2007):

Y=Basal+Em-BasalAKB+αβAB+τBABEC50nAKB+αβAB+τBABEC50n+EC50nKB+Bn

where Em is the maximal response of the system, Basal is the basal level of response in the absence of agonist, [A] and [B] are the concentrations of agonist and the PAM, respectively. KB denotes the equilibrium dissociation constant of the PAM. τB is the index of the coupling efficiency/efficacy of the PAM. α denotes the binding cooperativity between the agonist and PAM, whereas β denotes a scaling factor that quantifies the allosteric effect of the PAM on the orthosteric ligand efficacy. This model assumes that the agonists (A) are full agonists. n is the transducer slope that describes the stimulus–response coupling of the ligand–receptor to the signaling pathway and was constrained to 1. The allosteric binding affinity (KB) was constrained to the value determined from radioligand binding experiments. All potency, affinity, efficacy, and cooperativity parameters were estimated as logarithms.

Transduction coupling coefficients were calculated as follows Kenakin et al., 2012:

Y=logτAKA=logτA-logKA
ΔY=((ΔlogτA)2+(ΔKA)2)1/2

Efficacy modulation (β) was calculated as follows:

Y=logβ=logαβlogα
ΔY=((Δlog αβ)2+(Δlog α)2)1/2

Global analysis for sharing τB between interaction functional experiments using two different agonists:

Define:

Ag1=agonist 1 (e.g. ACh)

Ag2=agonist 2 (e.g. Ipx)

ab=αβ

Data for agonist 1 in columns <A:F>

Data for agonist 2 in columns <G:L>

Prism equation:

KB=10LogKB
tauB=10LogtauB
A=10X
B=ConcMod
EC50Ag1=10LogEC50Ag1
EC50Ag2=10LogEC50Ag2
abAg1=10LogAlphaBetaAg1
abAg2=10LogAlphaBetaAg2
Part1Ag1=(A(KB+abAg1B)+tauBBEC50Ag1)n
Part1Ag2=(A(KB+abAg2B)+tauBBEC50Ag2)n
Part2Ag1=(EC50Ag1n)((KB+B)n)
Part2Ag2=(EC50Ag2n)((KB+B)n)
SpanAg1=EmAg1BasalAg1
SpanAg2=EmAg2BasalAg2
<A:F>Y=BasalAg1+(SpanAg1Part1Ag1)/(Part1Ag1+Part2Ag1)
<G:L>Y=BasalAg2+(SpanAg2Part1Ag2)/(Part1Ag2+Part2Ag2)

Correcting efficacy τX values using receptor expression (Leach et al., 2011; Nawaratne et al., 2010):

Bmax–WT = maximal number of receptors for the WT receptor

Bmax–X = maximal number of receptors for the receptor construct X

τWT = efficacy at the WT receptor

τX = efficacy at receptor X

logτXcorr=logτX(log BmaxXlog BmaxWT)
ΔlogτXcorr=((ΔlogτX)2+(0.432×(ΔBmaxXBmaxX))2+(0.432×(ΔBmaxWTBmaxWT))2)1/2
Appendix 1—key resources table
Reagent type (species) or resourceDesignationSource or referenceIdentifiersAdditional information
AntibodyAnti-FLAG M1 (mouse polyclonal IgG2a)Gift from Prof. Brian
Kobilka (PMID::17962520)
Antibody was used to make anti-FLAG mAb resin that was used for the purification of FLAG-tagged M4 mAChR
Strain, strain background (Escherichia coli)DH5αNew England BiolabsC2987H
Strain, strain background (E. coli)DH10bacThermo Fisher Scientific10361012
Cell line (Spodoptera frugiperda)Sf9Expression Systems94-001S
Cell line (Trichoplusia ni)TniExpression Systems94-002S
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR WTPMID:26958838
Cell line (Chinese hamster ovary)CHO K1 mouse M4 mAChR WTPMID:21198541
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR D432EThis studypEF5-FRT-V5-DEST plasmid
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR T433RThis studypEF5-FRT-V5-DEST plasmid
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR
V91L, D432E, T433R
This studypEF5-FRT-V5-DEST plasmid
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR Y89APMID:26958838
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR Q184APMID:26958838
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR F186APMID:20406819
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR W435APMID:26958838
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR W439APMID:20406819
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR S85APMID:26958838
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR Y89APMID:26958838
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR Y92APMID:26958838
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR I93T, I94V, K95IPMID:20406819
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR I93TPMID:20406819
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR I94VPMID:20406819
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR K95IPMID:20406819
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR Y97APMID:26958838
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR W98APMID:26958838
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR G101APMID:26958838
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR D106APMID:21300722
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR W108APMID:21300722
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR L109APMID:21300722
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR D112EPMID:21300722
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR D112NPMID:21300722
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR S116APMID:21300722
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR N117APMID:21300722
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR V120APMID:21300722
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR D129EPMID:21300722
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR D129NPMID:21300722
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR W164APMID:26958838
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR F170APMID:26958838
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR W171APMID:26958838
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR Q172APMID:26958838
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR F173APMID:26958838
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR Q184APMID:26958838
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR F186APMID:20406819
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR I187APMID:26958838
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR Q188APMID:26958838
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR F189APMID:26958838
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR L190APMID:26958838
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR W413APMID:26958838
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR Y416APMID:26958838
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR N423APMID:26958838
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR Q427APMID:26958838
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR S428PPMID:20406819
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR D432NPMID:20406819
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR W435APMID:26958838
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR Y439APMID:20406819
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR Y440APMID:26958838
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR C442APMID:20406819
Cell line (Chinese hamster ovary)Flp-In CHO human M4 mAChR Y443APMID:20406819
Cell line (Chinese hamster ovary)Flp-In CHO cell lineThermo Fisher ScientificR75807
Cell line (Homo sapiens)293A cell lineThermo Fisher ScientificR70507
Recombinant DNA reagentHuman FLAG-M4∆i3-HisThis studypVL1392 vector
Recombinant DNA reagentHuman Gαi1 dominant negative mutantPMID:32193322pFastBac vector
Recombinant DNA reagentHuman Gβ1γ2Gift from Prof. Brian
Kobilka (PMID:24256733)
pVL1392 vector
Recombinant DNA reagentscFv16PMID:30213947pFastBac vector
Recombinant DNA reagentpcDNA5/FRT/TO-GAlphai1-RLuc8Gift from Prof. Bryan
Roth (PMID:32367019)
TRUPATH assay
Recombinant DNA reagentpcDNA3.1-Beta3Gift from Prof. Bryan
Roth (PMID:32367019)
TRUPATH assay
Recombinant DNA reagentpcDNA3.1-GGamma9-GFP2Gift from Prof. Bryan
Roth (PMID:32367019)
TRUPATH assay
Chemical compound, drugAcetylcholineSigma-Aldrich
Chemical compound, drugIperoxoSigma-Aldrich
Chemical compound, drugLY2033298Sigma-Aldrich
Chemical compound, drugVU0467154Gift from Prof. Craig Lindsley (PMID:25137629)PMID:25137629
Chemical compound, drugProlume PurpleNanolight Technology
Chemical compound, drug[3H]-NMSPerkinElmer
Chemical compound, drugPolyethylenimine (PEI)Sigma-Aldrich
Chemical compound, drugAtropineSigma-Aldrich
Commercial assay or kitAlphaScreen SureFire pERK 1/2
(Thr202/Tyr204) Assay Kits
PerkinElmer
Software, algorithmPrism 8.0GraphPad
Software, algorithmPyMOLSchrödinger
Software, algorithmGaMDPMID:26300708
Software, algorithmAMBER20https://ambermd.org
Software, algorithmCPPTRAJPMID:26300708
Software, algorithmPyReweightingPMID:25061441
Software, algorithmPhenix suitePMID:20124702
Software, algorithmCootPMID:31730249
Software, algorithmChimeraPMID:15264254
Software, algorithmChimera XPMID:32881101
Software, algorithmcryoSPARCPMID:28165473
Software, algorithmRelion 3.1PMID:30412051
Software, algorithmMotioncor2PMID:28250466
Software, algorithmGCTFPMID:26592709
Software, algorithmcrYOLOPMID:31240256
Software, algorithmISOLDEPMID:29872003
Software, algorithmMolProbityPMID:29067766
Software, algorithmDAQ scorePMID:35953671

Data availability

All data generated or analysed during this study are included in the manuscript. Structural data has been deposited in the Protein Data Bank (PDB) and Electron Microscopy Data Bank (EMDB) under the following codes:(1) M4R-Gi1-Ipx PDB: 7TRK and EMD-26099 (2) M4R-Gi1-Ipx-LY298 PDB: 7TRP and EMD-26100 (3) M4R-Gi1-Ipx-VU154 PDB: 7TRQ and EMD-26101 (4) M4R-Gi1-ACh PDB: 7TRS and EMD-26102.

The following data sets were generated
    1. Vuckovic Z
    2. Mobbs JI
    3. Belousoff MJ
    4. Glukhova A
    5. Sexton PM
    6. Danev R
    7. Thal DM
    (2023) RCSB Protein Data Bank
    ID 7TRK. Human M4 muscarinic acetylcholine receptor complex with Gi1 and the agonist iperoxo.
    1. Vuckovic Z
    2. Mobbs JI
    3. Belousoff MJ
    4. Glukhova A
    5. Sexton PM
    6. Danev R
    7. Thal DM
    (2023) RCSB Protein Data Bank
    ID 7TRP. Human M4 muscarinic acetylcholine receptor complex with Gi1 and the agonist iperoxo and positive allosteric modulator LY2033298.
    1. Vuckovic Z
    2. Mobbs JI
    3. Belousoff MJ
    4. Glukhova A
    5. Sexton PM
    6. Danev R
    7. Thal DM
    (2023) RCSB Protein Data Bank
    ID 7TRQ. Human M4 muscarinic acetylcholine receptor complex with Gi1 and the agonist iperoxo and positive allosteric modulator VU0467154.
    1. Vuckovic Z
    2. Mobbs JI
    3. Belousoff MJ
    4. Glukhova A
    5. Sexton PM
    6. Danev R
    7. Thal DM
    (2023) RCSB Protein Data Bank
    ID 7TRS. Human M4 muscarinic acetylcholine receptor complex with Gi1 and the endogenous agonist acetylcholine.
    1. Vuckovic Z
    2. Mobbs JI
    3. Belousoff MJ
    4. Glukhova A
    5. Sexton PM
    6. Danev R
    7. Thal DM
    (2023) Electron Microscopy Data Bank
    ID EMD-26099. Human M4 muscarinic acetylcholine receptor complex with Gi1 and the agonist iperoxo.
    1. Vuckovic Z
    2. Mobbs JI
    3. Belousoff MJ
    4. Glukhova A
    5. Sexton PM
    6. Danev R
    7. Thal DM
    (2023) Electron Microscopy Data Bank
    ID EMD-26100. Human M4 muscarinic acetylcholine receptor complex with Gi1 and the agonist iperoxo and positive allosteric modulator LY2033298.
    1. Vuckovic Z
    2. Mobbs JI
    3. Belousoff MJ
    4. Glukhova A
    5. Sexton PM
    6. Danev R
    7. Thanl DM
    (2023) Electron Microscopy Data Bank
    ID EMD-26101. Human M4 muscarinic acetylcholine receptor complex with Gi1 and the agonist iperoxo and positive allosteric modulator VU0467154.
    1. Vuckovic Z
    2. Mobbs JI
    3. Belousoff MJ
    4. Glukhova A
    5. Sexton PM
    6. Danev R
    7. Thal DM
    (2023) Electron Microscopy Data Bank
    ID EMD-26102. Human M4 muscarinic acetylcholine receptor complex with Gi1 and the endogenous agonist acetylcholine.

References

    1. Black JW
    2. Leff P
    (1983) Operational models of pharmacological agonism
    Proceedings of the Royal Society of London. Series B, Biological Sciences 220:141–162.
    https://doi.org/10.1098/rspb.1983.0093
    1. Ehlert FJ
    (1988)
    Estimation of the affinities of allosteric ligands using radioligand binding and pharmacological null methods
    Molecular Pharmacology 33:187–194.

Article and author information

Author details

  1. Ziva Vuckovic

    Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
    Contribution
    Investigation, Methodology, Writing – review and editing
    Contributed equally with
    Jinan Wang, Vi Pham and Jesse I Mobbs
    Competing interests
    No competing interests declared
  2. Jinan Wang

    Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, United States
    Contribution
    Software, Formal analysis, Validation, Investigation, Visualization, Methodology, Writing – review and editing
    Contributed equally with
    Ziva Vuckovic, Vi Pham and Jesse I Mobbs
    Competing interests
    No competing interests declared
  3. Vi Pham

    Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
    Contribution
    Formal analysis, Investigation, Writing – review and editing
    Contributed equally with
    Ziva Vuckovic, Jinan Wang and Jesse I Mobbs
    Competing interests
    No competing interests declared
  4. Jesse I Mobbs

    1. Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
    2. ARC Centre for Cryo-electron Microscopy of Membrane Proteins, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
    Contribution
    Validation, Investigation, Visualization, Writing – review and editing
    Contributed equally with
    Ziva Vuckovic, Jinan Wang and Vi Pham
    Competing interests
    No competing interests declared
  5. Matthew J Belousoff

    1. Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
    2. ARC Centre for Cryo-electron Microscopy of Membrane Proteins, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
    Contribution
    Resources, Formal analysis
    Competing interests
    No competing interests declared
  6. Apurba Bhattarai

    Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, United States
    Contribution
    Investigation, Methodology, Writing – review and editing
    Competing interests
    No competing interests declared
  7. Wessel AC Burger

    1. Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
    2. ARC Centre for Cryo-electron Microscopy of Membrane Proteins, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
    Contribution
    Formal analysis, Investigation, Writing – review and editing
    Competing interests
    No competing interests declared
  8. Geoff Thompson

    Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
    Contribution
    Investigation
    Competing interests
    No competing interests declared
  9. Mahmuda Yeasmin

    Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
    Contribution
    Investigation
    Competing interests
    No competing interests declared
  10. Vindhya Nawaratne

    Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
    Contribution
    Investigation
    Competing interests
    No competing interests declared
  11. Katie Leach

    1. Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
    2. ARC Centre for Cryo-electron Microscopy of Membrane Proteins, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
    Contribution
    Writing – review and editing
    Competing interests
    No competing interests declared
  12. Emma T van der Westhuizen

    Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
    Contribution
    Investigation, Writing – review and editing
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-9165-8526
  13. Elham Khajehali

    Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
    Contribution
    Investigation, Writing – review and editing
    Competing interests
    No competing interests declared
  14. Yi-Lynn Liang

    Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
    Contribution
    Investigation
    Competing interests
    No competing interests declared
  15. Alisa Glukhova

    1. Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
    2. ARC Centre for Cryo-electron Microscopy of Membrane Proteins, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
    Contribution
    Investigation, Writing – review and editing
    Competing interests
    No competing interests declared
  16. Denise Wootten

    1. Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
    2. ARC Centre for Cryo-electron Microscopy of Membrane Proteins, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
    Contribution
    Resources, Supervision, Project administration, Writing – review and editing
    Competing interests
    P.M.S, D.W., and A.C. are shareholders of Septerna Inc
  17. Craig W Lindsley

    Department of Pharmacology, Warren Center for Neuroscience Drug Discovery and Department of Chemistry, Warren Center for Neuroscience Drug Discovery, Vanderbilt University, Nashville, United States
    Contribution
    Resources, Writing – review and editing
    Competing interests
    No competing interests declared
  18. Andrew Tobin

    The Centre for Translational Pharmacology, Advanced Research Centre (ARC), College of Medical, Veterinary and Life Sciences, University of Glasgow, Glasgow, United Kingdom
    Contribution
    Conceptualization, Resources, Funding acquisition, Project administration, Writing – review and editing
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-1807-3123
  19. Patrick Sexton

    1. Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
    2. ARC Centre for Cryo-electron Microscopy of Membrane Proteins, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
    Contribution
    Conceptualization, Resources, Supervision, Funding acquisition, Project administration, Writing – review and editing
    Competing interests
    P.M.S, D.W., and A.C. are shareholders of Septerna Inc
  20. Radostin Danev

    Graduate School of Medicine, University of Tokyo, Tokyo, Japan
    Contribution
    Resources, Funding acquisition, Investigation, Writing – review and editing
    Competing interests
    No competing interests declared
  21. Celine Valant

    Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
    Contribution
    Conceptualization, Formal analysis, Supervision, Validation, Visualization, Writing – original draft, Project administration, Writing – review and editing
    For correspondence
    celine.valant@monash.edu
    Competing interests
    No competing interests declared
  22. Yinglong Miao

    Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, United States
    Contribution
    Conceptualization, Resources, Software, Formal analysis, Supervision, Funding acquisition, Validation, Visualization, Methodology, Writing – original draft, Project administration, Writing – review and editing
    For correspondence
    miao@ku.edu
    Competing interests
    No competing interests declared
  23. Arthur Christopoulos

    1. Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
    2. ARC Centre for Cryo-electron Microscopy of Membrane Proteins, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
    3. Neuromedicines Discovery Centre, Monash University, Parkville, Australia
    Contribution
    Conceptualization, Resources, Supervision, Funding acquisition, Writing – original draft, Project administration, Writing – review and editing
    For correspondence
    Arthur.Christopoulos@monash.edu
    Competing interests
    P.M.S, D.W., and A.C. are shareholders of Septerna Inc
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-4442-3294
  24. David M Thal

    1. Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
    2. ARC Centre for Cryo-electron Microscopy of Membrane Proteins, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
    Contribution
    Conceptualization, Resources, Software, Formal analysis, Supervision, Funding acquisition, Validation, Investigation, Visualization, Methodology, Writing – original draft, Project administration, Writing – review and editing
    For correspondence
    david.thal@monash.edu
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-0325-2524

Funding

Wellcome Trust (201529/Z/16/Z)

  • Andrew Tobin
  • Patrick Sexton
  • Arthur Christopoulos

National Health and Medical Research Council (1055134)

  • Arthur Christopoulos

National Health and Medical Research Council (1150083)

  • Arthur Christopoulos

National Health and Medical Research Council (1138448)

  • David M Thal

Australian Research Council (DE170100152)

  • David M Thal

Australian Research Council (DP190102950)

  • Celine Valant

Australian Research Council (IC200100052)

  • Denise Wootten

National Institutes of Health (GM132572)

  • Yinglong Miao

National Health and Medical Research Council (1154434)

  • Patrick Sexton

National Health and Medical Research Council (1155302)

  • Denise Wootten

National Health and Medical Research Council (1196951)

  • David M Thal

Australian Research Council (160100075)

  • Katie Leach

Takeda Science Foundation (2019 Medical Research Grant)

  • Radostin Danev

Japan Science and Technology Agency (18069571)

  • Radostin Danev

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication. For the purpose of Open Access, the authors have applied a CC BY public copyright license to any Author Accepted Manuscript version arising from this submission.

Acknowledgements

This work was supported by a Wellcome Trust Collaborative Award (201529/Z/16/Z; PMS, ABT, AC), the National Health and Medical Research Council of Australia (1055134, 1150083, and 1138448), the Australian Research Council (DE170100152, DP190102950, and IC200100052), and the National Institutes of Health (GM132572). PMS is a Senior Principal Research Fellow (1154434), DW a Senior Research Fellow (1155302), DMT an Early Career Research Fellow (1196951), and KL a Future Fellow (160100075). RD was supported by Takeda Science Foundation 2019 Medical Research Grant and Japan Science and Technology Agency PRESTO (18069571). This work was partially supported by the Monash University Ramaciotti Centre for cryo-electron microscopy and the Monash University MASSIVE high-performance computing facility and supercomputing resources with the XSEDE allocation award TG-MCB180049, BIO220137 from the Advanced Cyberinfrastructure Coordination Ecosystem: Services & Support (ACCESS) program, and NERSC project M2874. We thank John Tesmer for discussion and the suggestion of calculating DAQ scores.

Copyright

© 2023, Vuckovic, Wang, Pham 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

  • 3,105
    views
  • 478
    downloads
  • 13
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

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)

  1. Ziva Vuckovic
  2. Jinan Wang
  3. Vi Pham
  4. Jesse I Mobbs
  5. Matthew J Belousoff
  6. Apurba Bhattarai
  7. Wessel AC Burger
  8. Geoff Thompson
  9. Mahmuda Yeasmin
  10. Vindhya Nawaratne
  11. Katie Leach
  12. Emma T van der Westhuizen
  13. Elham Khajehali
  14. Yi-Lynn Liang
  15. Alisa Glukhova
  16. Denise Wootten
  17. Craig W Lindsley
  18. Andrew Tobin
  19. Patrick Sexton
  20. Radostin Danev
  21. Celine Valant
  22. Yinglong Miao
  23. Arthur Christopoulos
  24. David M Thal
(2023)
Pharmacological hallmarks of allostery at the M4 muscarinic receptor elucidated through structure and dynamics
eLife 12:e83477.
https://doi.org/10.7554/eLife.83477

Share this article

https://doi.org/10.7554/eLife.83477

Further reading

    1. Biochemistry and Chemical Biology
    2. Structural Biology and Molecular Biophysics
    Lingzhi Gao, Dian Chen, Yu Liu
    Research Article

    Riboswitches represent a class of non-coding RNA that possess the unique ability to specifically bind ligands and, in response, regulate gene expression. A recent report unveiled a type of riboswitch, known as the guanidine-IV riboswitch, which responds to guanidine levels to regulate downstream genetic transcription. However, the precise molecular mechanism through which the riboswitch senses its target ligand and undergoes conformational changes remain elusive. This gap in understanding has impeded the potential applications of this riboswitch. To bridge this knowledge gap, our study investigated the conformational dynamics of the guanidine-IV riboswitch RNA upon ligand binding. We employed single-molecule fluorescence resonance energy transfer (smFRET) to dissect the behaviors of the aptamer, terminator, and full-length riboswitch. Our findings indicated that the aptamer portion exhibited higher sensitivity to guanidine compared to the terminator and full-length constructs. Additionally, we utilized Position-specific Labelling of RNA (PLOR) combined with smFRET to observe, at the single-nucleotide and single-molecule level, the structural transitions experienced by the guanidine-IV riboswitch during transcription. Notably, we discovered that the influence of guanidine on the riboswitch RNA’s conformations was significantly reduced after the transcription of 88 nucleotides. Furthermore, we proposed a folding model for the guanidine-IV riboswitch in the absence and presence of guanidine, thereby providing insights into its ligand-response mechanism.

    1. Biochemistry and Chemical Biology
    Marius Landau, Sherif Elsabbagh ... Joachim E Schultz
    Research Article

    The biosynthesis of cyclic 3′,5′-adenosine monophosphate (cAMP) by mammalian membrane-bound adenylyl cyclases (mACs) is predominantly regulated by G-protein-coupled receptors (GPCRs). Up to now the two hexahelical transmembrane domains of mACs were considered to fix the enzyme to membranes. Here, we show that the transmembrane domains serve in addition as signal receptors and transmitters of lipid signals that control Gsα-stimulated mAC activities. We identify aliphatic fatty acids and anandamide as receptor ligands of mAC isoforms 1–7 and 9. The ligands enhance (mAC isoforms 2, 3, 7, and 9) or attenuate (isoforms 1, 4, 5, and 6) Gsα-stimulated mAC activities in vitro and in vivo. Substitution of the stimulatory membrane receptor of mAC3 by the inhibitory receptor of mAC5 results in a ligand inhibited mAC5–mAC3 chimera. Thus, we discovered a new class of membrane receptors in which two signaling modalities are at a crossing, direct tonic lipid and indirect phasic GPCR–Gsα signaling regulating the biosynthesis of cAMP.