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

The role of G protein conformation in receptor-G protein selectivity

G protein-coupled receptors (GPCRs) selectively activate at least one of the four families of heterotrimeric G proteins, but the mechanism of coupling selectivity remains unclear. Structural studies emphasize structural complementarity of GPCRs and nucleotide-free G proteins, but selectivity is likely to be determined by transient intermediate-state complexes that exist before nucleotide release. Here we study coupling to nucleotide-decoupled G protein variants that can adopt conformations similar to receptor-bound G proteins without releasing nucleotide, and are therefore able to bypass intermediate-state complexes. We find that selectivity is degraded when nucleotide release is not required for GPCR-G protein complex formation, to the extent that most GPCRs interact with most nucleotide-decoupled G proteins. These findings demonstrate the absence of absolute structural incompatibility between noncognate receptor-G protein pairs, and are consistent with the hypothesis that transient intermediate states are partly responsible for coupling selectivity.

Structure and function of the metagenomic plastic-degrading polyester hydrolase PHL7 bound to its product

The recently discovered metagenomic-derived polyester hydrolase PHL7 is able to efficiently degrade amorphous polyethylene terephthalate (PET) in post-consumer plastic waste. We present the cocrystal structure of this hydrolase with its hydrolysis product terephthalic acid and elucidate the influence of 17 single mutations on the PET-hydrolytic activity and thermal stability of PHL7. The substrate-binding mode of terephthalic acid is similar to that of the thermophilic polyester hydrolase LCC and deviates from the mesophilic IsPETase. The subsite I modifications L93F and Q95Y, derived from LCC, increased the thermal stability, while exchange of H185S, derived from IsPETase, reduced the stability of PHL7. The subsite II residue H130 is suggested to represent an adaptation for high thermal stability, whereas L210 emerged as the main contributor to the observed high PET-hydrolytic activity. Variant L210T showed significantly higher activity, achieving a degradation rate of 20 µm h^-1 with amorphous PET films.

Low Carbon Footprint Recycling of Post-Consumer PET Plastic with a Metagenomic Polyester Hydrolase

Earth is flooded with plastics and the need for sustainable recycling strategies for polymers has become increasingly urgent. Enzyme-based hydrolysis of post-consumer plastic is an emerging strategy for closed-loop recycling of polyethylene terephthalate (PET). The polyester hydrolase PHL7, isolated from a compost metagenome, completely hydrolyzes amorphous PET films, releasing 91 mg of terephthalic acid per hour and mg of enzyme. Vertical scanning interferometry shows degradation rates of the PET film of 6.8 μm h-1 . Structural analysis indicates the importance of leucine at position 210 for the extraordinarily high PET-hydrolyzing activity of PHL7. Within 24 h, 0.6 mgenzyme  gPET -1 completely degrades post-consumer thermoform PET packaging in an aqueous buffer at 70 °C without any energy-intensive pretreatments. Terephthalic acid recovered from the enzymatic hydrolysate is then used to synthesize virgin PET, demonstrating the potential of polyester hydrolases as catalysts in sustainable PET recycling processes with a low carbon footprint.

Predicted Structure of Fully Activated Tas1R3/1R3' Homodimer Bound to G Protein and Natural Sugars: Structural Insights into G Protein Activation by a Class C Sweet Taste Homodimer with Natural Sugars

The Tas1R3 G protein-coupled receptor constitutes the main component of sweet taste sensory response in humans via forming a heterodimer with Tas1R2 or a homodimer with Tas1R3. The Tas1R3/1R3' homodimer serves as a low-affinity sweet taste receptor, stimulating gustducin G protein (GGust) signaling in the presence of a high concentration of natural sugars. This provides an additional means to detect the taste of natural sugars, thereby differentiating the flavors between natural sugars and artificial sweeteners. We report here the predicted 3D structure of active state Tas1R3/1R3' homodimer complexed with heterotrimeric GGust and sucrose. We discovered that the GGust makes ionic anchors to intracellular loops 1 and 2 of Tas1R3 while the Gα-α5 helix engages the cytoplasmic region extensively through salt bridge and hydrophobic interactions. We show that in the activation of this complex the Venus flytrap domains of the homodimer undergo a remarkable twist up to ∼100° rotation around the vertical axis to adopt a closed-closed conformation while the intracellular region relaxes to an open-open conformation. We find that binding of sucrose to the homodimer stabilizes a preactivated conformation with a largely open intracellular region that recruits and activates the GGust. Upon activation, the Gα subunit spontaneously opens up the nucleotide-binding site, making nucleotide exchange facile for signaling. This activation of GGust promotes the interdomain twist of the Venus flytrap domains. These structures and transformations could potentially be a basis for the design of new sweeteners with higher activity and less unpleasant flavors.

Conformational switch that induces GDP release from Gi

Heterotrimeric guanine nucleotide-binding proteins (G proteins) are composed of α, β, and γ subunits. Gα switches between guanosine diphosphate (GDP)-bound inactive and guanosine triphosphate (GTP)-bound active states, and Gβγ interacts with the GDP-bound state. The GDP-binding regions are composed of two sites: the phosphate-binding and guanine-binding regions. The turnover of GDP and GTP is induced by guanine nucleotide-exchange factors (GEFs), including G protein-coupled receptors (GPCRs), Ric8A, and GIV/Girdin. However, the key structural factors for stabilizing the GDP-bound state of G proteins and the direct structural event for GDP release remain unclear. In this study, we investigated structural factors affecting GDP release by introducing point mutations in selected, conserved residues in Gαi3. We examined the effects of these mutations on the GDP/GTP turnover rate and the overall conformation of Gαi3 as well as the binding free energy between Gαi3 and GDP. We found that dynamic changes in the phosphate-binding regions are an immediate factor for the release of GDP.

G protein-coupled receptor-G protein interactions: a single-molecule perspective

