Structural insights into agonist-induced activation of G-protein-coupled receptors

Differential Behavior of Conformational Dynamics in Active and Inactive States of Cannabinoid Receptor 1

Cannabinoid receptor 1 (CB1) is a G protein-coupled receptor that regulates critical physiological processes including pain, appetite, and cognition. Understanding the conformational dynamics of CB1 associated with transitions between inactive and active signaling states is imperative for developing targeted modulators. Using microsecond-level all-atom molecular dynamics simulations, we identified marked differences in the conformational ensembles of inactive and active CB1 in apo. The inactive state exhibited substantially increased structural heterogeneity and plasticity compared to the more rigidified active state in the absence of stabilizing ligands. Transmembrane helices TM3 and TM7 were identified as distinguishing factors modulating the state-dependent dynamics. TM7 displayed amplified fluctuations selectively in the inactive state simulations attributed to disruption of conserved electrostatic contacts anchoring it to surrounding helices in the active state. Additionally, we identified significant reorganizations in key salt bridge and hydrogen bond networks contributing to the CB1 activation/inactivation. For instance, D213-Y224 hydrogen bond and D184-K192 salt bridge showed marked rearrangements between the states. Collectively, these findings reveal the specialized role of TM7 in directing state-dependent CB1 dynamics through electrostatic switch mechanisms. By elucidating the intrinsic enhanced flexibility of inactive CB1, this study provides valuable insights into the conformational landscape enabling functional transitions. Our perspective advances understanding of CB1 activation mechanisms and offers opportunities for structure-based drug discovery targeting the state-specific conformational dynamics of this receptor.

The complex nature of CXCR4 mutations in WHIM syndrome

Heterozygous autosomal dominant mutations in the CXCR4 gene cause WHIM syndrome, a severe combined immunodeficiency disorder. The mutations primarily affect the C-terminal region of the CXCR4 chemokine receptor, specifically several potential phosphorylation sites critical for agonist (CXCL12)-mediated receptor internalization and desensitization. Mutant receptors have a prolonged residence time on the cell surface, leading to hyperactive signaling that is responsible for some of the symptoms of WHIM syndrome. Recent studies have shown that the situation is more complex than originally thought, as mutant WHIM receptors and CXCR4 exhibit different dynamics at the cell membrane, which also influences their respective cellular functions. This review examines the functional mechanisms of CXCR4 and the impact of WHIM mutations in both physiological and pathological conditions.

GPER binding site detection and description: A flavonoid-based docking and molecular dynamics simulations study

Flavonoids, a phenolic compounds class widely distributed in the plant kingdom, have attracted much interest for their implications on several health and disease processes. Usually, the consumption of this type of compounds is approximately 1 g/d, primarily obtained from cereals, chocolate, and dry legumes ensuring its beneficial role in maintaining the homeostasis of the human body. In this context, in cancer disease prominent data points to the role of flavonoids as adjuvant treatment aimed at the regression of the disease. GPER, an estrogen receptor on the cell surface, has been postulated as a probable orchestrator of the beneficial effects of several flavonoids through modulation/inhibition of various mechanisms that lead to cancer progression. Therefore, applying pocket and cavity protein detection and docking and molecular dynamics simulations (MD), we generate, from a cluster composed of 39 flavonoids, crucial insights into the potential role as GPER ligands, of Puerarin, Isoquercetin, Kaempferol 3-O-glucoside and Petunidin 3-O-glucoside, aglycones whose sugar moiety delimits a new described sugar-acceptor sub-cavity into the cavity binding site on the receptor, as well as of the probable activation mechanism of the receptor and the pivotal residues involved in it. Altogether, our results shed light on the potential use of the aforementioned flavonoids as GPER ligands and for further evaluations in in vitro and in vivo assays to elucidate their probable anti-cancer activity.

Differential Behavior of Conformational Dynamics in Active and Inactive States of Cannabinoid Receptor 1

The cannabinoid receptor CB1 is a G protein-coupled receptor that regulates critical physiological processes including pain, appetite, and cognition. Understanding the conformational dynamics of CB1 associated with transitions between inactive and active signaling states is imperative for developing targeted modulators. Using microsecond-level all-atom molecular dynamics simulations, we identified marked differences in the conformational ensembles of inactive and active CB1 states in apo conditions. The inactive state exhibited substantially increased structural heterogeneity and plasticity compared to the more rigidified active state in the absence of stabilizing ligands. Transmembrane helices TM3 and TM7 were identified as distinguishing factors modulating the state-dependent dynamics. TM7 displayed amplified fluctuations selectively in the inactive state simulations attributed to disruption of conserved electrostatic contacts anchoring it to surrounding helices in the active state. Additionally, we identified significant reorganization of key salt bridge and hydrogen bond networks known to control CB1 activation between states. For instance, a conserved D213-Y224 hydrogen bond and D184-K192 salt bridge interactions showed marked rearrangements between the states. Collectively, these findings reveal the specialized role of TM7 in directing state-dependent CB1 dynamics through electrostatic switch mechanisms. By elucidating the intrinsic enhanced flexibility of inactive CB1, this study provides valuable insights into the conformational landscape enabling functional transitions. Our perspective advances understanding of CB1 activation mechanisms and offers opportunities for structure-based drug discovery targeting the state-specific conformational dynamics of this receptor.

Glu1022.53-Mediated Early Conformational Changes in the Process of Light-Induced Green Cone Pigment Activation

Ligand-triggered activation of G protein-coupled receptors (GPCRs) relies on the phenomenon of loose allosteric coupling, which involves conformational alterations spanning from the extracellular ligand-binding domain to the cytoplasmic region, where interactions with G proteins occur. During the GPCR activation process, several intermediate and equilibrium states orchestrate the movement of the flexible and rigid transmembrane (TM) segments of the GPCR. Monitoring early conformational changes is important in unraveling the structural intricacies of the loose allosteric coupling. Here, we focus on the lumi intermediate formed by thermal relaxation from the initial photointermediate, batho in primate green cone pigment (MG), a light-sensitive GPCR responsible for color vision. Our findings from light-induced Fourier transform infrared difference spectroscopy reveal its similarity with rhodopsin, which mediates twilight vision, specifically involving the flip motion of the β-ionone ring, the relaxation of the torsional structure of the retinal, and local perturbations in the α-helix upon lumi intermediate formation. Conversely, we observe a hydrogen bond modification specific to MG's protonated carboxylic acid, identifying its origin as Glu1022.53 situated in TM2. The weakening of the hydrogen bond strength at Glu1022.53 during the transition from the batho to the lumi intermediates corresponds to a slight outward movement of TM2. Additionally, within the X-ray crystal structure of the rhodopsin lumi intermediate, we note the relocation of the Met862.53 side chain in TM2, expanding the volume of the retinal binding pocket. Consequently, the position of 2.53 emerges as the early step in the conformational shift toward light-induced activation. Moreover, given the prevalence of IR-insensitive hydrophobic amino acids at position 2.53 in many rhodopsin-like GPCRs, including rhodopsin, the hydrogen bond alteration in the C═O stretching band at Glu1022.53 of MG can be used as a probe for tracing conformational changes during the GPCR activation process.

