Helix dipole movement and conformational variability contribute to allosteric GDP release in Galphai subunits

Recent advances in computational studies of GPCR-G protein interactions

Protein-protein interactions are key in cellular signaling. G protein-coupled receptors (GPCRs), the largest superfamily of human membrane proteins, are able to transduce extracellular signals (e.g., hormones and neurotransmitters) to intracellular proteins, in particular the G proteins. Since GPCRs serve as primary targets of ~1/3 of currently marketed drugs, it is important to understand mechanisms of GPCR signaling in order to design selective and potent drug molecules. This chapter focuses on recent advances in computational studies of the GPCR-G protein interactions using bioinformatics, protein-protein docking and molecular dynamics simulation approaches.

Dynamic regulation of GDP binding to G proteins revealed by magnetic field-dependent NMR relaxation analyses

Heterotrimeric guanine-nucleotide-binding proteins (G proteins) serve as molecular switches in signalling pathways, by coupling the activation of cell surface receptors to intracellular responses. Mutations in the G protein α-subunit (Gα) that accelerate guanosine diphosphate (GDP) dissociation cause hyperactivation of the downstream effector proteins, leading to oncogenesis. However, the structural mechanism of the accelerated GDP dissociation has remained unclear. Here, we use magnetic field-dependent nuclear magnetic resonance relaxation analyses to investigate the structural and dynamic properties of GDP bound Gα on a microsecond timescale. We show that Gα rapidly exchanges between a ground-state conformation, which tightly binds to GDP and an excited conformation with reduced GDP affinity. The oncogenic D150N mutation accelerates GDP dissociation by shifting the equilibrium towards the excited conformation.

Heterotrimeric guanine-nucleotide-binding proteins (G proteins) act as molecular switches. Here the authors use NMR relaxation analyses, which reveal the dynamics of G protein alpha subunit binding to GDP on a microsecond timescale.

Computational Simulation of the Activation Cycle of Gα Subunit in the G Protein Cycle Using an Elastic Network Model

Agonist-activated G protein-coupled receptors (GPCRs) interact with GDP-bound G protein heterotrimers (Gαβγ) promoting GDP/GTP exchange, which results in dissociation of Gα from the receptor and Gβγ. The GTPase activity of Gα hydrolyzes GTP to GDP, and the GDP-bound Gα interacts with Gβγ, forming a GDP-bound G protein heterotrimer. The G protein cycle is allosterically modulated by conformational changes of the Gα subunit. Although biochemical and biophysical methods have elucidated the structure and dynamics of Gα, the precise conformational mechanisms underlying the G protein cycle are not fully understood yet. Simulation methods could help to provide additional details to gain further insight into G protein signal transduction mechanisms. In this study, using the available X-ray crystal structures of Gα, we simulated the entire G protein cycle and described not only the steric features of the Gα structure, but also conformational changes at each step. Each reference structure in the G protein cycle was modeled as an elastic network model and subjected to normal mode analysis. Our simulation data suggests that activated receptors trigger conformational changes of the Gα subunit that are thermodynamically favorable for opening of the nucleotide-binding pocket and GDP release. Furthermore, the effects of GTP binding and hydrolysis on mobility changes of the C and N termini and switch regions are elucidated. In summary, our simulation results enabled us to provide detailed descriptions of the structural and dynamic features of the G protein cycle.

Dynamic Coupling and Allosteric Networks in the α Subunit of Heterotrimeric G Proteins

G protein α subunits cycle between active and inactive conformations to regulate a multitude of intracellular signaling cascades. Important structural transitions occurring during this cycle have been characterized from extensive crystallographic studies. However, the link between observed conformations and the allosteric regulation of binding events at distal sites critical for signaling through G proteins remain unclear. Here we describe molecular dynamics simulations, bioinformatics analysis, and experimental mutagenesis that identifies residues involved in mediating the allosteric coupling of receptor, nucleotide, and helical domain interfaces of Gαi. Most notably, we predict and characterize novel allosteric decoupling mutants, which display enhanced helical domain opening, increased rates of nucleotide exchange, and constitutive activity in the absence of receptor activation. Collectively, our results provide a framework for explaining how binding events and mutations can alter internal dynamic couplings critical for G protein function.

A conserved phenylalanine as a relay between the α5 helix and the GDP binding region of heterotrimeric Gi protein α subunit

