Computational approaches for modeling GPCR dimerization

Integration and Spatial Organization of Signaling by G Protein-Coupled Receptor Homo- and Heterodimers

Information flow from a source to a receiver becomes informative when the recipient can process the signal into a meaningful form. Information exchange and interpretation is essential in biology and understanding how cells integrate signals from a variety of information-coding molecules into complex orchestrated responses is a major challenge for modern cell biology. In complex organisms, cell to cell communication occurs mostly through neurotransmitters and hormones, and receptors are responsible for signal recognition at the membrane level and information transduction inside the cell. The G protein-coupled receptors (GPCRs) are the largest family of membrane receptors, with nearly 800 genes coding for these proteins. The recognition that GPCRs may physically interact with each other has led to the hypothesis that their dimeric state can provide the framework for temporal coincidence in signaling pathways. Furthermore, the formation of GPCRs higher order oligomers provides the structural basis for organizing distinct cell compartments along the plasma membrane where confined increases in second messengers may be perceived and discriminated. Here, we summarize evidence that supports these conjectures, fostering new ideas about the physiological role played by receptor homo- and hetero-oligomerization in cell biology.

Structure-Based Design and Development of Chemical Probes Targeting Putative MOR-CCR5 Heterodimers to Inhibit Opioid Exacerbated HIV-1 Infectivity

Crystal structures of ligand-bound G-protein-coupled receptors provide tangible templates for rationally designing molecular probes. Herein, we report the structure-based design, chemical synthesis, and biological investigations of bivalent ligands targeting putative mu opioid receptor C-C motif chemokine ligand 5 (MOR-CCR5) heterodimers. The bivalent ligand VZMC013 possessed nanomolar level binding affinities for both the MOR and CCR5, inhibited CCL5-stimulated calcium mobilization, and remarkably improved anti-HIV-1BaL activity over previously reported bivalent ligands. VZMC013 inhibited viral infection in TZM-bl cells coexpressing CCR5 and MOR to a greater degree than cells expressing CCR5 alone. Furthermore, VZMC013 blocked human immunodeficiency virus (HIV)-1 entry in peripheral blood mononuclear cells (PBMC) cells in a concentration-dependent manner and inhibited opioid-accelerated HIV-1 entry more effectively in phytohemagglutinin-stimulated PBMC cells than in the absence of opioids. A three-dimensional molecular model of VZMC013 binding to the MOR-CCR5 heterodimer complex is constructed to elucidate its mechanism of action. VZMC013 is a potent chemical probe targeting MOR-CCR5 heterodimers and may serve as a pharmacological agent to inhibit opioid-exacerbated HIV-1 entry.

Recent advances in dopaminergic strategies for the treatment of Parkinson's disease

Parkinson's disease (PD) is the second most common progressive neurodegenerative disease worldwide. However, there is no available therapy reversing the neurodegenerative process of PD. Based on the loss of dopamine or dopaminergic dysfunction in PD patients, most of the current therapies focus on symptomatic relief to improve patient quality of life. As dopamine replacement treatment remains the most effective symptomatic pharmacotherapy for PD, herein we provide an overview of the current pharmacotherapies, summarize the clinical development status of novel dopaminergic agents, and highlight the challenge and opportunity of emerging preclinical dopaminergic approaches aimed at managing the features and progression of PD.

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.

The effect of dimerization on the activation and conformational dynamics of adenosine A1 receptor

The adenosine A1 receptor (A1R) is one of four adenosine receptors in humans, which are involved in the function of the cardiovascular, respiratory and central nervous systems. Experimental results indicate that A1R can form a homodimer and that the protomer-protomer interaction in the A1R dimer is related to certain pharmacological characteristics of A1R activation. In this work, we performed docking, metadynamics simulation, conventional molecular dynamics simulations, Gaussian-accelerated molecular dynamics simulations, potential of mean force calculations, dynamic cross-correlation motions analysis and community network analysis to study the binding mode of 5'-N-ethylcarboxamidoadenosine (NECA) to A1R and the effect of dimerization on the activation of A1R. Our results show that NECA binds to A1R in a similar mode to adenosine in the A1R crystal structure and NECA in the A2AR crystal structure. The A1R homodimer can be activated by one or two agonists with NECA occupying its orthosteric pockets in one (which we call the NECA-A1R system) or both protomers (which we call the dNECA-A1R system). In the NECA-A1R system, activation is predicated in the protomer without NECA bound. In the dNECA-A1R system, only one protomer achieves the active state. These findings suggest an asymmetrical activation mechanism of the homodimer and a negative cooperativity between the two protomers. We envision that our results may further facilitate the drug development of A1R.

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.

Computational Modeling of Realistic Cell Membranes

