An alpha-carbon template for the transmembrane helices in the rhodopsin family of G-protein-coupled receptors

PDB-BRE: A ligand-protein interaction binding residue extractor based on Protein Data Bank

Proteins typically exert their biological functions by interacting with other biomolecules or ligands. The study of ligand-protein interactions is crucial in elucidating the biological mechanisms of proteins. Most existing studies have focused on analyzing ligand-protein interactions, and they ignore the additional situational of inserted and modified residues. Besides, the resources often support only a single ligand type and cannot obtain satisfied results in analyzing novel complexes. Therefore, it is important to develop a general analytical tool to extract the binding residues of ligand-protein interactions in complexes fully. In this study, we propose a ligand-protein interaction binding residue extractor (PDB-BRE), which can be used to automatically extract interacting ligand or protein-binding residues from complex three-dimensional (3D) structures based on the RCSB Protein Data Bank (RCSB PDB). PDB-BRE offers a notable advantage in its comprehensive support for analyzing six distinct types of ligands, including proteins, peptides, DNA, RNA, mixed DNA and RNA entities, and non-polymeric entities. Moreover, it takes into account the consideration of inserted and modified residues within complexes. Compared to other state-of-the-art methods, PDB-BRE is more suitable for massively parallel batch analysis, and can be directly applied for downstream tasks, such as predicting binding residues of novel complexes. PDB-BRE is freely available at http://bliulab.net/PDB-BRE.

Structural aspects of rod opsin and their implication in genetic diseases

Vision in dim-light conditions is triggered by photoactivation of rhodopsin, the visual pigment of rod photoreceptor cells. Rhodopsin is made of a protein, the G protein coupled receptor (GPCR) opsin, and the chromophore 11-cis-retinal. Vertebrate rod opsin is the GPCR best characterized at the atomic level of detail. Since the release of the first crystal structure 20 years ago, a huge number of structures have been released that, in combination with valuable spectroscopic determinations, unveiled most aspects of the photobleaching process. A number of spontaneous mutations of rod opsin have been found linked to vision-impairing diseases like autosomal dominant or autosomal recessive retinitis pigmentosa (adRP or arRP, respectively) and autosomal congenital stationary night blindness (adCSNB). While adCSNB is mainly caused by constitutive activation of rod opsin, RP shows more variegate determinants affecting different aspects of rod opsin function. The vast majority of missense rod opsin mutations affects folding and trafficking and is linked to adRP, an incurable disease that awaits light on its molecular structure determinants. This review article summarizes all major structural information available on vertebrate rod opsin conformational states and the insights gained so far into the structural determinants of adCSNB and adRP linked to rod opsin mutations. Strategies to design small chaperones with therapeutic potential for selected adRP rod opsin mutants will be discussed as well.

Revealing the Positive Binding Cooperativity Mechanism between the Orthosteric and the Allosteric Antagonists of CCR2 by Metadynamics and Gaussian Accelerated Molecular Dynamics Simulations

CC chemokine receptor 2 (CCR2) and its endogenous CC chemokine ligands are associated with numerous inflammatory, neurodegenerative diseases, and cancer. CCR2 is becoming an attractive target in the treatment of autoimmune disease and neurodegenerative diseases. The orthosteric antagonist BMS-681 and allosteric antagonist CCR2-RA-[R] of CCR2 show positive binding cooperativity. We performed well-tempered metadynamics simulations and Gaussian accelerated MD simulations to reveal the influence of the orthosteric antagonist on the unbinding of allosteric antagonist of CCR2. We revealed different unbinding pathways of CCR2-RA-[R] in binary complex CCR2-VT5 and ternary complex CCR2-73R-VT5. The different unbinding pathways of CCR2-RA-[R] are due to the conformational dynamics of TM6. We obtained the significant conformational differences of the intracellular side of TM6 upon CCR2 binding to different ligands by GaMD simulation. The conformational dynamics of TM6 are consistent with the unbinding pathway analysis. GaMD simulations indicate that BMS-681 binding restricts the bend of intracellular side of TM6 by stabilizing the extracellular sides of TM6 and TM7. The charged residues Arg206^5.43 of TM5 and Glu291^7.39 of TM7 play key roles in stabling TM7 and TM6. TM6 and TM7 are crucial components in the orthosteric and allosteric binding sites. Our results illustrate the conformational details about the effect of the orthosteric antagonist on the allosteric antagonist of CCR2. The conformational dynamics of CCR2 upon binding to different ligands can provide a rational basis for development of allosteric ligands of CCR2.

Extracellular loops 2 and 3 of the calcitonin receptor selectively modify agonist binding and efficacy

Class B peptide hormone GPCRs are targets for the treatment of major chronic disease. Peptide ligands of these receptors display biased agonism and this may provide future therapeutic advantage. Recent active structures of the calcitonin (CT) and glucagon-like peptide-1 (GLP-1) receptors reveal distinct engagement of peptides with extracellular loops (ECLs) 2 and 3, and mutagenesis of the GLP-1R has implicated these loops in dynamics of receptor activation. In the current study, we have mutated ECLs 2 and 3 of the human CT receptor (CTR), to interrogate receptor expression, peptide affinity and efficacy. Integration of these data with insights from the CTR and GLP-1R active structures, revealed marked diversity in mechanisms of peptide engagement and receptor activation between the CTR and GLP-1R. While the CTR ECL2 played a key role in conformational propagation linked to Gs/cAMP signalling this was mechanistically distinct from that of GLP-1R ECL2. Moreover, ECL3 was a hotspot for distinct ligand- and pathway-specific effects, and this has implications for the future design of biased agonists of class B GPCRs.

A Structural Framework for GPCR Chemogenomics: What's In a Residue Number?

The recent surge of crystal structures of G protein-coupled receptors (GPCRs), as well as comprehensive collections of sequence, structural, ligand bioactivity, and mutation data, has enabled the development of integrated chemogenomics workflows for this important target family. This chapter will focus on cross-family and cross-class studies of GPCRs that have pinpointed the need for, and the implementation of, a generic numbering scheme for referring to specific structural elements of GPCRs. Sequence- and structure-based numbering schemes for different receptor classes will be introduced and the remaining caveats will be discussed. The use of these numbering schemes has facilitated many chemogenomics studies such as consensus binding site definition, binding site comparison, ligand repurposing (e.g. for orphan receptors), sequence-based pharmacophore generation for homology modeling or virtual screening, and class-wide chemogenomics studies of GPCRs.

