Ligand- and mutation-induced conformational selection in the CCR5 chemokine G protein-coupled receptor

The chemokine receptor CCR5: multi-faceted hook for HIV-1

Chemokines are cytokines whose primary role is cellular activation and stimulation of leukocyte migration. They perform their various functions by interacting with G protein-coupled cell surface receptors (GPCRs) and are involved in the regulation of many biological processes such as apoptosis, proliferation, angiogenesis, hematopoiesis or organogenesis. They contribute to the maintenance of the homeostasis of lymphocytes and coordinate the function of the immune system. However, chemokines and their receptors are sometimes hijacked by some pathogens to infect the host organism. For a given chemokine receptor, there is a wide structural, organizational and conformational diversity. In this review, we describe the evidence for structural variety reported for the chemokine receptor CCR5, how this variability can be exploited by HIV-1 to infect its target cells and what therapeutic solutions are currently being developed to overcome this problem.

Intermediate-state-trapped mutants pinpoint G protein-coupled receptor conformational allostery

Understanding the roles of intermediate states in signaling is pivotal to unraveling the activation processes of G protein-coupled receptors (GPCRs). However, the field is still struggling to define these conformational states with sufficient resolution to study their individual functions. Here, we demonstrate the feasibility of enriching the populations of discrete states via conformation-biased mutants. These mutants adopt distinct distributions among five states that lie along the activation pathway of adenosine A2A receptor (A2AR), a class A GPCR. Our study reveals a structurally conserved cation-π lock between transmembrane helix VI (TM6) and Helix8 that regulates cytoplasmic cavity opening as a "gatekeeper" for G protein penetration. A GPCR activation process based on the well-discerned conformational states is thus proposed, allosterically micro-modulated by the cation-π lock and a previously well-defined ionic interaction between TM3 and TM6. Intermediate-state-trapped mutants will also provide useful information in relation to receptor-G protein signal transduction.

Tracking receptor motions at the plasma membrane reveals distinct effects of ligands on CCR5 dynamics depending on its dimerization status

G-protein-coupled receptors (GPCR) are present at the cell surface in different conformational and oligomeric states. However, how these states impact GPCRs biological function and therapeutic targeting remains incompletely known. Here, we investigated this issue in living cells for the CC chemokine receptor 5 (CCR5), a major receptor in inflammation and the principal entry co-receptor for Human Immunodeficiency Viruses type 1 (HIV-1). We used TIRF microscopy and a statistical method to track and classify the motion of different receptor subpopulations. We showed a diversity of ligand-free forms of CCR5 at the cell surface constituted of various oligomeric states and exhibiting transient Brownian and restricted motions. These forms were stabilized differently by distinct ligands. In particular, agonist stimulation restricted the mobility of CCR5 and led to its clustering, a feature depending on β-arrestin, while inverse agonist stimulation exhibited the opposite effect. These results suggest a link between receptor activation and immobilization. Applied to HIV-1 envelope glycoproteins gp120, our quantitative analysis revealed agonist-like properties of gp120s. Distinct gp120s influenced CCR5 dynamics differently, suggesting that they stabilize different CCR5 conformations. Then, using a dimerization-compromized mutant, we showed that dimerization (i) impacts CCR5 precoupling to G proteins, (ii) is a pre-requisite for the immobilization and clustering of receptors upon activation, and (iii) regulates receptor endocytosis, thereby impacting the fate of activated receptors. This study demonstrates that tracking the dynamic behavior of a GPCR is an efficient way to link GPCR conformations to their functions, therefore improving the development of drugs targeting specific receptor conformations.

Development of enhanced conformational sampling methods to probe the activation landscape of GPCRs

G protein-coupled receptors (GPCRs) make up the largest superfamily of integral membrane proteins and play critical signal transduction roles in many physiological processes. Developments in molecular biology, biophysical, biochemical, pharmacological, and computational techniques aimed at these important therapeutic targets are beginning to provide unprecedented details on the structural as well as functional basis of their pleiotropic signaling mediated by G proteins, β arrestins, and other transducers. This pleiotropy presents a pharmacological challenge as the same ligand-receptor interaction can cause a therapeutic effect as well as an undesirable on-target side-effect through different downstream pathways. GPCRs don't function as simple binary on-off switches but as finely tuned shape-shifting machines described by conformational ensembles, where unique subsets of conformations may be responsible for specific signaling cascades. X-ray crystallography and more recently cryo-electron microscopy are providing snapshots of some of these functionally-important receptor conformations bound to ligands and/or transducers, which are being utilized by computational methods to describe the dynamic conformational energy landscape of GPCRs. In this chapter, we review the progress in computational conformational sampling methods based on molecular dynamics and discrete sampling approaches that have been successful in complementing biophysical and biochemical studies on these receptors in terms of their activation mechanisms, allosteric effects, actions of biased ligands, and effects of pathological mutations. Some of the sampled simulation time scales are beginning to approach receptor activation time scales. The list of conformational sampling methods and example uses discussed is not exhaustive but includes representative examples that have pushed the limits of classical molecular dynamics and discrete sampling methods to describe the activation energy landscape of GPCRs.

