G-protein-coupled receptor structures were not built in a day

Development of enhanced conformational sampling methods to probe the activation landscape of GPCRs

G protein-coupled receptors (GPCRs) make up the largest superfamily of integral membrane proteins and play critical signal transduction roles in many physiological processes. Developments in molecular biology, biophysical, biochemical, pharmacological, and computational techniques aimed at these important therapeutic targets are beginning to provide unprecedented details on the structural as well as functional basis of their pleiotropic signaling mediated by G proteins, β arrestins, and other transducers. This pleiotropy presents a pharmacological challenge as the same ligand-receptor interaction can cause a therapeutic effect as well as an undesirable on-target side-effect through different downstream pathways. GPCRs don't function as simple binary on-off switches but as finely tuned shape-shifting machines described by conformational ensembles, where unique subsets of conformations may be responsible for specific signaling cascades. X-ray crystallography and more recently cryo-electron microscopy are providing snapshots of some of these functionally-important receptor conformations bound to ligands and/or transducers, which are being utilized by computational methods to describe the dynamic conformational energy landscape of GPCRs. In this chapter, we review the progress in computational conformational sampling methods based on molecular dynamics and discrete sampling approaches that have been successful in complementing biophysical and biochemical studies on these receptors in terms of their activation mechanisms, allosteric effects, actions of biased ligands, and effects of pathological mutations. Some of the sampled simulation time scales are beginning to approach receptor activation time scales. The list of conformational sampling methods and example uses discussed is not exhaustive but includes representative examples that have pushed the limits of classical molecular dynamics and discrete sampling methods to describe the activation energy landscape of GPCRs.

Complex Pharmacology of Free Fatty Acid Receptors

G protein-coupled receptors (GPCRs) are historically the most successful family of drug targets. In recent times it has become clear that the pharmacology of these receptors is far more complex than previously imagined. Understanding of the pharmacological regulation of GPCRs now extends beyond simple competitive agonism or antagonism by ligands interacting with the orthosteric binding site of the receptor to incorporate concepts of allosteric agonism, allosteric modulation, signaling bias, constitutive activity, and inverse agonism. Herein, we consider how evolving concepts of GPCR pharmacology have shaped understanding of the complex pharmacology of receptors that recognize and are activated by nonesterified or "free" fatty acids (FFAs). The FFA family of receptors is a recently deorphanized set of GPCRs, the members of which are now receiving substantial interest as novel targets for the treatment of metabolic and inflammatory diseases. Further understanding of the complex pharmacology of these receptors will be critical to unlocking their ultimate therapeutic potential.

Stabilization of G protein-coupled receptors by point mutations

G protein-coupled receptors (GPCRs) are flexible integral membrane proteins involved in transmembrane signaling. Their involvement in many physiological processes makes them interesting targets for drug development. Determination of the structure of these receptors will help to design more specific drugs, however, their structural characterization has so far been hampered by the low expression and their inherent instability in detergents which made protein engineering indispensable for structural and biophysical characterization. Several approaches to stabilize the receptors in a particular conformation have led to breakthroughs in GPCR structure determination. These include truncations of the flexible regions, stabilization by antibodies and nanobodies, fusion partners, high affinity and covalently bound ligands as well as conformational stabilization by mutagenesis. In this review we focus on stabilization of GPCRs by insertion of point mutations, which lead to increased conformational and thermal stability as well as improved expression levels. We summarize existing mutagenesis strategies with different coverage of GPCR sequence space and depth of information, design and transferability of mutations and the molecular basis for stabilization. We also discuss whether mutations alter the structure and pharmacological properties of GPCRs.

SuperBiHelix method for predicting the pleiotropic ensemble of G-protein-coupled receptor conformations

