G protein-specific mechanisms in the serotonin 5-HT2A receptor regulate psychosis-related effects and memory deficits.

G protein-coupled receptors (GPCRs) are sophisticated signaling machines able to simultaneously elicit multiple intracellular signaling pathways upon activation. Complete (in)activation of all pathways can be counterproductive for specific therapeutic applications. This is the case for the serotonin 2 A receptor (5-HT2AR), a prominent target for the treatment of schizophrenia. In this study, we elucidate the complex 5-HT2AR coupling signature in response to different signaling probes, and its physiological consequences by combining computational modeling, in vitro and in vivo experiments with human postmortem brain studies. We show how chemical modification of the endogenous agonist serotonin dramatically impacts the G protein coupling profile of the 5-HT2AR and the associated behavioral responses. Importantly, among these responses, we demonstrate that memory deficits are regulated by Gαq protein activation, whereas psychosis-related behavior is modulated through Gαi1 stimulation. These findings emphasize the complexity of GPCR pharmacology and physiology and open the path to designing improved therapeutics for the treatment of stchizophrenia.

Do they make a good match? Molecular dynamics studies on CALB-catalyzed esterification of 3-phenylpropionic and cinnamic acids.

Lipases are versatile catalysts widely used in industrial biotransformations and laboratory-scale developed reactions with industrial potential. Despite the fact that lipase B from Candida antarctica (CALB) is one of the most widely used lipolytic enzymes, its substrate specificity is still poorly understood. One observed trend is that reactions carried out with carboxylic acids containing a double bond are less efficient on average. Here, we have utilized a combination of in vitro and in silico techniques, to better understand the negative impact of a double bond on CALB-mediated esterification. Then through extensive molecular dynamics (MD) simulations, we were able to map the entry pathway of cinnamic acid and its derivative into the CALB active site, and their interactions with catalytic residues. We observed a 2 step binding mechanism of studied compounds, where they first penetrate the enzyme pocket in a conformation where their carboxylic groups are extended towards the solvent. This is followed by further penetration of the acid into the enzymatic active pocket, and a full rotation within the active site, which orients the acid in a conformation that allows further steps of the esterification reaction. As acids containing a double bond are more rigid, their mobility and thus ability to rotate in the narrow CALB active site is hampered, which provides a structural explanation for the decreased efficiency of such acids. Our data provide insight into the substrate specificity of CALB-mediated esterification, providing important structural guidelines to better understand and potentially improve the efficiency of such reactions.

Ligand entry pathways control the chemical space recognized by GPR183.

The G protein-coupled receptor GPR183 is a chemotactic receptor with an important function in the immune system and association with a variety of diseases. It recognizes ligands with diverse physicochemical properties as both the endogenous oxysterol ligand 7α,25-OHC and synthetic molecules can activate the G protein pathway of the receptor. To better understand the ligand promiscuity of GPR183, we utilized both molecular dynamics simulations and cell-based validation experiments. Our work reveals that the receptor possesses two ligand entry channels: one lateral between transmembrane helices 4 and 5 facing the membrane, and one facing the extracellular environment. Using enhanced sampling, we provide a detailed structural model of 7α,25-OHC entry through the lateral membrane channel. Importantly, the first ligand recognition point at the receptor surface has been captured in diverse experimentally solved structures of different GPCRs. The proposed ligand binding pathway is supported by in vitro data employing GPR183 mutants with a sterically blocked lateral entrance, which display diminished binding and signaling. In addition, computer simulations and experimental validation confirm the existence of a polar water channel which might serve as an alternative entrance gate for less lipophilic ligands from the extracellular milieu. Our study reveals knowledge to understand GPR183 functionality and ligand recognition with implications for the development of drugs for this receptor. Beyond, our work provides insights into a general mechanism GPCRs may use to respond to chemically diverse ligands.

Phosphorylation barcodes direct biased chemokine signaling at CXCR3.

G protein-coupled receptor (GPCR)-biased agonism, selective activation of certain signaling pathways relative to others, is thought to be directed by differential GPCR phosphorylation "barcodes." At chemokine receptors, endogenous chemokines can act as "biased agonists", which may contribute to the limited success when pharmacologically targeting these receptors. Here, mass spectrometry-based global phosphoproteomics revealed that CXCR3 chemokines generate different phosphorylation barcodes associated with differential transducer activation. Chemokine stimulation resulted in distinct changes throughout the kinome in global phosphoproteomics studies. Mutation of CXCR3 phosphosites altered ß-arrestin 2 conformation in cellular assays and was consistent with conformational changes observed in molecular dynamics simulations. T cells expressing phosphorylation-deficient CXCR3 mutants resulted in agonist- and receptor-specific chemotactic profiles. Our results demonstrate that CXCR3 chemokines are non-redundant and act as biased agonists through differential encoding of phosphorylation barcodes, leading to distinct physiological processes.