Synergic Effects in the Activation of the Sweet Receptor GPCR Heterodimer for Various Sweeteners Predicted Using Molecular Metadynamics Simulations

The sweet taste is elicited by activation of the TAS1R2/1R3 heterodimer G protein-coupled receptor. This is a therapeutic target for treatment of obesity and metabolic dysfunctions. Sweetener blends provide attractive strategies to lower the sugar level while preserving the attractive taste of food. To understand the synergic effect of various sweetener blend combinations of artificial and natural sweeteners, we carried out our molecular dynamics studies using predicted structures of the TAS1R2/1R3 heterodimer and predicted structures for the sweeteners. We used as a measure of activation the intracellular ionic lock distance between transmembrane helices 3 and 6 of TAS1R3. We find that full synergic combinations [rebaudioside A (Reb-A)/acesulfame K and Reb-A/sucralose] and partial synergic combinations (sucralose/acesulfame K) show significantly more negative changes in the free energy compared to single-ligand cases, while a pair known to be suppressive (saccharin and acesulfame K) shows significantly less changes than for the single-ligand case. This study provides an atomistic understanding of the mechanism for synergy and identifies new combinations of sweeteners to reduce the caloric content for treating diseases.

Predicted Structure of Fully Activated Tas1R3/1R3' Homodimer Bound to G Protein and Natural Sugars: Structural Insights into G Protein Activation by a Class C Sweet Taste Homodimer with Natural Sugars

The Tas1R3 G protein-coupled receptor constitutes the main component of sweet taste sensory response in humans via forming a heterodimer with Tas1R2 or a homodimer with Tas1R3. The Tas1R3/1R3' homodimer serves as a low-affinity sweet taste receptor, stimulating gustducin G protein (GGust) signaling in the presence of a high concentration of natural sugars. This provides an additional means to detect the taste of natural sugars, thereby differentiating the flavors between natural sugars and artificial sweeteners. We report here the predicted 3D structure of active state Tas1R3/1R3' homodimer complexed with heterotrimeric GGust and sucrose. We discovered that the GGust makes ionic anchors to intracellular loops 1 and 2 of Tas1R3 while the Gα-α5 helix engages the cytoplasmic region extensively through salt bridge and hydrophobic interactions. We show that in the activation of this complex the Venus flytrap domains of the homodimer undergo a remarkable twist up to ∼100° rotation around the vertical axis to adopt a closed-closed conformation while the intracellular region relaxes to an open-open conformation. We find that binding of sucrose to the homodimer stabilizes a preactivated conformation with a largely open intracellular region that recruits and activates the GGust. Upon activation, the Gα subunit spontaneously opens up the nucleotide-binding site, making nucleotide exchange facile for signaling. This activation of GGust promotes the interdomain twist of the Venus flytrap domains. These structures and transformations could potentially be a basis for the design of new sweeteners with higher activity and less unpleasant flavors.

Structures and Agonist Binding Sites of Bitter Taste Receptor TAS2R5 Complexed with Gi Protein and Validated against Experiment

Bitter taste receptors (TAS2Rs) function in taste perception, but are also expressed in many extraoral tissues, presenting attractive therapeutic targets. TAS2R5s expressed on human airway smooth muscle cells can induce bronchodilation for treating asthma and other obstructive diseases. But TAS2R5s display low agonist affinity and the lack of a 3D structure has hindered efforts to design more active ligands. We report the structure of the activated TAS2R5 coupled to the Gi protein and bound to each of 19 agonists, using computational approaches. These agonists bind to two polar residues in TM3 that are unique for TAS2R5 among 25 TAS2R subtypes. Our predicted results correlate well with experimental results of agonist-receptor signaling coefficients, providing validation of the predicted structure. These results provide highly specific data on how agonists activate TAS2R5, how modifications of ligand structure alter receptor activation, and a guide to structure-based drug design.

Modeling of CCR5 Recognition by HIV-1 gp120: How the Viral Protein Exploits the Conformational Plasticity of the Coreceptor

The chemokine receptor CCR5 is a key player in HIV-1 infection. The cryo-EM 3D structure of HIV-1 envelope glycoprotein (Env) subunit gp120 in complex with CD4 and CCR5 has provided important structural insights into HIV-1/host cell interaction, yet it has not explained the signaling properties of Env nor the fact that CCR5 exists in distinct forms that show distinct Env binding properties. We used classical molecular dynamics and site-directed mutagenesis to characterize the CCR5 conformations stabilized by four gp120s, from laboratory-adapted and primary HIV-1 strains, and which were previously shown to bind differentially to distinct CCR5 forms and to exhibit distinct cellular tropisms. The comparative analysis of the simulated structures reveals that the different gp120s do indeed stabilize CCR5 in different conformational ensembles. They differentially reorient extracellular loops 2 and 3 of CCR5 and thus accessibility to the transmembrane binding cavity. They also reshape this cavity differently and give rise to different positions of intracellular ends of transmembrane helices 5, 6 and 7 of the receptor and of its third intracellular loop, which may in turn influence the G protein binding region differently. These results suggest that the binding of gp120s to CCR5 may have different functional outcomes, which could result in different properties for viruses.