G protein-coupled receptors (GPCRs) regulate many cellular and physiological processes, responding to a diverse range of extracellular stimuli including hormones, neurotransmitters, odorants, and light. Decades of biochemical and pharmacological studies have provided fundamental insights into the mechanisms of GPCR signaling. Thanks to recent advances in structural biology, we now possess an atomistic understanding of receptor activation and G protein coupling. However, how GPCRs and G proteins interact in living cells to confer signaling efficiency and specificity remains insufficiently understood. The development of advanced optical methods, including single-molecule microscopy, has provided the means to study receptors and G proteins in living cells with unprecedented spatio-temporal resolution. The results of these studies reveal an unexpected level of complexity, whereby GPCRs undergo transient interactions among themselves as well as with G proteins and structural elements of the plasma membrane to form short-lived signaling nanodomains that likely confer both rapidity and specificity to GPCR signaling. These findings may provide new strategies to pharmaceutically modulate GPCR function, which might eventually pave the way to innovative drugs for common diseases such as diabetes or heart failure.

Internal water channel formation in CXCR4 is crucial for Gi-protein coupling upon activation by CXCL12

Chemokine receptor CXCR4 is a major drug target for numerous diseases because of its involvement in the regulation of cell migration and the developmental process. In this study, atomic-level molecular dynamics simulations were used to determine the activation mechanism and internal water formation of CXCR4 in complex with chemokine CXCL12 and G_i-protein. The results indicated that CXCL12-bound CXCR4 underwent transmembrane 6 (TM6) outward movement and a decrease in tyrosine toggle switch by eliciting the breakage of hydrophobic layer to form a continuous internal water channel. In the GDP-bound G_αi-protein state, the rotation and translation of the α5-helix of G_αi-protein toward the cytoplasmic pocket of CXCR4 induced an increase in interdomain distance for GDP leaving. Finally, an internal water channel formation model was proposed based on our simulations for CXCL12-bound CXCR4 in complex with G_αi-protein upon activation for downstream signaling. This model could be useful in anticancer drug development.

Rosetta Machine Learning Models Accurately Classify Positional Effects of Thioamides on Proteolysis

Thioamide substitutions of the peptide backbone have been shown to stabilize therapeutic and imaging peptides toward proteolysis. In order to rationally design thioamide modifications, we have developed a novel Rosetta custom score function to classify thioamide positional effects on proteolysis in substrates of serine and cysteine proteases. Peptides of interest were docked into proteases using the FlexPepDock application in Rosetta. Docked complexes were modified to contain thioamides parametrized through the creation of custom atom types in Rosetta based on ab intio simulations. Thioamide complexes were simulated, and the resultant structural complexes provided features for machine learning classification as the decomposed values of the Rosetta score function. An ensemble, majority voting model was developed to be a robust predictor of previously unpublished thioamide proteolysis holdout data. Theoretical control simulations with pseudo-atoms that modulate only one physical characteristic of the thioamide show differential effects on prediction accuracy by the optimized voting classification model. These pseudo-atom model simulations, as well as statistical analyses of the full thioamide simulations, implicate steric effects on peptide binding as being primarily responsible for thioamide positional effects on proteolytic resistance.

Rapid Simulation of Unprocessed DEER Decay Data for Protein Fold Prediction

Despite advances in sampling and scoring strategies, Monte Carlo modeling methods still struggle to accurately predict de novo the structures of large proteins, membrane proteins, or proteins of complex topologies. Previous approaches have addressed these shortcomings by leveraging sparse distance data gathered using site-directed spin labeling and electron paramagnetic resonance spectroscopy to improve protein structure prediction and refinement outcomes. However, existing computational implementations entail compromises between coarse-grained models of the spin label that lower the resolution and explicit models that lead to resource-intense simulations. These methods are further limited by their reliance on distance distributions, which are calculated from a primary refocused echo decay signal and contain uncertainties that may require manual refinement. Here, we addressed these challenges by developing RosettaDEER, a scoring method within the Rosetta software suite capable of simulating double electron-electron resonance spectroscopy decay traces and distance distributions between spin labels fast enough to fold proteins de novo. We demonstrate that the accuracy of resulting distance distributions match or exceed those generated by more computationally intensive methods. Moreover, decay traces generated from these distributions recapitulate intermolecular background coupling parameters even when the time window of data collection is truncated. As a result, RosettaDEER can discriminate between poorly folded and native-like models by using decay traces that cannot be accurately converted into distance distributions using regularized fitting approaches. Finally, using two challenging test cases, we demonstrate that RosettaDEER leverages these experimental data for protein fold prediction more effectively than previous methods. These benchmarking results confirm that RosettaDEER can effectively leverage sparse experimental data for a wide array of modeling applications built into the Rosetta software suite.

Assembly of a GPCR-G Protein Complex

The activation of G proteins by G protein-coupled receptors (GPCRs) underlies the majority of transmembrane signaling by hormones and neurotransmitters. Recent structures of GPCR-G protein complexes obtained by crystallography and cryoelectron microscopy (cryo-EM) reveal similar interactions between GPCRs and the alpha subunit of different G protein isoforms. While some G protein subtype-specific differences are observed, there is no clear structural explanation for G protein subtype-selectivity. All of these complexes are stabilized in the nucleotide-free state, a condition that does not exist in living cells. In an effort to better understand the structural basis of coupling specificity, we used time-resolved structural mass spectrometry techniques to investigate GPCR-G protein complex formation and G-protein activation. Our results suggest that coupling specificity is determined by one or more transient intermediate states that serve as selectivity filters and precede the formation of the stable nucleotide-free GPCR-G protein complexes observed in crystal and cryo-EM structures.

Active Site Flexibility as a Hallmark for Efficient PET Degradation by I. sakaiensis PETase