Phosphorylation barcodes direct biased chemokine signaling at CXCR3.

G protein-coupled receptor (GPCR) biased agonism, the activation of some signaling pathways over others, is thought to largely be due to differential receptor phosphorylation, or "phosphorylation barcodes." At chemokine receptors, ligands act as "biased agonists" with complex signaling profiles, which contributes to the limited success in pharmacologically targeting these receptors. Here, mass spectrometry-based global phosphoproteomics revealed that CXCR3 chemokines generate different phosphorylation barcodes associated with differential transducer activation. Chemokine stimulation resulted in distinct changes throughout the kinome in global phosphoproteomic studies. Mutation of CXCR3 phosphosites altered ß-arrestin conformation in cellular assays and was confirmed by molecular dynamics simulations. T cells expressing phosphorylation-deficient CXCR3 mutants resulted in agonist- and receptor-specific chemotactic profiles. Our results demonstrate that CXCR3 chemokines are non-redundant and act as biased agonists through differential encoding of phosphorylation barcodes and lead to distinct physiological processes.

Allosteric modulation of GPCR-induced ß-arrestin trafficking and signaling by a synthetic intrabody.

Agonist-induced phosphorylation of G protein-coupled receptors (GPCRs) is a primary determinant of ß-arrestin (ßarr) recruitment and trafficking. For several GPCRs such as the vasopressin receptor subtype 2 (V2R), agonist-stimulation first drives the translocation of ßarrs to the plasma membrane, followed by endosomal trafficking, which is generally considered to be orchestrated by multiple phosphorylation sites. We have previously shown that mutation of a single phosphorylation site in the V2R (i.e., V2RT360A) results in near-complete loss of ßarr translocation to endosomes despite robust recruitment to the plasma membrane, and compromised ERK1/2 activation. Here, we discover that a synthetic intrabody (Ib30), which selectively recognizes activated ßarr1, efficiently rescues the endosomal trafficking of ßarr1 and ERK1/2 activation for V2RT360A. Molecular dynamics simulations reveal that Ib30 enriches active-like ßarr1 conformation with respect to the inter-domain rotation, and cellular assays demonstrate that it also enhances ßarr1-ß2-adaptin interaction. Our data provide an experimental framework to positively modulate the receptor-transducer-effector axis for GPCRs using intrabodies, which can be potentially integrated in the paradigm of GPCR-targeted drug discovery.

Entrectinib-A SARS-CoV-2 Inhibitor in Human Lung Tissue (HLT) Cells.

Since the start of the COVID-19 outbreak, pharmaceutical companies and research groups have focused on the development of vaccines and antiviral drugs against SARS-CoV-2. Here, we apply a drug repurposing strategy to identify drug candidates that are able to block the entrance of the virus into human cells. By combining virtual screening with in vitro pseudovirus assays and antiviral assays in Human Lung Tissue (HLT) cells, we identify entrectinib as a potential antiviral drug.

Mechanistic insights into dopaminergic and serotonergic neurotransmission - concerted interactions with helices 5 and 6 drive the functional outcome.

Brain functions rely on neurotransmitters that mediate communication between billions of neurons. Disruption of this communication can result in a plethora of psychiatric and neurological disorders. In this work, we combine molecular dynamics simulations, live-cell biosensor and electrophysiological assays to investigate the action of the neurotransmitter dopamine at the dopaminergic D2 receptor (D2R). The study of dopamine and closely related chemical probes reveals how neurotransmitter binding translates into the activation of distinct subsets of D2R effectors (i.e.: Gi2, GoB, Gz and ß-arrestin 2). Ligand interactions with key residues in TM5 (S5.42) and TM6 (H6.55) in the D2R binding pocket yield a dopamine-like coupling signature, whereas exclusive TM5 interaction is typically linked to preferential G protein coupling (in particular GoB) over ß-arrestin. Further experiments for serotonin receptors indicate that the reported molecular mechanism is shared by other monoaminergic neurotransmitter receptors. Ultimately, our study highlights how sequence variation in position 6.55 is used by nature to fine-tune ß-arrestin recruitment and in turn receptor signaling and internalization of neurotransmitter receptors.