Predicted structure of fully activated human bitter taste receptor TAS2R4 complexed with G protein and agonists

Bitter taste is sensed by bitter taste receptors (TAS2Rs) that belong to the G protein-coupled receptor (GPCR) superfamily. In addition to bitter taste perception, TAS2Rs have been reported recently to be expressed in many extraoral tissues and are now known to be involved in health and disease. Despite important roles of TAS2Rs in biological functions and diseases, no crystal structure is available to help understand the signal transduction mechanism or to help develop selective ligands as new therapeutic targets. We report here the three-dimensional structure of the fully activated TAS2R4 human bitter taste receptor predicted using the GEnSeMBLE complete sampling method. This TAS2R4 structure is coupled to the gustducin G protein and to each of several agonists. We find that the G protein couples to TAS2R4 by forming strong salt bridges to each of the three intracellular loops, orienting the activated Gα5 helix of the Gα subunit to interact extensively with the cytoplasmic region of the activated receptor. We find that the TAS2Rs exhibit unique motifs distinct from typical Class A GPCRs, leading to a distinct activation mechanism and a less stable inactive state. This fully activated bitter taste receptor complex structure provides insight into the signal transduction mechanism and into ligand binding to TAS2Rs.

Design, Synthesis, and Biological Evaluation of a Series of 5- and 7-Hydroxycoumarin Derivatives as 5-HT1A Serotonin Receptor Antagonists

We have designed and synthesized a series of 60 new 5- and 7-hydroxycoumarin derivatives bearing the piperazine moiety with the expected binding to 5-HT_1A and 5-HT_2A receptors. Molecular docking of all investigated compounds revealed subnanomolar estimates of 5-HT_1AR K_i for three ligands and 5-HT_2AR Ki for one ligand as well as numerous low nanomolar estimates of K_i for both receptors. Intrigued by these results we synthesized all 60 new derivatives using microwave-assisted protocols. We show that three new compounds show a relatively high antagonistic activity against the 5HT_1A receptor, although lower than the reference compound WAY-100635. These compounds also showed relatively low binding affinities to the 5-HT_2A receptor. We also provide a detailed structure-activity analysis of this series of compounds and compare it with previously obtained results for an exhaustive series of coumarin derivatives.

Non-full-length Water-Soluble CXCR4QTY and CCR5QTY Chemokine Receptors: Implication for Overlooked Truncated but Functional Membrane Receptors

It was posited that functionalities of GPCRs require full-length sequences that are negated by residue deletions. Here we report that significantly truncated nfCCR5^QTY and nfCXCR4^QTY still bind native ligands. Receptor-ligand interactions were discovered from yeast 2-hybrid screening and confirmed by mating selection. Two nfCCR5^QTY (SZ218a, SZ190b) and two nfCXCR4^QTY (SZ158a, SZ146a) were expressed in E. coli. Synthesized receptors exhibited α-helical structures and bound respective ligands with reduced affinities. SZ190b and SZ158a were reconverted into non-QTY forms and expressed in HEK293T cells. Reconverted receptors localized on cell membranes and functioned as negative regulators for ligand-induced signaling when co-expressed with full-length receptors. CCR5-SZ190b individually can perform signaling at a reduced level with higher ligand concentration. Our findings provide insight into essential structural components for CCR5 and CXCR4 functionality, while raising the possibility that non-full-length receptors may be resulted from alternative splicing and that pseudo-genes in genomes may be present and functional in living organisms.

•

Y2H screening reveals ligand interaction from truncated CXCR4 and CCR5 in QTY form

•

Truncated CCR5^QTY and CXCR4^QTY can be produced in E. coli and bind native ligands

•

Reconverted receptors localize on membranes and regulate cell signaling in HEK293

•

Our finding indicates potential presence and function for truncated receptors

Understanding G Protein Selectivity of Muscarinic Acetylcholine Receptors Using Computational Methods

