Molecular dynamics simulations and structure-based network analysis reveal structural and functional aspects of G-protein coupled receptor dimer interactions

Cholesterol Dependence of the Conformational Changes in Metabotropic Glutamate Receptor 1

Metabotropic glutamate receptors (mGluRs) are class C G protein-coupled receptors that function as obligate dimers in regulating neurotransmission and synaptic plasticity in the central nervous system. The mGluR1 subtype has been shown to be modulated by the membrane lipid environment, particularly cholesterol, though the molecular mechanisms remain elusive. In this study, we employed all-atom molecular dynamics simulations to investigate the effects of cholesterol on the conformational dynamics of the mGluR1 seven-transmembrane (7TM) domain in an inactive state model. Simulations were performed with three different cholesterol concentrations (0%, 10%, and 25%) in a palmitoyl-oleoyl phosphatidylcholine (POPC) lipid bilayer system. Our results demonstrate that cholesterol induces conformational changes in the mGluR1 dimer more significantly than in the individual protomers. Notably, cholesterol modulates the dynamics and conformations of the TM1 and TM2 helices at the dimer interface. Interestingly, an intermediate cholesterol concentration of 10% elicits more pronounced conformational changes compared to both cholesterol-depleted (0%) and cholesterol-enriched (25%) systems. Specific electrostatic interaction unique to the 10% cholesterol system further corroborate these conformational differences. Given the high sequence conservation of the 7TM domains across mGluR subtypes, the cholesterol-dependent effects observed in mGluR1 are likely applicable to other members of this receptor family. Our findings provide atomistic insights into how cholesterol modulates the conformational landscape of mGluRs, which could impact their function and signaling mechanisms.

The transmembrane domains of GPCR dimers as targets for drug development

G-protein-coupled receptors (GPCRs) can form homodimers or heterodimers that modulate specific signal transduction pathways to regulate a wide range of physiological and pathological functions. As such, GPCR dimers are novel drug targets for disorders including depression, hypertension, diabetes, and vascular dementia. The interaction between two receptors in a GPCR dimer involves a conformational change in the transmembrane domain (TMD). It has been demonstrated that the TMD has an important role in GPCR dimer formation and stability in vitro and in vivo. Moreover, increasing evidence shows that the TMD of GPCRs affects the function of dimers. Therefore, the TMD of GPCRs is an emerging target for the development of drugs to treat diseases that involve GPCR dimerization.

Structural Characterization of Receptor-Receptor Interactions in the Allosteric Modulation of G Protein-Coupled Receptor (GPCR) Dimers

G protein-coupled receptor (GPCR) oligomerization, while contentious, continues to attract the attention of researchers. Numerous experimental investigations have validated the presence of GPCR dimers, and the relevance of dimerization in the effectuation of physiological functions intensifies the attractiveness of this concept as a potential therapeutic target. GPCRs, as a single entity, have been the main source of scrutiny for drug design objectives for multiple diseases such as cancer, inflammation, cardiac, and respiratory diseases. The existence of dimers broadens the research scope of GPCR functions, revealing new signaling pathways that can be targeted for disease pathogenesis that have not previously been reported when GPCRs were only viewed in their monomeric form. This review will highlight several aspects of GPCR dimerization, which include a summary of the structural elucidation of the allosteric modulation of class C GPCR activation offered through recent solutions to the three-dimensional, full-length structures of metabotropic glutamate receptor and γ-aminobutyric acid B receptor as well as the role of dimerization in the modification of GPCR function and allostery. With the growing influence of computational methods in the study of GPCRs, we will also be reviewing recent computational tools that have been utilized to map protein-protein interactions (PPI).

FLIM-FRET-Based Structural Characterization of a Class-A GPCR Dimer in the Cell Membrane

Class-A G protein-coupled receptors (GPCRs) are known to homo-dimerize in the membrane. Yet, methods to characterize the structure of GPCR dimer in the native environment are lacking. Accordingly, the molecular basis and functional relevance of the class-A GPCR dimerization remain unclear. Here, we present the dimeric structural model of GPR17 in the cell membrane. The dimer mainly involves transmembrane helix 5 (TM5) at the interface, with F229 in TM5, a critical residue. An F229A mutation makes GPR17 monomeric regardless of the expression level of the receptor. Monomeric mutants of GPR17 display impaired ERK1/2 activation and cannot be properly internalized upon agonist treatment. Conversely, the F229C mutant is cross-linked as a dimer and behaves like wild-type. Importantly, the GPR17 dimer structure has been modeled using sparse inter-protomer FRET distance restraints obtained from fluorescence lifetime imaging microscopy. The same approach can be applied to characterizing the interactions of other important membrane proteins in the cell.

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.