Modeling and Deorphanization of Orphan GPCRs

Despite tremendous efforts, approximately 120 GPCRs remain orphan. Their physiological functions and their potential roles in diseases are poorly understood. Orphan GPCRs are extremely important because they may provide novel therapeutic targets for unmet medical needs. As a complement to experimental approaches, molecular modeling and virtual screening are efficient techniques to discover synthetic surrogate ligands which can help to elucidate the role of oGPCRs. Constitutively activated mutants and recently published active structures of GPCRs provide stimulating opportunities for building active molecular models for oGPCRs and identifying activators using virtual screening of compound libraries. We describe the molecular modeling and virtual screening process we have applied in the discovery of surrogate ligands, and provide examples for CCKA, a simulated oGPCR, and for two oGPCRs, GPR52 and GPR34.

Apelinergic System Structure and Function

Apelin and apela (ELABELA/ELA/Toddler) are two peptide ligands for a class A G-protein-coupled receptor named the apelin receptor (AR/APJ/APLNR). Ligand-AR interactions have been implicated in regulation of the adipoinsular axis, cardiovascular system, and central nervous system alongside pathological processes. Each ligand may be processed into a variety of bioactive isoforms endogenously, with apelin ranging from 13 to 55 amino acids and apela from 11 to 32, typically being cleaved C-terminal to dibasic proprotein convertase cleavage sites. The C-terminal region of the respective precursor protein is retained and is responsible for receptor binding and subsequent activation. Interestingly, both apelin and apela exhibit isoform-dependent variability in potency and efficacy under various physiological and pathological conditions, but most studies focus on a single isoform. Biophysical behavior and structural properties of apelin and apela isoforms show strong correlations with functional studies, with key motifs now well determined for apelin. Unlike its ligands, the AR has been relatively difficult to characterize by biophysical techniques, with most characterization to date being focused on effects of mutagenesis. This situation may improve following a recently reported AR crystal structure, but there are still barriers to overcome in terms of comprehensive biophysical study. In this review, we summarize the three components of the apelinergic system in terms of structure-function correlation, with a particular focus on isoform-dependent properties, underlining the potential for regulation of the system through multiple endogenous ligands and isoforms, isoform-dependent pharmacological properties, and biological membrane-mediated receptor interaction. © 2018 American Physiological Society. Compr Physiol 8:407-450, 2018.

Unconventional Roles of Opsins

Rhodopsin is the classical light sensor. Although rhodopsin has long been known to be important for image formation in the eye, the requirements for opsins in non-image formation and in extraocular light sensation were revealed much later. Most recent is the demonstration that an opsin in the fruit fly, Drosophila melanogaster, is expressed in pacemaker neurons in the brain and functions in light entrainment of circadian rhythms. However, the biggest surprise is that opsins have light-independent roles, countering more than a century of dogma that they function exclusively as light sensors. Through studies in Drosophila, light-independent roles of opsins have emerged in temperature sensation and hearing. Although these findings have been uncovered in the fruit fly, there are hints that opsins have light-independent roles in a wide array of animals, including mammals. Thus, despite the decades of focus on opsins as light detectors, they represent an important new class of polymodal sensory receptor.

Retention of duplicated long-wavelength opsins in mosquito lineages by positive selection and differential expression

Background

Opsins are light sensitive receptors associated with visual processes. Insects typically possess opsins that are stimulated by ultraviolet, short and long wavelength (LW) radiation. Six putative LW-sensitive opsins predicted in the yellow fever mosquito, Aedes aegypti and malaria mosquito, Anopheles gambiae, and eight in the southern house mosquito, Culex quinquefasciatus, suggest gene expansion in the Family Culicidae (mosquitoes) relative to other insects. Here we report the first detailed molecular and evolutionary analyses of LW opsins in three mosquito vectors, with a goal to understanding the molecular basis of opsin-mediated visual processes that could be exploited for mosquito control.

Results

Time of divergence estimates suggest that the mosquito LW opsins originated from 18 or 19 duplication events between 166.9/197.5 to 1.07/0.94 million years ago (MY) and that these likely occurred following the predicted divergence of the lineages Anophelinae and Culicinae 145–226 MY. Fitmodel analyses identified nine amino acid residues in the LW opsins that may be under positive selection. Of these, eight amino acids occur in the N and C termini and are shared among all three species, and one residue in TMIII was unique to culicine species. Alignment of 5′ non-coding regions revealed potential Conserved Non-coding Sequences (CNS) and transcription factor binding sites (TFBS) in seven pairs of LW opsin paralogs.

Conclusions

Our analyses suggest opsin gene duplication and residues possibly associated with spectral tuning of LW-sensitive photoreceptors. We explore two mechanisms - positive selection and differential expression mediated by regulatory units in CNS – that may have contributed to the retention of LW opsin genes in Culicinae and Anophelinae. We discuss the evolution of mosquito LW opsins in the context of major Earth events and possible adaptation of mosquitoes to LW-dominated photo environments, and implications for mosquito control strategies based on disrupting vision-mediated behaviors.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0910-6) contains supplementary material, which is available to authorized users.

Differential CCR7 Targeting in Dendritic Cells by Three Naturally Occurring CC-Chemokines