The neurotransmitter molecule acetylcholine is capable of activating five muscarinic acetylcholine receptors, M1 through M5, which belong to the superfamily of G-protein-coupled receptors (GPCRs). These five receptors share high sequence and structure homology; however, the M1, M3, and M5 receptor subtypes signal preferentially through the Gαq/11 subset of G proteins, whereas the M2 and M4 receptor subtypes signal through the Gαi/o subset of G proteins, resulting in very different intracellular signaling cascades and physiological effects. The structural basis for this innate ability of the M1/M3/M5 set of receptors and the highly homologous M2/M4 set of receptors to couple to different G proteins is poorly understood. In this study, we used molecular dynamics (MD) simulations coupled with thermodynamic analyses of M1 and M2 receptors coupled to both Gαi and Gαq to understand the structural basis of the M1 receptor’s preference for the Gαq protein and the M2 receptor’s preference for the Gαi protein. The MD studies showed that the M1 and M2 receptors can couple to both Gα proteins such that the M1 receptor engages with the two Gα proteins in slightly different orientations and the M2 receptor engages with the two Gα proteins in the same orientation. Thermodynamic studies of the free energy of binding of the receptors to the Gα proteins showed that the M1 and M2 receptors bind more strongly to their cognate Gα proteins compared to their non-cognate ones, which is in line with previous experimental studies on the M3 receptor. A detailed analysis of receptor–G protein interactions showed some cognate-complex-specific interactions for the M2:Gαi complex; however, G protein selectivity determinants are spread over a large overlapping subset of residues. Conserved interaction between transmembrane helices 5 and 6 far away from the G-protein-binding receptor interface was found only in the two cognate complexes and not in the non-cognate complexes. An analysis of residues implicated previously in G protein selectivity, in light of the cognate and non-cognate structures, shaded a more nuanced role of those residues in affecting G protein selectivity. The simulation of both cognate and non-cognate receptor–G protein complexes fills a structural gap due to difficulties in determining non-cognate complex structures and provides an enhanced framework to probe the mechanisms of G protein selectivity exhibited by most GPCRs.

Molecular binding mode of PF-232798, a clinical anti-HIV candidate, at chemokine receptor CCR5

The chemokine receptor CCR5 is an important anti-HIV (human immunodeficiency virus) drug target owning to its pivotal role in HIV-1 viral entry as a co-receptor. Here, we present a 2.9 Å resolution crystal structure of CCR5 bound to PF-232798, a second-generation oral CCR5 antagonist currently in phase II clinical trials. PF-232798 and the marketed HIV drug maraviroc share a similar tropane scaffold with different amino (N)- and carboxyl (C)- substituents. Comparison of the CCR5-PF-232798 structure with the previously determined structure of CCR5 in complex with maraviroc reveals different binding modes of the two allosteric antagonists and subsequent conformational changes of the receptor. Our results not only offer insights into the phenomenon that PF-232798 has higher affinity and alternative resistance profile to maraviroc, but also will facilitate the design of new anti-HIV drugs.

CCR5 structural plasticity shapes HIV-1 phenotypic properties

CCR5 plays immune functions and is the coreceptor for R5 HIV-1 strains. It exists in diverse conformations and oligomerization states. We interrogated the significance of the CCR5 structural diversity on HIV-1 infection. We show that envelope glycoproteins (gp120s) from different HIV-1 strains exhibit divergent binding levels to CCR5 on cell lines and primary cells, but not to CD4 or the CD4i monoclonal antibody E51. This owed to differential binding of the gp120s to different CCR5 populations, which exist in varying quantities at the cell surface and are differentially expressed between different cell types. Some, but not all, of these populations are antigenically distinct conformations of the coreceptor. The different binding levels of gp120s also correspond to differences in their capacity to bind CCR5 dimers/oligomers. Mutating the CCR5 dimerization interface changed conformation of the CCR5 homodimers and modulated differentially the binding of distinct gp120s. Env-pseudotyped viruses also use particular CCR5 conformations for entry, which may differ between different viruses and represent a subset of those binding gp120s. In particular, even if gp120s can bind both CCR5 monomers and oligomers, impairment of CCR5 oligomerization improved viral entry, suggesting that HIV-1 prefers monomers for entry. From a functional standpoint, we illustrate that the nature of the CCR5 molecules to which gp120/HIV-1 binds shapes sensitivity to inhibition by CCR5 ligands and cellular tropism. Differences exist in the CCR5 populations between T-cells and macrophages, and this is associated with differential capacity to bind gp120s and to support viral entry. In macrophages, CCR5 structural plasticity is critical for entry of blood-derived R5 isolates, which, in contrast to prototypical M-tropic strains from brain tissues, cannot benefit from enhanced affinity for CD4. Collectively, our results support a role for CCR5 heterogeneity in diversifying the phenotypic properties of HIV-1 isolates and provide new clues for development of CCR5-targeting drugs.

Modulation of allosteric coupling by mutations: from protein dynamics and packing to altered native ensembles and function

A large body of work has gone into understanding the effect of mutations on protein structure and function. Conventional treatments have involved quantifying the change in stability, activity and relaxation rates of the mutants with respect to the wild-type protein. However, it is now becoming increasingly apparent that mutational perturbations consistently modulate the packing and dynamics of a significant fraction of protein residues, even those that are located >10–15 Å from the mutated site. Such long-range modulation of protein features can distinctly tune protein stability and the native conformational ensemble contributing to allosteric modulation of function. In this review, I summarize a series of experimental and computational observations that highlight the incredibly pliable nature of proteins and their response to mutational perturbations manifested via the intra-protein interaction network. I highlight how an intimate understanding of mutational effects could pave the way for integrating stability, folding, cooperativity and even allostery within a single physical framework.

CCR5 Revisited: How Mechanisms of HIV Entry Govern AIDS Pathogenesis