There is overwhelming evidence that G-protein-coupled receptors (GPCRs) exhibit several distinct low-energy conformations, each of which might favor binding to different ligands and/or lead to different downstream functions. Understanding the function of such proteins requires knowledge of the ensemble of low-energy configurations that might play a role in this pleiotropic functionality. We earlier reported the BiHelix method for efficiently sampling the (12)(7) = 35 million conformations resulting from 30° rotations about the axis (η) of all seven transmembrane helices (TMHs), showing that the experimental structure is reliably selected as the best conformation from this ensemble. However, various GPCRs differ sufficiently in the tilts of the TMHs that this method need not predict the optimum conformation starting from any other template. In this paper, we introduce the SuperBiHelix method in which the tilt angles (θ, ϕ) are optimized simultaneously with rotations (η) efficiently enough that it is practical and sufficient to sample (5 × 3 × 5)(7) = 13 trillion configurations. This method can correctly identify the optimum structure of a GPCR starting with the template from a different GPCR. We have validated this method by predicting known crystal structure conformations starting from the template of a different protein structure. We find that the SuperBiHelix conformational ensemble includes the higher energy conformations associated with the active protein in addition to those associated with the more stable inactive protein. This methodology was then applied to design and experimentally confirm structures of three mutants of the CB1 cannabinoid receptor associated with different functions.

Phasing statistics for alpha helical membrane protein structures

In this report we highlight the latest trends in phasing methods used to solve alpha helical membrane protein structures and analyze the use of heavy atom metals for the purpose of experimental phasing. Our results reveal that molecular replacement is emerging as the most successful method for phasing alpha helical membrane proteins, with the notable exception of the transporter family, where experimentally derived phase information still remains the most effective method. To facilitate selection of heavy atoms salts for experimental phasing an analysis of these was undertaken and indicates that organic mercury salts are still the most successful heavy atoms reagents. Interestingly the use of seleno-L-methionine incorporated protein has increased since earlier studies into membrane protein phasing, so too the use of SAD and MAD as techniques for phase determination. Taken together this study provides a brief snapshot of phasing methods for alpha helical membrane proteins and suggests possible routes for heavy atom selection and phasing methods based on currently available data.

Solid-state NMR analysis of membrane proteins and protein aggregates by proton detected spectroscopy

Solid-state NMR has emerged as an important tool for structural biology and chemistry, capable of solving atomic-resolution structures for proteins in membrane-bound and aggregated states. Proton detection methods have been recently realized under fast magic-angle spinning conditions, providing large sensitivity enhancements for efficient examination of uniformly labeled proteins. The first and often most challenging step of protein structure determination by NMR is the site-specific resonance assignment. Here we demonstrate resonance assignments based on high-sensitivity proton-detected three-dimensional experiments for samples of different physical states, including a fully-protonated small protein (GB1, 6 kDa), a deuterated microcrystalline protein (DsbA, 21 kDa), a membrane protein (DsbB, 20 kDa) prepared in a lipid environment, and the extended core of a fibrillar protein (α-synuclein, 14 kDa). In our implementation of these experiments, including CONH, CO(CA)NH, CANH, CA(CO)NH, CBCANH, and CBCA(CO)NH, dipolar-based polarization transfer methods have been chosen for optimal efficiency for relatively high protonation levels (full protonation or 100 % amide proton), fast magic-angle spinning conditions (40 kHz) and moderate proton decoupling power levels. Each H-N pair correlates exclusively to either intra- or inter-residue carbons, but not both, to maximize spectral resolution. Experiment time can be reduced by at least a factor of 10 by using proton detection in comparison to carbon detection. These high-sensitivity experiments are especially important for membrane proteins, which often have rather low expression yield. Proton-detection based experiments are expected to play an important role in accelerating protein structure elucidation by solid-state NMR with the improved sensitivity and resolution.

Lipid crystallization: from self-assembly to hierarchical and biological ordering

Lipid crystallization is ubiquitous in nature, observed in biological structures as well as in commercial products and applications. In a dehydrated state most of the lipids form well ordered crystals, whereas in an aqueous environment they self-assemble into various crystalline, liquid crystalline or sometimes macroscopically disordered phases. Lipid self-organization extends further to hierarchical levels including structured emulsions and nanostructured particles. Many consumer products including cosmetics, foods and medicines account for such lipid architectures. Cell membranes primarily consist of planar lipid bilayers; however sub-cellular biomembranes are more of a convoluted type. Some of the biological entities have lipids in truly crystalline form; yet liquid crystalline lipid phases are prevalent, in general. Crystallization of fats - triglyceride lipids - has been relatively well documented and reviewed more often, but this review features other areas where lipid organization is crucial and diverse. Some recent advances along with a few explicit examples of model lipid phases and biological evidences are also reported.

