The role of internal water molecules in the structure and function of the rhodopsin family of G protein-coupled receptors

Triptorelin-functionalized PEG-coated biosynthesized gold nanoparticles: Effects of receptor-ligand interactions on adhesion to triple negative breast cancer cells

This paper presents the results of an experimental and computational study of the adhesion of triptorelin-conjugated PEG-coated biosynthesized gold nanoparticles (GNP-PEG-TRP) to triple-negative breast cancer (TNBC) cells. The adhesion is studied at the nanoscale using a combination of atomic force microscopy (AFM) experiments and molecular dynamics (MD) simulations. The AFM measurements showed that the triptorelin-functionalized gold nanoparticles (GNP-TRP and GNP-PEG-TRP) have higher adhesion to triple-negative breast cancer cells (TNBC) than non-tumorigenic breast cells. The increased adhesion of GNP-TRP and GNP-PEG-TRP to TNBC is also attributed to the overexpression of LHRH receptors on the surfaces of both TNBC. Finally, the molecular dynamics model reveals insights into the effects of receptor density, molecular configuration, and receptor-ligand docking characteristics on the interactions of triptorelin-functionalized PEG-coated gold nanoparticles with TNBC. A three to nine-fold increase in the adhesion is predicted between triptorelin-functionalized PEG-coated gold nanoparticles and TNBC cells. The implications of the results are then discussed for the specific targeting of TNBC.

Gonadotropin-Releasing Hormone and GnRH Receptor: Structure, Function and Drug Development

BACKGROUND

Gonadotropin-Releasing Hormone (GnRH) is a key element in sexual maturation and regulation of the reproductive cycle in the human organism. GnRH interacts with the pituitary cells through the activation of the Gonadotropin Releasing Hormone Receptors (GnRHR). Any impairments/dysfunctions of the GnRH-GnRHR complex lead to the development of various cancer types and disorders. Furthermore, the identification of GnRHR as a potential drug target has led to the development of agonist and antagonist molecules implemented in various treatment protocols. The development of these drugs was based on the information derived from the functional studies of GnRH and GnRHR.

OBJECTIVE

This review aims at shedding light on the versatile function of GnRH and GnRH receptor and offers an apprehensive summary regarding the development of different agonists, antagonists and non-peptide GnRH analogues.

CONCLUSION

The information derived from these studies can enhance our understanding of the GnRH-GnRHR versatile nature and offer valuable insight into the design of new more potent molecules.

Experimental and computational analysis of biased agonism on full-length and a C-terminally truncated adenosine A2A receptor

Biased agonism, the ability of agonists to differentially activate downstream signaling pathways by stabilizing specific receptor conformations, is a key issue for G protein-coupled receptor (GPCR) signaling. The C-terminal domain might influence this functional selectivity of GPCRs as it engages G proteins, GPCR kinases, β-arrestins, and several other proteins. Thus, the aim of this paper is to compare the agonist-dependent selectivity for intracellular pathways in a heterologous system expressing the full-length (A_2AR) and a C-tail truncated (A_2A^Δ40R lacking the last 40 amino acids) adenosine A_2A receptor, a GPCR that is already targeted in Parkinson’s disease using a first-in-class drug. Experimental data such as ligand binding, cAMP production, β-arrestin recruitment, ERK1/2 phosphorylation and dynamic mass redistribution assays, which correspond to different aspects of signal transduction, were measured upon the action of structurally diverse compounds (the agonists adenosine, NECA, CGS-21680, PSB-0777 and LUF-5834 and the SCH-58261 antagonist) in cells expressing A_2AR and A_2A^Δ40R. The results show that taking cAMP levels and the endogenous adenosine agonist as references, the main difference in bias was obtained with PSB-0777 and LUF-5834. The C-terminus is dispensable for both G-protein and β-arrestin recruitment and also for MAPK activation. Unrestrained molecular dynamics simulations, at the μs timescale, were used to understand the structural arrangements of the binding cavity, triggered by these chemically different agonists, facilitating G protein binding with different efficacy.

HomolWat: a web server tool to incorporate 'homologous' water molecules into GPCR structures

Internal water molecules play an essential role in the structure and function of membrane proteins including G protein-coupled receptors (GPCRs). However, technical limitations severely influence the number and certainty of observed water molecules in 3D structures. This may compromise the accuracy of further structural studies such as docking calculations or molecular dynamics simulations. Here we present HomolWat, a web application for incorporating water molecules into GPCR structures by using template-based modelling of homologous water molecules obtained from high-resolution structures. While there are various tools available to predict the positions of internal waters using energy-based methods, the approach of borrowing lacking water molecules from homologous GPCR structures makes HomolWat unique. The tool can incorporate water molecules into a protein structure in about a minute with around 85% of water recovery. The web server is freely available at http://lmc.uab.es/homolwat.

QTY code designed thermostable and water-soluble chimeric chemokine receptors with tunable ligand affinity

Chemokine receptors are of great interest as they play a critical role in many immunological and pathological processes. The ability to study chemokine receptors in aqueous solution without detergent would be significant because natural receptors require detergents to become soluble. We previously reported using the QTY code to design detergent-free chemokine receptors. We here report the design of 2 detergent-free chimeric chemokine receptors that were experimentally unattainable in detergent solution. We designed chimeric receptors by switching the N terminus and 3 extracellular (EC) loops between different receptors. Specifically, we replaced the N terminus and 3 EC loops of CCR5^QTY with the N terminus and 3 EC loops of CXCR4. The ligand for CXCR4; namely CXCL12, binds to the chimeric receptor CCR5^QTY (7TM)-CXCR4 (N terminus+3 EC loops), but with lower affinity compared to CXCR4; the CCL5 ligand of CCR5 binds the chimeric receptor with ∼20× lower affinity. The chimeric design helps to elucidate the mechanism of native receptor-ligand interaction. We also show that all detergent-free QTY-designed chemokine receptors, expressed in Escherichia coli, bind to their respective chemokines with affinities in the nanomolar (nM) range, similar to the affinities of native receptors and SF9-produced QTY variants. These QTY-designed receptors exhibit remarkable thermostability in the presence of arginine and retain ligand-binding activity after heat treatment at 60 °C for 4 h and 24 h, and at 100 °C for 10 min. Our design approach enables affordable scale-up production of detergent-free QTY variant chemokine receptors with tunable functionality for various uses.

The role of water molecules in phototransduction of retinal proteins and G protein-coupled receptors

