Effect of macrocyclization and tetramethylrhodamine labeling on chemokine binding peptides.

Receptor-derived peptides have played an important role in elucidating chemokine-receptor interactions. For the inflammatory chemokine CXC-class chemokine ligand 8 (CXCL8), a site II-mimetic peptide has been derived from parts of extracellular loops 2 and 3 and adjacent transmembrane helices of its receptor CXC-class chemokine receptor 1 (Helmer et al., RSC Adv., 2015, 5, 25657). The peptide sequence with a C-terminal glutamine did not bind to CXCL8, whereas one with a C-terminal glutamate did but with low micromolar affinity. We sought to improve the affinity and protease stability of the latter peptide through cyclization while also cyclizing the former for control purposes. To identify a cyclization strategy that permits a receptor-like interaction, we conducted a molecular dynamics simulation of CXCL8 in complex with full-length CXC-class chemokine receptor 1. We introduced a linker to provide an appropriate spacing between the termini and used an on-resin side-chain-to-tail cyclization strategy. Upon chemokine binding, the fluorescence intensity of the tetramethylrhodamine (TAMRA)-labeled cyclic peptides increased whereas the fluorescence anisotropy decreased. Additional molecular dynamics simulations indicated that the fluorophore interacts with the peptide macrocycle so that chemokine binding leads to its displacement and observed changes in fluorescence. Macrocyclization of both 18-amino acid-long peptides led to the same low micromolar affinity for CXCL8. Likewise, both TAMRA-labeled linear peptides interacted with CXCL8 with similar affinities. Interestingly, the linear TAMRA-labeled peptides were more resistant to tryptic digestion than the unlabeled counterparts, whereas the cyclized peptides were not degraded at all. We conclude that the TAMRA fluorophore tends to interact with peptides altering their protease stability and behavior in fluorescence-based assays.

Discovery of a true bivalent dopamine D2 receptor agonist.

Employing two different alkyne-modified dopamine agonists to construct bivalent compounds via click chemistry resulted in the identification of a bivalent ligand (11c) for dopamine D2 receptor homodimer, which, compared to its parent monomeric alkyne, showed a 16-fold higher binding affinity for the dopamine D2 receptor and a 5-fold higher potency in a cAMP assay in HEK 293T cells stably expressing D2R. Molecular modeling revealed that 11c can indeed bridge the orthosteric binding sites of a D2R homodimer in a relaxed conformation via the TM5-TM6 interface and allows to largely rationalize the results of the receptor assays.

Structural Assessment of Agonist Efficacy in the µ-Opioid Receptor: Morphine and Fentanyl Elicit Different Activation Patterns.

Over the past two decades, the opioid epidemic in the United States and Canada has evidenced the need for a better understanding of the molecular mechanisms of medications used to fight pain. Morphine and fentanyl are widely used in opiate-mediated analgesia for the treatment of chronic pain. These compounds target the µ-opioid receptor (MOR), a class A G protein-coupled receptor (GPCR). In light of described higher efficacy of fentanyl with respect to morphine, we have performed independent µs-length unbiased molecular dynamics (MD) simulations of MOR complexes with each of these ligands, including the MOR antagonist naltrexone as a negative control. Consequently, MD simulations totaling 58 µs have been conducted to elucidate at the atomic level ligand-specific receptor activity and signal transmission in the MOR. In particular, we have identified stable binding poses of morphine and fentanyl, which interact differently with the MOR. Different ligand-receptor interaction landscapes directly induce sidechain conformational changes of orthosteric pocket residues: Asp1493.32, Tyr1503.33, Gln1262.60, and Lys2355.39. The induced conformations determine Asp1493.32-Tyr3287.43 sidechain-sidechain interactions and Trp2956.48-Ala2425.46 sidechain-backbone H-bond formations, as well as Met1533.36 conformational changes. In addition to differences in ligand binding, different intracellular receptor conformational changes are observed as morphine preferentially activates transmembrane (TM) helices: TM3 and TM5, while fentanyl preferentially activates TM6 and TM7. As conformational changes in TM6 and TM7 are widely described as being the most crucial aspect in GPCR activation, this may contribute to the greater efficacy of fentanyl over morphine. These computationally observed functional differences between fentanyl and morphine may provide new avenues for the design of safer but not weaker opioid drugs because it is desirable to increase the safety of medicines without sacrificing their efficacy.

