Lessons from free energy simulations of delta-opioid receptor homodimers involving the fourth transmembrane helix.
Ontology highlight
ABSTRACT:
SUBMITTER: Provasi D
PROVIDER: S-EPMC2914489 | biostudies-literature | 2010 Aug
REPOSITORIES: biostudies-literature
Publications
Lessons from free energy simulations of delta-opioid receptor homodimers involving the fourth transmembrane helix.
Provasi Davide D Johnston Jennifer M JM Filizola Marta M
Biochemistry 20100801 31
Several G protein-coupled receptors (GPCRs), including opioid receptors deltaOR, muOR, and kappaOR, have been reported to form stable dimers or oligomers in lipid bilayers and cell membranes. This notion has been recently challenged by imaging data supporting a transient nature of GPCR association. Here we use umbrella sampling reconstructed free energies of deltaOR homodimers involving the fourth transmembrane helix to predict their association constant. The results of these simulations, combin ...[more]
Similar Datasets
Project description:Opioid receptors, like other members of the G protein-coupled receptor (GPCR) family, have been shown to associate to form dimers and/or oligomers at the plasma membrane. Whether this association is stable or transient is not known. Recent compelling evidence suggests that at least some GPCRs rapidly associate and dissociate. We have recently calculated binding affinities from free energy estimates to predict transient association between mouse delta opioid receptor (DOR) protomers at a symmetric interface involving the fourth transmembrane (TM4) helix (herein termed "4" dimer). Here we present disulfide cross-linking experiments with DOR constructs with cysteines substituted at the extracellular ends of TM4 or TM5 that confirm the formation of DOR complexes involving these helices. Our results are consistent with the involvement of TM4 and/or TM5 at the DOR homodimer interface, but possibly with differing association propensities. Coarse-grained (CG) well-tempered metadynamics simulations of two different dimeric arrangements of DOR involving TM4 alone or with TM5 (herein termed "4/5" dimer) in an explicit lipid-water environment confirmed the presence of two structurally and energetically similar configurations of the 4 dimer, as previously assessed by umbrella sampling calculations, and revealed a single energetic minimum of the 4/5 dimer. Additional CG umbrella sampling simulations of the 4/5 dimer indicated that the strength of association between DOR protomers varies depending on the protein region at the interface, with the 4 dimer being more stable than the 4/5 dimer.
| S-EPMC3050604 | biostudies-literature
Project description:Opioid receptors (ORs) have been observed as homo- and heterodimers, but it is unclear if the dimers are stable under physiological conditions, and whether monomers or dimers comprise the predominant fraction in a cell. Here, we use three live-cell imaging approaches to assess dimerization of ORs at expression levels that are 10-100 × smaller than in classical biochemical assays. At membrane densities around 25/µm2, a split-GFP assay reveals that κOR dimerizes, while µOR and δOR stay monomeric. At receptor densities < 5/µm2, single-molecule imaging showed no κOR dimers, supporting the concept that dimer formation depends on receptor membrane density. To directly observe the transition from monomers to dimers, we used a single-molecule assay to assess membrane protein interactions at densities up to 100 × higher than conventional single-molecule imaging. We observe that κOR is monomeric at densities < 10/µm2 and forms dimers at densities that are considered physiological. In contrast, µOR and δOR stay monomeric even at the highest densities covered by our approach. The observation of long-lasting co-localization of red and green κOR spots suggests that it is a specific effect based on OR dimerization and not an artefact of coincidental encounters.
| S-EPMC8629795 | biostudies-literature
Project description:The gating of pentameric ligand-gated ion channels is sensitive to a variety of allosteric modulators that act on structures peripheral to those involved in the allosteric pathway leading from the agonist site to the channel gate. One such structure, the lipid-exposed transmembrane α helix, M4, is the target of lipids, neurosteroids, and disease-causing mutations. Here we show that M4 interactions with the adjacent transmembrane α helices, M1 and M3, modulate pLGIC function. Enhanced M4 interactions promote channel function while ineffective interactions reduce channel function. The interface chemistry governs the intrinsic strength of M4-M1/M3 inter-helical interactions, both influencing channel gating and imparting distinct susceptibilities to the potentiating effects of a lipid-facing M4 congenital myasthenic syndrome mutation. Through aromatic substitutions, functional studies, and molecular dynamics simulations, we elucidate a mechanism by which M4 modulates channel function.
| S-EPMC4824752 | biostudies-literature
Project description:Protein-lipid interaction and bilayer regulation of membrane protein functions are largely controlled by the hydrophobic match between the transmembrane (TM) domain of membrane proteins and the surrounding lipid bilayer. To systematically characterize responses of a TM helix and lipid adaptations to a hydrophobic mismatch, we have performed a total of 5.8-mus umbrella sampling simulations and calculated the potentials of mean force (PMFs) as a function of TM helix tilt angle under various mismatch conditions. Single-pass TM peptides called WALPn (n = 16, 19, 23, and 27) were used in two lipid bilayers with different hydrophobic thicknesses to consider hydrophobic mismatch caused by either the TM length or the bilayer thickness. In addition, different flanking residues, such as alanine, lysine, and arginine, instead of tryptophan in WALP23 were used to examine their influence. The PMFs, their decomposition, and trajectory analysis demonstrate that 1), tilting of a single-pass TM helix is the major response to a hydrophobic mismatch; 2), TM helix tilting up to approximately 10 degrees is inherent due to the intrinsic entropic contribution arising from helix precession around the membrane normal even under a negative mismatch; 3), the favorable helix-lipid interaction provides additional driving forces for TM helix tilting under a positive mismatch; 4), the minimum-PMF tilt angle is generally located where there is the hydrophobic match and little lipid perturbation; 5), TM helix rotation is dependent on the specific helix-lipid interaction; and 6), anchoring residues at the hydrophilic/hydrophobic interface can be an important determinant of TM helix orientation.
| S-EPMC2895360 | biostudies-literature
Project description:Orientation, dynamics, and packing of transmembrane helical peptides are important determinants of membrane protein structure, dynamics, and function. Because it is difficult to investigate these aspects by studying real membrane proteins, model transmembrane helical peptides are widely used. NMR experiments provide information on both orientation and dynamics of peptides, but they require that motional models be interpreted. Different motional models yield different interpretations of quadrupolar splittings (QS) in terms of helix orientation and dynamics. Here, we use coarse-grained (CG) molecular dynamics (MD) simulations to investigate the behavior of a well-known model transmembrane peptide, WALP23, under different hydrophobic matching/mismatching conditions. We compare experimental (2)H-NMR QS (directly measured in experiments), as well as helix tilt angle and azimuthal rotation (not directly measured), with CG MD simulation results. For QS, the agreement is significantly better than previously obtained with atomistic simulations, indicating that equilibrium sampling is more important than atomistic details for reproducing experimental QS. Calculations of helix orientation confirm that the interpretation of QS depends on the motional model used. Our simulations suggest that WALP23 can form dimers, which are more stable in an antiparallel arrangement. The origin of the preference for the antiparallel orientation lies not only in electrostatic interactions but also in better surface complementarity. In most cases, a mixture of monomers and antiparallel dimers provides better agreement with NMR data compared to the monomer and the parallel dimer. CG MD simulations allow predictions of helix orientation and dynamics and interpretation of QS data without requiring any assumption about the motional model.
| S-EPMC2931731 | biostudies-literature
Project description:Structural information of a transmembrane (TM) helix dimer is useful in understanding molecular mechanisms of important biological phenomena such as signal transduction across the cell membrane. Here, we describe an umbrella sampling (US) scheme for predicting the structure of a TM helix dimer in implicit membrane using the interhelical crossing angle and the TM-TM relative rotation angles as the reaction coordinates. This scheme conducts an efficient conformational search on TM-TM contact interfaces, and its robustness is tested by predicting the structures of glycophorin A (GpA) and receptor tyrosine kinase EphA1 (EphA1) TM dimers. The nuclear magnetic resonance (NMR) structures of both proteins correspond to the global free-energy minimum states in their free-energy landscapes. In addition, using the landscape of GpA as a reference, we also examine the protocols of temperature replica-exchange molecular dynamics (REMD) simulations for structure prediction of TM helix dimers in implicit membrane. A wide temperature range in REMD simulations, for example, 250-1000 K, is required to efficiently obtain a free-energy landscape consistent with the US simulations. The interhelical crossing angle and the TM-TM relative rotation angles can be used as reaction coordinates in multidimensional US and be good measures for conformational sampling of REMD simulations.
| S-EPMC4083740 | biostudies-literature
Project description:Methods that predict membrane helices have become increasingly useful in the context of analyzing entire proteomes, as well as in everyday sequence analysis. Here, we analyzed 27 advanced and simple methods in detail. To resolve contradictions in previous works and to reevaluate transmembrane helix prediction algorithms, we introduced an analysis that distinguished between performance on redundancy-reduced high- and low-resolution data sets, established thresholds for significant differences in performance, and implemented both per-segment and per-residue analysis of membrane helix predictions. Although some of the advanced methods performed better than others, we showed in a thorough bootstrapping experiment based on various measures of accuracy that no method performed consistently best. In contrast, most simple hydrophobicity scale-based methods were significantly less accurate than any advanced method as they overpredicted membrane helices and confused membrane helices with hydrophobic regions outside of membranes. In contrast, the advanced methods usually distinguished correctly between membrane-helical and other proteins. Nonetheless, few methods reliably distinguished between signal peptides and membrane helices. We could not verify a significant difference in performance between eukaryotic and prokaryotic proteins. Surprisingly, we found that proteins with more than five helices were predicted at a significantly lower accuracy than proteins with five or fewer. The important implication is that structurally unsolved multispanning membrane proteins, which are often important drug targets, will remain problematic for transmembrane helix prediction algorithms. Overall, by establishing a standardized methodology for transmembrane helix prediction evaluation, we have resolved differences among previous works and presented novel trends that may impact the analysis of entire proteomes.
| S-EPMC2373751 | biostudies-literature
Project description:The pharmacological effect of opioids originates, at the cellular level, by their interaction with the ?-opioid receptor (mOR) resulting in the regulation of voltage-gated Ca2+ channels and inwardly rectifying K+ channels that ultimately modulate the synaptic transmission. Recently, an alternative six trans-membrane helix isoform of mOR, (6TM-mOR) has been identified, but its function and signaling are still largely unknown. Here, we present the structural and functional mechanisms of 6TM-mOR signaling activity upon binding to morphine. Our data suggest that despite the similarity of binding modes of the alternative 6TM-mOR and the dominant seven trans-membrane helix variant (7TM-mOR), the interaction with morphine generates different dynamic responses in the two receptors, thus, promoting the activation of different mOR-specific signaling pathways. We characterize a series of 6TM-mOR-specific cellular responses, and observed that they are significantly different from those for 7TM-mOR. Morphine stimulation of 6TM-mOR does not promote a cellular cAMP response, while it increases the intracellular Ca2+ concentration and reduces the cellular K+ conductance. Our findings indicate that 6TM-mOR has a unique contribution to the cellular opioid responses. Therefore, it should be considered as a relevant target for the development of novel pharmacological tools and medical protocols involving the use of opioids.
| S-EPMC4640872 | biostudies-literature
Project description:Human Connexin26 gene mutations cause hearing loss. These hereditary mutations are the leading cause of childhood deafness worldwide. Mutations in gap junction proteins (connexins) can impair intercellular communication by eliminating protein synthesis, mis-trafficking, or inducing channels that fail to dock or have aberrant function. We previously identified a new class of mutants that form non-functional gap junction channels and hemichannels (connexons) by disrupting packing and inter-helix interactions. Here we analyzed fourteen point mutations in the fourth transmembrane helix of connexin26 (Cx26) that cause non-syndromic hearing loss. Eight mutations caused mis-trafficking (K188R, F191L, V198M, S199F, G200R, I203K, L205P, T208P). Of the remaining six that formed gap junctions in mammalian cells, M195T and A197S formed stable hemichannels after isolation with a baculovirus/Sf9 protein purification system, while C202F, I203T, L205V and N206S formed hemichannels with varying degrees of instability. The function of all six gap junction-forming mutants was further assessed through measurement of dye coupling in mammalian cells and junctional conductance in paired Xenopus oocytes. Dye coupling between cell pairs was reduced by varying degrees for all six mutants. In homotypic oocyte pairings, only A197S induced measurable conductance. In heterotypic pairings with wild-type Cx26, five of the six mutants formed functional gap junction channels, albeit with reduced efficiency. None of the mutants displayed significant alterations in sensitivity to transjunctional voltage or induced conductive hemichannels in single oocytes. Intra-hemichannel interactions between mutant and wild-type proteins were assessed in rescue experiments using baculovirus expression in Sf9 insect cells. Of the four unstable mutations (C202F, I203T, L205V, N206S) only C202F and N206S formed stable hemichannels when co-expressed with wild-type Cx26. Stable M195T hemichannels displayed an increased tendency to aggregate. Thus, mutations in TM4 cause a range of phenotypes of dysfunctional gap junction channels that are discussed within the context of the X-ray crystallographic structure.
| S-EPMC3744544 | biostudies-literature
Project description:We have previously studied the unfolding equilibrium of bacterioopsin in a single phase solvent, using Förster mechanism fluorescence resonance energy transfer (FRET) as a probe, from tryptophan donors to a dansyl acceptor. We observed an apparent unfolding transition in bacterioopsin perturbed by increasing ethanol concentrations [Nannepaga et al. (2004) Biochemistry 43, 50-59]. We have further investigated this transition and find that the unfolding is pH-dependent. We have now measured the apparent pK of acid-induced unfolding of bacterioopsin in 90% ethanol. When the acceptor is on helix B (Lys 41), the apparent pK for unfolding is 4.75; on the EF connecting loop (Cys 163), 5.15; and on helix G (Cys 222), 5.75. Five-helix proteolytic fragments are less stable. The apparent unfolding pKs are 5.46 for residues 72-248 (Cys 163) and 7.36 for residues 1-166 (Lys 41). When interpreted in terms of a simple equilibrium model for unfolding, the apparent pKs give relative free energies of unfolding in the range of -0.54 to -3.5 kcal/mol. The results suggest that the C-terminal helix of bacterioopsin is less stably folded than the N-terminal helices. We analyzed the pairwise helix-helix interaction surfaces of bacteriorhodopsin and three other seven-transmembrane-helix proteins on the basis of crystal structures. The results show that the interaction surfaces are smoother and the helix axis separations are closer in the amino-terminal two-thirds of the proteins compared with the carboxyl-terminal one-third. However, the F helix is important in stabilizing the folded structure, as shown by the instability of the 1-166 fragment. Considering the high-resolution crystal structure of bacteriorhodopsin, there are no obvious helix-helix interactions involving protein side chains which would be destabilized by protonation at the estimated pH of the unfolding transitions. However, a number of helix-bridging water molecules could become protonated, thereby weakening the helix-helix interactions.
| S-EPMC2519877 | biostudies-literature