Energy landscapes as a tool to integrate GPCR structure, dynamics, and function

Computational modeling of membrane proteins

The determination of membrane protein (MP) structures has always trailed that of soluble proteins due to difficulties in their overexpression, reconstitution into membrane mimetics, and subsequent structure determination. The percentage of MP structures in the protein databank (PDB) has been at a constant 1-2% for the last decade. In contrast, over half of all drugs target MPs, only highlighting how little we understand about drug-specific effects in the human body. To reduce this gap, researchers have attempted to predict structural features of MPs even before the first structure was experimentally elucidated. In this review, we present current computational methods to predict MP structure, starting with secondary structure prediction, prediction of trans-membrane spans, and topology. Even though these methods generate reliable predictions, challenges such as predicting kinks or precise beginnings and ends of secondary structure elements are still waiting to be addressed. We describe recent developments in the prediction of 3D structures of both α-helical MPs as well as β-barrels using comparative modeling techniques, de novo methods, and molecular dynamics (MD) simulations. The increase of MP structures has (1) facilitated comparative modeling due to availability of more and better templates, and (2) improved the statistics for knowledge-based scoring functions. Moreover, de novo methods have benefited from the use of correlated mutations as restraints. Finally, we outline current advances that will likely shape the field in the forthcoming decade.

Correlation of chemical shifts predicted by molecular dynamics simulations for partially disordered proteins

There has been a longstanding interest in being able to accurately predict NMR chemical shifts from structural data. Recent studies have focused on using molecular dynamics (MD) simulation data as input for improved prediction. Here we examine the accuracy of chemical shift prediction for intein systems, which have regions of intrinsic disorder. We find that using MD simulation data as input for chemical shift prediction does not consistently improve prediction accuracy over use of a static X-ray crystal structure. This appears to result from the complex conformational ensemble of the disordered protein segments. We show that using accelerated molecular dynamics (aMD) simulations improves chemical shift prediction, suggesting that methods which better sample the conformational ensemble like aMD are more appropriate tools for use in chemical shift prediction for proteins with disordered regions. Moreover, our study suggests that data accurately reflecting protein dynamics must be used as input for chemical shift prediction in order to correctly predict chemical shifts in systems with disorder.

Principles of allosteric interactions in cell signaling

Linking cell signaling events to the fundamental physicochemical basis of the conformational behavior of single molecules and ultimately to cellular function is a key challenge facing the life sciences. Here we outline the emerging principles of allosteric interactions in cell signaling, with emphasis on the following points. (1) Allosteric efficacy is not a function of the chemical composition of the allosteric pocket but reflects the extent of the population shift between the inactive and active states. That is, the allosteric effect is determined by the extent of preferred binding, not by the overall binding affinity. (2) Coupling between the allosteric and active sites does not decide the allosteric effect; however, it does define the propagation pathways, the allosteric binding sites, and key on-path residues. (3) Atoms of allosteric effectors can act as "driver" or "anchor" and create attractive "pulling" or repulsive "pushing" interactions. Deciphering, quantifying, and integrating the multiple co-occurring events present daunting challenges to our scientific community.

Molecular insights into the dynamics of pharmacogenetically important N-terminal variants of the human β2-adrenergic receptor

The human β₂-adrenergic receptor (β₂AR), a member of the G-protein coupled receptor (GPCR) family, is expressed in bronchial smooth muscle cells. Upon activation by agonists, β₂AR causes bronchodilation and relief in asthma patients. The N-terminal polymorphism of β₂AR at the 16^th position, Arg16Gly, has warranted a lot of attention since it is linked to variations in response to albuterol (agonist) treatment. Although the β₂AR is one of the well-studied GPCRs, the N-terminus which harbors this mutation, is absent in all available experimental structures. The goal of this work was to study the molecular level differences between the N-terminal variants using structural modeling and atomistic molecular dynamics simulations. Our simulations reveal that the N-terminal region of the Arg variant shows greater dynamics than the Gly variant, leading to differential placement. Further, the position and dynamics of the N-terminal region, further, affects the ligand binding-site accessibility. Interestingly, long-range effects are also seen at the ligand binding site, which is marginally larger in the Gly as compared to the Arg variant resulting in the preferential docking of albuterol to the Gly variant. This study thus reveals key differences between the variants providing a molecular framework towards understanding the variable drug response in asthma patients.