The CCR7 ligands CCL19 and CCL21 are increasingly recognized as functionally different (biased). Using mature human dendritic cells (DCs), we show that CCL19 is more potent than CCL21 in inducing 3D chemotaxis. Intriguingly, CCL21 induces prolonged and more efficient ERK1/2 activation compared with CCL19 and a C-terminal truncated (tailless) CCL21 in DCs. In contrast, tailless-CCL21 displays increased potency in DC chemotaxis compared with native CCL21. Using a CCL21-specific antibody, we show that CCL21, but not tailless-CCL21, accumulates at the cell surface. In addition, removal of sialic acid from the cell surface by neuraminidase treatment impairs ERK1/2 activation by CCL21, but not by CCL19 or tailless-CCL21. Using standard laboratory cell lines, we observe low potency of both CCL21 and tailless-CCL21 in G protein activation and β-arrestin recruitment compared with CCL19, indicating that the tail itself does not improve receptor interaction. Chemokines interact with their receptors in a stepwise manner with ultimate docking of their N-terminus into the main binding pocket. Employing site-directed mutagenesis we identify residues in this pocket of selective CCL21 importance. We also identify a molecular switch in the top of TM7 important for keeping CCR7 in an inactive conformation (Tyr312), as introduction of the chemokine receptor-conserved Glu (or Ala) induces high constitutive activity. Summarized, we show that the interaction of the tail of CCL21 with polysialic acid is needed for strong ERK signaling, whereas it impairs CCL21-mediated chemotaxis and has no impact on receptor docking consistent with the current model of chemokine:receptor interaction. This indicates that future selective pharmacological targeting of CCL19 versus CCL21 should focus on a differential targeting of the main receptor pocket, while selective targeting of tailless-CCL21 versus CCL21 and CCL19 requires targeting of the glycosaminoglycan (GAG) interaction.

Structural Insights into 5-HT1A/D4 Selectivity of WAY-100635 Analogues: Molecular Modeling, Synthesis, and in Vitro Binding

The resurgence of interest in 5-HT1A receptors as a therapeutic target requires the existence of highly selective 5-HT1A ligands. To date, WAY-100635 has been the prototypical antagonist of these receptors. However, this compound also has significant affinity for and activity at D4 dopamine receptors. In this context, this work was aimed at better understanding the 5-HT1A/D4 selectivity of WAY-100635 and analogues from a structural point of view. In silico investigations revealed two key interactions for the 5-HT1A/D4 selectivity of WAY-100635 and analogues. First, a hydrogen bond only found with the Ser 7.36 of D4 receptor appeared to be the key for a higher D4 affinity for newly synthesized aza analogues. The role of Ser 7.36 was confirmed as the affinity of aza analogues for the mutant D4 receptor S7.36A was reduced. Then, the formation of another hydrogen bond with the conserved Ser 5.42 residue appeared to be also critical for D4 binding.

Receptor Activity-modifying Proteins 2 and 3 Generate Adrenomedullin Receptor Subtypes with Distinct Molecular Properties

Adrenomedullin (AM) is a peptide hormone with numerous effects in the vascular systems. AM signals through the AM1 and AM2 receptors formed by the obligate heterodimerization of a G protein-coupled receptor, the calcitonin receptor-like receptor (CLR), and receptor activity-modifying proteins 2 and 3 (RAMP2 and RAMP3), respectively. These different CLR-RAMP interactions yield discrete receptor pharmacology and physiological effects. The effective design of therapeutics that target the individual AM receptors is dependent on understanding the molecular details of the effects of RAMPs on CLR. To understand the role of RAMP2 and -3 on the activation and conformation of the CLR subunit of AM receptors, we mutated 68 individual amino acids in the juxtamembrane region of CLR, a key region for activation of AM receptors, and determined the effects on cAMP signaling. Sixteen CLR mutations had differential effects between the AM1 and AM2 receptors. Accompanying this, independent molecular modeling of the full-length AM-bound AM1 and AM2 receptors predicted differences in the binding pocket and differences in the electrostatic potential of the two AM receptors. Druggability analysis indicated unique features that could be used to develop selective small molecule ligands for each receptor. The interaction of RAMP2 or RAMP3 with CLR induces conformational variation in the juxtamembrane region, yielding distinct binding pockets, probably via an allosteric mechanism. These subtype-specific differences have implications for the design of therapeutics aimed at specific AM receptors and for understanding the mechanisms by which accessory proteins affect G protein-coupled receptor function.

Bitter taste receptors: Novel insights into the biochemistry and pharmacology

Bitter taste receptors (T2Rs) belong to the super family of G protein-coupled receptors (GPCRs). There are 25 T2Rs expressed in humans, and these interact with a large and diverse group of bitter ligands. T2Rs are expressed in many extra-oral tissues and can perform diverse physiological roles. Structure-function studies led to the identification of similarities and dissimilarities between T2Rs and Class A GPCRs including amino acid conservation and novel motifs. However, the efficacy of most of the T2R ligands is not yet elucidated and the biochemical pharmacology of T2Rs is poorly understood. Recent studies on T2Rs characterized novel ligands including blockers for these receptors that include inverse agonist and antagonists. In this review we discuss the techniques used for elucidating bitter blockers, concept of ligand bias, generic amino acid numbering, the role of cholesterol, and conserved water molecules in the biochemistry and pharmacology of T2Rs.

Synthesis and evaluation of aporphine analogs containing C1 allyl isosteres at the h5-HT(2A) receptor

A series of C1 aporphine analogs related to compound 5 and that contain substituted allylic, alkynyl, nitrile, ester and benzyl groups was synthesized and evaluated for affinity at h5HT2A and α1A receptors in functional activity assays that measure calcium release. The presence of branched allylic substituent groups diminished affinity for the h5HT2A receptor. Likewise, the alkynyl, nitrile and ester derivatives evaluated displayed lower 5-HT2A receptor affinity as compared to 5. Hydrophobic, steric and electronic effects impact the affinity of p-substituted benzyl derivatives 8i-8k but in different ways. High hydrophobicity and size favor 5-HT2A affinity whereas, high electronegativity disfavors 5-HT2A affinity. p-Bromobenzyl analog 8k was identified as a 5-HT2A receptor selective ligand, with the highest 5-HT2A receptor affinity of any aporphine known to date. Most of the other analogs were selective for the 5-HT2A versus the α1A receptor. ChemScore binding energies from docking studies correlated qualitatively with the observed trends in affinity for 8i-8k, although the binding energies were not well differentiated quantitatively.

Determination of the binding mode for the cyclopentapeptide CXCR4 antagonist FC131 using a dual approach of ligand modifications and receptor mutagenesis

BACKGROUND AND PURPOSE

The cyclopentapeptide FC131 (cyclo(-L-Arg(1) -L-Arg(2) -L-2-Nal(3) -Gly(4) -D-Tyr(5) -)) is an antagonist at the CXC chemokine receptor CXCR4, which plays a role in human immunodeficiency virus infection, cancer and stem cell recruitment. Binding modes for FC131 in CXCR4 have previously been suggested based on molecular docking guided by structure-activity relationship (SAR) data; however, none of these have been verified by in vitro experiments.

