Pharmacological hallmarks of allostery at the M4 muscarinic receptor elucidated through structure and dynamics
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.sa0Introduction
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.
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).
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).
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).
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.
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.
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.
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.
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).
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).
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.
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.
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).
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 protocolDH5α (New England Biolabs) and DH10bac (Thermo Fisher Scientific) Escherichia coli cells were grown in LB at 37°C.
Cell culture
Request a detailed protocolTni 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 protocolFlp-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 protocolUpon 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 protocolThe 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 protocolTni 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 protocolThe 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 protocolSamples (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 protocolSpecific 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 protocolAn 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αi1Gβ1Gγ2 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 protocol3D 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 protocolGaMD 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 is lower than a reference energy E, the modified potential of the system is calculated as
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 and found on the original energy surface, if , should be a monotonic function that does not change the relative order of the biased potential values; that is, . Second, if , the potential difference observed on the smoothened energy surface should be smaller than that of the original; i.e., . By combining the first two criteria and plugging in the formula of and , we obtain
where and are the system minimum and maximum potential energies. To ensure that Equation 2 is valid, k has to satisfy . Let us define , then . Third, the standard deviation (SD) of needs to be small enough (i.e. narrow distribution) to ensure accurate reweighting using cumulant expansion to the second order: , where and are the average and SD of with as a user-specified upper limit (e.g. ) for accurate reweighting. When E is set to the lower bound according to Equation 2, can be calculated as
Alternatively, when the threshold energy E is set to its upper bound , is set to
If is calculated between 0 and 1. Otherwise, is calculated using Equation 3.
Energetic reweighting of GaMD simulations
Request a detailed protocolFor energetic reweighting of GaMD simulations to calculate potential of mean force (PMF), the probability distribution along a reaction coordinate is written as . Given the boost potential of each frame, can be reweighted to recover the canonical ensemble distribution , as
where M is the number of bins, , and is the ensemble-averaged Boltzmann factor of for simulation frames found in the jth bin. The ensemble-averaged reweighting factor can be approximated using cumulant expansion:
where the first two cumulants are given by
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 is calculated as
where is the modified free energy obtained from GaMD simulation and is a constant.
System setup
Request a detailed protocolThe 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 protocolThe 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 protocolCPPTRAJ (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 protocolAll 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):
[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:
PAM binding affinity for the agonist-bound state:
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:
Pharmacological parameters related to function
To determine efficacy values () of ACh and Ipx in TruPath and pERK1/2 assays, data were directly fit to the operational model of agonism (Black and Leff, 1983):
where 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 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 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):
where 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. 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 () 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:
Efficacy modulation () was calculated as follows:
Global analysis for sharing between interaction functional experiments using two different agonists:
Define:
Ag1=agonist 1 (e.g. ACh)
Ag2=agonist 2 (e.g. Ipx)
Data for agonist 1 in columns <A:F>
Data for agonist 2 in columns <G:L>
Prism equation:
Correcting efficacy 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
= efficacy at the WT receptor
= efficacy at receptor X
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.
-
RCSB Protein Data BankID 7TRK. Human M4 muscarinic acetylcholine receptor complex with Gi1 and the agonist iperoxo.
-
RCSB Protein Data BankID 7TRP. Human M4 muscarinic acetylcholine receptor complex with Gi1 and the agonist iperoxo and positive allosteric modulator LY2033298.
-
RCSB Protein Data BankID 7TRQ. Human M4 muscarinic acetylcholine receptor complex with Gi1 and the agonist iperoxo and positive allosteric modulator VU0467154.
-
RCSB Protein Data BankID 7TRS. Human M4 muscarinic acetylcholine receptor complex with Gi1 and the endogenous agonist acetylcholine.
-
Electron Microscopy Data BankID EMD-26099. Human M4 muscarinic acetylcholine receptor complex with Gi1 and the agonist iperoxo.
-
Electron Microscopy Data BankID EMD-26100. Human M4 muscarinic acetylcholine receptor complex with Gi1 and the agonist iperoxo and positive allosteric modulator LY2033298.
-
Electron Microscopy Data BankID EMD-26101. Human M4 muscarinic acetylcholine receptor complex with Gi1 and the agonist iperoxo and positive allosteric modulator VU0467154.
-
Electron Microscopy Data BankID EMD-26102. Human M4 muscarinic acetylcholine receptor complex with Gi1 and the endogenous agonist acetylcholine.
References
-
BookIntegrated methods for the construction of three-dimensional models and computational probing of structure-function relations in G protein-coupled receptorsIn: Sealfon SC, editors. Methods in Neurosciences, Receptor Molecular Biology. Academic Press. pp. 366–428.https://doi.org/10.1016/S1043-9471(05)80049-7
-
Activation of the beta 2-adrenergic receptor involves disruption of an ionic lock between the cytoplasmic ends of transmembrane segments 3 and 6The Journal of Biological Chemistry 276:29171–29177.https://doi.org/10.1074/jbc.M103747200
-
Molecular mechanisms of action of M5 muscarinic acetylcholine receptor allosteric modulatorsMolecular Pharmacology 90:427–436.https://doi.org/10.1124/mol.116.104182
-
Announcing the worldwide protein data bankNature Structural Biology 10:980.https://doi.org/10.1038/nsb1203-980
-
Gaussian accelerated molecular dynamics for elucidation of drug pathwaysExpert Opinion on Drug Discovery 13:1055–1065.https://doi.org/10.1080/17460441.2018.1538207
-
Operational models of pharmacological agonismProceedings of the Royal Society of London. Series B, Biological Sciences 220:141–162.https://doi.org/10.1098/rspb.1983.0093
-
Conformational changes in tyrosine 11 of neurotensin are required to activate the neurotensin receptor 1ACS Pharmacology & Translational Science 3:690–705.https://doi.org/10.1021/acsptsci.0c00026
-
Toward an understanding of the structural basis of allostery in muscarinic acetylcholine receptorsThe Journal of General Physiology 150:1360–1372.https://doi.org/10.1085/jgp.201711979
-
Allostery in GPCRs: `` MWC'' revisitedTrends in Biochemical Sciences 36:663–672.https://doi.org/10.1016/j.tibs.2011.08.005
-
Structural and functional diversity among agonist-bound states of the GLP-1 receptorNature Chemical Biology 18:256–263.https://doi.org/10.1038/s41589-021-00945-w
-
G protein-coupled receptor allosterism and complexingPharmacological Reviews 54:323–374.https://doi.org/10.1124/pr.54.2.323
-
Advances in G protein-coupled receptor allostery: from function to structureMolecular Pharmacology 86:463–478.https://doi.org/10.1124/mol.114.094342
-
Subtype-Selective allosteric modulators of muscarinic receptors for the treatment of CNS disordersTrends in Pharmacological Sciences 30:148–155.https://doi.org/10.1016/j.tips.2008.12.002
-
ISOLDE: a physically realistic environment for model building into low-resolution electron-density mapsActa Crystallographica. Section D, Structural Biology 74:519–530.https://doi.org/10.1107/S2059798318002425
-
Routine sub-2.5 Å cryo-EM structure determination of GPCRsNature Communications 12:4333.https://doi.org/10.1038/s41467-021-24650-3
-
Particle mesh ewald: an N⋅log (N) method for ewald sums in large systemsThe Journal of Chemical Physics 98:10089–10092.https://doi.org/10.1063/1.464397
-
Estimation of the affinities of allosteric ligands using radioligand binding and pharmacological null methodsMolecular Pharmacology 33:187–194.
-
Novel allosteric modulators of G protein-coupled receptorsThe Journal of Biological Chemistry 290:19478–19488.https://doi.org/10.1074/jbc.R115.662759
-
Identification of orthosteric and allosteric site mutations in M2 muscarinic acetylcholine receptors that contribute to ligand-selective signaling biasThe Journal of Biological Chemistry 285:7459–7474.https://doi.org/10.1074/jbc.M109.094011
-
Gpcr activation: a mutagenic spotlight on crystal structuresTrends in Pharmacological Sciences 34:67–84.https://doi.org/10.1016/j.tips.2012.11.002
-
VMD: visual molecular dynamicsJournal of Molecular Graphics 14:33–38.https://doi.org/10.1016/0263-7855(96)00018-5
-
The ligand paradox between affinity and efficacy: can you be there and not make a difference?Trends in Pharmacological Sciences 23:275–280.https://doi.org/10.1016/S0165-6147(02)02036-9
-
New concepts in drug discovery: collateral efficacy and permissive antagonismNature Reviews. Drug Discovery 4:919–927.https://doi.org/10.1038/nrd1875
-
Biased signalling and allosteric machines: new vistas and challenges for drug discoveryBritish Journal of Pharmacology 165:1659–1669.https://doi.org/10.1111/j.1476-5381.2011.01749.x
-
A simple method for quantifying functional selectivity and agonist biasACS Chemical Neuroscience 3:193–203.https://doi.org/10.1021/cn200111m
-
Signalling bias in new drug discovery: detection, quantification and therapeutic impactNature Reviews. Drug Discovery 12:205–216.https://doi.org/10.1038/nrd3954
-
Update of the CHARMM all-atom additive force field for lipids: validation on six lipid typesThe Journal of Physical Chemistry. B 114:7830–7843.https://doi.org/10.1021/jp101759q
-
FFParam: standalone package for CHARMM additive and drude polarizable force field parametrization of small moleculesJournal of Computational Chemistry 41:958–970.https://doi.org/10.1002/jcc.26138
-
Allosteric GPCR modulators: taking advantage of permissive receptor pharmacologyTrends in Pharmacological Sciences 28:382–389.https://doi.org/10.1016/j.tips.2007.06.004
-
Dominant negative G proteins enhance formation and purification of agonist-GPCR-G protein complexes for structure determinationACS Pharmacology & Translational Science 1:12–20.https://doi.org/10.1021/acsptsci.8b00017
-
Macromolecular structure determination using x-rays, neutrons and electrons: recent developments in phenixActa Crystallographica Section D Structural Biology 75:861–877.https://doi.org/10.1107/S2059798319011471
-
Structural insights into binding specificity, efficacy and bias of a β2AR partial agonistNature Chemical Biology 14:1059–1066.https://doi.org/10.1038/s41589-018-0145-x
-
Structural determinants of 5-HT2B receptor activation and biased agonismNature Structural & Molecular Biology 25:787–796.https://doi.org/10.1038/s41594-018-0116-7
-
Improved reweighting of accelerated molecular dynamics simulations for free energy calculationJournal of Chemical Theory and Computation 10:2677–2689.https://doi.org/10.1021/ct500090q
-
Gaussian accelerated molecular dynamics: unconstrained enhanced sampling and free energy calculationJournal of Chemical Theory and Computation 11:3584–3595.https://doi.org/10.1021/acs.jctc.5b00436
-
Gaussian accelerated molecular dynamics: theory, implementation, and applicationsAnnual Reports in Computational Chemistry 13:231–278.https://doi.org/10.1016/bs.arcc.2017.06.005
-
Ligand binding and micro-switches in 7TM receptor structuresTrends in Pharmacological Sciences 30:249–259.https://doi.org/10.1016/j.tips.2009.02.006
-
TRUPATH, an open-source biosensor platform for interrogating the GPCR transduceromeNature Chemical Biology 16:841–849.https://doi.org/10.1038/s41589-020-0535-8
-
UCSF chimera? A visualization system for exploratory research and analysisJournal of Computational Chemistry 25:1605–1612.https://doi.org/10.1002/jcc.20084
-
Scalable molecular dynamics with NAMDJournal of Computational Chemistry 26:1781–1802.https://doi.org/10.1002/jcc.20289
-
3D variability analysis: resolving continuous flexibility and discrete heterogeneity from single particle cryo-EMJournal of Structural Biology 213:107702.https://doi.org/10.1016/j.jsb.2021.107702
-
PTRAJ and CPPTRAJ: software for processing and analysis of molecular dynamics trajectory dataJournal of Chemical Theory and Computation 9:3084–3095.https://doi.org/10.1021/ct400341p
-
Numerical integration of the Cartesian equations of motion of a system with constraints: molecular dynamics of n-alkanesJournal of Computational Physics 23:327–341.https://doi.org/10.1016/0021-9991(77)90098-5
-
Agonists with supraphysiological efficacy at the muscarinic M2 ACh receptorBritish Journal of Pharmacology 169:357–370.https://doi.org/10.1111/bph.12003
-
Structural basis of GPCR coupling to distinct signal transducers: implications for biased signalingTrends in Biochemical Sciences 47:570–581.https://doi.org/10.1016/j.tibs.2022.03.009
-
Beta2 adrenergic receptor activation. Modulation of the proline kink in transmembrane 6 by a rotamer toggle switchThe Journal of Biological Chemistry 277:40989–40996.https://doi.org/10.1074/jbc.M206801200
-
Biased allosteric modulators: new frontiers in GPCR drug discoveryTrends in Pharmacological Sciences 42:283–299.https://doi.org/10.1016/j.tips.2020.12.005
-
CASTp 3.0: computed atlas of surface topography of proteinsNucleic Acids Research 46:W363–W367.https://doi.org/10.1093/nar/gky473
-
Automation of the CHARMM General force field (cgenff) I: bond perception and atom typingJournal of Chemical Information and Modeling 52:3144–3154.https://doi.org/10.1021/ci300363c
-
Automation of the CHARMM General force field (cgenff) II: assignment of bonded parameters and partial atomic chargesJournal of Chemical Information and Modeling 52:3155–3168.https://doi.org/10.1021/ci3003649
-
CHARMM additive and polarizable force fields for biophysics and computer-aided drug designBiochimica et Biophysica Acta 1850:861–871.https://doi.org/10.1016/j.bbagen.2014.08.004
-
Gaussian accelerated molecular dynamics (gamd): principles and applicationsWiley Interdisciplinary Reviews. Computational Molecular Science 11:e1521.https://doi.org/10.1002/wcms.1521
-
Discovery of vu0467485/az13713945: an M4 PAM evaluated as a preclinical candidate for the treatment of schizophreniaACS Medicinal Chemistry Letters 8:233–238.https://doi.org/10.1021/acsmedchemlett.6b00461
-
Challenges in the development of an M4 PAM in vivo tool compound: the discovery of VU0467154 and unexpected DMPK profiles of close analogsBioorganic & Medicinal Chemistry Letters 27:171–175.https://doi.org/10.1016/j.bmcl.2016.11.086
-
Emerging paradigms in GPCR allostery: implications for drug discoveryNature Reviews. Drug Discovery 12:630–644.https://doi.org/10.1038/nrd4052
-
Gctf: real-time CTF determination and correctionJournal of Structural Biology 193:1–12.https://doi.org/10.1016/j.jsb.2015.11.003
Article and author information
Author details
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
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Biochemistry and Chemical Biology
- Structural Biology and Molecular Biophysics
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.
-
- Biochemistry and Chemical Biology
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.