Cell membranes contain a large variety of lipid types and are crowded with proteins, endowing them with the plasticity needed to fulfill their key roles in cell functioning. The compositional complexity of cellular membranes gives rise to a heterogeneous lateral organization, which is still poorly understood. Computational models, in particular molecular dynamics simulations and related techniques, have provided important insight into the organizational principles of cell membranes over the past decades. Now, we are witnessing a transition from simulations of simpler membrane models to multicomponent systems, culminating in realistic models of an increasing variety of cell types and organelles. Here, we review the state of the art in the field of realistic membrane simulations and discuss the current limitations and challenges ahead.

1988-2018: Thirty years of drug smuggling at the nano scale. Challenges and opportunities of cell-penetrating peptides in biomedical research

In 1988, two unrelated papers reported the discovery of peptide vectors with innate cell translocation properties, setting the ground for a new area of research that over the years has grown into considerable therapeutic potential. The vectors, named cell-penetrating peptides (CPPs), constitute a now large and diversified family, sharing the extraordinary ability to diffuse unaltered across cell membranes while ferrying diverse associated cargos. Such properties have made CPPs ideal tools for delivery of nucleic acids, proteins and other therapeutic/diagnostic molecules to cells and tissues via covalent conjugation or complexation. This year 2018 marks the 30^th anniversary of a peptide research landmark opening new perspectives in drug delivery. Given its vastness, exhaustive coverage of the main features and accomplishments in the CPP field is virtually impossible. Hence this manuscript, after saluting the above 30^th jubilee, focuses by necessity on the most recent contributions, providing a comprehensive list of recognized CPPs and their latest-reported applications over the last two years. In addition, it thoroughly reviews three areas of peptide vector research of particular interest to us, namely (i) efficient transport of low-bioavailability drugs into the brain; (ii) CPP-delivered disruptors of G protein-coupled receptor (GPCRs) heteromers related to several disorders, and (iii) CPP-mediated delivery of useful but poorly internalized drugs into parasites.

Mapping the Interface of a GPCR Dimer: A Structural Model of the A2A Adenosine and D2 Dopamine Receptor Heteromer

The A_2A adenosine (A_2AR) and D₂ dopamine (D₂R) receptors form oligomers in the cell membrane and allosteric interactions across the A_2AR–D₂R heteromer represent a target for development of drugs against central nervous system disorders. However, understanding of the molecular determinants of A_2AR–D₂R heteromerization and the allosteric antagonistic interactions between the receptor protomers is still limited. In this work, a structural model of the A_2AR–D₂R heterodimer was generated using a combined experimental and computational approach. Regions involved in the heteromer interface were modeled based on the effects of peptides derived from the transmembrane (TM) helices on A_2AR–D₂R receptor–receptor interactions in bioluminescence resonance energy transfer (BRET) and proximity ligation assays. Peptides corresponding to TM-IV and TM-V of the A_2AR blocked heterodimer interactions and disrupted the allosteric effect of A_2AR activation on D₂R agonist binding. Protein–protein docking was used to construct a model of the A_2AR–D₂R heterodimer with a TM-IV/V interface, which was refined using molecular dynamics simulations. Mutations in the predicted interface reduced A_2AR–D₂R interactions in BRET experiments and altered the allosteric modulation. The heterodimer model provided insights into the structural basis of allosteric modulation and the technique developed to characterize the A_2AR–D₂R interface can be extended to study the many other G protein-coupled receptors that engage in heteroreceptor complexes.

Ghrelin Through GHSR1a and OX1R Heterodimers Reveals a Gαs-cAMP-cAMP Response Element Binding Protein Signaling Pathway in Vitro

Growth hormone secretagogue receptor 1α (GHSR1a) and Orexin 1 receptor (OX1R) are involved in various important physiological processes, and have many similar characteristics in function and distribution in peripheral tissues and the central nervous system. We explored the possibility of heterodimerization between GHSR1a and OX1R and revealed a signal transduction pathway mechanism. In this study, bioluminescence and fluorescence resonance energy transfer and co-immunoprecipitation (Co-IP) analyses were performed to demonstrate the formation of functional GHSR1a/OX1R heterodimers. This showed that a peptide corresponding to the 5-transmembrane domain of OX1R impaired heterodimer construction. We found that ghrelin stimulated GHSR1a/OX1R heterodimer cells to increase the activation of Gαs protein, compared to the cells that express GHSR1a. Stimulation of GHSR1a/OX1R heterodimers with orexin-A did not alter GPCR interactions with Gα protein subunits. GHSR1a/OX1R heterodimers induced Gαs and downstream signaling pathway activity, including increase of cAMP-response element luciferase reporter activity and cAMP levels. In addition, ghrelin induced a higher proliferation rate in SH-SY5Y cells than in controls. This suggests that ghrelin GHSR1a/OX1R heterodimers promotes an upregulation of a Gαs-cAMP-cAMP-responsive element signaling pathway in vitro and an increase in neuroblastoma cell proliferation.

Methods used to study the oligomeric structure of G-protein-coupled receptors