The chemokine receptor CCR5 has been the focus of intensive studies since its role as a coreceptor for HIV entry was discovered in 1996. These studies lead to the development of small molecular drugs targeting CCR5, with maraviroc becoming in 2007 the first clinically approved chemokine receptor inhibitor. More recently, the apparent HIV cure in a patient transplanted with hematopoietic stem cells devoid of functional CCR5 rekindled the interest for inactivating CCR5 through gene therapy and pharmacological approaches. Fundamental research on CCR5 has also been boosted by key advances in the field of G-protein coupled receptor research, with the realization that CCR5 adopts a variety of conformations, and that only a subset of these conformations may be targeted by chemokine ligands. In addition, recent genetic and pathogenesis studies have emphasized the central role of CCR5 expression levels in determining the risk of HIV and SIV acquisition and disease progression. In this article, we propose to review the key properties of CCR5 that account for its central role in HIV pathogenesis, with a focus on mechanisms that regulate CCR5 expression, conformation, and interaction with HIV envelope glycoproteins.

The 3D Structure of Human DP Prostaglandin G-Protein-Coupled Receptor Bound to Cyclopentanoindole Antagonist, Predicted Using the DuplexBiHelix Modification of the GEnSeMBLE Method

Prostaglandins play a critical physiological role in both cardiovascular and immune systems, acting through their interactions with 9 prostanoid G protein-coupled receptors (GPCRs). These receptors are important therapeutic targets for a variety of diseases including arthritis, allergies, type 2 diabetes, and cancer. The DP prostaglandin receptor is of interest because it has unique structural and physiological properties. Most notably, DP does not have the 3-6 ionic lock common to Class A GPCRs. However, the lack of X-ray structures for any of the 9 prostaglandin GPCRs hampers the application of structure-based drug design methods to develop more selective and active medications to specific receptors. We predict here 3D structures for the DP prostaglandin GPCR, based on the GEnSeMBLE complete sampling with hierarchical scoring (CS-HS) methodology. This involves evaluating the energy of 13 trillion packings to finally select the best 20 that are stable enough to be relevant for binding to antagonists, agonists, and modulators. To validate the predicted structures, we predict the binding site for the Merck cyclopentanoindole (CPI) selective antagonist docked to DP. We find that the CPI binds vertically in the 1-2-7 binding pocket, interacting favorably with residues R310^7.40 and K76^2.54 with additional interactions with S313^7.43, S316^7.46, S19^1.35, etc. This binding site differs significantly from that of antagonists to known Class A GPCRs where the ligand binds in the 3-4-5-6 region. We find that the predicted binding site leads to reasonable agreement with experimental Structure-Activity Relationship (SAR). We suggest additional mutation experiments including K76^2.54, E129^3.49, L123^3.43, M270^6.40, F274^6.44 to further validate the structure, function, and activation mechanism of receptors in the prostaglandin family. Our structures and binding sites are largely consistent and improve upon the predictions by Li et al. ( J. Am. Chem. Soc. 2007 , 129 ( 35 ), 10720 ) that used our earlier MembStruk prediction methodology.

Conformational and Thermodynamic Landscape of GPCR Activation from Theory and Computation

We present a hybrid computational methodology to predict multiple energetically accessible conformations for G protein-coupled receptors (GPCRs) that might play a role in binding to ligands and different signaling partners. To our knowledge, this method, termed ActiveGEnSeMBLE, enables the first quantitative energy profile for GPCR activation that is consistent with the qualitative profile deduced from experiments. ActiveGEnSeMBLE starts with a systematic coarse grid sampling of helix tilts/rotations (∼13 trillion transmembrane-domain conformations) and selects the conformational landscape based on energy. This profile identifies multiple potential active-state energy wells, with the TM3-TM6 intracellular distance as an approximate activation coordinate. These energy wells are then sampled locally using a finer grid to find locally minimized conformation in each energy well. We validate this strategy using the inactive and active experimental structures of β2 adrenergic receptor (hβ2AR) and M2 muscarinic acetylcholine receptor. Structures of membrane-embedded hβ2AR along its activation coordinate are subjected to molecular-dynamics simulations for relaxation and interaction energy analysis to generate a quantitative energy landscape for hβ2AR activation. This landscape reveals several metastable states along this coordinate, indicating that for hβ2AR, the agonist alone is not enough to stabilize the active state and that the G protein is necessary, consistent with experimental observations. The method's application to somatostatin receptor SSTR5 (no experimental structure available) shows that to predict an active conformation it is better to start from an inactive structure template based on a close homolog than to start from an active template based on a distant homolog. The energy landscape for hSSTR5 activation is consistent with hβ2AR in the role of the G protein. These results demonstrate the utility of the ActiveGEnSeMBLE method for predicting multiple conformations along the pathways for activating GPCRs and the corresponding energy landscapes, thereby providing detailed structural insights into the initial molecular events of GPCR function that are not easily accessible by experiments.

Activation mechanism of the G protein-coupled sweet receptor heterodimer with sweeteners and allosteric agonists

