Mechanism of μ-Opioid Receptor-Magnesium Interaction and Positive Allosteric Modulation

Functional dynamics of G protein-coupled receptors reveal new routes for drug discovery

G protein-coupled receptors (GPCRs) are the largest human membrane protein family that transduce extracellular signals into cellular responses. They are major pharmacological targets, with approximately 26% of marketed drugs targeting GPCRs, primarily at their orthosteric binding site. Despite their prominence, predicting the pharmacological effects of novel GPCR-targeting drugs remains challenging due to the complex functional dynamics of these receptors. Recent advances in X-ray crystallography, cryo-electron microscopy, spectroscopic techniques and molecular simulations have enhanced our understanding of receptor conformational dynamics and ligand interactions with GPCRs. These developments have revealed novel ligand-binding modes, mechanisms of action and druggable pockets. In this Review, we highlight such aspects for recently discovered small-molecule drugs and drug candidates targeting GPCRs, focusing on three categories: allosteric modulators, biased ligands, and bivalent and bitopic compounds. Although studies so far have largely been retrospective, integrating structural data on ligand-induced receptor functional dynamics into the drug discovery pipeline has the potential to guide the identification of drug candidates with specific abilities to modulate GPCR interactions with intracellular effector proteins such as G proteins and β-arrestins, enabling more tailored selectivity and efficacy profiles.

Positive allosteric modulation of µ-opioid receptor - A new possible approach in the pain management?

The antinociceptive effect of the opioid drugs is achieved through activation of the µ-opioid receptor (MOP). The orthosteric and allosteric sites of opioid receptors may be modulated, orthosteric site by endogenous i.e.β-endorphin and exogenous opioids (morphine, oxycodone, fentanyl); whereas BMS-986121, BMS-986122, Comp5, MS1, Ignavine or even oxytocin act on the allosteric site of the MOP. Opioid therapy is associated with numerous side effects, such as: respiratory depression, sedation, constipation, and importantly, prolonged therapy can influence the development of tolerance, overdose, and addiction. Opioid tolerance is a result of MOP internalization and desensitization, preceded by MOP phosphorylation, performed by protein kinases such as: PKA, PKC, GRKs or CaMKII. In vitro and in vivo data suggest that positive allosteric modulators may enhance antinociception triggered by orthosteric ligands and reduce side effects, which would allow the dose of opioids to be reduced and thus provide a more effective therapy. In this review, we present that positive modulation of the allosteric sites of MOP may constitute a new strategy for pain therapy.

Structural and dynamic insights into the activation of the μ-opioid receptor by an allosteric modulator

G-protein-coupled receptors (GPCRs) play pivotal roles in various physiological processes. These receptors are activated to different extents by diverse orthosteric ligands and allosteric modulators. However, the mechanisms underlying these variations in signaling activity by allosteric modulators remain largely elusive. Here, we determine the three-dimensional structure of the μ-opioid receptor (MOR), a class A GPCR, in complex with the Gi protein and an allosteric modulator, BMS-986122, using cryogenic electron microscopy. Our results reveal that BMS-986122 binding induces changes in the map densities corresponding to R1673.50 and Y2545.58, key residues in the structural motifs conserved among class A GPCRs. Nuclear magnetic resonance analyses of MOR in the absence of the Gi protein reveal that BMS-986122 binding enhances the formation of the interaction between R1673.50 and Y2545.58, thus stabilizing the fully-activated conformation, where the intracellular half of TM6 is outward-shifted to allow for interaction with the Gi protein. These findings illuminate that allosteric modulators like BMS-986122 can potentiate receptor activation through alterations in the conformational dynamics in the core region of GPCRs. Together, our results demonstrate the regulatory mechanisms of GPCRs, providing insights into the rational development of therapeutics targeting GPCRs.

Stereoselective recognition of morphine enantiomers by -opioid receptor

Stereospecific recognition of chiral molecules plays a crucial role in biological systems. The μ-opioid receptor (MOR) exhibits binding affinity towards (-)-morphine, a well-established gold standard in pain management, while it shows minimal binding affinity for the (+)-morphine enantiomer, resulting in a lack of analgesic activity. Understanding how MOR stereoselectively recognizes morphine enantiomers has remained a puzzle in neuroscience and pharmacology for over half-a-century due to the lack of direct observation techniques. To unravel this mystery, we constructed the binding and unbinding processes of morphine enantiomers with MOR via molecular dynamics simulations to investigate the thermodynamics and kinetics governing MOR's stereoselective recognition of morphine enantiomers. Our findings reveal that the binding of (-)-morphine stabilizes MOR in its activated state, exhibiting a deep energy well and a prolonged residence time. In contrast, (+)-morphine fails to sustain the activation state of MOR. Furthermore, the results suggest that specific residues, namely D1142.50 and D1473.32, are deprotonated in the active state of MOR bound to (-)-morphine. This work highlights that the selectivity in molecular recognition goes beyond binding affinities, extending into the realm of residence time.

