AlloViz: A tool for the calculation and visualisation of protein allosteric communication networks.

Allostery, the presence of functional interactions between distant parts of proteins, is a critical concept in the field of biochemistry and molecular biology, particularly in the context of protein function and regulation. Understanding the principles of allosteric regulation is essential for advancing our knowledge of biology and developing new therapeutic strategies. This paper presents AlloViz, an open-source Python package designed to quantitatively determine, analyse, and visually represent allosteric communication networks on the basis of molecular dynamics (MD) simulation data. The software integrates well-known techniques for understanding allosteric properties simplifying the process of accessing, rationalising, and representing protein allostery and communication routes. It overcomes the inefficiency of having multiple methods with heterogeneous implementations and showcases the advantages of using MD simulations and multiple replicas to obtain statistically sound information on protein dynamics; it also enables the calculation of "consensus-like" scores aggregating methods that consider multiple structural aspects of allosteric networks. We demonstrate the features of AlloViz on two proteins: ß-arrestin 1, a key player for regulating G protein-coupled receptor (GPCR) signalling, and the protein tyrosine phosphatase 1B, an important pharmaceutical target for allosteric inhibitors. The software includes comprehensive documentation and examples, tutorials, and a user-friendly graphical interface.

Pharmacological fingerprint of antipsychotic drugs at the serotonin 5-HT2A receptor.

The intricate involvement of the serotonin 5-HT2A receptor (5-HT2AR) both in schizophrenia and in the activity of antipsychotic drugs is widely acknowledged. The currently marketed antipsychotic drugs, although effective in managing the symptoms of schizophrenia to a certain extent, are not without their repertoire of serious side effects. There is a need for better therapeutics to treat schizophrenia for which understanding the mechanism of action of the current antipsychotic drugs is imperative. With bioluminescence resonance energy transfer (BRET) assays, we trace the signaling signature of six antipsychotic drugs belonging to three generations at the 5-HT2AR for the entire spectrum of signaling pathways activated by serotonin (5-HT). The antipsychotic drugs display previously unidentified pathway preference at the level of the individual Gα subunits and ß-arrestins. In particular, risperidone, clozapine, olanzapine and haloperidol showed G protein-selective inverse agonist activity. In addition, G protein-selective partial agonism was found for aripiprazole and cariprazine. Pathway-specific apparent dissociation constants determined from functional analyses revealed distinct coupling-modulating capacities of the tested antipsychotics at the different 5-HT-activated pathways. Computational analyses of the pharmacological and structural fingerprints support a mechanistically based clustering that recapitulate the clinical classification (typical/first generation, atypical/second generation, third generation) of the antipsychotic drugs. The study provides a new framework to functionally classify antipsychotics that should represent a useful tool for the identification of better and safer neuropsychiatric drugs and allows formulating hypotheses on the links between specific signaling cascades and in the clinical outcomes of the existing drugs.

Molecular insights into atypical modes of ß-arrestin interaction with seven transmembrane receptors.

ß-arrestins (ßarrs) are multifunctional proteins involved in signaling and regulation of seven transmembrane receptors (7TMRs), and their interaction is driven primarily by agonist-induced receptor activation and phosphorylation. Here, we present seven cryo-electron microscopy structures of ßarrs either in the basal state, activated by the muscarinic receptor subtype 2 (M2R) through its third intracellular loop, or activated by the ßarr-biased decoy D6 receptor (D6R). Combined with biochemical, cellular, and biophysical experiments, these structural snapshots allow the visualization of atypical engagement of ßarrs with 7TMRs and also reveal a structural transition in the carboxyl terminus of ßarr2 from a ß strand to an α helix upon activation by D6R. Our study provides previously unanticipated molecular insights into the structural and functional diversity encoded in 7TMR-ßarr complexes with direct implications for exploring novel therapeutic avenues.

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.

On the construction of LIECE models for the serotonin receptor 5-HT[Formula: see text]R.

Computer-aided approaches to ligand design need to balance accuracy with speed. This is particularly true for one of the key parameters to be optimized during ligand development, the free energy of binding ([Formula: see text]G[Formula: see text]). Here, we developed simple models based on the Linear Interaction Energy approximation to free energy calculation for a G protein-coupled receptor, the serotonin receptor 2A, and critically evaluated their accuracy. Several lessons can be taken from our calculations, providing information on the influence of the docking software used, the conformational state of the receptor, the cocrystallized ligand, and its comparability to the training/test ligands.

