Skip to main content
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian
  • Log in
  • My Cart

Main menu

  • Home
  • Articles
    • Current
    • Latest Articles
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • Archive
  • Front Matter
  • News
    • For the Press
    • Highlights from Latest Articles
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Purpose and Scope
    • Editorial and Journal Policies
    • Submission Procedures
    • For Reviewers
    • Author FAQ
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian

User menu

  • Log in
  • My Cart

Search

  • Advanced search
Home
Home

Advanced Search

  • Home
  • Articles
    • Current
    • Latest Articles
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • Archive
  • Front Matter
  • News
    • For the Press
    • Highlights from Latest Articles
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Purpose and Scope
    • Editorial and Journal Policies
    • Submission Procedures
    • For Reviewers
    • Author FAQ

New Research In

Physical Sciences

Featured Portals

  • Physics
  • Chemistry
  • Sustainability Science

Articles by Topic

  • Applied Mathematics
  • Applied Physical Sciences
  • Astronomy
  • Computer Sciences
  • Earth, Atmospheric, and Planetary Sciences
  • Engineering
  • Environmental Sciences
  • Mathematics
  • Statistics

Social Sciences

Featured Portals

  • Anthropology
  • Sustainability Science

Articles by Topic

  • Economic Sciences
  • Environmental Sciences
  • Political Sciences
  • Psychological and Cognitive Sciences
  • Social Sciences

Biological Sciences

Featured Portals

  • Sustainability Science

Articles by Topic

  • Agricultural Sciences
  • Anthropology
  • Applied Biological Sciences
  • Biochemistry
  • Biophysics and Computational Biology
  • Cell Biology
  • Developmental Biology
  • Ecology
  • Environmental Sciences
  • Evolution
  • Genetics
  • Immunology and Inflammation
  • Medical Sciences
  • Microbiology
  • Neuroscience
  • Pharmacology
  • Physiology
  • Plant Biology
  • Population Biology
  • Psychological and Cognitive Sciences
  • Sustainability Science
  • Systems Biology

Diverse GPCRs exhibit conserved water networks for stabilization and activation

A. J. Venkatakrishnan, Anthony K. Ma, Rasmus Fonseca, Naomi R. Latorraca, Brendan Kelly, Robin M. Betz, Chaitanya Asawa, Brian K. Kobilka, and Ron O. Dror
PNAS February 19, 2019 116 (8) 3288-3293; published ahead of print February 19, 2019 https://doi.org/10.1073/pnas.1809251116
A. J. Venkatakrishnan
aDepartment of Computer Science, Stanford University, Stanford, CA 94305;bInstitute for Computational and Mathematical Engineering, Stanford University, Stanford, CA 94305;cDepartment of Molecular and Cellular Physiology, Stanford University School of Medicine, Stanford, CA 94305;dDepartment of Structural Biology, Stanford University School of Medicine, Stanford, CA 94305;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Anthony K. Ma
aDepartment of Computer Science, Stanford University, Stanford, CA 94305;bInstitute for Computational and Mathematical Engineering, Stanford University, Stanford, CA 94305;cDepartment of Molecular and Cellular Physiology, Stanford University School of Medicine, Stanford, CA 94305;dDepartment of Structural Biology, Stanford University School of Medicine, Stanford, CA 94305;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Rasmus Fonseca
cDepartment of Molecular and Cellular Physiology, Stanford University School of Medicine, Stanford, CA 94305;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Naomi R. Latorraca
aDepartment of Computer Science, Stanford University, Stanford, CA 94305;bInstitute for Computational and Mathematical Engineering, Stanford University, Stanford, CA 94305;cDepartment of Molecular and Cellular Physiology, Stanford University School of Medicine, Stanford, CA 94305;dDepartment of Structural Biology, Stanford University School of Medicine, Stanford, CA 94305;eBiophysics Program, Stanford University, Stanford, CA 94305
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Brendan Kelly
aDepartment of Computer Science, Stanford University, Stanford, CA 94305;bInstitute for Computational and Mathematical Engineering, Stanford University, Stanford, CA 94305;cDepartment of Molecular and Cellular Physiology, Stanford University School of Medicine, Stanford, CA 94305;dDepartment of Structural Biology, Stanford University School of Medicine, Stanford, CA 94305;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Robin M. Betz
aDepartment of Computer Science, Stanford University, Stanford, CA 94305;bInstitute for Computational and Mathematical Engineering, Stanford University, Stanford, CA 94305;cDepartment of Molecular and Cellular Physiology, Stanford University School of Medicine, Stanford, CA 94305;dDepartment of Structural Biology, Stanford University School of Medicine, Stanford, CA 94305;eBiophysics Program, Stanford University, Stanford, CA 94305
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Chaitanya Asawa
aDepartment of Computer Science, Stanford University, Stanford, CA 94305;bInstitute for Computational and Mathematical Engineering, Stanford University, Stanford, CA 94305;cDepartment of Molecular and Cellular Physiology, Stanford University School of Medicine, Stanford, CA 94305;dDepartment of Structural Biology, Stanford University School of Medicine, Stanford, CA 94305;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Brian K. Kobilka
cDepartment of Molecular and Cellular Physiology, Stanford University School of Medicine, Stanford, CA 94305;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: kobilka@stanford.eduron.dror@stanford.edu
Ron O. Dror
aDepartment of Computer Science, Stanford University, Stanford, CA 94305;bInstitute for Computational and Mathematical Engineering, Stanford University, Stanford, CA 94305;cDepartment of Molecular and Cellular Physiology, Stanford University School of Medicine, Stanford, CA 94305;dDepartment of Structural Biology, Stanford University School of Medicine, Stanford, CA 94305;eBiophysics Program, Stanford University, Stanford, CA 94305
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Ron O. Dror
  • For correspondence: kobilka@stanford.eduron.dror@stanford.edu
  1. Contributed by Brian K. Kobilka, December 18, 2018 (sent for review June 5, 2018; reviewed by David E. Gloriam and Leonardo Pardo)

  • Article
  • Figures & SI
  • Info & Metrics
  • PDF