Structure prediction of G protein-coupled receptors and their ensemble of functionally important conformations

G protein-coupled receptors (GPCRs) are integral membrane proteins whose "pleiotropic" nature enables transmembrane (TM) signal transduction, amplification, and diversification via G protein-coupled and β arrestin-coupled pathways. GPCRs appear to enable this by being structurally flexible and by existing in different conformational states with potentially different signaling and functional consequences. We describe a method for the prediction of the three-dimensional structures of these different conformations of GPCRs starting from their amino acid sequence. It combines a unique protocol of computational methods that first predict the TM regions of these receptors and TM helix shapes based on those regions, which is followed by a locally complete sampling of TM helix packings and their scoring that results in a few (~10-20) lowest energy conformations likely to play a role in binding to different ligands and signaling events. Prediction of the structures for multiple conformations of a GPCR is starting to enable the testing of multiple hypotheses related to GPCR activation and binding to ligands with different signaling profiles.

Bihelix: Towards de novo structure prediction of an ensemble of G-protein coupled receptor conformations

G-Protein Coupled Receptors (GPCRs) play a critical role in cellular signal transduction pathways and are prominent therapeutic targets. Recently there has been major progress in obtaining experimental structures for a few GPCRs. Each GPCR, however, exhibits multiple conformations that play a role in their function and we have been developing methods aimed at predicting structures for all these conformations. Analysis of available structures shows that these conformations differ in relative helix tilts and rotations. The essential issue is, determining how to orient each of the seven helices about its axis since this determines how it interacts with the other six helices. Considering all possible helix rotations to ensure that no important packings are overlooked, and using rotation angle increments of 30° about the helical axis would still lead to 12(7) or 35 million possible conformations each with optimal residue positions. We show in this paper how to accomplish this. The fundamental idea is to optimize the interactions between each pair of contacting helices while ignoring the other 5 and then to estimate the energies of all 35 million combinations using these pair-wise interactions. This BiHelix approach dramatically reduces the effort to examine the complete set of conformations and correctly identifies the crystal packing for the experimental structures plus other near-native packings we believe may play an important role in activation. This approach also enables a detailed structural analysis of functionally distinct conformations using helix-helix interaction energy landscapes and should be useful for other helical transmembrane proteins as well.

Characterizing and predicting the functional and conformational diversity of seven-transmembrane proteins

The activation of seven-transmembrane receptors (7TMRs) allows cells to sense their environment and convert extracellular signals (like hormone binding) into intracellular signals (through G protein-coupled and/or β arrestin-coupled pathways). A single 7TMR is capable of transducing a wide spectrum of physiological responses inside a cell by coupling to these pathways. This intracellular pleiotropic action is enabled by multiple conformations exhibited by these receptors. Developments in membrane protein structure determination technologies have led to a rapid increase in crystal structures for many 7TMRs. Majority of these receptors have been crystallized in their inactive conformation and, for some, one of the many active conformations has also been crystallized. Given the topological constraints of a lipid bilayer that results in a single fold of seven almost parallel TM helices connected by mostly unstructured loops, these structures exhibit a diversity of conformations not only across the receptors but also across the different functional forms for receptors with structures for one of the functionally active conformations. Here we present a method to characterize this conformational diversity in terms of transmembrane helix topology (TMHTOP) parameters and how to use these helix orientation parameters to predict functionally-distinct multiple conformations for these receptors. The TMHTOP parameters enable a quantification of the structural changes that underlie 7TMR activation and also sheds a unique mechanistic light on the pleiotropic nature of these receptors. It provides a common language to describe the 7TMR activation mechanisms as well as differences across many receptors in terms of visually intuitive structural parameters. Protein structure prediction methods can use these parameters to describe 7TMR conformational ensembles, which coupled to experimental data can be used to develop testable hypotheses for the structural basis of 7TMR functions.

Crystallization chaperone strategies for membrane proteins

