Structure and dynamics of G-protein coupled receptors

Genetic variation in the chemokine receptor 5 gene and course of HIV infection; review on genetics and immunological aspect

Chemokines are small protein molecules associated with various physiological events precisely in immune modulation via chemokine receptors. The chemokine receptors are G-protein coupled receptors express mainly on the cell surface of immune cells. Retroviruses, including HIV in the early stage of infection, primarily target chemokines receptors and get internalized easily into immune cells; T cell and escape from immune surveillance. HIV glycoprotein selectively develops an affinity for the extracellular domain of chemokines receptors and allows the pathogen to internalize via CCR-5. Now, CCR-5 remains a crucial signaling pathway that can be translated into the therapeutic target by changing the receptor protein environment. Many populations have a mutation in coding and promoter regions of CCR-5, tuning a resistance for HIV infection. Natively, there are several mechanisms where the human genome remains in the dynamic state by changing its composition and acquiring variations. Single nucleotide polymorphism is spontaneous phenomenon responsible for precise and point mutation at the genome. Several studies have demonstrated that European and African American populations are enriched in significant CCR5 promoter SNP (CCR5Δ32) in the coding and promoter region as well. Now, such SNP can be an early-stage biomarker in studying HIV and other similar infections. Here, in this study, we have elucidated the role of SNP (both the promoter and coding region) and the fate of HIV infections. We also empathized with the genetics of such SNPs, mostly frequency and its immunological impact.

Characterizing Interhelical Interactions of G-Protein Coupled Receptors with the Fragment Molecular Orbital Method

G-protein coupled receptors (GPCRs) are the largest superfamily of membrane proteins, regulating almost every aspect of cellular activity and serving as key targets for drug discovery. We have identified an accurate and reliable computational method to characterize the strength and chemical nature of the interhelical interactions between the residues of transmembrane (TM) domains during different receptor activation states, something that cannot be characterized solely by visual inspection of structural information. Using the fragment molecular orbital (FMO) quantum mechanics method to analyze 35 crystal structures representing different branches of the class A GPCR family, we have identified 69 topologically equivalent TM residues that form a consensus network of 51 inter-TM interactions, providing novel results that are consistent with and help to rationalize experimental data. This discovery establishes a comprehensive picture of how defined molecular forces govern specific interhelical interactions which, in turn, support the structural stability, ligand binding, and activation of GPCRs.

Binding Interactions of Ergotamine and Dihydroergotamine to 5-Hydroxytryptamine Receptor 1B (5-HT1b) Using Molecular Dynamics Simulations and Dynamic Network Analysis

Ergotamine (ERG) and dihydroergotamine (DHE), common migraine drugs, have small structural differences but lead to clinically important distinctions in their pharmacological profiles. For example, DHE is less potent than ERG by about 10-fold at the 5-hydroxytrptamine receptor 1B (5-HT_1B). Although the high-resolution crystal structures of the 5-HT_1B receptor with both ligands have been solved, the high similarity between these two complex structures does not sufficiently explain their activity differences and the activation mechanism of the receptor. Hence, an examination of the dynamic motion of both drugs with the receptor is required. In this study, we ran a total of 6.0 μs molecular dynamics simulations on each system. Our simulation data show the subtle variations between the two systems in terms of the ligand-receptor interactions and receptor secondary structures. More importantly, the ligand and protein root-mean-square fluctuations (RMSFs) for the two systems were distinct, with ERG having a trend of lower RMSF values, indicating it to be bound tighter to 5-HT_1B with less fluctuations. The molecular mechanism-general born surface area (MM-GBSA) binding energies illustrate this further, proving ERG has an overall stronger MM-GBSA binding energy. Analysis of several different microswitches has shown that the 5-HT_1B-ERG complex is in a more active conformation state than 5-HT_1B-DHE, which is further supported by the dynamic network model, with reference to mutagenesis data with the critical nodes and the first three low-energy modes from the normal mode analysis. We also identify Trp327^6.48 and Phe331^6.52 as key residues involved in the active state 5-HT_1B for both ligands. Using the detailed dynamic information from our analysis, we made predictions for possible modifications to DHE and ERG that yielded five derivatives that might have more favorable binding energies and reduced structural fluctuations.

Identification of functional divergence sites in dopamine receptors of vertebrates