Loading

Significance

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.

Abstract

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.

  • GPCR dynamics
  • water molecules
  • activation
  • molecular dynamics simulation
  • polar network

G protein–coupled receptors (GPCRs) are membrane proteins that act as a control panel for our cells: they recognize the presence of extracellular molecules such as hormones or neurotransmitters and stimulate intracellular signaling pathways in response. Humans have over 800 GPCRs, which share a common architecture comprising seven transmembrane (TM) helices (1).

GPCRs transmit signals across the cell membrane by transitioning from an inactive conformational state to an active conformational state. The active conformational state, which is favored by binding of certain ligands to the extracellular side of the GPCR, couples to and stimulates G proteins. Available structures indicate that diverse GPCRs undergo similar conformational changes upon activation (2, 3). However, these GPCRs have evolved to recognize dramatically different ligands: human GPCRs recognize molecules ranging from amines, nucleosides, lipids, and steroids to entire proteins, in most cases with high specificity. This ability to stimulate cellular pathways in response to the presence of specific extracellular ligands makes GPCRs excellent drug targets, and in fact roughly one-third of all drugs act by binding to GPCRs (4). The process by which diverse ligands bind to GPCRs and favor activation is of great interest in molecular biology and drug discovery but remains incompletely understood.

Recent high-resolution crystal structures reveal dozens of water molecules buried in the transmembrane region of GPCRs (5⇓⇓⇓⇓–10). Water has been shown to play a critical role both in the binding of drugs to their targets (11) and in the function of many proteins (12, 13), raising the question of how the water molecules in GPCRs behave and what function they might serve. Indeed, this question was raised even before the high-resolution structures were available (14, 15).

Unfortunately, answering this question from crystal structures alone is difficult, because each crystal structure represents a static snapshot. In reality, water molecules may be highly mobile, and GPCRs themselves are constantly changing conformation (16). Moreover, most of the available GPCR crystal structures were solved at cryogenic temperatures at which water molecules are expected to be much less mobile than at body temperature. Even recovering average positions of water molecules from crystallographic data can be difficult when the corresponding electron density is weak.

Here, we investigate the structural dynamics of water molecules in a diverse set of GPCRs using molecular dynamics (MD) simulations, which allow us to follow the trajectory of each individual water molecule at physiological temperature and in a native-like membrane environment. We compare water dynamics not only between different GPCRs but also between the active and inactive conformational states of the same GPCR. Our results show that the water molecules observed in GPCR crystal structures are not all equal: a few are stable in their crystallographically observed position, but most are highly mobile, spending only a small fraction of their time near that position.

Remarkably, stable water molecules near the G protein-binding site form a network of polar interactions (hydrogen bonds) that rearranges upon receptor activation. This network was almost perfectly conserved across the GPCRs that we examined, in sharp contrast to water-mediated interactions with ligands, which vary widely from one GPCR to another. This conserved water network appears to play a key role in the GPCR activation mechanism that was preserved as GPCRs evolved to recognize different ligands.

Results

We focused on class A GPCRs, as they represent the largest and most diverse class of GPCRs. We chose functionally distinct and evolutionarily distant ligand-activated class A GPCRs for which both inactive-state and active-state structures are available: β2 adrenergic receptor (β2AR) (17, 18), M2 muscarinic receptor (M2R) (19, 20), μ-opioid receptor (MOR) (5, 21), and adenosine A2A receptor (A2AR) (6, 22). We also included δ-opioid receptor (DOR) as it has a particularly high-resolution inactive-state crystal structure (1.8 Å) (7). Collectively, these GPCRs bind to different types of ligands (amines, peptides, nucleosides) and different G proteins (Gs, Gi). All simulations included an explicitly represented lipid bilayer, water, salt ions, and the cocrystallized ligand (SI Appendix, Table S1).

Only a Minority of Crystallographic Water Molecules Are Stable.

In simulations at physiological temperature, only a few of the waters resolved in high-resolution GPCR crystal structures remain near their crystallographically observed positions (Fig. 1A and SI Appendix, Fig. S1). The position of a given molecule was tracked using its 3D coordinates after aligning the simulation frames on the residues proximal to the water molecule (SI Appendix, Materials and Methods). Our simulations are consistent with crystallographic data in that the same regions of the receptor are populated by waters in both crystal structures and simulations (SI Appendix, Fig. S1), but in simulation, most of the waters are highly mobile, showing little preference for the specific positions occupied by waters in the crystal structures. A few crystallographic water positions, however, are consistently occupied by water molecules in simulation. We refer to these crystallographic water molecules as “stable,” even though an individual water molecule at such a position is generally replaced by a different water molecule at the same position at least every few hundred nanoseconds. In DOR simulations, for example, the water in the “kink” of transmembrane helix 6 (TM6) was highly stable with a water molecule observed within 1 Å of the crystallographic position for 99% of the time. On the other hand, the waters within the transmembrane helical bundle were highly mobile; for some of these, a water molecule was found within 1 Å of the crystallographic position less than 10% of the time.

Fig. 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 1.

