Molecular modeling of the second extracellular loop of G-protein coupled receptors and its implication on structure-based virtual screening

The current study describes the validation of high-throughput modeling procedures for the construction of the second extracellular loop (ecl2) of all nonolfactory human G Protein-coupled receptors. Our modeling flowchart is based on the alignment of essential residues determining the particular ecl2 fold observed in the bovine rhodopsin (bRho) crystal structure. For a set of GPCR targets, the dopamine D2 receptor (DRD2), adenosine A3 receptor (AA3R), and the thromboxane A2 receptor (TA2R), the implications of including ecl2 atomic coordinates is evaluated in terms of structure-based virtual screening accuracy: the suitability of the 3D models to distinguish between known antagonists and randomly chosen decoys using automated docking approaches. The virtual screening results of different models describing increasingly exhaustive receptor representations (seven helices only, seven helices and ecl2 loop, full model) have been compared. Explicit modeling of the ecl2 loop was found to be important in only one of three test cases whereas a loopless model was shown to be accurate enough in the two other receptors. An exhaustive comparison of ecl2 loops of 365 receptors to that of bRho suggests that explicit ecl2 loop modeling should be reserved to receptors where loop building can be guided by experimental restraints.     

3D structural model of the G-protein-coupled cannabinoid CB2 receptor

The potential for therapeutic specificity in regulating diseases and for reduced side effects has made cannabinoid (CB) receptors one of the most important G-protein-coupled receptor (GPCR) targets for drug discovery. The cannabinoid (CB) receptor subtype CB2 is of particular interest due to its involvement in signal transduction in the immune system and its increased characterization by mutational and other studies. However, our understanding of their mode of action has been limited by the absence of an experimental receptor structure. In this study, we have developed a 3D model of the CB2 receptor based on the recent crystal structure of a related GPCR, bovine rhodopsin. The model was developed using multiple sequence alignment of homologous receptor sub-types in humans and mammals, and compared with other GPCRs. Alignments were analyzed with mutation scores, pairwise hydrophobicity profiles and Kyte-Doolittle plots. The 3D model of the transmembrane segment was generated by mapping the CB2 sequence onto the homologous residues of the rhodopsin structure. The extra- and intracellular loop regions of the CB2 were generated by searching for homologous C(alpha) backbone sequences in published structures in the Brookhaven Protein Databank (PDB). Residue side chains were positioned through a combination of rotamer library searches, simulated annealing and minimization. Intermediate models of the 7TM helix bundles were analyzed in terms of helix tilt angles, hydrogen-bond networks, conserved residues and motifs, possible disulfide bonds. The amphipathic cytoplasmic helix domain was also correlated with biological and site-directed mutagenesis data. Finally, the model receptor-binding cavity was characterized using solvent-accessible surface approach.     

A transmembrane helix-bundle from G-protein coupled receptor CB2: biosynthesis, purification, and NMR characterization

The cannabinoid receptor subtype 2 (CB2) is a member of the G-protein coupled receptor (GPCR) superfamily. As the relationship between structure and function for this receptor remains poorly understood, the present study was undertaken to characterize the structure of a segment including the first and second transmembrane helix (TM1 and TM2) domains of CB2. To accomplish this, a transmembrane double-helix bundle from this region was expressed, purified, and characterized by NMR. Milligrams of this hydrophobic fragment of the receptor were biosynthesized using a fusion protein overexpression strategy and purified by affinity chromatography combined with reverse phase HPLC. Chemical and enzymatic cleavage methods were implemented to remove the fusion tag. The resultant recombinant protein samples were analyzed and confirmed by HPLC, mass spectrometry, and circular dichroism (CD). The CD analyses of HPLC-purified protein in solution and in DPC micelle preparations suggested predominant alpha-helical structures under both conditions. The 13C/15N double-labeled protein CB2(27-101) was further verified and analyzed by NMR spectroscopy. Sequential assignment was accomplished for more than 80% of residues. The 15N HSQC NMR results show a clear chemical shift dispersion of the amide nitrogen-proton correlation indicative of a pure double-labeled polypeptide molecule. The results suggest that this method is capable of generating transmembrane helical bundles from GPCRs in quantity and purity sufficient for NMR and other biophysical studies. Therefore, the biosynthesis of GPCR transmembrane helix bundles represents a satisfactory alternative strategy to obtain and assemble NMR structures from recombinant "building blocks."     