A molecular sensor for cholesterol in the human serotonin1A receptor.

The function of several G protein-coupled receptors (GPCRs) exhibits cholesterol sensitivity. Cholesterol sensitivity of GPCRs could be attributed to specific sequence and structural features, such as the cholesterol recognition/interaction amino acid consensus (CRAC) motif, that facilitate their cholesterol-receptor interaction. In this work, we explored the molecular basis of cholesterol sensitivity exhibited by the serotonin1A receptor, the most studied GPCR in the context of cholesterol sensitivity, by generating mutants of key residues in CRAC motifs in transmembrane helix 2 (TM2) and TM5 of the receptor. Our results show that a lysine residue (K101) in one of the CRAC motifs is crucial for sensing altered membrane cholesterol levels. Insights from all-atom molecular dynamics simulations showed that cholesterol-sensitive functional states of the serotonin1A receptor are associated with reduced conformational dynamics of extracellular loops of the receptor. These results constitute one of the first reports on the molecular mechanism underlying cholesterol sensitivity of GPCRs.

Dopamine D2 Receptor Agonist Binding Kinetics-Role of a Conserved Serine Residue.

The forward (kon) and reverse (koff) rate constants of drug-target interactions have important implications for therapeutic efficacy. Hence, time-resolved assays capable of measuring these binding rate constants may be informative to drug discovery efforts. Here, we used an ion channel activation assay to estimate the kons and koffs of four dopamine D2 receptor (D2R) agonists; dopamine (DA), p-tyramine, (R)- and (S)-5-OH-dipropylaminotetralin (DPAT). We further probed the role of the conserved serine S1935.42 by mutagenesis, taking advantage of the preferential interaction of (S)-, but not (R)-5-OH-DPAT with this residue. Results suggested similar koffs for the two 5-OH-DPAT enantiomers at wild-type (WT) D2R, both being slower than the koffs of DA and p-tyramine. Conversely, the kon of (S)-5-OH-DPAT was estimated to be higher than that of (R)-5-OH-DPAT, in agreement with the higher potency of the (S)-enantiomer. Furthermore, S1935.42A mutation lowered the kon of (S)-5-OH-DPAT and reduced the potency difference between the two 5-OH-DPAT enantiomers. Kinetic Kds derived from the koff and kon estimates correlated well with EC50 values for all four compounds across four orders of magnitude, strengthening the notion that our assay captured meaningful information about binding kinetics. The approach presented here may thus prove valuable for characterizing D2R agonist candidate drugs.

Ligand with Two Modes of Interaction with the Dopamine D2 Receptor-An Induced-Fit Mechanism of Insurmountable Antagonism.

A solid understanding of the mechanisms governing ligand binding is crucial for rational design of therapeutics targeting the dopamine D2 receptor (D2R). Here, we use G protein-coupled inward rectifier potassium (GIRK) channel activation in Xenopus oocytes to measure the kinetics of D2R antagonism by a series of aripiprazole analogues, as well as the recovery of dopamine (DA) responsivity upon washout. The aripiprazole analogues comprise an orthosteric and a secondary pharmacophore and differ by the length of the saturated carbon linker joining these two pharmacophores. Two compounds containing 3- and 5-carbon linkers allowed for a similar extent of recovery from antagonism in the presence of 1 or 100 µM DA (>25 and >90% of control, respectively), whereas recovery was less prominent (∼20%) upon washout of the 4-carbon linker compound, SV-III-130, both with 1 and 100 µM DA. Prolonging the coincubation time with SV-III-130 further diminished recovery. Curve-shift experiments were consistent with competition between SV-III-130 and DA. Two mutations in the secondary binding pocket (V91A and E95A) of D2R decreased antagonistic potency and increased recovery from SV-III-130 antagonism, whereas a third mutation (L94A) only increased recovery. Our results suggest that the secondary binding pocket influences recovery from inhibition by the studied aripiprazole analogues. We propose a mechanism, supported by in silico modeling, whereby SV-III-130 initially binds reversibly to the D2R, after which the drug-receptor complex undergoes a slow transition to a second ligand-bound state, which is dependent on secondary binding pocket integrity and irreversible during the time frame of our experiments.

Distinct phosphorylation sites in a prototypical GPCR differently orchestrate ß-arrestin interaction, trafficking, and signaling.