Statistics for the analysis of molecular dynamics simulations: providing P values for agonist-dependent GPCR activation.

Molecular dynamics (MD) is the common computational technique for assessing efficacy of GPCR-bound ligands. Agonist efficacy measures the capability of the ligand-bound receptor of reaching the active state in comparison with the free receptor. In this respect, agonists, neutral antagonists and inverse agonists can be considered. A collection of MD simulations of both the ligand-bound and the free receptor are needed to provide reliable conclusions. Variability in the trajectories needs quantification and proper statistical tools for meaningful and non-subjective conclusions. Multiple-factor (time, ligand, lipid) ANOVA with repeated measurements on the time factor is proposed as a suitable statistical method for the analysis of agonist-dependent GPCR activation MD simulations. Inclusion of time factor in the ANOVA model is consistent with the time-dependent nature of MD. Ligand and lipid factors measure agonist and lipid influence on receptor activation. Previously reported MD simulations of adenosine A2a receptor (A2aR) are reanalyzed with this statistical method. TM6-TM3 and TM7-TM3 distances are selected as dependent variables in the ANOVA model. The ligand factor includes the presence or absence of adenosine whereas the lipid factor considers DOPC or DOPG lipids. Statistical analysis of MD simulations shows the efficacy of adenosine and the effect of the membrane lipid composition. Subsequent application of the statistical methodology to NECA A2aR agonist, with resulting P values in consistency with its pharmacological profile, suggests that the method is useful for ligand comparison and potentially for dynamic structure-based virtual screening.

Insights into adenosine A2A receptor activation through cooperative modulation of agonist and allosteric lipid interactions.

The activation process of G protein-coupled receptors (GPCRs) has been extensively studied, both experimentally and computationally. In particular, Molecular Dynamics (MD) simulations have proven useful in exploring GPCR conformational space. The typical behaviour of class A GPCRs, when subjected to unbiased MD simulations from their crystallized inactive state, is to fluctuate between inactive and intermediate(s) conformations, even with bound agonist. Fully active conformation(s) are rarely stabilized unless a G protein is also bound. Despite several crystal structures of the adenosine A2a receptor (A2aR) having been resolved in complex with co-crystallized agonists and Gs protein, its agonist-mediated activation process is still not completely understood. In order to thoroughly examine the conformational landscape of A2aR activation, we performed unbiased microsecond-length MD simulations in quadruplicate, starting from the inactive conformation either in apo or with bound agonists: endogenous adenosine or synthetic NECA, embedded in two homogeneous phospholipid membranes: 1,2-dioleoyl-sn-glycerol-3-phosphoglycerol (DOPG) or 1,2-dioleoyl-sn-glycerol-3-phosphocholine (DOPC). In DOPC with bound adenosine or NECA, we observe transition to an intermediate receptor conformation consistent with the known adenosine-bound crystal state. In apo state in DOPG, two different intermediate conformations are obtained. One is similar to that observed with bound adenosine in DOPC, while the other is closer to the active state but not yet fully active. Exclusively, in DOPG with bound adenosine or NECA, we reproducibly identify receptor conformations with fully active features, which are able to dock Gs protein. These different receptor conformations can be attributed to the action/absence of agonist and phospholipid-mediated allosteric effects on the intracellular side of the receptor.

Exploring the Activation Mechanism of the mGlu5 Transmembrane Domain.

As a class C GPCR and regulator of synaptic activity, mGlu5 is an attractive drug target, potentially offering treatment for several neurologic and psychiatric disorders. As little is known about the activation mechanism of mGlu5 at a structural level, potential of mean force calculations linked to molecular dynamics simulations were performed on the mGlu5 transmembrane domain crystal structure to explore various internal mechanisms responsible for its activation. Our results suggest that the hydrophilic interactions between intracellular loop 1 and the intracellular side of TM6 have to be disrupted to reach a theoretically active-like conformation. In addition, interactions between residues that are key for mGlu5 activation (Tyr6593.44 and Ile7515.51) and mGlu5 inactivation (Tyr6593.44 and Ser8097.39) have been identified. Inasmuch as mGlu5 receptor signaling is poorly understood, potentially showing a complex network of micro-switches and subtle structure-activity relationships, the present study represents a step forward in the understanding of mGlu5 transmembrane domain activation.