Trifluoroethanol and binding to model membranes stabilize a predicted turn in a peptide corresponding to the first extracellular loop of the angiotensin II AT(1A) receptor

Deuterium oxide solutions of a triple-helical polysaccharide schizophyllan, undergoing an order-disorder transition centered around 17 degrees C, were studied by the time-domain reflectometry (TDR) to obtain dielectric dispersions in the solution and frozen states. In the solution state, the dispersion below the transition temperature is resolved in three dispersions (relaxation times at 0 degrees C) ascribed to side chain glucose residue (1; 102 ns), structured water (s; 2.0 ns) and bulk water (h), respectively, from low to high frequencies. Bulk water is divided into slow water (h2; 0.04 ns) and free or pure water (h1; 0.02 ns). Above the transition temperature structured water almost disappears and is compensated by slow water. Structured water is similar to bound water for proteins but different from it because of this transition behavior. Another dispersion (l) seen at the lowest frequency is assigned to the rotation of side-chain glucose residue coupled with hydrated water. Parts of this dispersion and structured water are suggested to constitute bound water. In the frozen state were observed a major dispersion (h; 0.14 ns) and a minor one (m; 28 ns), which were ascribed to considerably mobile and less mobile waters. They are similar to but not exactly the same as that for unfreezable water in bovine serum albumin solutions argued by Miura et al. (Biopolymers, 1995, Vol. 36, p. 9). Water is molded into different structures by the triple helix.     

Involvement of the second extracellular loop and transmembrane residues of CCR5 in inhibitor binding and HIV-1 fusion: insights into the mechanism of allosteric inhibition

C-C chemokine receptor 5 (CCR5), a member of G-protein-coupled receptors, serves as a coreceptor for human immunodeficiency virus type 1 (HIV-1). In the present study, we examined the interactions between CCR5 and novel CCR5 inhibitors containing the spirodiketopiperazine scaffolds AK530 and AK317, both of which were lodged in the hydrophobic cavity located between the upper transmembrane domain and the second extracellular loop (ECL2) of CCR5. Although substantial differences existed between the two inhibitors—AK530 had 10-fold-greater CCR5-binding affinity (K_d = 1.4 nM) than AK317 (16.7 nM)—their antiviral potencies were virtually identical (IC₅₀ = 2.1 nM and 1.5 nM, respectively). Molecular dynamics simulations for unbound CCR5 showed hydrogen bond interactions among transmembrane residues Y108, E283, and Y251, which were crucial for HIV-1-gp120/sCD4 complex binding and HIV-1 fusion. Indeed, AK530 and AK317, when bound to CCR5, disrupted these interhelix hydrogen bond interactions, a salient molecular mechanism enabling allosteric inhibition. Mutagenesis and structural analysis showed that ECL2 consists of a part of the hydrophobic cavity for both inhibitors, although AK317 is more tightly engaged with ECL2 than AK530, explaining their similar anti-HIV-1 potencies despite the difference in K_d values. We also found that amino acid residues in the β-hairpin structural motif of ECL2 are critical for HIV-1-elicited fusion and binding of the spirodiketopiperazine-based inhibitors to CCR5. The direct ECL2-engaging property of the inhibitors likely produces an ECL2 conformation, which HIV-1 gp120 cannot bind to, but also prohibits HIV-1 from utilizing the “inhibitor-bound” CCR5 for cellular entry—a mechanism of HIV-1's resistance to CCR5 inhibitors. The data should not only help delineate the dynamics of CCR5 following inhibitor binding but also aid in designing CCR5 inhibitors that are more potent against HIV-1 and prevent or delay the emergence of resistant HIV-1 variants.     

Role of the transmembrane domain 4/extracellular loop 2 junction of the human gonadotropin-releasing hormone receptor in ligand binding and receptor conformational selection

