Time-resolved cryo-EM of G-protein activation by a GPCR

Relevance of G protein-coupled receptor (GPCR) dynamics for receptor activation, signalling bias and allosteric modulation

G protein-coupled receptors (GPCRs) are one of the major drug targets. In recent years, computational drug design for GPCRs has mainly focused on static structures obtained through X-ray crystallography, cryogenic electron microscopy (cryo-EM) or in silico modelling as a starting point for virtual screening campaigns. However, GPCRs are highly flexible entities with the ability to adopt different conformational states that elicit different physiological responses. Including this knowledge in the drug discovery pipeline can help to tailor novel conformation-specific drugs with an improved therapeutic profile. In this review, we outline our current knowledge about GPCR dynamics that is relevant for receptor activation, signalling bias and allosteric modulation. Ultimately, we highlight new technological implementations such as time-resolved X-ray crystallography and cryo-EM as well as computational algorithms that can contribute to a more comprehensive understanding of receptor dynamics and its relevance for GPCR functionality.

Ligand-Receptor Interactions and Structure-Function Relationships in Off-Target Binding of the β3-Adrenergic Agonist Mirabegron to α1A-Adrenergic Receptors

The β3-adrenoceptor agonist mirabegron is available for the treatment of storage symptoms of overactive bladder, including frequency, urgency, and incontinence. The off-target effects of mirabegron include binding to α1-adrenoceptors, which are central in the treatment of voiding symptoms. Here, we examined the structure-function relationships in the binding of mirabegron to a cryo-electron microscopy structure of α1A. The binding was simulated by docking mirabegron to a 3D structure of a human α1A-adrenoceptor (7YMH) using Autodock Vina. The simulations identified two binding states: slope orientation involving 10 positions and horizontal binding to the receptor surface involving 4 positions. No interactions occurred with positions constituting the α1A binding pocket, including Asp-106, Ser-188, or Phe-312, despite the positioning of the phenylethanolamine moiety in transmembrane regions close to the binding pocket by contact with Phe-288, -289, and Val-107. Contact with the unique positions of α1A included the transmembrane Met-292 during slope binding and exosite Phe-86 during horizontal binding. Exosite binding in slope orientation involved contact of the anilino part, rather than the aminothiazol end, to Ile-178, Ala-103, and Asn-179. In conclusion, contact with Met-292 and Phe-86, which are unique positions of α1A, accounts for mirabegron binding to α1A. Because of its lack of interactions with the binding pocket, mirabegron has lower affinity compared to α1A-blockers and no effects on voiding symptoms.

Structural and Computational Insights into Dynamics and Intermediate States of Orexin 2 Receptor Signaling

Orexin 2 receptor (OX2R) is a G protein-coupled receptor (GPCR) whose activation is crucial to regulation of the sleep-wake cycle. Recently, inactive and active state structures were determined from X-ray crystallography and cryo-electron microscopy single particle analysis, and the activation mechanisms have been discussed based on these static data. GPCRs have multiscale intermediate states during activation, and insights into these dynamics and intermediate states may aid the precise control of intracellular signaling by ligands in drug discovery. Molecular dynamics (MD) simulations are used to investigate dynamics induced in response to thermal perturbations, such as structural fluctuations of main and side chains. In this study, we proposed collective motions of the TM domain during activation by performing 30 independent microsecond-scale MD simulations for various OX2R systems and applying relaxation mode analysis. The analysis results suggested that TM3 had a vertical structural movement relative to the membrane surface during activation. In addition, we extracted three characteristic amino acid residues on TM3, i.e., Q1343.32, V1423.40, and R1523.50, which exhibited large conformational fluctuations. We quantitatively evaluated the changes in their equilibrium during activation in relation to the movement of TM3. We also discuss the regulation of ligand binding recognition and intracellular signal selectivity by changes in the equilibrium of Q1343.32 and R1523.50, respectively, according to MD simulations and GPCR database. Additionally, the OX2R-Gi signaling complex is stabilized in the conformation resembling a non-canonical (NC) state, which was previously proposed as an intermediate state during activation of neurotensin 1 receptor. Insights into the dynamics and intermediate states during activation gained from this study may be useful for developing biased agonists for OX2R.

Multiple recent HCAR2 structures demonstrate a highly dynamic ligand binding and G protein activation mode

A surprisingly clear picture of the allosteric mechanism connecting G protein-coupled receptor agonists with G protein binding-and back - is revealed by a puzzle of thirty novel 3D structures of the hydroxycarboxylic acid receptor 2 (HCAR2) in complex with eight different orthosteric and a single allosteric agonist. HCAR2 is a sensor of β-hydroxybutyrate, niacin and certain anti-inflammatory drugs. Surprisingly, agonists with and without on-target side effects bound very similarly and in a completely occluded orthosteric binding site. Thus, despite the many structures we are still left with a pertinent need to understand the molecular dynamics of this and similar systems.

Mechanistic insights into G-protein coupling with an agonist-bound G-protein-coupled receptor

G-protein-coupled receptors (GPCRs) activate heterotrimeric G proteins by promoting guanine nucleotide exchange. Here, we investigate the coupling of G proteins with GPCRs and describe the events that ultimately lead to the ejection of GDP from its binding pocket in the Gα subunit, the rate-limiting step during G-protein activation. Using molecular dynamics simulations, we investigate the temporal progression of structural rearrangements of GDP-bound Gs protein (Gs·GDP; hereafter GsGDP) upon coupling to the β2-adrenergic receptor (β2AR) in atomic detail. The binding of GsGDP to the β2AR is followed by long-range allosteric effects that significantly reduce the energy needed for GDP release: the opening of α1-αF helices, the displacement of the αG helix and the opening of the α-helical domain. Signal propagation to the Gs occurs through an extended receptor interface, including a lysine-rich motif at the intracellular end of a kinked transmembrane helix 6, which was confirmed by site-directed mutagenesis and functional assays. From this β2AR-GsGDP intermediate, Gs undergoes an in-plane rotation along the receptor axis to approach the β2AR-Gsempty state. The simulations shed light on how the structural elements at the receptor-G-protein interface may interact to transmit the signal over 30 Å to the nucleotide-binding site. Our analysis extends the current limited view of nucleotide-free snapshots to include additional states and structural features responsible for signaling and G-protein coupling specificity.

Structure and function of an intermediate GPCR-Gαβγ complex

Unraveling the signaling roles of intermediate complexes is pivotal for G protein-coupled receptor (GPCR) drug development. Despite hundreds of GPCR-Gαβγ structures, these snapshots primarily capture the fully activated end-state complex. Consequently, a comprehensive understanding of the conformational transitions during GPCR activation and the roles of intermediate GPCR-G protein complexes in signaling remain elusive. Guided by a conformational landscape profiled by 19 F quantitative NMR ( 19 F-qNMR) and Molecular Dynamics (MD) simulations, we resolved the structure of an unliganded GPCR-G protein intermediate complex by blocking its transition to the fully activated end-state complex. More importantly, we presented direct evidence that the intermediate GPCR-Gαsβγ complex initiates a rate-limited nucleotide exchange without progressing to the fully activated end-state complex, thereby bridging a significant gap in our understanding the complexity of GPCR signaling. Understanding the roles of individual conformational states and their complexes in signaling efficacy and bias will help us to design drugs that discriminately target a disease-related conformation.