Artificial Intelligence: A Novel Approach for Drug Discovery.

Molecular dynamics (MD) simulations can mechanistically explain receptor function. However, the enormous data sets that they may imply can be a hurdle. Plante and colleagues (Molecules, 2019) recently described a machine learning approach to the analysis of MD simulations. The approach successfully classified ligands and identified functional receptor motifs and thus it seems promising for mechanism-based drug discovery.

Revealing the Mechanism of Agonist-Mediated Cannabinoid Receptor 1 (CB1) Activation and Phospholipid-Mediated Allosteric Modulation.

Cannabinoid receptor 1 (CB1) mediates the functional responses of Δ9-tetrahydrocannabinol. Although progress has been made in understanding cannabinoid binding and receptor activation, detailed knowledge of the dynamics involved in the activation mechanism of CB1 is lacking. Here, we use recently determined CB1 crystal structures to analyze its transition from inactive to active state by performing unbiased microsecond-length molecular dynamics (MD) simulations, totaling 32 µs, with and without bound potent cannabinoid agonist CP-55940. CB1 activation is characterized by an upward axial movement of transmembrane (TM) helix 3, inward movement of TM7, and outward movement of TM6. These conformational changes collectively allow Gi protein docking, although fully active states of the receptor occur only transiently during MD simulations. Additionally, positive allosteric modulation of CB1 by anionic phospholipids is found to increase action of the bound agonist. Specifically, this involves protein-lipid interactions at intracellular loop 3, TM6, and ionic lock residue Arg2143.50.

Synthesis toward Bivalent Ligands for the Dopamine D2 and Metabotropic Glutamate 5 Receptors.

In this study, we designed and synthesized heterobivalent ligands targeting heteromers consisting of the metabotropic glutamate 5 receptor (mGluR5) and the dopamine D2 receptor (D2R). Bivalent ligand 22a with a linker consisting of 20 atoms showed 4-fold increase in affinity for cells coexpressing D2R and mGluR5 compared to cells solely expressing D2R. Likewise, the affinity of 22a for mGluR5 increased 2-fold in the coexpressing cells. Additionally, 22a exhibited a 5-fold higher mGluR5 affinity than its monovalent precursor 21a in cells coexpressing D2R and mGluR5. These results indicate that 22a is able to bridge binding sites on both receptors constituting the heterodimer. Likewise, cAMP assays revealed that 22a had a 4-fold higher potency in stable D2R and mGluR5 coexpressing cell lines than 1. Furthermore, molecular modeling reveals that 22a is able to simultaneously bind both receptors by passing between the TM5-TM6 interface and establishing six protein-ligand H-bonds.

Structural insights into positive and negative allosteric regulation of a G protein-coupled receptor through protein-lipid interactions.

Lipids are becoming known as essential allosteric modulators of G protein-coupled receptor (GPCRs). However, how they exert their effects on GPCR conformation at the atomic level is still unclear. In light of recent experimental data, we have performed several long-timescale molecular dynamics (MD) simulations, totalling 24 µs, to rigorously map allosteric modulation and conformational changes in the ß2 adrenergic receptor (ß2AR) that occur as a result of interactions with three different phospholipids. In particular, we identify different sequential mechanisms behind receptor activation and deactivation, respectively, mediated by specific lipid interactions with key receptor regions. We show that net negatively charged lipids stabilize an active-like state of ß2AR that is able to dock Gsα protein. Clustering of anionic lipids around the receptor with local distortion of membrane thickness is also apparent. On the other hand, net-neutral zwitterionic lipids inactivate the receptor, generating either fully inactive or intermediate states, with kinetics depending on lipid headgroup charge distribution and hydrophobicity. These chemical differences alter membrane thickness and density, which differentially destabilize the ß2AR active state through lateral compression effects.