EXPERIMENTAL APPROACH

Heterologous (125) I-12G5-competition binding and functional assays (inhibition of CXCL12-mediated activation) of FC131 and three analogues were performed on wild-type CXCR4 and 25 receptor mutants. Computational modelling was used to rationalize the experimental data.

KEY RESULTS

The Arg(2) and 2-Nal(3) side chains of FC131 interact with residues in TM-3 (His(113) , Asp(171) ) and TM-5 (hydrophobic pocket) respectively. Arg(1) forms charge-charge interactions with Asp(187) in ECL-2, while D-Tyr(5) points to the extracellular side of CXCR4. Furthermore, the backbone of FC131 interacts with the chemokine receptor-conserved Glu(288) via two water molecules. Intriguingly, Tyr(116) and Glu(288) form a H-bond in CXCR4 crystal structures and mutation of either residue to Ala abolishes CXCR4 activity.

CONCLUSIONS AND IMPLICATIONS

Ligand modification, receptor mutagenesis and computational modelling approaches were used to identify the binding mode of FC131 in CXCR4, which was in agreement with binding modes suggested from previous SAR studies. Furthermore, insights into the mechanism for CXCR4 activation by CXCL12 were gained. The combined findings will facilitate future design of novel CXCR4 antagonists.

In silico characterization of binding mode of CCR8 inhibitor: homology modeling, docking and membrane based MD simulation study

Human CC-chemokine receptor 8 (CCR8) is a crucial drug target in asthma that belongs to G-protein-coupled receptor superfamily, which is characterized by seven transmembrane helices. To date, there is no X-ray crystal structure available for CCR8; this hampers active research on the target. Molecular basis of interaction mechanism of antagonist with CCR8 remains unclear. In order to provide binding site information and stable binding mode, we performed modeling, docking and molecular dynamics (MD) simulation of CCR8. Docking study of biaryl-ether-piperidine derivative (13C) was performed inside predefined CCR8 binding site to get the representative conformation of 13C. Further, MD simulations of receptor and complex (13C-CCR8) inside dipalmitoylphosphatidylcholine lipid bilayers were performed to explore the effect of lipids. Results analyses showed that the Gln91, Tyr94, Cys106, Val109, Tyr113, Cys183, Tyr184, Ser185, Lys195, Thr198, Asn199, Met202, Phe254, and Glu286 were conserved in both docking and MD simulations. This indicated possible role of these residues in CCR8 antagonism. However, experimental mutational studies on these identified residues could be effective to confirm their importance in CCR8 antagonism. Furthermore, calculated Coulombic interactions represented the crucial roles of Glu286, Lys195, and Tyr113 in CCR8 antagonism. Important residues identified in this study overlap with the previous non-peptide agonist (LMD-009) binding site. Though, the non-peptide agonist and currently studied inhibitor (13C) share common substructure, but they differ in their effects on CCR8. So, to get more insight into their agonist and antagonist effects, further side-by-side experimental studies on both agonist (LMD-009) and antagonist (13C) are suggested.

Generic GPCR residue numbers - aligning topology maps while minding the gaps

Generic residue numbers facilitate comparisons of, for example, mutational effects, ligand interactions, and structural motifs. The numbering scheme by Ballesteros and Weinstein for residues within the class A GPCRs (G protein-coupled receptors) has more than 1100 citations, and the recent crystal structures for classes B, C, and F now call for a community consensus in residue numbering within and across these classes. Furthermore, the structural era has uncovered helix bulges and constrictions that offset the generic residue numbers. The use of generic residue numbers depends on convenient access by pharmacologists, chemists, and structural biologists. We review the generic residue numbering schemes for each GPCR class, as well as a complementary structure-based scheme, and provide illustrative examples and GPCR database (GPCRDB) web tools to number any receptor sequence or structure.

Improving homology modeling of G-protein coupled receptors through multiple-template derived conserved inter-residue interactions

Evidenced by the three-rounds of G-protein coupled receptors (GPCR) Dock competitions, improving homology modeling methods of helical transmembrane proteins including the GPCRs, based on templates of low sequence identity, remains an eminent challenge. Current approaches addressing this challenge adopt the philosophy of "modeling first, refinement next". In the present work, we developed an alternative modeling approach through the novel application of available multiple templates. First, conserved inter-residue interactions are derived from each additional template through conservation analysis of each template-target pairwise alignment. Then, these interactions are converted into distance restraints and incorporated in the homology modeling process. This approach was applied to modeling of the human β2 adrenergic receptor using the bovin rhodopsin and the human protease-activated receptor 1 as templates and improved model quality was demonstrated compared to the homology model generated by standard single-template and multiple-template methods. This method of "refined restraints first, modeling next", provides a fast and complementary way to the current modeling approaches. It allows rational identification and implementation of additional conserved distance restraints extracted from multiple templates and/or experimental data, and has the potential to be applicable to modeling of all helical transmembrane proteins.

Identification and functional comparison of seven-transmembrane G-protein-coupled BILF1 receptors in recently discovered nonhuman primate lymphocryptoviruses

Coevolution of herpesviruses with their respective host has resulted in a delicate balance between virus-encoded immune evasion mechanisms and host antiviral immunity. BILF1 encoded by human Epstein-Barr virus (EBV) is a 7-transmembrane (7TM) G-protein-coupled receptor (GPCR) with multiple immunomodulatory functions, including attenuation of PKR phosphorylation, activation of G-protein signaling, and downregulation of major histocompatibility complex (MHC) class I surface expression. In this study, we explored the evolutionary and functional relationships between BILF1 receptor family members from EBV and 12 previously uncharacterized nonhuman primate (NHP) lymphocryptoviruses (LCVs). Phylogenetic analysis defined 3 BILF1 clades, corresponding to LCVs of New World monkeys (clade A) or Old World monkeys and great apes (clades B and C). Common functional properties were suggested by a high degree of sequence conservation in functionally important regions of the BILF1 molecules. A subset of BILF1 receptors from EBV and LCVs from NHPs (chimpanzee, orangutan, marmoset, and siamang) were selected for multifunctional analysis. All receptors exhibited constitutive signaling activity via G protein Gαi and induced activation of the NF-κB transcription factor. In contrast, only 3 of 5 were able to activate NFAT (nuclear factor of activated T cells); chimpanzee and orangutan BILF1 molecules were unable to activate NFAT. Similarly, although all receptors were internalized, BILF1 from the chimpanzee and orangutan displayed an altered cellular localization pattern with predominant cell surface expression. This study shows how biochemical characterization of functionally important orthologous viral proteins can be used to complement phylogenetic analysis to provide further insight into diverse microbial evolutionary relationships and immune evasion function.

