Structural organization of G-protein-coupled receptors

Structural Space of the Duffy Antigen/Receptor for Chemokines' Intrinsically Disordered Ectodomain 1 Explored by Temperature Replica-Exchange Molecular Dynamics Simulations

Plasmodium vivax malaria affects 14 million people each year. Its invasion requires interactions between the parasitic Duffy-binding protein (PvDBP) and the N-terminal extracellular domain (ECD1) of the host's Duffy antigen/receptor for chemokines (DARC). ECD1 is highly flexible and intrinsically disordered, therefore it can adopt different conformations. We computationally modeled the challenging ECD1 local structure. With T-REMD simulations, we sampled its dynamic behavior and collected its most representative conformations. Our results suggest that most of the DARC ECD1 domain remains in a disordered state during the simulated time. Globular local conformations are found in the analyzed local free-energy minima. These globular conformations share an α-helix spanning residues Ser18 to Ser29 and in many cases they comprise an antiparallel β-sheet, whose β-strands are formed around residues Leu10 and Ala49. The formation of a parallel β-sheet is almost negligible. So far, progress in understanding the mechanisms forming the basis of the P. vivax malaria infection of reticulocytes has been hampered by experimental difficulties, along with a lack of DARC structural information. Our collection of the most probable ECD1 structural conformations will help to advance modeling of the DARC structure and to explore DARC-ECD1 interactions with a range of physiological and pathological ligands.

Bioinformatics and evolution of vertebrate nociceptin and opioid receptors

G protein-coupled receptors (GPCRs) are ancestrally related membrane proteins on cells that mediate the pharmacological effect of most drugs and neurotransmitters. GPCRs are the largest group of membrane receptor proteins encoded in the human genome. One of the most famous types of GPCRs is the opioid receptors. Opioid family receptors consist of four closely related proteins expressed in all vertebrate brains and spinal cords examined to date. The three classical types of opioid receptors shown unequivocally to mediate analgesia in animal models and in humans are the mu- (MOR), delta- (DOR), and kappa-(KOR) opioid receptor proteins. The fourth and most recent member of the opioid receptor family discovered is the nociceptin or orphanin FQ receptor (ORL). The role of ORL and its ligands in producing analgesia is not as clear, with both analgesic and hyperalgesic effects reported. All four opioid family receptor genes were cloned from expressed mRNA in a number of vertebrate species, and there are enough sequences presently available to carry out bioinformatic analysis. This chapter presents the results of a comparative analysis of vertebrate opioid receptors using pharmacological studies, bioinformatics, and the latest data from human whole-genome studies. Results confirm our initial hypotheses that the four opioid receptor genes most likely arose by whole-genome duplication, that there is an evolutionary vector of opioid receptor type divergence in sequence and function, and that the hMOR gene shows evidence of positive selection or adaptive evolution in Homo sapiens.

High- and low-affinity sites for sodium in δ-OR-Gi1α (Cys (351)-Ile (351)) fusion protein stably expressed in HEK293 cells; functional significance and correlation with biophysical state of plasma membrane

The effect of sodium, potassium, and lithium on δ-opioid receptor ligand binding parameters and coupling with the cognate G proteins was compared in model HEK293 cell line stably expressing PTX-insensitive δ-OR-Gi1α (Cys(351)-Ile(351)) fusion protein. Agonist [(3)H]DADLE binding was decreased in the order Na(+) ≫ Li(+) > K(+) > (+)NMDG. When plotted as a function of increasing NaCl concentrations, the binding was best-fitted with a two-phase exponential decay considering two Na(+)-responsive sites (r (2) = 0.99). High-affinity Na(+)-sites were characterized by Kd = 7.9 mM and represented 25 % of the basal level determined in the absence of ions. The remaining 75 % represented the low-affinity sites (Kd = 463 mM). Inhibition of [(3)H]DADLE binding by lithium, potassium, and (+)-NMDG proceeded in low-affinity manner only. Surprisingly, the affinity/potency of DADLE-stimulated [(35)S]GTPγS binding was increased in a reverse order: Na(+) < K(+) < Li(+). This result was demonstrated in PTX-treated as well as PTX-untreated cells. Therefore, it is not restricted to Gi1α(Cys(351)-Ile(351)) within the δ-OR-Gi1α fusion protein, but is also valid for stimulation of endogenous G proteins of Gi/Go family in HEK293 cells. Biophysical studies of interaction of ions with polar head-group region of lipids using Laurdan generalized polarization indicated the low-affinity type of interaction only proceeding in the order: Cs(+) < K(+) < Na(+) < Li(+). The results are discussed in terms of interaction of Na(+), K(+) and Li(+) with the high- and low-affinity sites located in water-accessible part of δ-OR binding pocket. We also consider the role of negatively charged Cl(-), Br(-), and I(-) counter anions in inhibition of both [(3)H]DADLE and [(35)S]GTPγS binding.

Comparative domain modeling of human EGF-like module EMR2 and study of interaction of the fourth domain of EGF with chondroitin 4-sulphate

EMR2 is an EGF-like module containing mucin-like hormone receptor-2 precursor, a G-protein coupled receptor (G-PCR). Mutation in EMR2 causes complicated disorders like polycystic kidney disease (PKD). The structure of EMR2 shows that the fifth domain is comprised of EGF-TM7 helices. Functional assignment of EMR2 by support vector machine (SVM) revealed that along with transporter activity, several novel functions are predicted. A twenty amino acid sequence "MGGRVFLVFLAFCVWLTLPG" acts as the signal peptide responsible for posttranslational transport. Eight amino acids are involved in N-glycosylation sites and two cleavage sites are Leu517 and Ser518 in EMR2. The residue Arg241 is responsible for interaction with glycosaminoglycan and chondroitin sulfate. On the basis of structure, function and ligand binding sites, competitive EMR2 inhibitors designed may decrease the rate of human diseases like Usher's syndrome, bilateral frontoparietal polymicrogyria and PKD.