G protein activation by G protein-coupled receptors is one of the critical steps for many cellular signal transduction pathways. Previously, we and other groups reported that the α5 helix in the G protein α subunit plays a major role during this activation process. However, the precise signaling pathway between the α5 helix and the guanosine diphosphate (GDP) binding pocket remains elusive. Here, using structural, biochemical, and computational techniques, we probed different residues around the α5 helix for their role in signaling. Our data showed that perturbing the Phe-336 residue disturbs hydrophobic interactions with the β2-β3 strands and α1 helix, leading to high basal nucleotide exchange. However, mutations in β strands β5 and β6 do not perturb G protein activation. We have highlighted critical residues that leverage Phe-336 as a relay. Conformational changes are transmitted starting from Phe-336 via β2-β3/α1 to Switch I and the phosphate binding loop, decreasing the stability of the GDP binding pocket and triggering nucleotide release. When the α1 and α5 helices were cross-linked, inhibiting the receptor-mediated displacement of the C-terminal α5 helix, mutation of Phe-336 still leads to high basal exchange rates. This suggests that unlike receptor-mediated activation, helix 5 rotation and translocation are not necessary for GDP release from the α subunit. Rather, destabilization of the backdoor region of the Gα subunit is sufficient for triggering the activation process.

A transient interaction between the phosphate binding loop and switch I contributes to the allosteric network between receptor and nucleotide in Gαi1

Receptor-mediated activation of the Gα subunit of heterotrimeric G proteins requires allosteric communication between the receptor binding site and the guanine nucleotide binding site, which are separated by >30 Å. Structural changes in the allosteric network connecting these sites are predicted to be transient in the wild-type Gα subunit, making studies of these connections challenging. In the current work, site-directed mutants that alter the energy barriers between the activation states are used as tools to better understand the transient features of allosteric signaling in the Gα subunit. The observed differences in relative receptor affinity for intact Gαi1 subunits versus C-terminal Gαi1 peptides harboring the K345L mutation are consistent with this mutation modulating the allosteric network in the protein subunit. Measurement of nucleotide exchange rates, affinity for metarhodopsin II, and thermostability suggest that the K345L Gαi1 variant has reduced stability in both the GDP-bound and nucleotide-free states as compared with wild type but similar stability in the GTPγS-bound state. High resolution x-ray crystal structures reveal conformational changes accompanying the destabilization of the GDP-bound state. Of these, the conformation for Switch I was stabilized by an ionic interaction with the phosphate binding loop. Further site-directed mutagenesis suggests that this interaction between Switch I and the phosphate binding loop is important for receptor-mediated nucleotide exchange in the wild-type Gαi1 subunit.

Energetic analysis of the rhodopsin-G-protein complex links the α5 helix to GDP release

We present a model of interaction of Gi protein with the activated receptor (R*) rhodopsin, which pinpoints energetic contributions to activation and reconciles the β2 adrenergic receptor-Gs crystal structure with new and previously published experimental data. In silico analysis demonstrated energetic changes when the Gα C-terminal helix (α5) interacts with the R* cytoplasmic pocket, thus leading to displacement of the helical domain and GDP release. The model features a less dramatic domain opening compared with the crystal structure. The α5 helix undergoes a 63° rotation, accompanied by a 5.7-Å translation, that reorganizes interfaces between α5 and α1 helices and between α5 and β6-α5. Changes in the β6-α5 loop displace αG. All of these movements lead to opening of the GDP-binding pocket. The model creates a roadmap for experimental studies of receptor-mediated G-protein activation.

Molecular basis of cannabinoid CB1 receptor coupling to the G protein heterotrimer Gαiβγ: identification of key CB1 contacts with the C-terminal helix α5 of Gαi