GPU-specific algorithms for improved solute sampling in grand canonical Monte Carlo simulations

The Grand Canonical Monte Carlo (GCMC) ensemble defined by the excess chemical potential, μex , volume, and temperature, in the context of molecular simulations allows for variations in the number of particles in the system. In practice, GCMC simulations have been widely applied for the sampling of rare gasses and water, but limited in the context of larger molecules. To overcome this limitation, the oscillating μex GCMC method was introduced and shown to be of utility for sampling small solutes, such as formamide, propane, and benzene, as well as for ionic species such as monocations, acetate, and methylammonium. However, the acceptance of GCMC insertions is low, and the method is computationally demanding. In the present study, we improved the sampling efficiency of the GCMC method using known cavity-bias and configurational-bias algorithms in the context of GPU architecture. Specifically, for GCMC simulations of aqueous solution systems, the configurational-bias algorithm was extended by applying system partitioning in conjunction with a random interval extraction algorithm, thereby improving the efficiency in a highly parallel computing environment. The method is parallelized on the GPU using CUDA and OpenCL, allowing for the code to run on both Nvidia and AMD GPUs, respectively. Notably, the method is particularly well suited for GPU computing as the large number of threads allows for simultaneous sampling of a large number of configurations during insertion attempts without additional computational overhead. In addition, the partitioning scheme allows for simultaneous insertion attempts for large systems, offering considerable efficiency. Calculations on the BK Channel, a transporter, including a lipid bilayer with over 760,000 atoms, show a speed up of ~53-fold through the use of system partitioning. The improved algorithm is then combined with an enhanced μex oscillation protocol and shown to be of utility in the context of the site-identification by ligand competitive saturation (SILCS) co-solvent sampling approach as illustrated through application to the protein CDK2.

Molecular Dynamics Simulations of Ion Permeation in Human Voltage-Gated Sodium Channels

The recent determination of cryo-EM structures of voltage-gated sodium (Na_v) channels has revealed many details of these proteins. However, knowledge of ionic permeation through the Na_v pore remains limited. In this work, we performed atomistic molecular dynamics (MD) simulations to study the structural features of various neuronal Na_v channels based on homology modeling of the cryo-EM structure of the human Na_v1.4 channel and, in addition, on the recently resolved configuration for Na_v1.2. In particular, single Na⁺ permeation events during standard MD runs suggest that the ion resides in the inner part of the Na_v selectivity filter (SF). On-the-fly free energy parametrization (OTFP) temperature-accelerated molecular dynamics (TAMD) was also used to calculate two-dimensional free energy surfaces (FESs) related to single/double Na⁺ translocation through the SF of the homology-based Na_v1.2 model and the cryo-EM Na_v1.2 structure, with different realizations of the DEKA filter domain. These additional simulations revealed distinct mechanisms for single and double Na⁺ permeation through the wild-type SF, which has a charged lysine in the DEKA ring. Moreover, the configurations of the ions in the SF corresponding to the metastable states of the FESs are specific for each SF motif. Overall, the description of these mechanisms gives us new insights into ion conduction in human Na_v cryo-EM-based and cryo-EM configurations that could advance understanding of these systems and how they differ from potassium and bacterial Na_v channels.

Three-Dimensional Structural Insights Have Revealed the Distinct Binding Interactions of Agonists, Partial Agonists, and Antagonists with the µ Opioid Receptor

The United States is experiencing the most profound and devastating opioid crisis in history, with the number of deaths involving opioids, including prescription and illegal opioids, continuing to climb over the past two decades. This severe public health issue is difficult to combat as opioids remain a crucial treatment for pain, and at the same time, they are also highly addictive. Opioids act on the opioid receptor, which in turn activates its downstream signaling pathway that eventually leads to an analgesic effect. Among the four types of opioid receptors, the µ subtype is primarily responsible for the analgesic cascade. This review describes available 3D structures of the µ opioid receptor in the protein data bank and provides structural insights for the binding of agonists and antagonists to the receptor. Comparative analysis on the atomic details of the binding site in these structures was conducted and distinct binding interactions for agonists, partial agonists, and antagonists were observed. The findings in this article deepen our understanding of the ligand binding activity and shed some light on the development of novel opioid analgesics which may improve the risk benefit balance of existing opioids.