Polyethylene terephthalate (PET) is one of the most-consumed synthetic polymers, with an annual production of 50 million tons. Unfortunately, PET accumulates as waste and is highly resistant to biodegradation. Recently, fungal and bacterial thermophilic hydrolases were found to catalyze PET hydrolysis with optimal activities at high temperatures. Strikingly, an enzyme from Ideonella sakaiensis, termed PETase, was described to efficiently degrade PET at room temperature, but the molecular basis of its activity is not currently understood. Here, a crystal structure of PETase was determined at 2.02 Å resolution and employed in molecular dynamics simulations showing that the active site of PETase has higher flexibility at room temperature than its thermophilic counterparts. This flexibility is controlled by a novel disulfide bond in its active site, with its removal leading to destabilization of the catalytic triad and reduction of the hydrolase activity. Molecular docking of a model substrate predicts that PET binds to PETase in a unique and energetically favorable conformation facilitated by several residue substitutions within its active site when compared to other enzymes. These computational predictions are in excellent agreement with recent mutagenesis and PET film degradation analyses. Finally, we rationalize the increased catalytic activity of PETase at room temperature through molecular dynamics simulations of enzyme-ligand complexes for PETase and other thermophilic PET-degrading enzymes at 298, 323, and 353 K. Our results reveal that both the binding pose and residue substitutions within PETase favor proximity between the catalytic residues and the labile carbonyl of the substrate at room temperature, suggesting a more favorable hydrolytic reaction. These results are valuable for enabling detailed evolutionary analysis of PET-degrading enzymes and for rational design endeavors aiming at increasing the efficiency of PETase and similar enzymes toward plastic degradation.

A New Multisystem Disorder Caused by the Gαs Mutation p.F376V

CONTEXT

The α subunit of the stimulatory G protein (Gαs) links numerous receptors to adenylyl cyclase. Gαs, encoded by GNAS, is expressed predominantly from the maternal allele in certain tissues. Thus, maternal heterozygous loss-of-function mutations cause hormonal resistance, as in pseudohypoparathyroidism type Ia, whereas somatic gain-of-function mutations cause hormone-independent endocrine stimulation, as in McCune-Albright syndrome.

OBJECTIVE

We report two unrelated boys presenting with a new combination of clinical findings that suggest both gain and loss of Gαs function.

DESIGN AND SETTING

Clinical features were studied and sequencing of GNAS was performed. Signaling capacities of wild-type and mutant Gαs were determined in the presence of different G protein-coupled receptors (GPCRs) under basal and agonist-stimulated conditions.

RESULTS

Both unrelated patients presented with unexplained hyponatremia in infancy, followed by severe early onset gonadotrophin-independent precocious puberty and skeletal abnormalities. An identical heterozygous de novo variant (c.1136T>G; p.F376V) was found on the maternal GNAS allele in both patients; this resulted in a clinical phenotype that differed from known Gαs-related diseases and suggested gain of function at the vasopressin 2 receptor (V2R) and lutropin/choriogonadotropin receptor (LHCGR), yet increased serum PTH concentrations indicative of impaired proximal tubular PTH1 receptor (PTH1R) function. In vitro studies demonstrated that Gαs-F376V enhanced ligand-independent signaling at the PTH1R, LHCGR, and V2R and, at the same time, blunted ligand-dependent responses. Structural homology modeling suggested mutation-induced modifications at the C-terminal α5 helix of Gαs that are relevant for interaction with GPCRs and signal transduction.

CONCLUSIONS

The Gαs p.F376V mutation causes a previously unrecognized multisystem disorder.

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.

Comparative structural dynamic analysis of GTPases

GTPases regulate a multitude of essential cellular processes ranging from movement and division to differentiation and neuronal activity. These ubiquitous enzymes operate by hydrolyzing GTP to GDP with associated conformational changes that modulate affinity for family-specific binding partners. There are three major GTPase superfamilies: Ras-like GTPases, heterotrimeric G proteins and protein-synthesizing GTPases. Although they contain similar nucleotide-binding sites, the detailed mechanisms by which these structurally and functionally diverse superfamilies operate remain unclear. Here we compare and contrast the structural dynamic mechanisms of each superfamily using extensive molecular dynamics (MD) simulations and subsequent network analysis approaches. In particular, dissection of the cross-correlations of atomic displacements in both the GTP and GDP-bound states of Ras, transducin and elongation factor EF-Tu reveals analogous dynamic features. This includes similar dynamic communities and subdomain structures (termed lobes). For all three proteins the GTP-bound state has stronger couplings between equivalent lobes. Network analysis further identifies common and family-specific residues mediating the state-specific coupling of distal functional sites. Mutational simulations demonstrate how disrupting these couplings leads to distal dynamic effects at the nucleotide-binding site of each family. Collectively our studies extend current understanding of GTPase allosteric mechanisms and highlight previously unappreciated similarities across functionally diverse families.

GTPases are a large superfamily of essential enzymes that regulate a variety of cellular processes. They share a common core structure supporting nucleotide binding and hydrolysis, and are potentially descended from the same ancestor. Yet their biological functions diverge dramatically, ranging from cell division and movement to signal transduction and translation. It has been shown that conformational changes through binding to different substrates underlie the regulation of their activities. Here we investigate the conformational dynamics of three typical GTPases by in silico simulation. We find that these three GTPases possess overall similar substrate-associated dynamic features, beyond their distinct functions. Further identification of key common and family-specific elements in these three families helps us understand how enzymes are adapted to acquire distinct functions from a common core structure. Our results provide unprecedented insights into the functional mechanism of GTPases in general, which potentially facilitates novel protein design in the future.

Tracing the evolution of the heterotrimeric G protein α subunit in Metazoa

BACKGROUND

Heterotrimeric G proteins are fundamental signaling proteins composed of three subunits, Gα and a Gβγ dimer. The role of Gα as a molecular switch is critical for transmitting and amplifying intracellular signaling cascades initiated by an activated G protein Coupled Receptor (GPCR). Despite their biochemical and therapeutic importance, the study of G protein evolution has been limited to the scope of a few model organisms. Furthermore, of the five primary Gα subfamilies, the underlying gene structure of only two families has been thoroughly investigated outside of Mammalia evolution. Therefore our understanding of Gα emergence and evolution across phylogeny remains incomplete.