Illuminating the brain-genetically encoded single wavelength fluorescent biosensors to unravel neurotransmitter dynamics

Understanding how neuronal networks generate complex behavior is one of the major goals of Neuroscience. Neurotransmitter and Neuromodulators are crucial for information flow between neurons and understanding their dynamics is the key to unravel their role in behavior. To understand how the brain transmits information and how brain states arise, it is essential to visualize the dynamics of neurotransmitters, neuromodulators and neurochemicals. In the last five years, an increasing number of single-wavelength biosensors either based on periplasmic binding proteins (PBPs) or on G-protein-coupled receptors (GPCR) have been published that are able to detect neurotransmitter release in vitro and in vivo with high spatial and temporal resolution. Here we review and discuss recent progress in the development of these sensors, their limitations and future directions.

Structure, function and drug discovery of GPCR signaling

G protein-coupled receptors (GPCRs) are versatile and vital proteins involved in a wide array of physiological processes and responses, such as sensory perception (e.g., vision, taste, and smell), immune response, hormone regulation, and neurotransmission. Their diverse and essential roles in the body make them a significant focus for pharmaceutical research and drug development. Currently, approximately 35% of marketed drugs directly target GPCRs, underscoring their prominence as therapeutic targets. Recent advances in structural biology have substantially deepened our understanding of GPCR activation mechanisms and interactions with G-protein and arrestin signaling pathways. This review offers an in-depth exploration of both traditional and recent methods in GPCR structure analysis. It presents structure-based insights into ligand recognition and receptor activation mechanisms and delves deeper into the mechanisms of canonical and noncanonical signaling pathways downstream of GPCRs. Furthermore, it highlights recent advancements in GPCR-related drug discovery and development. Particular emphasis is placed on GPCR selective drugs, allosteric and biased signaling, polyphamarcology, and antibody drugs. Our goal is to provide researchers with a thorough and updated understanding of GPCR structure determination, signaling pathway investigation, and drug development. This foundation aims to propel forward-thinking therapeutic approaches that target GPCRs, drawing upon the latest insights into GPCR ligand selectivity, activation, and biased signaling mechanisms.

Structural basis of G protein-Coupled receptor CMKLR1 activation and signaling induced by a chemerin-derived agonist

Chemokine-like receptor 1 (CMKLR1), also known as chemerin receptor 23 (ChemR23) or chemerin receptor 1, is a chemoattractant G protein-coupled receptor (GPCR) that responds to the adipokine chemerin and is highly expressed in innate immune cells, including macrophages and neutrophils. The signaling pathways of CMKLR1 can lead to both pro- and anti-inflammatory effects depending on the ligands and physiological contexts. To understand the molecular mechanisms of CMKLR1 signaling, we determined a high-resolution cryo-electron microscopy (cryo-EM) structure of the CMKLR1-Gi signaling complex with chemerin9, a nanopeptide agonist derived from chemerin, which induced complex phenotypic changes of macrophages in our assays. The cryo-EM structure, together with molecular dynamics simulations and mutagenesis studies, revealed the molecular basis of CMKLR1 signaling by elucidating the interactions at the ligand-binding pocket and the agonist-induced conformational changes. Our results are expected to facilitate the development of small molecule CMKLR1 agonists that mimic the action of chemerin9 to promote the resolution of inflammation.

Pro-phagocytic function and structural basis of GPR84 signaling

GPR84 is a unique orphan G protein-coupled receptor (GPCR) that can be activated by endogenous medium-chain fatty acids (MCFAs). The signaling of GPR84 is largely pro-inflammatory, which can augment inflammatory response, and GPR84 also functions as a pro-phagocytic receptor to enhance phagocytic activities of macrophages. In this study, we show that the activation of GPR84 by the synthetic agonist 6-OAU can synergize with the blockade of CD47 on cancer cells to induce phagocytosis of cancer cells by macrophages. We also determine a high-resolution structure of the GPR84-Gi signaling complex with 6-OAU. This structure reveals an occluded binding pocket for 6-OAU, the molecular basis of receptor activation involving non-conserved structural motifs of GPR84, and an unusual Gi-coupling interface. Together with computational docking and simulations studies, this structure also suggests a mechanism for the high selectivity of GPR84 for MCFAs and a potential routes of ligand binding and dissociation. These results provide a framework for understanding GPR84 signaling and developing new drugs targeting GPR84.

Structural basis for the ligand recognition and signaling of free fatty acid receptors

Free fatty acid receptors 1-4 (FFA1-4) are class A G protein-coupled receptors (GPCRs). FFA1-3 share substantial sequence similarity whereas FFA4 is unrelated. Despite this FFA1 and FFA4 are activated by the same range of long chain fatty acids (LCFAs) whilst FFA2 and FFA3 are instead activated by short chain fatty acids (SCFAs) generated by the intestinal microbiota. Each of FFA1, 2 and 4 are promising targets for novel drug development in metabolic and inflammatory conditions. To gain insights into the basis of ligand interactions with, and molecular mechanisms underlying activation of, FFAs by LCFAs and SCFAs, we determined the active structures of FFA1 and FFA4 bound to the polyunsaturated LCFA docosahexaenoic acid (DHA), FFA4 bound to the synthetic agonist TUG-891, as well as SCFA butyrate-bound FFA2, each complexed with an engineered heterotrimeric Gq protein (miniGq), by cryo-electron microscopy. Together with computational simulations and mutagenesis studies, we elucidated the similarities and differences in the binding modes of fatty acid ligands with varying chain lengths to their respective GPCRs. Our findings unveil distinct mechanisms of receptor activation and G protein coupling. We anticipate that these outcomes will facilitate structure-based drug development and underpin future research to understand allosteric modulation and biased signaling of this group of GPCRs.

Structural basis of CMKLR1 signaling induced by chemerin9

