Marta Filizola

G Protein–Coupled Receptor Oligomerization Revisited: Functional and Pharmacological Perspectives

Most evidence indicates that, as for family C G protein–coupled receptors (GPCRs), family A GPCRs form homo- and heteromers. Homodimers seem to be a predominant species, with potential dynamic formation of higher-order oligomers, particularly tetramers. Although monomeric GPCRs can activate G proteins, the pentameric structure constituted by one GPCR homodimer and one heterotrimeric G protein may provide a main functional unit, and oligomeric entities can be viewed as multiples of dimers. It still needs to be resolved if GPCR heteromers are preferentially heterodimers or if they are mostly constituted by heteromers of homodimers. Allosteric mechanisms determine a multiplicity of possible unique pharmacological properties of GPCR homomers and heteromers. Some general mechanisms seem to apply, particularly at the level of ligand-binding properties. In the frame of the dimer-cooperativity model, the two-state dimer model provides the most practical method to analyze ligand–GPCR interactions when considering receptor homomers. In addition to ligand-binding properties, unique properties for each GPCR oligomer emerge in relation to different intrinsic efficacy of ligands for different signaling pathways (functional selectivity). This gives a rationale for the use of GPCR oligomers, and particularly heteromers, as novel targets for drug development. Herein, we review the functional and pharmacological properties of GPCR oligomers and provide some guidelines for the application of discrete direct screening and high-throughput screening approaches to the discovery of receptor-heteromer selective compounds.

Making Structural Sense of Dimerization Interfaces of Delta Opioid Receptor Homodimers

Opioid receptors, like other members of the G protein-coupled receptor (GPCR) family, have been shown to associate to form dimers and/or oligomers at the plasma membrane. Whether this association is stable or transient is not known. Recent compelling evidence suggests that at least some GPCRs rapidly associate and dissociate. We have recently calculated binding affinities from free energy estimates to predict transient association between mouse delta opioid receptor (DOR) protomers at a symmetric interface involving the fourth transmembrane (TM4) helix (herein termed “4” dimer). Here we present disulfide cross-linking experiments with DOR constructs with cysteines substituted at the extracellular ends of TM4 or TM5 that confirm the formation of DOR complexes involving these helices. Our results are consistent with the involvement of TM4 and/or TM5 at the DOR homodimer interface, but possibly with differing association propensities. Coarse-grained (CG) well-tempered metadynamics simulations of two different dimeric arrangements of DOR involving TM4 alone or with TM5 (herein termed “4/5” dimer) in an explicit lipid−water environment confirmed the presence of two structurally and energetically similar configurations of the 4 dimer, as previously assessed by umbrella sampling calculations, and revealed a single energetic minimum of the 4/5 dimer. Additional CG umbrella sampling simulations of the 4/5 dimer indicated that the strength of association between DOR protomers varies depending on the protein region at the interface, with the 4 dimer being more stable than the 4/5 dimer.

Mechanistic insights into the allosteric modulation of opioid receptors by sodium ions.

The idea of sodium ions altering G-protein-coupled receptor (GPCR) ligand binding and signaling was first suggested for opioid receptors (ORs) in the 1970s and subsequently extended to other GPCRs. Recently published ultra-high-resolution crystal structures of GPCRs, including that of the δ-OR subtype, have started to shed light on the mechanism underlying sodium control in GPCR signaling by revealing details of the sodium binding site. Whether sodium accesses different receptor subtypes from the extra- or intracellular sides, following similar or different pathways, is still an open question. Earlier experiments in brain homogenates suggested a differential sodium regulation of ligand binding to the three major OR subtypes, in spite of their high degree of sequence similarity. Intrigued by this possibility, we explored the dynamic nature of sodium binding to δ-OR, μ-OR, and κ-OR by means of microsecond-scale, all-atom molecular dynamics (MD) simulations. Rapid sodium permeation was observed exclusively from the extracellular milieu, and following similar binding pathways in all three ligand-free OR systems, notwithstanding extra densities of sodium observed near nonconserved residues of κ-OR and δ-OR, but not in μ-OR. We speculate that these differences may be responsible for the differential increase in antagonist binding affinity of μ-OR by sodium resulting from specific ligand binding experiments in transfected cells. On the other hand, sodium reduced the level of binding of subtype-specific agonists to all OR subtypes. Additional biased and unbiased MD simulations were conducted using the δ-OR ultra-high-resolution crystal structure as a model system to provide a mechanistic explanation for this experimental observation.

Preferred supramolecular organization and dimer interfaces of opioid receptors from simulated self-association.

