Two distinct conformations of helix 6 observed in antagonist-bound structures of a beta1-adrenergic receptor

A mutagenesis and screening strategy to generate optimally thermostabilized membrane proteins for structural studies

The thermostability of an integral membrane protein (MP) in detergent solution is a key parameter that dictates the likelihood of obtaining well-diffracting crystals that are suitable for structure determination. However, many mammalian MPs are too unstable for crystallization. We developed a thermostabilization strategy based on systematic mutagenesis coupled to a radioligand-binding thermostability assay that can be applied to receptors, ion channels and transporters. It takes ∼6-12 months to thermostabilize a G-protein-coupled receptor (GPCR) containing 300 amino acid (aa) residues. The resulting thermostabilized MPs are more easily crystallized and result in high-quality structures. This methodology has facilitated structure-based drug design applied to GPCRs because it is possible to determine multiple structures of the thermostabilized receptors bound to low-affinity ligands. Protocols and advice are given on how to develop thermostability assays for MPs and how to combine mutations to make an optimally stable mutant suitable for structural studies. The steps in the procedure include the generation of ∼300 site-directed mutants by Ala/Leu scanning mutagenesis, the expression of each mutant in mammalian cells by transient transfection and the identification of thermostable mutants using a thermostability assay that is based on binding of an (125)I-labeled radioligand to the unpurified, detergent-solubilized MP. Individual thermostabilizing point mutations are then combined to make an optimally stable MP that is suitable for structural biology and other biophysical studies.

Allosteric nanobodies reveal the dynamic range and diverse mechanisms of G-protein-coupled receptor activation

G-protein coupled receptors (GPCRs) modulate many physiological processes by transducing a variety of extracellular cues into intracellular responses. Ligand binding to an extracellular orthosteric pocket propagates conformational change to the receptor cytosolic region to promote binding and activation of downstream signaling effectors such as G proteins and β-arrestins. It is widely appreciated that different agonists can share the same binding pocket but evoke unique receptor conformations leading to a wide range of downstream responses (i.e., ‘efficacy’)¹. Furthermore, mounting biophysical evidence, primarily using the β-adrenergic receptor (β₂AR) as a model system, supports the existence of multiple active and inactive conformational states^2–5. However, how agonists with varying efficacy modulate these receptor states to initiate cellular responses is not well understood. Here we report stabilization of two distinct β₂AR conformations using single domain camelid antibodies (nanobodies): a previously described positive allosteric nanobody (Nb80) and a newly identified negative allosteric nanobody (Nb60)^6,7. We show that Nb60 stabilizes a previously unappreciated low affinity receptor state which corresponds to one of two inactive receptor conformations as delineated by X-ray crystallography and NMR spectroscopy. We find that the agonist isoproterenol has a 15,000-fold higher affinity for the β₂AR in the presence of Nb80 compared to Nb60, highlighting the full allosteric range of a GPCR. Assessing the binding of 17 ligands of varying efficacy to the β₂AR in the absence and presence of Nb60 or Nb80 reveals large ligand-specific effects that can only be explained using an allosteric model which assumes equilibrium amongst at least three receptor states. Agonists generally exert efficacy by stabilizing the active Nb80-stabilized receptor state (R₈₀). In contrast, for a number of partial agonists, both stabilization of R₈₀ and destabilization of the inactive, Nb60-bound state (R₆₀) contribute to their ability to modulate receptor activation. These data demonstrate that ligands can initiate a wide range of cellular responses by differentially stabilizing multiple receptor states.

In silico Exploration of the Conformational Universe of GPCRs

The structural plasticity of G protein coupled receptors (GPCRs) leads to a conformational universe going from inactive to active receptor states with several intermediate states. Many of them have not been captured yet and their role for GPCR activation is not well understood. The study of this conformational space and the transition dynamics between different receptor populations is a major challenge in molecular biophysics. The rational design of effector molecules that target such receptor populations allows fine-tuning receptor signalling with higher specificity to produce drugs with safer therapeutic profiles. In this minireview, we outline highly conserved receptor regions which are considered determinant for the establishment of distinct receptor states. We then discuss in-silico approaches such as dimensionality reduction methods and Markov State Models to explore the GPCR conformational universe and exploit the obtained conformations through structure-based drug design.

Application of advanced X-ray methods in life sciences

BACKGROUND

Synchrotron radiation (SR) sources provide diverse X-ray methods for the investigation of structure-function relationships in biological macromolecules.

SCOPE OF REVIEW

Recent developments in SR sources and in the X-ray tools they offer for life sciences are reviewed. Specifically, advances in macromolecular crystallography, small angle X-ray solution scattering, X-ray absorption and fluorescence spectroscopy, and imaging are discussed with examples.

MAJOR CONCLUSIONS

SR sources offer a range of X-ray techniques that can be used in a complementary fashion in studies of biological systems at a wide range of resolutions from atomic to cellular scale. Emerging applications of X-ray techniques include the characterization of disordered proteins, noncrystalline and nonequilibrium systems, elemental imaging of tissues, cells and organs, and detection of time-resolved changes in molecular structures.

GENERAL SIGNIFICANCE

X-ray techniques are in the center of hybrid approaches that are used to gain insight into complex problems relating to biomolecular mechanisms, disease and possible therapeutic solutions. This article is part of a Special Issue entitled "Science for Life". Guest Editors: Dr. Austen Angell, Dr. Salvatore Magazù and Dr. Federica Migliardo.

In vitro, molecular modeling and behavioral studies of 3-{[4-(5-methoxy-1H-indol-3-yl)-1,2,3,6-tetrahydropyridin-1-yl]methyl}-1,2-dihydroquinolin-2-one (D2AAK1) as a potential antipsychotic