FoldGPCR: structure prediction protocol for the transmembrane domain of G protein-coupled receptors from class A

Building reliable structural models of G protein-coupled receptors (GPCRs) is a difficult task because of the paucity of suitable templates, low sequence identity, and the wide variety of ligand specificities within the superfamily. Template-based modeling is known to be the most successful method for protein structure prediction. However, refinement of homology models within 1-3 A C alpha RMSD of the native structure remains a major challenge. Here, we address this problem by developing a novel protocol (foldGPCR) for modeling the transmembrane (TM) region of GPCRs in complex with a ligand, aimed to accurately model the structural divergence between the template and target in the TM helices. The protocol is based on predicted conserved inter-residue contacts between the template and target, and exploits an all-atom implicit membrane force field. The placement of the ligand in the binding pocket is guided by biochemical data. The foldGPCR protocol is implemented by a stepwise hierarchical approach, in which the TM helical bundle and the ligand are assembled by simulated annealing trials in the first step, and the receptor-ligand complex is refined with replica exchange sampling in the second step. The protocol is applied to model the human beta(2)-adrenergic receptor (beta(2)AR) bound to carazolol, using contacts derived from the template structure of bovine rhodopsin. Comparison with the X-ray crystal structure of the beta(2)AR shows that our protocol is particularly successful in accurately capturing helix backbone irregularities and helix-helix packing interactions that distinguish rhodopsin from beta(2)AR.

Human Cytomegalovirus US28: a functionally selective chemokine binding receptor

Chemokines are small cytokines that are part of a large family of molecules that bind to G-protein coupled receptors, which, as a family, are the most widely targeted group of molecules in the treatment of disease. Chemokines are critical for recruiting and activating the cells of the immune system during inflammation especially during viral infections. However, a number of viruses including the large herpes virus human cytomegalovirus (HCMV) encode mechanisms to impede the effects of chemokines or has gained the ability to use these molecules to its own advantage. The Human Cytomegalovirus (HCMV)-encoded chemokine receptor US28 is the best characterized of the four unique chemokine receptor-like molecules found in the HCMV genome. US28 has been studied as an important virulence factor for HCMV-mediated vascular disease and, more recently, in models of HCMV-associated malignancy. US28 is a rare multi-chemokine family binding receptor with the ability to bind ligands from two distinct chemokine classes. Ligand binding to US28 activates cell-type and ligand-specific signaling pathways leading to cellular migration, which is an important example of receptor functional selectivity. Additionally, US28 has been demonstrated to constitutively activate phospholipase C (PLC) and NF-kB signaling pathways. Understanding the structure/function relationships between US28, its ligands and intracellular signaling molecules will provide essential clues for effective pharmacological targeting of this multifunctional chemokine receptor.

Comparative docking study of anibamine as the first natural product CCR5 antagonist in CCR5 homology models

Anibamine, a novel pyridine quaternary alkaloid recently isolated from Aniba sp., has been found to effectively bind to the chemokine receptor CCR5 with an IC(50) at 1 microM in competition with (125)I-gp120, an HIV viral envelope protein binding to CCR5 with high affinity. Since CCR5, a G-protein-coupled receptor, is an essential coreceptor for the human immunodeficiency virus type I (HIV-1) entry to host cells, a CCR5 antagonist that inhibits the cellular entry of HIV-1 provides a new therapy choice for the treatment of HIV. Anibamine provides a novel structural skeleton that is remarkably different from all lead compounds previously identified as CCR5 antagonists. Here, we report comparative docking studies of anibamine with several other known CCR5 antagonists in two CCR5 homology models built based on the crystal structures of bovine rhodopsin and human beta(2)-adrenergic receptor. The binding pocket of anibamine has some common features shared with other high affinity CCR5 antagonists, suggesting that they may bind in similar binding sites and/or modes. At the same time, several unique binding features of anibamine were identified, and it will likely prove beneficial in future molecular design of novel CCR5 antagonists based on the anibamine scaffold.

The evolution of vertebrate opioid receptors

The proteins that mediate the analgesic and other effects of opioid drugs and endogenous opioid peptides are known as opioid receptors. Opioid receptors consist of a family of four closely-related proteins belonging to the large superfamily of G-protein coupled receptors. The three types of opioid receptors shown unequivocally to mediate analgesia in animal models are the mu (MOR), delta (DOR), and kappa (KOR) opioid receptor proteins. The role of the fourth member of the opioid receptor family, the nociceptin or orphanin FQ receptor (ORL), is not as clear as hyperalgesia, analgesia, and no effect was reported after administration of ORL agonists. There are now cDNA sequences for all four types of opioid receptors that are expressed in the brain of six species from three different classes of vertebrates. This review presents a comparative analysis of vertebrate opioid receptors using bioinformatics and data from recent human genome studies. Results indicate that opioid receptors arose by gene duplication, that there is a vector of opioid receptor divergence, and that MOR shows evidence of rapid evolution.

Modelling the structures of G protein-coupled receptors aided by three-dimensional validation

Background