Insights into divalent cation regulation and G13-coupling of orphan receptor GPR35

Endogenous ions play important roles in the function and pharmacology of G protein-coupled receptors (GPCRs) with limited atomic evidence. In addition, compared with G protein subtypes G_s, G_i/o, and G_q/11, insufficient structural evidence is accessible to understand the coupling mechanism of G_12/13 protein by GPCRs. Orphan receptor GPR35, which is predominantly expressed in the gastrointestinal tract and is closely related to inflammatory bowel diseases (IBDs), stands out as a prototypical receptor for investigating ionic modulation and G₁₃ coupling. Here we report a cryo-electron microscopy structure of G₁₃-coupled GPR35 bound to an anti-allergic drug, lodoxamide. This structure reveals a novel divalent cation coordination site and a unique ionic regulatory mode of GPR35 and also presents a highly positively charged binding pocket and the complementary electrostatic ligand recognition mode, which explain the promiscuity of acidic ligand binding by GPR35. Structural comparison of the GPR35-G₁₃ complex with other G protein subtypes-coupled GPCRs reveals a notable movement of the C-terminus of α5 helix of the Gα₁₃ subunit towards the receptor core and the least outward displacement of the cytoplasmic end of GPR35 TM6. A featured 'methionine pocket' contributes to the G₁₃ coupling by GPR35. Together, our findings provide a structural basis for divalent cation modulation, ligand recognition, and subsequent G₁₃ protein coupling of GPR35 and offer a new opportunity for designing GPR35-targeted drugs for the treatment of IBDs.

Bias-inducing allosteric binding site in mu-opioid receptor signaling

Abstract

G-protein-biased agonism of the mu-opioid receptor (μ-OR) is emerging as a promising strategy in analgesia. A deep understanding of how biased agonists modulate and differentiate G-protein-coupled receptors (GPCR) signaling pathways and how this is transferred into the cell are open questions. Here, using extensive all-atom molecular dynamics simulations, we analyzed the binding recognition process and signaling effects of three prototype μ-OR agonists. Our suggested structural mechanism of biased signaling in μ-OR involves an allosteric sodium ion site, water networks, conformational rearrangements in conserved motifs and collective motions of loops and transmembrane helices. These analyses led us to highlight the relevance of a bias-inducing allosteric binding site in the understanding of μ-OR’s G-protein-biased signaling. These results also suggest a competitive equilibrium between the agonists and the allosteric sodium ion, where the bias-inducing allosteric binding site can be modulated by this ion or an agonist such as herkinorin. Notably, herkinorin arises as the archetype modulator of μ-OR and its interactive pattern could be used for screening efforts via protein–ligand interaction fingerprint (PLIF) studies.

Article highlights

The Effects of Sodium Ions on Ligand Binding and Conformational States of G Protein-Coupled Receptors-Insights from Mass Spectrometry

The use of mass spectrometry to investigate proteins is now well established and provides invaluable information for both soluble and membrane protein assemblies. Maintaining transient noncovalent interactions under physiological conditions, however, remains challenging. Here, using nanoscale electrospray ionization emitters, we establish conditions that enable mass spectrometry of two G protein-coupled receptors (GPCR) from buffers containing high concentrations of sodium ions. For the Class A GPCR, the adenosine 2A receptor, we observe ligand-induced changes to sodium binding of the receptor at the level of individual sodium ions. We find that antagonists promote sodium binding while agonists attenuate sodium binding. These findings are in line with high-resolution X-ray crystallography wherein only inactive conformations retain sodium ions in allosteric binding pockets. For the glucagon receptor (a Class B GPCR) we observed enhanced ligand binding in electrospray buffers containing high concentrations of sodium, as opposed to ammonium acetate buffers. A combination of native and -omics mass spectrometry revealed the presence of a lipophilic negative allosteric modulator. These experiments highlight the advantages of implementing native mass spectrometry, from electrospray buffers containing high concentrations of physiologically relevant salts, to inform on allosteric ions or ligands with the potential to define their roles on GPCR function.

Slc1a3-2A-CreERT2 mice reveal unique features of Bergmann glia and augment a growing collection of Cre drivers and effectors in the 129S4 genetic background