Antipsychotics currently available to treat schizophrenia suffer several limitations: (1) they are efficient against positive but not negative and cognitive symptoms of the disease; (2) they help only a half of patients; (3) they have severe side effects including neurological and metabolic side effects. Thus, novel drugs to treat schizophrenia are highly demanded. We identified a novel dopamine D2 receptor antagonist, D2AAK1, with Ki of 58 nM using structure-based virtual screening. D2AAK1 possesses additional nanomolar or low micromolar affinity to D1, D3, 5-HT1A and 5-HT2A receptors, making it an ideal candidate for a multi-target drug. Here we present homology modeling, molecular docking and molecular dynamics of D2AAK1 and its molecular targets and animal studies of D2AAK1 as a potential antipsychotic. The main contact of D2AAK1 and all the receptors studied is the electrostatic interaction between the protonable nitrogen atom of the ligand and the conserved Asp(3.32) as typical for orthosteric ligands of aminergic GPCRs. We confirmed antagonistic/partial agonistic properties of D2AAK1 towards all the receptors in in vitro essays and in in silico studies as the ligand stabilizes the ionic lock interaction. We also demonstrated neuroleptic, anxiolytic and, importantly, procognitive properties of D2AAK1 in mouse models.

Comparative MD Simulations Indicate a Dual Role for Arg1323.50 in Dopamine-Dependent D2R Activation

Residue Arg^3.50 belongs to the highly conserved DRY-motif of class A GPCRs, which is located at the bottom of TM3. On the one hand, Arg^3.50 has been reported to help stabilize the inactive state of GPCRs, but on the other hand has also been shown to be crucial for stabilizing active receptor conformations and mediating receptor-G protein coupling. The combined results of these studies suggest that the exact function of Arg^3.50 is likely to be receptor-dependent and must be characterized independently for every GPCR. Consequently, we now present comparative molecular-dynamics simulations that use our recently described inactive-state and Gα-bound active-state homology models of the dopamine D2 receptor (D2R), which are either bound to dopamine or ligand-free, performed to identify the function of Arg132^3.50 in D2R. Our results are consistent with a dynamic model of D2R activation in which Arg132^3.50 adopts a dual role, both by stabilizing the inactive-state receptor conformation and enhancing dopamine-dependent D2R-G protein coupling.

Identification of biosynthetic intermediates of teaghrelins and teaghrelin-like compounds in oolong teas, and their molecular docking to the ghrelin receptor

Teaghrelins are unique acylated flavonoid tetraglycosides found in Chin-shin oolong tea, and have been demonstrated to be promising oral ghrelin analogues. The biosynthetic pathway of teaghrelins from quercetin-3-O-rutinoside (rutin) or kaempferol-3-O-rutinoside (nicotiflorin) was proposed to comprise three enzymatic steps according to the identification of putative intermediates in Chin-shin oolong tea. In addition to the two known teaghrelins in Chin-shin oolong tea, four teaghrelin-like compounds with different attachments of glycosides were identified in various oolong teas. Molecular modeling and docking were used to evaluate theoretically whether the putative biosynthetic intermediates of teaghrelins and the four teaghrelin-like compounds could be potential candidates of ghrelin analogues. The results showed that the attachment of a coumaroyl group was crucial for these tea compounds to bind to the ghrelin receptor. However, the additional attachment of a rhamnosyl glycoside to the flavonoid backbone of teaghrelin-like compounds at C-7 significantly reduced their binding affinity with the ghrelin receptor.

How Can Mutations Thermostabilize G-Protein-Coupled Receptors?

Structures of over 30 different G-protein-coupled receptors (GPCRs) have advanced our understanding of cell signaling and have provided a foundation for structure-guided drug design. This exciting progress has required the development of three complementary methods to facilitate GPCR crystallization, one of which is the thermostabilization of receptors by systematic mutagenesis. However, the reason why a particular mutation, or combination of mutations, stabilizes the receptor is not always evident from a static crystal structure. Molecular dynamics (MD) simulations have been used to identify and estimate the energetic factors that affect thermostability through comparing the dynamics of the thermostabilized receptors with structure-based models of the wild-type receptor. The data indicate that receptors are stabilized through a combination of factors, including an increase in receptor rigidity, a decrease in collective motion, reduced stress at specific residues, and the presence of ordered water molecules. Predicting thermostabilizing mutations computationally represents a major challenge for the field.

Towards predictive docking at aminergic G-protein coupled receptors

G protein-coupled receptors (GPCRs) are hard to crystallize. However, attempts to predict their structure have boomed as a result of advancements in crystallographic techniques. This trend has allowed computer-aided molecular modeling of GPCRs. We analyzed the performance of four molecular modeling programs in pose evaluation of re-docked antagonists / inverse agonists to 11 original crystal structures of aminergic GPCRs using an induced fit-docking procedure. AutoDock and Glide were used for docking. AutoDock binding energy function, GlideXP, Prime MM-GB/SA, and YASARA binding function were used for pose scoring. Root mean square deviation (RMSD) of the best pose ranged from 0.09 to 1.58 Å, and median RMSD of the top 60 poses ranged from 1.47 to 3.83 Å. However, RMSD of the top pose ranged from 0.13 to 7.33 Å and ranking of the best pose ranged from the 1st to 60th out of 60 poses. Moreover, analysis of ligand-receptor interactions of top poses revealed substantial differences from interactions found in crystallographic structures. Bad ranking of top poses and discrepancies between top docked poses and crystal structures render current simple docking methods unsuitable for predictive modeling of receptor-ligand interactions. Prime MM-GB/SA optimized for 3NY9 by multiple linear regression did not work well at 3NY8 and 3NYA, structures of the same receptor with different ligands. However, 9 of 11 trajectories of molecular dynamics simulations by Desmond of top poses converged with trajectories of crystal structures. Key interactions were properly detected for all structures. This procedure also worked well for cross-docking of tested β2-adrenergic antagonists. Thus, this procedure represents a possible way to predict interactions of antagonists with aminergic GPCRs.

Re-Designing of Existing Pharmaceuticals for Environmental Biodegradability: A Tiered Approach with β-Blocker Propranolol as an Example