The human β₂-adrenergic receptor (β₂AR) is an important member of the GPCR family and a mutation at the 16^th position, Arg16Gly, is commonly found in the population. This variation in asthma patients is linked to differential (good/bad) response to the drug albuterol, an agonist of the β₂AR. To date, the coordinates of the N-terminal residues harboring the 16^th position mutation have not been resolved. In our study we sought to glean insights into the dynamics of the variants that could address the differential response to albuterol. We used knowledge from class A GPCRs to build the N-terminal region of β₂AR variants in conjunction with the available structure of the inactive receptor. This was followed by atomistic simulations in triplicate totaling to a sampling of 6 µs. We observe that the N-terminal region of the Arg variant is more dynamic than the Gly variant. Amongst the various differences between the variants, we observe long-range effects at the binding site leading to preferential docking of albuterol to the Gly variant. Our work is a first step to unravel the molecular mechanism linking the Arg16Gly variation to the differential response to albuterol in asthma patients.

Conformational selection and equilibrium governs the ability of retinals to bind opsin

Despite extensive study, how retinal enters and exits the visual G protein-coupled receptor rhodopsin remains unclear. One clue may lie in two openings between transmembrane helix 1 (TM1) and TM7 and between TM5 and TM6 in the active receptor structure. Recently, retinal has been proposed to enter the inactive apoprotein opsin (ops) through these holes when the receptor transiently adopts the active opsin conformation (ops*). Here, we directly test this "transient activation" hypothesis using a fluorescence-based approach to measure rates of retinal binding to samples containing differing relative fractions of ops and ops*. In contrast to what the transient activation hypothesis model would predict, we found that binding for the inverse agonist, 11-cis-retinal (11CR), slowed when the sample contained more ops* (produced using M257Y, a constitutively activating mutation). Interestingly, the increased presence of ops* allowed for binding of the agonist, all-trans-retinal (ATR), whereas WT opsin showed no binding. Shifting the conformational equilibrium toward even more ops* using a G protein peptide mimic (either free in solution or fused to the receptor) accelerated the rate of ATR binding and slowed 11CR binding. An arrestin peptide mimic showed little effect on 11CR binding; however, it stabilized opsin · ATR complexes. The TM5/TM6 hole is apparently not involved in this conformational selection. Increasing its size by mutagenesis did not enable ATR binding but instead slowed 11CR binding, suggesting that it may play a role in trapping 11CR. In summary, our results indicate that conformational selection dictates stable retinal binding, which we propose involves ATR and 11CR binding to different states, the latter a previously unidentified, open-but-inactive conformation.

Biological redundancy of endogenous GPCR ligands in the gut and the potential for endogenous functional selectivity

This review focuses on the existence and function of multiple endogenous agonists of the somatostatin and opioid receptors with an emphasis on their expression in the gastrointestinal tract. These agonists generally arise from the proteolytic cleavage of prepropeptides during peptide maturation or from degradation of peptides by extracellular or intracellular endopeptidases. In other examples, endogenous peptide agonists for the same G protein-coupled receptors can be products of distinct genes but contain high sequence homology. This apparent biological redundancy has recently been challenged by the realization that different ligands may engender distinct receptor conformations linked to different intracellular signaling profiles and, as such the existence of distinct ligands may underlie mechanisms to finely tune physiological responses. We propose that further characterization of signaling pathways activated by these endogenous ligands will provide invaluable insight into the mechanisms governing biased agonism. Moreover, these ligands may prove useful in the design of novel therapeutic tools to target distinct signaling pathways, thereby favoring desirable effects and limiting detrimental on-target effects. Finally we will discuss the limitations of this area of research and we will highlight the difficulties that need to be addressed when examining endogenous bias in tissues and in animals.

Cholesterol modulates the dimer interface of the β₂-adrenergic receptor via cholesterol occupancy sites

The β2-adrenergic receptor is an important member of the G-protein-coupled receptor (GPCR) superfamily, whose stability and function are modulated by membrane cholesterol. The recent high-resolution crystal structure of the β2-adrenergic receptor revealed the presence of possible cholesterol-binding sites in the receptor. However, the functional relevance of cholesterol binding to the receptor remains unexplored. We used MARTINI coarse-grained molecular-dynamics simulations to explore dimerization of the β2-adrenergic receptor in lipid bilayers containing cholesterol. A novel (to our knowledge) aspect of our results is that receptor dimerization is modulated by membrane cholesterol. We show that cholesterol binds to transmembrane helix IV, and cholesterol occupancy at this site restricts its involvement at the dimer interface. With increasing cholesterol concentration, an increased presence of transmembrane helices I and II, but a reduced presence of transmembrane helix IV, is observed at the dimer interface. To our knowledge, this study is one of the first to explore the correlation between cholesterol occupancy and GPCR organization. Our results indicate that dimer plasticity is relevant not just as an organizational principle but also as a subtle regulatory principle for GPCR function. We believe these results constitute an important step toward designing better drugs for GPCR dimer targets.

A functional selectivity mechanism at the serotonin-2A GPCR involves ligand-dependent conformations of intracellular loop 2