Genetic variation is a primary determinant of phenotypic diversity. In laboratory mice, genetic variation can be a serious experimental confounder, and thus minimized through inbreeding. However, generalizations of results obtained with inbred strains must be made with caution, especially when working with complex phenotypes and disease models. Here we compared behavioral characteristics of C57Bl/6—the strain most widely used in biomedical research—with those of 129S4. In contrast to 129S4, C57Bl/6 demonstrated high within-strain and intra-litter behavioral hyperactivity. Although high consistency would be advantageous, the majority of disease models and transgenic tools are in C57Bl/6. We recently established six Cre driver lines and two Cre effector lines in 129S4. To augment this collection, we genetically engineered a Cre line to study astrocytes in 129S4. It was validated with two Cre effector lines: calcium indicator gCaMP5g-tdTomato and RiboTag—a tool widely used to study cell type-specific translatomes. These reporters are in different genomic loci, and in both the Cre was functional and astrocyte-specific. We found that calcium signals lasted longer and had a higher amplitude in cortical compared to hippocampal astrocytes, genes linked to a single neurodegenerative disease have highly divergent expression patterns, and that ribosome proteins are non-uniformly expressed across brain regions and cell types.

Lipid nanoparticle technologies for the study of G protein-coupled receptors in lipid environments

G protein-coupled receptors (GPCRs) are a large family of integral membrane proteins which conduct a wide range of biological roles and represent significant drug targets. Most biophysical and structural studies of GPCRs have been conducted on detergent-solubilised receptors, and it is clear that detergents can have detrimental effects on GPCR function. Simultaneously, there is increasing appreciation of roles for specific lipids in modulation of GPCR function. Lipid nanoparticles such as nanodiscs and styrene maleic acid lipid particles (SMALPs) offer opportunities to study integral membrane proteins in lipid environments, in a form that is soluble and amenable to structural and biophysical experiments. Here, we review the application of lipid nanoparticle technologies to the study of GPCRs, assessing the relative merits and limitations of each system. We highlight how these technologies can provide superior platforms to detergents for structural and biophysical studies of GPCRs and inform on roles for protein-lipid interactions in GPCR function.

Contributions and competition of Mg2+ and K+ in folding and stabilization of the Twister ribozyme

Native folded and compact intermediate states of RNA typically involve tertiary structures in the presence of divalent ions such as Mg²⁺ in a background of monovalent ions. In a recent study, we have shown how the presence of Mg²⁺ impacts the transition from partially unfolded to folded states through a "push-pull" mechanism where the ion both favors and disfavors the sampling of specific phosphate-phosphate interactions. To further understand the ion atmosphere of RNA in folded and partially folded states results from atomistic umbrella sampling and oscillating chemical potential grand canonical Monte Carlo/molecular dynamics (GCMC/MD) simulations are used to obtain atomic-level details of the distributions of Mg²⁺ and K⁺ ions around Twister RNA. Results show the presence of 100 mM Mg²⁺ to lead to increased charge neutralization over that predicted by counterion condensation theory. Upon going from partially unfolded to folded states, overall charge neutralization increases at all studied ion concentrations that, while associated with an increase in the number of direct ion-phosphate interactions, is fully accounted for by the monovalent K⁺ ions. Furthermore, K⁺ preferentially interacts with purine N7 atoms of helical regions in partially unfolded states, thereby potentially stabilizing the helical regions. Thus, both secondary helical structures and formation of tertiary structures leads to increased counterion condensation, thereby stabilizing those structural features of Twister. Notably, it is shown that K⁺ can act as a surrogate for Mg²⁺ by participating in specific interactions with nonsequential phosphate pairs that occur in the folded state, explaining the ability of Twister to self-cleave at submillimolar Mg²⁺ concentrations.

Tuning Down the Pain - An Overview of Allosteric Modulation of Opioid Receptors: Mechanisms of Modulation, Allosteric Sites, Modulator Syntheses

Opioid signaling plays a central role in pain perception. As such, it remains the main target in the development of antinociceptive agents, despite serious side effects involved. In recent years, hopes for improved opioid painkillers are rising, together with our understanding of allosterism and biased signaling mechanisms. In this review, we focus on recently discovered allosteric modulators of opioid receptors, insights into phenomena underlying their action, as well as on how they extend our understanding of mechanisms of previously known compounds. A brief overlook of their synthesis is also presented.

Impact of electronic polarizability on protein-functional group interactions