RESULTS

We have computationally identified the presence and absence of every Gα gene (GNA-) across all major branches of Deuterostomia and evaluated the conservation of the underlying exon-intron structures across these phylogenetic groups. We provide evidence of mutually exclusive exon inclusion through alternative splicing in specific lineages. Variations of splice site conservation and isoforms were found for several paralogs which coincide with conserved, putative motifs of DNA-/RNA-binding proteins. In addition to our curated gene annotations, within Primates, we identified 15 retrotranspositions, many of which have undergone pseudogenization. Most importantly, we find numerous deviations from previous findings regarding the presence and absence of individual GNA- genes, nuanced differences in phyla-specific gene copy numbers, novel paralog duplications and subsequent intron gain and loss events.

CONCLUSIONS

Our curated annotations allow us to draw more accurate inferences regarding the emergence of all Gα family members across Metazoa and to present a new, updated theory of Gα evolution. Leveraging this, our results are critical for gaining new insights into the co-evolution of the Gα subunit and its many protein binding partners, especially therapeutically relevant G protein - GPCR signaling pathways which radiated in Vertebrata evolution.

Integrated Structural Biology for α-Helical Membrane Protein Structure Determination

While great progress has been made, only 10% of the nearly 1,000 integral, α-helical, multi-span membrane protein families are represented by at least one experimentally determined structure in the PDB. Previously, we developed the algorithm BCL::MP-Fold, which samples the large conformational space of membrane proteins de novo by assembling predicted secondary structure elements guided by knowledge-based potentials. Here, we present a case study of rhodopsin fold determination by integrating sparse and/or low-resolution restraints from multiple experimental techniques including electron microscopy, electron paramagnetic resonance spectroscopy, and nuclear magnetic resonance spectroscopy. Simultaneous incorporation of orthogonal experimental restraints not only significantly improved the sampling accuracy but also allowed identification of the correct fold, which is demonstrated by a protein size-normalized transmembrane root-mean-square deviation as low as 1.2 Å. The protocol developed in this case study can be used for the determination of unknown membrane protein folds when limited experimental restraints are available.

Improved in Vitro Folding of the Y2 G Protein-Coupled Receptor into Bicelles

Prerequisite for structural studies on G protein-coupled receptors is the preparation of highly concentrated, stable, and biologically active receptor samples in milligram amounts of protein. Here, we present an improved protocol for Escherichia coli expression, functional refolding, and reconstitution into bicelles of the human neuropeptide Y receptor type 2 (Y₂R) for solution and solid-state NMR experiments. The isotopically labeled receptor is expressed in inclusion bodies and purified using SDS. We studied the details of an improved preparation protocol including the in vitro folding of the receptor, e.g., the native disulfide bridge formation, the exchange of the denaturating detergent SDS, and the functional reconstitution into bicelle environments of varying size. Full pharmacological functionality of the Y₂R preparation was shown by a ligand affinity of 4 nM and G-protein activation. Further, simple NMR experiments are used to test sample quality in high micromolar concentration.

Structure and dynamics of GPCR signaling complexes

G-protein-coupled receptors (GPCRs) relay numerous extracellular signals by triggering intracellular signaling through coupling with G proteins and arrestins. Recent breakthroughs in the structural determination of GPCRs and GPCR-transducer complexes represent important steps toward deciphering GPCR signal transduction at a molecular level. A full understanding of the molecular basis of GPCR-mediated signaling requires elucidation of the dynamics of receptors and their transducer complexes as well as their energy landscapes and conformational transition rates. Here, we summarize current insights into the structural plasticity of GPCR-G-protein and GPCR-arrestin complexes that underlies the regulation of the receptor's intracellular signaling profile.

Molecular Mechanisms of GPCR Signaling: A Structural Perspective

G protein-coupled receptors (GPCRs) are cell surface receptors that respond to a wide variety of stimuli, from light, odorants, hormones, and neurotransmitters to proteins and extracellular calcium. GPCRs represent the largest family of signaling proteins targeted by many clinically used drugs. Recent studies shed light on the conformational changes that accompany GPCR activation and the structural state of the receptor necessary for the interactions with the three classes of proteins that preferentially bind active GPCRs, G proteins, G protein-coupled receptor kinases (GRKs), and arrestins. Importantly, structural and biophysical studies also revealed activation-related conformational changes in these three types of signal transducers. Here, we summarize what is already known and point out questions that still need to be answered. Clear understanding of the structural basis of signaling by GPCRs and their interaction partners would pave the way to designing signaling-biased proteins with scientific and therapeutic potential.

Structural dynamics of Giα protein revealed by single molecule FRET

The heterotrimeric G proteins (Gαβγ) act as molecular switches to mediate signal transduction from G protein-coupled receptors to downstream effectors. Upon interaction with an activated receptor, G protein exchanges its bound GDP with GTP, stimulating downstream signal transmission. Release of GDP requires a structural rearrangement between the GTPase domain and helical domain of the Gα subunit. Here, we used single molecule fluorescence resonance energy transfer (smFRET) technique to study the conformational dynamics of these two domains in the apo state and in the binding of different ligands. Direct imaging of individual molecules showed that the G_iα subunit is highly dynamic, and at least three major conformations of G_iα could be observed in the apo state. Upon binding of GDP, G_iα becomes dramatically less dynamic, resulting in a closed conformation between the two domains. We postulate that changes between the three conformations are sequential, and the three conformations appear to have distinct affinities toward GDP.

Molecular mechanism of Gαi activation by non-GPCR proteins with a Gα-Binding and Activating motif