The sweet taste in humans is mediated by the TAS1R2/TAS1R3 G protein-coupled receptor (GPCR), which belongs to the class C family that also includes the metabotropic glutamate and γ-aminobutyric acid receptors. We report here the predicted 3D structure of the full-length TAS1R2/TAS1R3 heterodimer, including the Venus Flytrap Domains (VFDs) [in the closed-open (co) active conformation], the cysteine-rich domains (CRDs), and the transmembrane domains (TMDs) at the TM56/TM56 interface. We observe that binding of agonists to VFD2 of TAS1R2 leads to major conformational changes to form a TM6/TM6 interface between TMDs of TAS1R2 and TAS1R3, which is consistent with the activation process observed biophysically on the metabotropic glutamate receptor 2 homodimer. We find that the initial effect of the agonist is to pull the bottom part of VFD3/TAS1R3 toward the bottom part of VFD2/TAS1R2 by ∼6 Å and that these changes get transmitted from VFD2 of TAS1R2 (where agonists bind) through the VFD3 and the CRD3 to the TMD3 of TAS1R3 (which couples to the G protein). These structural transformations provide a detailed atomistic mechanism for the activation process in GPCR, providing insights and structural details that can now be validated through mutation experiments.

Molecular-Docking-Based Drug Design and Discovery

Currently 30-50% of drug targets are G Protein-Coupled Receptors (GPCRs). However, the clinical useful drugs for targeting GPCR have been limited by the lack of subtype selectivity or efficacy, leading to undesirable side effects. To develop subtype-selective GPCR ligands with desired molecular properties, better understanding is needed of the pharmacophore elements and of the binding mechanism required for subtype selectivity. To illustrate these issues, we describe here three successful applications to understand the binding mechanism associated with subtype selectivity: 5-HT2B (5-Hydroxytryptamine, 5-HT) serotonin receptor (HT2BR), H3 histamine receptor (H3HR) and A3 adenosine receptor (A3AR). The understanding of structure-function relationships among individual types and subtypes of GPCRs gained from such computational predictions combined with experimental validation and testing is expected the development of new highly selective and effective ligands to address such diseases while minimizing side-effects.

mutLBSgeneDB: mutated ligand binding site gene DataBase

Mutations at the ligand binding sites (LBSs) can influence protein structure stability, binding affinity with small molecules, and drug resistance in cancer patients. Our recent analysis revealed that ligand binding residues had a significantly higher mutation rate than other parts of the protein. Here, we built mutLBSgeneDB (mutated Ligand Binding Site gene DataBase) available at http://zhaobioinfo.org/mutLBSgeneDB. We collected and curated over 2300 genes (mutLBSgenes) having ∼12 000 somatic mutations at ∼10 000 LBSs across 16 cancer types and selected 744 drug targetable genes (targetable_mutLBSgenes) by incorporating kinases, transcription factors, pharmacological genes, and cancer driver genes. We analyzed LBS mutation information, differential gene expression network, drug response correlation with gene expression, and protein stability changes for all mutLBSgenes using integrated genetic, genomic, transcriptomic, proteomic, network and functional information. We calculated and compared the binding affinities of 20 carefully selected genes with their drugs in wild type and mutant forms. mutLBSgeneDB provides a user-friendly web interface for searching and browsing through seven categories of annotations: Gene summary, Mutated information, Protein structure related information, Differential gene expression and gene-gene network, Phenotype information, Pharmacological information, and Conservation information. mutLBSgeneDB provides a useful resource for functional genomics, protein structure, drug and disease research communities.

Short Communication: Limited Anti-HIV-1 Activity of Maraviroc in Mucosal Tissues

The potential of maraviroc (MVC), a small-molecule CCR5 antagonist, as a candidate to prevent HIV-1 sexual transmission by oral or topical dosing has not yet been completely established. Using relevant cellular and mucosal tissue explant models, we show partial antiviral activity of MVC when tested in multiple preclinical dosing strategies.

Structure-Based Sequence Alignment of the Transmembrane Domains of All Human GPCRs: Phylogenetic, Structural and Functional Implications

The understanding of G-protein coupled receptors (GPCRs) is undergoing a revolution due to increased information about their signaling and the experimental determination of structures for more than 25 receptors. The availability of at least one receptor structure for each of the GPCR classes, well separated in sequence space, enables an integrated superfamily-wide analysis to identify signatures involving the role of conserved residues, conserved contacts, and downstream signaling in the context of receptor structures. In this study, we align the transmembrane (TM) domains of all experimental GPCR structures to maximize the conserved inter-helical contacts. The resulting superfamily-wide GpcR Sequence-Structure (GRoSS) alignment of the TM domains for all human GPCR sequences is sufficient to generate a phylogenetic tree that correctly distinguishes all different GPCR classes, suggesting that the class-level differences in the GPCR superfamily are encoded at least partly in the TM domains. The inter-helical contacts conserved across all GPCR classes describe the evolutionarily conserved GPCR structural fold. The corresponding structural alignment of the inactive and active conformations, available for a few GPCRs, identifies activation hot-spot residues in the TM domains that get rewired upon activation. Many GPCR mutations, known to alter receptor signaling and cause disease, are located at these conserved contact and activation hot-spot residue positions. The GRoSS alignment places the chemosensory receptor subfamilies for bitter taste (TAS2R) and pheromones (Vomeronasal, VN1R) in the rhodopsin family, known to contain the chemosensory olfactory receptor subfamily. The GRoSS alignment also enables the quantification of the structural variability in the TM regions of experimental structures, useful for homology modeling and structure prediction of receptors. Furthermore, this alignment identifies structurally and functionally important residues in all human GPCRs. These residues can be used to make testable hypotheses about the structural basis of receptor function and about the molecular basis of disease-associated single nucleotide polymorphisms.