Chemokine-like receptor 1 (CMKLR1), also known as chemerin receptor 23 (ChemR23) or chemerin receptor 1, is a chemoattractant G protein-coupled receptor (GPCR) that responds to the adipokine chemerin and is highly expressed in innate immune cells, including macrophages and neutrophils. The signaling pathways of CMKLR1 can lead to both pro- and anti-inflammatory effects depending on the ligands and physiological contexts. To understand the molecular mechanisms of CMKLR1 signaling, we determined a high-resolution cryo-electron microscopy (cryo-EM) structure of the CMKLR1-Gi signaling complex with chemerin9, a nanopeptide agonist derived from chemerin, which induced complex phenotypic changes of macrophages in our assays. The cryo-EM structure, together with molecular dynamics simulations and mutagenesis studies, revealed the molecular basis of CMKLR1 signaling by elucidating the interactions at the ligand-binding pocket and the agonist-induced conformational changes. Our results are expected to facilitate the development of small molecule CMKLR1 agonists that mimic the action of chemerin9 to promote the resolution of inflammation.

Inhibiting a promiscuous GPCR: iterative discovery of bitter taste receptor ligands

The human GPCR family comprises circa 800 members, activated by hundreds of thousands of compounds. Bitter taste receptors, TAS2Rs, constitute a large and distinct subfamily, expressed orally and extra-orally and involved in physiological and pathological conditions. TAS2R14 is the most promiscuous member, with over 150 agonists and 3 antagonists known prior to this study. Due to the scarcity of inhibitors and to the importance of chemical probes for exploring TAS2R14 functions, we aimed to discover new ligands for this receptor, with emphasis on antagonists. To cope with the lack of experimental structure of the receptor, we used a mixed experimental/computational methodology which iteratively improved the performance of the predicted structure. The increasing number of active compounds, obtained here through experimental screening of FDA-approved drug library, and through chemically synthesized flufenamic acid derivatives, enabled the refinement of the binding pocket, which in turn improved the structure-based virtual screening reliability. This mixed approach led to the identification of 10 new antagonists and 200 new agonists of TAS2R14, illustrating the untapped potential of rigorous medicinal chemistry for TAS2Rs. 9% of the ~ 1800 pharmaceutical drugs here tested activate TAS2R14, nine of them at sub-micromolar concentrations. The iterative framework suggested residues involved in the activation process, is suitable for expanding bitter and bitter-masking chemical space, and is applicable to other promiscuous GPCRs lacking experimental structures.

Role of Monomer/Tetramer Equilibrium of Rod Visual Arrestin in the Interaction with Phosphorylated Rhodopsin

The phototransduction cascade in vertebrate rod visual cells is initiated by the photoactivation of rhodopsin, which enables the activation of the visual G protein transducin. It is terminated by the phosphorylation of rhodopsin, followed by the binding of arrestin. Here we measured the solution X-ray scattering of nanodiscs containing rhodopsin in the presence of rod arrestin to directly observe the formation of the rhodopsin/arrestin complex. Although arrestin self-associates to form a tetramer at physiological concentrations, it was found that arrestin binds to phosphorylated and photoactivated rhodopsin at 1:1 stoichiometry. In contrast, no complex formation was observed for unphosphorylated rhodopsin upon photoactivation, even at physiological arrestin concentrations, suggesting that the constitutive activity of rod arrestin is sufficiently low. UV-visible spectroscopy demonstrated that the rate of the formation of the rhodopsin/arrestin complex well correlates with the concentration of arrestin monomer rather than the tetramer. These findings indicate that arrestin monomer, whose concentration is almost constant due to the equilibrium with the tetramer, binds to phosphorylated rhodopsin. The arrestin tetramer would act as a reservoir of monomer to compensate for the large changes in arrestin concentration in rod cells caused by intense light or adaptation.

Structural basis of lysophosphatidylserine receptor GPR174 ligand recognition and activation

Lysophosphatidylserine (LysoPS) is a lipid mediator that induces multiple cellular responses through binding to GPR174. Here, we present the cryo-electron microscopy (cryo-EM) structure of LysoPS-bound human GPR174 in complex with G_s protein. The structure reveals a ligand recognition mode, including the negatively charged head group of LysoPS forms extensive polar interactions with surrounding key residues of the ligand binding pocket, and the L-serine moiety buries deeply into a positive charged cavity in the pocket. In addition, the structure unveils a partially open pocket on transmembrane domain helix (TM) 4 and 5 for a lateral entry of ligand. Finally, the structure reveals a G_s engaging mode featured by a deep insertion of a helix 5 (αH5) and extensive polar interactions between receptor and αH5. Taken together, the information revealed by our structural study provides a framework for understanding LysoPS signaling and a rational basis for designing LysoPS receptor-targeting drugs.

Superconserved receptors expressed in the brain: Expression, function, motifs and evolution of an orphan receptor family

GPR27, GPR85 and GPR173 constitute a small family of G protein-coupled receptors (GPCR) that share the distinctive characteristics of being highly conserved throughout vertebrate evolution and predominantly expressed in the brain. Accordingly, they have been coined as "Superconserved Receptors Expressed in the Brain" (SREB), although their expression profile is more complex than what was originally thought. SREBs have no known validated endogenous ligands and are thus labeled as "orphan" receptors. The investigation of this particular category of uncharacterized receptors holds great promise both in terms of physiology and drug development. In the largest GPCR family, the Rhodopsin-like or Class A, around 100 receptors are considered orphans. Because GPCRs are the most successful source of drug targets, the discovery of a novel function or ligand most likely will lead to significant breakthroughs for the discovery of innovative therapies. The high level of conservation is one of the characteristic features of the SREBs. We propose herein a detailed analysis of the putative evolutionary origin of this family. We highlight the properties that distinguish SREBs from other rhodopsin-like GPCRs. We present the current evidence for these receptors downstream signaling pathways and functions. We discuss the pharmacological challenge for the identification of natural or synthetic ligands of orphan receptors like SREBs. The different SREB-related scientific questions are presented with a highlight on what should be addressed in the near future, including the confirmation of published evidence and their validation as drug targets. In particular, we discuss in which pathological conditions these receptors may be of great relevance to solve unmet medical needs.

Biased Activation Mechanism Induced by GPCR Heterodimerization: Observations from μOR/δOR Dimers