Recent crystal structures of G protein-coupled receptors (GPCRs) show the remarkable structural diversity of extracellular loop 2 (ECL2), implying its potential role in ligand binding and ligand-induced receptor conformational selectivity. Here we have applied molecular modeling and mutagenesis studies to the TM4/ECL2 junction (residues Pro(174(4.59))-Met(180(4.66))) of the human gonadotropin-releasing hormone (GnRH) receptor, which uniquely has one functional type of receptor but two endogenous ligands in humans. We suggest that the above residues assume an α-helical extension of TM4 in which the side chains of Gln(174(4.60)) and Phe(178(4.64)) face toward the central ligand binding pocket to make H-bond and aromatic contacts with pGlu(1) and Trp(3) of both GnRH I and GnRH II, respectively. The interaction between the side chains of Phe(178(4.64)) of the receptor and Trp(3) of the GnRHs was supported by reciprocal mutations of the interacting residues. Interestingly, alanine mutations of Leu(175(4.61)), Ile(177(4.63)), and Met(180(4.66)) decreased mutant receptor affinity for GnRH I but, in contrast, increased affinity for GnRH II. This suggests that these residues make intramolecular or intermolecular contacts with residues of transmembrane (TM) domain 3, TM5, or the phospholipid bilayer, which couple the ligand structure to specific receptor conformational switches. The marked decrease in signaling efficacy of I177A and F178A also indicates that IIe(177(4.63)) and Phe(178(4.64)) are important in stabilizing receptor-active conformations. These findings suggest that the TM4/ECL2 junction is crucial for peptide ligand binding and, consequently, for ligand-induced receptor conformational selection.     

Heteromultimerization of cannabinoid CB(1) receptor and orexin OX(1) receptor generates a unique complex in which both protomers are regulated by orexin A

Agonist-induced internalization was observed for both inducible and constitutively expressed forms of the cannabinoid CB(1) receptor. These were also internalized by the peptide orexin A, which has no direct affinity for the cannabinoid CB(1) receptor, but only when the orexin OX(1) receptor was co-expressed along with the cannabinoid CB(1) receptor. This effect of orexin A was concentration-dependent and blocked by OX(1) receptor antagonists. Moreover, the ability of orexin A to internalize the CB(1) receptor was also blocked by CB(1) receptor antagonists. Remarkably, orexin A was substantially more potent in producing internalization of the CB(1) receptor than in causing internalization of the bulk OX(1) receptor population, and this was true in cells in which the CB(1) receptor was maintained at a constant level, whereas levels of OX(1) could be varied and vice versa. Both co-immunoprecipitation and cell surface, homogenous time-resolved fluorescence resonance energy transfer based on covalent labeling of N-terminal "SNAP" and "CLIP" tags present in the extracellular N-terminal domain of the receptors confirmed the capacity of these two receptors to heteromultimerize. These studies confirm the capacity of the CB(1) and OX(1) receptors to interact directly and demonstrate that this complex has unique regulatory characteristics. The higher potency of the agonist orexin A to regulate the CB(1)-OX(1) heteromer compared with the OX(1)-OX(1) homomer present in the same cells and the effects of CB(1) receptor antagonists on the function of orexin A suggest an interplay between these two systems that may modulate appetite, feeding, and wakefulness.     

Homology modeling and molecular dynamics simulations of the active state of the nociceptin receptor reveal new insights into agonist binding and activation

The opioid receptor-like receptor, also known as the nociceptin receptor (NOP), is a class A G protein-coupled receptor (GPCR) in the opioid receptor family. Although NOP shares a significant homology with the other opioid receptors, it does not bind known opioid ligands and has been shown to have a distinct mechanism of activation compared to the closely related opioid receptors mu, delta, and kappa. Previously reported homology models of the NOP receptor, based on the inactive-state GPCR crystal structures, give limited information on the activation and selectivity features of this fourth member of the opioid receptor family. We report here the first active-state homology model of the NOP receptor based on the opsin GPCR crystal structure. An inactive-state homology model of NOP was also built using a multiple template approach. Molecular dynamics simulation of the active-state NOP model and comparison to the inactive-state model suggest that NOP activation involves movements of transmembrane (TM)3 and TM6 and several activation microswitches, consistent with GPCR activation. Docking of the selective nonpeptidic NOP agonist ligand Ro 64-6198 into the active-state model reveals active-site residues in NOP that play a role in the high selectivity of this ligand for NOP over the other opioid receptors. Docking the shortest active fragment of endogenous agonist nociceptin/orphaninFQ (residues 1-13) shows that the NOP extracellular loop 2 (EL2) loop interacts with the positively charged residues (8-13) of N/OFQ. Both agonists show extensive polar interactions with residues at the extracellular end of the TM domain and EL2 loop, suggesting agonist-induced reorganization of polar networks, during receptor activation.     