Dopamine is one of the major neurotransmitters in the brain and body, and regulates a wide variety of functions via its binding with dopamine receptors. Abnormalities in dopamine receptors have also been found to be related to various neurological disorders. For such reason, dopamine receptors are among the key components to understanding the molecular mechanisms of many diseases, they are also the potential drug targets for the treatment of many diseases. Till now, five different dopamine receptors (D1-D5) have been identified in mammals, which are assumed to be evolved from a common ancestor after multiple gene duplication events and functional divergence. Thus, identifying the specific features of each dopamine receptor, will not only provide clues for understanding the functional differences between the receptors, but also help us to design drugs specific for a certain subtype of receptor. In this study, we investigated the functional divergence in dopamine receptors in representative vertebrate species by analyzing their molecular evolution features. Our results showed that the coefficients for type I functional divergence (θ_I) were significantly greater than 0 for all the pairwise comparisons between the five dopamine receptors, suggesting that type I functional divergence, i.e., altered functional constraints or different evolutionary rates, may have taken place at some amino acids in the receptors. We further identified 84 potential type I functional divergence peptide sites for the pairwise comparisons between the D1-like and D2-like are identified in total. When these sites were mapped to the 3D structure of dopamine receptors, most of them were included in ICL3, M6 and M7 domains. Especially, sixteen of these sites may be the major sites associated with the changes of properties between D1-like and D2-like receptors. These sites provide molecular basis for further studies such as dopamine receptor function exploration and subtype specific drug design and screening.

In Silico Workflow for the Discovery of Natural Products Activating the G Protein-Coupled Bile Acid Receptor 1

The G protein-coupled bile acid receptor (GPBAR1) has been recognized as a promising new target for the treatment of diverse diseases, including obesity, type 2 diabetes, fatty liver disease and atherosclerosis. The identification of novel and potent GPBAR1 agonists is highly relevant, as these diseases are on the rise and pharmacological unmet therapeutic needs are pervasive. Therefore, the aim of this study was to develop a proficient workflow for the in silico prediction of GPBAR1 activating compounds, primarily from natural sources. A protocol was set up, starting with a comprehensive collection of structural information of known ligands. This information was used to generate ligand-based pharmacophore models in LigandScout 4.08 Advanced. After theoretical validation, the two most promising models, namely BAMS22 and TTM8, were employed as queries for the virtual screening of natural product and synthetic small molecule databases. Virtual hits were progressed to shape matching experiments and physicochemical clustering. Out of 33 diverse virtual hits subjected to experimental testing using a reporter gene-based assay, two natural products, farnesiferol B (27) and microlobidene (28), were confirmed as GPBAR1 activators reaching more than 50% receptor activation at 20 μM with EC50s of 13.53 μM and 13.88 μM, respectively. This activity is comparable to that of the endogenous ligand lithocholic acid (1). Seven further virtual hits showed activity reaching at least 15% receptor activation either at 5 or 20 μM, including new scaffolds from natural and synthetic origin.

Distinctive binding properties of the negative allosteric modulator, [3H]SB269,652, at recombinant dopamine D3 receptors

Recently, employing radioligand displacement and functional coupling studies, we demonstrated that SB269,652 (N-[(1r,4r)-4-[2-(7-cyano-1,2,3,4-tetrahydroisoquinolin-2-yl)ethyl]cyclohexyl]-1H-indole-2-carboxamide) interacts in an atypical manner with dopamine D₃ receptor displaying a unique profile reminiscent of a negative allosteric ligand. Here, we characterized the binding of radiolabelled [³H]SB269,652 to human dopamine D₃ receptor stably expressed in Chinese Hamster Ovary cells. Under saturating conditions, SB269,652 showed a KD value of ≈ 1nM. Consistent with high selectivity for human dopamine D₃ receptor, [³H]SB269,652 binding was undetectable in cells expressing human dopamine D₁, D_2L or D₄ receptors and absent in synaptosomes from dopamine D₃ receptor knockout vs. wild-type mice. In contrast to saturation binding experiments, the dissociation kinetics of [³H]SB269,652 from human dopamine D₃ receptors initiated with an excess of unlabelled ligand were best fitted by a bi-exponential binding model. Supporting the kinetic data, competition experiments with haloperidol, S33084 (a dopamine D₃ receptor antagonist) or dopamine, were best described by a two-site model. In co-transfection experiments binding of SB269,652 to dopamine D₃ receptor was able to influence the functional coupling of dopamine D₂ receptor, supporting the notion that SB269,652 is a negative allosteric modulator across receptor dimers. However, because SB269,652 decreases the rate of [³H]nemonapride dissociation, the present data suggest that SB269,652 behaves as a bitopic antagonist at unoccupied dopamine D₃ receptor, binding simultaneously to both orthosteric and allosteric sites, and as a pure negative allosteric modulator when receptors are occupied and it can solely bind to the allosteric site.