Prediction and targeting of GPCR oligomer interfaces

GPCR oligomerization has emerged as a hot topic in the GPCR field in the last years. Receptors that are part of these oligomers can influence each other's function, although it is not yet entirely understood how these interactions work. The existence of such a highly complex network of interactions between GPCRs generates the possibility of alternative targets for new therapeutic approaches. However, challenges still exist in the characterization of these complexes, especially at the interface level. Different experimental approaches, such as FRET or BRET, are usually combined to study GPCR oligomer interactions. Computational methods have been applied as a useful tool for retrieving information from GPCR sequences and the few X-ray-resolved oligomeric structures that are accessible, as well as for predicting new and trustworthy GPCR oligomeric interfaces. Machine-learning (ML) approaches have recently helped with some hindrances of other methods. By joining and evaluating multiple structure-, sequence- and co-evolution-based features on the same algorithm, it is possible to dilute the issues of particular structures and residues that arise from the experimental methodology into all-encompassing algorithms capable of accurately predict GPCR-GPCR interfaces. All these methods used as a single or a combined approach provide useful information about GPCR oligomerization and its role in GPCR function and dynamics. Altogether, we present experimental, computational and machine-learning methods used to study oligomers interfaces, as well as strategies that have been used to target these dynamic complexes.

Computational prediction of GPCR oligomerization

There has been a recent and prolific expansion in the number of GPCR crystal structures being solved: in both active and inactive forms and in complex with ligand, with G protein and with each other. Despite this, there is relatively little experimental information about the precise configuration of GPCR oligomers during these different biologically relevant states. While it may be possible to identify the experimental conditions necessary to crystallize a GPCR preferentially in a specific structural conformation, computational approaches afford a potentially more tractable means of describing the probability of formation of receptor dimers and higher order oligomers. Ensemble-based computational methods based on structurally determined dimers, coupled with a computational workflow that uses quantum mechanical methods to analyze the chemical nature of the molecular interactions at a GPCR dimer interface, will generate the reproducible and accurate predictions needed to predict previously unidentified GPCR dimers and to inform future advances in our ability to understand and begin to precisely manipulate GPCR oligomers in biological systems. It may also provide information needed to achieve an increase in the number of experimentally determined oligomeric GPCR structures.

Effects of Pro1266Leu mutation on structure and function of glycoprotein Ib binding domain of von Willebrand factor

Glycoprotein Ibα (GpIbα) binding ability of A1 domain of von Willebrand factor (vWF) facilitates platelet adhesion that plays a crucial role in maintaining hemostasis and thrombosis at the site of vascular damage. There are both "loss as well as gain of function" mutations observed in this domain. Naturally occurring "gain of function" mutations leave self-activating impacts on the A1 domain which turns the normal binding to characteristic constitutive binding with GPIbα. These "gain of function" mutations are associated with the von Willebrand disease type 2B. In recent years, studies focused on understanding the mechanism and conformational patterns attached to these phenomena have been conducted, but the conformational pathways leading to such binding patterns are poorly understood as of now. To obtain a microscopic picture of such events for the better understanding of pathways, we used molecular dynamics (MD) simulations along with principal component analysis and normal mode analysis to study the effects of Pro1266Leu (Pro503Leu in structural context) mutation on the structure and function of A1 domain of vWF. MD simulations have provided atomic-level details of intermolecular motions as a function of time to understand the dynamic behavior of A1 domain of vWF. Comparative analysis of the trajectories obtained from MD simulations of both the wild type and Pro503Leu mutant suggesting appreciable conformational changes in the structure of mutant which might provide a basis for assuming the "gain of function" effects of these mutations on the A1 domain of vWF, resulting in the constitutive binding with GpIbα.

The Endocannabinoid System as a Target in Cancer Diseases: Are We There Yet?