Plasma membrane preassociation drives ß-arrestin coupling to receptors and activation.

ß-arrestin plays a key role in G protein-coupled receptor (GPCR) signaling and desensitization. Despite recent structural advances, the mechanisms that govern receptor-ß-arrestin interactions at the plasma membrane of living cells remain elusive. Here, we combine single-molecule microscopy with molecular dynamics simulations to dissect the complex sequence of events involved in ß-arrestin interactions with both receptors and the lipid bilayer. Unexpectedly, our results reveal that ß-arrestin spontaneously inserts into the lipid bilayer and transiently interacts with receptors via lateral diffusion on the plasma membrane. Moreover, they indicate that, following receptor interaction, the plasma membrane stabilizes ß-arrestin in a longer-lived, membrane-bound state, allowing it to diffuse to clathrin-coated pits separately from the activating receptor. These results expand our current understanding of ß-arrestin function at the plasma membrane, revealing a critical role for ß-arrestin preassociation with the lipid bilayer in facilitating its interactions with receptors and subsequent activation.

Migration mediated by the oxysterol receptor GPR183 depends on arrestin coupling but not receptor internalization.

The chemotactic G protein-coupled receptor GPR183 and its most potent endogenous oxysterol ligand 7α,25-dihydroxycholesterol (7α,25-OHC) are important for immune cell positioning in secondary lymphoid tissues. This receptor-ligand pair is associated with various diseases, in some cases contributing favorably and in other cases adversely, making GPR183 an attractive target for therapeutic intervention. We investigated the mechanisms underlying GPR183 internalization and the role of internalization in the main biological function of the receptor, chemotaxis. We found that the C terminus of the receptor was important for ligand-induced internalization but less so for constitutive (ligand-independent) internalization. ß-arrestin potentiated ligand-induced internalization but was not required for ligand-induced or constitutive internalization. Caveolin and dynamin were the main mediators of both constitutive and ligand-induced receptor internalization in a mechanism independent of G protein activation. Clathrin-mediated endocytosis also contributed to constitutive GPR183 internalization in a ß-arrestin-independent manner, suggesting the existence of different pools of surface-localized GPR183. Chemotaxis mediated by GPR183 depended on receptor desensitization by ß-arrestins but could be uncoupled from internalization, highlighting an important biological role for the recruitment of ß-arrestin to GPR183. The role of distinct pathways in internalization and chemotaxis may aid in the development of GPR183-targeting drugs for specific disease contexts.

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.

Impact of membrane lipid polyunsaturation on dopamine D2 receptor ligand binding and signaling.

Increasing evidence supports a relationship between lipid metabolism and mental health. In particular, the biostatus of polyunsaturated fatty acids (PUFAs) correlates with some symptoms of psychiatric disorders, as well as the efficacy of pharmacological treatments. Recent findings highlight a direct association between brain PUFA levels and dopamine transmission, a major neuromodulatory system implicated in the etiology of psychiatric symptoms. However, the mechanisms underlying this relationship are still unknown. Here we demonstrate that membrane enrichment in the n-3 PUFA docosahexaenoic acid (DHA), potentiates ligand binding to the dopamine D2 receptor (D2R), suggesting that DHA acts as an allosteric modulator of this receptor. Molecular dynamics simulations confirm that DHA has a high preference for interaction with the D2R and show that membrane unsaturation selectively enhances the conformational dynamics of the receptor around its second intracellular loop. We find that membrane unsaturation spares G protein activity but potentiates the recruitment of ß-arrestin in cells. Furthermore, in vivo n-3 PUFA deficiency blunts the behavioral effects of two D2R ligands, quinpirole and aripiprazole. These results highlight the importance of membrane unsaturation for D2R activity and provide a putative mechanism for the ability of PUFAs to enhance antipsychotic efficacy.

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.

In Silico Study of Allosteric Communication Networks in GPCR Signaling Bias.