Worldwide, contamination of aquatic systems with micropollutants, including pharmaceuticals, is one of the challenges for sustainable management of water resources. Although micropollutants are present at low concentrations, many of them raise considerable toxicological concerns, particularly when present as components of complex mixtures. Recent research has shown that this problem cannot be sustainably solved with advanced effluent treatment. Therefore, an alternative that might overcome these environmental problems is the design of new pharmaceutical molecules or the redesign of existing pharmaceutical molecules that present the functionality needed for their application and have improved environmental biodegradability. Such redesigning can be performed by small molecular changes in the drug molecule with intact drug moiety which could incorporate the additional attribute such as biodegradability while retaining its pharmacological potency. This proof of concept study provides an approach for the rational redesign of a given pharmaceutical (Propranolol as an example). New derivatives with small molecular changes as compared to propranolol molecule were generated by a nontargeted photolysis process. Generated derivatives with intact drug moieties (an aromatic ring and a β-ethanolamine moiety) were further screened for aerobic biodegradability and pharmacological potency. The feasibility of the approach of redesigning an existing pharmaceutical through nontargeted generation of new derivatives with intact drug moiety and through subsequent screening was demonstrated in this study. Application of such approaches in turn might contribute to the protection of water resources in a truly sustainable manner.

Pharmacological Analysis and Structure Determination of 7-Methylcyanopindolol-Bound β1-Adrenergic Receptor

Comparisons between structures of the β1-adrenergic receptor (AR) bound to either agonists, partial agonists, or weak partial agonists led to the proposal that rotamer changes of Ser(5.46), coupled to a contraction of the binding pocket, are sufficient to increase the probability of receptor activation. (RS)-4-[3-(tert-butylamino)-2-hydroxypropoxy]-1H-indole-2-carbonitrile (cyanopindolol) is a weak partial agonist of β1AR and, based on the hypothesis above, we predicted that the addition of a methyl group to form 4-[(2S)-3-(tert-butylamino)-2-hydroxypropoxy]-7-methyl-1H-indole-2-carbonitrile (7-methylcyanopindolol) would dramatically reduce its efficacy. An eight-step synthesis of 7-methylcyanopindolol was developed and its pharmacology was analyzed. 7-Methylcyanopindolol bound with similar affinity to cyanopindolol to both β1AR and β2AR. As predicted, the efficacy of 7-methylcyanopindolol was reduced significantly compared with cyanopindolol, acting as a very weak partial agonist of turkey β1AR and an inverse agonist of human β2AR. The structure of 7-methylcyanopindolol-bound β1AR was determined to 2.4-Å resolution and found to be virtually identical to the structure of cyanopindolol-bound β1AR. The major differences in the orthosteric binding pocket are that it has expanded by 0.3 Å in 7-methylcyanopindolol-bound β1AR and the hydroxyl group of Ser(5.46) is positioned 0.8 Å further from the ligand, with respect to the position of the Ser(5.46) side chain in cyanopindolol-bound β1AR. Thus, the molecular basis for the reduction in efficacy of 7-methylcyanopindolol compared with cyanopindolol may be regarded as the opposite of the mechanism proposed for the increase in efficacy of agonists compared with antagonists.

Structural Insights into the Dynamic Process of β2-Adrenergic Receptor Signaling

G-protein-coupled receptors (GPCRs) transduce signals from the extracellular environment to intracellular proteins. To gain structural insight into the regulation of receptor cytoplasmic conformations by extracellular ligands during signaling, we examine the structural dynamics of the cytoplasmic domain of the β2-adrenergic receptor (β2AR) using (19)F-fluorine NMR and double electron-electron resonance spectroscopy. These studies show that unliganded and inverse-agonist-bound β2AR exists predominantly in two inactive conformations that exchange within hundreds of microseconds. Although agonists shift the equilibrium toward a conformation capable of engaging cytoplasmic G proteins, they do so incompletely, resulting in increased conformational heterogeneity and the coexistence of inactive, intermediate, and active states. Complete transition to the active conformation requires subsequent interaction with a G protein or an intracellular G protein mimetic. These studies demonstrate a loose allosteric coupling of the agonist-binding site and G-protein-coupling interface that may generally be responsible for the complex signaling behavior observed for many GPCRs.

Quality control in eukaryotic membrane protein overproduction

The overexpression of authentically folded eukaryotic membrane proteins in milligramme quantities is a fundamental prerequisite for structural studies. One of the most commonly used expression systems for the production of mammalian membrane proteins is the baculovirus expression system in insect cells. However, a detailed analysis by radioligand binding and comparative Western blotting of G protein-coupled receptors and a transporter produced in insect cells showed that a considerable proportion of the expressed protein was misfolded and incapable of ligand binding. In contrast, production of the same membrane proteins in stable inducible mammalian cell lines suggested that the majority was folded correctly. It was noted that detergent solubilisation of the misfolded membrane proteins using either digitonin or dodecylmaltoside was considerablyless efficient than using sodium dodecyl sulfate or foscholine-12, whilst these detergents were equally efficient at solubilising correctly folded membrane proteins. This provides a simple and rapid test to suggest whether heterologously expressed mammalian membrane proteins are indeed correctly folded, without requiring radioligand binding assays. This will greatly facilitate the high-throughput production of fully functional membrane proteins for structural studies.

3D structure prediction of human β1-adrenergic receptor via threading-based homology modeling for implications in structure-based drug designing

Dilated cardiomyopathy is a disease of left ventricular dysfunction accompanied by impairment of the β₁-adrenergic receptor (β₁-AR) signal cascade. The disturbed β₁-AR function may be based on an elevated sympathetic tone observed in patients with heart failure. Prolonged adrenergic stimulation may induce metabolic and electrophysiological disturbances in the myocardium, resulting in tachyarrhythmia that leads to the development of heart failure in human and sudden death. Hence, β₁-AR is considered as a promising drug target but attempts to develop effective and specific drug against this tempting pharmaceutical target is slowed down due to the lack of 3D structure of Homo sapiens β₁-AR (hsβADR1). This study encompasses elucidation of 3D structural and physicochemical properties of hsβADR1 via threading-based homology modeling. Furthermore, the docking performance of several docking programs including Surflex-Dock, FRED, and GOLD were validated by re-docking and cross-docking experiments. GOLD and Surflex-Dock performed best in re-docking and cross docking experiments, respectively. Consequently, Surflex-Dock was used to predict the binding modes of four hsβADR1 agonists. This study provides clear understanding of hsβADR1 structure and its binding mechanism, thus help in providing the remedial solutions of cardiovascular, effective treatment of asthma and other diseases caused by malfunctioning of the target protein.

