Dual Role of NMDAR Containing NR2A and NR2B Subunits in Alzheimer's Disease

Alzheimer's disease (AD) is the main cause of dementia worldwide. Given that learning and memory are impaired in this pathology, NMDA receptors (NMDARs) appear as key players in the onset and progression of the disease. NMDARs are glutamate receptors, mainly located at the post-synapse, which regulate voltage-dependent influx of calcium into the neurons. They are heterotetramers, and there are different subunits that can be part of the receptors, which are usually composed of two obligatory GluN1 subunits plus either two NR2A or two NR2B subunits. NR2A are mostly located at the synapse, and their activation is involved in the expression of pro-survival genes. Conversely, NR2B are mainly extrasynaptic, and their activation has been related to cell death and neurodegeneration. Thus, activation of NR2A and/or inactivation of NR2B-containing NMDARS has been proposed as a therapeutic strategy to treat AD. Here, we wanted to investigate the main differences between both subunits signalling in neuronal primary cultures of the cortex and hippocampus. It has been observed that Aβ induces a significant increase in calcium release and also in MAPK phosphorylation signalling in NR2B-containing NMDAR in cortical and hippocampal neurons. However, while NR2A-containing NMDAR decreases neuronal death and favours cell viability after Aβ treatment, NR2B-containing NMDAR shows higher levels of cytotoxicity and low levels of neuronal survival. Finally, it has been detected that NMDAR has no effect on pTau axonal transport. The present results demonstrate a different role between GluNA and GluNB subunits in neurodegenerative diseases such as Alzheimer's.

Discovery of a macromolecular complex mediating the hunger suppressive actions of cocaine: Structural and functional properties

Cocaine not only increases brain dopamine levels but also activates the sigma1 receptor (σ1 R) that in turn regulates orexigenic receptor function. Identification of interactions involving dopamine D1 (D1 R), ghrelin (GHS-R1a ), and σ1 receptors have been addressed by biophysical techniques and a complementation approach using interfering peptides. The effect of cocaine on receptor functionality was assayed by measuring second messenger, cAMP and Ca2+ , levels. The effect of acute or chronic cocaine administration on receptor complex expression was assayed by in situ proximity ligation assay. In silico procedures were used for molecular model building. σ1 R KO mice were used for confirming involvement of this receptor. Upon identification of protomer interaction and receptor functionality, a unique structural model for the macromolecular complex formed by σ1 R, D1 R, and GHS-R1a is proposed. The functionality of the complex, able to couple to both Gs and Gq proteins, is affected by cocaine binding to the σ1 R, as confirmed using samples from σ1 R-/- mice. The expression of the macromolecular complex was differentially affected upon acute and chronic cocaine administration to rats. The constructed 3D model is consistent with biochemical, biophysical, and available structural data. The σ1 R, D1 R, and GHS-R1a complex constitutes a functional unit that is altered upon cocaine binding to the σ1 R. Remarkably, the heteromer can simultaneously couple to two G proteins, thus allowing dopamine to signal via Ca2+ and ghrelin via cAMP. The anorexic action of cocaine is mediated by such complex whose expression is higher after acute than after chronic administration regimens.

Design of Negative and Positive Allosteric Modulators of the Cannabinoid CB2 Receptor Derived from the Natural Product Cannabidiol

Cannabidiol (CBD), the second most abundant of the active compounds found in the Cannabis sativa plant, is of increasing interest because it is approved for human use and is neither euphorizing nor addictive. Here, we design and synthesize novel compounds taking into account that CBD is both a partial agonist, when it binds to the orthosteric site, and a negative allosteric modulator, when it binds to the allosteric site of the cannabinoid CB₂ receptor. Molecular dynamic simulations and site-directed mutagenesis studies have identified the allosteric site near the receptor entrance. This knowledge has permitted to perform structure-guided design of negative and positive allosteric modulators of the CB₂ receptor with potential therapeutic utility.

Structure and function of adenosine receptor heteromers

Adenosine is one of the most ancient signaling molecules and has receptors in both animals and plants. In mammals there are four specific receptors, A1, A2A, A2B, and A3, which belong to the superfamily of G-protein-coupled receptors (GPCRs). Evidence accumulated in the last 20 years indicates that GPCRs are often expressed as oligomeric complexes formed by a number of equal (homomers) or different (heteromers) receptors. This review presents the data showing the occurrence of heteromers formed by A1 and A2A, A2A and A2B, and A2A and A3 receptors highlighting (i) their tetrameric structural arrangements, and (ii) the functional diversity that those heteromers provide to adenosinergic signaling.

TMSNP: a web server to predict pathogenesis of missense mutations in the transmembrane region of membrane proteins