The endocannabinoid system (ECS) has been placed in the anti-cancer spotlight in the last decade. The immense data load published on its dual role in both tumorigenesis and inhibition of tumor growth and metastatic spread has transformed the cannabinoid receptors CB1 (CB1R) and CB2 (CB2R), and other members of the endocannabinoid-like system, into attractive new targets for the treatment of various cancer subtypes. Although the clinical use of cannabinoids has been extensively documented in the palliative setting, clinical trials on their application as anti-cancer drugs are still ongoing. As drug repurposing is significantly faster and more economical than de novo introduction of a new drug into the clinic, there is hope that the existing pharmacokinetic and safety data on the ECS ligands will contribute to their successful translation into oncological healthcare. CB1R and CB2R are members of a large family of membrane proteins called G protein-coupled receptors (GPCR). GPCRs can form homodimers, heterodimers and higher order oligomers with other GPCRs or non-GPCRs. Currently, several CB1R and CB2R-containing heteromers have been reported and, in cancer cells, CB2R form heteromers with the G protein-coupled chemokine receptor CXCR4, the G protein-coupled receptor 55 (GPR55) and the tyrosine kinase receptor (TKR) human V-Erb-B2 Avian Erythroblastic Leukemia Viral Oncogene Homolog 2 (HER2). These protein complexes possess unique pharmacological and signaling properties, and their modulation might affect the antitumoral activity of the ECS. This review will explore the potential of the endocannabinoid network in the anti-cancer setting as well as the clinical and ethical pitfalls behind it, and will develop on the value of cannabinoid receptor heteromers as potential new targets for anti-cancer therapies and as prognostic biomarkers.

Multiscale simulations on human Frizzled and Taste2 GPCRs

Recently, molecular dynamics simulations, from all atom and coarse grained to hybrid methods bridging the two scales, have provided exciting functional insights into class F (Frizzled and Taste2) GPCRs (about 40 members in humans). Findings include: (i) The activation of one member of the Frizzled receptors (FZD4) involves a bending of transmembrane helix TM7 far larger than that in class A GPCRs. (ii) The affinity of an anticancer drug targeting another member (Smoothened receptor) decreases in a specific drug-resistant variant, because the mutation ultimately disrupts the binding cavity and affects TM6. (iii) A novel two-state recognition mechanism explains the very large agonist diversity for at least one member of the Taste2 GPCRs, hTAS2R46.

Interaction of POPC, DPPC, and POPE with the μ opioid receptor: A coarse-grained molecular dynamics study

The μ opioid receptor (μOR), which is part of the G protein-coupled receptors family, is a membrane protein that is modulated by its lipid environment. In the present work, we model μOR in three different membrane systems: POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine), POPE (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine), and DPPC (1, 2-dipalmitoyl-sn-glycero-3-phosphocholine) through 45 μs molecular dynamics (MD) simulations at the coarse-grained level. Our theoretical studies provide new insights to the lipid-induced modulation of the receptor. Particularly, to characterize how μOR interacts with each lipid, we analyze the tilt of the protein, the number of contacts occurring between the lipids and each amino acid of the receptor, and the μOR-lipid interface described as a network graph. We also analyze the variations in the number and the nature of the protein contacts that are induced by the lipid structure. We show that POPC interacts preferentially with helix 1 (H1) and helices H5-H6, POPE, with H5-H6 and H6-H7, and DPPC, with H4 and H6. We demonstrate how each of the three lipids shape the structure of the μOR.

The heterotetrameric structure of the adenosine A1-dopamine D1 receptor complex: Pharmacological implication for restless legs syndrome