G protein-coupled receptors (GPCRs) are abundant, activate complex signalling and represent the targets for up to ~60% of pharmaceuticals but there is a paucity of structural data. Bovine rhodopsin is the first GPCR for which high-resolution structures have been completed but significant variations in structure are likely to exist among the GPCRs. Because of this, considerable effort has been expended on developing in silico tools for refining structures of individual GPCRs. We have developed REPIMPS, a modification of the inverse-folding software Profiles-3D, to assess and predict the rotational orientation and vertical position of helices within the helix bundle of individual GPCRs. We highlight the value of the method by applying it to the Baldwin GPCR template but the method can, in principle, be applied to any low- or high-resolution membrane protein template or structure.

Results

3D models were built for transmembrane helical segments of 493 GPCRs based on the Baldwin template, and the models were then scored using REPIMPS and Profiles-3D. The compatibility scores increased significantly using REPIMPS because it takes into account the physicochemical properties of the (lipid) environment surrounding the helix bundle. The arrangement of helices in the helix bundle of the 493 models was then altered systematically by rotating the individual helices. For most GPCRs in the set, changes in the rotational position of one or more helices resulted in significant improvement in the compatibility scores. In particular, for most GPCRs, a rotation of helix VII by 240–300° resulted in improved scores. Bovine rhodopsin modelled using this method showed 3.31 Å RMSD to its crystal structure for 198 C^α atom pairs, suggesting the utility of the method even when starting with idealised structures such as the Baldwin template.

Conclusion

We have developed an in silico tool which can be used to test the validity of, and refine, models of GPCRs with respect to helix rotation and vertical position based on the physicochemical properties of amino acids and the surrounding environment. The method can be applied to any multi-pass membrane protein and potentially can be used in combination with other high-throughput methodologies to generate and refine models of membrane proteins.

Evidence for a strong selenium-aromatic interaction derived from crystallographic data and ab initio quantum chemical calculations

Increasing attention is being paid to the role of selenium, both as an essential component required for the activity of many enzymes and in the context of selenium-based pharmaceutical agents. A wide range of therapeutics that include selenium are on the market and under development, such as antihypertensive, anticancerogenic, antiviral, and immunosuppressive agents. Computer-aided drug design (CADD) has proven to be an important tool for the development of new drugs. Many CADD techniques, including docking, molecular dynamics simulation, and other receptor-based approaches, require an accurate understanding of the nature of the intermolecular forces that act to stabilize protein-ligand complexes; moreover, a quantitative assessment of these interactions furthers our efforts to rationalize the drug design process. In this paper, we consider one class of interaction involving selenium, that between Se and aromatic rings. Prior work has shown that interactions between divalent sulfur and aromatic rings are observed much more frequently than would be expected on the basis of chance, both in protein structures and the crystal structures of organic compounds that include these moieties. Recent studies on the optimization of inhibitor-protein binding also suggest that sulfur-aromatic interactions are important in stabilizing these complexes and may be crucial focal point for CADD. Given that selenium and sulfur have similar chemistry, and that selenium is significantly more polarizable, we propose that Se-aromatic interactions may also play an important stabilizing role in the structure of folded proteins and in drug-protein complexes. We have tested this hypothesis against data from the Cambridge Crystallographic Database and ab initio quantum chemical calculations. We have found evidence that selenium does interact strongly with aromatic rings and may play a role analogous to sulfur in stabilizing protein folds. In addition, selenium should be considered along with sulfur in rational drug design strategies that seek to improve binding to target protein sites that include aromatic rings.

Location of the hydrophobic pocket in the binding site of fentanyl analogs in the µ-opioid receptor

Fentanyl is a highly potent and clinically widely used narcotic analgesic. The synthesis of its analogs remains a challenge in an attempt to develop highly selective ?-opioid receptor agonists with specific pharmacological properties. In this paper, the use of flexible molecular docking of several specific fentanyl analogs to the ?-opioid receptor model, in order to test the hypothesis that the hydrophobic pocket accommodates alkyl groups at position 3 of the fentanyl skeleton, is described. The stereoisomers of the following compounds were studied: cis- and trans-3-methylfentanyl, 3,3-dimethylfentanyl, cis- and trans-3-ethylfentanyl, cis- and trans-3-propylfentanyl, cis-3-isopropylfentanyl and cis-3-benzylfentanyl. The optimal position and orientation of these fentanyl analogs in the binding pocket of the ?-receptor, explaining their enantiospecific potency, were determined. It was found that the 3-alkyl group of cis-3R,4S and trans-3S,4S stereoisomers of all the active compounds occupies the hydrophobic pocket between TM5, TM6 and TM7, made up of the amino acids Trp318 (TM7), Ile322 (TM7), Ile301 (TM6) and Phe237 (TM5). However, the fact that this hydrophobic pocket can also accommodate the bulky 3-alkyl substituents of the two inactive compounds: cis-3-isopropylfentanyl, and cis-3-benzylfentanyl, indicates that this hydrophobic pocket in the employed receptor model is probably too large. .

The chemokine receptor homologue encoded by US27 of human cytomegalovirus is heavily glycosylated and is present in infected human foreskin fibroblasts and enveloped virus particles