Stability of water molecules and water-mediated interactions in GPCRs. (A) Positions of water molecules (red spheres) in one region of a GPCR [the δ-opioid receptor (DOR)] in three simulation snapshots. Receptor transmembrane helices are shown as orange ribbons. Most waters are mobile, but the water at the TM6 kink (circled) is not. (B, Left) Water-mediated interactions present in the crystal structure of DOR. (B, Right) Water-mediated interactions formed over 60% of the time in simulation. Hydrogen bonds forming water-mediated interactions are shown as black dashed lines, and the residues involved in these interactions are shown with orange sticks.

To further characterize the behavior of water molecules in the structural framework of the receptor, we investigated the polar interactions (hydrogen bonds) in the receptor formed by the water molecules (Fig. 1B). We focused on water-mediated interactions in which a pair of receptor residues is connected by one or two water molecules in the network of hydrogen bonds. We estimated the stability of a given water-mediated interaction based on the frequency of the interaction’s presence in simulation (SI Appendix, Materials and Methods).

Only a fraction of the crystallographically observed water-mediated interactions are stable in simulation because a majority of the water molecules are highly mobile, forming many transient, alternative polar interactions. For example, of the 42 water-mediated interactions in the DOR crystal structure, only 19 are formed in 60% or more of simulation frames (Fig. 1B and SI Appendix, Fig. S2).

Conserved Networks of State-Independent and State-Dependent Water-Mediated Interactions in GPCR Activation.

Might water molecules play a role in GPCR activation? To address this question, we first compared water-mediated interactions in simulations of inactive and active GPCR states. We also compared the water-mediated interactions across different GPCRs (β2AR, A2AR, MOR, and M2R). For the comparisons, we analyzed interactions between residues in structurally equivalent positions at different receptors. Structurally equivalent residues were identified using GPCRdb numbers (23), which are receptor-independent generic numbers for referencing structurally equivalent positions in GPCRs (SI Appendix, Materials and Methods).

We found that a surprisingly large set of stable water-mediated interactions (i.e., those formed over 60% of the time in simulation) was shared across all receptors examined in a given activation state (Fig. 2 and SI Appendix, Fig. S3): six interactions mediated generally by four waters in the inactive state and 12 interactions mediated generally by five waters in the active state.

Fig. 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 2.

Conserved water-mediated interactions in inactive and active states of GPCRs. Comparison of the stability of structurally equivalent water-mediated interactions across different GPCRs in (A) inactive states and (B) active states. Four functionally distinct and evolutionarily distant class A GPCRs were compared: A2AR, β2AR, M2R, and MOR. The stability of the water-mediated interactions is estimated based on the frequency of formation in simulations. (A and B, Left) Frames from simulations of four functionally distinct GPCRs. Receptors are shown as ribbons and water molecules are shown as spheres. (A and B, Middle) Heatmap comparing the stability of water-mediated interactions. Columns indicate receptors and rows indicate pairs of structurally equivalent positions, shown using GPCRdb numbers. The cells indicate the measure of stability of water-mediated interactions. Stability is shown as percentage values using corresponding colors. Lighter shades indicate water-mediated interactions that have low stability and darker shades indicate water-mediated interactions that have high stability. Rows with dark cells in all columns indicate interactions that have high stability in all of the receptors being compared. (A and B, Right) Water-mediated interactions that are stable (frequency >60%) across all four receptors mapped onto a GPCR structure. Hydrogen bonds forming water-mediated interactions are shown as black dashed lines, with the residues involved shown as sticks. Inactive state is shown in orange and active state is shown in green.

These conserved interactions represent 20–30% of the stable water-mediated interactions observed in each receptor individually. In contrast, a majority of the remaining stable interactions are found in just one of A2AR, β2AR, M2R, and MOR. For example, of 32 stable water-mediated interactions in active-state M2R, 12 are shared across all four active-state receptors, and 14 of the remaining 20 are unique to M2R.

Among the conserved, stable, water-mediated interactions, three interactions are maintained in a state-independent manner; i.e., they are present in both the inactive states and the active states of A2AR, β2AR, M2R, and MOR. While one of these interactions stabilizes a helical kink of TM6, the other two water molecules are proximal to the G protein-coupling region and mediate interactions between residues on TM1 and TM2 and between residues on TM7 (Fig. 3A).

Fig. 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 3.

State-independent water-mediated interaction network: conserved and stable water-mediated interactions maintained across the inactive and the active states of diverse GPCRs. (A) Stable water-mediated interactions (frequency >60%) present in both the inactive (orange) and the active (green) states of diverse GPCRs. (A, Top) The conserved stable water-mediated interactions mapped onto either an inactive (orange) or active (green) structure. Residues are shown as sticks, water molecules are shown as spheres, and hydrogen bonds forming water-mediated interactions are shown as black dashed lines. (A, Bottom) The conserved stable water-mediated interactions are shown using “flareplots.” In flareplots, the amino acid residues in the transmembrane helices and helix 8 are shown as points on the circle. The water-mediated interactions between the residues are shown as chords connecting the points on the circle. (B) Water-mediated interactions common across the high-resolution structures (resolution of 2.1 Å or better) of diverse GPCRs: inactive Orexin-2 (PDB ID code 5WQC), D4 dopamine (PDB ID code 5WIU), δ-opioid (PDB ID code 4N6H), A2A adenosine (PDB ID code 5IU4), and β1 adrenergic (PDB ID code 4BVN) receptors, and active μ-opioid receptor (PDB ID: 5C1M). Residues are shown as sticks, water molecules are shown as spheres, and hydrogen bonds forming water-mediated interactions are shown as black dashed lines.