Correction: Rational design of a peptide capture agent for CXCL8 based on a model of the CXCL8:CXCR1 complex.

[This corrects the article DOI: 10.1039/C4RA13749C.].

Analysis of positive and negative allosteric modulation in metabotropic glutamate receptors 4 and 5 with a dual ligand.

As class C GPCRs and regulators of synaptic activity, human metabotropic glutamate receptors (mGluRs) 4 and 5 are prime targets for allosteric modulation, with mGlu5 inhibition or mGlu4 stimulation potentially treating conditions like chronic pain and Parkinson's disease. As an allosteric modulator that can bind both receptors, 2-Methyl-6-(phenylethynyl)pyridine (MPEP) is able to negatively modulate mGlu5 or positively modulate mGlu4. At a structural level, how it elicits these responses and how mGluRs undergo activation is unclear. Here, we employ homology modelling and 30 µs of atomistic molecular dynamics (MD) simulations to probe allosteric conformational change in mGlu4 and mGlu5, with and without docked MPEP. Our results identify several structural differences between mGlu4 and mGlu5, as well as key differences responsible for MPEP-mediated positive and negative allosteric modulation, respectively. A novel mechanism of mGlu4 activation is revealed, which may apply to all mGluRs in general. This involves conformational changes in TM3, TM4 and TM5, separation of intracellular loop 2 (ICL2) from ICL1/ICL3, and destabilization of the ionic-lock. On the other hand, mGlu5 experiences little disturbance when MPEP binds, maintaining its inactive state with reduced conformational fluctuation. In addition, when MPEP is absent, a lipid molecule can enter the mGlu5 allosteric pocket.

Angiotensin II type 1/adenosine A 2A receptor oligomers: a novel target for tardive dyskinesia.

Tardive dyskinesia (TD) is a serious motor side effect that may appear after long-term treatment with neuroleptics and mostly mediated by dopamine D2 receptors (D2Rs). Striatal D2R functioning may be finely regulated by either adenosine A2A receptor (A2AR) or angiotensin receptor type 1 (AT1R) through putative receptor heteromers. Here, we examined whether A2AR and AT1R may oligomerize in the striatum to synergistically modulate dopaminergic transmission. First, by using bioluminescence resonance energy transfer, we demonstrated a physical AT1R-A2AR interaction in cultured cells. Interestingly, by protein-protein docking and molecular dynamics simulations, we described that a stable heterotetrameric interaction may exist between AT1R and A2AR bound to antagonists (i.e. losartan and istradefylline, respectively). Accordingly, we subsequently ascertained the existence of AT1R/A2AR heteromers in the striatum by proximity ligation in situ assay. Finally, we took advantage of a TD animal model, namely the reserpine-induced vacuous chewing movement (VCM), to evaluate a novel multimodal pharmacological TD treatment approach based on targeting the AT1R/A2AR complex. Thus, reserpinized mice were co-treated with sub-effective losartan and istradefylline doses, which prompted a synergistic reduction in VCM. Overall, our results demonstrated the existence of striatal AT1R/A2AR oligomers with potential usefulness for the therapeutic management of TD.

Illuminating Phenylazopyridines To Photoswitch Metabotropic Glutamate Receptors: From the Flask to the Animals.

Phenylazopyridines are photoisomerizable compounds with high potential to control biological functions with light. We have obtained a series of phenylazopyridines with light dependent activity as negative allosteric modulators (NAM) of metabotropic glutamate receptor subtype 5 (mGlu5). Here we describe the factors needed to achieve an operational molecular photoisomerization and its effective translation into in vitro and in vivo receptor photoswitching, which includes zebrafish larva motility and the regulation of the antinociceptive effects in mice. The combination of light and some specific phenylazopyridine ligands displays atypical pharmacological profiles, including light-dependent receptor overactivation, which can be observed both in vitro and in vivo. Remarkably, the localized administration of light and a photoswitchable compound in the peripheral tissues of rodents or in the brain amygdalae results in an illumination-dependent analgesic effect. The results reveal a robust translation of the phenylazopyridine photoisomerization to a precise photoregulation of biological activity.