Human cytomegalovirus (HCMV), a member of the beta-herpesvirus family, encodes four homologues of cellular G protein-coupled receptors (GPCRs). One of these, the protein product of HCMV open reading frame (ORF) UL33, has been identified in HCMV-infected cells and virus particles and shown to be heat-aggregatable and N-glycosylated. Another, the product of ORF US28, has been functionally characterized as a beta-chemokine receptor. Here we report the use of RT-PCR, coupled in vitro transcription-translation, immunoprecipitation, and Western immunoassays to (i) show that RNA from the open reading frame US27 appears predominantly during the late phase of replication; (ii) identify the protein encoded by HCMV US27 in infected cells and enveloped virus particles; (iii) demonstrate that the US27-encoded protein is heterogeneously N-glycosylated and resolves as two species following treatment with peptide N-glycosidase F; and (iv) show that both the recombinant and deglycoylated infected cell US27 protein aggregate when heated in the presence of SDS prior to electrophoresis in polyacrylamide gels, a property which is abrogated with the addition of urea to sample buffer.

Homology modeling of opioid receptor-ligand complexes using experimental constraints

Opioid receptors interact with a variety of ligands, including endogenous peptides, opiates, and thousands of synthetic compounds with different structural scaffolds. In the absence of experimental structures of opioid receptors, theoretical modeling remains an important tool for structure-function analysis. The combination of experimental studies and modeling approaches allows development of realistic models of ligand-receptor complexes helpful for elucidation of the molecular determinants of ligand affinity and selectivity and for understanding mechanisms of functional agonism or antagonism. In this review we provide a brief critical assessment of the status of such theoretical modeling and describe some common problems and their possible solutions. Currently, there are no reliable theoretical methods to generate the models in a completely automatic fashion. Models of higher accuracy can be produced if homology modeling, based on the rhodopsin X-ray template, is supplemented by experimental structural constraints appropriate for the active or inactive receptor conformations, together with receptor-specific and ligand-specific interactions. The experimental constraints can be derived from mutagenesis and cross-linking studies, correlative replacements of ligand and receptor groups, and incorporation of metal binding sites between residues of receptors or receptors and ligands. This review focuses on the analysis of similarity and differences of the refined homology models of mu, delta, and kappa-opioid receptors in active and inactive states, emphasizing the molecular details of interaction of the receptors with some representative peptide and nonpeptide ligands, underlying the multiple modes of binding of small opiates, and the differences in binding modes of agonists and antagonists, and of peptides and alkaloids.

Location and nature of the residues important for ligand recognition in G-protein coupled receptors

The overall structure of the biogenic amine subclass of the G-protein-coupled receptors, and of their ligand binding sites, is discussed with the aim of highlighting the major structural features of these receptors that are responsible for ligand recognition. A comparison is made between biogenic amine receptors, peptide receptors of the rhodopsin class, and the secretin receptors which all have peptide ligands. The question of where the peptide ligands bind, whether at extracellular sites or within the transmembrane helix bundle, is discussed. The suitability of the rhodopsin crystal structure as a template for construction of homology models is discussed and it is concluded that there are many reasons why a caution should be issued against using it uncritically.

High-throughput modeling of human G-protein coupled receptors: amino acid sequence alignment, three-dimensional model building, and receptor library screening

The current study describes the development of a computer package (GPCRmod) aimed at the high-throughput modeling of the therapeutically important family of human G-protein coupled receptors (GPCRs). GPCRmod first proposes a reliable alignment of the seven transmembrane domains (7 TMs) of most druggable human GPCRs based on pattern/motif recognition for each of the 7 TMs that are considered independently. It then converts the alignment into knowledge-based three-dimensional (3-D) models starting from a set of 3-D backbone templates and two separate rotamer libraries for side chain positioning. The 7 TMs of 277 human GPCRs have been accurately aligned, unambiguously clustered in three different classes (rhodopsin-like, secretin-like, metabotropic glutamate-like), and converted into high-quality 3-D models at a remarkable throughput (ca. 3s/model). A 3-D GPCR target library of 277 receptors has consequently been setup. Its utility for "in silico" inverse screening purpose has been demonstrated by recovering among top scorers the receptor of a selective GPCR antagonist as well as the receptors of a promiscuous antagonist. The current GPCR target library thus constitutes a 3-D database of choice to address as soon as possible the "virtual selectivity" profile of any GPCR antagonist or inverse agonist in an early hit optimization process.

Membrane model for the G-protein-coupled receptor rhodopsin: hydrophobic interface and dynamical structure

Rhodopsin is the only member of the pharmacologically important superfamily of G-protein-coupled receptors with a known structure at atomic resolution. A molecular dynamics model of rhodopsin in a POPC phospholipid bilayer was simulated for 15 ns, revealing a conformation significantly different from the recent crystal structures. The structure of the bilayer compared with a protein-free POPC control indicated hydrophobic matching with the nonpolar interface of the receptor, in agreement with deuterium NMR experiments. A new generalized molecular surface method, based on a three-dimensional Voronoi cell construction for atoms with different radii, was developed to quantify cross-sectional area profiles for the protein, lipid acyl chains and headgroups, and water. Thus, it was possible to investigate the bilayer deformation due to curvature of the individual lipid monolayers. Moreover, the generalized molecular surface derived hydrophobic interface allowed benchmarking of the hydropathy sequence analysis, an important structural genomics tool. Five water molecules diffused into internal hydration sites during the simulation, yielding a total of 12 internal waters. The cytoplasmic loops and the C-terminal tail, containing the G-protein recognition and protein sorting sequences, exhibited a high mobility, in marked contrast to the extracellular and transmembrane domains. The proposed functional coupling of the highly conserved ERY motif to the lipid-water interface via the cytoplasmic loops provides insight into lipid effects on G-protein-coupled receptor activation in terms of a flexible surface model, involving the spontaneous monolayer curvature.