Second extracellular loop of human glucagon-like peptide-1 receptor (GLP-1R) has a critical role in GLP-1 peptide binding and receptor activation

The glucagon-like peptide-1 receptor (GLP-1R) is a therapeutically important family B G protein-coupled receptor (GPCR) that is pleiotropically coupled to multiple signaling effectors and, with actions including regulation of insulin biosynthesis and secretion, is one of the key targets in the management of type II diabetes mellitus. However, there is limited understanding of the role of the receptor core in orthosteric ligand binding and biological activity. To assess involvement of the extracellular loop (ECL) 2 in ligand-receptor interactions and receptor activation, we performed alanine scanning mutagenesis of loop residues and assessed the impact on receptor expression and GLP-1(1-36)-NH(2) or GLP-1(7-36)-NH(2) binding and activation of three physiologically relevant signaling pathways as follows: cAMP formation, intracellular Ca(2+) (Ca(2+)(i)) mobilization, and phosphorylation of extracellular signal-regulated kinases 1 and 2 (pERK1/2). Although antagonist peptide binding was unaltered, almost all mutations affected GLP-1 peptide agonist binding and/or coupling efficacy, indicating an important role in receptor activation. However, mutation of several residues displayed distinct pathway responses with respect to wild type receptor, including Arg-299 and Tyr-305, where mutation significantly enhanced both GLP-1(1-36)-NH(2)- and GLP-1(7-36)-NH(2)-mediated signaling bias for pERK1/2. In addition, mutation of Cys-296, Trp-297, Asn-300, Asn-302, and Leu-307 significantly increased GLP-1(7-36)-NH(2)-mediated signaling bias toward pERK1/2. Of all mutants studied, only mutation of Trp-306 to alanine abolished all biological activity. These data suggest a critical role of ECL2 of the GLP-1R in the activation transition(s) of the receptor and the importance of this region in the determination of both GLP-1 peptide- and pathway-specific effects.     

Second extracellular loop of human glucagon-like peptide-1 receptor (GLP-1R) differentially regulates orthosteric but not allosteric agonist binding and function

The glucagon-like peptide-1 receptor (GLP-1R) is a prototypical family B G protein-coupled receptor that exhibits physiologically important pleiotropic coupling and ligand-dependent signal bias. In our accompanying article (Koole, C., Wootten, D., Simms, J., Miller, L. J., Christopoulos, A., and Sexton, P. M. (2012) J. Biol. Chem. 287, 3642-3658), we demonstrate, through alanine-scanning mutagenesis, a key role for extracellular loop (ECL) 2 of the receptor in propagating activation transition mediated by GLP-1 peptides that occurs in a peptide- and pathway-dependent manner for cAMP formation, intracellular (Ca(2+)(i)) mobilization, and phosphorylation of extracellular signal-regulated kinases 1 and 2 (pERK1/2). In this study, we examine the effect of ECL2 mutations on the binding and signaling of the peptide mimetics, exendin-4 and oxyntomodulin, as well as small molecule allosteric agonist 6,7-dichloro-2-methylsulfonyl-3-tert-butylaminoquinoxaline (compound 2). Lys-288, Cys-296, Trp-297, and Asn-300 were globally important for peptide signaling and also had critical roles in governing signal bias of the receptor. Peptide-specific effects on relative efficacy and signal bias were most commonly observed for residues 301-305, although R299A mutation also caused significantly different effects for individual peptides. Met-303 was more important for exendin-4 and oxyntomodulin action than those of GLP-1 peptides. Globally, ECL2 mutation was more detrimental to exendin-4-mediated Ca(2+)i release than GLP-1(7-36)-NH(2), providing additional evidence for subtle differences in receptor activation by these two peptides. Unlike peptide activation of the GLP-1R, ECL2 mutations had only limited impact on compound 2 mediated cAMP and pERK responses, consistent with this ligand having a distinct mechanism for receptor activation. These data suggest a critical role of ECL2 of the GLP-1R in the activation transition of the receptor by peptide agonists.     

Computationally‐predicted CB1 cannabinoid receptor mutants show distinct patterns of salt‐bridges that correlate with their level of constitutive activity reflected in G protein coupling levels,thermal stability,and ligand binding