The cannabinoid (CB1) receptor is a member of the rhodopsin-like G protein-coupled receptor superfamily. The human CB1 receptor, which is among the most expressed receptors in the brain, has been implicated in several disease states, including drug addiction, anxiety, depression, obesity, and chronic pain. Different classes of CB1 agonists evoke signaling pathways through the activation of specific subtypes of G proteins. The molecular basis of CB1 receptor coupling to its cognate G protein is unknown. As a first step toward understanding CB1 receptor-mediated G protein signaling, we have constructed a ternary complex structural model of the CB1 receptor and Gi heterotrimer (CB1-Gi), guided by the x-ray structure of β2-adrenergic receptor (β2AR) in complex with Gs (β2AR-Gs), through 824-ns duration molecular dynamics simulations in a fully hydrated 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine bilayer environment. We identified a group of residues at the juxtamembrane regions of the intracellular loops 2 and 3 (IC2 and IC3) of the CB1 receptor, including Ile-218(3.54), Tyr-224(IC2), Asp-338(6.30), Arg-340(6.32), Leu-341(6.33), and Thr-344(6.36), as potential key contacts with the extreme C-terminal helix α5 of Gαi. Ala mutations of these residues at the receptor-Gi interface resulted in little G protein coupling activity, consistent with the present model of the CB1-Gi complex, which suggests tight interactions between CB1 and the extreme C-terminal helix α5 of Gαi. The model also suggests that unique conformational changes in the extreme C-terminal helix α5 of Gα play a crucial role in the receptor-mediated G protein activation.

Structural features of the G-protein/GPCR interactions

BACKGROUND

The details of the functional interaction between G proteins and the G protein coupled receptors (GPCRs) have long been subjected to extensive investigations with structural and functional assays and a large number of computational studies.

SCOPE OF REVIEW

The nature and sites of interaction in the G-protein/GPCR complexes, and the specificities of these interactions selecting coupling partners among the large number of families of GPCRs and G protein forms, are still poorly defined.

MAJOR CONCLUSIONS

Many of the contact sites between the two proteins in specific complexes have been identified, but the three dimensional molecular architecture of a receptor-Gα interface is only known for one pair. Consequently, many fundamental questions regarding this macromolecular assembly and its mechanism remain unanswered.

GENERAL SIGNIFICANCE

In the context of current structural data we review the structural details of the interfaces and recognition sites in complexes of sub-family A GPCRs with cognate G-proteins, with special emphasis on the consequences of activation on GPCR structure, the prevalence of preassembled GPCR/G-protein complexes, the key structural determinants for selective coupling and the possible involvement of GPCR oligomerization in this process.

Conformational flexibility and structural dynamics in GPCR-mediated G protein activation: a perspective

Structure and dynamics of G proteins and their cognate receptors, both alone and in complex, are becoming increasingly accessible to experimental techniques. Understanding the conformational changes and timelines that govern these changes can lead to new insights into the processes of ligand binding and associated G protein activation. Experimental systems may involve the use of, or otherwise stabilize, non-native environments. This can complicate our understanding of structural and dynamic features of processes such as the ionic lock, tryptophan toggle, and G protein flexibility. While elements in the receptor's transmembrane helices and the C-terminal α5 helix of Gα undergo well-defined structural changes, regions subject to conformational flexibility may be important in fine-tuning the interactions between activated receptors and G proteins. The pairing of computational and experimental approaches will continue to provide powerful tools to probe the conformation and dynamics of receptor-mediated G protein activation.

Linking receptor activation to changes in Sw I and II of Gα proteins

G-protein coupled receptors catalyze nucleotide exchange on G proteins, which results in subunit dissociation and effector activation. In the recent β2AR-Gs structure, portions of Switch I and II of Gα are not fully elucidated. We paired fluorescence studies of receptor-Gαi interactions with the β2AR-Gs and other Gi structures to investigate changes in Switch I and II during receptor activation and GTP binding. The β2/β3 loop containing Leu194 of Gαi is located between Switches I and II, in close proximity to IC2 of the receptor and the C-terminus of Gα, thus providing an allosteric connection between these Switches and receptor activation. We compared the environment of residues in myristoylated Gαi proteins in the heterotrimer to that upon receptor activation and subsequent GTP binding. Upon receptor activation, residues in both Switch regions are less solvent-exposed, as compared to the heterotrimer. Upon GTPγS binding, the environment of several residues in Switch I resemble the receptor-bound state, while Switch II residues display effects on their environment which are consistent with their role in GTP binding and Gβγ dissociation. The ability to merge available crystal structures with solution studies is a powerful tool to gain insight into conformational changes associated with receptor-mediated Gi protein activation.

Myristoylation exerts direct and allosteric effects on Gα conformation and dynamics in solution