With recent progress in determination of G protein-coupled receptor (GPCR) structure with crystallography, a variety of other experimental approaches (e.g., NMR spectroscopy, fluorescent-based assays, mass spectrometry techniques) are also being used to characterize state-specific and ligand-specific conformational states. MD simulations offer a powerful complementary approach to elucidate the dynamic features associated with ligand-specific GPCR conformations. To shed light on the conformational elements and dynamics of the important aspect of GPCR functional selectivity, we carried out unbiased microsecond-length MD simulations of the human serotonin 2A receptor (5-HT(2A)R) in the absence of ligand and bound to four distinct serotonergic agonists. The 5-HT(2A)R is a suitable system to study the structural features involved in the ligand-dependent conformational heterogeneity of GPCRs because it is well-characterized experimentally and exhibits a strong agonist-specific phenotype in that some 5-HT(2A)R agonists induce LSD-like hallucinations, while others lack this psychoactive property entirely. Here we report evidence for structural and dynamic differences in 5-HT(2A)R interacting with such pharmacologically distinct ligands, hallucinogens, and nonhallucinogens obtained from all-atom MD simulations. Differential ligand binding contacts were identified for structurally similar hallucinogens and nonhallucinogens and found to correspond to different conformations in the intracellular loop 2 (ICL2). From the different ICL2 conformations, functional selective phenotypes are suggested through effects on dimerization and/or distinct direct interaction with effector proteins. The findings are presented in the context of currently proposed hallucinogenesis mechanisms, and ICL2 is proposed as a fine-tuning selective switch that can differentiates modes of 5-HT(2A)R activation.

Structure-bias relationships for fenoterol stereoisomers in six molecular and cellular assays at the β2-adrenoceptor

Functional selectivity is well established as an underlying concept of ligand-specific signaling via G protein-coupled receptors (GPCRs). Functionally, selective drugs could show greater therapeutic efficacy and fewer adverse effects. Dual coupling of the β2-adrenoceptor (β2AR) triggers a signal transduction via Gsα and Giα proteins. Here, we examined 12 fenoterol stereoisomers in six molecular and cellular assays. Using β2AR-Gsα and β2AR-Giα fusion proteins, (R,S')- and (S,S')-isomers of 4'-methoxy-1-naphthyl-fenoterol were identified as biased ligands with preference for Gs. G protein-independent signaling via β-arrestin-2 was disfavored by these ligands. Isolated human neutrophils constituted an ex vivo model of β2AR signaling and demonstrated functional selectivity through the dissociation of cAMP accumulation and the inhibition of formyl peptide-stimulated production of reactive oxygen species. Ligand bias was calculated using an operational model of agonism and revealed that the fenoterol scaffold constitutes a promising lead structure for the development of Gs-biased β2AR agonists.

Agonists with supraphysiological efficacy at the muscarinic M2 ACh receptor

BACKGROUND AND PURPOSE

Artificial agonists may have higher efficacy for receptor activation than the physiological agonist. Until now, such 'superagonism' has rarely been reported for GPCRs. Iperoxo is an extremely potent muscarinic receptor agonist. We hypothesized that iperoxo is a 'superagonist'.

EXPERIMENTAL APPROACH

Signalling of iperoxo and newly synthesized structural analogues was compared with that of ACh at label-free M2 muscarinic receptors applying whole cell dynamic mass redistribution, measurement of G-protein activation, evaluation of cell surface agonist binding and computation of operational efficacies.

KEY RESULTS

In CHO-hM2 cells, iperoxo significantly exceeds ACh in Gi /Gs signalling competence. In the orthosteric loss-of-function mutant M2 -Y104(3.33) A, the maximum effect of iperoxo is hardly compromised in contrast to ACh. 'Superagonism' is preserved in the physiological cellular context of MRC-5 human lung fibroblasts. Structure-signalling relationships including iperoxo derivatives with either modified positively charged head group or altered tail suggest that 'superagonism' of iperoxo is mechanistically based on parallel activation of the receptor protein via two orthosteric interaction points.

CONCLUSION AND IMPLICATIONS

Supraphysiological agonist efficacy at muscarinic M2 ACh receptors is demonstrated for the first time. In addition, a possible underlying molecular mechanism of GPCR 'superagonism' is provided. We suggest that iperoxo-like orthosteric GPCR activation is a new avenue towards a novel class of receptor activators.

GPCRs: Lipid-Dependent Membrane Receptors That Act as Drug Targets