First principles predictions of the structure and function of g-protein-coupled receptors: validation for bovine rhodopsin

G-protein-coupled receptors (GPCRs) are involved in cell communication processes and with mediating such senses as vision, smell, taste, and pain. They constitute a prominent superfamily of drug targets, but an atomic-level structure is available for only one GPCR, bovine rhodopsin, making it difficult to use structure-based methods to design receptor-specific drugs. We have developed the MembStruk first principles computational method for predicting the three-dimensional structure of GPCRs. In this article we validate the MembStruk procedure by comparing its predictions with the high-resolution crystal structure of bovine rhodopsin. The crystal structure of bovine rhodopsin has the second extracellular (EC-II) loop closed over the transmembrane regions by making a disulfide linkage between Cys-110 and Cys-187, but we speculate that opening this loop may play a role in the activation process of the receptor through the cysteine linkage with helix 3. Consequently we predicted two structures for bovine rhodopsin from the primary sequence (with no input from the crystal structure)-one with the EC-II loop closed as in the crystal structure, and the other with the EC-II loop open. The MembStruk-predicted structure of bovine rhodopsin with the closed EC-II loop deviates from the crystal by 2.84 A coordinate root mean-square (CRMS) in the transmembrane region main-chain atoms. The predicted three-dimensional structures for other GPCRs can be validated only by predicting binding sites and energies for various ligands. For such predictions we developed the HierDock first principles computational method. We validate HierDock by predicting the binding site of 11-cis-retinal in the crystal structure of bovine rhodopsin. Scanning the whole protein without using any prior knowledge of the binding site, we find that the best scoring conformation in rhodopsin is 1.1 A CRMS from the crystal structure for the ligand atoms. This predicted conformation has the carbonyl O only 2.82 A from the N of Lys-296. Making this Schiff base bond and minimizing leads to a final conformation only 0.62 A CRMS from the crystal structure. We also used HierDock to predict the binding site of 11-cis-retinal in the MembStruk-predicted structure of bovine rhodopsin (closed loop). Scanning the whole protein structure leads to a structure in which the carbonyl O is only 2.85 A from the N of Lys-296. Making this Schiff base bond and minimizing leads to a final conformation only 2.92 A CRMS from the crystal structure. The good agreement of the ab initio-predicted protein structures and ligand binding site with experiment validates the use of the MembStruk and HierDock first principles' methods. Since these methods are generic and applicable to any GPCR, they should be useful in predicting the structures of other GPCRs and the binding site of ligands to these proteins.

Homology modelling and binding site mapping of the human histamine H1 receptor

Three-dimensional model of the human histamine H1 receptor was developed by homology modelling using the high resolution structure of bovine rhodopsin as template. Genetic algorithm based docking calculations were used to identify the role of several amino acids having an effect on agonist or antagonist binding. Binding mode analyses of mepyramine, desloratidine, loratidine and acrivastine allowed us to rationalise their binding affinity. Binding site mapping resulted in seven new potential aromatic interaction points (Tyr 108, Phe 184, Phe 190, Phe 199, Phe 424, Trp 428, Tyr 431), that took part in forming the lipophilic pocket of the antagonist binding cavity.

Novel approaches for modeling of the A1adenosine receptor and its agonist binding site

The present work describes the building of a human A(1) adenosine receptor (hA(1)AR) model, based on the X-ray crystal structure of bovine rhodopsin, and its use as a basis for the investigation of some important structural characteristics of the receptor. One of the issues investigated was the protonation position of two histidine residues known to influence ligand binding, with protonation of His251 (6.52) in epsilon position and His278 (7.43) in delta position showing the best agreement with experimental evidence. The model was also used to study the position and structural role of water molecules present in the helical bundle. Finally, the binding site location and the ligand docking were investigated using an objective strategy. A suitable site for the binding of the ribose moiety of adenosine was first postulated and further confirmed by means of a novel chemometric strategy based on GRIND descriptors. Using this position as an anchor point, the binding of adenosine was studied by docking and molecular dynamics simulations obtaining two putative binding positions in good agreement with experimental data.

Development and validation of opioid ligand-receptor interaction models: the structural basis of mu vs delta selectivity

Opioid receptor binding conformations for two structurally related, conformationally constrained tetrapeptides, JOM-6 ( micro receptor selective) and JOM-13 (delta receptor selective), were deduced using conformational analysis of these ligands and analogs with additional conformational restrictions. Docking of these ligands in their binding conformations to opioid receptor structural models, based upon the published rhodopsin X-ray structure, implicates specific structural features of the micro and delta receptor ligand binding sites as forming the basis for the micro selectivity of JOM-6 and the delta selectivity of JOM-13. In particular, the presence of E229 in the micro receptor (in place of the corresponding D210 of the delta receptor) causes an adverse electrostatic interaction with C-terminal carboxylate-containing ligands, resulting in the observed preference of ligands with an uncharged C-terminus for the micro receptor. In addition, the requirement that the Phe3 side chain of JOM-13 assume a gauche orientation for optimal delta binding, whereas the Phe3 side chain of JOM-6 must be in a trans orientation for high-affinity micro binding can be largely attributed to the steric effect of replacement of L300 of the delta receptor by W318 of the micro receptor. Testing this hypothesis by examining the binding of JOM-6 and several of its key analogs with specific micro receptor mutants is described. Our initial results are consistent with the proposed ligand-receptor interaction models.