From G protein-coupled receptors to ion channels, membrane proteins represent over half of known drug targets. Yet, structure-based drug discovery is hampered by the dearth of available three-dimensional models for this large category of proteins. Other than efforts to improve membrane protein expression and stability, current strategies to improve the ability of membrane proteins to crystallize involve examining many orthologs and DNA constructs, testing the effects of different detergents for purification and crystallization, creating a lipidic environment during crystallization, and cocrystallizing with covalent or non-covalent soluble protein chaperones with an intrinsic high propensity to crystallize. In this review, we focus on this last category, highlighting successes of crystallization chaperones in membrane protein structure determination and recent developments in crystal chaperone engineering, including molecular display to enhance chaperone crystallizability, and end with a novel generic approach in development to target any membrane protein of interest.

From molecular details of the interplay between transmembrane helices of the thyrotropin receptor to general aspects of signal transduction in family a G-protein-coupled receptors (GPCRs)

Transmembrane helices (TMHs) 5 and 6 are known to be important for signal transduction by G-protein-coupled receptors (GPCRs). Our aim was to characterize the interface between TMH5 and TMH6 of the thyrotropin receptor (TSHR) to gain molecular insights into aspects of signal transduction and regulation. A proline at TMH5 position 5.50 is highly conserved in family A GPCRs and causes a twist in the helix structure. Mutation of the TSHR-specific alanine (Ala-593⁵·⁵⁰) at this position to proline resulted in a 20-fold reduction of cell surface expression. This indicates that TMH5 in the TSHR might have a conformation different from most other family A GPCRs by forming a regular α-helix. Furthermore, linking our own and previous data from directed mutagenesis with structural information led to suggestions of distinct pairs of interacting residues between TMH5 and TMH6 that are responsible for stabilizing either the basal or the active state. Our insights suggest that the inactive state conformation is constrained by a core set of polar interactions among TMHs 2, 3, 6, and 7 and in contrast that the active state conformation is stabilized mainly by non-polar interactions between TMHs 5 and 6. Our findings might be relevant for all family A GPCRs as supported by a statistical analysis of residue properties between the TMHs of a vast number of GPCR sequences.

Somatostatin receptor-3 mediated intracellular signaling and apoptosis is regulated by its cytoplasmic terminal

In the present study, we describe the role of cytoplasmic terminal (C-tail) domain in regulating coupling to adenylyl cyclase, signaling, and apoptosis in human embryonic kidney (HEK-293) cells transfected with wild type (wt)-hSSTR3 and C-tail deleted mutants. Cells transfected with wt-hSSTR3 and C-tail mutants show comparable membrane expression; however, display decreased expression in presence of agonist. wt-hSSTR3 exists as preformed homodimer at cell surface in basal conditions and decreases in response to agonist. Cells expressing C-tail mutants also show evidence of homodimerization with the same intensity as wt-hSSTR3. The agonist-dependent inhibition of cyclic adenosine monophosphate (cAMP) was lost in cells expressing C-tail mutants. Agonist treatment in cells expressing wt-hSSTR3 resulted in inhibition of cell proliferation, increased expression of PARP-1, and TUNEL positivity in proliferating cell nuclear antigen (PCNA)-positive cells. The agonist mediated increase in membrane expression of protein tyrosine phosphatase (PTP) seen with wt-hSSTR3 was diminished in C-tail mutants, which was accompanied with the loss of receptor's ability to induce apoptosis. Taken together, our data provide new insights into C-tail-dependent regulation of cell signaling and apoptosis by hSSTR3.

Crystallizing transmembrane peptides in lipidic mesophases

Structure determination of membrane proteins by crystallographic means has been facilitated by crystallization in lipidic mesophases. It has been suggested, however, that this so-called in meso method, as originally implemented, would not apply to small protein targets having </=4 transmembrane crossings. In our study, the hypothesis that the inherent flexibility of the mesophase would enable crystallogenesis of small proteins was tested using a transmembrane pentadecapeptide, linear gramicidin, which produced structure-grade crystals. This result suggests that the in meso method should be considered as a viable means for high-resolution structure determination of integral membrane peptides, many of which are predicted to be coded for in the human genome.