G protein-coupled receptors (GPCRs) are the largest class of molecules involved in signal transduction across cell membranes and represent major targets in the development of novel drug candidates in all clinical areas. Although there have been some recent leads, structural information on GPCRs is relatively rare due to the difficulty associated with crystallization. A specific reason for this is the intrinsic flexibility displayed by GPCRs, which is necessary for their functional diversity. Since GPCRs are integral membrane proteins, interaction of membrane lipids with them constitutes an important area of research in GPCR biology. In particular, membrane cholesterol has been reported to have a modulatory role in the function of a number of GPCRs. The role of membrane cholesterol in GPCR function is discussed with specific example of the serotonin1A receptor. Recent results show that GPCRs are characterized with structural motifs that preferentially associate with cholesterol. An emerging and important concept is oligomerization of GPCRs and its role in GPCR function and signaling. The role of membrane cholesterol in GPCR oligomerization is highlighted. Future research in GPCR biology would offer novel insight in basic biology and provide new avenues for drug discovery.

Model for coupled insertion and folding of membrane-spanning proteins

Current understanding of the forces directing the folding of integral membrane proteins is very limited compared to the detailed picture available for water-soluble proteins. While mechanistic studies of the folding process in vitro have been conducted for only a small number of membrane proteins, the available evidence indicates that their folding process is thermodynamically driven like that of soluble proteins. In vivo, however, the majority of integral membrane proteins are installed in membranes by dedicated machinery, suggesting that the cellular systems may act to facilitate and regulate the spontaneous physical process of folding. Both the in vitro folding process and the in vivo pathway must navigate an energy landscape dominated by the energetically favorable burial of hydrophobic segments in the membrane interior and the opposition to folding due to the need for passage of polar segments across the membrane. This manuscript describes a simple, exactly solvable model which incorporates these essential features of membrane protein folding. The model is used to compare the folding time under conditions which depict both the in vitro and in vivo pathways. It is proposed that the cellular complexes responsible for insertion of membrane proteins act by lowering the energy barrier for passage of polar regions through the membrane, thereby allowing the chain to more rapidly achieve the folded state.

Novel GPCR paradigms at the μ-opioid receptor

Opioids, such as morphine, are the most clinically useful class of analgesic drugs for the treatment of acute and chronic pain. However, the use of opioids is greatly limited by the development of severe adverse side effects. Consequently, drug discovery efforts have been directed towards improving the therapeutic profile of opioid-based treatments. Opioid receptors are members of the family of GPCRs. As such, the recent GPCR paradigms of biased agonism and allosterism may provide novel avenues for more effective analgesics. Biased agonism (or functional selectivity) has been described for all the opioid receptor family members. Furthermore, the first allosteric modulators of opioid receptors have very recently been described. However, identification and quantification of biased agonism in a manner that is informative to medicinal chemists and drug discovery programmes still remains a challenge. In this review, we examine the progress, to date, towards identification and quantification of biased agonism and allosterism at the μ-opioid receptor in the context of its implications for the discovery of better and safer analgesics.

LINKED ARTICLES

This article is part of a themed section on Opioids: New Pathways to Functional Selectivity. To view the other articles in this section visit http://dx.doi.org/10.1111/bph.2015.172.issue-2.

Regulation of phospholipase C-β(1) GTPase-activating protein (GAP) function and relationship to G(q) efficacy

How cells regulate Gq efficacy (initiation and termination of Gq signaling) to effect response remains a central question in pharmacology and drug discovery. Phospholipase C-β1 (PLC-β1) is an effector and a GTPase activating protein (GAP) specific to Gαq. The physiological function of PLC-β1 GAP remains unclear and controversial. GAPs are generally thought to function in deactivation of Gq signaling. However, PLC-β1 GAP has also been shown to increase signaling efficiency through kinetic coupling with the ligand-activated GPCR. GPCRs function as guanine nucleotide exchange factors (GEF) on the G protein activation cycle. This article sets forth a new hypothesis that could unify these conflicting paradigms as it pertains to physiological signaling and native levels of protein. It is proposed that the physiological function of PLC-β1 GAP is context-dependent and regulated by phosphatidic acid (PA). PA stimulates PLC-β1 GAP activity. In the absence of ligand, PLC-β1 GAP does indeed deactivate Gq signaling, limiting leaky activation to set the threshold for stimulation to sharpen signal kinetics. However in the presence of activating ligand, the increase in levels of PA would stimulate PLC-β1 GAP to kinetically couple with GPCR GEF to increase signaling efficiency. We found that PA-increased Gq efficiency is dependent on signaling via the unique PLC-β1 PA binding domain.

Unraveling structural mechanisms of allosteric drug action

Orthosteric drugs block the active site to obstruct function; allosteric drugs modify the population of the active state, to modulate function. Available data lead us to propose that allosteric drugs can constitute anchors and drivers. The anchor docks into an allosteric pocket. The conformation with which it interacts is unchanged during the transition between the inactive and active states. The anchor provides the foundation that allows the driver to exert a 'pull' and/or 'push' action that shifts the receptor population from the inactive to the active state. The presence or absence of driver atom in an allosteric drug can exert opposite agonism. We map a strategy for driver identification and expect the allosteric trigger concept to transform agonist/antagonist drug discovery.