Positional isomers of bispyridine benzene derivatives induce efficacy changes on mGlu5 negative allosteric modulation.

Modulation of metabotropic glutamate receptor 5 (mGlu5) with partial allosteric antagonists has received increased interest due to their favourable in vivo activity profiles compared to the unfavourable side-effects of full inverse agonists. Here we report on a series of bispyridine benzene derivatives with a functional molecular switch affecting antagonistic efficacy, shifting from inverse agonism to partial antagonism with only a single change in the substitution pattern of the benzene ring. These efficacy changes are explained through computational docking, revealing two different receptor conformations of different energetic stability and different positional isomer binding preferences.

Shining Light on an mGlu5 Photoswitchable NAM: A Theoretical Perspective.

Metabotropic glutamate receptors (mGluRs) are important drug targets because of their involvement in several neurological diseases. Among mGluRs, mGlu5 is a particularly high-profile target because its positive or negative allosteric modulation can potentially treat schizophrenia or anxiety and chronic pain, respectively. Here, we computationally and experimentally probe the functional binding of a novel photoswitchable mGlu5 NAM, termed alloswitch-1, which loses its NAM functionality under violet light. We show alloswitch-1 binds deep in the allosteric pocket in a similar fashion to mavoglurant, the co-crystallized NAM in the mGlu5 transmembrane domain crystal structure. Alloswitch-1, like NAM 2-Methyl-6-(phenylethynyl)pyridine (MPEP), is significantly affected by P655M mutation deep in the allosteric pocket, eradicating its functionality. In MD simulations, we show alloswitch-1 and MPEP stabilize the co-crystallized water molecule located at the bottom of the allosteric site that is seemingly characteristic of the inactive receptor state. Furthermore, both NAMs form H-bonds with S809 on helix 7, which may constitute an important stabilizing interaction for NAM-induced mGlu5 inactivation. Alloswitch-1, through isomerization of its amide group from trans to cis is able to form an additional interaction with N747 on helix 5. This may be an important interaction for amide-containing mGlu5 NAMs, helping to stabilize their binding in a potentially unusual cis-amide state. Simulated conformational switching of alloswitch-1 in silico suggests photoisomerization of its azo group from trans to cis may be possible within the allosteric pocket. However, photoexcited alloswitch-1 binds in an unstable fashion, breaking H-bonds with the protein and destabilizing the co-crystallized water molecule. This suggests photoswitching may have destabilizing effects on mGlu5 binding and functionality.

Helix 3 acts as a conformational hinge in Class A GPCR activation: An analysis of interhelical interaction energies in crystal structures.

A collection of crystal structures of rhodopsin, ß2-adrenergic and adenosine A2A receptors in active, intermediate and inactive states were selected for structural and energetic analyses to identify the changes involved in the activation/deactivation of Class A GPCRs. A set of helix interactions exclusive to either inactive or active/intermediate states were identified. The analysis of these interactions distinguished some local conformational changes involved in receptor activation, in particular, a packing between the intracellular domains of transmembrane helices H3 and H7 and a separation between those of H2 and H6. Also, differential movements of the extracellular and intracellular domains of these helices are apparent. Moreover, a segment of residues in helix H3, including residues L/I3.40 to L3.43, is identified as a key component of the activation mechanism, acting as a conformational hinge between extracellular and intracellular regions. Remarkably, the influence on the activation process of some glutamic and aspartic acidic residues and, as a consequence, the influence of variations on local pH is highlighted. Structural hypotheses that arose from the analysis of rhodopsin, ß2-adrenergic and adenosine A2A receptors were tested on the active and inactive M2 muscarinic acetylcholine receptor structures and further discussed in the context of the new mechanistic insights provided by the recently determined active and inactive crystal structures of the µ-opioid receptor. Overall, the structural and energetic analyses of the interhelical interactions present in this collection of Class A GPCRs suggests the existence of a common general activation mechanism featuring a chemical space useful for drug discovery exploration.

Selective Protonation of Acidic Residues Triggers Opsin Activation.