Expression of G-protein coupled receptors in Escherichia coli for structural studies

To elaborate a high-performance system for expression of genes of G-protein coupled receptors (GPCR), methods of direct and hybrid expression of 17 GPCR genes in Escherichia coli and selection of strains and bacteria cultivation conditions were investigated. It was established that expression of most of the target GPCR fused with the N-terminal fragment of OmpF or Mistic using media for autoinduction provides high output (up to 50 mg/liter).

Structural insights into conformational stability of wild-type and mutant beta1-adrenergic receptor

Recent experiments to derive a thermally stable mutant of turkey beta-1-adrenergic receptor (beta1AR) have shown that a combination of six single point mutations resulted in a 20 degrees C increase in thermal stability in mutant beta1AR. Here we have used the all-atom force-field energy function to calculate a stability score to detect stabilizing point mutations in G-protein coupled receptors. The calculated stability score shows good correlation with the measured thermal stability for 76 single point mutations and 22 multiple mutants in beta1AR. We have demonstrated that conformational sampling of the receptor for various mutants improve the prediction of thermal stability by 50%. Point mutations Y227A5.58, V230A5.61, and F338M7.48 in the thermally stable mutant m23-beta1AR stabilizes key microdomains of the receptor in the inactive conformation. The Y227A5.58 and V230A5.61 mutations stabilize the ionic lock between R139(3.50) on transmembrane helix3 and E285(6.30) on transmembrane helix6. The mutation F338M7.48 on TM7 alters the interaction of the conserved motif NPxxY(x)5,6F with helix8 and hence modulates the interaction of TM2-TM7-helix8 microdomain. The D186-R317 salt bridge (in extracellular loops 2 and 3) is stabilized in the cyanopindolol-bound wild-type beta1AR, whereas the salt bridge between D184-R317 is preferred in the mutant m23. We propose that this could be the surrogate to a similar salt bridge found between the extracellular loop 2 and TM7 in beta2AR reported recently. We show that the binding energy difference between the inactive and active states is less in m23 compared to the wild-type, which explains the activation of m23 at higher norepinephrine concentration compared to the wild-type. Results from this work throw light into the mechanism behind stabilizing mutations. The computational scheme proposed in this work could be used to design stabilizing mutations for other G-protein coupled receptors.

Independent phenotype of binuclear hepatocytes and cellular localization of UbD

Mice fed DDC (0.1%) for 10 weeks, and then withdrawn from the drug for 1 month, retain the ability to form Mallory-Denk bodies (MDBs) when the drug is refed for 7 days. The number of liver cells that form MDBs increased and partially replaced normal liver cells, at the end of 7 days of refeeding DDC. The MDBs that formed were associated with increased expression of UbD (also called FAT10) in the Mallory-Denk body forming cells. UbD is over expressed in 70% of human HCCs, but its cellular localization is not well established. UbD belongs to the UbL family (ubiquitin-like), and can be linked to others proteins with their 2 C-terminal glycine to lysine. By Western Blot, UbD was found to be covalently linked with proteins. We performed immunohistochemistry on tissue from mouse liver and found that UbD was located in the cytoplasm and in one or two nuclei of the same hepatocyte. However, in primary cell culture, UbD formed speckles within the cytoplasm of the liver cell. A similar pattern of cytoplasmic localization was observed in the Hepa 1-6 cell lines, which over expressed UbD fused with GFP at the C-Terminal. The localization and the control of UbD localization remain unclear. The identification of proteins that interact with UbD and the post translational modification of UbD would help to determine the regulation of this localization and function.

Beyond rhodopsin: G protein-coupled receptor structure and modeling incorporating the beta2-adrenergic and adenosine A(2A) crystal structures

For quite some time, the majority of GPCR models have been based on a single template structure: dark-adapted bovine rhodopsin. The recent solution of β2AR, β1AR and adenosine A(2A) receptor crystal structures has dramatically expanded the GPCR structural landscape and provided many new insights into receptor conformation and ligand binding. They will serve as templates for the next generation of GPCR models, but also allow direct validation of previous models and computational techniques. This review summarizes key findings from the new structures, comparison of existing models to these structures and highlights new models constructed from these templates.