The cannabinoid receptor 1 (CB1), a member of the class A G‐protein‐coupled receptor (GPCR) family, possesses an observable level of constitutive activity. Its activation mechanism, however, has yet to be elucidated. Previously we discovered dramatic changes in CB1 activity due to single mutations; T3.46A, which made the receptor inactive, and T3.46I and L3.43A, which made it essentially fully constitutively active. Our subsequent prediction of the structures of these mutant receptors indicated that these changes in activity are explained in terms of the pattern of salt‐bridges in the receptor region involving transmembrane domains 2, 3, 5, and 6. Here we identified key salt‐bridges, R2.37 + D6.30 and D2.63 + K3.28, critical for CB1 inactive and active states, respectively, and generated new mutant receptors that we predicted would change CB1 activity by either precluding or promoting these interactions. We find that breaking the R2.37 + D6.30 salt‐bridge resulted in substantial increase in G‐protein coupling activity and reduced thermal stability relative to the wild‐type reflecting the changes in constitutive activity from inactive to active. In contrast, breaking the D2.63 + K3.28 salt‐bridge produced the opposite profile suggesting this interaction is critical for the receptor activation. Thus, we demonstrate an excellent correlation with the predicted pattern of key salt‐bridges and experimental levels of activity and conformational flexibility. These results are also consistent with the extended ternary complex model with respect to shifts in agonist and inverse agonist affinity and provide a powerful framework for understanding the molecular basis for the multiple stages of CB1 activation and that of other GPCRs in general. Proteins 2013; 81:1304–1317. © 2013 Wiley Periodicals, Inc.     

Understanding interactions of gastric inhibitory polypeptide (GIP) with its G-protein coupled receptor through NMR and molecular modeling.