IMPORTANCE

Epstein-Barr virus (EBV), known as an oncovirus, is the only human herpesvirus in the genus Lymphocryptovirus (LCV). EBV uses multiple strategies to hijack infected host cells, establish persistent infection in B cells, and evade antiviral immune responses. As part of EBV's immune evasion strategy, the virus encodes a multifunctional 7-transmembrane (7TM) G-protein-coupled receptor (GPCR), EBV BILF1. In addition to multiple immune evasion-associated functions, EBV BILF1 has transforming properties, which are linked to its high constitutive activity. We identified BILF1 receptor orthologues in 12 previously uncharacterized LCVs from nonhuman primates (NHPs) of Old and New World origin. As 7TM receptors are excellent drug targets, our unique insight into the molecular mechanism of action of the BILF1 family and into the evolution of primate LCVs may enable validation of EBV BILF1 as a drug target for EBV-mediated diseases, as well as facilitating the design of drugs targeting EBV BILF1.

Insights into the binding modes of CC chemokine receptor 4 (CCR4) inhibitors: a combined approach involving homology modelling, docking, and molecular dynamics simulation studies

CC chemokine receptor 4 (CCR4), a G protein-coupled receptor (GPCR), plays a vital role in the progression of asthma, T-cell lymphoma, inflammation, and Alzheimer's disease. To date, the structure of CCR4 has not been determined. Therefore, the nature of the interactions between inhibitors and CCR4 is not well known. In this study, we used CCR5 as a template to model the structure of CCR4. Docking studies were performed for four naphthalene-sulphonamide derivatives and crucial ligand-protein interactions were analysed. Molecular dynamics (MD) simulations of these complexes (100 ns each) were carried out to gain insights into the interactions between ligands and CCR4. MD simulations revealed that the residues identified by the docking were displaced and new residues were inserted near the ligands. Results of a principal component analysis (PCA) suggested that CCR4 unfolds at the extracellular site surrounding the ligands. Our simulations identified crucial residues involved in CCR4 antagonism, which were supported by previous mutational studies. Additionally, we identified Ser3.29, Leu3.33, Ser5.39, Phe6.47, Ile7.35, Thr7.38, Thr7.40, and Ala7.42 as residues that play crucial roles in CCR4 antagonism. Mutational studies will help elucidate the significance of these residues in CCR4 antagonism. An understanding of ligand-CCR4 interactions might aid in the design of novel CCR4 inhibitors.

Molecular basis of cannabinoid CB1 receptor coupling to the G protein heterotrimer Gαiβγ: identification of key CB1 contacts with the C-terminal helix α5 of Gαi

The cannabinoid (CB1) receptor is a member of the rhodopsin-like G protein-coupled receptor superfamily. The human CB1 receptor, which is among the most expressed receptors in the brain, has been implicated in several disease states, including drug addiction, anxiety, depression, obesity, and chronic pain. Different classes of CB1 agonists evoke signaling pathways through the activation of specific subtypes of G proteins. The molecular basis of CB1 receptor coupling to its cognate G protein is unknown. As a first step toward understanding CB1 receptor-mediated G protein signaling, we have constructed a ternary complex structural model of the CB1 receptor and Gi heterotrimer (CB1-Gi), guided by the x-ray structure of β2-adrenergic receptor (β2AR) in complex with Gs (β2AR-Gs), through 824-ns duration molecular dynamics simulations in a fully hydrated 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine bilayer environment. We identified a group of residues at the juxtamembrane regions of the intracellular loops 2 and 3 (IC2 and IC3) of the CB1 receptor, including Ile-218(3.54), Tyr-224(IC2), Asp-338(6.30), Arg-340(6.32), Leu-341(6.33), and Thr-344(6.36), as potential key contacts with the extreme C-terminal helix α5 of Gαi. Ala mutations of these residues at the receptor-Gi interface resulted in little G protein coupling activity, consistent with the present model of the CB1-Gi complex, which suggests tight interactions between CB1 and the extreme C-terminal helix α5 of Gαi. The model also suggests that unique conformational changes in the extreme C-terminal helix α5 of Gα play a crucial role in the receptor-mediated G protein activation.

Structural insight into a novel human CCR5-V130I variant associated with resistance to HIV-1 infection

We report the identification of a novel CC chemokine receptor 5 (CCR5) variant that seems associated with resistance to HIV-1 infection. The V130I mutation of the CCR5 receptor is located in the intracellular loop ICL2 known as DRY box and described in the literature as a nonsynonymous mutation present in nonhuman primates group. Extensive molecular modeling and dynamics simulations were performed to elucidate the mechanism by which the V130I mutation may induce conformational change of the CCR5 folding protein and prevent the interaction with the β-arrestin protein. Our study provides new mechanistic insight into how a specific mutation in the regulatory domain of CCR5 might alter the structural folding of the DRY box and the possible ICL2 loop binding with the β-arrestin protein, as described in our previous computational study. The results from our large-scale simulations complement recent experimental results and clinical features and offer useful insights into the mechanism behind CCR5 protein folding and signal transduction. In order for HIV, the entry of the virus to the cells must fuse with the CCR5 receptor that sits on the surface of T-helper immune cells. The described V130I mutation in the gene encoding the CCR5 protein may results in a defective CCR5-Arrestin binding complex that blocks entry of the virus.

On homology modeling of the M₂ muscarinic acetylcholine receptor subtype