The eye photoreceptor protein rhodopsin. Structural implications for retinal disease

Rhodopsin is the membrane receptor responsible for photoreception in the vertebrate retina. Its characteristic seven-transmembrane helical structural motif is today widely recognised as a paradigm in signal transduction. Rhodopsin and the phototransduction system are frequently used as structural and mechanistic models for the G-protein coupled receptor superfamily. Recent advances in the activation mechanism (as derived from the structural available data) and the implications for normal and pathological - in retinal disorders - visual function will be reviewed.

Interatomic potentials and solvation parameters from protein engineering data for buried residues

Van der Waals (vdW) interaction energies between different atom types, energies of hydrogen bonds (H-bonds), and atomic solvation parameters (ASPs) have been derived from the published thermodynamic stabilities of 106 mutants with available crystal structures by use of an originally designed model for the calculation of free-energy differences. The set of mutants included substitutions of uncharged, inflexible, water-inaccessible residues in alpha-helices and beta-sheets of T4, human, and hen lysozymes and HI ribonuclease. The determined energies of vdW interactions and H-bonds were smaller than in molecular mechanics and followed the "like dissolves like" rule, as expected in condensed media but not in vacuum. The depths of modified Lennard-Jones potentials were -0.34, -0.12, and -0.06 kcal/mole for similar atom types (polar-polar, aromatic-aromatic, and aliphatic-aliphatic interactions, respectively) and -0.10, -0.08, -0.06, -0.02, and nearly 0 kcal/mole for different types (sulfur-polar, sulfur-aromatic, sulfur-aliphatic, aliphatic-aromatic, and carbon-polar, respectively), whereas the depths of H-bond potentials were -1.5 to -1.8 kcal/mole. The obtained solvation parameters, that is, transfer energies from water to the protein interior, were 19, 7, -1, -21, and -66 cal/moleA(2) for aliphatic carbon, aromatic carbon, sulfur, nitrogen, and oxygen, respectively, which is close to the cyclohexane scale for aliphatic and aromatic groups but intermediate between octanol and cyclohexane for others. An analysis of additional replacements at the water-protein interface indicates that vdW interactions between protein atoms are reduced when they occur across water.

Electrostatic properties of membrane lipids coupled to metarhodopsin II formation in visual transduction

Changes in lipid composition have recently been shown to exert appreciable influences on the activities of membrane-bound proteins and peptides. We tested the hypothesis that the conformational states of rhodopsin linked to visual signal transduction are related to biophysical properties of the membrane lipid bilayer. For bovine rhodopsin, the meta I-meta II conformational transition was studied in egg phosphatidylcholine (PC) recombinants versus the native rod outer segment (ROS) membranes by means of flash photolysis. Formation of metarhodopsin II was observed by the change in absorbance at 478 nm after a single actinic flash was delivered to the sample. The meta I/meta II ratio was investigated as a function of both temperature and pH. The data clearly demonstrated thermodynamic reversibility of the transition for both the egg PC recombinants and the native ROS membranes. A significant shift of the apparent pK(a) for the acid-base equilibrium to lower values was evident in the egg PC recombinant, with little meta II produced under physiological conditions. Calculations of the membrane surface pH using a Poisson-Boltzmann model suggested the free energies of the meta I and meta II states were significantly affected by electrostatic properties of the bilayer lipids. In the ROS membranes, phosphatidylserine (PS) is needed for full formation of meta II, in combination with phosphatidylethanolamine (PE) and polyunsaturated docosahexaenoic acid (DHA; 22:6omega3) chains. We propose that the PS surface potential leads to an accumulation of hydronium ions, H(3)O(+), in the electrical double layer, which drive the reaction together with the large negative spontaneous curvature (H(0)) conferred by PE plus DHA chains. The elastic stress/strain of the bilayer arises from an interplay of the approximately zero H(0) from PS and the negative H(0) due to the PE headgroups and polyunsaturated chains. The lipid influences are further explained in terms of matching of the bilayer spontaneous curvature to the curvature at the lipid/rhodopsin interface, as formulated by the Helfrich bending energy. These new findings guide current ideas as to how bilayer properties govern the conformational energetics of integral membrane proteins. Moreover, they yield knowledge of how membrane lipid-protein interactions involving acidic phospholipids such as PS and neutral polyunsaturated DHA chains are implicated in key biological functions such as vision.

Chemokines and chemokine receptors in leukocyte trafficking

Chemokines regulate inflammation, leukocyte trafficking, and immune cell differentiation. The role of chemokines in homing of naive T lymphocytes to secondary lymphatic organs is probably the best understood of these processes, and information on chemokines in inflammation, asthma, and neurological diseases is rapidly increasing. Over the past 15 years, understanding of the size and functional complexity of the chemokine family of peptide chemoattractants has grown substantially. In this review, we first present information regarding the structure, expression, and signaling properties of chemokines and their receptors. The second part is a systems physiology-based overview of the roles that chemokines play in tissue-specific homing of lymphocyte subsets and in trafficking of inflammatory cells. This review draws on recent experimental findings as well as current models proposed by experts in the chemokine field.

NMR structure of the second intracellular loop of the alpha 2A adrenergic receptor: evidence for a novel cytoplasmic helix