G protein coupled receptors (GPCRs) are a key family of membrane proteins in all eukaryotes and also very important drug targets for medical intervention. The extensively studied visual pigment rhodopsin is a prime example of a family A GPCR. Its chromophore ligand retinal is covalently linked to a lysine in helix seven forming a protonated Schiff base. Interestingly, this is the same situation in other-non-GPCR-retinal proteins, like the prototype light-driven microbial proton pump bacteriorhodopsin, albeit there is no (or only a very remote) phylogenetical link. Close to the retinal ligand, several water molecules help to organise a functionally important hydrogen bond network that undergoes significant changes during photo-activation. Such water-mediated networks are also critical for ligand binding of other GPCRs and they are becoming increasingly important in drug discovery. GPCRs also contain a partially conserved water mediated hydrogen bond network that stabilises the ground state of the receptor, and rearrangement of this network leads to the stabilization of the active state. Some water molecules have a specific role in this process to appropriately orient specific residues relative to the Schiff base, and to modulate the fine structure of the transmembrane bundle, particularly near the intracellular G protein binding site. While the atomic details of these mechanisms are still missing, the recent developments in free electron lasers (FELs) are enabling us to begin to observe the changes in waters and relevant side chains shortly after photo activation at an unprecedented level of spatial and temporal resolution.

Understanding Ligand Binding to G-Protein Coupled Receptors Using Multiscale Simulations

Human G-protein coupled receptors (GPCRs) convey a wide variety of extracellular signals inside the cell and they are one of the main targets for pharmaceutical intervention. Rational drug design requires structural information on these receptors; however, the number of experimental structures is scarce. This gap can be filled by computational models, based on homology modeling and docking techniques. Nonetheless, the low sequence identity across GPCRs and the chemical diversity of their ligands may limit the quality of these models and hence refinement using molecular dynamics simulations is recommended. This is the case for olfactory and bitter taste receptors, which constitute the first and third largest GPCR groups and show sequence identities with the available GPCR templates below 20%. We have developed a molecular dynamics approach, based on the combination of molecular mechanics and coarse grained (MM/CG), tailored to study ligand binding in GPCRs. This approach has been applied so far to bitter taste receptor complexes, showing significant predictive power. The protein/ligand interactions observed in the simulations were consistent with extensive mutagenesis and functional data. Moreover, the simulations predicted several binding residues not previously tested, which were subsequently verified by carrying out additional experiments. Comparison of the simulations of two bitter taste receptors with different ligand selectivity also provided some insights into the binding determinants of bitter taste receptors. Although the MM/CG approach has been applied so far to a limited number of GPCR/ligand complexes, the excellent agreement of the computational models with the mutagenesis and functional data supports the applicability of this method to other GPCRs for which experimental structures are missing. This is particularly important for the challenging case of GPCRs with low sequence identity with available templates, for which molecular docking shows limited predictive power.

Heterogeneity of water structure and dynamics at the protein-water interface

In this molecular dynamics simulation study, we analyze the local structural and dynamic properties of water hydrating the protein ubiquitin on a spatial grid with 1 Å resolution. This allows for insights into the spatial distribution of water number densities, molecular orientations, translations, and rotations as a function of distance from the protein surface. Water molecule orientations follow a heterogeneous distribution with preferred local orientations of water dipoles and O-H bond vectors up to 10-15 Å distances from the protein, while local variations of the water number density converge to homogeneous bulk-like values within less than 8 Å. Interestingly, we find that the long-ranged orientational structure of water does not impact either the translational or rotational dynamics of water. Instead, heterogeneous distributions of local dynamical parameters and averaged dynamical retardation factors are only found close to the protein surface and follow a distance dependence comparable to heterogeneities in the local water number density. This study shows that the formation of nanodomains of preferred water orientations far from the protein does not significantly impact dynamical processes probed as a non-local average in most experiments.

Mechanisms Underlying Allosteric Molecular Switches of Metabotropic Glutamate Receptor 5

The metabotropic glutamate 5 (mGlu₅) receptor is a class C G protein-coupled receptor (GPCR) that is implicated in several CNS disorders making it a popular drug discovery target. Years of research have revealed allosteric mGlu₅ ligands showing an unexpected complete switch in functional activity despite only small changes in their chemical structure, resulting in positive allosteric modulators (PAM) or negative allosteric modulators (NAM) for the same scaffold. Up to now, the origins of this effect are not understood, causing difficulties in a drug discovery context. In this work, experimental data was gathered and analyzed alongside docking and Molecular Dynamics (MD) calculations for three sets of PAM and NAM pairs. The results consistently show the role of specific interactions formed between ligand substituents and amino acid side chains that block or promote local movements associated with receptor activation. The work provides an explanation for how such small structural changes lead to remarkable differences in functional activity. While this work can greatly help drug discovery programs avoid these switches, it also provides valuable insight into the mechanisms of class C GPCR allosteric activation. Furthermore, the approach shows the value of applying MD to understand functional activity in drug design programs, even for such close structural analogues.

Diverse GPCRs exhibit conserved water networks for stabilization and activation

G protein-coupled receptors (GPCRs) represent both the largest class of drug targets and the largest family of human membrane proteins. Recent three-dimensional structures have revealed the presence of many water molecules buried inside GPCRs, but the functional role of these water molecules has been unclear. Using extensive atomic-level computer simulations, we find that, although most of these waters are highly mobile, a few are stable. These stable water molecules mediate state-dependent and state-independent polar networks that are conserved across diverse GPCRs. These findings promise to help guide the rational design of GPCR-targeted drugs.

G protein-coupled receptors (GPCRs) have evolved to recognize incredibly diverse extracellular ligands while sharing a common architecture and structurally conserved intracellular signaling partners. It remains unclear how binding of diverse ligands brings about GPCR activation, the common structural change that enables intracellular signaling. Here, we identify highly conserved networks of water-mediated interactions that play a central role in activation. Using atomic-level simulations of diverse GPCRs, we show that most of the water molecules in GPCR crystal structures are highly mobile. Several water molecules near the G protein-coupling interface, however, are stable. These water molecules form two kinds of polar networks that are conserved across diverse GPCRs: (i) a network that is maintained across the inactive and the active states and (ii) a network that rearranges upon activation. Comparative analysis of GPCR crystal structures independently confirms the striking conservation of water-mediated interaction networks. These conserved water-mediated interactions near the G protein-coupling region, along with diverse water-mediated interactions with extracellular ligands, have direct implications for structure-based drug design and GPCR engineering.