Twelve homology models of the human M2 muscarinic receptor using different sets of templates have been designed using the Prime program or the modeller program and compared to crystallographic structure (PDB:3UON). The best models were obtained using single template of the closest published structure, the M3 muscarinic receptor (PDB:4DAJ). Adding more (structurally distant) templates led to worse models. Data document a key role of the template in homology modeling. The models differ substantially. The quality checks built into the programs do not correlate with the RMSDs to the crystallographic structure and cannot be used to select the best model. Re-docking of the antagonists present in crystallographic structure and relative binding energy estimation by calculating MM/GBSA in Prime and the binding energy function in YASARA suggested it could be possible to evaluate the quality of the orthosteric binding site based on the prediction of relative binding energies. Although estimation of relative binding energies distinguishes between relatively good and bad models it does not indicate the best one. On the other hand, visual inspection of the models for known features and knowledge-based analysis of the intramolecular interactions allows an experimenter to select overall best models manually.

Functional residues essential for the activation of the CB1 cannabinoid receptor

Recently developed X-ray crystal structures of active state G-protein-coupled receptors have opened the way for detailed examination of the movement of the transmembrane (TM) helices and the specific residues involved in the receptor activation upon ligand binding. Previous modeling studies have indicated that the brain cannabinoid (CB1) receptor binds with a ligand at least in part through a hydrophobic tail on the ligand. This interaction is believed to be similar to the rotameric toggle switch proposed for the β2 adrenergic receptor (β2AR). In the present study, an active state model for the CB1 receptor, guided by the X-ray structure of the active state for β2AR, was constructed with HU210 bound as a ligand. Molecular dynamics (MD) simulations were employed to provide a smooth progression between inactive and active states of the receptor. This model was compared with our previously published CB1 receptor model to identify the functional residues that play key roles in triggering receptor conformational changes upon agonist binding. Movements in TM helices and functional residues are discussed. W279(5.43), contributing to an inward movement of the fifth TM helix (TM5) to the helical core, could serve as another rotameric switch for receptor activation. V282(5.46), interacting with the ligand's hydrophobic C3 alkyl chain, appears to play a key role in TM5 inward movement centered at L286(5.50) and subsequent coupling to V204(3.40). V204(3.40), closely interacting with the TM5 and TM6 hydrophobic patch residues in the middle of the receptor, particularly I290(5.54) and L352(6.44), appears to facilitate helical rearrangements, leading to the breakage of the ionic lock and the rotameric change of Y397(7.53), which are key features of the active state.

Structural insight into the binding complex: β-arrestin/CCR5 complex

The chemokine receptor 5 (CCR5) belongs to the superfamily of serpentine G protein-coupled receptors (GPCRs). The DRY motif (Asp, Arg, Tyr) of the intracellular loop 2 (ICL2), which is highly conserved in the GPCRs has been shown to be essential for the stability of folding of CCR5 and the interaction with β-arrestin. But the molecular mechanism by which it recognizes and interacts with β-arrestin has not been elucidated. In the present study, we described the active state of the β-arrestin structure using normal mode analysis and characterized the binding cleft of CCR5-ICL2 with β-arrestin using SABRE© docking tool and molecular dynamics simulation. Based on our computational results, we proposed a mode of binding between the ICL2 loop of CCR5 and β-arrestin structure, and modeled the energetically stable β-arrestin/CCR5 complex. In view of CCR5's importance as a therapeutic target for the treatment of HIV, this observation provides novel insight into the β-arrestin/CCR5 pathway. As a result, the current computational study of the detailed β-arrestin/CCR5 binding complex could provide the rationale for the development of next generation of HIV peptide inhibitors as therapeutic agents.

Identifying subset errors in multiple sequence alignments

Multiple sequence alignment (MSA) accuracy is important, but there is no widely accepted method of judging the accuracy that different alignment algorithms give. We present a simple approach to detecting two types of error, namely block shifts and the misplacement of residues within a gap. Given a MSA, subsets of very similar sequences are generated through the use of a redundancy filter, typically using a 70-90% sequence identity cut-off. Subsets thus produced are typically small and degenerate, and errors can be easily detected even by manual examination. The errors, albeit minor, are inevitably associated with gaps in the alignment, and so the procedure is particularly relevant to homology modelling of protein loop regions. The usefulness of the approach is illustrated in the context of the universal but little known [K/R]KLH motif that occurs in intracellular loop 1 of G protein coupled receptors (GPCR); other issues relevant to GPCR modelling are also discussed.

Structural features of the apelin receptor N-terminal tail and first transmembrane segment implicated in ligand binding and receptor trafficking

G-protein coupled receptors (GPCRs) comprise a large family of membrane proteins with rich functional diversity. Signaling through the apelin receptor (AR or APJ) influences the cardiovascular system, central nervous system and glucose regulation. Pathophysiological involvement of apelin has been shown in atherosclerosis, cancer, human immunodeficiency virus-1 (HIV-1) infection and obesity. Here, we present the high-resolution nuclear magnetic resonance (NMR) spectroscopy-based structure of the N-terminus and first transmembrane (TM) segment of AR (residues 1-55, AR55) in dodecylphosphocholine micelles. AR55 consists of two disrupted helices, spanning residues D14-K25 and A29-R55(1.59). Molecular dynamics (MD) simulations of AR built from a hybrid of experimental NMR and homology model-based restraints allowed validation of the AR55 structure in the context of the full-length receptor in a hydrated bilayer. AR55 structural features were functionally probed using mutagenesis in full-length AR through monitoring of apelin-induced extracellular signal-regulated kinase (ERK) phosphorylation in transiently transfected human embryonic kidney (HEK) 293A cells. Residues E20 and D23 form an extracellular anionic face and interact with lipid headgroups during MD simulations in the absence of ligand, producing an ideal binding site for a cationic apelin ligand proximal to the membrane-water interface, lending credence to membrane-catalyzed apelin-AR binding. In the TM region of AR55, N46(1.50) is central to a disruption in helical character. G42(1.46), G45(1.49) and N46(1.50), which are all involved in the TM helical disruption, are essential for proper trafficking of AR. In summary, we introduce a new correlative NMR spectroscopy and computational biochemistry methodology and demonstrate its utility in providing some of the first high-resolution structural information for a peptide-activated GPCR TM domain.