Gastric inhibitory polypeptide (GIP, or glucose-dependent insulinotropic polypeptide) is a 42-amino acid incretin hormone moderating glucose-induced insulin secretion. Antidiabetic therapy based on GIP holds great promise because of the fact that its insulinotropic action is highly dependent on the level of glucose, overcoming the sideeffects of hypoglycemia associated with the current therapy of Type 2 diabetes. The truncated peptide, GIP(1-30)NH2, has the same activity as the full length native peptide. We have studied the structure of GIP(1-30)NH2 and built a model of its G-protein coupled receptor (GPCR). The structure of GIP(1-30)NH2 in DMSO-d6 and H2O has been studied using 2D NMR (total correlation spectroscopy (TOCSY), nuclear overhauser effect spectroscopy (NOESY), double quantum filtered-COSY (DQF-COSY), 13C-heteronuclear single quantum correlation (HSQC) experiments, and its conformation built by MD simulations with the NMR data as constraints. The peptide in DMSO-d6 exhibits an alpha-helix between residues Ile12 and Lys30 with a discontinuity at residues Gln19 and Gln20. In H2O, the alpha-helix starts at Ile7, breaks off at Gln19, and then continues right through to Lys30. GIP(1-30)NH2 has all the structural features of peptides belonging to family B1 GPCRs, which are characterized by a coil at the N-terminal and a long C-terminal alpha-helix with or without a break. A model of the seven transmembrane (TM) helices of the GIP receptor (GIPR) has been built on the principles of comparative protein modeling, using the crystal structure of bovine rhodopsin as a template. The N-terminal domain of GIPR has been constructed from the NMR structure of the N-terminal of corticoptropin releasing factor receptor (CRFR), a family B1 GCPR. The intra and extra cellular loops and the C-terminal have been modeled from fragments retrieved from the PDB. On the basis of the experimental data available for some members of family B1 GPCRs, four pairs of constraints between GIP(1-30)NH2 and its receptor were used in the FTDOCK program, to build the complete model of the GIP(1-30)NH2:GIPR complex. The model can rationalize the various experimental observations including the potency of the truncated GIP peptide. This work is the first complete model at the atomic level of GIP(1-30)NH2 and of the complex with its GPCR.     

Investigation on the binding mechanism of loratinib with the c-ros oncogene 1 (ROS1) receptor tyrosine kinase via molecular dynamics simulation and binding free energy calculations

The c-ros oncogene 1 (ROS1) has proven to be an important cancer target for the treatment of various human cancers. The anaplastic lymphoma kinase inhibitor crizotinib has been granted approval for the treatment of patients with ROS1 positive metastatic non-small-cell lung cancer by the Food and Drug Administration on 2016. However, serious resistance due to the secondary mutation of glycine 2032 to arginine (G2032R) was developed in clinical studies. Loratinib (PF-06463922), a macrocyclic analog of crizotinib, showed significantly improved inhibitory activity against wild–type (WT) ROS1 and ROS1^G2032R mutant. To provide insights into the inhibition mechanism, molecular dynamics simulations and free energy calculations were carried out for the complexes of loratinib with WT and G2032R mutated ROS1. The apo-ROS1^WT and apo-ROS1^G2032R systems showed similar RMSF distributions, while ROS1^G2032R-loratinib showed significantly higher than that of WT ROS1-loratinib, which revealed that the binding of loratinib to ROS1^G2032R significantly interfered the ?uctuation of protein. Calculations of binding free energies indicate that G2032R mutation significantly reduces the binding affinity of loratinib for ROS1, which arose mostly from the increase of conformation entropy and the decrease of solvation energy. Furthermore, detailed per-residue binding free energies highlighted the increased and decreased contributions of some residues in the G2032R mutated systems. The present study revealed the detailed inhibitory mechanism of loratinib as potent WT and G2032R mutated ROS1 inhibitor, which was expected to provide a basis for rational drug design.     

Computer modeling and molecular dynamics simulations of ligand bound complexes of bovine angiogenin: dinucleotide topology at the active site of RNase a family proteins

We have undertaken the modeling of substrate-bound structures of angiogenin. In our recent study, we modeled the dinucleotide ligand binding to human angiogenin. In the present study, the substrates CpG, UpG, and CpA were docked onto bovine angiogenin. This was achieved by overcoming the problem of an obstruction to the B1 site by the C-terminus and identifying residues that bind to the second base. The modeled complexes retain biochemically important interactions. The docked models were subjected to 1 ns of molecular dynamics, and structures from the simulation were refined by using simulated annealing. Our models explained the enzyme's specificity for both B1 and B2 bases as observed experimentally. The nature of binding of the dinucleotide substrate was compared with that of the mononucleotide product. The models of these complexes were also compared with those obtained earlier with human angiogenin. On the basis of the simulations and annealed structures, we came up with a consensus topology of dinucleotide ligands that binds to human and bovine angiogenins. This dinucleotide conformation can serve as a starting model for ligand-bound complex structures for RNase A family of proteins. We demonstrated this capability by generating the complex structure of CpA bound to eosinophil-derived neurotoxin (EDN) by fitting the consensus topology of CpA to the crystal structure of native EDN.     

Extracellular loops 1 and 3 and their associated transmembrane regions of the calcitonin receptor-like receptor are needed for CGRP receptor function

The first and third extracellular loops (ECL) of G protein-coupled receptors (GPCRs) have been implicated in ligand binding and receptor function. This study describes the results of an alanine/leucine scan of ECLs 1 and 3 and loop-associated transmembrane (TM) domains of the secretin-like GPCR calcitonin receptor-like receptor which associates with receptor activity modifying protein 1 to form the CGRP receptor. Leu195Ala, Val198Ala and Ala199Leu at the top of TM2 all reduced αCGRP-mediated cAMP production and internalization; Leu195Ala and Ala199Leu also reduced αCGRP binding. These residues form a hydrophobic cluster within an area defined as the “minor groove” of rhodopsin-like GPCRs. Within ECL1, Ala203Leu and Ala206Leu influenced the ability of αCGRP to stimulate adenylate cyclase. In TM3, His219Ala, Leu220Ala and Leu222Ala have influences on αCGRP binding and cAMP production; they are likely to indirectly influence the binding site for αCGRP as well as having an involvement in signal transduction. On the exofacial surfaces of TMs 6 and 7, a number of residues were identified that reduced cell surface receptor expression, most noticeably Leu351Ala and Glu357Ala in TM6. The residues may contribute to the RAMP1 binding interface. Ile360Ala impaired αCGRP-mediated cAMP production. Ile360 is predicted to be located close to ECL2 and may facilitate receptor activation. Identification of several crucial functional loci gives further insight into the activation mechanism of this complex receptor system and may aid rational drug design.     

New insights into the GABA(A) receptor structure and orthosteric ligand binding: receptor modeling guided by experimental data

GABA(A) receptors (GABA(A)Rs) are ligand gated chloride ion channels that mediate overall inhibitory signaling in the CNS. A detailed understanding of their structure is important to gain insights in, e.g., ligand binding and functional properties of this pharmaceutically important target. Homology modeling is a necessary tool in this regard because experimentally determined structures are lacking. Here we present an exhaustive approach for creating a high quality model of the α(1)β(2)γ(2) subtype of the GABA(A)R ligand binding domain, and we demonstrate its usefulness in understanding details of orthosteric ligand binding. The model was constructed by using multiple templates and by incorporation of knowledge from biochemical/pharmacological experiments. It was validated on the basis of objective energy functions, its ability to account for available residue specific information, and its stability in molecular dynamics (MD) compared with that of the two homologous crystal structures. We then combined the model with extensive structure-activity relationships available from two homologous series of orthosteric GABA(A)R antagonists to create a detailed hypothesis for their binding modes. Excellent agreement with key experimental data was found, including the ability of the model to accommodate and explain a previously developed pharmacophore model. A coupling to agonist binding was thereby established and discussed in relation to activation mechanisms. Our results highlight the importance of critical evaluation and optimization of each step in the homology modeling process. The approach taken here can greatly aid in increasing the understanding of GABA(A)Rs and related receptors where structural insight is limited and reliable models are difficult to obtain.     

The conformation of the extracellular binding domain of Death Receptor 5 in the presence and absence of the activating ligand TRAIL: a molecular dynamics study

The Death Receptor 5 (DR5), a member of tumor necrosis factor receptor (TNFR) superfamily of receptors, triggers apoptosis (programmed cell death) when stimulated by its tridentate ligand TRAIL. Until recently it was generally assumed that the activation of DR5 resulted from the recruitment of three independent receptor units, leading to the trimerization of intracellular domains. However, there is mounting evidence to suggest that, in the absence of ligand, such cytokine receptors primarily reside as preformed complexes. In this work, molecular dynamics simulations of the TRAIL-DR5 complex, the unbound receptor trimer and individual receptor monomers are compared to gain insight in the mechanism of activation. The results suggest that, in the absence of TRAIL, DR5 has a strong propensity to self-associate and that this is primarily mediated through interactions of the membrane proximal domains. The association of the free receptors leads to a loss of the threefold symmetry found within the receptor-ligand complex. The simulations suggest that the primary role of TRAIL is to induce threefold-symmetry within the DR5 complex and to constrain the receptor to a specific conformation. The implications of this in terms of the mechanism by which the receptor switches from an inactive to an active state are discussed.     

Identification of Distinct Conformations of the Angiotensin-II Type 1 Receptor Associated with the Gq/11 Protein Pathway and the β-Arrestin Pathway Using Molecular Dynamics Simulations

Biased signaling represents the ability of G protein-coupled receptors to engage distinct pathways with various efficacies depending on the ligand used or on mutations in the receptor. The angiotensin-II type 1 (AT₁) receptor, a prototypical class A G protein-coupled receptor, can activate various effectors upon stimulation with the endogenous ligand angiotensin-II (AngII), including the G_q/11 protein and β-arrestins. It is believed that the activation of those two pathways can be associated with distinct conformations of the AT₁ receptor. To verify this hypothesis, microseconds of molecular dynamics simulations were computed to explore the conformational landscape sampled by the WT-AT₁ receptor, the N111G-AT₁ receptor (constitutively active and biased for the G_q/11 pathway), and the D74N-AT₁ receptor (biased for the β-arrestin1 and -2 pathways) in their apo-forms and in complex with AngII. The molecular dynamics simulations of the AngII-WT-AT₁, N111G-AT₁, and AngII-N111G-AT₁ receptors revealed specific structural rearrangements compared with the initial and ground state of the receptor. Simulations of the D74N-AT₁ receptor revealed that the mutation stabilizes the receptor in the initial ground state. The presence of AngII further stabilized the ground state of the D74N-AT₁ receptor. The biased agonist [Sar¹,Ile⁸]AngII also showed a preference for the ground state of the WT-AT₁ receptor compared with AngII. These results suggest that activation of the G_q/11 pathway is associated with a specific conformational transition stabilized by the agonist, whereas the activation of the β-arrestin pathway is linked to the stabilization of the ground state of the receptor.