Heterotrimeric G proteins are quintessential signalling switches activated by nucleotide exchange on Gα. Although activation is predominantly carried out by G-protein-coupled receptors (GPCRs), non-receptor guanine-nucleotide exchange factors (GEFs) have emerged as critical signalling molecules and therapeutic targets. Here we characterize the molecular mechanism of G-protein activation by a family of non-receptor GEFs containing a Gα-binding and -activating (GBA) motif. We combine NMR spectroscopy, computational modelling and biochemistry to map changes in Gα caused by binding of GBA proteins with residue-level resolution. We find that the GBA motif binds to the SwitchII/α3 cleft of Gα and induces changes in the G-1/P-loop and G-2 boxes (involved in phosphate binding), but not in the G-4/G-5 boxes (guanine binding). Our findings reveal that G-protein-binding and activation mechanisms are fundamentally different between GBA proteins and GPCRs, and that GEF-mediated perturbation of nucleotide phosphate binding is sufficient for Gα activation.

Nonreceptor guanine-nucleotide exchange factors (GEFs) are emerging as important regulators of heterotrimeric G proteins. Here, the authors present structural and mechanistic insights into how a class of nonreceptor GEFs containing the Ga-Binding and Activating motif interact and modulate G proteins.

A new approach to distance measurements between two spin labels in the >10 nm range

ESR spectroscopy can be efficiently used to acquire the distance between two spin labels placed on a macromolecule by measuring their mutual dipolar interaction frequency, as long as the distance is not greater than ∼10 nm. Any hope to significantly increase this figure is hampered by the fact that all available spin labels have a phase memory time (T_m), restricted to the microseconds range, which provides a limited window during which the dipolar interaction frequency can be measured. Thus, due to the inverse cubic dependence of the dipolar frequency over the labels' separation distance, evaluating much larger distances, e.g. 20 nm, would require to have a T_m that is ∼200 microsecond, clearly beyond any hope. Here we propose a new approach to greatly enhancing the maximum measured distance available by relying on another type of dipole interaction-mediated mechanism called spin diffusion. This mechanism operates and can be evaluated during the spin lattice relaxation time, T₁ (commonly in the milliseconds range), rather than only during T_m. Up until recently, the observation of spin diffusion in solid electron spin systems was considered experimentally impractical. However, recent developments have enabled its direct measurement by means of high sensitivity pulsed ESR that employs intense short magnetic field gradients, thus opening the door to the subsequent utilization of these capabilities. The manuscript presents the subject of spin diffusion, the ways it can be directly measured, and a theoretical discussion on how intramolecular spin-pair distance, even in the range of 20-30 nm, could be accurately extracted from spin diffusion measurements.

Recent Progress in Understanding the Conformational Mechanism of Heterotrimeric G Protein Activation

Heterotrimeric G proteins are key intracellular coordinators that receive signals from cells through activation of cognate G protein-coupled receptors (GPCRs). The details of their atomic interactions and structural mechanisms have been described by many biochemical and biophysical studies. Specifically, a framework for understanding conformational changes in the receptor upon ligand binding and associated G protein activation was provided by description of the crystal structure of the β2-adrenoceptor-Gs complex in 2011. This review focused on recent findings in the conformational dynamics of G proteins and GPCRs during activation processes.

Interplay of G Protein-Coupled Receptors with the Membrane: Insights from Supra-Atomic Coarse Grain Molecular Dynamics Simulations

G protein-coupled receptors (GPCRs) are central to many fundamental cellular signaling pathways. They transduce signals from the outside to the inside of cells in physiological processes ranging from vision to immune response. It is extremely challenging to look at them individually using conventional experimental techniques. Recently, a pseudo atomistic molecular model has emerged as a valuable tool to access information on GPCRs, more specifically on their interactions with their environment in their native cell membrane and the consequences on their supramolecular organization. This approach uses the Martini coarse grain (CG) model to describe the receptors, lipids, and solvent in molecular dynamics (MD) simulations and in enough detail to allow conserving the chemical specificity of the different molecules. The elimination of unnecessary degrees of freedom has opened up large-scale simulations of the lipid-mediated supramolecular organization of GPCRs. Here, after introducing the Martini CGMD method, we review these studies carried out on various members of the GPCR family, including rhodopsin (visual receptor), opioid receptors, adrenergic receptors, adenosine receptors, dopamine receptor, and sphingosine 1-phosphate receptor. These studies have brought to light an interesting set of novel biophysical principles. The insights range from revealing localized and heterogeneous deformations of the membrane bilayer at the surface of the protein, specific interactions of lipid molecules with individual GPCRs, to the effect of the membrane matrix on global GPCR self-assembly. The review ends with an overview of the lessons learned from the use of the CGMD method, the biophysical-chemical findings on lipid-protein interplay.

Mechanistic insights into GPCR-G protein interactions

G protein-coupled receptors (GPCRs) respond to extracellular stimuli and interact with several intracellular binding partners to elicit cellular responses, including heterotrimeric G proteins. Recent structural and biophysical studies have highlighted the dynamic nature of GPCRs and G proteins and have identified specific conformational changes important for receptor-mediated nucleotide exchange on Gα. While domain separation within Gα is necessary for GDP release, opening the inter-domain interface is insufficient to stimulate nucleotide exchange. Rather, an activated receptor promotes GDP release by allosterically disrupting the nucleotide-binding site via interactions with the Gα N-termini and C-termini. Highlighting the allosteric nature of GPCRs, recent studies suggest that agonist binding alone poorly stabilizes an active conformation of several receptors. Rather, full stabilization of the receptor in an active state requires formation of the agonist-receptor-G protein ternary complex. In turn, nucleotide-free Gα is able to stabilize conformational changes around the receptor's agonist-binding site to enhance agonist affinity.

Engineering a minimal G protein to facilitate crystallisation of G protein-coupled receptors in their active conformation