Two additional surprising observations link the conserved, stable, water-mediated interactions to the GPCR activation process. First, the remaining conserved, stable water-mediated interactions are in the lower half of the GPCR, near the G protein-binding interface, where the largest structural rearrangements occur upon activation. Second, although these interactions are conserved across diverse GPCRs in a given activation state, they differ dramatically between the active and inactive states. Every conserved stable water molecule near the G protein-binding interface changes interaction partners upon activation, and a number of water molecules form conserved, stable water-mediated interactions in the active state but not in the inactive state. Water-mediated interactions between TM1, TM2, and TM7 are stably formed in the inactive states (Fig. 4A). This network is disrupted upon activation and a new stable network of water-mediated interactions is formed between TM3, TM5, TM6, and TM7 (Fig. 4A). Tyr7x53 (GPCRdb number) of the NPxxY motif acts as a switch between the inactive and active water-mediated networks. In the inactive state, Tyr7x53 interacts with a water molecule at the TM1–TM2 interface, and upon activation, Tyr7x53 switches conformation and interacts with water molecules at the interface of TM3, TM5, and TM6. Tyr7x53 has previously been shown to be an activation-associated microswitch in a water-independent role as well (2, 24).

Fig. 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 4.

State-dependent water-mediated interaction network: conserved and stable water-mediated interactions maintained exclusively in the inactive states or the active states of diverse GPCRs. (A) Stable water-mediated interactions (frequency >60%) present exclusively in all of the inactive (orange) or all of the active (green) states of diverse GPCRs. (A, Top) The conserved stable water-mediated interactions of inactive states and active states mapped onto inactive (orange) and active (green) structures, respectively. Residues are shown as sticks, water molecules are shown as spheres, and hydrogen bonds forming water-mediated interactions are shown as black dashed lines. (A, Bottom) The conserved stable water-mediated interactions are shown using flareplots. (B) Water-mediated interactions present exclusively in the crystal structures of GPCRs in inactive state or active state. The high-resolution inactive state structures are of the following GPCRs: Orexin-2 (PDB ID code 5WQC), D4 dopamine (PDB ID code 5WIU), δ-opioid (PDB ID code 4N6H), and A2A adenosine (PDB ID code 5IU4) receptors. The active-state crystal structures of the following GPCRs are considered: rhodopsin (PDB ID code 2X72), viral chemokine receptor US28 (PDB ID code 4XT1), and μ-opioid receptor (PDB ID code 5C1M). Residues are shown as sticks, water molecules are shown as spheres, and hydrogen bonds forming water-mediated interactions are shown as black dashed lines.

All of the side chains involved in this rewiring of interactions with stable waters upon activation are highly conserved as polar residues across class A GPCRs (>89%) (SI Appendix, Table S2). Taken together, our observations indicate that stable waters near the G protein-coupling interface form a conserved polar interaction network that rearranges upon receptor activation. The fact that this water-mediated activation switch was preserved as the GPCR family evolved suggests that it is functionally important.

Comparative Analysis of Crystal Structures Confirms Conservation of State-Dependent and State-Independent Water-Mediated Interactions.

We postulated that the conserved stable water-mediated interactions identified from the MD simulations would correspond to water-mediated interactions that are maintained across crystal structures of different receptors. To investigate this, we compared the water-mediated interactions from all available high-resolution crystal structures (SI Appendix, Materials and Methods): inactive Orexin-2, D4 dopamine, δ-opioid, A2A adenosine, and β1 adrenergic receptors and active μ-opioid receptor (Fig. 3B). We compared the water-mediated interactions across structurally equivalent positions (SI Appendix, Materials and Methods). A network of seven water-mediated interactions involving three water molecules is maintained in all these crystal structures (Fig. 3B). Interestingly, this network of water-mediated interactions is highly similar to the network of conserved, state-independent, stable interactions observed in the MD simulations. All of the three water-mediated interactions found to be conserved and stable across the simulations (Fig. 3A) are maintained consistently across the crystal structures (Fig. 3B). The remaining interactions that are consistently maintained in the crystal structures involve the same three waters as the conserved stable water-mediated interactions in the simulations.

Similarly, we compared the water-mediated interactions across high-resolution crystal structures of inactive states and crystal structures of active states. We found a network of two water-mediated interactions formed by Tyr7x53 to be maintained exclusively in the inactive state across Orexin-2, D4 dopamine, δ-opioid, and adenosine A2A receptors (Fig. 4B). The only exception in which this network is absent is in the case of the inactive β1 adrenergic receptor, and this is explained by the presence of a Tyr7x53Leu thermostabilizing mutation, which would prevent formation of any polar interactions. Interestingly, the conserved inactive state-dependent network identified from crystal structures is identical to the conserved, stable water-mediated interaction network found in the MD simulations of inactive states (Fig. 4 A and B). Similarly, we also found a network of three water-mediated interactions formed by Tyr7x53 to be maintained exclusively in the crystal structures of active-state rhodopsin, chemokine receptor US28, and μ-opioid receptor (Fig. 4B). All of the three water-mediated interactions identified from the crystal structures of active states are present in the network of conserved, stable water-mediated interactions found in the MD simulations of the active states (Fig. 4 A and B).

Taken together, the comparative analysis of crystal structures supports the presence of conserved state-dependent and state-independent water-mediated interactions in GPCRs.

Stable Water Molecules in the Ligand-Binding Pocket Vary Between Different GPCRs.