Effect of a toggle switch mutation in TM6 of the human adenosine A₃ receptor on Gi protein-dependent signalling and Gi-independent receptor internalization

BACKGROUND AND PURPOSE

The highly conserved tryptophan (W6.48) in transmembrane domain 6 of GPCRs has been shown to play a central role in forming an active conformation in response to agonist binding. We set out to characterize the effect of this mutation on the efficacy of two agonists at multiple signalling pathways downstream of the adenosine A₃ receptor.

EXPERIMENTAL APPROACH

Residue W6.48 in the human adenosine A₃ receptor fused to yellow fluorescent protein was mutated to phenylalanine and expressed in CHO-K1 cells containing a cAMP response element reporter gene. The effects on agonist-mediated receptor internalization were monitored by automated confocal microscopy and image analysis. Further experiments were carried out to investigate agonist-mediated ERK1/2 phosphorylation, inhibition of [³H]-cAMP accumulation and β-arrestin2 binding.

KEY RESULTS

NECA was able to stimulate agonist-mediated internalization of the W6.48F mutant receptor, while the agonist HEMADO was inactive. Investigation of other downstream signalling pathways indicated that G-protein coupling was impaired for both agonists tested. Mutation of W6.48F therefore resulted in differential effects on agonist efficacy, and introduced signalling pathway bias for HEMADO at the adenosine A₃ receptor.

CONCLUSIONS AND IMPLICATIONS

Investigation of the pharmacology of the W6.48F mutant of the adenosine A₃ receptor confirms that this region is important in forming the active conformation of the receptor for stimulating a number of different signalling pathways and that mutations in this residue can lead to changes in agonist efficacy and signalling bias.

Large-scale production and protein engineering of G protein-coupled receptors for structural studies

Structural studies of G protein-coupled receptors (GPCRs) gave insights into molecular mechanisms of their action and contributed significantly to molecular pharmacology. This is primarily due to technical advances in protein engineering, production and crystallization of these important receptor targets. On the other hand, NMR spectroscopy of GPCRs, which can provide information about their dynamics, still remains challenging due to difficulties in preparation of isotopically labeled receptors and their low long-term stabilities. In this review, we discuss methods used for expression and purification of GPCRs for crystallographic and NMR studies. We also summarize protein engineering methods that played a crucial role in obtaining GPCR crystal structures.

Molecular Determinants of CGS21680 Binding to the Human Adenosine A2A Receptor

The adenosine A2A receptor (A(2A)R) plays a key role in transmembrane signaling mediated by the endogenous agonist adenosine. Here, we describe the crystal structure of human A2AR thermostabilized in an active-like conformation bound to the selective agonist 2-[p-(2-carboxyethyl)phenylethyl-amino]-5'-N-ethylcarboxamido adenosine (CGS21680) at a resolution of 2.6 Å. Comparison of A(2A)R structures bound to either CGS21680, 5'-N-ethylcarboxamido adenosine (NECA), UK432097 [6-(2,2-diphenylethylamino)-9-[(2R,3R,4S,5S)-5-(ethylcarbamoyl)-3,4-dihydroxy-tetrahydrofuran-2-yl]-N-[2-[[1-(2-pyridyl)-4-piperidyl]carbamoylamino]ethyl]purine-2-carboxamide], or adenosine shows that the adenosine moiety of the ligands binds to the receptor in an identical fashion. However, an extension in CGS21680 compared with adenosine, the (2-carboxyethyl)phenylethylamino group, binds in an extended vestibule formed from transmembrane regions 2 and 7 (TM2 and TM7) and extracellular loops 2 and 3 (EL2 and EL3). The (2-carboxyethyl)phenylethylamino group makes van der Waals contacts with side chains of amino acid residues Glu169(EL2), His264(EL3), Leu267(7.32), and Ile274(7.39), and the amine group forms a hydrogen bond with the side chain of Ser67(2.65). Of these residues, only Ile274(7.39) is absolutely conserved across the human adenosine receptor subfamily. The major difference between the structures of A(2A)R bound to either adenosine or CGS21680 is that the binding pocket narrows at the extracellular surface when CGS21680 is bound, due to an inward tilt of TM2 in that region. This conformation is stabilized by hydrogen bonds formed by the side chain of Ser67(2.65) to CGS21680, either directly or via an ordered water molecule. Mutation of amino acid residues Ser67(2.65), Glu169(EL2), and His264(EL3), and analysis of receptor activation either in the presence or absence of ligands implicates this region in modulating the level of basal activity of A(2A)R.

Characterization of G protein coupling mediated by the conserved D134(3.49) of DRY motif, M241(6.34), and F251(6.44) residues on human CXCR1

•

Mutations in the DRY motif of CXCR1 abolish ligand binding and receptor activation.

•

Point mutations between TM6 and i3 loop result in constitutive activity of CXCR1.

•

Constitutive activity of mutant CXCR1 occurs via Gα15 signaling activation.

•

The highly conserved DRY motifs have distinct roles in CXCR1 and CXCR2.