Structural Connection between Activation Microswitch and Allosteric Sodium Site in GPCR Signaling

Sodium ions are endogenous allosteric modulators of many G-protein-coupled receptors (GPCRs). Mutation of key residues in the sodium binding motif causes a striking effect on G-protein signaling. We report the crystal structures of agonist complexes for two variants in the first sodium coordination shell of the human A_2A adenosine receptor, D52^2.50N and S91^3.39A. Both structures present an overall active-like conformation; however, the variants show key changes in the activation motif NPxxY. Changes in the hydrogen bonding network in this microswitch suggest a possible mechanism for modified G-protein signaling and enhanced thermal stability. These structures, signaling data, and thermal stability analysis with a panel of pharmacological ligands provide a basis for understanding the role of the sodium-coordinating residues on stability and G-protein signaling. Utilizing the D^2.50N variant is a promising method for stabilizing class A GPCRs to accelerate structural efforts and drug discovery.

Intracellular Transfer of Na+ in an Active-State G-Protein-Coupled Receptor

Playing a central role in cell signaling, G-protein-coupled receptors (GPCRs) are the largest superfamily of membrane proteins and form the majority of drug targets in humans. How extracellular agonist binding triggers the activation of GPCRs and associated intracellular effector proteins remains, however, poorly understood. Structural studies have revealed that inactive class A GPCRs harbor a conserved binding site for Na+ ions in the center of their transmembrane domain, accessible from the extracellular space. Here, we show that the opening of a conserved hydrated channel in the activated state receptors allows the Na+ ion to egress from its binding site into the cytosol. Coupled with protonation changes, this ion movement occurs without significant energy barriers, and can be driven by physiological transmembrane ion and voltage gradients. We propose that Na+ ion exchange with the cytosol is a key step in GPCR activation. Further, we hypothesize that this transition locks receptors in long-lived active-state conformations.

Evolutionary action and structural basis of the allosteric switch controlling β2AR functional selectivity

Functional selectivity of G-protein-coupled receptors is believed to originate from ligand-specific conformations that activate only subsets of signaling effectors. In this study, to identify molecular motifs playing important roles in transducing ligand binding into distinct signaling responses, we combined in silico evolutionary lineage analysis and structure-guided site-directed mutagenesis with large-scale functional signaling characterization and non-negative matrix factorization clustering of signaling profiles. Clustering based on the signaling profiles of 28 variants of the β₂-adrenergic receptor reveals three clearly distinct phenotypical clusters, showing selective impairments of either the Gi or βarrestin/endocytosis pathways with no effect on Gs activation. Robustness of the results is confirmed using simulation-based error propagation. The structural changes resulting from functionally biasing mutations centered around the DRY, NPxxY, and PIF motifs, selectively linking these micro-switches to unique signaling profiles. Our data identify different receptor regions that are important for the stabilization of distinct conformations underlying functional selectivity.

Ligand-induced biased signaling is thought to result in part from ligand-specific receptor conformations that cause the engagement of distinct effectors. Here the authors trace and evaluate the impact of mutations of the β2–adrenergic receptor on multiple signaling outputs to provide structural-level insight into the determinants of GPCR functional selectivity.

Agonist Binding to Chemosensory Receptors: A Systematic Bioinformatics Analysis

Human G-protein coupled receptors (hGPCRs) constitute a large and highly pharmaceutically relevant membrane receptor superfamily. About half of the hGPCRs' family members are chemosensory receptors, involved in bitter taste and olfaction, along with a variety of other physiological processes. Hence these receptors constitute promising targets for pharmaceutical intervention. Molecular modeling has been so far the most important tool to get insights on agonist binding and receptor activation. Here we investigate both aspects by bioinformatics-based predictions across all bitter taste and odorant receptors for which site-directed mutagenesis data are available. First, we observe that state-of-the-art homology modeling combined with previously used docking procedures turned out to reproduce only a limited fraction of ligand/receptor interactions inferred by experiments. This is most probably caused by the low sequence identity with available structural templates, which limits the accuracy of the protein model and in particular of the side-chains' orientations. Methods which transcend the limited sampling of the conformational space of docking may improve the predictions. As an example corroborating this, we review here multi-scale simulations from our lab and show that, for the three complexes studied so far, they significantly enhance the predictive power of the computational approach. Second, our bioinformatics analysis provides support to previous claims that several residues, including those at positions 1.50, 2.50, and 7.52, are involved in receptor activation.

Modeling of Polarization Transfer Kinetics in Protein Hydration Using Hyperpolarized Water

Water-protein interactions play a central role in protein structure, dynamics, and function. These interactions, traditionally, have been studied using nuclear magnetic resonance (NMR) by measuring chemical exchange and nuclear Overhauser effect (NOE). Polarization transferred from hyperpolarized water can result in substantial transient signal enhancements of protein resonances due to these processes. Here, we use dissolution dynamic nuclear polarization and flow-NMR for measuring the pH dependence of transferred signals to the protein trypsin. A maximum enhancement of 20 is visible in the amide proton region of the spectrum at pH 6.0, and of 47 at pH 7.5. The aliphatic region is enhanced up to 2.3 times at pH 6.0 and up to 2.5 times at pH 7.5. The time dependence of these observed signals can be modeled quantitatively using rate equations incorporating chemical exchange to amide sites and, optionally, intramolecular NOE to aliphatic protons. On the basis of these two- and three-site models, average exchange (k_ex) and cross-relaxation rates (σ) obtained were k_ex = 12 s^-1, σ = -0.33 s^-1 for pH 7.5 and k_ex = 1.8 s^-1, σ = -0.72 s^-1 for pH 6.0 at a temperature of 304 K. These values were validated using conventional EXSY and NOESY measurements. In general, a rapid measurement of exchange and cross-relaxation rates may be of interest for the study of structural changes of the protein occurring on the same time scale. Besides protein-water interactions, interactions with cosolvent or solutes can further be investigated using the same methods.

Evolution of nonspectral rhodopsin function at high altitudes