A number of the water molecules in high-resolution crystal structures of GPCRs form direct interactions with the ligand, and previous studies have pointed out the importance of water in determining ligand affinities and binding rates (25). We thus also examined the occupancy of waters in the binding pocket and the stability of water-mediated interactions between ligands and receptor residues in simulation.

Most waters in the binding pocket are highly mobile. Certain waters and water-mediated interactions, however, are stable. For example, in simulations of inactive-state DOR, two waters in the binding pocket are highly stable near their crystallographic positions (SI Appendix, Fig. S4A). These two stable waters mediate interactions between a hydroxyl group on the ligand and polar residues on TM5 and TM6.

No stable water-mediated interactions between the ligand and the receptor appear to be conserved across diverse GPCRs. Our simulations showed no such interactions conserved across all of the receptors examined, and only very limited structurally equivalent interactions between subsets of these receptors (SI Appendix, Fig. S4B). This variability in water-mediated interactions across ligand-binding pockets suggests that, unlike the G protein-coupling region, the binding pocket has been adapted to exploit different water-mediated interactions in different GPCRs as the GPCR family evolved.

Discussion

The human genome includes hundreds of class A GPCRs, and these receptors have evolved to recognize extremely diverse extracellular ligands—ranging from nucleosides and amines to peptides and even entire proteins—often with very high specificity. On the other hand, all of these GPCRs couple to the same set of intracellular G proteins, which share very similar structures. One might thus imagine that GPCRs would have evolved a common activation mechanism. Indeed, previous studies have identified conserved networks of amino acid contacts and ionic interactions that appear to facilitate such a mechanism (2, 26, 27).

Our results suggest that certain water molecules also play an important, conserved role in GPCR activation. In particular, we discovered a network of highly stable water molecules near the G protein-binding site that is conserved almost perfectly across diverse class A GPCRs and that rearranges in a conserved manner upon GPCR activation. Waters also play a critical role in mediating ligand binding in GPCRs, but the water-mediated interactions in the ligand-binding pocket vary widely among class A GPCRs, likely reflecting the fact that the binding pockets of these different GPCRs have evolved to recognize different ligands.

Our results complement those of previous experimental and computational studies. Recent high-resolution crystal structures have revealed the presence of dozens of waters within GPCRs. Analyses of earlier, lower-resolution crystal structures have suggested the presence of certain conserved structural water molecules (14, 15). Radiolytic foot-printing has provided insights into tightly bound waters and dynamic waters in rhodopsin (28). Simulation-based studies have shown the presence of activation-associated internal water pathways in GPCRs (28, 29) and indicated a role for waters in ligand binding (30⇓–32).

A few caveats are in order. First, the force fields used for MD simulation are not perfect, potentially introducing artifacts into our results; simulations, however, have proven useful in elucidating water dynamics that are difficult to characterize experimentally (33, 34). Second, the strengths of hydrogen bonds in proteins can vary depending on factors such as the local structural environment (35), but our current analysis does not take this into account. Third, we lack a direct method to determine the functional importance of a water molecule in a particular position. We thus use stability and conservation across receptors as rough proxies for importance. Fourth, our analysis of stable water molecules and water-mediated interactions does not cover all of the ways in which water may contribute to GPCR function; for example, displacement of certain water molecules from the binding pocket upon ligand binding is a determinant of ligand-binding affinity (36). Finally, our analysis has focused on class A GPCRs as they are the largest and most diverse class. Future work will be necessary to understand the role of water molecules in the activation of other GPCR classes.

Our findings have several implications for GPCR drug discovery and GPCR engineering. First, understanding the behavior of water molecules in the binding pocket can aid structure-based drug discovery (25). Our results show that crystallographically resolved waters cannot all be treated as equal. One needs to characterize the dynamics of water molecules (e.g., through simulations) to identify the stable waters, which can then be exploited for drug design (37).

Second, the conserved stable waters and water-mediated interactions identified in this study can aid in structural modeling. Homology models of GPCRs are widely used for drug discovery, both when no crystal structure of the target is available and when the available structures are not in the desired activation state. The conserved water molecules identified in this study for active-state and inactive-state class A GPCRs can be incorporated into homology models and used to refine those models (see templates provided in Datasets S1 and S2). In both the inactive-state and the active-state models, the conserved water molecules at the TM6 helical kink, between positions 1x50 and 2x50 and between positions 7x45 and 7x49, should typically be included (Fig. 3). In addition, in the inactive-state models, in the water-mediated interaction of Tyr7x53 with positions 1x50 and 2x50, an additional water molecule may be included unless Tyr7x53 can directly interact with the water molecule interacting with 1x50 and 2x50 (Fig. 4). On the other hand, in the active-state models, the water molecules mediating interactions between positions 3x43, 5x58, and 7x53 and between 6x40 and 7x53 should typically be included (Fig. 4). An additional water molecule may be included between 6x40 and 7x53, as seen in the case of US28 (Fig. 4). While including these waters in the model, one needs to ensure that the chemical nature of the interacting amino acids and their side-chain geometry allow for the formation of a hydrogen bond with the water. Similar to inclusion of conserved water molecules in homology models, these water molecules can also be incorporated into low-resolution experimental GPCR structures, such as cryo-EM structures (38), in which few waters are resolved.

Third, mutating GPCR residues to stabilize either the active state or the inactive state has proven useful in GPCR crystallography as well as in functional investigations. The distinct water-mediated interaction networks associated with different activation states provide new avenues for biochemically perturbing the conformational landscape of GPCRs. In particular, residues interacting with the conserved, stable waters identified in this study are excellent candidate sites for introducing mutations to shift the conformational ensemble toward the active or inactive state.