Identification of a Different Agonist-Binding Site and Activation Mechanism of the Human P2Y1 Receptor

The human P2Y1 receptor (P2Y1R) is a purinergic G-protein-coupled receptor (GPCR) that functions as a receptor for adenosine 5'-diphosphate (ADP). An antagonist of P2Y1R might potentially have antithrombotic effects, whereas agonists might serve as antidiabetic agents. On the basis of the antagonist-bound MRS2500-P2Y1R crystal structure, we constructed computational models of apo-P2Y1R and the agonist-receptor complex 2MeSADP-P2Y1R. We then performed conventional molecular dynamics (cMD) and accelerated molecular dynamics (aMD) simulations to study the conformational dynamics after binding with agonist/antagonist as well as the P2Y1R activation mechanism. We identified a new agonist-binding site of P2Y1R that is consistent with previous mutagenesis data. This new site is deeper than those of the agonist ADP in the recently simulated ADP-P2Y1R structure and the antagonist MRS2500 in the MRS2500-P2Y1R crystal structure. During P2Y1R activation, the cytoplasmic end of helix VI shifts outward 9.1 Å, the Ser1463.47-Tyr2375.58 hydrogen bond breaks, a Tyr2375.58-Val2626.37 hydrogen bond forms, and the conformation of the χ1 rotamer of Phe2696.44 changes from parallel to perpendicular to helix VI. The apo-P2Y1R system and the MRS2500-P2Y1R system remain inactive. The newly identified agonist binding site and activation mechanism revealed in this study may aid in the design of P2Y1R antagonists/agonists as antithrombotic/antidiabetic agents, respectively.

CXCL10 Acts as a Bifunctional Antimicrobial Molecule against Bacillus anthracis

Bacillus anthracis is killed by the interferon-inducible, ELR(−) CXC chemokine CXCL10. Previous studies showed that disruption of the gene encoding FtsX, a conserved membrane component of the ATP-binding cassette transporter-like complex FtsE/X, resulted in resistance to CXCL10. FtsX exhibits some sequence similarity to the mammalian CXCL10 receptor, CXCR3, suggesting that the CXCL10 N-terminal region that interacts with CXCR3 may also interact with FtsX. A C-terminal truncated CXCL10 was tested to determine if the FtsX-dependent antimicrobial activity is associated with the CXCR3-interacting N terminus. The truncated CXCL10 exhibited antimicrobial activity against the B. anthracis parent strain but not the ΔftsX mutant, which supports a key role for the CXCL10 N terminus. Mutations in FtsE, the conserved ATP-binding protein of the FtsE/X complex, resulted in resistance to both CXCL10 and truncated CXCL10, indicating that both FtsX and FtsE are important. Higher concentrations of CXCL10 overcame the resistance of the ΔftsX mutant to CXCL10, suggesting an FtsX-independent killing mechanism, likely involving its C-terminal α-helix, which resembles a cationic antimicrobial peptide. Membrane depolarization studies revealed that CXCL10 disrupted membranes of the B. anthracis parent strain and the ΔftsX mutant, but only the parent strain underwent depolarization with truncated CXCL10. These findings suggest that CXCL10 is a bifunctional molecule that kills B. anthracis by two mechanisms. FtsE/X-dependent killing is mediated through an N-terminal portion of CXCL10 and is not reliant upon the C-terminal α-helix. The FtsE/X-independent mechanism involves membrane depolarization by CXCL10, likely because of its α-helix. These findings present a new paradigm for understanding mechanisms by which CXCL10 and related chemokines kill bacteria.