GPCRs regulate multiple intracellular signaling cascades. Biasedly activating one signaling pathway over the others provides additional clinical utility to optimize GPCR-based therapies. GPCR heterodimers possess different functions from their monomeric states, including their selectivity to different transducers. However, the biased signaling mechanism induced by the heterodimerization remains unclear. Motivated by the issue, we select an important GPCR heterodimer (μOR/δOR heterodimer) as a case and use microsecond Gaussian accelerated molecular dynamics simulation coupled with potential of mean force and protein structure network (PSN) to probe mechanisms regarding the heterodimerization-induced constitutive β-arrestin activity and efficacy change of the agonist DAMGO. The results show that only the lowest energy state of the μOR/δOR heterodimer, which adopts a slightly outward shift of TM6 and an ICL2 conformation close to the receptor core, can selectively accommodate β-arrestins. PSN further reveals important roles of H8, ICL1, and ICL2 in regulating the constitutive β-arrestin-biased activity for the apo μOR/δOR heterodimer. In addition, the heterodimerization can allosterically alter the binding mode of DAMGO mainly by means of W7.35. Consequently, DAMGO transmits the structural signal mainly through TM6 and TM7 in the dimer, rather than TM3 similar to the μOR monomer, thus changing the efficacy of DAMGO from a balanced agonist to the β-arrestin-biased one. On the other side, the binding of DAMGO to the heterodimer can stabilize μOR/δOR heterodimers through a stronger interaction of TM1/TM1 and H8/H8, accordingly enhancing the interaction of μOR with δOR and the binding affinity of the dimer to the β-arrestin. The agonist DAMGO does not change main compositions of the regulation network from the dimer interface to the transducer binding pocket of the μOR protomer, but induces an increase in the structural communication of the network, which should contribute to the enhanced β-arrestin coupling. Our observations, for the first time, reveal the molecular mechanism of the biased signaling induced by the heterodimerization for GPCRs, which should be beneficial to more comprehensively understand the GPCR bias signaling.

Formylpeptide receptor 2: Nomenclature, structure, signalling and translational perspectives: IUPHAR review 35

We discuss the fascinating pharmacology of formylpeptide receptor 2 (FPR2; often referred to as FPR2/ALX since it binds lipoxin A₄ ). Initially identified as a low-affinity 'relative' of FPR1, FPR2 presents complex and diverse biology. For instance, it is activated by several classes of agonists (from peptides to proteins and lipid mediators) and displays diverse expression patterns on myeloid cells as well as epithelial cells and endothelial cells, to name a few. Over the last decade, the pharmacology of FPR2 has progressed from being considered a weak chemotactic receptor to a master-regulator of the resolution of inflammation, the second phase of the acute inflammatory response. We propose that exploitation of the biology of FPR2 offers innovative ways to rectify chronic inflammatory states and represents a viable avenue to develop novel therapies. Recent elucidation of FPR2 structure will facilitate development of the anti-inflammatory and pro-resolving drugs of next decade.

Computational Methods for Understanding the Selectivity and Signal Transduction Mechanism of Aminomethyl Tetrahydronaphthalene to Opioid Receptors

Opioid receptors are members of the group of G protein-couple receptors, which have been proven to be effective targets for treating severe pain. The interactions between the opioid receptors and corresponding ligands and the receptor’s activation by different agonists have been among the most important fields in opioid research. In this study, with compound M1, an active metabolite of tramadol, as the clue compound, several aminomethyl tetrahydronaphthalenes were designed, synthesized and assayed upon opioid receptors. With the resultant compounds FW-AII-OH-1 (Ki = 141.2 nM for the κ opioid receptor), FW-AII-OH-2 (Ki = 4.64 nM for the δ opioid receptor), FW-DI-OH-2 (Ki = 8.65 nM for the δ opioid receptor) and FW-DIII-OH-2 (Ki = 228.45 nM for the δ opioid receptor) as probe molecules, the structural determinants responsible for the subtype selectivity and activation mechanisms were further investigated by molecular modeling and molecular dynamics simulations. It was shown that Y^7.43 was a key residue in determining the selectivity of the three opioid receptors, and W^6.58 was essential for the selectivity of the δ opioid receptor. A detailed stepwise discovered agonist-induced signal transduction mechanism of three opioid receptors by aminomethyl tetrahydronaphthalene compounds was proposed: the 3–7 lock between TM3 and TM7, the DRG lock between TM3 and TM6 and rearrangement of I^3.40, P^5.50 and F^6.44, which resulted in the cooperative movement in 7 TMs. Then, the structural relaxation left room for the binding of the G protein at the intracellular site, and finally the opioid receptors were activated.

Molecular recognition of formylpeptides and diverse agonists by the formylpeptide receptors FPR1 and FPR2

The formylpeptide receptors (FPRs) mediate pattern recognition of formylated peptides derived from invading pathogens or mitochondria from dead host cells. They can also sense other structurally distinct native peptides and even lipid mediators to either promote or resolve inflammation. Pharmacological targeting of FPRs represents a novel therapeutic approach in treating inflammatory diseases. However, the molecular mechanisms underlying FPR ligand recognition are elusive. We report cryo-EM structures of G_i-coupled FPR1 and FPR2 bound to a formylpeptide and G_i-coupled FPR2 bound to two synthetic peptide and small-molecule agonists. Together with mutagenesis data, our structures reveal the molecular mechanism of formylpeptide recognition by FPRs and structural variations of FPR1 and FPR2 leading to their different ligand preferences. Structural analysis also suggests that diverse FPR agonists sample a conserved activation chamber at the bottom of ligand-binding pockets to activate FPRs. Our results provide a basis for rational drug design on FPRs.

Dimeric Rhodopsin R135L Mutant-Transducin-like Complex Sheds Light on Retinitis Pigmentosa Misfunctions

Rhodopsin (RHO) is a light-sensitive pigment in the retina and the main prototypical protein of the G-protein-coupled receptor (GCPR) family. After receiving a light stimulus, RHO and its cofactor retinylidene undergo a series of structural changes that initiate an intricate transduction mechanism. Along with RHO, other partner proteins play key roles in the signaling pathway. These include transducin, a GTPase, kinases that phosphorylate RHO, and arrestin (Arr), which ultimately stops the signaling process and promotes RHO regeneration. A large number of RHO genetic mutations may lead to very severe retinal dysfunction and eventually to impaired dark adaptation disease called autosomal dominant retinitis pigmentosa (adRP). In this study, we used molecular dynamics (MD) simulations to evaluate the different behaviors of the dimeric form of wild-type RHO (WT dRHO) and its mutant at position 135 of arginine to leucine (dR135L), both in the free (noncomplexed) and in complex with the transducin-like protein (Gtl). Gtl is a heterotrimeric model composed of a mixture of human and bovine G proteins. Our calculations allow us to explain how the mutation causes structural changes in the RHO dimer and how this can affect the signal that transducin generates when it is bound to RHO. Moreover, the structural modifications induced by the R135L mutation can also account for other misfunctions observed in the up- and downstream signaling pathways. The mechanism of these dysfunctions, together with the transducin activity reduction, provides structure-based explanations of the impairment of some key processes that lead to adRP.