Free energy landscape of G-protein coupled receptors, explored by accelerated molecular dynamics

G-protein coupled receptors (GPCRs) mediate cellular responses to various hormones and neurotransmitters and are important targets for treating a wide spectrum of diseases. They are known to adopt multiple conformational states (e.g., inactive, intermediate and active) during their modulation of various cell signaling pathways. Here, the free energy landscape of GPCRs is explored using accelerated molecular dynamics (aMD) simulations as demonstrated on the M2 muscarinic receptor, a key GPCR that regulates human heart rate and contractile forces of cardiomyocytes. Free energy profiles of important structural motifs that undergo conformational transitions upon GPCR activation and allosteric signaling are analyzed in detail, including the Arg(3.50)-Glu(6.30) ionic lock, the Trp(6.48) toggle switch and the hydrogen interactions between Tyr(5.58)-Tyr(7.53).

Recent developments in biased agonism

The classic paradigm of G protein-coupled receptor (GPCR) activation was based on the understanding that agonist binding to a receptor induces or stabilizes a conformational change to an 'active' conformation. In the past decade, however, it has been appreciated that ligands can induce distinct 'active' receptor conformations with unique downstream functional signaling profiles. Building on the initial recognition of the existence of such 'biased ligands', recent years have witnessed significant developments in several areas of GPCR biology. These include increased understanding of structural and biophysical mechanisms underlying biased agonism, improvements in characterization and quantification of ligand efficacy, as well as clinical development of these novel ligands. Here we review recent major developments in these areas over the past several years.

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

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

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.

β-adrenoceptors as molecular targets in the treatment of hypertension

Regulation of sympathoadrenal activity has been a long-time target in the management of hypertension. Regulation of β-adrenoceptor (βAR) function has been the most therapeutically important of these targets. The development of effective antihypertensive treatments based on βAR antagonism paralleled the elucidation of the molecular basis of β-adrenergic effects by the family of βARs, which are members of the G-protein-coupled receptor (GPCR) superfamily. βARs serve as the extracellular face of the transmembrane signalling pathway that results in the consequent activation of heterotrimeric G-proteins and the activation of several other newly appreciated signalling molecules that include β-arrestins and GPCR kinases (GRKs). The aggregate effect of the activation of these signalling pathways mediates the response to βAR activation. Paradoxically, the hypertensive state is characterized by impaired βAR responsiveness. This defect is common to many other receptor systems linked to the stimulator G protein (Gs) and adenylyl cyclase activation. This impairment is principally mediated by receptor-G-protein uncoupling, which has been linked to increased expression and activity of GRK2.

Mapping of allosteric druggable sites in activation-associated conformers of the M2 muscarinic receptor

G-protein-coupled receptors (GPCRs) are key cellular signaling proteins and have been targeted by approximately 30-40% of marketed drugs for treating many human diseases including cancer and heart failure. Recently, we directly observed activation of the M2 muscarinic receptor through long-timescale accelerated molecular dynamics (aMD) simulation, which revealed distinct inactive, intermediate and active conformers of the receptor. Here, FTMAP is applied to search for 'hot spots' in these activation-associated conformers using a library of 16 organic probe molecules that represent fragments of potential drugs. Seven allosteric (non-orthosteric) binding sites are identified in the M2 receptor through the FTMAP analysis. These sites are distributed in the solvent-exposed extracellular and intracellular mouth regions, as well as the lipid-exposed pockets formed by the transmembrane α helices TM3-TM4, TM5-TM6 and TM7-TM1/TM2. They serve as promising target sites for designing novel allosteric modulators as receptor-selective drugs.

A unified view of "how allostery works"

The question of how allostery works was posed almost 50 years ago. Since then it has been the focus of much effort. This is for two reasons: first, the intellectual curiosity of basic science and the desire to understand fundamental phenomena, and second, its vast practical importance. Allostery is at play in all processes in the living cell, and increasingly in drug discovery. Many models have been successfully formulated, and are able to describe allostery even in the absence of a detailed structural mechanism. However, conceptual schemes designed to qualitatively explain allosteric mechanisms usually lack a quantitative mathematical model, and are unable to link its thermodynamic and structural foundations. This hampers insight into oncogenic mutations in cancer progression and biased agonists' actions. Here, we describe how allostery works from three different standpoints: thermodynamics, free energy landscape of population shift, and structure; all with exactly the same allosteric descriptors. This results in a unified view which not only clarifies the elusive allosteric mechanism but also provides structural grasp of agonist-mediated signaling pathways, and guides allosteric drug discovery. Of note, the unified view reasons that allosteric coupling (or communication) does not determine the allosteric efficacy; however, a communication channel is what makes potential binding sites allosteric.