High-altitude environments present a range of biochemical and physiological challenges for organisms through decreases in oxygen, pressure, and temperature relative to lowland habitats. Protein-level adaptations to hypoxic high-altitude conditions have been identified in multiple terrestrial endotherms; however, comparable adaptations in aquatic ectotherms, such as fishes, have not been as extensively characterized. In enzyme proteins, cold adaptation is attained through functional trade-offs between stability and activity, often mediated by substitutions outside the active site. Little is known whether signaling proteins [e.g., G protein-coupled receptors (GPCRs)] exhibit natural variation in response to cold temperatures. Rhodopsin (RH1), the temperature-sensitive visual pigment mediating dim-light vision, offers an opportunity to enhance our understanding of thermal adaptation in a model GPCR. Here, we investigate the evolution of rhodopsin function in an Andean mountain catfish system spanning a range of elevations. Using molecular evolutionary analyses and site-directed mutagenesis experiments, we provide evidence for cold adaptation in RH1. We find that unique amino acid substitutions occur at sites under positive selection in high-altitude catfishes, located at opposite ends of the RH1 intramolecular hydrogen-bonding network. Natural high-altitude variants introduced into these sites via mutagenesis have limited effects on spectral tuning, yet decrease the stability of dark-state and light-activated rhodopsin, accelerating the decay of ligand-bound forms. As found in cold-adapted enzymes, this phenotype likely compensates for a cold-induced decrease in kinetic rates-properties of rhodopsin that mediate rod sensitivity and visual performance. Our results support a role for natural variation in enhancing the performance of GPCRs in response to cold temperatures.

Water permeation through the internal water pathway in activated GPCR rhodopsin

Rhodopsin is a light-driven G-protein-coupled receptor that mediates signal transduction in eyes. Internal water molecules mediate activation of the receptor in a rhodopsin cascade reaction and contribute to conformational stability of the receptor. However, it remains unclear how internal water molecules exchange between the bulk and protein inside, in particular through a putative solvent pore on the cytoplasmic. Using all-atom molecular dynamics simulations, we identified the solvent pore on cytoplasmic side in both the Meta II state and the Opsin. On the other hand, the solvent pore does not exist in the dark-adapted rhodopsin. We revealed two characteristic narrow regions located within the solvent pore in the Meta II state. The narrow regions distinguish bulk and the internal hydration sites, one of which is adjacent to the conserved structural motif "NPxxY". Water molecules in the solvent pore diffuse by pushing or sometimes jumping a preceding water molecule due to the geometry of the solvent pore. These findings revealed a total water flux between the bulk and the protein inside in the Meta II state, and suggested that these pathways provide water molecules to the crucial sites of the activated rhodopsin.

The Roles of Water in the Protein Matrix: A Largely Untapped Resource for Drug Discovery

The value of thoroughly understanding the thermodynamics specific to a drug discovery/design study is well known. Over the past decade, the crucial roles of water molecules in protein structure, function, and dynamics have also become increasingly appreciated. This Perspective explores water in the biological environment by adopting its point of view in such phenomena. The prevailing thermodynamic models of the past, where water was seen largely in terms of an entropic gain after its displacement by a ligand, are now known to be much too simplistic. We adopt a set of terminology that describes water molecules as being "hot" and "cold", which we have defined as being easy and difficult to displace, respectively. The basis of these designations, which involve both enthalpic and entropic water contributions, are explored in several classes of biomolecules and structural motifs. The hallmarks for characterizing water molecules are examined, and computational tools for evaluating water-centric thermodynamics are reviewed. This Perspective's summary features guidelines for exploiting water molecules in drug discovery.

Gonadotropin-Releasing Hormone (GnRH) Receptor Structure and GnRH Binding

Gonadotropin-releasing hormone (GnRH) regulates reproduction. The human GnRH receptor lacks a cytoplasmic carboxy-terminal tail but has amino acid sequence motifs characteristic of rhodopsin-like, class A, G protein-coupled receptors (GPCRs). This review will consider how recent descriptions of X-ray crystallographic structures of GPCRs in inactive and active conformations may contribute to understanding GnRH receptor structure, mechanism of activation and ligand binding. The structures confirmed that ligands bind to variable extracellular surfaces, whereas the seven membrane-spanning α-helices convey the activation signal to the cytoplasmic receptor surface, which binds and activates heterotrimeric G proteins. Forty non-covalent interactions that bridge topologically equivalent residues in different transmembrane (TM) helices are conserved in class A GPCR structures, regardless of activation state. Conformation-independent interhelical contacts account for a conserved receptor protein structure and their importance in the GnRH receptor structure is supported by decreased expression of receptors with mutations of residues in the network. Many of the GnRH receptor mutations associated with congenital hypogonadotropic hypogonadism, including the Glu^2.53(90) Lys mutation, involve amino acids that constitute the conserved network. Half of the ~250 intramolecular interactions in GPCRs differ between inactive and active structures. Conformation-specific interhelical contacts depend on amino acids changing partners during activation. Conserved inactive conformation-specific contacts prevent receptor activation by stabilizing proximity of TM helices 3 and 6 and a closed G protein-binding site. Mutations of GnRH receptor residues involved in these interactions, such as Arg^3.50(139) of the DRY/S motif or Tyr^7.53(323) of the N/DPxxY motif, increase or decrease receptor expression and efficiency of receptor coupling to G protein signaling, consistent with the native residues stabilizing the inactive GnRH receptor structure. Active conformation-specific interhelical contacts stabilize an open G protein-binding site. Progress in defining the GnRH-binding site has recently slowed, with evidence that Tyr^6.58(290) contacts Tyr⁵ of GnRH, whereas other residues affect recognition of Trp³ and Gly¹⁰NH₂. The surprisingly consistent observations that GnRH receptor mutations that disrupt GnRH binding have less effect on "conformationally constrained" GnRH peptides may now be explained by crystal structures of agonist-bound peptide receptors. Analysis of GPCR structures provides insight into GnRH receptor function.

Membrane potentials regulating GPCRs: insights from experiments and molecular dynamics simulations

G-protein coupled receptors (GPCRs) form the largest class of membrane proteins in humans and the targets of most present drugs. Membrane potential is one of the defining characteristics of living cells. Recent work has shown that the membrane voltage, and changes thereof, modulates signal transduction and ligand binding in GPCRs. As it may allow differential signalling patterns depending on tissue, cell type, and the excitation status of excitable cells, GPCR voltage sensitivity could have important implications for their pharmacology. This review summarises recent experimental insights on GPCR voltage regulation and the role of molecular dynamics simulations in identifying the structural basis of GPCR voltage-sensing. We discuss the potential significance for drug design on GPCR targets from excitable and non-excitable cells.