Coupling of heterotrimeric G proteins to activated G protein-coupled receptors results in nucleotide exchange on the Gα subunit, which in turn decreases its affinity for both Gβγ and activated receptors. N-Terminal myristoylation of Gα subunits aids in membrane localization of inactive G proteins. Despite the presence of the covalently attached myristoyl group, Gα proteins are highly soluble after GTP binding. This study investigated factors facilitating the solubility of the activated, myristoylated protein. In doing so, we also identified myristoylation-dependent differences in regions of Gα known to play important roles in interactions with receptors, effectors, and nucleotide binding. Amide hydrogen-deuterium exchange and site-directed fluorescence of activated proteins revealed a solvent-protected amino terminus that was enhanced by myristoylation. Furthermore, fluorescence quenching confirmed that the myristoylated amino terminus is in proximity to the Switch II region in the activated protein. Myristoylation also stabilized the interaction between the guanine ring and the base of the α5 helix that contacts the bound nucleotide. The allosteric effects of myristoylation on protein structure, function, and localization indicate that the myristoylated amino terminus of Gα(i) functions as a myristoyl switch, with implications for myristoylation in the stabilization of nucleotide binding and in the spatial regulation of G protein signaling.

Allosteric mechanisms of G protein-Coupled Receptor signaling: a structural perspective

G protein-Coupled Receptors (GPCRs) use a complex series of intramolecular conformational changes to couple agonist binding to the binding and activation of cognate heterotrimeric G protein (Gαβγ). The mechanisms underlying this long-range activation have been identified using a variety of biochemical and structural approaches and have primarily used visual signal transduction via the GPCR rhodopsin and cognate heterotrimeric G protein transducin (G(t)) as a model system. In this chapter, we review the methods that have revealed allosteric signaling through rhodopsin and transducin. These methods can be applied to a variety of GPCR-mediated signaling pathways.

Interaction of a G protein with an activated receptor opens the interdomain interface in the alpha subunit

In G-protein signaling, an activated receptor catalyzes GDP/GTP exchange on the G(α) subunit of a heterotrimeric G protein. In an initial step, receptor interaction with G(α) acts to allosterically trigger GDP release from a binding site located between the nucleotide binding domain and a helical domain, but the molecular mechanism is unknown. In this study, site-directed spin labeling and double electron-electron resonance spectroscopy are employed to reveal a large-scale separation of the domains that provides a direct pathway for nucleotide escape. Cross-linking studies show that the domain separation is required for receptor enhancement of nucleotide exchange rates. The interdomain opening is coupled to receptor binding via the C-terminal helix of G(α), the extension of which is a high-affinity receptor binding element.

Coupling efficiency of rhodopsin and transducin in bicelles

G protein coupled receptors (GPCRs) can be activated by various extracellular stimuli, including hormones, peptides, odorants, neurotransmitters, nucleotides, or light. After activation, receptors interact with heterotrimeric G proteins and catalyze GDP release from the Gα subunit, the rate limiting step in G protein activation, to form a high affinity nucleotide-free GPCR-G protein complex. In vivo, subsequent GTP binding reduces affinity of the Gα protein for the activated receptor. In this study, we investigated the biochemical and structural characteristics of the prototypical GPCR, rhodopsin, and its signaling partner, transducin (G(t)), in bicelles to better understand the effects of membrane composition on high affinity complex formation, stability, and receptor mediated nucleotide release. Our results demonstrate that the high-affinity complex (rhodopsin-G(t)(empty)) forms more readily and has dramatically increased stability when rhodopsin is integrated into bicelles of a defined composition. We increased the half-life of functional complex to 1 week in the presence of negatively charged phospholipids. These data suggest that a membrane-like structure is an important contributor to the formation and stability of functional receptor-G protein complexes and can extend the range of studies that investigate properties of these complexes.

Trp fluorescence reveals an activation-dependent cation-pi interaction in the Switch II region of Galphai proteins

Crystal structures of Galpha(i) (and closely related family member Galpha(t)) reveal much of what we currently know about G protein structure, including changes which occur in Switch regions. Galpha(t) exhibits a low rate of basal (uncatalyzed) nucleotide exchange and an ordered Switch II region in the GDP-bound state, unlike Galpha(i), which exhibits higher basal exchange and a disordered Switch II region in Galpha(i)GDP structures. Using purified Galpha(i) and Galpha(t), we examined the intrinsic tryptophan fluorescence of these proteins, which reports conformational changes associated with activation and deactivation of Galpha proteins. In addition to the expected enhancement in tryptophan fluorescence intensity, activation of GalphaGDP proteins was accompanied by a modest but notable red shift in tryptophan emission maxima. We identified a cation-pi interaction between tryptophan and arginine residues in the Switch II of Galpha(i) family proteins that mediates the observed red shift in emission maxima. Furthermore, amino-terminal myristoylation of Galpha(i) resulted in a less polar environment for tryptophan residues in the GTPase domain, consistent with an interaction between the myristoylated amino terminus and the GTPase domain of Galpha proteins. These results reveal unique insights into conformational changes which occur upon activation and deactivation of G proteins in solution.