G protein-coupled receptors (GPCRs) modulate cytoplasmic signalling in response to extracellular stimuli, and are important therapeutic targets in a wide range of diseases. Structure determination of GPCRs in all activation states is important to elucidate the precise mechanism of signal transduction and to facilitate optimal drug design. However, due to their inherent instability, crystallisation of GPCRs in complex with cytoplasmic signalling proteins, such as heterotrimeric G proteins and β-arrestins, has proved challenging. Here, we describe the design of a minimal G protein, mini-G_s, which is composed solely of the GTPase domain from the adenylate cyclase stimulating G protein G_s. Mini-G_s is a small, soluble protein, which efficiently couples GPCRs in the absence of Gβγ subunits. We engineered mini-G_s, using rational design mutagenesis, to form a stable complex with detergent-solubilised β₁-adrenergic receptor (β₁AR). Mini G proteins induce similar pharmacological and structural changes in GPCRs as heterotrimeric G proteins, but eliminate many of the problems associated with crystallisation of these complexes, specifically their large size, conformational dynamics and instability in detergent. They are therefore novel tools, which will facilitate the biochemical and structural characterisation of GPCRs in their active conformation.

G Protein-selective GPCR Conformations Measured Using FRET Sensors in a Live Cell Suspension Fluorometer Assay

Fӧrster resonance energy transfer (FRET)-based studies have become increasingly common in the investigation of GPCR signaling. Our research group developed an intra-molecular FRET sensor to detect the interaction between Gα subunits and GPCRs in live cells following agonist stimulation. Here, we detail the protocol for detecting changes in FRET between the β2-adrenergic receptor and the Gαs C-terminus peptide upon treatment with 100 µM isoproterenol hydrochloride as previously characterized(1). Our FRET sensor is a single polypeptide consisting serially of a full-length GPCR, a FRET acceptor fluorophore (mCitrine), an ER/K SPASM (systematic protein affinity strength modulation) linker, a FRET donor fluorophore (mCerulean), and a Gα C-terminal peptide. This protocol will detail cell preparation, transfection conditions, equipment setup, assay execution, and data analysis. This experimental design detects small changes in FRET indicative of protein-protein interactions, and can also be used to compare the strength of interaction across ligands and GPCR-G protein pairings. To enhance the signal-to-noise in our measurements, this protocol requires heightened precision in all steps, and is presented here to enable reproducible execution.

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.

Invited review: Activation of G proteins by GTP and the mechanism of Gα-catalyzed GTP hydrolysis

This review addresses the regulatory consequences of the binding of GTP to the alpha subunits (Gα) of heterotrimeric G proteins, the reaction mechanism of GTP hydrolysis catalyzed by Gα and the means by which GTPase activating proteins (GAPs) stimulate the GTPase activity of Gα. The high energy of GTP binding is used to restrain and stabilize the conformation of the Gα switch segments, particularly switch II, to afford stable complementary to the surfaces of Gα effectors, while excluding interaction with Gβγ, the regulatory binding partner of GDP-bound Gα. Upon GTP hydrolysis, the energy of these conformational restraints is dissipated and the two switch segments, particularly switch II, become flexible and are able to adopt a conformation suitable for tight binding to Gβγ. Catalytic site pre-organization presents a significant activation energy barrier to Gα GTPase activity. The glutamine residue near the N-terminus of switch II (Glncat ) must adopt a conformation in which it orients and stabilizes the γ phosphate and the water nucleophile for an in-line attack. The transition state is probably loose with dissociative character; phosphoryl transfer may be concerted. The catalytic arginine in switch I (Argcat ), together with amide hydrogen bonds from the phosphate binding loop, stabilize charge at the β-γ bridge oxygen of the leaving group. GAPs that harbor "regulator of protein signaling" (RGS) domains, or structurally unrelated domains within G protein effectors that function as GAPs, accelerate catalysis by stabilizing the pre-transition state for Gα-catalyzed GTP hydrolysis, primarily by restraining Argcat and Glncat to their catalytic conformations. © 2016 Wiley Periodicals, Inc. Biopolymers 105: 449-462, 2016.

A Conserved Hydrophobic Core in Gαi1 Regulates G Protein Activation and Release from Activated Receptor

G protein-coupled receptor-mediated heterotrimeric G protein activation is a major mode of signal transduction in the cell. Previously, we and other groups reported that the α5 helix of Gαi1, especially the hydrophobic interactions in this region, plays a key role during nucleotide release and G protein activation. To further investigate the effect of this hydrophobic core, we disrupted it in Gαi1 by inserting 4 alanine amino acids into the α5 helix between residues Gln(333) and Phe(334) (Ins4A). This extends the length of the α5 helix without disturbing the β6-α5 loop interactions. This mutant has high basal nucleotide exchange activity yet no receptor-mediated activation of nucleotide exchange. By using structural approaches, we show that this mutant loses critical hydrophobic interactions, leading to significant rearrangements of side chain residues His(57), Phe(189), Phe(191), and Phe(336); it also disturbs the rotation of the α5 helix and the π-π interaction between His(57) and Phe(189) In addition, the insertion mutant abolishes G protein release from the activated receptor after nucleotide binding. Our biochemical and computational data indicate that the interactions between α5, α1, and β2-β3 are not only vital for GDP release during G protein activation, but they are also necessary for proper GTP binding (or GDP rebinding). Thus, our studies suggest that this hydrophobic interface is critical for accurate rearrangement of the α5 helix for G protein release from the receptor after GTP binding.

Conformational dynamics of a G-protein α subunit is tightly regulated by nucleotide binding