A Strategy Combining Differential Low-Throughput Screening and Virtual Screening (DLS-VS) Accelerating the Discovery of new Modulators for the Orphan GPR34 Receptor

The DLS-VS strategy was developed as an integrated method for identifying chemical modulators for orphan GPCRs. It combines differential low-throughput screening (DLS) and virtual screening (VS). The two cascaded techniques offer complementary advantages and allow the experimental testing of a minimal number of compounds. First, DLS identifies modulators specific for the considered receptor among a set of receptors, through the screening of a small library with diverse chemical compounds. Then, an active molecular model of the receptor is built by homology to a validated template, and it is progressively refined by rotamers modification for key side-chains, by VS of the already screened library, and by iterative selection of the model generating the best enrichment. The refined active model is finally used for the VS of a large chemical library and the selection of a small set of compounds for experimental testing. Applied to the orphan receptor GPR34, the DLS-VS strategy combined the experimental screening of 20 000 compounds and the virtual screening of 1 250 000 compounds. It identified one agonist and eight inverse agonists, showing a high chemical diversity. We describe the method. The strategy can be applied to other GPCRs.

Asymmetry of the rhodopsin dimer in complex with transducin

A large body of evidence for G-protein-coupled receptor (GPCR) oligomerization has accumulated over the past 2 decades. The smallest of these oligomers in vivo most likely is a dimer that buries 1000-Å(2) intramolecular surfaces and on stimulation forms a complex with heterotrimeric G protein in 2:1 stoichiometry. However, it is unclear whether each of the monomers adopts the same or a different conformation and function after activation of this dimer. With bovine rhodopsin (Rho) and its cognate bovine G-protein transducin (Gt) as a model system, we used the retinoid chromophores 11-cis-retinal and 9-cis-retinal to monitor each monomer of the dimeric GPCR within a stable complex with nucleotide-free Gt. We found that only 50% of Rho* in the Rho*-Gt complex is trapped in a Meta II conformation, while 50% evolves toward an opsin conformation and can be regenerated with 9-cis-retinal. We also found that all-trans-retinal can regenerate chromophore-depleted Rho*e complexed with Gt and FAK*TSA peptide containing Lys(296) with the attached all-trans retinoid (m/z of 934.5[MH](+)) was identified by mass spectrometry. Thus, our study shows that each of the monomers contributes unequally to the pentameric (2:1:1:1) complex of Rho dimer and Gt heterotrimer, validating the oligomeric structure of the complex and the asymmetry of the GPCR dimer, and revealing its structural/functional signature. This study provides a clear functional distinction between monomers of family A GPCRs in their oligomeric form.

Similarity between class A and class B G-protein-coupled receptors exemplified through calcitonin gene-related peptide receptor modelling and mutagenesis studies

Modelling class B G-protein-coupled receptors (GPCRs) using class A GPCR structural templates is difficult due to lack of homology. The plant GPCR, GCR1, has homology to both class A and class B GPCRs. We have used this to generate a class A–class B alignment, and by incorporating maximum lagged correlation of entropy and hydrophobicity into a consensus score, we have been able to align receptor transmembrane regions. We have applied this analysis to generate active and inactive homology models of the class B calcitonin gene-related peptide (CGRP) receptor, and have supported it with site-directed mutagenesis data using 122 CGRP receptor residues and 144 published mutagenesis results on other class B GPCRs. The variation of sequence variability with structure, the analysis of polarity violations, the alignment of group-conserved residues and the mutagenesis results at 27 key positions were particularly informative in distinguishing between the proposed and plausible alternative alignments. Furthermore, we have been able to associate the key molecular features of the class B GPCR signalling machinery with their class A counterparts for the first time. These include the [K/R]KLH motif in intracellular loop 1, [I/L]xxxL and KxxK at the intracellular end of TM5 and TM6, the NPXXY/VAVLY motif on TM7 and small group-conserved residues in TM1, TM2, TM3 and TM7. The equivalent of the class A DRY motif is proposed to involve Arg^2.39, His^2.43 and Glu^3.46, which makes a polar lock with T^6.37. These alignments and models provide useful tools for understanding class B GPCR function.

Restructuring G-protein- coupled receptor activation

G-protein-coupled receptors serve as key signal transduction conduits, linking extracellular inputs with diverse cellular responses. These receptors eluded structural characterization for decades following their identification. A landmark structure of rhodopsin provided a basis for structure-function studies and homology modeling, but advances in receptor biology suffered from a lack of receptor-specific structural insights. The recent explosion in GPCR structures confirms some features predicted by rhodopsin-based models, and more importantly, it reveals unexpected ligand-binding modes and critical aspects of the receptor activation process. The new structures also promise to foster studies testing emerging models for GPCR function such as receptor dimerization and ligand-biased signaling.

Modeling of mammalian olfactory receptors and docking of odorants

Olfactory receptors (ORs) belong to the superfamily of G protein-coupled receptors (GPCRs), the second largest class of genes after those related to immunity, and account for about 3 % of mammalian genomes. ORs are present in all multicellular organisms and represent more than half the GPCRs in mammalian species (e.g., the mouse OR repertoire contains >1,000 functional genes). ORs are mainly expressed in the olfactory epithelium where they detect odorant molecules, but they are also expressed in a number of other cells, such as sperm cells, although their functions in these cells remain mostly unknown. It has recently been reported that ORs are present in tumoral tissues where they are expressed at different levels than in healthy tissues. A specific OR is over-expressed in prostate cancer cells, and activation of this OR has been shown to inhibit the proliferation of these cells. Odorant stimulation of some of these receptors results in inhibition of cell proliferation. Even though their biological role has not yet been elucidated, these receptors might constitute new targets for diagnosis and therapeutics. It is important to understand the activation mechanism of these receptors at the molecular level, in particular to be able to predict which ligands are likely to activate a particular receptor ('deorphanization') or to design antagonists for a given receptor. In this review, we describe the in silico methodologies used to model the three-dimensional (3D) structure of ORs (in the more general framework of GPCR modeling) and to dock ligands into these 3D structures.

Leu128(3.43) (l128) and Val247(6.40) (V247) of CXCR1 are critical amino acid residues for g protein coupling and receptor activation