Interactions of proteins with functional groups are key to their biological functions, making it essential that they be accurately modeled. To investigate the impact of the inclusion of explicit treatment of electronic polarizability in force fields on protein-functional group interactions, the additive CHARMM and Drude polarizable force field are compared in the context of the Site-Identification by Ligand Competitive Saturation (SILCS) simulation methodology from which functional group interaction patterns with five proteins for which experimental binding affinities of multiple ligands are available, were obtained. The explicit treatment of polarizability produces significant differences in the functional group interactions in the ligand binding sites including overall enhanced binding of functional groups to the proteins. This is associated with variations of the dipole moments of solutes representative of functional groups in the binding sites relative to aqueous solution with higher dipole moments systematically occurring in the latter, though exceptions occur with positively charged methylammonium. Such variation indicates the complex, heterogeneous nature of the electronic environments of ligand binding sites and emphasizes the inherent limitation of fixed charged, additive force fields for modeling ligand-protein interactions. These effects yield more defined orientation of the functional groups in the binding pockets and a small, but systematic improvement in the ability of the SILCS method to predict the binding orientation and relative affinities of ligands to their target proteins. Overall, these results indicate that the physical model associated with the explicit treatment of polarizability along with the presence of lone pairs in a force field leads to changes in the nature of the interactions of functional groups with proteins versus that occurring with additive force fields, suggesting the utility of polarizable force fields in obtaining a more realistic understanding of protein-ligand interactions.

Mg2+ Impacts the Twister Ribozyme through Push-Pull Stabilization of Nonsequential Phosphate Pairs

RNA molecules perform a variety of biological functions for which the correct three-dimensional structure is essential, including as ribozymes where they catalyze chemical reactions. Metal ions, especially Mg2+, neutralize these negatively charged nucleic acids and specifically stabilize RNA tertiary structures as well as impact the folding landscape of RNAs as they assume their tertiary structures. Specific binding sites of Mg2+ in folded conformations of RNA have been studied extensively; however, the full range of interactions of the ion with compact intermediates and unfolded states of RNA is challenging to investigate, and the atomic details of the mechanism by which the ion facilitates tertiary structure formation is not fully known. Here, umbrella sampling combined with oscillating chemical potential Grand Canonical Monte Carlo/molecular dynamics simulations are used to capture the energetics and atomic-level details of Mg2+-RNA interactions that occur along an unfolding pathway of the Twister ribozyme. The free energy profiles reveal stabilization of partially unfolded states by Mg2+, as observed in unfolding experiments, with this stabilization being due to increased sampling of simultaneous interactions of Mg2+ with two or more nonsequential phosphate groups. Notably, these results indicate a push-pull mechanism in which the Mg2+-RNA interactions actually lead to destabilization of specific nonsequential phosphate-phosphate interactions (i.e., pushed apart), whereas other interactions are stabilized (i.e., pulled together), a balance that stabilizes unfolded states and facilitates the folding of Twister, including the formation of hydrogen bonds associated with the tertiary structure. This study establishes a better understanding of how Mg2+-ion interactions contribute to RNA structural properties and stability.

Opioid Pharmacology under the Microscope

The powerful analgesic effects of opioid drugs have captivated the interest of physicians and scientists for millennia, and the ability of opioid drugs to produce serious undesired effects has been recognized for a similar period of time (Kieffer and Evans, 2009). Many of these develop progressively with prolonged or repeated drug use and then persist, motivating particular interest in understanding how opioid drugs initiate adaptive or maladaptive modifications in neural function or regulation. Exciting advances have been made over the past several years in elucidating drug-induced changes at molecular, cellular, and physiologic scales of analysis. The present review will highlight some recent cellular studies that we believe bridge across scales and will focus on optical imaging approaches that put opioid drug action "under the microscope." SIGNIFICANCE STATEMENT: Opioid receptors are major pharmacological targets, but their signaling at the cellular level results from a complex interplay between pharmacology, regulation, subcellular localization, and membrane trafficking. This minireview discusses recent advances in understanding the cellular biology of opioid receptors, emphasizing particular topics discussed at the 50th anniversary of the International Narcotics Research Conference. Our goal is to highlight distinct signaling and regulatory properties emerging from the cellular biology of opioid receptors and discuss potential relevance to therapeutics.

Ions Everywhere? Mg2+ in the μ-Opioid GPCR and Atomic Details of Their Impact on Function

Filizola and co-workers have applied a combination of long-time molecular dynamics and oscillating chemical potential grand canonical Monte Carlo/molecular dynamics to investigate the distribution of Mg²⁺ and Na⁺ in the μ-opioid receptor and their impact on its function. Results indicate atomic details of potential mechanisms by which Mg²⁺ leads to increased efficacy of opioid analgesics. The presence of information flow between the extracellular loops and the intracellular region of the G-protein-coupled receptors that interacts with G-proteins in the presence of Mg²⁺ may be a phenomenon occurring in other G-protein-coupled receptors and, therefore, potentially of broad impact.