Divergent label-free cell phenotypic pharmacology of ligands at the overexpressed β₂-adrenergic receptors

We present subclone sensitive cell phenotypic pharmacology of ligands at the β₂-adrenergic receptor (β₂-AR) stably expressed in HEK-293 cells. The parental cell line was transfected with green fluorescent protein (GFP)-tagged β₂-AR. Four stable subclones were established and used to profile a library of sixty-nine AR ligands. Dynamic mass redistribution (DMR) profiling resulted in a pharmacological activity map suggesting that HEK293 endogenously expresses functional G_i-coupled α₂-AR and G_s-coupled β₂-AR, and the label-free cell phenotypic activity of AR ligands are subclone dependent. Pathway deconvolution revealed that the DMR of epinephrine is originated mostly from the remodeling of actin microfilaments and adhesion complexes, to less extent from the microtubule networks and receptor trafficking, and certain agonists displayed different efficacy towards the cAMP-Epac pathway. We demonstrate that receptor signaling and ligand pharmacology is sensitive to the receptor expression level, and the organization of the receptor and its signaling circuitry.

Beyond standard molecular dynamics: investigating the molecular mechanisms of G protein-coupled receptors with enhanced molecular dynamics methods

The majority of biological processes mediated by G Protein-Coupled Receptors (GPCRs) take place on timescales that are not conveniently accessible to standard molecular dynamics (MD) approaches, notwithstanding the current availability of specialized parallel computer architectures, and efficient simulation algorithms. Enhanced MD-based methods have started to assume an important role in the study of the rugged energy landscape of GPCRs by providing mechanistic details of complex receptor processes such as ligand recognition, activation, and oligomerization. We provide here an overview of these methods in their most recent application to the field.

Constitutive activities in the thyrotropin receptor: regulation and significance

The thyroid-stimulating hormone receptor (TSHR, or thyrotropin receptor) is a family A G protein-coupled receptor. It not only binds thyroid-stimulating hormone (TSH, or thyrotropin) but also interacts with autoantibodies under pathological conditions. The TSHR and TSH are essential for thyroid growth and function and thus for all thyroid hormone-associated physiological superordinated processes, including metabolism and development of the central nervous system. In vitro studies have found that the TSHR permanently stimulates ligand-independent (constitutive) activation of Gs, which ultimately leads to intracellular cAMP accumulation. Furthermore, a vast variety of constitutively activating mutations of TSHR-at more than 50 different amino acid positions-have been reported to enhance basal signaling. These lead in vivo to a "gain-of-function" phenotype of nonautoimmune hyperthyroidism or toxic adenomas. Moreover, many naturally occurring inactivating mutations are known to cause a "loss-of-function" phenotype, resulting in resistance to thyroid hormone or hyperthyrotropinemia. Several of these mutations are also characterized by impaired basal signaling, and these are designated here as "constitutively inactivating mutations" (CIMs). More than 30 amino acid positions with CIMs have been identified so far. Moreover, the permanent TSHR signaling capacity can also be blocked by inverse agonistic antibodies or small drug-like molecules, which both have a potential for clinical usage. In this chapter, information on constitutive activity in the TSHR is described, including up- and downregulation, linked protein conformations, physiological and pathophysiological conditions, and related intracellular signaling.

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.

Opsin, a structural model for olfactory receptors?

Receptor-ligand interaction: Olfactory receptors (ORs) are G-protein-coupled receptors (GPCRs), which detect signaling molecules such as hormones and odorants. The structure of opsin, the GPCR employed in vision, with a detergent molecule bound deep in its orthosteric ligand-binding pocket provides a template for OR homology modeling, thus enabling investigation of the structural basis of the mechanism of odorant-receptor recognition.

A constitutively activating mutation alters the dynamics and energetics of a key conformational change in a ligand-free G protein-coupled receptor

G protein-coupled receptors (GPCRs) undergo dynamic transitions between active and inactive conformations. Usually, these conversions are triggered when the receptor detects an external signal, but some so-called constitutively activating mutations, or CAMs, induce a GPCR to bind and activate G proteins in the absence of external stimulation, in ways still not fully understood. Here, we investigated how a CAM alters the structure of a GPCR and the dynamics involved as the receptor transitions between different conformations. Our approach used site-directed fluorescence labeling (SDFL) spectroscopy to compare opsin, the ligand-free form of the GPCR rhodopsin, with opsin containing the CAM M257Y, focusing specifically on key movements that occur in the sixth transmembrane helix (TM6) during GPCR activation. The site-directed fluorescence labeling data indicate opsin is constrained to an inactive conformation both in detergent micelles and lipid membranes, but when it contains the M257Y CAM, opsin is more dynamic and can interact with a G protein mimetic. Further study of these receptors using tryptophan-induced quenching (TrIQ) methods indicates that in detergent, the CAM significantly increases the population of receptors in the active state, but not in lipids. Subsequent Arrhenius analysis of the TrIQ data suggests that, both in detergent and lipids, the CAM lowers the energy barrier for TM6 movement, a key transition required for conversion between the inactive and active conformations. Together, these data suggest that the lowered energy barrier is a primary effect of the CAM on the receptor dynamics and energetics.