Substantial evidence in support of the formation of opioid receptor (OR) di-/oligomers suggests previously unknown mechanisms used by these proteins to exert their biological functions. In an attempt to guide experimental assessment of the identity of the minimal signaling unit for ORs, we conducted extensive coarse-grained (CG) molecular dynamics (MD) simulations of different combinations of the three major OR subtypes, i.e., μ-OR, δ-OR, and κ-OR, in an explicit lipid bilayer. Specifically, we ran multiple, independent MD simulations of each homomeric μ-OR/μ-OR, δ-OR/δ-OR, and κ-OR/κ-OR complex, as well as two of the most studied heteromeric complexes, i.e., δ-OR/μ-OR and δ-OR/κ-OR, to derive the preferred supramolecular organization and dimer interfaces of ORs in a cell membrane model. These simulations yielded over 250 microseconds of accumulated data, which correspond to approximately 1 millisecond of effective simulated dynamics according to established scaling factors of the CG model we employed. Analysis of these data indicates similar preferred supramolecular organization and dimer interfaces of ORs across the different receptor subtypes, but also important differences in the kinetics of receptor association at specific dimer interfaces. We also investigated the kinetic properties of interfacial lipids, and explored their possible role in modulating the rate of receptor association and in promoting the formation of filiform aggregates, thus supporting a distinctive role of the membrane in OR oligomerization and, possibly, signaling.

Synthetic and Receptor Signaling Explorations of the Mitragyna Alkaloids: Mitragynine as an Atypical Molecular Framework for Opioid Receptor Modulators.

Mu-opioid receptor agonists represent mainstays of pain management. However, the therapeutic use of these agents is associated with serious side effects, including potentially lethal respiratory depression. Accordingly, there is a longstanding interest in the development of new opioid analgesics with improved therapeutic profiles. The alkaloids of the Southeast Asian plant Mitragyna speciosa, represented by the prototypical member mitragynine, are an unusual class of opioid receptor modulators with distinct pharmacological properties. Here we describe the first receptor-level functional characterization of mitragynine and related natural alkaloids at the human mu-, kappa-, and delta-opioid receptors. These results show that mitragynine and the oxidized analogue 7-hydroxymitragynine, are partial agonists of the human mu-opioid receptor and competitive antagonists at the kappa- and delta-opioid receptors. We also show that mitragynine and 7-hydroxymitragynine are G-protein-biased agonists of the mu-opioid receptor, which do not recruit β-arrestin following receptor activation. Therefore, the Mitragyna alkaloid scaffold represents a novel framework for the development of functionally biased opioid modulators, which may exhibit improved therapeutic profiles. Also presented is an enantioselective total synthesis of both (-)-mitragynine and its unnatural enantiomer, (+)-mitragynine, employing a proline-catalyzed Mannich-Michael reaction sequence as the key transformation. Pharmacological evaluation of (+)-mitragynine revealed its much weaker opioid activity. Likewise, the intermediates and chemical transformations developed in the total synthesis allowed the elucidation of previously unexplored structure-activity relationships (SAR) within the Mitragyna scaffold. Molecular docking studies, in combination with the observed chemical SAR, suggest that Mitragyna alkaloids adopt a binding pose at the mu-opioid receptor that is distinct from that of classical opioids.

Molecular determinants of non-specific recognition of δ, μ, and κ opioid receptors

Abstract Identification of the molecular determinants of recognition common to all three opioid receptors embedded in a single three-dimensional (3D) non-specific recognition pharmacophore has been carried out. The working hypothesis that underlies the computational study reported here is that ligands that bind with significant affinity to all three cloned opioid receptors, δ, μ, and κ, but with different combinations of activation and inhibition properties at these receptors, could be promising behaviorally selective analgesics with diminished side effects. The study presented here represents the first step towards the rational design of such therapeutic agents. The common 3D pharmacophore developed for recognition of δ, μ, and κ opioid receptors was based on the receptor affinities determined for 23 different opioid ligands that display no specificity for any of the receptor subtypes. The pharmacophore centers identified are a protonated amine, two hydrophobic groups, and the centroid of an aromatic group in a geometric arrangement common to all 23, non-specific, opioid ligands studied. Using this three-dimensional pharmacophore as a query for searching 3D structural databases, novel compounds potentially involved in non-specific recognition of δ, μ, and κ opioid receptors were retrieved. These compounds can be valuable candidates for novel behaviorally selective analgesics with diminished or no side effects, and thus with potential therapeutic usefulness.

Discovery, synthesis, and molecular pharmacology of selective positive allosteric modulators of the δ-opioid receptor