Structural Mechanisms of Voltage Sensing in G Protein-Coupled Receptors

G-protein-coupled receptors (GPCRs) form the largest superfamily of membrane proteins and one-third of all drug targets in humans. A number of recent studies have reported evidence for substantial voltage regulation of GPCRs. However, the structural basis of GPCR voltage sensing has remained enigmatic. Here, we present atomistic simulations on the δ-opioid and M2 muscarinic receptors, which suggest a structural and mechanistic explanation for the observed voltage-induced functional effects. The simulations reveal that the position of an internal Na(+) ion, recently detected to bind to a highly conserved aqueous pocket in receptor crystal structures, strongly responds to voltage changes. The movements give rise to gating charges in excellent agreement with previous experimental recordings. Furthermore, free energy calculations show that these rearrangements of Na(+) can be induced by physiological membrane voltages. Due to its role in receptor function and signal bias, the repositioning of Na(+) has important general implications for signal transduction in GPCRs.

Bitter taste receptors: Novel insights into the biochemistry and pharmacology

Bitter taste receptors (T2Rs) belong to the super family of G protein-coupled receptors (GPCRs). There are 25 T2Rs expressed in humans, and these interact with a large and diverse group of bitter ligands. T2Rs are expressed in many extra-oral tissues and can perform diverse physiological roles. Structure-function studies led to the identification of similarities and dissimilarities between T2Rs and Class A GPCRs including amino acid conservation and novel motifs. However, the efficacy of most of the T2R ligands is not yet elucidated and the biochemical pharmacology of T2Rs is poorly understood. Recent studies on T2Rs characterized novel ligands including blockers for these receptors that include inverse agonist and antagonists. In this review we discuss the techniques used for elucidating bitter blockers, concept of ligand bias, generic amino acid numbering, the role of cholesterol, and conserved water molecules in the biochemistry and pharmacology of T2Rs.

Characterization of a novel water pocket inside the human Cx26 hemichannel structure

Connexins (Cxs) are a family of vertebrate proteins constituents of gap junction channels (GJCs) that connect the cytoplasm of adjacent cells by the end-to-end docking of two Cx hemichannels. The intercellular transfer through GJCs occurs by passive diffusion allowing the exchange of water, ions, and small molecules. Despite the broad interest to understand, at the molecular level, the functional state of Cx-based channels, there are still many unanswered questions regarding structure-function relationships, perm-selectivity, and gating mechanisms. In particular, the ordering, structure, and dynamics of water inside Cx GJCs and hemichannels remains largely unexplored. In this work, we describe the identification and characterization of a believed novel water pocket-termed the IC pocket-located in-between the four transmembrane helices of each human Cx26 (hCx26) monomer at the intracellular (IC) side. Using molecular dynamics (MD) simulations to characterize hCx26 internal water structure and dynamics, six IC pockets were identified per hemichannel. A detailed characterization of the dynamics and ordering of water including conformational variability of residues forming the IC pockets, together with multiple sequence alignments, allowed us to propose a functional role for this cavity. An in vitro assessment of tracer uptake suggests that the IC pocket residue Arg-143 plays an essential role on the modulation of the hCx26 hemichannel permeability.

Functional elements of the gastric inhibitory polypeptide receptor: Comparison between secretin- and rhodopsin-like G protein-coupled receptors

Innovative crystallographic techniques have resulted in an exponential growth in the number of solved G-protein coupled receptor (GPCR) structures and a better understanding of the mechanisms of class A receptor activation and G protein binding. The recent release of the type 1 receptor for the corticotropin-releasing factor and the glucagon receptor structures, two members of the secretin-like family, gives the opportunity to understand these mechanisms of activation in this family of GPCRs. Here, we addressed the comparison of the functional elements of class A and secretin-like GPCRs, using the glucose-dependent insulinotropic polypeptide receptor (GIPR) as a model receptor. Inactive and active models of GIPR permitted to select, by structural homology with class A GPCRs, several residues that may form key interactions presumably involved in receptor activation and Gs coupling, for pharmacological evaluation. Mutants on these amino acids were expressed in HEKT 293 cells and characterized in terms of GIP-induced cAMP production. We identified various functional domains spanning from the peptide-binding to the G protein pockets: including: a network linking the extracellular part of transmembrane (TM) 6 with TMs 2 and 7; a polar lock that resembles the ionic-lock in class A GPCRs; an interaction between TMs 3 and 7 that favors activation; and two clusters of polar/charged and of hydrophobic residues that interact with the C-terminus of the Gα. The results show that despite the low degree of sequence similarity between rhodopsin- and secretin-like GPCRs, the two families share conserved elements in their mechanisms of activation and G protein binding.

Structural dynamics and thermostabilization of neurotensin receptor 1

The neurotensin receptor NTSR1 binds the peptide agonist neurotensin (NTS) and signals preferentially via the Gq protein. Recently, Grisshammer and co-workers reported the crystal structure of a thermostable mutant NTSR1-GW5 with NTS bound. Understanding how the mutations thermostabilize the structure would allow efficient design of thermostable mutant GPCRs for protein purification, and subsequent biophysical studies. Using microsecond scale molecular dynamics simulations (4 μs) of the thermostable mutant NTSR1-GW5 and wild type NTSR1, we have elucidated the structural and energetic factors that affect the thermostability and dynamics of NTSR1. The thermostable mutant NTSR1-GW5 is found to be less flexible and less dynamic than the wild type NTSR1. The point mutations confer thermostability by improving the interhelical hydrogen bonds, hydrophobic packing, and receptor interactions with the lipid bilayer, especially in the intracellular regions. During MD, NTSR1-GW5 becomes more hydrated compared to wild type NTSR1, with tight hydrogen bonded water clusters within the transmembrane core of the receptor, thus providing evidence that water plays an important role in improving helical packing in the thermostable mutant. Our studies provide valuable insights into the stability and functioning of NTSR1 that will be useful in future design of thermostable mutants of other peptide GPCRs.

Predicted structures for kappa opioid G-protein coupled receptor bound to selective agonists