A G protein-coupled receptor at work: the rhodopsin model

G protein-coupled receptors (GPCRs) are ubiquitous signal transducers in cell membranes, as well as important drug targets. Interaction with extracellular agonists turns the seven transmembrane helix (7TM) scaffold of a GPCR into a catalyst for GDP and GTP exchange in heterotrimeric Galphabetagamma proteins. Activation of the model GPCR, rhodopsin, is triggered by photoisomerization of its retinal ligand. From the augmentation of biochemical and biophysical studies by recent high-resolution 3D structures, its activation intermediates can now be interpreted as the stepwise engagement of protein domains. Rearrangement of TM5-TM6 opens a crevice at the cytoplasmic side of the receptor into which the C terminus of the Galpha subunit can bind. The Galpha C-terminal helix is used as a transmission rod to the nucleotide binding site. The mechanism relies on dynamic interactions between conserved residues and could therefore be common to other GPCRs.

Structural evidence for a sequential release mechanism for activation of heterotrimeric G proteins

Heptahelical G-protein (heterotrimeric guanine nucleotide-binding protein)-coupled receptors couple to heterotrimeric G proteins to relay extracellular signals to intracellular signaling networks, but the molecular mechanism underlying guanosine 5'-diphosphate (GDP) release by the G protein alpha-subunit is not well understood. Amino acid substitutions in the conserved alpha5 helix of G(i), which extends from the C-terminal region to the nucleotide-binding pocket, cause dramatic increases in basal (receptor-independent) GDP release rates. For example, mutant Galpha(i1)-T329A shows an 18-fold increase in basal GDP release rate and, when expressed in culture, it causes a significant decrease in forskolin-stimulated cAMP accumulation. The crystal structure of Galpha(i1)-T329A.GDP shows substantial conformational rearrangement of the switch I region and additional striking alterations of side chains lining the catalytic pocket that disrupt the Mg(+2) coordination sphere and dislodge bound Mg(+2). We propose a "sequential release" mechanism whereby a transient conformational change in the alpha5 helix alters switch I to induce GDP release. Interestingly, this mechanistic model for heterotrimeric G protein activation is similar to that suggested for the activation of the plant small G protein Rop4 by RopGEF8.

Structural and kinetic modeling of an activating helix switch in the rhodopsin-transducin interface

Extracellular signals prompt G protein-coupled receptors (GPCRs) to adopt an active conformation (R*) and catalyze GDP/GTP exchange in the alpha-subunit of intracellular G proteins (Galphabetagamma). Kinetic analysis of transducin (G(t)alphabetagamma) activation shows that an intermediary R*xG(t)alphabetagamma.GDP complex is formed that precedes GDP release and formation of the nucleotide-free R*xG protein complex. Based on this reaction sequence, we explore the dynamic interface between the proteins during formation of these complexes. We start from the R* conformation stabilized by a G(t)alpha C-terminal peptide (GalphaCT) obtained from crystal structures of the GPCR opsin. Molecular modeling allows reconstruction of the fully elongated C-terminal alpha-helix of G(t)alpha (alpha5) and shows how alpha5 can be docked to the open binding site of R*. Two modes of interaction are found. One of them--termed stable or S-interaction--matches the position of the GalphaCT peptide in the crystal structure and reproduces the hydrogen-bonding networks between the C-terminal reverse turn of GalphaCT and conserved E(D)RY and NPxxY(x)(5,6)F regions of the GPCR. The alternative fit--termed intermediary or I-interaction--is distinguished by a tilt (42 degrees ) and rotation (90 degrees ) of alpha5 relative to the S-interaction and shows different alpha5 contacts with the NPxxY(x)(5,6)F region and the second cytoplasmic loop of R*. From the 2 alpha5 interactions, we derive a "helix switch" mechanism for the transition of R*xG(t)alphabetagamma.GDP to the nucleotide-free R*xG protein complex that illustrates how alpha5 might act as a transmission rod to propagate the conformational change from the receptor-G protein interface to the nucleotide binding site.