Precision vs flexibility in GPCR signaling

The G protein coupled receptor (GPCR) rhodopsin activates the heterotrimeric G protein transducin (Gt) to transmit the light signal into retinal rod cells. The rhodopsin activity is virtually zero in the dark and jumps by more than one billion fold after photon capture. Such perfect switching implies both high fidelity and speed of rhodopsin/Gt coupling. We employed Fourier transform infrared (FTIR) spectroscopy and supporting all-atom molecular dynamics (MD) simulations to study the conformational diversity of rhodopsin in membrane environment and extend the static picture provided by the available crystal structures. The FTIR results show how the equilibria of inactive and active protein states of the receptor (so-called metarhodopsin states) are regulated by the highly conserved E(D)RY and Yx7K(R) motives. The MD data identify an intrinsically unstructured cytoplasmic loop region connecting transmembrane helices 5 and 6 (CL3) and show how each protein state is split into conformational substates. The C-termini of the Gtγ- and Gtα-subunits (GαCT and GγCT), prepared as synthetic peptides, are likely to bind sequentially and at different sites of the active receptor. The peptides have different effects on the receptor conformation. While GγCT stabilizes the active states but preserves CL3 flexibility, GαCT selectively stabilizes a single conformational substate with largely helical CL3, as it is found in crystal structures. Based on these results we propose a mechanism for the fast and precise signal transfer from rhodopsin to Gt, which assumes a stepwise and mutual reduction of their conformational space. The mechanism relies on conserved amino acids and may therefore underlie GPCR/G protein coupling in general.

Functional selectivity of G-protein-coupled receptors: from recombinant systems to native human cells

In the mid 1990s, it was assumed that a two-state model, postulating an inactive (R) state and an active (R*) state provides the molecular basis for GPCR activation. However, it became clear that this model could not accommodate many experimental observations. Accordingly, the two-state model was superseded by a multi-state model according to which any given ligand stabilizes a unique receptor conformation with distinct capabilities of activating down-stream G-proteins and β-arrestin. Much of this research was conducted with the β2-adrenoceptor in recombinant systems. At the molecular level, there is now no doubt anymore that ligand-specific receptor conformations, also referred to as functional selectivity, exist. This concept holds great potential for drug discovery in terms of developing drugs with higher selectivity for specific cells and/or cell functions and fewer side effects. A major challenge is the analysis for functional selectivity in native cells. Here, I discuss our current knowledge on functional selectivity of three representative GPCRs, the β2-adrenoceptor and the histamine H2- and H4-receptors, in recombinant systems and native human cells. Studies with human neutrophils and eosinophils support the concept of functional selectivity. A major strategy for the analysis of functional selectivity in native cells is to generate complete concentration/response curves with a large set of structurally diverse ligands for multiple parameters. Next, correlations of potencies and efficacies are analyzed, and deviations of the correlations from linearity are indicative for functional selectivity. Additionally, pharmacological inhibitors are used to dissect cell functions from each other.

Simulations of biased agonists in the β(2) adrenergic receptor with accelerated molecular dynamics

The biased agonism of the G protein-coupled receptors (GPCRs), where in addition to a traditional G protein-signaling pathway a GPCR promotes intracellular signals though β-arrestin, is a novel paradigm in pharmacology. Biochemical and biophysical studies have suggested that a GPCR forms a distinct ensemble of conformations signaling through the G protein and β-arrestin. Here we report on the dynamics of the β2 adrenergic receptor bound to the β-arrestin and G protein-biased agonists and the empty receptor to further characterize the receptor conformational changes caused by biased agonists. We use conventional and accelerated molecular dynamics (aMD) simulations to explore the conformational transitions of the GPCR from the active state to the inactive state. We found that aMD simulations enable monitoring of the transition within the nanosecond time scale while capturing the known microscopic characteristics of the inactive states, such as the ionic lock, the inward position of F6.44, and water clusters. Distinct conformational states are shown to be stabilized by each biased agonist. In particular, in simulations of the receptor with the β-arrestin-biased agonist N-cyclopentylbutanepherine, we observe a different pattern of motions in helix 7 when compared to simulations with the G protein-biased agonist salbutamol that involves perturbations of the network of interactions within the NPxxY motif. Understanding the network of interactions induced by biased ligands and the subsequent receptor conformational shifts will lead to development of more efficient drugs.