The massive amount of data generated from genome sequencing brings tons of newly identified mutations, whose pathogenic/non-pathogenic effects need to be evaluated. This has given rise to several mutation predictor tools that, in general, do not consider the specificities of the various protein groups. We aimed to develop a predictor tool dedicated to membrane proteins, under the premise that their specific structural features and environment would give different responses to mutations compared to globular proteins. For this purpose, we created TMSNP, a database that currently contains information from 2624 pathogenic and 196 705 non-pathogenic reported mutations located in the transmembrane region of membrane proteins. By computing various conservation parameters on these mutations in combination with annotations, we trained a machine-learning model able to classify mutations as pathogenic or not. TMSNP (freely available at http://lmc.uab.es/tmsnp/) improves considerably the prediction power of commonly used mutation predictors trained with globular proteins.

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.

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.

DIMERBOW: exploring possible GPCR dimer interfaces

Motivation

G protein-coupled receptors (GPCRs) can form homo-, heterodimers and larger order oligomers that exert different functions than monomers. The pharmacological potential of such complexes is hampered by the limited information available on the type of complex formed and its quaternary structure. Several GPCR structures in the Protein Data Bank display crystallographic interfaces potentially compatible with physiological interactions.

Results

Here, we present DIMERBOW, a database and web application aimed to visually browse the complete repertoire of potential GPCR dimers present in solved structures. The tool is suited to help finding the best possible structural template to model GPCR homomers.

Availability and implementation

DIMERBOW is available at http://lmc.uab.es/dimerbow/.

Supplementary information

Supplementary data are available at Bioinformatics online.

Control of glutamate release by complexes of adenosine and cannabinoid receptors

Background

It has been hypothesized that heteromers of adenosine A_2A receptors (A2AR) and cannabinoid CB₁ receptors (CB1R) localized in glutamatergic nerve terminals mediate the integration of adenosine and endocannabinoid signaling involved in the modulation of striatal excitatory neurotransmission. Previous studies have demonstrated the existence of A2AR-CB1R heteromers in artificial cell systems. A dependence of A2AR signaling for the Gi protein-mediated CB1R signaling was described as one of its main biochemical characteristics. However, recent studies have questioned the localization of functionally significant A2AR-CB1R heteromers in striatal glutamatergic terminals.

Results

Using a peptide-interfering approach combined with biophysical and biochemical techniques in mammalian transfected cells and computational modeling, we could establish a tetrameric quaternary structure of the A2AR-CB1R heterotetramer. This quaternary structure was different to the also tetrameric structure of heteromers of A2AR with adenosine A₁ receptors or dopamine D₂ receptors, with different heteromeric or homomeric interfaces. The specific quaternary structure of the A2A-CB1R, which depended on intermolecular interactions involving the long C-terminus of the A2AR, determined a significant A2AR and Gs protein-mediated constitutive activation of adenylyl cyclase. Using heteromer-interfering peptides in experiments with striatal glutamatergic terminals, we could then demonstrate the presence of functionally significant A2AR-CB1R heteromers with the same biochemical characteristics of those studied in mammalian transfected cells. First, either an A2AR agonist or an A2AR antagonist allosterically counteracted Gi-mediated CB1R agonist-induced inhibition of depolarization-induced glutamate release. Second, co-application of both an A2AR agonist and an antagonist cancelled each other effects. Finally, a CB1R agonist inhibited glutamate release dependent on a constitutive activation of A2AR by a canonical Gs-Gi antagonistic interaction at the adenylyl cyclase level.

Conclusions

We demonstrate that the well-established cannabinoid-induced inhibition of striatal glutamate release can mostly be explained by a CB1R-mediated counteraction of the A2AR-mediated constitutive activation of adenylyl cyclase in the A2AR-CB1R heteromer.

The size matters? A computational tool to design bivalent ligands

Motivation

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.

Results

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.

Availability and implementation

Freely accessible from the Molecular Operating Environment svl exchange server (https://svl.chemcomp.com/).

Supplementary information

Supplementary data are available at Bioinformatics online.

Inter-residue interactions in alpha-helical transmembrane proteins

MOTIVATION

The number of available membrane protein structures has markedly increased in the last years and, in parallel, the reliability of the methods to detect transmembrane (TM) segments. In the present report, we characterized inter-residue interactions in α-helical membrane proteins using a dataset of 3462 TM helices from 430 proteins. This is by far the largest analysis published to date.

RESULTS

Our analysis of residue-residue interactions in TM segments of membrane proteins shows that almost all interactions involve aliphatic residues and Phe. There is lack of polar-polar, polar-charged and charged-charged interactions except for those between Thr or Ser sidechains and the backbone carbonyl of aliphatic and Phe residues. The results are discussed in the context of the preferences of amino acids to be in the protein core or exposed to the lipid bilayer and to occupy specific positions along the TM segment. Comparison to datasets of β-barrel membrane proteins and of α-helical globular proteins unveils the specific patterns of interactions and residue composition characteristic of α-helical membrane proteins that are the clue to understanding their structure.

AVAILABILITY AND IMPLEMENTATION

Results data and datasets used are available at http://lmc.uab.cat/TMalphaDB/interactions.php.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

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 implications.

Cocaine Blocks Effects of Hunger Hormone, Ghrelin, Via Interaction with Neuronal Sigma-1 Receptors

Despite ancient knowledge on cocaine appetite-suppressant action, the molecular basis of such fact remains unknown. Addiction/eating disorders (e.g., binge eating, anorexia, bulimia) share a central control involving reward circuits. However, we here show that the sigma-1 receptor (σ1R) mediates cocaine anorectic effects by interacting in neurons with growth/hormone/secretagogue (ghrelin) receptors. Cocaine increases colocalization of σ1R and GHS-R1a at the cell surface. Moreover, in transfected HEK-293T and neuroblastoma SH-SY5Y cells, and in primary neuronal cultures, pretreatment with cocaine or a σ1R agonist inhibited ghrelin-mediated signaling, in a similar manner as the GHS-R1a antagonist YIL-781. Results were similar in G protein-dependent (cAMP accumulation and calcium release) and in partly dependent or independent (ERK1/2 phosphorylation and label-free) assays. We provide solid evidence for direct interaction between receptors and the functional consequences, as well as a reliable structural model of the macromolecular σ1R-GHS-R1a complex, which arises as a key piece in the puzzle of the events linking cocaine consumption and appetitive/consummatory behaviors.

Design of a True Bivalent Ligand with Picomolar Binding Affinity for a G Protein-Coupled Receptor Homodimer

Bivalent ligands have emerged as chemical tools to study G protein-coupled receptor dimers. Using a combination of computational, chemical, and biochemical tools, here we describe the design of bivalent ligand 13 with high affinity ( K_DB1 = 21 pM) for the dopamine D₂ receptor (D₂R) homodimer. Bivalent ligand 13 enhances the binding affinity relative to monovalent compound 15 by 37-fold, indicating simultaneous binding at both protomers. Using synthetic peptides with amino acid sequences of transmembrane (TM) domains of D₂R, we provide evidence that TM6 forms the interface of the homodimer. Notably, the disturber peptide TAT-TM6 decreased the binding of bivalent ligand 13 by 52-fold and had no effect on monovalent compound 15, confirming the D₂R homodimer through TM6 ex vivo. In conclusion, by using a versatile multivalent chemical platform, we have developed a precise strategy to generate a true bivalent ligand that simultaneously targets both orthosteric sites of the D₂R homodimer.

Ligand-Triggered Structural Changes in the M2 Muscarinic Acetylcholine Receptor

The muscarinic M₂ acetylcholine receptor, one of the few G-protein coupled receptors that has not only been crystallized in both active and inactive conformations but also in the presence of a positive allosteric modulator, is an interesting system to study the molecular mechanisms of GPCR activation and ligand allosterism. Here, we have employed molecular dynamics (MD) simulations (adding to 14 μs in total) to study conformational changes triggered by the inverse agonist R-(-)-3-quinuclidinyl-benzilate (QNB) in the structure of the active M₂ receptor (PBD ID 4MQS ) after replacement of the agonist iperoxo by the inverse agonist QNB. This permitted us to identify the sequence of events in the deactivation mechanism of the M₂ acetylcholine receptor, which results first in the rearrangement of the transmission switch, the subsequent opening of the extracellular portion of the receptor and finally, the closure of the intracellular part. We also evaluate the effect of the positive allosteric modulator LY2119620 when bound simultaneously with the orthosteric agonist iperoxo and find that it restricts the conformation of Trp422^7.35 in a position that modulates the orientation of the Tyr426^7.39 at the orthosteric-binding pocket.

Essential Control of the Function of the Striatopallidal Neuron by Pre-coupled Complexes of Adenosine A2A-Dopamine D2 Receptor Heterotetramers and Adenylyl Cyclase

The central adenosine system and adenosine receptors play a fundamental role in the modulation of dopaminergic neurotransmission. This is mostly achieved by the strategic co-localization of different adenosine and dopamine receptor subtypes in the two populations of striatal efferent neurons, striatonigral and striatopallidal, that give rise to the direct and indirect striatal efferent pathways, respectively. With optogenetic techniques it has been possible to dissect a differential role of the direct and indirect pathways in mediating "Go" responses upon exposure to reward-related stimuli and "NoGo" responses upon exposure to non-rewarded or aversive-related stimuli, respectively, which depends on their different connecting output structures and their differential expression of dopamine and adenosine receptor subtypes. The striatopallidal neuron selectively expresses dopamine D₂ receptors (D2R) and adenosine A_2A receptors (A2AR), and numerous experiments using multiple genetic and pharmacological in vitro, in situ and in vivo approaches, demonstrate they can form A2AR-D2R heteromers. It was initially assumed that different pharmacological interactions between dopamine and adenosine receptor ligands indicated the existence of different subpopulations of A2AR and D2R in the striatopallidal neuron. However, as elaborated in the present essay, most evidence now indicates that all interactions can be explained with a predominant population of striatal A2AR-D2R heteromers forming complexes with adenylyl cyclase subtype 5 (AC5). The A2AR-D2R heteromer has a tetrameric structure, with two homodimers, which allows not only multiple allosteric interactions between different orthosteric ligands, agonists, and antagonists, but also the canonical Gs-Gi antagonistic interaction at the level of AC5. We present a model of the function of the A2AR-D2R heterotetramer-AC5 complex, which acts as an integrative device of adenosine and dopamine signals that determine the excitability and gene expression of the striatopallidal neurons. The model can explain most behavioral effects of A2AR and D2R ligands, including the psychostimulant effects of caffeine. The model is also discussed in the context of different functional striatal compartments, mainly the dorsal and the ventral striatum. The current accumulated knowledge of the biochemical properties of the A2AR-D2R heterotetramer-AC5 complex offers new therapeutic possibilities for Parkinson's disease, schizophrenia, SUD and other neuropsychiatric disorders with dysfunction of dorsal or ventral striatopallidal neurons.

Human Blue Cone Opsin Regeneration Involves Secondary Retinal Binding with Analog Specificity

Human color vision is mediated by the red, green, and blue cone visual pigments. Cone opsins are G-protein-coupled receptors consisting of an opsin apoprotein covalently linked to the 11-cis-retinal chromophore. All visual pigments share a common evolutionary origin, and red and green cone opsins exhibit a higher homology, whereas blue cone opsin shows more resemblance to the dim light receptor rhodopsin. Here we show that chromophore regeneration in photoactivated blue cone opsin exhibits intermediate transient conformations and a secondary retinoid binding event with slower binding kinetics. We also detected a fine-tuning of the conformational change in the photoactivated blue cone opsin binding site that alters the retinal isomer binding specificity. Furthermore, the molecular models of active and inactive blue cone opsins show specific molecular interactions in the retinal binding site that are not present in other opsins. These findings highlight the differential conformational versatility of human cone opsin pigments in the chromophore regeneration process, particularly compared to rhodopsin, and point to relevant functional, unexpected roles other than spectral tuning for the cone visual pigments.

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 (A₁R or A_2AR) can form A₁R-A_2AR heteromers (A₁-A_2AHet), and their activation leads to canonical G-protein-dependent (adenylate cyclase mediated) and -independent (β-arrestin mediated) signaling. Adenosine has different affinities for A₁R and A_2AR, allowing the heteromeric receptor to detect its concentration by integrating the downstream G_i- and G_s-dependent signals. cAMP accumulation and β-arrestin recruitment assays have shown that, within the complex, activation of A_2AR impedes signaling via A₁R.

Results

We examined the mechanism by which A₁-A_2AHet integrates G_i- and G_s-dependent signals. A₁R blockade by A_2AR in the A₁-A_2AHet is not observed in the absence of A_2AR activation by agonists, in the absence of the C-terminal domain of A_2AR, or in the presence of synthetic peptides that disrupt the heteromer interface of A₁-A_2AHet, indicating that signaling mediated by A₁R and A_2AR is controlled by both G_i and G_s proteins.

Conclusions

We identified a new mechanism of signal transduction that implies a cross-communication between G_i and G_s proteins guided by the C-terminal tail of the A_2AR. This mechanism provides the molecular basis for the operation of the A₁-A_2AHet as an adenosine concentration-sensing device that modulates the signals originating at both A₁R and A_2AR.

Electronic supplementary material

The online version of this article (10.1186/s12915-018-0491-x) contains supplementary material, which is available to authorized users.

Analysis of the interactions of sulfur-containing amino acids in membrane proteins

The interactions of Met and Cys with other amino acid side chains have received little attention, in contrast to aromatic-aromatic, aromatic-aliphatic or/and aliphatic-aliphatic interactions. Precisely, these are the only amino acids that contain a sulfur atom, which is highly polarizable and, thus, likely to participate in strong Van der Waals interactions. Analysis of the interactions present in membrane protein crystal structures, together with the characterization of their strength in small-molecule model systems at the ab-initio level, predicts that Met-Met interactions are stronger than Met-Cys ≈ Met-Phe ≈ Cys-Phe interactions, stronger than Phe-Phe ≈ Phe-Leu interactions, stronger than the Met-Leu interaction, and stronger than Leu-Leu ≈ Cys-Leu interactions. These results show that sulfur-containing amino acids form stronger interactions than aromatic or aliphatic amino acids. Thus, these amino acids may provide additional driving forces for maintaining the 3D structure of membrane proteins and may provide functional specificity.

Beyond spectral tuning: human cone visual pigments adopt different transient conformations for chromophore regeneration

Human red and green visual pigments are seven transmembrane receptors of cone photoreceptor cells of the retina that mediate color vision. These pigments share a very high degree of homology and have been assumed to feature analogous structural and functional properties. We report on a different regeneration mechanism among red and green cone opsins with retinal analogs using UV-Vis/fluorescence spectroscopic analyses, molecular modeling and site-directed mutagenesis. We find that photoactivated green cone opsin adopts a transient conformation which regenerates via an unprotonated Schiff base linkage with its natural chromophore, whereas red cone opsin forms a typical protonated Schiff base. The chromophore regeneration kinetics is consistent with a secondary retinal uptake by the cone pigments. Overall, our findings reveal, for the first time, structural differences in the photoactivated conformation between red and green cone pigments that may be linked to their molecular evolution, and support the proposal of secondary retinal binding to visual pigments, in addition to binding to the canonical primary site, which may serve as a regulatory mechanism of dark adaptation in the phototransduction process.

Effect of Force Field Parameters on Sodium and Potassium Ion Binding to Dipalmitoyl Phosphatidylcholine Bilayers

The behavior of electrolytes in molecular dynamics simulations of zwitterionic phospholipid bilayers is very sensitive to the force field parameters used. Here, several 200 ns molecular dynamics of simulations of dipalmitoyl phosphotidylcholine (PC) bilayers in 0.2 M sodium or potassium chloride using various common force field parameters for the cations are presented. All employed parameter sets give a larger number of Na(+) ions than K(+) ions that bind to the lipid heads, but depending on the parameter choice quite different results are seen. A wide range of coordination numbers for the Na(+) and K(+) ions is also observed. These findings have been analyzed and compared to published experimental data. Some simulations produce aggregates of potassium chloride, indicating (in accordance with published simulations) that these force fields do not reproduce the delicate balance between salt and solvated ions. The differences between the force fields can be characterized by one single parameter, the electrostatic radius of the ion, which is correlated to σMO (M represents Na(+)/K(+)), the Lennard-Jones radius. When this parameter exceeds a certain threshold, binding to the lipid heads is no longer observed. One would, however, need more accurate experimental data to judge or rank the different force fields precisely. Still, reasons for the poor performance of some of the parameter sets are clearly demonstrated, and a quality control procedure is provided.

TMalphaDB and TMbetaDB: web servers to study the structural role of sequence motifs in α-helix and β-barrel domains of membrane proteins

Background

Membrane proteins represent over 25 % of human protein genes and account for more than 60 % of drug targets due to their accessibility from the extracellular environment. The increasing number of available crystal structures of these proteins in the Protein Data Bank permits an initial estimation of their structural properties.

Description

We have developed two web servers—TMalphaDB for α-helix bundles and TMbetaDB for β-barrels—to analyse the growing repertoire of available crystal structures of membrane proteins. TMalphaDB and TMbetaDB permit to search for these specific sequence motifs in a non-redundant structure database of transmembrane segments and quantify structural parameters such as ϕ and ψ backbone dihedral angles, χ₁ side chain torsion angle, unit bend and unit twist.

Conclusions

The structural information offered by TMalphaDB and TMbetaDB permits to quantify structural distortions induced by specific sequence motifs, and to elucidate their role in the 3D structure. This specific structural information has direct implications in homology modeling of the growing sequences of membrane proteins lacking experimental structure. TMalphaDB and TMbetaDB are freely available at http://lmc.uab.cat/TMalphaDB and http://lmc.uab.cat/TMbetaDB.

Cognitive Impairment Induced by Delta9-tetrahydrocannabinol Occurs through Heteromers between Cannabinoid CB1 and Serotonin 5-HT2A Receptors

Activation of cannabinoid CB1 receptors (CB1R) by delta9-tetrahydrocannabinol (THC) produces a variety of negative effects with major consequences in cannabis users that constitute important drawbacks for the use of cannabinoids as therapeutic agents. For this reason, there is a tremendous medical interest in harnessing the beneficial effects of THC. Behavioral studies carried out in mice lacking 5-HT2A receptors (5-HT2AR) revealed a remarkable 5-HT2AR-dependent dissociation in the beneficial antinociceptive effects of THC and its detrimental amnesic properties. We found that specific effects of THC such as memory deficits, anxiolytic-like effects, and social interaction are under the control of 5-HT2AR, but its acute hypolocomotor, hypothermic, anxiogenic, and antinociceptive effects are not. In biochemical studies, we show that CB1R and 5-HT2AR form heteromers that are expressed and functionally active in specific brain regions involved in memory impairment. Remarkably, our functional data shows that costimulation of both receptors by agonists reduces cell signaling, antagonist binding to one receptor blocks signaling of the interacting receptor, and heteromer formation leads to a switch in G-protein coupling for 5-HT2AR from Gq to Gi proteins. Synthetic peptides with the sequence of transmembrane helices 5 and 6 of CB1R, fused to a cell-penetrating peptide, were able to disrupt receptor heteromerization in vivo, leading to a selective abrogation of memory impairments caused by exposure to THC. These data reveal a novel molecular mechanism for the functional interaction between CB1R and 5-HT2AR mediating cognitive impairment. CB1R-5-HT2AR heteromers are thus good targets to dissociate the cognitive deficits induced by THC from its beneficial antinociceptive properties.

Gastrin-stimulated Gα13 Activation of Rgnef Protein (ArhGEF28) in DLD-1 Colon Carcinoma Cells

The guanine nucleotide exchange factor Rgnef (also known as ArhGEF28 or p190RhoGEF) promotes colon carcinoma cell motility and tumor progression via interaction with focal adhesion kinase (FAK). Mechanisms of Rgnef activation downstream of integrin or G protein-coupled receptors remain undefined. In the absence of a recognized G protein signaling homology domain in Rgnef, no proximal linkage to G proteins was known. Utilizing multiple methods, we have identified Rgnef as a new effector for Gα13 downstream of gastrin and the type 2 cholecystokinin receptor. In DLD-1 colon carcinoma cells depleted of Gα13, gastrin-induced FAK Tyr(P)-397 and paxillin Tyr(P)-31 phosphorylation were reduced. RhoA GTP binding and promoter activity were increased by Rgnef in combination with active Gα13. Rgnef co-immunoprecipitated with activated Gα13Q226L but not Gα12Q229L. The Rgnef C-terminal (CT, 1279-1582) region was sufficient for co-immunoprecipitation, and Rgnef-CT exogenous expression prevented Gα13-stimulated SRE activity. A domain at the C terminus of the protein close to the FAK binding domain is necessary to bind to Gα13. Point mutations of Rgnef-CT residues disrupt association with active Gα13 but not Gαq. These results show that Rgnef functions as an effector of Gα13 signaling and that this linkage may mediate FAK activation in DLD-1 colon carcinoma cells.

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.

Binding specificity of retinal analogs to photoactivated visual pigments suggest mechanism for fine-tuning GPCR-ligand interactions

11-cis-retinal acts as an inverse agonist stabilizing the inactive conformation of visual pigments, and upon photoactivation, it isomerizes to all-trans-retinal, initiating signal transduction. We have analyzed opsin regeneration with retinal analogs for rhodopsin and red cone opsin. We find differential binding of the analogs to the receptors after photobleaching and a dependence of the binding kinetics on the oligomerization state of the protein. The results outline the sensitivity of retinal entry to the binding pocket of visual receptors to the specific conformation adopted by the receptor and by the molecular architecture defined by specific amino acids in the binding pocket and the retinal entry site, as well as the topology of the retinal analog. Overall, our findings highlight the specificity of the ligand-opsin interactions, a feature that can be shared by other G-protein-coupled receptors.

Ligand binding determinants for angiotensin II type 1 receptor from computer simulations

The ligand binding determinants for the angiotensin II type 1 receptor (AT1R), a G protein-coupled receptor (GPCR), have been characterized by means of computer simulations. As a first step, a pharmacophore model of various known AT1R ligands exhibiting a wide range of binding affinities was generated. Second, a structural model of AT1R was built making use of the growing set of crystal structures of GPCRs, which was further used for the docking of the AT1R ligands based on the devised pharmacophore model. Next, ligand-receptor-lipid bilayer systems were studied by means of molecular dynamics (MD) simulations. Overall, the present study has permitted, combining the pharmacophore model with binding free energy calculations obtained from the MD simulations, to propose the molecular mechanisms by which sartans interact with AT1R.

Sulfur-containing amino acids in 7TMRs: molecular gears for pharmacology and function

Seven-transmembrane receptors (7TMRs) mediate the majority of physiological responses to hormones and neurotransmitters in higher organisms. Tertiary structure stability and activation of these versatile membrane proteins require formation or disruption of complex networks of well-recognized interactions (such as H-bonds, ionic, or aromatic-aromatic) but also of other type of interactions which have been less studied. In this review, we compile evidence from crystal structure, biophysical, and site-directed mutagenesis data that indicate or support the importance of interactions involving Met and Cys in 7TMRs in terms of pharmacology and function. We show examples of Met/Cys-aromatic and Met-Met interactions participating in ligand binding, in tuning the orientation of functionally important aromatic residues during activation or even in modulating the type of signaling response. Collectively, data presented enlarge the repertoire of interactions governing 7TMR functioning.

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.

Membrane Protein Simulations Using AMBER Force Field and Berger Lipid Parameters

AMBER force fields are among the most commonly used in molecular dynamics (MD) simulations of proteins. Unfortunately, they lack a specific set of lipid parameters, thus limiting its use in membrane protein simulations. In order to overcome this limitation we assessed whether the widely used united-atom lipid parameters described by Berger and co-workers could be used in conjunction with AMBER force fields in simulations of membrane proteins. Thus, free energies of solvation in water and in cyclohexane, and free energies of water to cyclohexane transfer, were computed by thermodynamic integration procedures for neutral amino acid side-chains employing AMBER99, AMBER03, and OPLS-AA amino acid force fields. In addition, MD simulations of three membrane proteins in a POPC lipid bilayer, the β2 adrenergic G protein-coupled receptor, Aquaporin-1, and the outer membrane protein Omp32, were performed with the aim of comparing the AMBER99SB/Berger combination of force fields with the OPLS-AA/Berger combination. We have shown that AMBER99SB and Berger force fields are compatible, they provide reliable free energy estimations relative to experimental values, and their combination properly describes both membrane and protein structural properties. We then suggest that the AMBER99SB/Berger combination is a reliable choice for the simulation of membrane proteins, which links the easiness of ligand parametrization and the ability to reproduce secondary structure of AMBER99SB force field with the largely validated Berger lipid parameters.

Molecular mechanisms of disease for mutations at Gly-90 in rhodopsin

Two different mutations at Gly-90 in the second transmembrane helix of the photoreceptor protein rhodopsin have been proposed to lead to different phenotypes. G90D has been classically associated with congenital night blindness, whereas the newly reported G90V substitution was linked to a retinitis pigmentosa phenotype. Here, we used Val/Asp replacements of the native Gly at position 90 to unravel the structure/function divergences caused by these mutations and the potential molecular mechanisms of inherited retinal disease. The G90V and G90D mutants have a similar conformation around the Schiff base linkage region in the dark state and same regeneration kinetics with 11-cis-retinal, but G90V has dramatically reduced thermal stability when compared with the G90D mutant rhodopsin. The G90V mutant also shows, like G90D, an altered photobleaching pattern and capacity to activate Gt in the opsin state. Furthermore, the regeneration of the G90V mutant with 9-cis-retinal was improved, achieving the same A(280)/A(500) as wild type isorhodopsin. Hydroxylamine resistance was also recovered, indicating a compact structure around the Schiff base linkage, and the thermal stability was substantially improved when compared with the 11-cis-regenerated mutant. These results support the role of thermal instability and/or abnormal photoproduct formation in eliciting a retinitis pigmentosa phenotype. The improved stability and more compact structure of the G90V mutant when it was regenerated with 9-cis-retinal brings about the possibility that this isomer or other modified retinoid analogues might be used in potential treatment strategies for mutants showing the same structural features.

Interactions of fatty acids with phosphatidylethanolamine membranes: X-ray diffraction and molecular dynamics studies

An experimental and theoretical study on 1,2-dielaidoyl-sn-glycero-3-phosphoethanolamine (DEPE) membranes containing fatty acids (FAs) was performed by means of X-ray diffraction analysis and molecular dynamics (MD) simulations. The study was aimed at understanding the interactions of several structurally related FAs with biomembranes, which is necessary for further rational lipid drug design in membrane-lipid therapy. The main effect of FAs was to promote the formation of a H(II) phase, despite a stabilization of the coexisting L(alpha) + H(II) phases. Derivatives of OA exhibited a specific density profile in the direction perpendicular to the bilayer that reflects differences in the relative localization of the carboxylate group within the polar region of the membrane as well as in the degree of membrane penetration of the FA acyl chain. Hydroxyl and methyl substituents at carbon-2 in the FA acyl chain were identified as effective modulators of the position of carboxylate group in the lipid bilayer. Our data highlight the specific potential of each FA in modulating the membrane structure properties.

Effect of ions on a dipalmitoyl phosphatidylcholine bilayer. a molecular dynamics simulation study

The effect of physiological concentrations of different chlorides on the structure of a dipalmitoyl phosphatidylcholine (DPPC) bilayer has been investigated through atomistic molecular dynamics simulations. These calculations provide support to the concept that Li+, Na+, Ca2+, Mg2+, Sr2+, Ba2+, and Ac3+, but not K+, bind to the lipid-head oxygens. Ion binding exhibits an influence on lipid order, area per lipid, orientation of the lipid head dipole, the charge distribution in the system, and therefore the electrostatic potential across the head-group region of the bilayer. These structural effects are sensitive to the specific characteristics of each cation, i.e., radius, charge, and coordination properties. These results provide evidence aimed at shedding some light into the apparent contradictions among different studies reported recently regarding the ordering effect of ions on zwitterionic phosphatidylcholine lipid bilayers.

Molecular Dynamics Simulations of Rhodopsin in Different One-Component Lipid Bilayers

Four 20 ns molecular dynamic simulations of rhodopsin embedded in different one-component lipid bilayers have been carried out to ascertain the importance of membrane lipids on the protein structure. Specifically, dimyristoyl phosphatidylcholine (DMPC), dipalmitoyl phosphatidylcholine (DPPC), palmitoyl oleoyl phosphatidylcholine (POPC), and palmitoyl linoleyl phosphatidylcholine (PLPC) lipid bilayers have been considered for the present work. The results reported here provide information on the hydrophobic matching between the protein and the bilayer and about the differential effects of the protein on the thickness of the different membranes. Furthermore, a careful analysis of the individual protein-lipid interactions permits the identification of residues that exhibit permanent interactions with atoms of the lipid environment that may putatively act as hooks of the protein to the membrane. The analysis of the trajectories also provides information about the effect of the bilayer on the protein structure, including secondary structural elements, salt bridges, and rigid-body motions.

Critical role of electrostatic interactions of amino acids at the cytoplasmic region of helices 3 and 6 in rhodopsin conformational properties and activation

The cytoplasmic sides of transmembrane helices 3 and 6 of G-protein-coupled receptors are connected by a network of ionic interactions that play an important role in maintaining its inactive conformation. To investigate the role of such a network in rhodopsin structure and function, we have constructed single mutants at position 134 in helix 3 and at positions 247 and 251 in helix 6, as well as combinations of these to obtain double mutants involving the two helices. These mutants have been expressed in COS-1 cells, immunopurified using the rho-1D4 antibody, and studied by UV-visible spectrophotometry. Most of the single mutations did not affect chromophore formation, but double mutants, especially those involving the T251K mutant, resulted in low yield of protein and impaired 11-cis-retinal binding. Single mutants E134Q, E247Q, and E247A showed the ability to activate transducin in the dark, and E134Q and E247A enhanced activation upon illumination, with regard to wild-type rhodopsin. Mutations E247A and T251A (in E134Q/E247A and E134Q/T251A double mutants) resulted in enhanced activation compared with the single E134Q mutant in the dark. A role for Thr(251) in this network is proposed for the first time in rhodopsin. As a result of these mutations, alterations in the hydrogen bond interactions between the amino acid side chains at the cytoplasmic region of transmembrane helices 3 and 6 have been observed using molecular dynamics simulations. Our combined experimental and modeling results provide new insights into the details of the structural determinants of the conformational change ensuing photoactivation of rhodopsin.

X-ray absorption and molecular dynamics study of cation binding sites in the purple membrane

The present work describes the results of a study aimed at identifying candidate cation binding sites on the extracellular region of bacteriorhodopsin, including a site near the retinal pocket. The approach used is a combined effort involving computational chemistry methods (computation of cation affinity maps and molecular dynamics) together with the Extended X-Ray Absorption Fine Structure (EXAFS) technique to obtain relevant information about the local structure of the protein in the neighborhood of Mn(2+) ions in different affinity binding sites. The results permit the identification of a high-affinity binding site where the ion is coordinated simultaneously to Asp212(-) and Asp85(-). Comparison of EXAFS data of the wild type protein with the quadruple mutant E9Q/E74Q/E194Q/E204Q at pH 7.0 and 10.0 demonstrate that extracellular glutamic acid residues are involved in cation binding.

Effect of different treatments of long-range interactions and sampling conditions in molecular dynamic simulations of rhodopsin embedded in a dipalmitoyl phosphatidylcholine bilayer

The present study analyzes the effect of the simulation conditions on the results of molecular dynamics simulations of G-protein coupled receptors (GPCRs) performed with an explicit lipid bilayer. Accordingly, the present work reports the analysis of different simulations of bovine rhodopsin embedded in a dipalmitoyl phosphatidylcholine (DPPC) lipid bilayer using two different sampling conditions and two different approaches for the treatment of long-range electrostatic interactions. Specifically, sampling was carried out either by using the statistical ensembles NVT or NPT (constant number of atoms, a pressure of 1 atm in all directions and fixed temperature), and the electrostatic interactions were treated either by using a twin-cutoff, or the particle mesh Ewald summation method (PME). The results of the present study suggest that the use of the NPT ensemble in combination with the PME method provide more realistic simulations. The use of NPT during the equilibration avoids the need of an a priori estimation of the box dimensions, giving the correct area per lipid. However, once the system is equilibrated, the simulations are irrespective of the sampling conditions used. The use of an electrostatic cutoff induces artifacts on both lipid thickness and the ion distribution, but has no direct effect on the protein and water molecules.

A methyl group at C7 of 11-cis-retinal allows chromophore formation but affects rhodopsin activation

The newly synthesized 11-cis-7-methylretinal can form an artificial visual pigment with kinetic and spectroscopic properties similar to the native pigment in the dark-state. However, its photobleaching behavior is altered, showing a Meta I-like photoproduct. This behavior reflects a steric constraint imposed by the 7-methyl group that affects the conformational change in the binding pocket as a result of retinal photoisomerization. Transducin activation is reduced, when compared to the native pigment with 11-cis-retinal. Molecular dynamics simulations suggest coupling of the C7 methyl group and the beta-ionone ring with Met207 in transmembrane helix 5 in agreement with recent experimental results.