Heterotrimeric G proteins play a pivotal role in the signal-transduction pathways initiated by G-protein-coupled receptor (GPCR) activation. Agonist-receptor binding causes GDP-to-GTP exchange and dissociation of the Gα subunit from the heterotrimeric G protein, leading to downstream signaling. Here, we studied the internal mobility of a G-protein α subunit in its apo and nucleotide-bound forms and characterized their dynamical features at multiple time scales using solution NMR, small-angle X-ray scattering, and molecular dynamics simulations. We find that binding of GTP analogs leads to a rigid and closed arrangement of the Gα subdomain, whereas the apo and GDP-bound forms are considerably more open and dynamic. Furthermore, we were able to detect two conformational states of the Gα Ras domain in slow exchange whose populations are regulated by binding to nucleotides and a GPCR. One of these conformational states, the open state, binds to the GPCR; the second conformation, the closed state, shows no interaction with the receptor. Binding to the GPCR stabilizes the open state. This study provides an in-depth analysis of the conformational landscape and the switching function of a G-protein α subunit and the influence of a GPCR in that landscape.

A GTPase chimera illustrates an uncoupled nucleotide affinity and release rate, providing insight into the activation mechanism

The release of GDP from GTPases signals the initiation of a GTPase cycle, where the association of GTP triggers conformational changes promoting binding of downstream effector molecules. Studies have implicated the nucleotide-binding G5 loop to be involved in the GDP release mechanism. For example, biophysical studies on both the eukaryotic Gα proteins and the GTPase domain (NFeoB) of prokaryotic FeoB proteins have revealed conformational changes in the G5 loop that accompany nucleotide binding and release. However, it is unclear whether this conformational change in the G5 loop is a prerequisite for GDP release, or, alternatively, the movement is a consequence of release. To gain additional insight into the sequence of events leading to GDP release, we have created a chimeric protein comprised of Escherichia coli NFeoB and the G5 loop from the human Giα1 protein. The protein chimera retains GTPase activity at a similar level to wild-type NFeoB, and structural analyses of the nucleotide-free and GDP-bound proteins show that the G5 loop adopts conformations analogous to that of the human nucleotide-bound Giα1 protein in both states. Interestingly, isothermal titration calorimetry and stopped-flow kinetic analyses reveal uncoupled nucleotide affinity and release rates, supporting a model where G5 loop movement promotes nucleotide release.

Hitchhiking on the heptahelical highway: structure and function of 7TM receptor complexes

A revolution in the analysis of seven transmembrane domain (7TM) receptors has provided detailed information about how these physiologically important signalling proteins interact with extracellular cues. However, it has proved much more challenging to understand how 7TM receptors convey information to their principal intracellular targets: heterotrimeric G proteins, G protein-coupled receptor kinases and arrestins. Recent structures now suggest a common mechanism that enables these structurally diverse cytoplasmic proteins to 'hitch a ride' on hundreds of different activated 7TM receptors in order to instigate physiological change.

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.

Characterization of the Domain Orientations of E. coli 5'-Nucleotidase by Fitting an Ensemble of Conformers to DEER Distance Distributions

Escherichia coli 5'-nucleotidase is a two-domain enzyme exhibiting a unique 96° domain motion that is required for catalysis. Here we present an integrated structural biology study that combines DEER distance distributions with structural information from X-ray crystallography and computational biology to describe the population of presumably almost isoenergetic open and closed states in solution. Ensembles of models that best represent the experimental distance distributions are determined by a Monte Carlo search algorithm. As a result, predominantly open conformations are observed in the unliganded state indicating that the majority of enzyme molecules await substrate binding for the catalytic cycle. The addition of a substrate analog yields ensembles with an almost equal mixture of open and closed states. Thus, in the presence of substrate, efficient catalysis is provided by the simultaneous appearance of open conformers (binding substrate or releasing product) and closed conformers (enabling the turnover of the substrate).

The experimental power of FR900359 to study Gq-regulated biological processes

Despite the discovery of heterotrimeric αβγ G proteins ∼25 years ago, their selective perturbation by cell-permeable inhibitors remains a fundamental challenge. Here we report that the plant-derived depsipeptide FR900359 (FR) is ideally suited to this task. Using a multifaceted approach we systematically characterize FR as a selective inhibitor of Gq/11/14 over all other mammalian Gα isoforms and elaborate its molecular mechanism of action. We also use FR to investigate whether inhibition of Gq proteins is an effective post-receptor strategy to target oncogenic signalling, using melanoma as a model system. FR suppresses many of the hallmark features that are central to the malignancy of melanoma cells, thereby providing new opportunities for therapeutic intervention. Just as pertussis toxin is used extensively to probe and inhibit the signalling of Gi/o proteins, we anticipate that FR will at least be its equivalent for investigating the biological relevance of Gq.

Role of Structural Dynamics at the Receptor G Protein Interface for Signal Transduction

GPCRs catalyze GDP/GTP exchange in the α-subunit of heterotrimeric G proteins (Gαßγ) through displacement of the Gα C-terminal α5 helix, which directly connects the interface of the active receptor (R*) to the nucleotide binding pocket of G. Hydrogen-deuterium exchange mass spectrometry and kinetic analysis of R* catalysed G protein activation have suggested that displacement of α5 starts from an intermediate GDP bound complex (R*•GGDP). To elucidate the structural basis of receptor-catalysed displacement of α5, we modelled the structure of R*•GGDP. A flexible docking protocol yielded an intermediate R*•GGDP complex, with a similar overall arrangement as in the X-ray structure of the nucleotide free complex (R*•Gempty), however with the α5 C-terminus (GαCT) forming different polar contacts with R*. Starting molecular dynamics simulations of GαCT bound to R* in the intermediate position, we observe a screw-like motion, which restores the specific interactions of α5 with R* in R*•Gempty. The observed rotation of α5 by 60° is in line with experimental data. Reformation of hydrogen bonds, water expulsion and formation of hydrophobic interactions are driving forces of the α5 displacement. We conclude that the identified interactions between R* and G protein define a structural framework in which the α5 displacement promotes direct transmission of the signal from R* to the GDP binding pocket.

BCL::MP-fold: Membrane protein structure prediction guided by EPR restraints