More generally, our results reveal a “hidden” structural feature of GPCRs, a network of stable water-mediated interactions that is not evident from crystal structures alone. Our study also highlights the untapped potential of comparing structural dynamics across a drug target family to identify shared and distinct molecular features.

Materials and Methods

The MD simulation trajectories of inactive β2AR, M2R, and MOR and active β2AR were obtained from previously published studies (SI Appendix, Table S1). New MD simulations were performed for inactive DOR and A2AR and active A2AR, M2R, and MOR (SI Appendix, Materials and Methods and Table S1). The simulations were all performed using CHARMM force fields (SI Appendix, Table S1). We verified that the simulations initiated from the inactive-state and active-state structures were maintained in their respective states. Water-mediated interactions were computed for every frame of the MD simulations based on geometric criteria (SI Appendix, Materials and Methods). The stability of a water-mediated interaction between a pair of residues was defined as the fraction of frames in which either a direct or extended water-mediated interaction is formed. The occupancy of water molecules near crystallographic waters was calculated as the fraction of simulation frames in which a water molecule is within 1 Å of the crystallographic water position (SI Appendix, Materials and Methods).

Acknowledgments

This work was supported by National Institutes of Health (NIH) Grant R01GM127359 (to R.O.D.); NIH Grants R37DA036246 and R01GM083118 (to B.K.K.); a Stanford ChEM-H seed grant (to A.J.V.); Grant NNF15OC0015268 from Novo Nordisk Foundation and Stanford Bio-X programme (to R.F.); and National Science Foundation Graduate Research Fellowships (to N.R.L. and R.M.B.). B.K.K. is a Chan–Zuckerberg Biohub Investigator.

Footnotes

  • ↵1A.J.V. and A.K.M. contributed equally to this work.

  • ↵2To whom correspondence may be addressed. Email: kobilka{at}stanford.edu or ron.dror{at}stanford.edu.
  • Author contributions: A.J.V. and R.O.D. designed research; A.J.V., A.K.M., R.F., N.R.L., B.K., R.M.B., C.A., and R.O.D. performed research; A.J.V. and R.F. contributed new reagents/analytic tools; A.J.V., A.K.M., B.K.K., and R.O.D. analyzed data; and A.J.V., A.K.M., B.K.K., and R.O.D. wrote the paper.

  • Reviewers: D.E.G., University of Copenhagen; and L.P., Universitat Autonoma de Barcelona.

  • Conflict of interest statement: B.K.K. is a cofounder of and consultant for ConfometRx, Inc.

  • This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1809251116/-/DCSupplemental.

  • Copyright © 2019 the Author(s). Published by PNAS.

This open access article is distributed under Creative Commons Attribution-NonCommercial-NoDerivatives License 4.0 (CC BY-NC-ND).