Deconstructing the transmembrane core of class A G protein-coupled receptors

Class A G protein-coupled receptors have evolved to recognize ligands ranging from small-molecule odorants to proteins. Although they are among the most diverse membrane receptors in eukaryotic organisms, they possess a highly conserved core within their seven-transmembrane helix framework. The conservation of the transmembrane core has led to the idea of a common mechanism by which ligand binding is coupled to the outward rotation of helix H6, the hallmark of an active receptor. Nevertheless, there is still no consensus on the mechanism of coupling or on the roles of specific residues within the core. Recent insights from crystallography and NMR spectroscopy provide a way to decompose the core into its essential structural and functional elements that shed new light on this important region.

Functional molecular switches of mammalian G protein-coupled bitter-taste receptors

Bitter taste receptors (TAS2Rs) are a poorly understood subgroup of G protein-coupled receptors (GPCRs). The experimental structure of these receptors has yet to be determined, and key-residues controlling their function remain mostly unknown. We designed an integrative approach to improve comparative modeling of TAS2Rs. Using current knowledge on class A GPCRs and existing experimental data in the literature as constraints, we pinpointed conserved motifs to entirely re-align the amino-acid sequences of TAS2Rs. We constructed accurate homology models of human TAS2Rs. As a test case, we examined the accuracy of the TAS2R16 model with site-directed mutagenesis and in vitro functional assays. This combination of in silico and in vitro results clarifies sequence-function relationships and proposes functional molecular switches that encode agonist sensing and downstream signaling mechanisms within mammalian TAS2Rs sequences.

Understanding the role of the CB1 toggle switch in interaction networks using molecular dynamics simulation

The cannabinoid receptor 1 (CB1) is a class A G-protein coupled receptor (GPCR) that can exert various effects on the human body through the endocannabinoid system. Understanding CB1 activation has many benefits for the medical use of cannabinoids. A previous study reported that CB1 has two notable residues referred to as the toggle switch, F3.36 and W6.48, which are important for its activation mechanism. We performed a molecular dynamics simulation with a mutation in the toggle switch to examine its role in active and inactive states. We also examined structural changes, the residue-residue interaction network, and the interaction network among helices and loops of wildtype and mutant CB1 for both activation states. As a result, we found that the energetic changes in the hydrogen-bond network of the Na⁺ pocket, extracellular N-terminus-TM2-ECL1-TM3 interface including D2.63-K3.28 salt-bridge, and extracellular ECL2-TM5-ECL3-TM6 interface directly linked to the toggle switch contribute to the stability of CB1 by the broken aromatic interaction of the toggle switch. It makes the conformation of inactive CB1 receptor to be unstable. Our study explained the role of the toggle switch regarding the energetic interactions related to the Na⁺ pocket and extracellular loop interfaces, which could contribute to a better understanding of the activation mechanism of CB1.

Probing Allosteric Regulation Mechanism of W7.35 on Agonist-Induced Activity for μOR by Mutation Simulation

The residue located at 15 positions before the most conserved residue in TM7 (7.35 of Ballesteros-Weinstein number) plays important roles in ligand binding and the receptor activity for class A GPCRs. Nevertheless, its regulation mechanism has not been clearly clarified in experiments, and some controversies also exist for its impact on μ-opioid receptors (μOR) bound by agonists. Thus, we chose the μ-opioid receptor (μOR) of class A GPCRs as a representative and utilized a microsecond accelerated molecular dynamics simulation (aMD) coupled with a protein structure network (PSN) to explore the effect of W318^7.35 on its functional activity induced by the agonist endomorphin2 mainly by a comparison of the wild system and its W7.35A mutant. When endomorphin2 binds to the wild-type μOR, TM6 in μOR moves outward to form an open intracellular conformation that is beneficial to accommodating the β-arrestin transducer, rather than the G-protein transducer due to the clash with the α5 helix of G-protein, thus acting as a β-arrestin biased agonist. However, the W318A mutation induces the intracellular part of μOR to form a closed state, which disfavors coupling with either G-protein or β-arrestin. The allosteric pathway analysis further reveals that the binding of endomorphin2 to the wild-type μOR transmits more activation signals to the β-arrestin binding site while the W318A mutation induces more structural signals to transmit to common binding residues of the G protein and β-arrestin. More interestingly, the residue at the 7.35 position regulates the shortest allosteric pathway in indirect ways by influencing the interactions between other ligand-binding residues and endomorphin2. W293^6.48 and F289^6.44 are important for regulating the different activities of μOR induced either by the agonist or by the mutation. Y336^7.53, F343^8.50, and D340^8.47 play crucial roles in activating the β-arrestin biased signal induced by the agonist endomorphin2, while L158^3.43 and V286^6.41 devote important contributions to the change in the activity of endomorphin2 from the β-arrestin biased agonist to the antagonist upon the W318A mutation.

Molecular basis for lipid recognition by the prostaglandin D2 receptor CRTH2

Prostaglandin D₂ (PGD₂) signals through the G protein-coupled receptor (GPCR) CRTH2 to mediate various inflammatory responses. CRTH2 is the only member of the prostanoid receptor family that is phylogenetically distant from others, implying a nonconserved mechanism of lipid action on CRTH2. Here, we report a crystal structure of human CRTH2 bound to a PGD₂ derivative, 15R-methyl-PGD₂ (15mPGD₂), by serial femtosecond crystallography. The structure revealed a "polar group in"-binding mode of 15mPGD₂ contrasting the "polar group out"-binding mode of PGE₂ in its receptor EP3. Structural comparison analysis suggested that these two lipid-binding modes, associated with distinct charge distributions of ligand-binding pockets, may apply to other lipid GPCRs. Molecular dynamics simulations together with mutagenesis studies also identified charged residues at the ligand entry port that function to capture lipid ligands of CRTH2 from the lipid bilayer. Together, our studies suggest critical roles of charge environment in lipid recognition by GPCRs.

Valine-279 Deletion-Mutation on Arginine Vasopressin Receptor 2 Causes Obstruction in G-Protein Binding Site: A Clinical Nephrogenic Diabetes Insipidus Case and Its Sub-Molecular Pathogenic Analysis