Human kappa opioid receptor (κ-OR), a G protein-coupled receptor (GPCR), has been identified as a drug target for treatment of such human disorders as pain perception, neuroendocrine physiology, affective behavior, and cognition. In order to find more selective and active agonists, one would like to do structure based drug design. Indeed, there is an X-ray structure for an antagonist bound to κ-OR, but structures for activated GPCRs are quite different from those for the inactive GPCRs. Here we predict the ensemble of 24 low-energy structures of human kappa opioid receptor (κ-OR), obtained by application of the GEnSeMBLE (GPCR Ensemble of Structures in Membrane Bilayer Environment) complete sampling method, which evaluates 13 trillion combinations of tilt and rotation angles for κ-OR to select the best 24. To validate these structures, we used the DarwinDock complete sampling method to predict the binding sites for five known agonists (ethylketocyclazocine, bremazocine, pentazocine, nalorphine, and morphine) bound to all 24 κ-OR conformations. We find that some agonists bind selectively to receptor conformations that lack the salt bridge between transmembrane domains 3 and 6 as expected for active conformations. These 3D structures for κ-OR provide a structural basis for understanding ligand binding and activation of κ-OR, which should be useful for guiding subtype specific drug design.

Structure and function of serotonin G protein-coupled receptors

Serotonin receptors are prevalent throughout the nervous system and the periphery, and remain one of the most lucrative and promising drug discovery targets for disorders ranging from migraine headaches to neuropsychiatric disorders such as schizophrenia and depression. There are 14 distinct serotonin receptors, of which 13 are G protein-coupled receptors (GPCRs), which are targets for approximately 40% of the approved medicines. Recent crystallographic and biochemical evidence has provided a converging understanding of the basic structure and functional mechanics of GPCR activation. Currently, two GPCR crystal structures exist for the serotonin family, the 5-HT1B and 5-HT2B receptor, with the antimigraine and valvulopathic drug ergotamine bound. The first serotonin crystal structures not only provide the first evidence of serotonin receptor topography but also provide mechanistic explanations into functional selectivity or biased agonism. This review will detail the findings of these crystal structures from a molecular and mutagenesis perspective for driving rational drug design for novel therapeutics incorporating biased signaling.

Rare platelet GPCR variants: what can we learn?

Platelet-expressed GPCRs are critical regulators of platelet function. Pharmacological blockade of these receptors forms a powerful therapeutic tool in the treatment and prevention of arterial thrombosis associated with coronary atherosclerosis and ischaemic stroke. However, anti-thrombotic drug therapy is associated with high inter-patient variability in therapeutic response and adverse bleeding side effects. In order to optimize the use of existing anti-platelet drugs and to develop new therapies, more detailed knowledge is required relating to the molecular mechanisms that regulate GPCR and therefore platelet function. One approach has been to identify rare, function-disrupting mutations within key platelet proteins in patients with bleeding disorders. In this review, we describe how an integrated functional genomics strategy has contributed important structure-function information about platelet GPCRs with specific emphasis upon purinergic and thromboxane A2 receptors. We also discuss the potential implications these findings have for pharmacotherapy and for understanding the molecular basis of mild bleeding disorders.

Differential light-induced responses in sectorial inherited retinal degeneration

Retinitis pigmentosa (RP) is a group of genetically and clinically heterogeneous inherited degenerative retinopathies caused by abnormalities of photoreceptors or retinal pigment epithelium in the retina leading to progressive sight loss. Rhodopsin is the prototypical G-protein-coupled receptor located in the vertebrate retina and is responsible for dim light vision. Here, novel M39R and N55K variants were identified as causing an intriguing sector phenotype of RP in affected patients, with selective degeneration in the inferior retina. To gain insights into the molecular aspects associated with this sector RP phenotype, whose molecular mechanism remains elusive, the mutations were constructed by site-directed mutagenesis, expressed in heterologous systems, and studied by biochemical, spectroscopic, and functional assays. M39R and N55K opsins had variable degrees of chromophore regeneration when compared with WT opsin but showed no gross structural misfolding or altered trafficking. M39R showed a faster rate for transducin activation than WT rhodopsin with a faster metarhodopsinII decay, whereas N55K presented a reduced activation rate and an altered photobleaching pattern. N55K also showed an altered retinal release from the opsin binding pocket upon light exposure, affecting its optimal functional response. Our data suggest that these sector RP mutations cause different protein phenotypes that may be related to their different clinical progression. Overall, these findings illuminate the molecular mechanisms of sector RP associated with rhodopsin mutations.

Perturbation of fluid dynamics properties of water molecules during G protein-coupled receptor-ligand recognition: the human A2A adenosine receptor as a key study

Recent advances in structural biology revealed that water molecules play a crucial structural role in the protein architecture and ligand binding of G protein-coupled receptors. In this work, we present an alternative approach to monitor the time-dependent organization of water molecules during the final stage of the ligand-receptor recognition process by means of membrane molecular dynamics simulations. We inspect the variation of fluid dynamics properties of water molecules upon ligand binding with the aim to correlate the results with the binding affinities. The outcomes of this analysis are transferred into a bidimensional graph called water fluid dynamics maps, that allow a fast graphical identification of protein "hot-spots" characterized by peculiar shape and electrostatic properties that can play a critical role in ligand binding. We hopefully believe that the proposed approach might represent a valuable tool for structure-based drug discovery that can be extended to cases where crystal structures are not yet available, or have not been solved at high resolution.

Unifying family A GPCR theories of activation

Several new pairs of active and inactive GPCR structures have recently been solved enabling detailed structural insight into the activation process, not only of rhodopsin but now also of the β2 adrenergic, M2 muscarinic and adenosine A2A receptors. Combined with structural analyses they have enabled us to examine the different recent theories proposed for GPCR activation and show that they are all indeed parts of the same process, and are intrinsically related through their effect on the central hydrophobic core of GPCRs. This new unifying general process of activation is consistent with the identification of known constitutively active mutants and an in-depth conservational analysis of significant residues implicated in the process.

Synthesis and structure evaluation of new complex butylarylpiperazin-1-yl derivatives

A series of arylpiperazine derivatives of 1,16-diphenyl-19-azahexacyclo-[14.5.1.0^2,15.0^3,8.0^9,14.0^17,21]docosa-2,3,5,7,8,9,11,13,14-nonaene-18,20,22-trione and 4,10-diphenyl-1H,2H,3H,5H-indeno[1,2-f]isoindole-1,3,5-trione was synthesized. The pharmacological profile of compound 4 at the 5-HT_1A receptor was measured by binding assay. The title compounds were tested in cell-based assay against the human immunodeficiency virus type-1. The X-ray crystallographic studies of derivatives 2, 6, 7, 11, 19, and 20 were presented.