G-protein-coupled receptors (GPCRs), which constitute the largest family of cell surface receptors, were originally thought to function as monomers, but are now recognized as being able to act in a wide range of oligomeric states and indeed, it is known that the oligomerization state of a GPCR can modulate its pharmacology and function. A number of experimental techniques have been devised to study GPCR oligomerization including those based upon traditional biochemistry such as blue-native PAGE (BN-PAGE), co-immunoprecipitation (Co-IP) and protein-fragment complementation assays (PCAs), those based upon resonance energy transfer, FRET, time-resolved FRET (TR-FRET), FRET spectrometry and bioluminescence resonance energy transfer (BRET). Those based upon microscopy such as FRAP, total internal reflection fluorescence microscopy (TIRFM), spatial intensity distribution analysis (SpIDA) and various single molecule imaging techniques. Finally with the solution of a growing number of crystal structures, X-ray crystallography must be acknowledged as an important source of discovery in this field. A different, but in many ways complementary approach to the use of more traditional experimental techniques, are those involving computational methods that possess obvious merit in the study of the dynamics of oligomer formation and function. Here, we summarize the latest developments that have been made in the methods used to study GPCR oligomerization and give an overview of their application.

Apelin receptor homodimer-oligomers revealed by single-molecule imaging and novel G protein-dependent signaling

The apelin receptor (APJ) belongs to family A of the G protein-coupled receptors (GPCRs) and is a potential pharmacotherapeutic target for heart failure, hypertension, and other cardiovascular diseases. There is evidence APJ heterodimerizes with other GPCRs; however, the existence of APJ homodimers and oligomers remains to be investigated. Here, we measured APJ monomer-homodimer-oligomer interconversion by monitoring APJ dynamically on cells and compared their proportions, spatial arrangement, and mobility using total internal reflection fluorescence microscopy, resonance energy transfer, and proximity biotinylation. In cells with <0.3 receptor particles/μm², approximately 60% of APJ molecules were present as dimers or oligomers. APJ dimers were present on the cell surface in a dynamic equilibrium with constant formation and dissociation of receptor complexes. Furthermore, we applied interference peptides and MALDI-TOF mass spectrometry to confirm APJ homo-dimer and explore the dimer-interfaces. Peptides corresponding to transmembrane domain (TMD)1, 2, 3, and 4, but not TMD5, 6, and 7, disrupted APJ dimerization. APJ mutants in TMD1 and TMD2 also decreased bioluminescence resonance energy transfer of APJ dimer. APJ dimerization resulted in novel functional characteristics, such as a distinct G-protein binding profile and cell responses after agonist stimulation. Thus, dimerization may serve as a unique mechanism for fine-tuning APJ-mediated functions.

Membrane-Mediated Oligomerization of G Protein Coupled Receptors and Its Implications for GPCR Function

The dimerization or even oligomerization of G protein coupled receptors (GPCRs) causes ongoing, controversial debates about its functional role and the coupled biophysical, biochemical or biomedical implications. A continously growing number of studies hints to a relation between oligomerization and function of GPCRs and strengthens the assumption that receptor assembly plays a key role in the regulation of protein function. Additionally, progress in the structural analysis of GPCR-G protein and GPCR-ligand interactions allows to distinguish between actively functional and non-signaling complexes. Recent findings further suggest that the surrounding membrane, i.e., its lipid composition may modulate the preferred dimerization interface and as a result the abundance of distinct dimeric conformations. In this review, the association of GPCRs and the role of the membrane in oligomerization will be discussed. An overview of the different reported oligomeric interfaces is provided and their capability for signaling discussed. The currently available data is summarized with regard to the formation of GPCR oligomers, their structures and dependency on the membrane microenvironment as well as the coupling of oligomerization to receptor function.

Allosteric signaling through an mGlu2 and 5-HT2A heteromeric receptor complex and its potential contribution to schizophrenia

Heterotrimeric guanine nucleotide-binding protein (G protein)-coupled receptors (GPCRs) can form multiprotein complexes (heteromers), which can alter the pharmacology and functions of the constituent receptors. Previous findings demonstrated that the Gq/11-coupled serotonin 5-HT2A receptor and the Gi/o-coupled metabotropic glutamate 2 (mGlu2) receptor-GPCRs that are involved in signaling alterations associated with psychosis-assemble into a heteromeric complex in the mammalian brain. In single-cell experiments with various mutant versions of the mGlu2 receptor, we showed that stimulation of cells expressing mGlu2-5-HT2A heteromers with an mGlu2 agonist led to activation of Gq/11 proteins by the 5-HT2A receptors. For this crosstalk to occur, one of the mGlu2 subunits had to couple to Gi/o proteins, and we determined the relative location of the Gi/o-contacting subunit within the mGlu2 homodimer of the heteromeric complex. Additionally, mGlu2-dependent activation of Gq/11, but not Gi/o, was reduced in the frontal cortex of 5-HT2A knockout mice and was reduced in the frontal cortex of postmortem brains from schizophrenic patients. These findings offer structural insights into this important target in molecular psychiatry.