The role of Cysteine 6.47 in class A GPCRs

Background: The CWxP motif of transmembrane helix 6 (x: any residue) is highly conserved in class A GPCRs. Within this motif, W6.48 is a big star in the theory of the global "toggle switch" because of its key role in the activation mechanism of GPCRs upon ligand binding. With all footlights focused on W6.48, the reason why the preceding residue, C6.47, is largely conserved is still unknown. The present study is aimed to fill up this lack of knowledge by characterizing the role of C6.47 of the CWxP motif. Results: A complete analysis of available crystal structures has been made alongside with molecular dynamics simulations of model peptides to explore a possible structural role for C6.47. Conclusions: We conclude that C6.47 does not modulate the conformation of the TM6 proline kink and propose that C6.47 participates in the rearrangement of the TM6 and TM7 interface accompanying activation

Inter-residue interactions in membrane proteins

Knowing the precise 3D-structure of a protein is crucial to understand its functional mechanism at the molecular level and to develop new pharmacological agents to targeting it. Nowadays only a few hundred integral membrane protein structures have been solved at high resolution due to the associated technical difficulties. In the present study we aim to characterize the main interactions in alpha and beta membrane proteins that are responsible of the maintenance of the overall structure. With this purpose, two nonredundant databases of alpha and beta transmembrane segments were constructed and analysed. The interactions that stabilize the structure of alpha and beta membrane proteins were quantified. The results reveal important differences in inter-residues interactions between alpha and beta membrane proteins. This novel structural information may be useful in predicting 3D models of proteins lacking structural information or in refining initial models of alpha and beta membrane proteins

Ligand binding mechanisms in human cone visual pigments

Vertebrate vision starts with light absorption by visual pigments in rod and cone photoreceptor cells of the retina. Rhodopsin, in rod cells, responds to dim light, whereas three types of cone opsins (red, green, and blue) function under bright light and mediate color vision. Cone opsins regenerate with retinal much faster than rhodopsin, but the molecular mechanism of regeneration is still unclear. Recent advances in the area pinpoint transient intermediate opsin conformations, and a possible secondary retinal-binding site, as determinant factors for regeneration. In this Review, we compile previous and recent findings to discuss possible mechanisms of ligand entry in cone opsins, involving a secondary binding site, which may have relevant functional and evolutionary implicationsPeer ReviewedPostprint (author's final draft

The size matters? A computational tool to design bivalent ligands

Bivalent ligands are increasingly important such as for targeting G protein-coupled receptor (GPCR) dimers or proteolysis targeting chimeras (PROTACs). They contain two pharmacophoric units that simultaneously bind in their corresponding binding sites, connected with a spacer chain. Here, we report a molecular modelling tool that links the pharmacophore units via the shortest pathway along the receptors van der Waals surface and then scores the solutions providing prioritization for the design of new bivalent ligands. Bivalent ligands of known dimers of GPCRs, PROTACs and a model bivalent antibody/antigen system were analysed. The tool could rapidly assess the preferred linker length for the different systems and recapitulated the best reported results. In the case of GPCR dimers the results suggest that in some cases these ligands might bind to a secondary binding site at the extracellular entrance (vestibule or allosteric site) instead of the orthosteric binding site

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

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

Funding: This work was partially supported by grants from the Spanish Ministry of Economy and Competitiveness (BFU2015-64405-R, SAF2017-84117-R, RTI2018-094204-B-I00 and PID2019- 109240RB-I00; they may include FEDER funds), the Alzheimer's Association (AARFD-17-503612) and by a grant from Fundacio "la Marato" de TV3 (201413-30).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 (AR) and a C-tail truncated (A Δ40 R lacking the last 40 amino acids) adenosine A 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 AR and A Δ40 R. 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

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 (A2AR) and a C-tail truncated (A2AΔ40R lacking the last 40 amino acids) adenosine A2A 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 A2AR and A2AΔ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

Discovery of Homobivalent Bitopic Ligands of the Cannabinoid CB2 Receptor**

Single chemical entities with potential to simultaneously interact with two binding sites are emerging strategies in medicinal chemistry. We have designed, synthesized and functionally characterized the first bitopic ligands for the CB2 receptor. These compounds selectively target CB2 versus CB1 receptors. Their binding mode was studied by molecular dynamic simulations and site-directed mutagenesis

Cross-communication between Gi and Gs in a G-protein-coupled receptor heterotetramer guided by a receptor C-terminal domain

BACKGROUND: G-protein-coupled receptor (GPCR) heteromeric complexes have distinct properties from homomeric GPCRs, giving rise to new receptor functionalities. Adenosine receptors (A1R or A2AR) can form A1R-A2AR heteromers (A1-A2AHet), and their activation leads to canonical G-protein-dependent (adenylate cyclase mediated) and -independent (β-arrestin mediated) signaling. Adenosine has different affinities for A1R and A2AR, allowing the heteromeric receptor to detect its concentration by integrating the downstream Gi- and Gs-dependent signals. cAMP accumulation and β-arrestin recruitment assays have shown that, within the complex, activation of A2AR impedes signaling via A1R. RESULTS: We examined the mechanism by which A1-A2AHet integrates Gi- and Gs-dependent signals. A1R blockade by A2AR in the A1-A2AHet is not observed in the absence of A2AR activation by agonists, in the absence of the C-terminal domain of A2AR, or in the presence of synthetic peptides that disrupt the heteromer interface of A1-A2AHet, indicating that signaling mediated by A1R and A2AR is controlled by both Gi and Gs proteins. CONCLUSIONS: We identified a new mechanism of signal transduction that implies a cross-communication between Gi and Gs proteins guided by the C-terminal tail of the A2AR. This mechanism provides the molecular basis for the operation of the A1-A2AHet as an adenosine concentration-sensing device that modulates the signals originating at both A1R and A2AR