Chemokines are a class of molecules known for their chemoattractant properties but more recently have been shown to possess antimicrobial activity against a wide range of Gram-positive and Gram-negative bacterial pathogens. The mechanism(s) by which these chemokines kill bacteria is not well understood, but it is generally thought to be due to the conserved amphipathic C-terminal α-helix that resembles cationic antimicrobial peptides in charge and secondary structure. Our present study indicates that the interferon-inducible, ELR(−) chemokine CXCL10 kills the Gram-positive pathogen Bacillus anthracis through multiple molecular mechanisms. One mechanism is mediated by interaction of CXCL10 with the bacterial FtsE/X complex and does not require the presence of the CXCL10 C-terminal α-helix. The second mechanism is FtsE/X receptor independent and kills through membrane disruption due to the C-terminal α-helix. This study represents a new paradigm for understanding how chemokines exert an antimicrobial effect that may prove applicable to other bacterial species.

Commercially available antibodies directed against α-adrenergic receptor subtypes and other G protein-coupled receptors with acceptable selectivity in flow cytometry experiments

Several previous reports suggested that many commercially available antibodies directed against G protein-coupled receptors (GPCR) lack sufficient selectivity. Accordingly, it has been proposed that receptor antibodies should be validated by at least one of several criteria, such as testing tissues or cells after knockout or silencing of the corresponding gene. Here, we tested whether 12 commercially available antibodies directed against α-adrenergic receptor (AR) subtypes (α1A/B/D, α2A/B/C), atypical chemokine receptor 3 (ACKR3), and vasopressin receptor 1A (AVPR1A) suffice these criteria. We detected in flow cytometry experiments with human vascular smooth muscle cells that the fluorescence signals from each of these antibodies were reduced by 46 ± 10 %-91 ± 2 % in cells treated with commercially available small interfering RNA (siRNA) specific for each receptor, as compared with cells that were incubated with non-targeting siRNA. The tested antibodies included anti-ACKR3 (R&D Systems, mab42273), for which specificity has previously been demonstrated. Staining with this antibody resulted in 72 ± 5 % reduction of the fluorescence signal after ACKR3 siRNA treatment. Furthermore, staining with anti-α1A-AR (Santa Cruz, sc1477) and anti-ACKR3 (Abcam, ab38089), which have previously been reported to be non-specific, resulted in 70 ± 19 % and 80 ± 4 % loss of the fluorescence signal after α1A-AR and ACKR3 siRNA treatment, respectively. Our findings demonstrate that the tested antibodies show reasonable selectivity for their receptor target under our experimental conditions. Furthermore, our observations suggest that the selectivity of GPCR antibodies depends on the method for which the antibody is employed, the species from which cells/tissues are obtained, and on the type of specimens (cell, tissue/cell homogenate, or section) tested.

Towards an understanding of the structural basis for insect olfaction by odorant receptors

Insects have co-opted a unique family of seven transmembrane proteins for odour sensing. Odorant receptors are believed to have evolved from gustatory receptors somewhere at the base of the Hexapoda and have expanded substantially to become the dominant class of odour recognition elements within the Insecta. These odorant receptors comprise an obligate co-receptor, Orco, and one of a family of highly divergent odorant "tuning" receptors. The two subunits are thought to come together at some as-yet unknown stoichiometry to form a functional complex that is capable of both ionotropic and metabotropic signalling. While there are still no 3D structures for these proteins, site-directed mutagenesis, resonance energy transfer, and structural modelling efforts, all mainly on Drosophila odorant receptors, are beginning to inform hypotheses of their structures and how such complexes function in odour detection. Some of the loops, especially the second extracellular loop that has been suggested to form a lid over the binding pocket, and the extracellular regions of some transmembrane helices, especially the third and to a less extent the sixth and seventh, have been implicated in ligand recognition in tuning receptors. The possible interaction between Orco and tuning receptor subunits through the final intracellular loop and the adjacent transmembrane helices is thought to be important for transducing ligand binding into receptor activation. Potential phosphorylation sites and a calmodulin binding site in the second intracellular loop of Orco are also thought to be involved in regulating channel gating. A number of new methods have recently been developed to express and purify insect odorant receptor subunits in recombinant expression systems. These approaches are enabling high throughput screening of receptors for agonists and antagonists in cell-based formats, as well as producing protein for the application of biophysical methods to resolve the 3D structure of the subunits and their complexes.

The Importance of Ligand-Receptor Conformational Pairs in Stabilization: Spotlight on the N/OFQ G Protein-Coupled Receptor