Rhodopsin, the visual pigment in the retina, is a Class A G protein-coupled receptor (GPCR) covalently bound to retinal chromophore. In dark conditions, retinal is in the cis-isomeric state, stabilizing the rhodopsin inactive state as an inverse agonist. After light absorption, retinal undergoes an isomerization photoreaction to trans-retinal, which includes a conformational change of the receptor to its active state. In the absence of retinal, the apoprotein opsin presents a low level of constitutive activity, which depends on pH, with higher propensity of activation at acidic pH. To examine the effect and the underlying mechanism that protonation may have on opsin activation, a number of MD simulations were run varying the number and identity of acidic residues selected for protonation. Results show that the combined protonation of D83, E113, and E247 is of special relevance for the induction of receptor activation. Subsequent conformational analysis of the MD trajectories provides a structural mechanistic insight into the opsin activation process. Furthermore, because protonation seems to be a determining step in the activation of other GPCRs, the methodology and rationale used herein can be extended to mechanistic studies of GPCRs in general.

Quantifying conformational changes in GPCRs: glimpse of a common functional mechanism.

BACKGROUND: G-protein-coupled receptors (GPCRs) are important drug targets and a better understanding of their molecular mechanisms would be desirable. The crystallization rate of GPCRs has accelerated in recent years as techniques have become more sophisticated, particularly with respect to Class A GPCRs interacting with G-proteins. These developments have made it possible for a quantitative analysis of GPCR geometrical features and binding-site conformations, including a statistical comparison between Class A GPCRs in active (agonist-bound) and inactive (antagonist-bound) states. RESULTS: Here we implement algorithms for the analysis of interhelical angles, distances, interactions and binding-site volumes in the transmembrane domains of 25 Class A GPCRs (7 active and 18 inactive). Two interhelical angles change in a statistically significant way between average inactive and active states: TM3-TM6 (by -9°) and TM6-TM7 (by +12°). A third interhelical angle: TM5-TM6 shows a trend, changing by -9°. In the transition from inactive to active states, average van der Waals interactions between TM3 and TM7 significantly increase as the average distance between them decreases by >2 Å. Average H-bonding between TM3 and TM6 decreases but is seemingly compensated by an increase in H-bonding between TM5 and TM6. In five Class A GPCRs, crystallized in both active and inactive states, increased H-bonding of agonists to TM6 and TM7, relative to antagonists, is observed. These protein-agonist interactions likely favour a change in the TM6-TM7 angle, which creates a narrowing in the binding pocket of activated receptors and an average ~200 Å(3) reduction in volume. CONCLUSIONS: In terms of similar conformational changes and agonist binding pattern, Class A GPCRs appear to share a common mechanism of activation, which can be exploited in future drug development.

Computational analysis of negative and positive allosteric modulator binding and function in metabotropic glutamate receptor 5 (in)activation.

Metabotropic glutamate receptors (mGluRs) are high-profile G-protein coupled receptors drug targets because of their involvement in several neurological disease states, and mGluR5 in particular is a subtype whose controlled allosteric modulation, both positive and negative, can potentially be useful for the treatment of schizophrenia and relief of chronic pain, respectively. Here we model mGluR5 with a collection of positive and negative allosteric modulators (PAMs and NAMs) in both active and inactive receptor states, in a manner that is consistent with experimental information, using a specialized protocol that includes homology to increase docking accuracy, and receptor relaxation to generate an individual induced fit with each allosteric modulator. Results implicate two residues in particular for NAM and PAM function: NAM interaction with W785 for receptor inactivation, and NAM/PAM H-bonding with S809 for receptor (in)activation. Models suggest the orientation of the H-bond between allosteric modulator and S809, controlled by PAM/NAM chemistry, influences the position of TM7, which in turn influences the shape of the allosteric site, and potentially the receptor state. NAM-bound and PAM-bound mGluR5 models also reveal that although PAMs and NAMs bind in the same pocket and share similar binding modes, they have distinct effects on the conformation of the receptor. Our models, together with the identification of a possible activation mechanism, may be useful in the rational design of new allosteric modulators for mGluR5.