Allosteric modulators of G protein-coupled receptors (GPCRs) have a number of potential advantages compared to agonists or antagonists that bind to the orthosteric site of the receptor. These include the potential for receptor selectivity, maintenance of the temporal and spatial fidelity of signaling in vivo, the ceiling effect of the allosteric cooperativity which may prevent overdose issues, and engendering bias by differentially modulating distinct signaling pathways. Here we describe the discovery, synthesis, and molecular pharmacology of δ-opioid receptor-selective positive allosteric modulators (δ PAMs). These δ PAMs increase the affinity and/or efficacy of the orthosteric agonists leu-enkephalin, SNC80 and TAN67, as measured by receptor binding, G protein activation, β-arrestin recruitment, adenylyl cyclase inhibition, and extracellular signal-regulated kinases (ERK) activation. As such, these compounds are useful pharmacological tools to probe the molecular pharmacology of the δ receptor and to explore the therapeutic potential of δ PAMs in diseases such as chronic pain and depression.

αIIbβ3 variants defined by next-generation sequencing: Predicting variants likely to cause Glanzmann thrombasthenia

Significance Next-generation sequencing is identifying millions of novel gene variants, presenting challenges to researchers and clinicians. Variations in the genes ITGA2B and ITGB3 affect integrin αIIbβ3, leading to the bleeding disorder Glanzmann thrombasthenia. We analyzed novel missense variants on ∼32,000 alleles of ITGA2B and ITGB3 and found missense variants affecting ∼10% of the amino acids in each protein in ∼1.3% of the population. Almost all variants are rare, indicating recent entry into the population. Two novel variants we predicted would be deleterious profoundly affected recombinant protein expression. At cut-off values that correctly predicted at least 69% of the known Glanzmann thrombasthenia mutations as deleterious, three variant prediction algorithms predicted that at least 27% of the novel variants are deleterious. Next-generation sequencing is transforming our understanding of human genetic variation but assessing the functional impact of novel variants presents challenges. We analyzed missense variants in the integrin αIIbβ3 receptor subunit genes ITGA2B and ITGB3 identified by whole-exome or -genome sequencing in the ThromboGenomics project, comprising ∼32,000 alleles from 16,108 individuals. We analyzed the results in comparison with 111 missense variants in these genes previously reported as being associated with Glanzmann thrombasthenia (GT), 20 associated with alloimmune thrombocytopenia, and 5 associated with aniso/macrothrombocytopenia. We identified 114 novel missense variants in ITGA2B (affecting ∼11% of the amino acids) and 68 novel missense variants in ITGB3 (affecting ∼9% of the amino acids). Of the variants, 96% had minor allele frequencies (MAF) < 0.1%, indicating their rarity. Based on sequence conservation, MAF, and location on a complete model of αIIbβ3, we selected three novel variants that affect amino acids previously associated with GT for expression in HEK293 cells. αIIb P176H and β3 C547G severely reduced αIIbβ3 expression, whereas αIIb P943A partially reduced αIIbβ3 expression and had no effect on fibrinogen binding. We used receiver operating characteristic curves of combined annotation-dependent depletion, Polyphen 2-HDIV, and sorting intolerant from tolerant to estimate the percentage of novel variants likely to be deleterious. At optimal cut-off values, which had 69–98% sensitivity in detecting GT mutations, between 27% and 71% of the novel αIIb or β3 missense variants were predicted to be deleterious. Our data have implications for understanding the evolutionary pressure on αIIbβ3 and highlight the challenges in predicting the clinical significance of novel missense variants.

Opioid receptors: Structural and mechanistic insights into pharmacology and signaling

Opioid receptors are important drug targets for pain management, addiction, and mood disorders. Although substantial research on these important subtypes of G protein-coupled receptors has been conducted over the past two decades to discover ligands with higher specificity and diminished side effects, currently used opioid therapeutics remain suboptimal. Luckily, recent advances in structural biology of opioid receptors provide unprecedented insights into opioid receptor pharmacology and signaling. We review here a few recent studies that have used the crystal structures of opioid receptors as a basis for revealing mechanistic details of signal transduction mediated by these receptors, and for the purpose of drug discovery.

Differential stability of the crystallographic interfaces of mu- and kappa-opioid receptors.

The recent mu-opioid receptor (MOPr) and kappa-opioid receptor (KOPr) crystal structures have inspired hypotheses of physiologically relevant dimerization contacts, specifically: a closely packed interface involving transmembrane (TM) helices TM5 and TM6, and a less compact interface, involving TM1, TM2, and helix 8 (H8). While the former was only found in MOPr crystals, similar arrangements of the latter were identified for both KOPr and MOPr. The relevance of these interfaces outside of a crystal lattice is called into question by the possibility that they might be influenced by the specific crystallization conditions. In this study, we have employed umbrella sampling molecular dynamics simulations of coarse-grained representations of the interacting MOPr or KOPr crystallographic structures, in the absence of the T4 lysozyme, and in an explicit lipid-water environment, to determine the strength of receptor dimerization at the different crystallographic interfaces. We note that the shape of the interface plays a dominant role in the strength of the interaction, and the pattern of contacting residues defines the shape of the potential of mean force. This information can be used to guide experiments aimed at exploring the role of dimerization in opioid receptor function.