A major, unresolved question in signal transduction by G protein coupled receptors (GPCRs) is to understand how, at atomic resolution, a GPCR activates a G protein. A step toward answering this question was made with the determination of the high-resolution structure of rhodopsin; we now know the intramolecular interactions that characterize the resting conformation of a GPCR. To what degree does this structure represent a structural paradigm for other GPCRs, especially at the cytoplasmic surface where GPCR-G protein interaction occurs and where the sequence homology is low among GPCRs? To address this question, we performed NMR studies on approximately 35-residue-long peptides including the critical second intracellular loop (i2) of the alpha 2A adrenergic receptor (AR) and of rhodopsin. To stabilize the secondary structure of the peptide termini, 4-12 residues from the adjacent transmembrane helices were included and structures determined in dodecylphosphocholine micelles. We also characterized the effects on an alpha 2A AR peptide of a D130I mutation in the conserved DRY motif. Our results show that in contrast to the L-shaped loop in the i2 of rhodopsin, the i2 of the alpha 2A AR is predominantly helical, supporting the hypothesis that there is structural diversity within GPCR intracellular loops. The D130I mutation subtly modulates the helical structure. The spacing of nonpolar residues in i2 with helical periodicity is a predictor of helical versus loop structure. These data should lead to more accurate models of the intracellular surface of GPCRs and of receptor-mediated G protein activation.

Molecular dynamics of 5-HT1A and 5-HT2A serotonin receptors with methylated buspirone analogues

In the present study experimentally determined ligand selectivity of three methylated buspirone analogues (denoted as MM2, MM5 and P55) towards 5-HT1A and 5-HT2A serotonin receptors was theoretically investigated on a molecular level. The relationships between the ligand structure and 5-HT1A and 5-HT2A receptor affinities were studied and the results were found to be in agreement with the available site-directed mutagenesis and binding affinity data. Molecular dynamics (MD) simulations of ligand-receptor complexes were performed for each investigated analogue, docked twice into the central cavity of 5-HT1A/5-HT2A, each time in a different orientation. Present results were compared with our previous theoretical results, obtained for buspirone and its non-methylated analogues. It was found that due to the presence of the methyl group in the piperazine ring the ligand position alters and the structure of the ligand-receptor complex is modified. Further, the positions of derivatives with pyrimidinyl aromatic moiety and quinolinyl moiety are significantly different at the 5-HT2A receptor. Thus, methylation of such derivatives alters the 3D structures of ligand-receptor complexes in different ways. The ligand-induced changes of the receptor structures were also analysed. The obtained results suggest, that helical domains of both receptors have different dynamical behaviour. Moreover, both location and topography of putative binding sites for buspirone analogues are different at 5-HT1A and 5-HT2A receptors.

The role of the hydrophilic Asn230 residue of the mu-opioid receptor in the potency of various opioid agonists

1. To investigate the effect of the hydrophilic Asn amino acid at position 230 of the human mu-opioid receptor (hMOR230) on the potency of various agonists, we mutated this residue to Thr and Leu (hMORN230T and hMORN230L respectively). 2. Taking advantage of the functional coupling of the opioid receptor with the heteromultimeric G-protein-coupled inwardly rectifying K(+) (GIRK1/GIRK2) channel, either the wild type hMOR or one of the mutated receptors (hMORN230L or hMORN230T) were functionally coexpressed with GIRK1/GIRK2 channels and a regulator of G-protein signalling (RGS4) in Xenopus laevis oocytes. 3. The two-microelectrode voltage clamp technique was used to measure the opioid receptor-activated GIRK1/GIRK2 channel responses. The potency of [D-Ala(2),N-MePhe(4),Gly(5)-ol]-enkephalin (DAMGO), remained unaffected as measured via hMORN230T and hMORN230L, while the potency of fentanyl and morphine significantly increased via these mutated receptors. 4. Our results are indicative for the existence of hydrophobic interactions between a methyl-group of the side chain of Thr or Leu on the one hand and the piperidine-ring of fentanyl and the hexene-ring of morphine on the other. The mutations also had no influence on the potency of morphine-6-glucuronide (M6G) and morphine-3-glucuronide (M3G). 5. We conclude that the hydrophilic side chain of Asn in position 230 is not involved in the formation of a H-bond with the aliphatic alcohol of morphine and that an enhancement of the potency of morphine and fentanyl can be explained by mutating this residue towards more hydrophobic amino acids.

A sequence and structural study of transmembrane helices

A comparison is made between the distribution of residue preferences, three dimensional nearest neighbour contacts, preferred rotamers, helix-helix crossover angles and peptide bond angles in three sets of proteins: a non-redundant set of accurately determined globular protein structures, a set of four-helix bundle structures and a set of membrane protein structures. Residue preferences for the latter two sets may reflect overall helix stabilising propensities but may also highlight differences arising out of the contrasting nature of the solvent environments in these two cases. The results bear out the expectation that there may be differences between residue type preferences in membrane proteins and in water soluble globular proteins. For example, the beta-branched residue types valine and isoleucine are considerably more frequently encountered in membrane helices. Likewise, glycine and proline. residue types normally associated with 'helix-breaking' propensity are found to be relatively more common in membrane helices. Three dimensional nearest neighbour contacts along the helix, preferred rotamers, and peptide bond angles are very similar in the three sets of proteins as far as can be ascertained within the limits of the relatively low resolution of the membrane proteins dataset. Crossing angles for helices in the membrane protein set resemble the four helix bundle set more than the general non-redundant set, but in contrast to both sets they have smaller crossing angles consistent with the dual requirements for the helices to form a compact structure while having to span the membrane. In addition to the pairwise packing of helices we investigate their global packing and consider the question of helix supercoiling in helix bundle proteins.

Retinylidene proteins: structures and functions from archaea to humans