Dopaminergic and purinergic signaling play a pivotal role in neurological diseases associated with motor symptoms, including Parkinson's disease (PD), multiple sclerosis, amyotrophic lateral sclerosis, Huntington disease, Restless Legs Syndrome (RLS), spinal cord injury (SCI), and ataxias. Extracellular dopamine and adenosine exert their functions interacting with specific dopamine (DR) or adenosine (AR) receptors, respectively, expressed on the surface of target cells. These receptors are members of the family A of G protein-coupled receptors (GPCRs), which is the largest protein superfamily in mammalian genomes. GPCRs are target of about 40% of all current marketed drugs, highlighting their importance in clinical medicine. The striatum receives the densest dopamine innervations and contains the highest density of dopamine receptors. The modulatory role of adenosine on dopaminergic transmission depends largely on the existence of antagonistic interactions mediated by specific subtypes of DRs and ARs, the so-called A_2AR-D₂R and A₁R-D₁R interactions. Due to the dopamine/adenosine antagonism in the CNS, it was proposed that ARs and DRs could form heteromers in the neuronal cell surface. Therefore, adenosine can affect dopaminergic signaling through receptor-receptor interactions and by modulations in their shared intracellular pathways in the striatum and spinal cord. In this work we describe the allosteric modulations between GPCR protomers, focusing in those of adenosine and dopamine within the A₁R-D₁R heteromeric complex, which is involved in RLS. We also propose that the knowledge about the intricate allosteric interactions within the A₁R-D₁R heterotetramer, may facilitate the treatment of motor alterations, not only when the dopamine pathway is hyperactivated (RLS, chorea, etc.) but also when motor function is decreased (SCI, aging, PD, etc.).

Signal Transduction and Pathogenic Modifications at the Melanocortin-4 Receptor: A Structural Perspective

The melanocortin-4 receptor (MC4R) can be endogenously activated by binding of melanocyte-stimulating hormones (MSH), which mediates anorexigenic effects. In contrast, the agouti-related peptide (AgRP) acts as an endogenous inverse agonist and suppresses ligand-independent basal signaling activity (orexigenic effects). Binding of ligands to MC4R leads to the activation of different G-protein subtypes or arrestin and concomitant signaling pathways. This receptor is a key protein in the hypothalamic regulation of food intake and energy expenditure and naturally-occurring inactivating MC4R variants are the most frequent cause of monogenic obesity. In general, obesity is a growing problem on a global scale and is of social, medical, and economic relevance. A significant goal is to develop optimized pharmacological tools targeting MC4R without adverse effects. To date, this has not been achieved because of inter alia non-selective ligands across the five functionally different MCR subtypes (MC1-5R). This motivates further investigation of (i) the three-dimensional MC4R structure, (ii) binding mechanisms of various ligands, and (iii) the molecular transfer process of signal transduction, with the aim of understanding how structural features are linked with functional-physiological aspects. Unfortunately, experimentally elucidated structural information is not yet available for the MC receptors, a group of class A G-protein coupled receptors (GPCRs). We, therefore, generated MC4R homology models and complexes with interacting partners to describe approximate structural properties associated with signaling mechanisms. In addition, molecular insights from pathogenic mutations were incorporated to discriminate more precisely their individual malfunction of the signal transfer mechanism.

Identifying G protein-coupled receptor dimers from crystal packings

Dimers of G protein-coupled receptors (GPCRs) are believed to be important for signaling with their associated G proteins. Low-resolution electron microscopy has shown rhodopsin dimers in native retinal membranes, and CXCR4 dimers have been found in several different crystal structures. Evidence for dimers of other GPCRs is more indirect. An alternative to computational modeling studies is to search for parallel dimers in the packing environments of the reported crystal structures of GPCRs. Two major structural types of GPCR dimers exist (as predicted by others), but there is considerable structural variation within each cluster. The different structural variants described here might reflect different functional properties and should provide a range of model structures for computational and experimental examination.

Oligomerization and cooperativity in GPCRs from the perspective of the angiotensin AT1 and dopamine D2 receptors

G Protein-Coupled Receptors (GPCRs) can form homo- and heterodimers or constitute higher oligomeric clusters with other heptahelical GPCRs. In this article, multiscale molecular modeling approaches as well as experimental techniques which are used to study oligomerization of GPCRs are reviewed. In particular, the effect of dimerization/oligomerization to the ligand binding affinity of individual protomers and also on the efficacy of the oligomer are discussed by including diverse examples from the literature. In addition, possible allosteric effects that may emerge upon interaction of GPCRs with membrane components, like cholesterol, is also discussed. Investigation of these above-mentioned interactions may greatly contribute to the candidate molecule screening studies and development of novel therapeutics with fewer adverse effects.

Probing the cooperative mechanism of the μ–δ opioid receptor heterodimer by multiscale simulation

The activation-cooperativity of the μ–δ opioid receptor heterodimer was probed by multiscale simulation coupled with a protein structure network.