Bridging molecular docking to membrane molecular dynamics to investigate GPCR-ligand recognition: the human A₂A adenosine receptor as a key study

G protein-coupled receptors (GPCRs) represent the largest family of cell-surface receptors and about one-third of the actual targets of clinically used drugs. Following the progress made in the field of GPCRs structural determination, docking-based screening for novel potent and selective ligands is becoming an increasingly adopted strategy in the drug discovery process. However, this methodology is not yet able to anticipate the "bioactive" binding mode and discern it among other conformations. In the present work, we present a novel approach consisting in the integration of molecular docking and membrane MD simulations with the aim to merge the rapid sampling of ligand poses into in the binding site, typical of docking algorithms, with the thermodynamic accuracy of MD simulations in describing, at the molecular level, the stability a GPCR-ligand complex embedded into explicit lipid-water environment. To validate our approach, we have chosen as a key study the human A(2A) adenosine receptor (hA(2A) AR) and selected four receptor-antagonist complexes and one receptor-agonist complex that have been recently crystallized. In light of the obtained results, we believe that our novel strategy can be extended to other GPCRs and might represent a valuable tool to anticipate the "bioactive" conformation of high-affinity ligands.

Constitutively active rhodopsin and retinal disease

Rhodopsin is the light receptor in rod photoreceptor cells of the retina that initiates scotopic vision. In the dark, rhodopsin is bound to the chromophore 11-cis retinal, which locks the receptor in an inactive state. The maintenance of an inactive rhodopsin in the dark is critical for rod photoreceptor cells to remain highly sensitive. Perturbations by mutation or the absence of 11-cis retinal can cause rhodopsin to become constitutively active, which leads to the desensitization of photoreceptor cells and, in some instances, retinal degeneration. Constitutive activity can arise in rhodopsin by various mechanisms and can cause a variety of inherited retinal diseases including Leber congenital amaurosis, congenital night blindness, and retinitis pigmentosa. In this review, the molecular and structural properties of different constitutively active forms of rhodopsin are overviewed, and the possibility that constitutive activity can arise from different active-state conformations is discussed.

The role of a sodium ion binding site in the allosteric modulation of the A(2A) adenosine G protein-coupled receptor

The function of G protein-coupled receptors (GPCRs) can be modulated by a number of endogenous allosteric molecules. In this study, we used molecular dynamics, radioligand binding, and thermostability experiments to elucidate the role of the recently discovered sodium ion binding site in the allosteric modulation of the human A(2A) adenosine receptor, conserved among class A GPCRs. While the binding of antagonists and sodium ions to the receptor was noncompetitive in nature, the binding of agonists and sodium ions appears to require mutually exclusive conformational states of the receptor. Amiloride analogs can also bind to the sodium binding pocket, showing distinct patterns of agonist and antagonist modulation. These findings suggest that physiological concentrations of sodium ions affect functionally relevant conformational states of GPCRs and can help to design novel synthetic allosteric modulators or bitopic ligands exploiting the sodium ion binding pocket.

MP:PD--a data base of internal packing densities, internal packing defects and internal waters of helical membrane proteins

