Determinants of Orexin Receptor Binding and Activation-A Molecular Dynamics Study.
Ontology highlight
ABSTRACT:
SUBMITTER: Karhu L
PROVIDER: S-EPMC6727383 | biostudies-literature | 2019 Mar
REPOSITORIES: biostudies-literature
Publications
Determinants of Orexin Receptor Binding and Activation-A Molecular Dynamics Study.
Karhu Lasse L Magarkar Aniket A Bunker Alex A Xhaard Henri H
The journal of physical chemistry. B 20190318 12
We assess the stability of two previously suggested binding modes for the neuropeptide orexin-A in the OX<sub>2</sub> receptor through extensive molecular dynamics simulations. As the activation determinants of the receptor remain unknown, we simulated an unliganded receptor and two small-molecular ligands, the antagonist suvorexant and the agonist Nag26 for comparison. Each system was simulated in pure POPC membrane as well as in the 25% cholesterol-POPC membrane. In total, we carried out 36 μs ...[more]
Similar Datasets
Project description:The epidermal growth factor receptor (EGFR) is a member of the receptor tyrosine kinase family that plays a role in multiple cellular processes. Activation of EGFR requires binding of a ligand on the extracellular domain to promote conformational changes leading to dimerization and transphosphorylation of intracellular kinase domains. Seven ligands are known to bind EGFR with affinities ranging from sub-nanomolar to near micromolar dissociation constants. In the case of EGFR, distinct conformational states assumed upon binding a ligand is thought to be a determining factor in activation of a downstream signaling network. Previous biochemical studies suggest the existence of both low affinity and high affinity EGFR ligands. While these studies have identified functional effects of ligand binding, high-resolution structural data are lacking. To gain a better understanding of the molecular basis of EGFR binding affinities, we docked each EGFR ligand to the putative active state extracellular domain dimer and 25.0 ns molecular dynamics simulations were performed. MM-PBSA/GBSA are efficient computational approaches to approximate free energies of protein-protein interactions and decompose the free energy at the amino acid level. We applied these methods to the last 6.0 ns of each ligand-receptor simulation. MM-PBSA calculations were able to successfully rank all seven of the EGFR ligands based on the two affinity classes: EGF>HB-EGF>TGF-?>BTC>EPR>EPG>AR. Results from energy decomposition identified several interactions that are common among binding ligands. These findings reveal that while several residues are conserved among the EGFR ligand family, no single set of residues determines the affinity class. Instead we found heterogeneous sets of interactions that were driven primarily by electrostatic and Van der Waals forces. These results not only illustrate the complexity of EGFR dynamics but also pave the way for structure-based design of therapeutics targeting EGF ligands or the receptor itself.
| S-EPMC3554757 | biostudies-literature
Project description:The endoplasmic reticulum (ER) Ca2+ sensors stromal interaction molecule 1 (STIM1) and STIM2, which connect ER Ca2+ depletion with extracellular Ca2+ influx, are crucial for the maintenance of Ca2+ homeostasis in mammalian cells. Despite the recent progress in unraveling the role of STIM2 in Ca2+ signaling, the mechanistic underpinnings of its activation remain underexplored. We use an engineering approach to direct ER-resident STIMs to the plasma membrane (PM) while maintaining their correct membrane topology, as well as Förster resonance energy transfer (FRET) sensors that enabled in cellulo real-time monitoring of STIM activities. This allowed us to determine the calcium affinities of STIM1 and STIM2 both in cellulo and in situ, explaining the current discrepancies in the literature. We also identified the key structural determinants, especially the corresponding G residue in STIM1, which define the distinct activation dynamics of STIM2. The chimeric E470G mutation could switch STIM2 from a slow and weak Orai channel activator into a fast and potent one like STIM1 and vice versa. The systemic dissection of STIM2 activation by protein engineering sets the stage for the elucidation of the regulation and function of STIM2-mediated signaling in mammals.
| S-EPMC6267984 | biostudies-literature
Project description:RDL receptors are GABA-activated inhibitory Cys-loop receptors found throughout the insect CNS. They are a key target for insecticides. Here, we characterize the GABA binding site in RDL receptors using computational and electrophysiological techniques. A homology model of the extracellular domain of RDL was generated and GABA docked into the binding site. Molecular dynamics simulations predicted critical GABA binding interactions with aromatic residues F206, Y254, and Y109 and hydrophilic residues E204, S176, R111, R166, S176, and T251. These residues were mutated, expressed in Xenopus oocytes, and their functions assessed using electrophysiology. The data support the binding mechanism provided by the simulations, which predict that GABA forms many interactions with binding site residues, the most significant of which are cation-? interactions with F206 and Y254, H-bonds with E204, S205, R111, S176, T251, and ionic interactions with R111 and E204. These findings clarify the roles of a range of residues in binding GABA in the RDL receptor, and also show that molecular dynamics simulations are a useful tool to identify specific interactions in Cys-loop receptors.
| S-EPMC3512037 | biostudies-literature
Project description:Neurosteroids positively modulate GABA-A receptor (GABAAR) channel activity by binding to a transmembrane domain intersubunit site. Understanding the interactions in this site that determine neurosteroid binding and its effect is essential for the design of neurosteroid-based therapeutics. Using photo-affinity labeling and an ELIC-?1GABAAR chimera, we investigated the impact of mutations (Q242L, Q242W and W246L) within the intersubunit site on neurosteroid binding. These mutations, which abolish the thermal stabilizing effect of allopregnanolone on the chimera, reduce neither photolabeling within the intersubunit site nor competitive prevention of labeling by allopregnanolone. Instead, these mutations change the orientation of neurosteroid photolabeling. Molecular docking of allopregnanolone in WT and Q242W receptors confirms that the mutation favors re-orientation of allopregnanolone within the binding pocket. Collectively, the data indicate that mutations at Gln242 or Trp246 that eliminate neurosteroid effects do not eliminate neurosteroid binding within the intersubunit site, but significantly alter the preferred orientation of the neurosteroid within the site. The interactions formed by Gln242 and Trp246 within this pocket play a vital role in determining the orientation of the neurosteroid that may be necessary for its functional effect.
| S-EPMC6708749 | biostudies-literature
Project description:The adenosine A2A receptor (A(2A)R) plays a key role in transmembrane signaling mediated by the endogenous agonist adenosine. Here, we describe the crystal structure of human A2AR thermostabilized in an active-like conformation bound to the selective agonist 2-[p-(2-carboxyethyl)phenylethyl-amino]-5'-N-ethylcarboxamido adenosine (CGS21680) at a resolution of 2.6 Å. Comparison of A(2A)R structures bound to either CGS21680, 5'-N-ethylcarboxamido adenosine (NECA), UK432097 [6-(2,2-diphenylethylamino)-9-[(2R,3R,4S,5S)-5-(ethylcarbamoyl)-3,4-dihydroxy-tetrahydrofuran-2-yl]-N-[2-[[1-(2-pyridyl)-4-piperidyl]carbamoylamino]ethyl]purine-2-carboxamide], or adenosine shows that the adenosine moiety of the ligands binds to the receptor in an identical fashion. However, an extension in CGS21680 compared with adenosine, the (2-carboxyethyl)phenylethylamino group, binds in an extended vestibule formed from transmembrane regions 2 and 7 (TM2 and TM7) and extracellular loops 2 and 3 (EL2 and EL3). The (2-carboxyethyl)phenylethylamino group makes van der Waals contacts with side chains of amino acid residues Glu169(EL2), His264(EL3), Leu267(7.32), and Ile274(7.39), and the amine group forms a hydrogen bond with the side chain of Ser67(2.65). Of these residues, only Ile274(7.39) is absolutely conserved across the human adenosine receptor subfamily. The major difference between the structures of A(2A)R bound to either adenosine or CGS21680 is that the binding pocket narrows at the extracellular surface when CGS21680 is bound, due to an inward tilt of TM2 in that region. This conformation is stabilized by hydrogen bonds formed by the side chain of Ser67(2.65) to CGS21680, either directly or via an ordered water molecule. Mutation of amino acid residues Ser67(2.65), Glu169(EL2), and His264(EL3), and analysis of receptor activation either in the presence or absence of ligands implicates this region in modulating the level of basal activity of A(2A)R.
| S-EPMC4720118 | biostudies-literature
Project description:The histamine H2 receptor (H2R) plays an important role in the regulation of gastric acid secretion. Therefore, it is a main drug target for the treatment of gastroesophageal reflux or peptic ulcer disease. However, there is as of yet no 3D-structural information available hampering a mechanistic understanding of H2R. Therefore, we created a model of the histamine-H2R-Gs complex based on the structure of the ternary complex of the β2-adrenoceptor and investigated the conformational stability of this active GPCR conformation. Since the physiologically relevant motions with respect to ligand binding and conformational changes of GPCRs can only partly be assessed on the timescale of conventional MD (cMD) simulations, we also applied metadynamics and Gaussian accelerated molecular dynamics (GaMD) simulations. A multiple walker metadynamics simulation in combination with cMD was applied for the determination of the histamine binding mode. The preferential binding pose detected is in good agreement with previous data from site directed mutagenesis and provides a basis for rational ligand design. Inspection of the H2R-Gs interface reveals a network of polar interactions that may contribute to H2R coupling selectivity. The cMD and GaMD simulations demonstrate that the active conformation is retained on a μs-timescale in the ternary histamine-H2R-Gs complex and in a truncated complex that contains only Gs helix α5 instead of the entire G protein. In contrast, histamine alone is unable to stabilize the active conformation, which is in line with previous studies of other GPCRs.
| S-EPMC7554837 | biostudies-literature
Project description:In many biological applications, we would like to be able to computationally predict mutational effects on affinity in protein-protein interactions. However, many commonly used methods to predict these effects perform poorly in important test cases. In particular, the effects of multiple mutations, non alanine substitutions, and flexible loops are difficult to predict with available tools and protocols. We present here an existing method applied in a novel way to a new test case; we interrogate affinity differences resulting from mutations in a host-virus protein-protein interface. We use steered molecular dynamics (SMD) to computationally pull the machupo virus (MACV) spike glycoprotein (GP1) away from the human transferrin receptor (hTfR1). We then approximate affinity using the maximum applied force of separation and the area under the force-versus-distance curve. We find, even without the rigor and planning required for free energy calculations, that these quantities can provide novel biophysical insight into the GP1/hTfR1 interaction. First, with no prior knowledge of the system we can differentiate among wild type and mutant complexes. Moreover, we show that this simple SMD scheme correlates well with relative free energy differences computed via free energy perturbation. Second, although the static co-crystal structure shows two large hydrogen-bonding networks in the GP1/hTfR1 interface, our simulations indicate that one of them may not be important for tight binding. Third, one viral site known to be critical for infection may mark an important evolutionary suppressor site for infection-resistant hTfR1 mutants. Finally, our approach provides a framework to compare the effects of multiple mutations, individually and jointly, on protein-protein interactions.
| S-EPMC3940602 | biostudies-literature
Project description:Orexin receptor antagonism represents a novel approach for the treatment of insomnia that directly targets sleep/wake regulation. Several such compounds have entered into clinical development, including the dual orexin receptor antagonists, suvorexant and almorexant. In this study, we have used equilibrium and kinetic binding studies with the orexin-2 (OX?) selective antagonist radioligand, [³H]-EMPA, to profile several orexin receptor antagonists. Furthermore, selected compounds were studied in cell-based assays of inositol phosphate accumulation and ERK-1/2 phosphorylation in CHO cells stably expressing the OX2 receptor that employ different agonist incubation times (30 and 5?min, respectively). EMPA, suvorexant, almorexant and TCS-OX-29 all bind to the OX? receptor with moderate to high affinity (pk(I) values ? 7.5), whereas the primarily OX1 selective antagonists SB-334867 and SB-408124 displayed low affinity (pK(I) values ca. 6). Competition kinetic analysis showed that the compounds displayed a range of dissociation rates from very fast (TCS-OX2-29, k(off) = 0.22?min?¹) to very slow (almorexant, k(off) = 0.005?min?¹). Notably, there was a clear correlation between association rate and affinity. In the cell-based assays, fast-offset antagonists EMPA and TCS-OX2-29 displayed surmountable antagonism of orexin-A agonist activity. However, both suvorexant and particularly almorexant cause concentration-dependent depression in the maximal orexin-A response, a profile that is more evident with a shorter agonist incubation time. Analysis according to a hemi-equilibrium model suggests that antagonist dissociation is slower in a cellular system than in membrane binding; under these conditions, almorexant effectively acts as a pseudo-irreversible antagonist.
| S-EPMC3904256 | biostudies-literature
Project description:Cells can sense mechanical force in regulating focal adhesion assembly. One vivid example is the force-induced recruitment of vinculin to reinforce initial contacts between a cell and the extracellular matrix. Crystal structures of the unbound proteins and bound complex between the vinculin head subdomain (Vh1) and the talin vinculin binding site 1 (VBS1) indicate that vinculin undergoes a conformational change upon binding to talin. However, the molecular basis for this event and the precise nature of the binding pathway remain elusive. In this article, molecular dynamics is used to investigate the binding mechanism of Vh1 and VBS1 under minimal constraints to facilitate binding. One simulation demonstrates binding of the two molecules in the complete absence of external force. VBS1 makes early hydrophobic contact with Vh1 by positioning the critical hydrophobic residues (L608, L615, and L622) in the groove formed by helices 1 and 2 of Vh1. The solvent-exposed hydrophobic residues (V619 and L623) then gradually penetrate the hydrophobic core of Vh1, thus further separating helix 1 from helix 2. These critical residues are highly conserved as large hydrophobic side groups in other vinculin binding sites; studies also have demonstrated that these residues are essential in Vh1-VBS1 binding. Similar binding mechanisms are also demonstrated in separate molecular dynamics simulations of Vh1 binding to other vinculin binding sites both in talin and alpha-actinin.
| S-EPMC2483755 | biostudies-literature
Project description:Background and purposeOne of the major responses upon orexin receptor activation is Ca(2+) influx, and this influx seems to amplify the other responses mediated by orexin receptors. However, the reduction in Ca(2+) , often used to assess the importance of Ca(2+) influx, might affect other properties, like ligand-receptor interactions, as suggested for some GPCR systems. Hence, we investigated the role of the ligand-receptor interaction and Ca(2+) signal cascades in the apparent Ca(2+) requirement of orexin-A signalling.Experimental approachReceptor binding was assessed in CHO cells expressing human OX1 receptors with [(125) I]-orexin-A by conventional ligand binding as well as scintillation proximity assays. PLC activity was determined by chromatography.Key resultsBoth orexin receptor binding and PLC activation were strongly dependent on the extracellular Ca(2+) concentration. The relationship between Ca(2+) concentration and receptor binding was the same as that for PLC activation. However, when Ca(2+) entry was reduced by depolarizing the cells or by inhibiting the receptor-operated Ca(2+) channels, orexin-A-stimulated PLC activity was much more strongly inhibited than orexin-A binding.Conclusions and implicationsCa(2+) plays a dual role in orexin signalling by being a prerequisite for both ligand-receptor interaction and amplifying orexin signals via Ca(2+) influx. Some previous results obtained utilizing Ca(2+) chelators have to be re-evaluated based on the results of the current study. From a drug discovery perspective, further experiments need to identify the target for Ca(2+) in orexin-A-OX1 receptor interaction and its mechanism of action.
| S-EPMC4290719 | biostudies-literature