For many membrane proteins, the determination of their topology remains a challenge for methods like X-ray crystallography and nuclear magnetic resonance (NMR) spectroscopy. Electron paramagnetic resonance (EPR) spectroscopy has evolved as an alternative technique to study structure and dynamics of membrane proteins. The present study demonstrates the feasibility of membrane protein topology determination using limited EPR distance and accessibility measurements. The BCL::MP-Fold (BioChemical Library membrane protein fold) algorithm assembles secondary structure elements (SSEs) in the membrane using a Monte Carlo Metropolis (MCM) approach. Sampled models are evaluated using knowledge-based potential functions and agreement with the EPR data and a knowledge-based energy function. Twenty-nine membrane proteins of up to 696 residues are used to test the algorithm. The RMSD100 value of the most accurate model is better than 8 Å for 27, better than 6 Å for 22, and better than 4 Å for 15 of the 29 proteins, demonstrating the algorithms' ability to sample the native topology. The average enrichment could be improved from 1.3 to 2.5, showing the improved discrimination power by using EPR data.

The High-Resolution Structure of Activated Opsin Reveals a Conserved Solvent Network in the Transmembrane Region Essential for Activation

Rhodopsin, a light-activated G protein coupled receptor (GPCR), has been the subject of numerous biochemical and structural investigations, serving as a model receptor for GPCRs and their activation. We present the 2.3-Å resolution structure of native source rhodopsin stabilized in a conformation competent for G protein binding. An extensive water-mediated hydrogen bond network linking the chromophore binding site to the site of G protein binding is observed, providing connections to conserved motifs essential for GPCR activation. Comparison of this extensive solvent-mediated hydrogen-bonding network with the positions of ordered solvent in earlier crystallographic structures of rhodopsin photointermediates reveals both static structural and dynamic functional water-protein interactions present during the activation process. When considered along with observations that solvent occupies similar positions in the structures of other GPCRs, these analyses strongly support an integral role for this dynamic ordered water network in both rhodopsin and GPCR activation.

Probing Gαi1 protein activation at single-amino acid resolution

We present comprehensive maps at single-amino acid resolution of the residues stabilizing the human Gαi1 subunit in nucleotide- and receptor-bound states. We generated these maps by measuring the effects of alanine mutations on the stability of Gαi1 and the rhodopsin-Gαi1 complex. We identified stabilization clusters in the GTPase and helical domains responsible for structural integrity and the conformational changes associated with activation. In activation cluster I, helices α1 and α5 pack against strands β1-β3 to stabilize the nucleotide-bound states. In the receptor-bound state, these interactions are replaced by interactions between α5 and strands β4-β6. Key residues in this cluster are Y320, which is crucial for the stabilization of the receptor-bound state, and F336, which stabilizes nucleotide-bound states. Destabilization of helix α1, caused by rearrangement of this activation cluster, leads to the weakening of the interdomain interface and release of GDP.

Universal allosteric mechanism for Gα activation by GPCRs

G protein-coupled receptors (GPCRs) allosterically activate heterotrimeric G proteins and trigger GDP release. Given that there are ∼800 human GPCRs and 16 different Gα genes, this raises the question of whether a universal allosteric mechanism governs Gα activation. Here we show that different GPCRs interact with and activate Gα proteins through a highly conserved mechanism. Comparison of Gα with the small G protein Ras reveals how the evolution of short segments that undergo disorder-to-order transitions can decouple regions important for allosteric activation from receptor binding specificity. This might explain how the GPCR-Gα system diversified rapidly, while conserving the allosteric activation mechanism.

SIGNAL TRANSDUCTION. Structural basis for nucleotide exchange in heterotrimeric G proteins

G protein-coupled receptors (GPCRs) relay diverse extracellular signals into cells by catalyzing nucleotide release from heterotrimeric G proteins, but the mechanism underlying this quintessential molecular signaling event has remained unclear. Here we use atomic-level simulations to elucidate the nucleotide-release mechanism. We find that the G protein α subunit Ras and helical domains-previously observed to separate widely upon receptor binding to expose the nucleotide-binding site-separate spontaneously and frequently even in the absence of a receptor. Domain separation is necessary but not sufficient for rapid nucleotide release. Rather, receptors catalyze nucleotide release by favoring an internal structural rearrangement of the Ras domain that weakens its nucleotide affinity. We use double electron-electron resonance spectroscopy and protein engineering to confirm predictions of our computationally determined mechanism.

Nullspace Sampling with Holonomic Constraints Reveals Molecular Mechanisms of Protein Gαs

Proteins perform their function or interact with partners by exchanging between conformational substates on a wide range of spatiotemporal scales. Structurally characterizing these exchanges is challenging, both experimentally and computationally. Large, diffusional motions are often on timescales that are difficult to access with molecular dynamics simulations, especially for large proteins and their complexes. The low frequency modes of normal mode analysis (NMA) report on molecular fluctuations associated with biological activity. However, NMA is limited to a second order expansion about a minimum of the potential energy function, which limits opportunities to observe diffusional motions. By contrast, kino-geometric conformational sampling (KGS) permits large perturbations while maintaining the exact geometry of explicit conformational constraints, such as hydrogen bonds. Here, we extend KGS and show that a conformational ensemble of the α subunit Gαs of heterotrimeric stimulatory protein Gs exhibits structural features implicated in its activation pathway. Activation of protein Gs by G protein-coupled receptors (GPCRs) is associated with GDP release and large conformational changes of its α-helical domain. Our method reveals a coupled α-helical domain opening motion while, simultaneously, Gαs helix α5 samples an activated conformation. These motions are moderated in the activated state. The motion centers on a dynamic hub near the nucleotide-binding site of Gαs, and radiates to helix α4. We find that comparative NMA-based ensembles underestimate the amplitudes of the motion. Additionally, the ensembles fall short in predicting the accepted direction of the full activation pathway. Taken together, our findings suggest that nullspace sampling with explicit, holonomic constraints yields ensembles that illuminate molecular mechanisms involved in GDP release and protein Gs activation, and further establish conformational coupling between key structural elements of Gαs.