The membrane protein packing database (MP:PD) (http://proteinformatics.charite.de/mppd) is a database of helical membrane proteins featuring internal atomic packing densities, cavities and waters. Membrane proteins are not tightly packed but contain a considerable number of internal cavities that differ in volume, polarity and solvent accessibility as well as in their filling with internal water. Internal cavities are supposed to be regions of high physical compressibility. By serving as mobile hydrogen bonding donors or acceptors, internal waters likely facilitate transition between different functional states. Despite these distinct functional roles, internal cavities of helical membrane proteins are not well characterized, mainly because most internal waters are not resolved by crystal structure analysis. Here we combined various computational biophysical techniques to characterize internal cavities, reassign positions of internal waters and calculate internal packing densities of all available helical membrane protein structures and stored them in MP:PD. The database can be searched using keywords and entries can be downloaded. Each entry can be visualized in Provi, a Jmol-based protein viewer that provides an integrated display of low energy waters alongside membrane planes, internal packing density, hydrophobic cavities and hydrogen bonds.

Minireview: Toward the establishment of a link between melatonin and glucose homeostasis: association of melatonin MT2 receptor variants with type 2 diabetes

The existence of interindividual variations in G protein-coupled receptor sequences has been recognized early on. Recent advances in large-scale exon sequencing techniques are expected to dramatically increase the number of variants identified in G protein-coupled receptors, giving rise to new challenges regarding their functional characterization. The current minireview will illustrate these challenges based on the MTNR1B gene, which encodes the melatonin MT2 receptor, for which exon sequencing revealed 40 rare nonsynonymous variants in the general population and in type 2 diabetes (T2D) cohorts. Functional characterization of these MT2 mutants revealed 14 mutants with loss of Gi protein activation that associate with increased risk of T2D development. This repertoire of disease-associated mutants is a rich source for structure-activity studies and will help to define the still poorly understood role of melatonin in glucose homeostasis and T2D development in humans. Defining the functional defects in carriers of rare MT2 mutations will help to provide personalized therapies to these patients in the future.

Computationally-predicted CB1 cannabinoid receptor mutants show distinct patterns of salt-bridges that correlate with their level of constitutive activity reflected in G protein coupling levels, thermal stability, and ligand binding

The cannabinoid receptor 1 (CB1), a member of the class A G-protein-coupled receptor (GPCR) family, possesses an observable level of constitutive activity. Its activation mechanism, however, has yet to be elucidated. Previously we discovered dramatic changes in CB1 activity due to single mutations; T3.46A, which made the receptor inactive, and T3.46I and L3.43A, which made it essentially fully constitutively active. Our subsequent prediction of the structures of these mutant receptors indicated that these changes in activity are explained in terms of the pattern of salt-bridges in the receptor region involving transmembrane domains 2, 3, 5, and 6. Here we identified key salt-bridges, R2.37 + D6.30 and D2.63 + K3.28, critical for CB1 inactive and active states, respectively, and generated new mutant receptors that we predicted would change CB1 activity by either precluding or promoting these interactions. We find that breaking the R2.37 + D6.30 salt-bridge resulted in substantial increase in G-protein coupling activity and reduced thermal stability relative to the wild-type reflecting the changes in constitutive activity from inactive to active. In contrast, breaking the D2.63 + K3.28 salt-bridge produced the opposite profile suggesting this interaction is critical for the receptor activation. Thus, we demonstrate an excellent correlation with the predicted pattern of key salt-bridges and experimental levels of activity and conformational flexibility. These results are also consistent with the extended ternary complex model with respect to shifts in agonist and inverse agonist affinity and provide a powerful framework for understanding the molecular basis for the multiple stages of CB1 activation and that of other GPCRs in general.

From atomic structures to neuronal functions of g protein-coupled receptors

G protein-coupled receptors (GPCRs) are essential mediators of signal transduction, neurotransmission, ion channel regulation, and other cellular events. GPCRs are activated by diverse stimuli, including light, enzymatic processing of their N-termini, and binding of proteins, peptides, or small molecules such as neurotransmitters. GPCR dysfunction caused by receptor mutations and environmental challenges contributes to many neurological diseases. Moreover, modern genetic technology has helped identify a rich array of mono- and multigenic defects in humans and animal models that connect such receptor dysfunction with disease affecting neuronal function. The visual system is especially suited to investigate GPCR structure and function because advanced imaging techniques permit structural studies of photoreceptor neurons at both macro and molecular levels that, together with biochemical and physiological assessment in animal models, provide a more complete understanding of GPCR signaling.

The dynamic process of β(2)-adrenergic receptor activation

G-protein-coupled receptors (GPCRs) can modulate diverse signaling pathways, often in a ligand-specific manner. The full range of functionally relevant GPCR conformations is poorly understood. Here, we use NMR spectroscopy to characterize the conformational dynamics of the transmembrane core of the β(2)-adrenergic receptor (β(2)AR), a prototypical GPCR. We labeled β(2)AR with (13)CH(3)ε-methionine and obtained HSQC spectra of unliganded receptor as well as receptor bound to an inverse agonist, an agonist, and a G-protein-mimetic nanobody. These studies provide evidence for conformational states not observed in crystal structures, as well as substantial conformational heterogeneity in agonist- and inverse-agonist-bound preparations. They also show that for β(2)AR, unlike rhodopsin, an agonist alone does not stabilize a fully active conformation, suggesting that the conformational link between the agonist-binding pocket and the G-protein-coupling surface is not rigid. The observed heterogeneity may be important for β(2)AR's ability to engage multiple signaling and regulatory proteins.

Platelet dysfunction associated with the novel Trp29Cys thromboxane A₂ receptor variant

BACKGROUND

Genetic variations that affect the structure of the thromboxane A2 receptor (TP receptor) provide insights into the function of this key platelet and vascular receptor, but are very rare in unselected populations.

OBJECTIVES

To determine the functional consequences of the TP receptor Trp29Cys (W29C) substitution.

PATIENTS/METHODS

We performed a detailed phenotypic analysis of an index case (P1) with reduced platelet aggregation and secretion responses to TP receptor pathway activators, and a heterozygous TP receptor W29C substitution. An analysis of the variant W29C TP receptor expressed in heterologous cells was performed.

RESULTS

Total TP receptor expression in platelets from P1 was similar to that of controls, but there was reduced maximum binding and reduced affinity of binding to the TP receptor antagonist [(3) H]SQ29548. HEK293 cells transfected with W29C TP receptor cDNA showed similar total TP receptor expression to wild-type (WT) controls. However, the TP receptor agonist U46619 was less potent at inducing rises in cytosolic free Ca(2+) in HEK293 cells expressing the W29C TP receptor than in WT controls, indicating reduced receptor function. Immunofluorescence microscopy and cell surface ELISA showed intracellular retention and reduced cell surface expression of the W29C TP receptor in HEK293 cells. Consistent with the platelet phenotype, both maximum binding and the affinity of binding of [(3) H]SQ29548 to the W29C TP receptor were reduced compared to WT controls.

CONCLUSION

These findings extend the phenotypic description of the very rare disorder TP receptor deficiency, and show that the W29C substitution reduces TP receptor function by reducing surface receptor expression and by disrupting ligand binding.

Characterization of the dynamic events of GPCRs by automated computational simulations

The recent advances in membrane protein crystallography have provided extremely valuable structural information of the superfamily of GPCRs (G-protein-coupled receptors). This has been particularly true for a few receptors whose structure was solved several times under different biochemical conditions. It follows that the mechanisms of receptor conformational equilibrium and related dynamic events can be explored by computational simulations. In the present article, we summarize our recent understanding of several dynamic features of GPCRs, accomplished through the use of MD (molecular dynamics) simulations. Our pipeline for the MD simulations of GPCRs, implemented in the web service http://gpcr.usc.es, is updated in the present paper and illustrated by recent applications. Special emphasis is put on the A2A adenosine receptor, one of the selected cases where crystal structures in several conformations and conditions exist, and on the dimerization process of the CXCR4 (CXC chemokine receptor 4).

Quantification of structural distortions in the transmembrane helices of GPCRs

A substantial part of the structural and much of the functional information about G protein-coupled receptors (GPCRs) comes from studies on rhodopsin. Thus, analysis tools for detailed structure comparison are key to see to what extent this information can be extended to other GPCRs. Among the methods to evaluate protein structures and, in particular, helix distortions, HELANAL has the advantage that it provides data (local bend and twist angles) that can be easily translated to structural effects, as a local opening/tightening of the helix.In this work I show how HELANAL can be used to extract detailed structural information of the transmembrane bundle of GPCRs, and I provide some examples on how these data can be interpreted to study basic principles of protein structure, to compare homologous proteins and to study mechanisms of receptor activation. Also, I show how in combination with the sequence analysis tools provided by the program GMoS, distortions in individual receptors can be put in the context of the whole Class A GPCR family. Specifically, quantification of the strong proline-induced distortions in the transmembrane bundle of rhodopsin shows that they are not standard proline kinks. Moreover, the helix distortions in transmembrane helix (TMH) 5 and TMH 6 of rhodopsin are also present in the rest of GPCR crystal structures obtained so far, and thus, rhodopsin-based homology models have modeled correctly these strongly distorted helices. While in some cases the inherent "rhodopsin bias" of many of the GPCR models to date has not been a disadvantage, the availability of more templates will clearly result in better homology models. This type of analysis can be, of course, applied to any protein, and it may be particularly useful for the structural analysis of other membrane proteins. A detailed knowledge of the local structural changes related to ligand binding and how they are translated into larger-scale movements of transmembrane domains is key to understand receptor activation.