Congenital nephrogenic diabetes insipidus (CNDI) is a genetic disorder caused by mutations in arginine vasopressin receptor 2 (AVPR2) or aquaporin 2 genes, rendering collecting duct cells insensitive to the peptide hormone arginine vasopressin stimulation for water reabsorption. This study reports a first identified AVPR2 mutation in Taiwan and demonstrates our effort to understand the pathogenesis caused by applying computational structural analysis tools. The CNDI condition of an 8-month-old male patient was confirmed according to symptoms, family history, and DNA sequence analysis. The patient was identified to have a valine 279 deletion-mutation in the AVPR2 gene. Cellular experiments using mutant protein transfected cells revealed that mutated AVPR2 is expressed successfully in cells and localized on cell surfaces. We further analyzed the pathogenesis of the mutation at sub-molecular levels via long-term molecular dynamics (MD) simulations and structural analysis. The MD simulations showed while the structure of the extracellular ligand-binding domain remains unchanged, the mutation alters the direction of dynamic motion of AVPR2 transmembrane helix 6 toward the center of the G-protein binding site, obstructing the binding of G-protein, thus likely disabling downstream signaling. This study demonstrated that the computational approaches can be powerful tools for obtaining valuable information on the pathogenesis induced by mutations in G-protein-coupled receptors. These methods can also be helpful in providing clues on potential therapeutic strategies for CNDI.

Capturing Peptide-GPCR Interactions and Their Dynamics

Many biological functions of peptides are mediated through G protein-coupled receptors (GPCRs). Upon ligand binding, GPCRs undergo conformational changes that facilitate the binding and activation of multiple effectors. GPCRs regulate nearly all physiological processes and are a favorite pharmacological target. In particular, drugs are sought after that elicit the recruitment of selected effectors only (biased ligands). Understanding how ligands bind to GPCRs and which conformational changes they induce is a fundamental step toward the development of more efficient and specific drugs. Moreover, it is emerging that the dynamic of the ligand–receptor interaction contributes to the specificity of both ligand recognition and effector recruitment, an aspect that is missing in structural snapshots from crystallography. We describe here biochemical and biophysical techniques to address ligand–receptor interactions in their structural and dynamic aspects, which include mutagenesis, crosslinking, spectroscopic techniques, and mass-spectrometry profiling. With a main focus on peptide receptors, we present methods to unveil the ligand–receptor contact interface and methods that address conformational changes both in the ligand and the GPCR. The presented studies highlight a wide structural heterogeneity among peptide receptors, reveal distinct structural changes occurring during ligand binding and a surprisingly high dynamics of the ligand–GPCR complexes.

Internal water channel formation in CXCR4 is crucial for Gi-protein coupling upon activation by CXCL12

Chemokine receptor CXCR4 is a major drug target for numerous diseases because of its involvement in the regulation of cell migration and the developmental process. In this study, atomic-level molecular dynamics simulations were used to determine the activation mechanism and internal water formation of CXCR4 in complex with chemokine CXCL12 and G_i-protein. The results indicated that CXCL12-bound CXCR4 underwent transmembrane 6 (TM6) outward movement and a decrease in tyrosine toggle switch by eliciting the breakage of hydrophobic layer to form a continuous internal water channel. In the GDP-bound G_αi-protein state, the rotation and translation of the α5-helix of G_αi-protein toward the cytoplasmic pocket of CXCR4 induced an increase in interdomain distance for GDP leaving. Finally, an internal water channel formation model was proposed based on our simulations for CXCL12-bound CXCR4 in complex with G_αi-protein upon activation for downstream signaling. This model could be useful in anticancer drug development.

Camelid-derived single-chain antibodies in hemostasis: Mechanistic, diagnostic, and therapeutic applications

Hemostasis is a complex process involving the concerted action of molecular and vascular components. Its basic understanding as well as diagnostic and therapeutic aspects have greatly benefited from the use of monoclonal antibodies. Interestingly, camelid-derived single-domain antibodies (sdAbs), also known as VHH or nanobodies, have become available during the previous 2 decades as alternative tools in this regard. Compared to classic antibodies, sdAbs are easier to produce and their small size facilitates their engineering and functionalization. It is not surprising, therefore, that sdAbs are increasingly used in hemostasis-related research. In addition, they have the capacity to recognize unique epitopes unavailable to full monoclonal antibodies. This property can be used to develop novel diagnostic tests identifying conformational variants of hemostatic proteins. Examples include sdAbs that bind active but not globular von Willebrand factor or free factor VIIa but not tissue factor-bound factor VIIa. Finally, sdAbs have a high therapeutic potential, exemplified by caplacizumab, a homodimeric sdAb targeting von Willebrand factor that is approved for the treatment of thrombotic thrombocytopenic purpura. In this review, the various applications of sdAbs in thrombosis and hemostasis-related research, diagnostics, and therapeutic strategies will be discussed.

Using modern approaches to sedimentation velocity to detect conformational changes in proteins

It has been known for decades that proteins undergo conformational changes in response to binding ligands. Such changes are usually accompanied by a loss of entropy by the protein, and thus conformational changes are integral to the thermodynamics of ligand association. Methods to detect these alterations are numerous; here, we focus on the sedimentation velocity (SV) mode of AUC, which has several advantages, including ease of use and rigorous data-selection criteria. In SV, it is assumed that conformational changes manifest primarily as differences in the sedimentation coefficient (the s-value). Two methods of determining s-value differences were assessed. The first method used the widely adopted c(s) distribution to gather statistics on the s-value differences to determine whether the observed changes were reliable. In the second method, a decades-old technique called "difference SV" was revived and updated to address its viability in this era of modern instrumentation. Both methods worked well to determine the extent of conformational changes to three model systems. Both simulations and experiments were used to explore the strengths and limitations of the methods. Finally, software incorporating these methodologies was produced.

Copper-mediated thiol potentiation and mutagenesis-guided modeling suggest a highly conserved copper-binding motif in human OR2M3