References

  1. ↵
    1. Fredriksson R,
    2. Lagerström MC,
    3. Lundin L-G,
    4. Schiöth HB
    (2003) The G-protein-coupled receptors in the human genome form five main families. Phylogenetic analysis, paralogon groups, and fingerprints. Mol Pharmacol 63:1256–1272.
    .
    OpenUrlAbstract/FREE Full Text
  2. ↵
    1. Venkatakrishnan AJ, et al.
    (2016) Diverse activation pathways in class A GPCRs converge near the G-protein-coupling region. Nature 536:484–487.
    .
    OpenUrlCrossRef
  3. ↵
    1. Venkatakrishnan AJ, et al.
    (2013) Molecular signatures of G-protein-coupled receptors. Nature 494:185–194.
    .
    OpenUrlCrossRefPubMed
  4. ↵
    1. Santos R, et al.
    (2017) A comprehensive map of molecular drug targets. Nat Rev Drug Discov 16:19–34.
    .
    OpenUrlCrossRefPubMed
  5. ↵
    1. Huang W, et al.
    (2015) Structural insights into µ-opioid receptor activation. Nature 524:315–321.
    .
    OpenUrlCrossRefPubMed
  6. ↵
    1. Segala E, et al.
    (2016) Controlling the dissociation of ligands from the adenosine A2A receptor through modulation of salt bridge strength. J Med Chem 59:6470–6479.
    .
    OpenUrl
  7. ↵
    1. Fenalti G, et al.
    (2014) Molecular control of δ-opioid receptor signalling. Nature 506:191–196.
    .
    OpenUrlCrossRefPubMed
  8. ↵
    1. Suno R, et al.
    (2018) Crystal structures of human orexin 2 receptor bound to the subtype-selective antagonist EMPA. Structure 26:7–19.e5.
    .
    OpenUrl
  9. ↵
    1. Wang S, et al.
    (2017) D4 dopamine receptor high-resolution structures enable the discovery of selective agonists. Science 358:381–386.
    .
    OpenUrlAbstract/FREE Full Text
  10. ↵
    1. Miller-Gallacher JL, et al.
    (2014) The 2.1 Å resolution structure of cyanopindolol-bound β1-adrenoceptor identifies an intramembrane Na+ ion that stabilises the ligand-free receptor. PLoS One 9:e92727.
    .
    OpenUrlCrossRefPubMed
  11. ↵
    1. Breiten B, et al.
    (2013) Water networks contribute to enthalpy/entropy compensation in protein-ligand binding. J Am Chem Soc 135:15579–15584.
    .
    OpenUrlCrossRefPubMed
  12. ↵
    1. Levy Y,
    2. Onuchic JN
    (2006) Water mediation in protein folding and molecular recognition. Annu Rev Biophys Biomol Struct 35:389–415.
    .
    OpenUrlCrossRefPubMed
  13. ↵
    1. Bellissent-Funel M-C, et al.
    (2016) Water determines the structure and dynamics of proteins. Chem Rev 116:7673–7697.
    .
    OpenUrlCrossRefPubMed
  14. ↵
    1. Pardo L,
    2. Deupi X,
    3. Dölker N,
    4. López-Rodríguez ML,
    5. Campillo M
    (2007) The role of internal water molecules in the structure and function of the rhodopsin family of G protein-coupled receptors. ChemBioChem 8:19–24.
    .
    OpenUrlCrossRefPubMed
  15. ↵
    1. Angel TE,
    2. Chance MR,
    3. Palczewski K
    (2009) Conserved waters mediate structural and functional activation of family A (rhodopsin-like) G protein-coupled receptors. Proc Natl Acad Sci USA 106:8555–8560.
    .
    OpenUrlAbstract/FREE Full Text
  16. ↵
    1. Latorraca NR,
    2. Venkatakrishnan AJ,
    3. Dror RO
    (2017) GPCR dynamics: Structures in motion. Chem Rev 117:139–155.
    .
    OpenUrl
  17. ↵
    1. Rasmussen SGF, et al.
    (2007) Crystal structure of the human β2 adrenergic G-protein-coupled receptor. Nature 450:383–387.
    .
    OpenUrlCrossRefPubMed
  18. ↵
    1. Rasmussen SGF, et al.
    (2011) Crystal structure of the β2 adrenergic receptor-Gs protein complex. Nature 477:549–555.
    .
    OpenUrlCrossRefPubMed
  19. ↵
    1. Kruse AC, et al.
    (2013) Activation and allosteric modulation of a muscarinic acetylcholine receptor. Nature 504:101–106.
    .
    OpenUrlCrossRefPubMed
  20. ↵
    1. Haga K, et al.
    (2012) Structure of the human M2 muscarinic acetylcholine receptor bound to an antagonist. Nature 482:547–551.
    .
    OpenUrlCrossRefPubMed
  21. ↵
    1. Manglik A, et al.
    (2012) Crystal structure of the µ-opioid receptor bound to a morphinan antagonist. Nature 485:321–326.
    .
    OpenUrlCrossRefPubMed
  22. ↵
    1. Carpenter B,
    2. Nehmé R,
    3. Warne T,
    4. Leslie AGW,
    5. Tate CG
    (2016) Structure of the adenosine A(2A) receptor bound to an engineered G protein. Nature 536:104–107.
    .
    OpenUrlCrossRefPubMed
  23. ↵
    1. Isberg V, et al.
    (2016) GPCRdb: An information system for G protein-coupled receptors. Nucleic Acids Res 44:D356–D364.
    .
    OpenUrlCrossRefPubMed
  24. ↵
    1. Carpenter B,
    2. Tate CG
    (2017) Active state structures of G protein-coupled receptors highlight the similarities and differences in the G protein and arrestin coupling interfaces. Curr Opin Struct Biol 45:124–132.
    .
    OpenUrl
  25. ↵
    1. Jazayeri A,
    2. Dias JM,
    3. Marshall FH
    (2015) From G protein-coupled receptor structure resolution to rational drug design. J Biol Chem 290:19489–19495.
    .
    OpenUrlAbstract/FREE Full Text
  26. ↵
    1. Flock T, et al.
    (2015) Universal allosteric mechanism for Gα activation by GPCRs. Nature 524:173–179.
    .
    OpenUrlCrossRefPubMed
  27. ↵
    1. Isom DG,
    2. Dohlman HG
    (2015) Buried ionizable networks are an ancient hallmark of G protein-coupled receptor activation. Proc Natl Acad Sci USA 112:5702–5707.
    .
    OpenUrlAbstract/FREE Full Text
  28. ↵
    1. Angel TE,
    2. Gupta S,
    3. Jastrzebska B,
    4. Palczewski K,
    5. Chance MR
    (2009) Structural waters define a functional channel mediating activation of the GPCR, rhodopsin. Proc Natl Acad Sci USA 106:14367–14372.
    .
    OpenUrlAbstract/FREE Full Text
  29. ↵
    1. Yuan S,
    2. Filipek S,
    3. Palczewski K,
    4. Vogel H
    (2014) Activation of G-protein-coupled receptors correlates with the formation of a continuous internal water pathway. Nat Commun 5:4733.
    .
    OpenUrlCrossRefPubMed
  30. ↵
    1. Yuan S, et al.
    (2015) The mechanism of ligand-induced activation or inhibition of μ- and κ-opioid receptors. Angew Chem Int Ed Engl 54:7560–7563.
    .
    OpenUrlCrossRef
  31. ↵
    1. Dror RO, et al.
    (2011) Pathway and mechanism of drug binding to G-protein-coupled receptors. Proc Natl Acad Sci USA 108:13118–13123.
    .
    OpenUrlAbstract/FREE Full Text
  32. ↵
    1. Bortolato A,
    2. Tehan BG,
    3. Bodnarchuk MS,
    4. Essex JW,
    5. Mason JS
    (2013) Water network perturbation in ligand binding: Adenosine A(2A) antagonists as a case study. J Chem Inf Model 53:1700–1713.
    .
    OpenUrlCrossRefPubMed
  33. ↵
    1. Tajkhorshid E, et al.
    (2002) Control of the selectivity of the aquaporin water channel family by global orientational tuning. Science 296:525–530.
    .
    OpenUrlAbstract/FREE Full Text
  34. ↵
    1. de Groot BL,
    2. Grubmüller H
    (2001) Water permeation across biological membranes: Mechanism and dynamics of aquaporin-1 and GlpF. Science 294:2353–2357.
    .
    OpenUrlAbstract/FREE Full Text
  35. ↵
    1. Pace CN
    (2009) Energetics of protein hydrogen bonds. Nat Struct Mol Biol 16:681–682.
    .
    OpenUrlCrossRefPubMed
  36. ↵
    1. Snyder PW, et al.
    (2011) Mechanism of the hydrophobic effect in the biomolecular recognition of arylsulfonamides by carbonic anhydrase. Proc Natl Acad Sci USA 108:17889–17894.
    .
    OpenUrlAbstract/FREE Full Text
  37. ↵
    1. Manglik A, et al.
    (2016) Structure-based discovery of opioid analgesics with reduced side effects. Nature 537:185–190.
    .
    OpenUrlCrossRefPubMed
  38. ↵
    1. Zhang Y, et al.
    (2017) Cryo-EM structure of the activated GLP-1 receptor in complex with a G protein. Nature 546:248–253.
    .
    OpenUrlCrossRefPubMed