Agonist-induced phosphorylation of G protein-coupled receptors (GPCRs) is a key determinant for their interaction with ß-arrestins (ßarrs) and subsequent functional responses. Therefore, it is important to decipher the contribution and interplay of different receptor phosphorylation sites in governing ßarr interaction and functional outcomes. Here, we find that several phosphorylation sites in the human vasopressin receptor (V2R), positioned either individually or in clusters, differentially contribute to ßarr recruitment, trafficking, and ERK1/2 activation. Even a single phosphorylation site in V2R, suitably positioned to cross-talk with a key residue in ßarrs, has a decisive contribution in ßarr recruitment, and its mutation results in strong G-protein bias. Molecular dynamics simulation provides mechanistic insights into the pivotal role of this key phosphorylation site in governing the stability of ßarr interaction and regulating the interdomain rotation in ßarrs. Our findings uncover important structural aspects to better understand the framework of GPCR-ßarr interaction and biased signaling.

How Do Molecular Dynamics Data Complement Static Structural Data of GPCRs.

G protein-coupled receptors (GPCRs) are implicated in nearly every physiological process in the human body and therefore represent an important drug targeting class. Advances in X-ray crystallography and cryo-electron microscopy (cryo-EM) have provided multiple static structures of GPCRs in complex with various signaling partners. However, GPCR functionality is largely determined by their flexibility and ability to transition between distinct structural conformations. Due to this dynamic nature, a static snapshot does not fully explain the complexity of GPCR signal transduction. Molecular dynamics (MD) simulations offer the opportunity to simulate the structural motions of biological processes at atomic resolution. Thus, this technique can incorporate the missing information on protein flexibility into experimentally solved structures. Here, we review the contribution of MD simulations to complement static structural data and to improve our understanding of GPCR physiology and pharmacology, as well as the challenges that still need to be overcome to reach the full potential of this technique.

GPCRmd uncovers the dynamics of the 3D-GPCRome.