Sulfur-containing compounds within a physiological relevant, natural odor space, such as the key food odorants, typically constitute the group of volatiles with the lowest odor thresholds. The observation that certain metals, such as copper, potentiate the smell of sulfur-containing, metal-coordinating odorants led to the hypothesis that their cognate receptors are metalloproteins. However, experimental evidence is sparse-so far, only one human odorant receptor, OR2T11, and a few mouse receptors, have been reported to be activated by sulfur-containing odorants in a copper-dependent way, while the activation of other receptors by sulfur-containing odorants did not depend on the presence of metals. Here we identified an evolutionary conserved putative copper interaction motif CC/CSSH, comprising two copper-binding sites in TMH5 and TMH6, together with the binding pocket for 3-mercapto-2-methylpentan-1-ol in the narrowly tuned human receptor OR2M3. To characterize the copper-binding motif, we combined homology modeling, docking studies, site-directed mutagenesis, and functional expression of recombinant ORs in a cell-based, real-time luminescence assay. Ligand activation of OR2M3 was potentiated in the presence of copper. This effect of copper was mimicked by ionic and colloidal silver. In two broadly tuned receptors, OR1A1 and OR2W1, which did not reveal a putative copper interaction motif, activation by their most potent, sulfur-containing key food odorants did not depend on the presence of copper. Our results suggest a highly conserved putative copper-binding motif to be necessary for a copper-modulated and thiol-specific function of members from three subfamilies of family 2 ORs.

Probing the correlation between ligand efficacy and conformational diversity at the α1A-adrenoreceptor reveals allosteric coupling of its microswitches

G protein-coupled receptors (GPCRs) use a series of conserved microswitches to transmit signals across the cell membrane via an allosteric network encompassing the ligand-binding site and the G protein-binding site. Crystal structures of GPCRs provide snapshots of their inactive and active states, but poorly describe the conformational dynamics of the allosteric network that underlies GPCR activation. Here, we analyzed the correlation between ligand binding and receptor conformation of the α1A-adrenoreceptor, a GPCR that stimulates smooth muscle contraction in response to binding noradrenaline. NMR of [13CϵH3]methionine-labeled α1A-adrenoreceptor variants, each exhibiting differing signaling capacities, revealed how different classes of ligands modulate the conformational equilibria of this receptor. [13CϵH3]Methionine residues near the microswitches exhibited distinct states that correlated with ligand efficacies, supporting a conformational selection mechanism. We propose that allosteric coupling among the microswitches controls the conformation of the α1A-adrenoreceptor and underlies the mechanism of ligand modulation of GPCR signaling in cells.

Structure of formylpeptide receptor 2-Gi complex reveals insights into ligand recognition and signaling

Formylpeptide receptors (FPRs) as G protein-coupled receptors (GPCRs) can recognize formylpeptides derived from pathogens or host cells to function in host defense and cell clearance. In addition, FPRs, especially FPR2, can also recognize other ligands with a large chemical diversity generated at different stages of inflammation to either promote or resolve inflammation in order to maintain a balanced inflammatory response. The mechanism underlying promiscuous ligand recognition and activation of FPRs is not clear. Here we report a cryo-EM structure of FPR2-Gi signaling complex with a peptide agonist. The structure reveals a widely open extracellular region with an amphiphilic environment for ligand binding. Together with computational docking and simulation, the structure suggests a molecular basis for the recognition of formylpeptides and a potential mechanism of receptor activation, and reveals conserved and divergent features in Gi coupling. Our results provide a basis for understanding the molecular mechanism of the functional promiscuity of FPRs.

Bridging non-overlapping reads illuminates high-order epistasis between distal protein sites in a GPCR

Epistasis emerges when the effects of an amino acid depend on the identities of interacting residues. This phenomenon shapes fitness landscapes, which have the power to reveal evolutionary paths and inform evolution of desired functions. However, there is a need for easily implemented, high-throughput methods to capture epistasis particularly at distal sites. Here, we combine deep mutational scanning (DMS) with a straightforward data processing step to bridge reads in distal sites within genes (BRIDGE). We use BRIDGE, which matches non-overlapping reads to their cognate templates, to uncover prevalent epistasis within the binding pocket of a human G protein-coupled receptor (GPCR) yielding variants with 4-fold greater affinity to a target ligand. The greatest functional improvements in our screen result from distal substitutions and substitutions that are deleterious alone. Our results corroborate findings of mutational tolerance in GPCRs, even in conserved motifs, but reveal inherent constraints restricting tolerated substitutions due to epistasis.

Epistasis effects among amino acids at distal sites within binding pockets can have important impacts on protein fitness landscapes. Here the authors present BRIDGE, which matches non-overlapping sequence reads with their cognate DNA templates.

Conformational Differences among Metarhodopsin I, Metarhodopsin II, and Opsin Probed by Wide-Angle X-ray Scattering

Among the photoproducts of vertebrate rhodopsin, only metarhodopsin II (Meta-II) preferentially adopts the active structure in which transmembrane helices are rearranged. Light-induced helical rearrangement of rhodopsin in membrane-embedded form was directly monitored by wide-angle X-ray scattering (WAXS) using nanodiscs. The change in the WAXS curve for the formation of Meta-II was characterized by a peak at 0.2 Å^-1 and a valley at 0.6 Å^-1, which were not observed in metarhodopsin I and opsin. However, acid-induced active opsin (Opsin*) showed a 0.2 Å^-1 peak, but no 0.6 Å^-1 valley. Analyses using the model structures based on the crystal structures of dark state and Meta-II suggest that the outward movement of helix VI occurred in Opsin*. However, the displaced helices III and V in Meta-II resulting from the disruption of cytoplasmic ionic lock were restored in Opsin*, which is likely to destabilize the G-protein-activating structure of opsin.

Toward Understanding the Impact of Dimerization Interfaces in Angiotensin II Type 1 Receptor

Angiotensin II type 1 receptor (AT1R) is a prototypical class A G protein-coupled receptor (GPCR) that has an important role in cardiovascular pathologies and blood pressure regulation as well as in the central nervous system. GPCRs may exist and function as monomers; however, they can assemble to form higher order structures, and as a result of oligomerization, their function and signaling profiles can be altered. In the case of AT1R, the classical Gα_q/11 pathway is initiated with endogenous agonist angiotensin II binding. A variety of cardiovascular pathologies such as heart failure, diabetic nephropathy, atherosclerosis, and hypertension are associated with this pathway. Recent findings reveal that AT1R can form homodimers and activate the noncanonical (β-arrestin-mediated) pathway. Nevertheless, the exact dimerization interface and atomic details of AT1R homodimerization have not been still elucidated. Here, six different symmetrical dimer interfaces of AT1R are considered, and homodimers were constructed using other published GPCR crystal dimer interfaces as template structures. These AT1R homodimers were then inserted into the model membrane bilayers and subjected to all-atom molecular dynamics simulations. Our simulation results along with the principal component analysis and water pathway analysis suggest four different interfaces as the most plausible: symmetrical transmembrane (TM)1,2,8; TM5; TM4; and TM4,5 AT1R dimer interfaces that consist of one inactive and one active protomer. Moreover, we identified ILE238^6.33 as a hub residue in the stabilization of the inactive state of AT1R.