CXCR1, a classic GPCR that binds IL-8, plays a key role in neutrophil activation and migration by activating phospholipase C (PLC)β through Gα₁₅ and Gα_i which generates diacylglycerol and inositol phosphates (IPs). In this study, two conserved amino acid residues of CXCR1 on the transmembrane domain (TM) 3 and TM6, Leu128^3.43 (L128) and Val247^6.40 (V247), respectively, were selectively substituted with other amino acids to investigate the role of these conserved residues in CXCR1 activation. Although two selective mutants on Leu128, Leu128Ala (L128A) and Leu128Arg (L128R), demonstrated high binding affinity to IL-8, they were not capable of coupling to G proteins and consequently lost the functional response of the receptors. By contrast, among the four mutants at residue Val247 (TM6.40), replacing Val247 with Ala (V247A) and Asn (V247N) led to constitutive activation of mutant receptors when cotransfected with Gα₁₅. The V247N mutant also constitutively activated the Gα_i protein. These results indicate that L128 on TM3.43 is involved in G protein coupling and receptor activation but is unimportant for ligand binding. On the other hand, V247 on TM6.40 plays a critical role in maintaining the receptor in the inactive state, and the substitution of V247 impaired the receptor constraint and stabilized an active conformation. Functionally, there was an increase in chemotaxis in response to IL-8 in cells expressing V247A and V247N. Our findings indicate that Leu128^3.43 and Val247^6.40 are critical for G protein coupling and activation of signaling effectors, providing a valuable insight into the mechanism of CXCR1 activation.

Identification of structural determinants of ligand selectivity in 5-HT₂ receptor subtypes on the basis of protein-ligand interactions

Drug selectivity is one of the most critical improvement steps in drug development. The 5-hydroxytryptamine 2 (5-HT₂) receptor has 3 subtypes that exhibit different pharmacological functions. Because of their high amino acid sequence similarity, designing small molecules that selectively activate only 1 receptor among the 3 subtypes is difficult. We performed homology modeling of the 5-HT₂ receptor subtypes using the β₂-adrenergic receptor as a template to identify differences in active sites that may influence 5-HT₂ receptor agonist selectivity. A subset of selective 5-HT₂ agonists was docked into the modeled protein structures to investigate their interactions with each receptor. Subtype-specific active site residues at positions xl2.54, 5.39, and 5.46 interacted differently with each ligand. Molecular dynamics simulations revealed that position 5.46 of the 5-HT(2A) receptor interacted more favorably with selective 5-HT(2A) agonists than with selective 5-HT(2B) agonists. These computationally obtained insights provided clues to improving agonist selectivity for specific pharmacological action at 5-HT₂ receptors.

The significance of G protein-coupled receptor crystallography for drug discovery

Crucial as molecular sensors for many vital physiological processes, seven-transmembrane domain G protein-coupled receptors (GPCRs) comprise the largest family of proteins targeted by drug discovery. Together with structures of the prototypical GPCR rhodopsin, solved structures of other liganded GPCRs promise to provide insights into the structural basis of the superfamily's biochemical functions and assist in the development of new therapeutic modalities and drugs. One of the greatest technical and theoretical challenges to elucidating and exploiting structure-function relationships in these systems is the emerging concept of GPCR conformational flexibility and its cause-effect relationship for receptor-receptor and receptor-effector interactions. Such conformational changes can be subtle and triggered by relatively small binding energy effects, leading to full or partial efficacy in the activation or inactivation of the receptor system at large. Pharmacological dogma generally dictates that these changes manifest themselves through kinetic modulation of the receptor's G protein partners. Atomic resolution information derived from increasingly available receptor structures provides an entrée to the understanding of these events and practically applying it to drug design. Supported by structure-activity relationship information arising from empirical screening, a unified structural model of GPCR activation/inactivation promises to both accelerate drug discovery in this field and improve our fundamental understanding of structure-based drug design in general. This review discusses fundamental problems that persist in drug design and GPCR structural determination.

In silico investigation of interactions between human cannabinoid receptor-1 and its antagonists

Cannabinoid receptor-1 (CB(1)) is widely expressed in the central nervous system and plays a vital role in regulating food intake and energy expenditure. CB(1) antagonists such as Rimonabant have been used in clinic to inhibit food intake, and therefore reduce body weight in obese animals and humans. To investigate the binding modes of CB(1) antagonists to the receptor, both receptor- and ligand-based methods were implemented in this study. At first, a pharmacophore model was generated based on 31 diverse CB(1) antagonists collected from literature. A test set validation and a simulated virtual screening evaluation were then performed to verify the reliability and discriminating ability of the pharmacophore. Meanwhile, the homology model of CB(1) receptor was constructed based on the crystal structure of human β (2) adrenergic receptor (β (2)-AR). Several classical antagonists were then docked into the optimized homology model with induced fit docking method. A hydrogen bond between the antagonists and Lys192 on the third transmembrane helix of the receptor was formed in the docking study, which has proven to be critical for receptor-ligand interaction by biological experiments. The structure obtained from induced fit docking was then confirmed to be a reliable model for molecular docking from the result of the simulated virtual screening. The consistency between the pharmacophore and the homology structure further proved the previous observation. The built receptor structure and antagonists' pharmacophore should be useful for the understanding of inhibitory mechanism and development of novel CB(1) antagonists.

5-HT₂c receptor selectivity and structure-activity relationship of N-methyl-N-(1-methylpiperidin-4-yl)benzenesulfonamide analogs

Agonists of the 5-HT(2C) receptor have attracted much attention as therapeutic agents for the treatment of obesity. Subtype selectivity against other 5-HT(2) receptors is one of the most important prerequisites for reducing side effects. We present the synthesis of N-methyl-N-(1-methylpiperidin-4-yl)benzenesulfonamide analogs and their structure-activity relationship studies on 5-HT(2A) and 5-HT(2C) receptors. Although the compounds showed nanomolar activity to the 5-HT(2C) receptor, their selectivity against the 5-HT(2A) receptor was modest to low. Molecular modeling studies using homology modeling and docking simulation revealed that selectivity originated from subtype specific residues. The observed binding modes and receptor-ligand interactions provided us a clue for optimizing the selectivity against the 5-HT(2A) receptor.