G-protein-coupled receptors (GPCRs) are involved in numerous physiological processes and are the most frequent targets of approved drugs. The explosion in the number of new three-dimensional (3D) molecular structures of GPCRs (3D-GPCRome) over the last decade has greatly advanced the mechanistic understanding and drug design opportunities for this protein family. Molecular dynamics (MD) simulations have become a widely established technique for exploring the conformational landscape of proteins at an atomic level. However, the analysis and visualization of MD simulations require efficient storage resources and specialized software. Here we present GPCRmd (http://gpcrmd.org/), an online platform that incorporates web-based visualization capabilities as well as a comprehensive and user-friendly analysis toolbox that allows scientists from different disciplines to visualize, analyze and share GPCR MD data. GPCRmd originates from a community-driven effort to create an open, interactive and standardized database of GPCR MD simulations.

Key phosphorylation sites in GPCRs orchestrate the contribution of ß-Arrestin 1 in ERK1/2 activation.

ß-arrestins (ßarrs) are key regulators of G protein-coupled receptor (GPCR) signaling and trafficking, and their knockdown typically leads to a decrease in agonist-induced ERK1/2 MAP kinase activation. Interestingly, for some GPCRs, knockdown of ßarr1 augments agonist-induced ERK1/2 phosphorylation although a mechanistic basis for this intriguing phenomenon is unclear. Here, we use selected GPCRs to explore a possible correlation between the spatial positioning of receptor phosphorylation sites and the contribution of ßarr1 in ERK1/2 activation. We discover that engineering a spatially positioned double-phosphorylation-site cluster in the bradykinin receptor (B2 R), analogous to that present in the vasopressin receptor (V2 R), reverses the contribution of ßarr1 in ERK1/2 activation from inhibitory to promotive. An intrabody sensor suggests a conformational mechanism for this role reversal of ßarr1, and molecular dynamics simulation reveals a bifurcated salt bridge between this double-phosphorylation site cluster and Lys294 in the lariat loop of ßarr1, which directs the orientation of the lariat loop. Our findings provide novel insights into the opposite roles of ßarr1 in ERK1/2 activation for different GPCRs with a direct relevance to biased agonism and novel therapeutics.

Publisher Correction: GPCRmd uncovers the dynamics of the 3D-GPCRome.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

A highly conserved δ-opioid receptor region determines RGS4 interaction.

The δ-opioid receptor (δ-OR) couples to Gi/Go proteins to modulate a variety of responses in the nervous system. Τhe regulator of G protein signalling 4 (RGS4) was previously shown to directly interact within the C-terminal region of δ-OR using its N-terminal domain to negatively modulate opioid receptor signalling. Herein, using molecular dynamics simulations and in vitro pull-down experiments we delimit this interaction to 12 helix 8 residues of δ-ΟR and to the first 17 N-terminal residues (NT) of RGS4. Monitoring the complex arrangement and stabilization between RGS4 and δ-OR by molecular dynamics simulations combined with mutagenesis studies, we defined that two critical interactions are formed: one between Phe329 of helix8 of δ-ΟR and Pro9 of the NT of RGS4 and the other a salt bridge between Glu323 of δ-ΟR and Lys17 of RGS4. Our observations allow drafting for the first time a structural model of a ternary complex including the δ-opioid receptor, a G protein and a RGS protein. Furthermore, the high degree of conservation among opioid receptors of the RGS4-binding region, points to a conserved interaction mode between opioid receptors and this important regulatory protein.

Synthesis, molecular modelling studies and biological evaluation of new oxoeicosanoid receptor 1 agonists.

The oxoeicosanoid receptor 1 (OXER1) is a member of the G-protein coupled receptors (GPCR) family, and is involved in inflammatory processes and oncogenesis. As such it is an attractive target for pharmacological intervention. The present study aimed to shed light on the molecular fundaments of OXER1 modulation using chemical probes structurally related to the natural agonist 5-oxo-ETE. In a first step, 5-oxo-ETE and its closely related derivatives (5-oxo-EPE and 4-oxo-DHA) were obtained by conducting concise and high-yielding syntheses. The biological activity of obtained compounds was assessed in terms of potency (EC50) and efficacy (Emax) for arrestin recruitment. Finally, molecular modelling and simulation were used to explore binding characteristics of 5-oxo-ETE and derivatives with the aim to rationalize biological activity. Our data suggest that the tested 5-oxo-ETE derivatives (i) insert quickly into the membrane, (ii) access the receptor via transmembrane helices (TMs) 5 and 6 from the membrane side and (iii) drive potency and efficacy by differential interaction with TM5 and 7. Most importantly, we found that the methyl ester of 5-oxo-ETE (1a) showed even a higher maximum response than the natural agonist (1). In contrast, shifting the 5-oxo group into position 4 results in inactive compounds (4-oxo DHA compounds (3) and (3a)). All in all, our study provides relevant structural data that help understanding better OXER1 functionality and its modulation. The structural information presented herein will be useful for designing new lead compounds with desired signalling profiles.

Protein-Protein Docking in Drug Design and Discovery.

Protein-protein interactions (PPIs) are responsible for a number of key physiological processes in the living cells and underlie the pathomechanism of many diseases. Nowadays, along with the concept of so-called "hot spots" in protein-protein interactions, which are well-defined interface regions responsible for most of the binding energy, these interfaces can be targeted with modulators. In order to apply structure-based design techniques to design PPIs modulators, a three-dimensional structure of protein complex has to be available. In this context in silico approaches, in particular protein-protein docking, are a valuable complement to experimental methods for elucidating 3D structure of protein complexes. Protein-protein docking is easy to use and does not require significant computer resources and time (in contrast to molecular dynamics) and it results in 3D structure of a protein complex (in contrast to sequence-based methods of predicting binding interfaces). However, protein-protein docking cannot address all the aspects of protein dynamics, in particular the global conformational changes during protein complex formation. In spite of this fact, protein-protein docking is widely used to model complexes of water-soluble proteins and less commonly to predict structures of transmembrane protein assemblies, including dimers and oligomers of G protein-coupled receptors (GPCRs). In this chapter we review the principles of protein-protein docking, available algorithms and software and discuss the recent examples, benefits, and drawbacks of protein-protein docking application to water-soluble proteins, membrane anchoring and transmembrane proteins, including GPCRs.

Challenges and Opportunities in Drug Discovery of Biased Ligands.

The observation of biased agonism in G protein-coupled receptors (GPCRs) has provided new approaches for the development of more efficacious and safer drugs. However, in order to rationally design biased drugs, one must understand the molecular basis of this phenomenon. Computational approaches can help in exploring the conformational universe of GPCRs and detecting conformational states with relevance for distinct functional outcomes. This information is extremely valuable for the development of new therapeutic agents that promote desired conformational receptor states and responses while avoiding the ones leading to undesired side-effects.This book chapter intends to introduce the reader to powerful computational approaches for sampling the conformational space of these receptors, focusing first on molecular dynamics and the analysis of the produced data through methods such as dimensionality reduction, Markov State Models and adaptive sampling. Then, we show how to seek for compounds that target distinct conformational states via docking and virtual screening. In addition, we describe how to detect receptor-ligand interactions that drive signaling bias and comment current challenges and opportunities of presented methods.