Molecular Determinants of the Intrinsic Efficacy of the Antipsychotic Aripiprazole

Partial agonists of the dopamine D₂ receptor (D₂R) have been developed to treat the symptoms of schizophrenia without causing the side effects elicited by antagonists. The receptor-ligand interactions that determine the intrinsic efficacy of such drugs, however, are poorly understood. Aripiprazole has an extended structure comprising a phenylpiperazine primary pharmacophore and a 1,2,3,4-tetrahydroquinolin-2-one secondary pharmacophore. We combined site-directed mutagenesis, analytical pharmacology, ligand fragments, and molecular dynamics simulations to identify the D₂R-aripiprazole interactions that contribute to affinity and efficacy. We reveal that an interaction between the secondary pharmacophore of aripiprazole and a secondary binding pocket defined by residues at the extracellular portions of transmembrane segments 1, 2, and 7 determines the intrinsic efficacy of aripiprazole. Our findings reveal a hitherto unappreciated mechanism for fine-tuning the intrinsic efficacy of D₂R agonists.

How Does Agonist and Antagonist Binding Lead to Different Conformational Ensemble Equilibria of the κ-Opioid Receptor: Insight from Long-Time Gaussian Accelerated Molecular Dynamics Simulation

The opioid receptors belong to the class A seven transmembrane-spanning (7TM) G protein-coupled receptors (GPCRs). The κ-opioid receptor (KOR) is a subfamily of four opioid receptors. The endogenous peptide and a variety of selective agonists and antagonists of KOR have been developed. The structurally similar ligands at the same site cause completely opposite biological functions and induce different conformational changes. To shed light on the conformation ensembles and conformational dynamics in activation and deactivation processes of KOR, we performed all-atom, long-time Gaussian accelerated molecular dynamics simulation (GaMD) on KOR binding with agonist epoxymorphinan MP1104 and antagonist JDTic, respectively. Our results revealed different conformation ensembles of KOR binding with agonist and with antagonist. Agonist binding stabilizes the active state of key motifs including DYYNM motif and CWxP motif, and biases the conformation equilibria toward the active state. Antagonist binding will not destroy inactive conformation equilibria, by keeping the stable inactive state of these crucial motifs. We found that the inactive apo form of KOR is the most stable state, while the active apo form relaxes readily to inactive state. Our results also revealed a stable intermediate (I), which is attributed to the hydrophobic interactions between Tyr246^5.58 and TM6, as well as the steric hindrance of them. Our results not only show the conformation equilibria bias of KOR by binding with agonist and antagonist, but also provide the structural information for the design and discovery of potential ligands with different functions.

Glu2.53(90) of the GnRH receptor is part of the conserved G protein-coupled receptor structure and does not form a salt-bridge with Lys3.32(121)

GnRH receptor mutations, Glu^2.53(90)Lys and Glu^2.53(90)Asp, cause congenital hypogonadotropic hypogonadism. The Glu^2.53(90) side-chain has been proposed to form an intramolecular salt-bridge with Lys^3.32(121), but conserved intramolecular interaction networks in G protein-coupled receptor crystal structures predict that it interacts with Ser^3.35(124) and Trp^6.48(280). We investigated interhelical interactions of Glu^2.53(90) that stabilise GnRH receptor folding using functional analyses and computational modelling of mutant receptors. The Glu^2.53(90)Asp mutant was non-functional, but mutants with hydrophobic amino acids or Arg substituted for Glu^2.53(90) were functional, excluding a salt-bridge interaction. The Glu^2.53(90)Arg and Trp^6.48(280)Arg mutants had decreased affinity for GnRH. Models showed that congenital Glu^2.53(90)Lys and Glu^2.53(90)Asp mutations disrupt interactions with Ser^3.35(124) and Trp^6.48(280) respectively, whereas the Glu^2.53(90)Arg and Trp^6.48(280)Arg mutations preserve intramolecular contacts, but increase distance between the transmembrane helices. Our results show that disruption of interhelical contacts that are conserved in G protein-coupled receptors accounts for the effects of some disease-associated GnRH receptor mutations.

Arrestin-1 engineering facilitates complex stabilization with native rhodopsin

Arrestin-1 desensitizes the activated and phosphorylated photoreceptor rhodopsin by forming transient rhodopsin−arrestin-1 complexes that eventually decay to opsin, retinal and arrestin-1. Via a multi-dimensional screening setup, we identified and combined arrestin-1 mutants that form lasting complexes with light-activated and phosphorylated rhodopsin in harsh conditions, such as high ionic salt concentration. Two quadruple mutants, D303A + T304A + E341A + F375A and R171A + T304A + E341A + F375A share similar heterologous expression and thermo-stability levels with wild type (WT) arrestin-1, but are able to stabilize complexes with rhodopsin with more than seven times higher half-maximal inhibitory concentration (IC₅₀) values for NaCl compared to the WT arrestin-1 protein. These quadruple mutants are also characterized by higher binding affinities to phosphorylated rhodopsin, light-activated rhodopsin and phosphorylated opsin, as compared with WT arrestin-1. Furthermore, the assessed arrestin-1 mutants are still specifically associating with phosphorylated or light-activated receptor states only, while binding to the inactive ground state of the receptor is not significantly altered. Additionally, we propose a novel functionality for R171 in stabilizing the inactive arrestin-1 conformation as well as the rhodopsin–arrestin-1 complex. The achieved stabilization of the active rhodopsin–arrestin-1 complex might be of great interest for future structure determination, antibody development studies as well as drug-screening efforts targeting G protein-coupled receptors (GPCRs).

Kinetic and thermodynamic insights into sodium ion translocation through the μ-opioid receptor from molecular dynamics and machine learning analysis

The differential modulation of agonist and antagonist binding to opioid receptors (ORs) by sodium (Na+) has been known for decades. To shed light on the molecular determinants, thermodynamics, and kinetics of Na+ translocation through the μ-OR (MOR), we used a multi-ensemble Markov model framework combining equilibrium and non-equilibrium atomistic molecular dynamics simulations of Na+ binding to MOR active or inactive crystal structures embedded in an explicit lipid bilayer. We identify an energetically favorable, continuous ion pathway through the MOR active conformation only, and provide, for the first time: i) estimates of the energy differences and required timescales of Na+ translocation in inactive and active MORs, ii) estimates of Na+-induced changes to agonist binding validated by radioligand measurements, and iii) testable hypotheses of molecular determinants and correlated motions involved in this translocation, which are likely to play a key role in MOR signaling.