G-protein coupled receptors (GPCRs) are a large superfamily of integral membrane proteins that share a characteristic 7 transmembrane helix fold. They detect various molecules outside of the cell and signal their presence to the inside of the cell. At least half of the 800 human GPCRs are potential drug targets, so understanding their structure and function is critical. Experimental structures are now available for at least one receptor from each GPCR class. The structure of the 7 helix fold is highly conserved even for receptors with very low sequence similarity. We analyze the available experimental structures and compare the common inter-helical contacts. Our analysis leads to a unified sequence-structure alignment of the GPCR superfamily that can then be used as the starting point for structure prediction of all other GPCRs. A key result of our analysis is a list of conserved contact residues and activation “hot-spots” residues that are critical for GPCR folding and function. We propose that mutations and natural variants of amino acids at these locations in the GPCRs can dramatically influence their activation state and alter intracellular signaling. This provides hypotheses for the molecular mechanisms underlying disease causing mutants for any GPCR.

Computational Prediction and Biochemical Analyses of New Inverse Agonists for the CB1 Receptor

Human cannabinoid type 1 (CB1) G-protein coupled receptor is a potential therapeutic target for obesity. The previously predicted and experimentally validated ensemble of ligand-free conformations of CB1 [Scott, C. E. et al. Protein Sci. 2013 , 22 , 101 - 113 ; Ahn, K. H. et al. Proteins 2013 , 81 , 1304 - 1317] are used here to predict the binding sites for known CB1-selective inverse agonists including rimonabant and its seven known derivatives. This binding pocket, which differs significantly from previously published models, is used to identify 16 novel compounds expected to be CB1 inverse agonists by exploiting potential new interactions. We show experimentally that two of these compounds exhibit inverse agonist properties including inhibition of basal and agonist-induced G-protein coupling activity, as well as an enhanced level of CB1 cell surface localization. This demonstrates the utility of using the predicted binding sites for an ensemble of CB1 receptor structures for designing new CB1 inverse agonists.

Molecular level understanding of resistance to nalidixic acid in Salmonella enteric serovar typhimurium associates with the S83F sequence type

Nalidixic acid is an antibiotic drug used for treatment of Salmonellosis, a gastrointestinal infection. DNA gyrase subunit A (GyrA) of Salmonella typhimurium is the drug target for nalidixic acid. Resistance of GyrA to nalidixic acid, because of a point mutation in S. typhimurium, was recently reported. Substitution of Phe in place of Ser at locus 83 in GyrA of S. typhimurium has been experimentally associated with nalidixic acid resistance. Despite recent efforts, the mechanism of this resistance is not well understood. In this investigation we used computational techniques to address this shortcoming. Our results showed that contact with residue Arg 91 is certainly important for efficient binding of nalidixic acid to the target protein, and that mutation of this residue results in 180° rotation of the antibiotic in its binding pocket, around its own long axis. It is hoped these findings may enable development of new antibiotics against resistant forms of Salmonella.

A single-residue change in the HIV-1 V3 loop associated with maraviroc resistance impairs CCR5 binding affinity while increasing replicative capacity

Background

Maraviroc (MVC) is an allosteric CCR5 inhibitor used against HIV-1 infection. While MVC-resistant viruses have been identified in patients, it still remains incompletely known how they adjust their CD4 and CCR5 binding properties to resist MVC inhibition while preserving their replicative capacity. It is thought that they maintain high efficiency of receptor binding. To date however, information about the binding affinities to receptors for inhibitor-resistant HIV-1 remains limited.

Results

Here, we show by means of viral envelope (gp120) binding experiments and virus-cell fusion kinetics that a MVC-resistant virus (MVC-Res) that had emerged as a dominant viral quasispecies in a patient displays reduced affinities for CD4 and CCR5 either free or bound to MVC, as compared to its MVC-sensitive counterpart isolated before MVC therapy. An alanine insertion within the GPG motif (G310_P311insA) of the MVC-resistant gp120 V3 loop is responsible for the decreased CCR5 binding affinity, while impaired binding to CD4 is due to sequence changes outside V3. Molecular dynamics simulations of gp120 binding to CCR5 further emphasize that the Ala insertion alters the structure of the V3 tip and weakens interaction with CCR5 ECL2. Paradoxically, infection experiments on cells expressing high levels of CCR5 also showed that Ala allows MVC-Res to use CCR5 efficiently, thereby improving viral fusion and replication efficiencies. Actually, although we found that the V3 loop of MVC-Res is required for high levels of MVC resistance, other regions outside V3 are sufficient to confer a moderate level of resistance. These sequence changes outside V3, however, come with a replication cost, which is compensated for by the Ala insertion in V3.