Signaling within Allosteric Machines: Signal Transmission Pathways Inside G Protein-Coupled Receptors

In recent years, our understanding of function of G protein-coupled receptors (GPCRs) has changed from a picture of simple signal relays, transmitting only a particular signal to a particular G protein heterotrimer, to versatile machines, capable of various responses to different stimuli and being modulated by various factors. Some recent reports provide not only the data on ligands/modulators and resultant signals induced by them, but also deeper insights into exact pathways of signal migration and mechanisms of signal transmission through receptor structure. Combination of these computational and experimental data sheds more light on underlying mechanisms of signal transmission and signaling bias in GPCRs. In this review we focus on available clues on allosteric pathways responsible for complex signal processing within GPCRs structures, with particular emphasis on linking compatible in silico- and in vitro-derived data on the most probable allosteric connections.

Structural-Functional Features of the Thyrotropin Receptor: A Class A G-Protein-Coupled Receptor at Work

The thyroid-stimulating hormone receptor (TSHR) is a member of the glycoprotein hormone receptors, a sub-group of class A G-protein-coupled receptors (GPCRs). TSHR and its endogenous ligand thyrotropin (TSH) are of essential importance for growth and function of the thyroid gland and proper function of the TSH/TSHR system is pivotal for production and release of thyroid hormones. This receptor is also important with respect to pathophysiology, such as autoimmune (including ophthalmopathy) or non-autoimmune thyroid dysfunctions and cancer development. Pharmacological interventions directly targeting the TSHR should provide benefits to disease treatment compared to currently available therapies of dysfunctions associated with the TSHR or the thyroid gland. Upon TSHR activation, the molecular events conveying conformational changes from the extra- to the intracellular side of the cell across the membrane comprise reception, conversion, and amplification of the signal. These steps are highly dependent on structural features of this receptor and its intermolecular interaction partners, e.g., TSH, antibodies, small molecules, G-proteins, or arrestin. For better understanding of signal transduction, pathogenic mechanisms such as autoantibody action and mutational modifications or for developing new pharmacological strategies, it is essential to combine available structural data with functional information to generate homology models of the entire receptor. Although so far these insights are fragmental, in the past few decades essential contributions have been made to investigate in-depth the involved determinants, such as by structure determination via X-ray crystallography. This review summarizes available knowledge (as of December 2016) concerning the TSHR protein structure, associated functional aspects, and based on these insights we suggest several receptor complex models. Moreover, distinct TSHR properties will be highlighted in comparison to other class A GPCRs to understand the molecular activation mechanisms of this receptor comprehensively. Finally, limitations of current knowledge and lack of information are discussed highlighting the need for intensified efforts toward TSHR structure elucidation.

Interplay of G Protein-Coupled Receptors with the Membrane: Insights from Supra-Atomic Coarse Grain Molecular Dynamics Simulations

G protein-coupled receptors (GPCRs) are central to many fundamental cellular signaling pathways. They transduce signals from the outside to the inside of cells in physiological processes ranging from vision to immune response. It is extremely challenging to look at them individually using conventional experimental techniques. Recently, a pseudo atomistic molecular model has emerged as a valuable tool to access information on GPCRs, more specifically on their interactions with their environment in their native cell membrane and the consequences on their supramolecular organization. This approach uses the Martini coarse grain (CG) model to describe the receptors, lipids, and solvent in molecular dynamics (MD) simulations and in enough detail to allow conserving the chemical specificity of the different molecules. The elimination of unnecessary degrees of freedom has opened up large-scale simulations of the lipid-mediated supramolecular organization of GPCRs. Here, after introducing the Martini CGMD method, we review these studies carried out on various members of the GPCR family, including rhodopsin (visual receptor), opioid receptors, adrenergic receptors, adenosine receptors, dopamine receptor, and sphingosine 1-phosphate receptor. These studies have brought to light an interesting set of novel biophysical principles. The insights range from revealing localized and heterogeneous deformations of the membrane bilayer at the surface of the protein, specific interactions of lipid molecules with individual GPCRs, to the effect of the membrane matrix on global GPCR self-assembly. The review ends with an overview of the lessons learned from the use of the CGMD method, the biophysical-chemical findings on lipid-protein interplay.