Signaling bias is a promising characteristic of G protein-coupled receptors (GPCRs) as it provides the opportunity to develop more efficacious and safer drugs. This is because biased ligands can avoid the activation of pathways linked to side effects whilst still producing the desired therapeutic effect. In this respect, a deeper understanding of receptor dynamics and implicated allosteric communication networks in signaling bias can accelerate the research on novel biased drug candidates. In this review, we aim to provide an overview of computational methods and techniques for studying allosteric communication and signaling bias in GPCRs. This includes (i) the detection of allosteric communication networks and (ii) the application of network theory for extracting relevant information pipelines and highly communicated sites in GPCRs. We focus on the most recent research and highlight structural insights obtained based on the framework of allosteric communication networks and network theory for GPCR signaling bias.

Community guidelines for GPCR ligand bias: IUPHAR review 32.

GPCRs modulate a plethora of physiological processes and mediate the effects of one-third of FDA-approved drugs. Depending on which ligand activates a receptor, it can engage different intracellular transducers. This 'biased signalling' paradigm requires that we now characterize physiological signalling not just by receptors but by ligand-receptor pairs. Ligands eliciting biased signalling may constitute better drugs with higher efficacy and fewer adverse effects. However, ligand bias is very complex, making reproducibility and description challenging. Here, we provide guidelines and terminology for any scientists to design and report ligand bias experiments. The guidelines will aid consistency and clarity, as the basic receptor research and drug discovery communities continue to advance our understanding and exploitation of ligand bias. Scientific insight, biosensors, and analytical methods are still evolving and should benefit from and contribute to the implementation of the guidelines, together improving translation from in vitro to disease-relevant in vivo models.

SCoV2-MD: a database for the dynamics of the SARS-CoV-2 proteome and variant impact predictions.

SCoV2-MD (www.scov2-md.org) is a new online resource that systematically organizes atomistic simulations of the SARS-CoV-2 proteome. The database includes simulations produced by leading groups using molecular dynamics (MD) methods to investigate the structure-dynamics-function relationships of viral proteins. SCoV2-MD cross-references the molecular data with the pandemic evolution by tracking all available variants sequenced during the pandemic and deposited in the GISAID resource. SCoV2-MD enables the interactive analysis of the deposited trajectories through a web interface, which enables users to search by viral protein, isolate, phylogenetic attributes, or specific point mutation. Each mutation can then be analyzed interactively combining static (e.g. a variety of amino acid substitution penalties) and dynamic (time-dependent data derived from the dynamics of the local geometry) scores. Dynamic scores can be computed on the basis of nine non-covalent interaction types, including steric properties, solvent accessibility, hydrogen bonding, and other types of chemical interactions. Where available, experimental data such as antibody escape and change in binding affinities from deep mutational scanning experiments are also made available. All metrics can be combined to build predefined or custom scores to interrogate the impact of evolving variants on protein structure and function.

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.

Discovery of a dual Ras and ARF6 inhibitor from a GPCR endocytosis screen.

Internalization and intracellular trafficking of G protein-coupled receptors (GPCRs) play pivotal roles in cell responsiveness. Dysregulation in receptor trafficking can lead to aberrant signaling and cell behavior. Here, using an endosomal BRET-based assay in a high-throughput screen with the prototypical GPCR angiotensin II type 1 receptor (AT1R), we sought to identify receptor trafficking inhibitors from a library of ~115,000 small molecules. We identified a novel dual Ras and ARF6 inhibitor, which we named Rasarfin, that blocks agonist-mediated internalization of AT1R and other GPCRs. Rasarfin also potently inhibits agonist-induced ERK1/2 signaling by GPCRs, and MAPK and Akt signaling by EGFR, as well as prevents cancer cell proliferation. In silico modeling and in vitro studies reveal a unique binding modality of Rasarfin within the SOS-binding domain of Ras. Our findings unveil a class of dual small G protein inhibitors for receptor trafficking and signaling, useful for the inhibition of oncogenic cellular responses.

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.

Structural dynamics bridge the gap between the genetic and functional levels of GPCRs.

G protein-coupled receptors (GPCRs) are implicated in nearly all physiological processes in the human body and represent an important drug targeting class. The genes encoding the different GPCR (sub)types determine their specific functionality, which can be altered by natural genetic variants and isoforms. Deciphering the molecular link between sequence diversity and its functional consequences is a current challenge and critical for the comprehension of the physiological response of GPCRs. It requires a global understanding of how protein sequence translates into protein structure, how this impacts the structural motions of the protein, and, finally, how all these factors determine the receptor functionality. Here, we discuss available resources and state-of-the-art computational approaches to address this question.