How ligands and signalling proteins affect G-protein-coupled receptors' conformational landscape

The dynamic character of GPCRs (G-protein-coupled receptors) is essential to their function. However, the details of how ligands and signalling proteins stabilize a receptor conformation to trigger the activation of a given signalling pathway remain largely unexplored. Multiple data, including recent results obtained with the purified ghrelin receptor, suggest a model where ligand efficacy and functional selectivity are directly related to different receptor conformations. Importantly, distinct effector proteins (G-proteins and arrestins) as well as ligands are likely to affect the conformational landscape of GPCRs in different manners, as we show with the isolated ghrelin receptor. Such modulation of the GPCR conformational landscape by pharmacologically distinct ligands and effector proteins has major implications for the design of new drugs that activate specific signalling pathways.

The concept of allosteric interaction and its consequences for the chemistry of the brain

Throughout this Reflections article, I have tried to follow up on the genesis in the 1960s and subsequent evolution of the concept of allosteric interaction and to examine its consequences within the past decades, essentially in the field of the neuroscience. The main conclusion is that allosteric mechanisms built on similar structural principles operate in bacterial regulatory enzymes, gene repressors (and the related nuclear receptors), rhodopsin, G-protein-coupled receptors, neurotransmitter receptors, ion channels, and so on from prokaryotes up to the human brain yet with important features of their own. Thus, future research on these basic cybernetic sensors is expected to develop in two major directions: at the elementary level, toward the atomic structure and molecular dynamics of the conformational changes involved in signal recognition and transduction, but also at a higher level of organization, the contribution of allosteric mechanisms to the modulation of brain functions.

Single-molecule force spectroscopy of G-protein-coupled receptors

The applicability of single-molecule force spectroscopy (SMFS) to characterize membrane proteins in vitro is developing rapidly and opening a wide range of fascinating possibilities to study how intra- and intermolecular interactions determine their structural stability and functional state. In particular, understanding how molecular interactions contribute to the functional state of G-protein-coupled receptors (GPCRs) is of importance because they mediate most of our physiological responses and act as therapeutic targets for a broad spectrum of diseases. In our review we focus on SMFS to characterize GPCRs embedded in their physiologically relevant membranes and exposed to physiologically relevant conditions. SMFS uses a molecularly sharp stylus to grasp the terminal end of a GPCR and to quickly unfold the receptor while recording interaction forces. The positional accuracy of SMFS localizes these interactions to structural segments of the GPCR whereas the sensitivity of SMFS enables their stabilizing interaction forces to be quantified. To further investigate the kinetic, energetic and mechanical properties of the structural segments, dynamic SMFS (DFS) probes their stability over a wide range of loading rates. These parameters provide insight into the energy landscape that provides information on the structural and functional properties of the GPCRs. Selected highlights exemplify the application of SMFS to characterize inter- and intramolecular interactions, which change the properties of GPCRs in relation to their functional state (e.g., ligand binding), diseased state (e.g., mutation), or lipid environment such as cholesterol.

The role of ligands on the equilibria between functional states of a G protein-coupled receptor

G protein-coupled receptors exhibit a wide variety of signaling behaviors in response to different ligands. When a small label was incorporated on the cytosolic interface of transmembrane helix 6 (Cys-265), (19)F NMR spectra of the β2 adrenergic receptor (β2AR) reconstituted in maltose/neopentyl glycol detergent micelles revealed two distinct inactive states, an activation intermediate state en route to activation, and, in the presence of a G protein mimic, a predominant active state. Analysis of the spectra as a function of temperature revealed that for all ligands, the activation intermediate is entropically favored and enthalpically disfavored. β2AR enthalpy changes toward activation are notably lower than those observed with rhodopsin, a likely consequence of basal activity and the fact that the ionic lock and other interactions stabilizing the inactive state of β2AR are weaker. Positive entropy changes toward activation likely reflect greater mobility (configurational entropy) in the cytoplasmic domain, as confirmed through an order parameter analysis. Ligands greatly influence the overall changes in enthalpy and entropy of the system and the corresponding changes in population and amplitude of motion of given states, suggesting a complex landscape of states and substates.

Production of GPCR and GPCR complexes for structure determination

Since the first high-resolution structure of the beta 2 adrenergic receptor (b2AR) in 2007, we have seen a growing number of G-protein-coupled receptor (GPCR) structures coming to the repertory, providing a significant progress in our understanding of the structural basis of their function. This has been achieved by the interdisciplinary collaborative work between scientists with various expertise and the development of new methodologies as well as combining and optimizing existing techniques.