Conclusion

These results indicate that changes in the V3 loop of MVC-resistant viruses can augment the efficiency of CCR5-dependent steps of viral entry other than gp120 binding, thereby compensating for their decreased affinity for entry receptors and improving their fusion and replication efficiencies. This study thus sheds light on unsuspected mechanisms whereby MVC-resistant HIV-1 could emerge and grow in treated patients.

Electronic supplementary material

The online version of this article (doi:10.1186/s12977-015-0177-1) contains supplementary material, which is available to authorized users.

CCR5 susceptibility to ligand-mediated down-modulation differs between human T lymphocytes and myeloid cells

CCR5 is a chemokine receptor expressed on leukocytes and a coreceptor used by HIV-1 to enter CD4(+) T lymphocytes and macrophages. Stimulation of CCR5 by chemokines triggers internalization of chemokine-bound CCR5 molecules in a process called down-modulation, which contributes to the anti-HIV activity of chemokines. Recent studies have shown that CCR5 conformational heterogeneity influences chemokine-CCR5 interactions and HIV-1 entry in transfected cells or activated CD4(+) T lymphocytes. However, the effect of CCR5 conformations on other cell types and on the process of down-modulation remains unclear. We used mAbs, some already shown to detect distinct CCR5 conformations, to compare the behavior of CCR5 on in vitro generated human T cell blasts, monocytes and MDMs and CHO-CCR5 transfectants. All human cells express distinct antigenic forms of CCR5 not detected on CHO-CCR5 cells. The recognizable populations of CCR5 receptors exhibit different patterns of down-modulation on T lymphocytes compared with myeloid cells. On T cell blasts, CCR5 is recognized by all antibodies and undergoes rapid chemokine-mediated internalization, whereas on monocytes and MDMs, a pool of CCR5 molecules is recognized by a subset of antibodies and is not removed from the cell surface. We demonstrate that this cell surface-retained form of CCR5 responds to prolonged treatment with more-potent chemokine analogs and acts as an HIV-1 coreceptor. Our findings indicate that the regulation of CCR5 is highly specific to cell type and provide a potential explanation for the observation that native chemokines are less-effective HIV-entry inhibitors on macrophages compared with T lymphocytes.

The predicted ensemble of low-energy conformations of human somatostatin receptor subtype 5 and the binding of antagonists

Human somatostatin receptor subtype 5 (hSSTR5) regulates cell proliferation and hormone secretion. However, the identification of effective therapeutic small-molecule ligands is impeded because experimental structures are not available for any SSTR subtypes. Here, we predict the ensemble of low-energy 3D structures of hSSTR5 using a modified GPCR Ensemble of Structures in Membrane BiLayer Environment (GEnSeMBLE) complete sampling computational method. We find that this conformational ensemble displays most interhelical interactions conserved in class A G protein-coupled receptors (GPCRs) plus seven additional interactions (e.g., Y2.43-D3.49, T3.38-S4.53, K5.64-Y3.51) likely conserved among SSTRs. We then predicted the binding sites for a series of five known antagonists, leading to predicted binding energies consistent with experimental results reported in the literature. Molecular dynamics (MD) simulation of 50 ns in explicit water and lipid retained the predicted ligand-bound structure and formed new interaction patterns (e.g. R3.50-T6.34) consistent with the inactive μ-opioid receptor X-ray structure. We suggest more than six mutations for experimental validation of our prediction. The final predicted receptor conformations and antagonist binding sites provide valuable insights for designing new small-molecule drugs targeting SSTRs.

Predicted structures for kappa opioid G-protein coupled receptor bound to selective agonists

Human kappa opioid receptor (κ-OR), a G protein-coupled receptor (GPCR), has been identified as a drug target for treatment of such human disorders as pain perception, neuroendocrine physiology, affective behavior, and cognition. In order to find more selective and active agonists, one would like to do structure based drug design. Indeed, there is an X-ray structure for an antagonist bound to κ-OR, but structures for activated GPCRs are quite different from those for the inactive GPCRs. Here we predict the ensemble of 24 low-energy structures of human kappa opioid receptor (κ-OR), obtained by application of the GEnSeMBLE (GPCR Ensemble of Structures in Membrane Bilayer Environment) complete sampling method, which evaluates 13 trillion combinations of tilt and rotation angles for κ-OR to select the best 24. To validate these structures, we used the DarwinDock complete sampling method to predict the binding sites for five known agonists (ethylketocyclazocine, bremazocine, pentazocine, nalorphine, and morphine) bound to all 24 κ-OR conformations. We find that some agonists bind selectively to receptor conformations that lack the salt bridge between transmembrane domains 3 and 6 as expected for active conformations. These 3D structures for κ-OR provide a structural basis for understanding ligand binding and activation of κ-OR, which should be useful for guiding subtype specific drug design.