Retinylidene proteins, containing seven membrane-embedded alpha-helices that form an internal pocket in which the chromophore retinal is bound, are ubiquitous in photoreceptor cells in eyes throughout the animal kingdom. They are also present in a diverse range of other organisms and locations, such as archaeal prokaryotes, unicellular eukaryotic microbes, the dermal tissue of frogs, the pineal glands of lizards and birds, the hypothalamus of toads, and the human brain. Their functions include light-driven ion transport and phototaxis signaling in microorganisms, and retinal isomerization and various types of photosignal transduction in higher animals. The aims of this review are to examine this group of photoactive proteins as a whole, to summarize our current understanding of structure/function relationships in the best-studied examples, and to report recent new developments.

Exploring the relationship between binding modes of 9-(aminomethyl)-9,10-dihydroanthracene and cyproheptadine analogues at the 5-HT2A serotonin receptor

Comparison of the serotonin 5-HT2A receptor affinities of a parallel series of structural analogues of the novel ligand 9-aminomethyl-9,10-dihydroanthracene (AMDA) and a structurally similar prototypical tricyclic amine cyproheptadine suggests that the two agents bind to the receptor in different fashions. Examination of ligand-receptor model complexes supports the experimental data and suggests a potential origin for the differences in binding modes.

Synthesis and biophysical analysis of transmembrane domains of a Saccharomyces cerevisiae G protein-coupled receptor

The Ste2p receptor for alpha-factor, a tridecapeptide mating pheromone of the yeast Saccharomyces cerevisiae, belongs to the G protein-coupled family of receptors. In this paper we report on the synthesis of peptides corresponding to five of the seven transmembrane domains (M1-M5) and two homologues of the sixth transmembrane domain corresponding to the wild-type sequence and a mutant sequence found in a constitutively active receptor. The secondary structures of all new transmembrane peptides and previously synthesized peptides corresponding to domains 6 and 7 were assessed using a detailed CD analysis in trifluoroethanol, trifluoroethanol-water mixtures, sodium dodecyl sulfate micelles, and dimyristoyl phosphatidyl choline bilayers. Tryptophan fluorescence quenching experiments were used to assess the penetration of the membrane peptides into lipid bilayers. All peptides were predominantly (40-80%) helical in trifluoroethanol and most trifluoroethanol-water mixtures. In contrast, two of the peptides M3-35 (KKKNIIQVLLVASIETSLVFQIKVIFTGDNFKKKG) and M6-31 (KQFDSFHILLINleSAQSLLVPSIIFILAYSLK) formed stable beta-sheet structures in both sodium dodecyl sulfate micelles and DMPC bilayers. Polyacrylamide gel electrophoresis showed that these two peptides formed high molecular aggregates in the presence of SDS whereas all other peptides moved as monomeric species. The peptide (KKKFDSFHILLIMSAQSLLVLSIIFILAYSLKKKS) corresponding to the sequence in the constitutive mutant was predominantly helical under a variety of conditions, whereas the homologous wild-type sequence (KKKFDSFHILLIMSAQSLLVPSIIFILAYSLKKKS) retained a tendency to form beta-structures. These results demonstrate a connection between a conformational shift in secondary structure, as detected by biophysical techniques, and receptor function. The aggregation of particular transmembrane domains may also reflect a tendency for intermolecular interactions that occur in the membrane environment facilitating formation of receptor dimers or multimers.

Complementarity of delta opioid ligand pharmacophore and receptor models

The elaboration of a pharmacophore model for the delta opioid receptor selective ligand JOM-13 (Tyr-c[D-Cys-Phe-D-Pen]OH) and the parallel, independent development of a structural model of the delta receptor are summarized. Although the backbone conformation of JOM-13's tripeptide cycle is well defined, considerable conformational lability is evident in the Tyr(1) residue and in the Phe(3) side chain, key pharmacophore elements of the ligand. Replacement of these flexible features of the ligand by more conformationally restricted analogues and subsequent correlation of receptor binding and conformational properties allowed the number of possible binding conformations of JOM-13 to be reduced to two. Of these, one was chosen as more likely, based on its better superposition with other conformationally constrained delta receptor ligands. Our model of the delta opioid receptor, constructed using a general approach that we have developed for all rhodopsin-like G protein-coupled receptors, contains a large cavity within the transmembrane domain that displays excellent complementarity in both shape and polarity to JOM-13 and other delta ligands. This binding pocket, however, cannot accommodate the conformer of JOM-13 preferred from analysis of ligands, alone. Rather, only the "alternate" allowed conformer, identified from analysis of the ligands but "disfavored" because it does not permit simultaneous superposition of all pharmacophore elements of JOM-13 with other delta ligands, fits the binding site. These results argue against a simple view of a single, common fit to a receptor binding site and suggest, instead, that at least some binding site interactions of different ligands may differ.

A putative chemoattractant receptor, C5L2, is expressed in granulocyte and immature dendritic cells, but not in mature dendritic cells

Leukocyte chemoattractants are involved in the inflammatory response act via G protein-coupled receptors. We report the cloning of a novel human gene encoding the putative orphan receptor, C5L2, belonging to a subfamily of C3a, C5a and formyl Met-Leu-Ph receptors that are related to the chemokine receptor family. C5L2 transcript levels were abundant in granulocytes and immature dendritic cells but not in mature dendritic cells. We speculate that this receptor may regulate the activation of immature dendritic cells and play a role in the chemoattraction of leukocytes to inflammatory regions.