CXCR1, a receptor for interleukin-8 (IL-8), plays an important role in defending against pathogen invasion during neutrophil-mediated innate immune response. Human CXCR1 is a G protein-coupled receptor (GPCR) with its characteristic seven transmembrane domains (TMs). Functional and structural analyses of several GPCRs have revealed that conserved residues on TM3 (including the highly conserved Asp-Arg-Tyr (DRY) motif) and TM6 near intracellular loops contain domains critical for G protein coupling as well as GPCR activation. The objective of this study was to elucidate the role of critical amino acid residues on TM3 near intracellular loop 2 (i2) and TM6 near intracellular loop 3 (i3), including S132^3.47 (Baldwin location), D134^3.49, M241^6.34, and F251^6.44, in G protein coupling and CXCR1 activation. The results demonstrate that mutations of D134^3.49 at DRY motif of CXCR1 (D134N and D134V) completely abolished the ligand binding and functional response of the receptor. Additionally, point mutations at positions 241 and 251 between TM6 and i3 loop generated mutant receptors with modest constitutive activity via Gα15 signaling activation. Our results show that D134^3.49 on the highly conserved DRY motif has a distinct role for CXCR1 compared to its homologues (CXCR2 and KSHV-GPCR) in G protein coupling and receptor activation. In addition, M241^6.34 and F251^6.44 along with our previously identified V247^6.40 on TM6 are spatially located in a “hot spot” likely essential for CXCR1 activation. Identification of these amino acid residues may be useful for elucidating mechanism of CXCR1 activation and designing specific antagonists for the treatment of CXCR1-mediated diseases.

Advances in receptor conformation research: the quest for functionally selective conformations focusing on the β2-adrenoceptor

Seven-transmembrane receptors, also called GPCRs, represent the largest class of drug targets. Upon ligand binding, a GPCR undergoes conformational rearrangement and thereby changes its interaction with effector proteins including the cognate G-proteins and the multifunctional adaptor proteins, β-arrestins. These proteins, by initiating distinct signal transduction mechanisms, mediate one or several functional responses. Recently, the concept of ligand-directed GPCR signalling, also called functional selectivity or biased agonism, has been proposed to explain the phenomenon that chemically diverse ligands exhibit different efficacies towards the different signalling pathways of a single GPCR, and thereby act as functionally selective or 'biased' ligands. Current concepts support the notion that ligand-specific GPCR conformations are the basis of ligand-directed signalling. Multiple studies using fluorescence spectroscopy, X-ray crystallography, mass spectroscopy, nuclear magnetic resonance spectroscopy, single-molecule force spectroscopy and other techniques have provided the evidence to support this notion. It is anticipated that these techniques will ultimately help elucidate the structural basis of ligand-directed GPCR signalling at a precision meaningful for structure-based drug design and how a specific ligand molecular structure induces a unique receptor conformation leading to biased signalling. In this review, we will summarize recent advances in experimental techniques applied in the study of functionally selective GPCR conformations and breakthrough data obtained in these studies particularly those of the β2-adrenoceptor.

Residues remote from the binding pocket control the antagonist selectivity towards the corticotropin-releasing factor receptor-1

The corticotropin releasing factors receptor-1 and receptor-2 (CRF1R and CRF2R) are therapeutic targets for treating neurological diseases. Antagonists targeting CRF1R have been developed for the potential treatment of anxiety disorders and alcohol addiction. It has been found that antagonists targeting CRF1R always show high selectivity, although CRF1R and CRF2R share a very high rate of sequence identity. This has inspired us to study the origin of the selectivity of the antagonists. We have therefore built a homology model for CRF2R and carried out unbiased molecular dynamics and well-tempered metadynamics simulations for systems with the antagonist CP-376395 in CRF1R or CRF2R to address this issue. We found that the side chain of Tyr(6.63) forms a hydrogen bond with the residue remote from the binding pocket, which allows Tyr(6.63) to adopt different conformations in the two receptors and results in the presence or absence of a bottleneck controlling the antagonist binding to or dissociation from the receptors. The rotameric switch of the side chain of Tyr356(6.63) allows the breaking down of the bottleneck and is a perquisite for the dissociation of CP-376395 from CRF1R.

Emoghrelin, a unique emodin derivative in Heshouwu, stimulates growth hormone secretion via activation of the ghrelin receptor

ETHNOPHARMACOLOGICAL RELEVANCE

Heshouwu, the root of Polygonum multiflorum, is an anti-aging Chinese traditional medicine. Fresh (raw) Heshouwu is commonly converted to processed Heshouwu by specialized heating to alleviate its side effects of diarrhea presumably caused by anthraquinones. However, raw Heshouwu has been noted to be better than processed Heshouwu regarding anti-aging effects. The therapeutic effects of raw Heshouwu on aging-related diseases were somehow similar to the anti-aging effects of growth hormone release induced by ghrelin

MATERIALS AND METHODS