Understanding the mechanism by which ligands affect receptor conformational equilibria is key in accelerating membrane protein structural biology. In the case of G protein-coupled receptors (GPCRs), we currently pursue a brute-force approach for identifying ligands that stabilize receptors and facilitate crystallogenesis. The nociceptin/orphanin FQ peptide receptor (NOP) is a member of the opioid receptor subfamily of GPCRs for which many structurally diverse ligands are available for screening. We observed that antagonist potency is correlated with a ligand's ability to induce receptor stability (Tm) and crystallogenesis. Using this screening strategy, we solved two structures of NOP in complex with top candidate ligands SB-612111 and C-35. Docking studies indicate that while potent, stabilizing antagonists strongly favor a single binding orientation, less potent ligands can adopt multiple binding modes, contributing to their low Tm values. These results suggest a mechanism for ligand-aided crystallogenesis whereby potent antagonists stabilize a single ligand-receptor conformational pair.

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

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

Insights into a defined secondary binding region on β-adrenoceptors and putative roles in ligand binding and drug design

Putative roles of a secondary binding region shared among beta-adrenoceptors.

Development of a highly selective allosteric antagonist radioligand for the type 1 cholecystokinin receptor and elucidation of its molecular basis of binding

Understanding the molecular basis of ligand binding to receptors provides insights useful for rational drug design. This work describes development of a new antagonist radioligand of the type 1 cholecystokinin receptor (CCK1R), (2-fluorophenyl)-2,3-dihydro-3-[(3-isoquinolinylcarbonyl)amino]-6-methoxy-2-oxo-l-H-indole-3-propanoate (T-0632), and exploration of the molecular basis of its binding. This radioligand bound specifically with high affinity within an allosteric pocket of CCK1R. T-0632 fully inhibited binding and action of CCK at this receptor, while exhibiting no saturable binding to the closely related type 2 cholecystokinin receptor (CCK2R). Chimeric CCK1R/CCK2R constructs were used to explore the molecular basis of T-0632 binding. Exchanging exonic regions revealed the functional importance of CCK1R exon 3, extending from the bottom of transmembrane segment (TM) 3 to the top of TM5, including portions of the intramembranous pocket as well as the second extracellular loop region (ECL2). However, CCK1R mutants in which each residue facing the pocket was changed to that present in CCK2R had no negative impact on T-0632 binding. Extending the chimeric approach to ECL2 established the importance of its C-terminal region, and site-directed mutagenesis of each nonconserved residue in this region revealed the importance of Ser(208) at the top of TM5. A molecular model of T-0632-occupied CCK1R was consistent with these experimental determinants, also identifying Met(121) in TM3 and Arg(336) in TM6 as important. Although these residues are conserved in CCK2R, mutating them had a distinct impact on the two closely related receptors, suggesting differential orientation. This establishes the molecular basis of binding of a highly selective nonpeptidyl allosteric antagonist of CCK1R, illustrating differences in docking that extend beyond determinants attributable to distinct residues lining the intramembranous pocket in the two receptor subtypes.

Constitutively active chemokine CXC receptors

Chemokines are low-molecular-weight, secreted proteins that act as leukocyte-specific chemoattractants. The chemokine family has more than 40 members. Based on the position of two conserved cysteines in the N-terminal domain, chemokines can be divided into the CXC, C, CC, and CX3C subfamilies. The interaction of chemokines with their receptors mediates signaling pathways that play critical roles in cell migration, differentiation, and proliferation. The receptors for chemokines are G protein-coupled receptors (GPCRs), and thus far, seven CXC receptors have been cloned and are designated CXCR1-7. Constitutively active GPCRs are present in several human immune-mediated diseases and in tumors, and they have provided valuable information in understanding the molecular mechanism of GPCR activation. Several constitutively active CXC chemokine receptors include the V6.40A and V6.40N mutants of CXCR1; the D3.49V variant of CXCR2; the N3.35A, N3.35S, and T2.56P mutants of CXCR3; the N3.35 mutation of CXCR4; and the naturally occurring KSHV-GPCR. Here, we review the regulation of CXC chemokine receptor signaling, with a particular focus on the constitutive activation of these receptors and the implications in physiological conditions and in pathogenesis. Understanding the mechanisms behind the constitutive activation of CXC chemokine receptors may aid in pharmaceutical design and the screening of inverse agonists and allosteric modulators for the treatment of autoimmune diseases and cancers.