View Abstract
PreviousNext
Back to top
Article Alerts
Email Article

Thank you for your interest in spreading the word on PNAS.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Diverse GPCRs exhibit conserved water networks for stabilization and activation
(Your Name) has sent you a message from PNAS
(Your Name) thought you would like to see the PNAS web site.
Citation Tools
Diverse GPCRs exhibit conserved water networks for stabilization and activation
A. J. Venkatakrishnan, Anthony K. Ma, Rasmus Fonseca, Naomi R. Latorraca, Brendan Kelly, Robin M. Betz, Chaitanya Asawa, Brian K. Kobilka, Ron O. Dror
Proceedings of the National Academy of Sciences Feb 2019, 116 (8) 3288-3293; DOI: 10.1073/pnas.1809251116

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Request Permissions
Share
Diverse GPCRs exhibit conserved water networks for stabilization and activation
A. J. Venkatakrishnan, Anthony K. Ma, Rasmus Fonseca, Naomi R. Latorraca, Brendan Kelly, Robin M. Betz, Chaitanya Asawa, Brian K. Kobilka, Ron O. Dror
Proceedings of the National Academy of Sciences Feb 2019, 116 (8) 3288-3293; DOI: 10.1073/pnas.1809251116
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Mendeley logo Mendeley
Proceedings of the National Academy of Sciences: 116 (8)
Current Issue

Submit

Sign up for Article Alerts

Jump to section

  • Article
    • Abstract
    • Results
    • Discussion
    • Materials and Methods
    • Acknowledgments
    • Footnotes
    • References
  • Figures & SI
  • Info & Metrics
  • PDF

You May Also be Interested in

News Feature: Cities serve as testbeds for evolutionary change
Urban living can pressure flora and fauna to adapt in intriguing ways. Biologists are starting to take advantage of this convenient laboratory of evolution.
Image credit: Kristin Winchell (Washington University in St. Louis, St. Louis).
Several aspects of the proposal, which aims to expand open access, require serious discussion and, in some cases, a rethink.
Opinion: “Plan S” falls short for society publishers—and for the researchers they serve
Several aspects of the proposal, which aims to expand open access, require serious discussion and, in some cases, a rethink.
Image credit: Dave Cutler (artist).
Featured Profile
PNAS Profile of NAS member and biochemist Hao Wu
 Nonmonogamous strawberry poison frog (Oophaga pumilio).  Image courtesy of Yusan Yang (University of Pittsburgh, Pittsburgh).
Putative signature of monogamy
A study suggests a putative gene-expression hallmark common to monogamous male vertebrates of some species, namely cichlid fishes, dendrobatid frogs, passeroid songbirds, common voles, and deer mice, and identifies 24 candidate genes potentially associated with monogamy.
Image courtesy of Yusan Yang (University of Pittsburgh, Pittsburgh).
Active lifestyles. Image courtesy of Pixabay/MabelAmber.
Meaningful life tied to healthy aging
Physical and social well-being in old age are linked to self-assessments of life worth, and a spectrum of behavioral, economic, health, and social variables may influence whether aging individuals believe they are leading meaningful lives.
Image courtesy of Pixabay/MabelAmber.

More Articles of This Classification

Biological Sciences

  • Major histocompatibility complex class I diversity limits the repertoire of T cell receptors
  • LINC00116 codes for a mitochondrial peptide linking respiration and lipid metabolism
  • Unnatural verticilide enantiomer inhibits type 2 ryanodine receptor-mediated calcium leak and is antiarrhythmic
Show more

Pharmacology

  • Ca2+ allostery in PTH-receptor signaling
  • Effective connectivity changes in LSD-induced altered states of consciousness in humans
Show more

Related Content

  • No related articles found.
  • PubMed
  • Google Scholar

Cited by...

  • No citing articles found.
  • Google Scholar

Similar Articles

Site Logo
Powered by HighWire
  • Submit Manuscript
  • Twitter
  • Facebook
  • RSS Feeds
  • Email Alerts

Articles

  • Current Issue
  • Latest Articles
  • Archive

PNAS Portals

  • Classics
  • Front Matter
  • Teaching Resources
  • Anthropology
  • Chemistry
  • Physics
  • Sustainability Science

Information

  • Authors
  • Editorial Board
  • Reviewers
  • Press
  • Site Map

Feedback    Privacy/Legal

Copyright © 2019 National Academy of Sciences. Online ISSN 1091-6490