Major ingredients in the methanol extract from raw Heshouwu were separated and identified. Emodin-8-O-(6'-O-malonyl)-glucoside, a unique anthraquinone glycoside known to be completely eliminated in the conversion process of Heshouwu was isolated. This emodin derivative, tentatively named emoghrelin, was examined for its cytotoxicity and capability of stimulating growth hormone release of rat primary anterior pituitary cells via activation of the ghrelin receptor. Moreover, molecular modeling of emoghrelin docking to the ghrelin receptor was exhibited to explore the possible interaction within the binding pocket.

RESULTS

No apparent cytotoxicity was observed for emoghrelin of 10(-7)-10(-4)M. Similar to growth hormone-releasing hormone-6 (GHRP-6), a synthetic analog of ghrelin, emoghrelin was demonstrated to stimulate growth hormone secretion of rat primary anterior pituitary cells in a dose dependent manner, and the stimulation was inhibited by [d-Arg(1), d-Phe(5), d-Trp(7,9), Leu(11)]-substance P, an antagonist of the ghrelin receptor. Molecular modeling and docking showed that emoghrelin as well as GHRP-6 could fit in and adequately interact with the binding pocket of the ghrelin receptor.

CONCLUSION

The results suggest that emoghrelin is a key ingredient accounting for the anti-aging effects of Heshouwu, and possesses great potential to be a promising non-peptidyl analog of ghrelin.

Insights into a defined secondary binding region on β-adrenoceptors and putative roles in ligand binding and drug design

Putative roles of a secondary binding region shared among beta-adrenoceptors.

Structural features of β2 adrenergic receptor: crystal structures and beyond

The beta2-adrenergic receptor (β2AR) family, which is the largest family of cell surface receptors in humans. Extra attention has been focused on the human GPCRs because they have been studied as important protein targets for pharmaceutical drug development. In fact, approximately 40% of marketed drugs directly work on GPCRs. GPCRs respond to various extracellular stimuli, such as sensory signals, neurotransmitters, chemokines, and hormones, to induce structural changes at the cytoplasmic surface, activating downstream signaling pathways, primarily through interactions with heterotrimeric G proteins or through G-protein independent pathways, such as arrestin. Most GPCRs, except for rhodhopsin, which contains covalently linked 11 cis-retinal, bind to diffusible ligands, having various conformational states between inactive and active structures. The first human GPCR structure was determined using an inverse agonist bound β2AR in 2007 and since then, more than 20 distinct GPCR structures have been solved. However, most GPCR structures were solved as inactive forms, and an agonist bound fully active structure is still hard to obtain. In a structural point of view, β2AR is relatively well studied since its fully active structure as a complex with G protein as well as several inactive structures are available. The structural comparison of inactive and active states gives an important clue in understanding the activation mechanism of β2AR. In this review, structural features of inactive and active states of β2AR, the interaction of β2AR with heterotrimeric G protein, and the comparison with β1AR will be discussed.

Mammalian olfactory receptors: molecular mechanisms of odorant detection, 3D-modeling, and structure-activity relationships

This chapter describes the main characteristics of olfactory receptor (OR) genes of vertebrates, including generation of this large multigenic family and pseudogenization. OR genes are compared in relation to evolution and among species. OR gene structure and selection of a given gene for expression in an olfactory sensory neuron (OSN) are tackled. The specificities of OR proteins, their expression, and their function are presented. The expression of OR proteins in locations other than the nasal cavity is regulated by different mechanisms, and ORs display various additional functions. A conventional olfactory signal transduction cascade is observed in OSNs, but individual ORs can also mediate different signaling pathways, through the involvement of other molecular partners and depending on the odorant ligand encountered. ORs are engaged in constitutive dimers. Ligand binding induces conformational changes in the ORs that regulate their level of activity depending on odorant dose. When present, odorant binding proteins induce an allosteric modulation of OR activity. Since no 3D structure of an OR has been yet resolved, modeling has to be performed using the closest G-protein-coupled receptor 3D structures available, to facilitate virtual ligand screening using the models. The study of odorant binding modes and affinities may infer best-bet OR ligands, to be subsequently checked experimentally. The relationship between spatial and steric features of odorants and their activity in terms of perceived odor quality are also fields of research that development of computing tools may enhance.

Structural determinants for binding, activation, and functional selectivity of the angiotensin AT1 receptor

The renin-angiotensin system (RAS) plays an important role in the pathophysiology of cardiovascular disorders. Pharmacologic interventions targeting the RAS cascade have led to the discovery of renin inhibitors, angiotensin-converting enzyme inhibitors, and AT(1) receptor blockers (ARBs) to treat hypertension and some cardiovascular and renal disorders. Mutagenesis and modeling studies have revealed that differential functional outcomes are the results of multiple active states conformed by the AT(1) receptor upon interaction with angiotensin II (Ang II). The binding of agonist is dependent on both extracellular and intramembrane regions of the receptor molecule, and as a consequence occupies more extensive area of the receptor than a non-peptide antagonist. Both agonist and antagonist bind to the same intramembrane regions to interfere with each other's binding to exhibit competitive, surmountable interaction. The nature of interactions with the amino acids in the receptor is different for each of the ARBs given the small differences in the molecular structure between drugs. AT(1) receptors attain different conformation states after binding various Ang II analogues, resulting in variable responses through activation of multiple signaling pathways. These include both classical and non-classical pathways mediated through growth factor receptor transactivations, and provide cross-communication between downstream signaling molecules. The structural requirements for AT(1) receptors to activate extracellular signal-regulated kinases 1 and 2 through G proteins, or G protein-independently through β-arrestin, are different. We review the structural and functional characteristics of Ang II and its analogs and antagonists, and their interaction with amino acid residues in the AT(1) receptor.

Molecular basis for dramatic changes in cannabinoid CB1 G protein-coupled receptor activation upon single and double point mutations

There is considerable interest in determining the activation mechanism of G protein-coupled receptors (GPCRs), one of the most important types of proteins for intercellular signaling. Recently, it was demonstrated for the cannabinoid CB1 GPCR, that a single mutation T210A could make CB1 completely inactive whereas T210I makes it essentially constitutively active. To obtain an understanding of this dramatic dependence of activity on mutation, we used first-principles-based methods to predict the ensemble of low-energy seven-helix conformations for the wild-type (WT) and mutants (T210A and T210I). We find that the transmembrane (TM) helix packings depend markedly on these mutations, leading for T210A to both TM3+TM6 and TM2+TM6 salt-bridge couplings in the cytoplasmic face that explains the inactivity of this mutant. In contrast T210I has no such couplings across the receptor explaining the ease in activating this mutant. WT has just the TM3+TM6 coupling, known to be broken upon GPCR activation. To test this hypothesis on activity, we predicted double mutants that would convert the inactive mutant to normal activity and then confirmed this experimentally. This CB1 activation mechanism, or one similar to it, is expected to play a role in other constitutively active GPCRs as well.

Strategies for improved modeling of GPCR-drug complexes: blind predictions of serotonin receptors bound to ergotamine

The recent increase in the number of atomic-resolution structures of G protein-coupled receptors (GPCRs) has contributed to a deeper understanding of ligand binding to several important drug targets. However, reliable modeling of GPCR-ligand complexes for the vast majority of receptors with unknown structure remains to be one of the most challenging goals for computer-aided drug design. The GPCR Dock 2013 assessment, in which researchers were challenged to predict the crystallographic structures of serotonin 5-HT(1B) and 5-HT(2B) receptors bound to ergotamine, provided an excellent opportunity to benchmark the current state of this field. Our contributions to GPCR Dock 2013 accurately predicted the binding mode of ergotamine with RMSDs below 1.8 Å for both receptors, which included the best submissions for the 5-HT(1B) complex. Our models also had the most accurate description of the binding sites and receptor-ligand contacts. These results were obtained using a ligand-guided homology modeling approach, which combines extensive molecular docking screening with incorporation of information from multiple crystal structures and experimentally derived restraints. In this work, we retrospectively analyzed thousands of structures that were generated during the assessment to evaluate our modeling strategies. Major contributors to accuracy were found to be improved modeling of extracellular loop two in combination with the use of molecular docking to optimize the binding site for ligand recognition. Our results suggest that modeling of GPCR-drug complexes has reached a level of accuracy at which structure-based drug design could be applied to a large number of pharmaceutically relevant targets.

Synthesis and structure-activity relationship studies in serotonin 5-HT4 receptor ligands based on a benzo[de][2,6]naphthridine scaffold

A small series of serotonin 5-HT4 receptor ligands has been designed from flexible 2-methoxyquinoline compounds 7a,b by applying the conformational constraint approach. Ligands 7a,b and the corresponding conformationally constrained analogues 8a-g were synthesized and their interactions with the 5-HT4 receptor were examined by measuring both binding affinity and the ability to promote or inhibit receptor-G protein coupling. Ester derivative 7a and conformationally constrained compound 8b were demonstrated to be the most interesting compounds showing a nanomolar 5-HT4R affinity similar to that shown by reference ligands cisapride (1) and RS-23,597-190 (4). The result was rationalized by docking studies in term of high similarity in the binding modalities of flexible 7a and conformationally constrained 8b. The intrinsic efficacy of some selected ligands was determined by evaluating the receptor-G protein coupling and the results obtained demonstrated that the nature and the position of substituents play a critical role in the interaction of these ligands with their receptor.

Allosteric sodium in class A GPCR signaling

Despite their functional and structural diversity, G-protein-coupled receptors (GPCRs) share a common mechanism of signal transduction via conformational changes in the seven-transmembrane (7TM) helical domain. New major insights into this mechanism come from the recent crystallographic discoveries of a partially hydrated sodium ion that is specifically bound in the middle of the 7TM bundle of multiple class A GPCRs. This review discusses the remarkable structural conservation and distinct features of the Na(+) pocket in this most populous GPCR class, as well as the conformational collapse of the pocket upon receptor activation. New insights help to explain allosteric effects of sodium on GPCR agonist binding and activation, and sodium's role as a potential co-factor in class A GPCR function.

Compound activity prediction using models of binding pockets or ligand properties in 3D

Transient interactions of endogenous and exogenous small molecules with flexible binding sites in proteins or macromolecular assemblies play a critical role in all biological processes. Current advances in high-resolution protein structure determination, database development, and docking methodology make it possible to design three-dimensional models for prediction of such interactions with increasing accuracy and specificity. Using the data collected in the Pocketome encyclopedia, we here provide an overview of two types of the three-dimensional ligand activity models, pocketbased and ligand property-based, for two important classes of proteins, nuclear and G-protein coupled receptors. For half the targets, the pocket models discriminate actives from property matched decoys with acceptable accuracy (the area under ROC curve, AUC, exceeding 84%) and for about one fifth of the targets with high accuracy (AUC > 95%). The 3D ligand property field models performed better than 95% in half of the cases. The high performance models can already become a basis of activity predictions for new chemicals. Family-wide benchmarking of the models highlights strengths of both approaches and helps identify their inherent bottlenecks and challenges.

Modeling the interactions between MC2R and ACTH models from human

Melanocortin system is composed of four peptide hormones namely α-, β-, -γ, and adrenocorticotropic hormone (ACTH), derived from post-translational cleavage of a polypeptide precursor 'proopiomelanocortin (POMC).' Among these hormones, ACTH, a 38 amino acid residue peptide fragment is an important hormone as it is involved in steroid secretion. In addition to this, to cite a few, this hormone is also known to induce variety of other effects, such as alterations in motor/sexual behavior, improvement in memory, and anti-inflammatory effects. To date, five melanocortin receptors (MC1R-MC5R) have been characterized with tissue-specific expression patterns and different binding affinities for each of the melanocortin hormones to regulate various biological functions. In the present work, three-dimensional (3D) models of MC2R and ACTH from human have been predicted, followed by docking and molecular dynamics simulation. While the 3D model of MC2R receptor has been predicted through threading approach, structure of ACTH was built based on ab initio technique. The MC2R model was later successfully docked onto the ACTH structure. Molecular dynamics (MD) simulation for 20 ns was used to compute the binding free energy of MC2R with ACTH model under implicit solvent conditions.

The 2.1 Å resolution structure of cyanopindolol-bound β1-adrenoceptor identifies an intramembrane Na+ ion that stabilises the ligand-free receptor

The β₁-adrenoceptor (β₁AR) is a G protein-coupled receptor (GPCR) that is activated by the endogenous agonists adrenaline and noradrenaline. We have determined the structure of an ultra-thermostable β₁AR mutant bound to the weak partial agonist cyanopindolol to 2.1 Å resolution. High-quality crystals (100 μm plates) were grown in lipidic cubic phase without the assistance of a T4 lysozyme or BRIL fusion in cytoplasmic loop 3, which is commonly employed for GPCR crystallisation. An intramembrane Na⁺ ion was identified co-ordinated to Asp87^2.50, Ser128^3.39 and 3 water molecules, which is part of a more extensive network of water molecules in a cavity formed between transmembrane helices 1, 2, 3, 6 and 7. Remarkably, this water network and Na⁺ ion is highly conserved between β₁AR and the adenosine A_2A receptor (rmsd of 0.3 Å), despite an overall rmsd of 2.4 Å for all Cα atoms and only 23% amino acid identity in the transmembrane regions. The affinity of agonist binding and nanobody Nb80 binding to β₁AR is unaffected by Na⁺ ions, but the stability of the receptor is decreased by 7.5°C in the absence of Na⁺. Mutation of amino acid side chains that are involved in the co-ordination of either Na⁺ or water molecules in the network decreases the stability of β₁AR by 5–10°C. The data suggest that the intramembrane Na⁺ and associated water network stabilise the ligand-free state of β₁AR, but still permits the receptor to form the activated state which involves the collapse of the Na⁺ binding pocket on agonist binding.

The role of protein dynamics in GPCR function: insights from the β2AR and rhodopsin

G protein-coupled receptors (GPCRs) are versatile signaling proteins that mediate complex cellular responses to hormones and neurotransmitters. Recent advances in GPCR crystallography have provided inactive and active state structures for rhodopsin and the β2 adrenergic receptor (β2AR). Although these structures suggest a two-state 'on-off' mechanism of receptor activation, other biophysical studies and observed signaling versatility suggest that GPCRs are highly dynamic and exist in a multitude of functionally distinct conformations. To fully understand how GPCRs work, we must characterize these conformations and determine how ligands affect their energetics and rates of interconversion. This brief review will compare and contrast the dynamic properties of rhodopsin and β2AR that shed light on the role of structural dynamics in their distinct signaling behaviors.

Constitutive activity in gonadotropin receptors

Constitutively active mutants (CAMs) of gonadotropin receptors are, in general, rare conditions. Luteinizing hormone-choriogonadotropin receptor (LHCGR) CAMs provoke the dramatic phenotype of familial gonadotropin-independent isosexual male-limited precocious puberty, whereas in females, there is not yet any identified phenotype. Only one isolated follicle-stimulating hormone receptor (FSHR) CAM (Asp567Gly) has so far been detected in a single male patient, besides other FSHR weak CAMs linked to pregnancy-associated ovarian hyperstimulation syndrome or to impaired desensitization and internalization. Several animal models have been developed for studying enhanced gonadotropin action; in addition to unraveling valuable new information about the possible phenotypes of isolated FSHR and LHCGR CAMs in women, the information obtained from these mouse models has served multiple translational goals, including the development of new diagnostic and therapeutic targets as well as the prediction of phenotypes for mutations not yet identified in humans. Mutagenesis and computational studies have shed important information on the physiopathogenic mechanisms leading to constitutive activity of gonadotropin receptors; a common feature in these receptor CAMs is the release of stabilizing interhelical interactions between transmembrane domains (TMDs) 3 and 6 leading to an increase, with respect to the wild-type receptor, in the solvent accessibility at the cytosolic extension of TMDs 3, 5, and 6, which involves the highly conserved Glu/Asp-Arg-Tyr/Trp sequence. In this chapter, we summarize the structural features, functional consequences, and mechanisms that lead to constitutive activation of gonadotropin receptor CAMs and provide information on pharmacological approaches that might potentially modulate gonadotropin receptor CAM function.

Recent structural advances of β1 and β2 adrenoceptors yield keys for ligand recognition and drug design

Because they represent attractive drug targets, adrenoceptors have been widely studied. Recent progress in structural data of β-adrenoceptors allows us to understand and predict key interactions in ligand recognition and receptor activation. Nevertheless, an important aspect of this process has only begun to be explored: the stabilization of a conformational state of these receptors upon contact with a ligand and the capacity of a ligand to influence receptor conformation through allosteric modulation, biased signaling, and selectivity. The aim of the present Perspective is to identify the well-defined orthosteric binding site and possible allosteric sites and to analyze the importance of the ligand-receptor interaction in the stabilization of certain receptor conformations. For this purpose, we have reviewed recent advances made through the use of X-ray data from ligand-β-adrenoceptor (including ADRB1 and ADRB2) crystal structures. Most importantly, implications in the medicinal chemistry field are explored in relation to drug design.

Relevance of rhodopsin studies for GPCR activation

Rhodopsin, the dim-light photoreceptor present in the rod cells of the retina, is both a retinal-binding protein and a G protein-coupled receptor (GPCR). Due to this conjunction, it benefits from an arsenal of spectroscopy techniques that can be used for its characterization, while being a model system for the important family of Class A (also referred to as "rhodopsin-like") GPCRs. For instance, rhodopsin has been a crucial player in the field of GPCR structural biology. Until 2007, it was the only GPCR for which a high-resolution crystal structure was available, so all structure-activity analyses on GPCRs, from structure-based drug discovery to studies of structural changes upon activation, were based on rhodopsin. At present, about a third of currently available GPCR structures are still from rhodopsin. In this review, I show some examples of how these structures can still be used to gain insight into general aspects of GPCR activation. First, the analysis of the third intracellular loop in rhodopsin structures allows us to gain an understanding of the structural and dynamic properties of this region, which is absent (due to protein engineering or poor electron density) in most of the currently available GPCR structures. Second, a detailed analysis of the structure of the transmembrane domains in inactive, intermediate and active rhodopsin structures allows us to detect early conformational changes in the process of ligand-induced GPCR activation. Finally, the analysis of a conserved ligand-activated transmission switch in the transmembrane bundle of GPCRs in the context of the rhodopsin activation cycle, allows us to suggest that the structures of many of the currently available agonist-bound GPCRs may correspond to intermediate active states. While the focus in GPCR structural biology is inevitably moving away from rhodopsin, in other aspects rhodopsin is still at the forefront. For instance, the first studies of the structural basis of disease mutants in GPCRs, or the most detailed analysis of cellular GPCR signal transduction networks using a systems biology approach, have been carried out in rhodopsin. Finally, due again to its unique properties among GPCRs, rhodopsin will likely play an important role in the application of X-ray free electron laser crystallography to time-resolved structural biology in membrane proteins. Rhodopsin, thus, still remains relevant as a model system to study the molecular mechanisms of GPCR activation. This article is part of a Special Issue entitled: Retinal Proteins-You can teach an old dog new tricks.

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.