In vitro dimerization of the bovine papillomavirus E5 protein transmembrane domain

The E5 protein from bovine papillomavirus is a type II membrane protein and the product of the smallest known oncogene. E5 causes cell transformation by binding and activating the platelet-derived growth factor beta receptor (PDGFbetaR). In order to productively interact with the receptor, it is thought that E5 binds as a dimer. However, wild-type E5 and various mutants have also been shown to form trimers, tetramers, and even higher order oligomers. The residues in E5 that drive and stabilize a dimeric state are also still in question. At present, two different models for the E5 dimer exist in the literature, one symmetric and one asymmetric. There is universal agreement, however, that the transmembrane (TM) domain plays a vital role in stabilizing the functional oligomer; indeed, mutation of various TM domain residues can abolish E5 function. In order to better resolve the role of the E5 TM domain in function, we have undertaken the first quantitative in vitro characterization of the E5 TM domain in detergent micelles and liposomes. Circular and linear dichroism analyses verify that the TM domain adopts a stable alpha-helical structure and is able to partition efficiently across lipid bilayers. SDS-PAGE and analytical ultracentrifugation demonstrate for the first time that the TM domain of E5 forms a strong dimer with a standard state free energy of dissociation of 5.0 kcal mol (-1). We have used our new results to interpret existing models of E5 dimer formation and provide a direct link between TM helix interactions and E5 function.     

Making structural sense of dimerization interfaces of delta opioid receptor homodimers

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.     

Energetics of ErbB1 Transmembrane Domain Dimerization in Lipid Bilayers

One of the most extensively studied receptor tyrosine kinases is EGFR/ErbB1. Although our knowledge of the role of the extracellular domains and ligands in ErbB1 activation has increased dramatically based on solved domain structures, the exact mechanism of signal transduction across the membrane remains unknown. The transmembrane domains are expected to play an important role in the dimerization process, but the contribution of ErbB1 TM domain to dimer stability is not known, with published results contradicting one another. We address this controversy by showing that ErbB1 TM domain dimerizes in lipid bilayers and by calculating its contribution to stability as −2.5 kcal/mol. The stability calculations use two different methods based on Förster resonance energy transfer, which give the same result. The ErbB1 TM domain contribution to stability exceeds the change in receptor tyrosine kinases dimerization propensities that can convert normal signaling processes into pathogenic processes, and is thus likely important for biological function.     

FGFR3 Transmembrane Domain Interactions Persist in the Presence of Its Extracellular Domain

Isolated receptor tyrosine kinase transmembrane (TM) domains have been shown to form sequence-specific dimers in membranes. Yet, it is not clear whether studies of isolated TM domains yield knowledge that is relevant to full-length receptors or whether the large glycosylated extracellular domains alter the interactions between the TM helices. Here, we address this question by quantifying the effect of the pathogenic A391E TM domain mutation on the stability of the fibroblast growth factor receptor 3 dimer in the presence of the extracellular domain and comparing these results to the case of the isolated TM fibroblast growth factor receptor 3 domains. We perform the measurements in plasma membrane-derived vesicles using a Förster-resonance-energy-transfer-based method. The effect of the mutation on dimer stability in both cases is the same (∼−1.5 kcal/mol), suggesting that the interactions observed in simple TM-peptide model systems are relevant in a biological context.     

Open and Closed Conformations of the Isolated Transmembrane Domain of Death Receptor 5 Support a New Model of Activation

It has long been presumed that activation of the apoptosis-initiating Death Receptor 5, as well as other structurally homologous members of the TNF-receptor superfamily, relies on ligand-stabilized trimerization of noninteracting receptor monomers. We and others have proposed an alternate model in which the TNF-receptor dimer—sitting at the vertices of a large supramolecular receptor network of ligand-bound receptor trimers—undergoes a closed-to-open transition, propagated through a scissorslike conformational change in a tightly bundled transmembrane (TM) domain dimer. Here we have combined electron paramagnetic resonance spectroscopy and potential-of-mean force calculations on the isolated TM domain of the long isoform of DR5. The experiments and calculations both independently validate that the opening transition is intrinsic to the physical character of the TM domain dimer, with a significant energy barrier separating the open and closed states.Death receptor 5 (DR5) is a member of the tumor necrosis factor receptor (TNFR) superfamily that mediates apoptosis when bound by its cognate ligand, TNF-related apoptosis-inducing ligand (1). Upregulated in cancer cells, DR5 is among the most actively pursued anticancer targets (2). TNF-related apoptosis-inducing ligand binds to preassembled DR5 trimers at their extracellular domains, causing the formation of oligomeric ligand-receptor networks that are held together by receptor dimers (3). In the long-isoform of DR5, this dimer is crosslinked via ligand-induced disulfide bond formation between two transmembrane (TM) domain α-helices at Cys-209, and is further stabilized by a GxxxG motif one helix-turn downstream (3).Our recent study of the structurally homologous TNFR1 showed that receptor activation involves a conformational change that propagates from the extracellular domain to the cytosolic domain through a separation (or opening) of the TM domains of the dimer (4). We have therefore hypothesized that the activation of DR5, and indeed all structurally homologous TNF-receptors, involves a scissorslike opening of the TM domain dimer (Fig. 1).Open in a separate windowFigure 1Activation model of the DR5-L TM dimer. The sequence and positions of the disulfide bond and TOAC spin label (top), along with our previously published model (bottom, left) are shown. We propose an activation model (bottom, right) in which the transmembrane dimer pivots at its disulfide bond to reach an active open conformation.Using electron paramagnetic resonance (EPR) spectroscopy, a technique that has been used previously to study TM helix architecture and dynamics (5,6), and potential-of-mean force (PMF) calculations (7,8), this study addresses the question of whether the isolated disulfide-linked DR5-L TM domain dimer occupies distinct open and closed states (Fig. 1), and how its dynamic behavior contributes to the free-energy landscape of the opening transition of the full-length receptor.The DR5-L TM domain was synthesized with TOAC, an amino acid with a nitroxide spin label rigidly fixed to the α-carbon (9), incorporated at position 32 (Fig. 1), with some minor modification to facilitate EPR measurements. Previous work confirmed that this peptide forms disulfide-linked dimers (e.g., via comparison to 2-ME treated sample) and a negligible population of higher-order oligomers (further supported by model fitting of the EPR data below). For peptide work, residues were renumbered such that Thr-204 corresponds to Thr-1, and so on. The cytosolic Cys-29 (which we previously showed does not participate in a disulfide bond in cells) was replaced with serine to prevent the formation of antiparallel disulfide-linked dimers, and Trp-34 was replaced with tyrosine to prevent intrinsic fluorescence in fluorescence studies (not published). Continuous-wave (CW) dipolar EPR (sensitive only to spin-spin distances <25 Å) was used to measure TOAC-TOAC distances within the TM dimers and revealed an ordered Gaussian distribution centered at 16 Å (full width half-maximum (FWHM) = 4 Å), corresponding to a closed state (Fig. 2 A). Double electron-electron resonance (DEER) (sensitive to spin-spin distances from 15 to 60 Å) also detected a short distance consistent with the dipolar EPR data, along with a longer, disordered component (32.9 Å, FWHM = 28 Å) (Fig. 2 B). Together, these measurements indicate the presence of a compact, ordered closed state and a broader, disordered open state. EPR on oriented membranes also indicated two structural states. Global fitting revealed two populations of spin-label tilt angles (orientation of the nitroxide principal axis relative to the membrane normal): a narrow conformation (24°, FWHM = 20°), and a disordered conformation (50°, FWHM = 48°) (Fig. 2 C). This bimodal orientational distribution (Fig. 2 C) is remarkably consistent with the bimodal distance distribution (Fig. 2 B).Open in a separate windowFigure 2EPR spectra (left) of 32-TOAC-DR5 in lipid, and resulting structural distributions (right). (A) CW dipolar EPR spectra (left) of dimer (1 mM diamide) and monomer (1 mM 2-mercaptoethanol). Best-fit spin-spin distance distribution was a single Gaussian centered at 16 ± 2 Å (right). (B) The DEER waveform (left) of 32-TOAC-DR5 dimer was best fit (right) to a two-Gaussian distribution. The short distance was constrained to agree with the CW data, because DEER has poor sensitivity for distances <20 Å. The long-distance distribution is centered at 32.9 Å and is much broader. (C) CW EPR spectra (left) of 32-TOAC-DR5, with the membrane-normal oriented parallel (red) and perpendicular (blue) to the field. Simultaneous (global) fitting of these spectra reveals narrow and broad components (right). (In panels B and C, the overall distribution is plotted as black, while the closed and open components are plotted as green and magenta, respectively.)We subsequently conducted a PMF calculation (10) using the DR5-L TM dimer starting configuration developed by our group previously (3), embedded in a DMPC bilayer, with the Leu-32/Leu-32 C_α distance as the reaction coordinate. Three calculations were run from independent starting configurations, each using 50 windows spaced in 0.5° increments, and run for 20 ns at each window (totaling 3 μs). Each of the calculations yielded a similar result, and the averaged free energy curve (Fig. 3 A) agrees remarkably well with our EPR measurements: a narrow distribution at the closed conformation (∼16 Å, Fig. 3 B) separated by an ∼3 kcal/mol energy barrier from a broad distribution of accessible open conformations at ∼27 Å, (Fig. 3 C). Each of the three individual PMF plots can be found in Fig. S1 in the Supporting Material.Open in a separate windowFigure 3(A) PMF calculation of the DR5 TM domain dimer along the Leu-32/Leu-32 distance reaction coordinate. The PMF calculation reveals a narrow closed state and a broader open state separated by a free energy barrier. Representative snapshots of the (B) closed state and (C) open state.In the closed state, the helices are tightly packed at the GxxxG interfacial motif and all the way down the juxtaposed helix faces at residues Ala-18, Leu-22, Ala-25, and Val-26. The tight packing is aided by kinking and twisting of the two helices around their common axis, increasing the interacting surface area. In the open conformations, the Ala-18, Leu-22, Ala-25, and Val-26 pairs are dissociated and, interestingly, the GxxxG motif at Gly-10 and Gly-14 remains tightly packed. The open state energy well is only slightly less favorable than the closed state (by ∼2 kcal/mol), and its free energy profile is relatively broad and flat. The increased crossing angle in the open state is facilitated by straightening of the helix kink and is not accommodated by a change in bilayer thickness (see Fig. S3, A and B).The observed change in helix-helix distance (11 Å between the two minima in the PMF) is extremely close to that observed previously in live-cell FRET studies of a constitutively active form of TNFR1 (∼8 Å change between states using large fluorescence probes at the cytosolic domains) (4). The change observed in the EPR data (17 Å) may be an overestimate because the measurement is made between TOAC spin labels that likely protrude from the two helices, depending on rotational orientation. These results collectively show that activation of these receptors requires a small, but clearly significant conformational opening of the TM domains. One important note is that our EPR experiments recapitulate the equilibrium distribution of the two states despite there being no driving force to traverse the barrier between them (∼3 kcal/mol in the closed-to-open transition and ∼1 kcal/mol in the open-to-closed transition, Fig. 3). We do not interpret the results to mean that the dimer necessarily traverses these barriers at 4°C. Rather, there likely exist multiple reaction paths for dimerization of the abstracted TM domains. Finally, in the context of the full-length receptor, how the ligand induces a conformational change capable of overcoming the closed-to-open barrier remains an important question.Whether the observed structural transition in the TM domain dimer of the long-isoform of DR5 is a ubiquitous conformational switch that acts over the entire TNFR superfamily remains unknown. Vilar et al. (11) first proposed a similar scissors-model for activation of p75 neurotrophin receptor, which has a cysteine at the center of its TM helix. The short isoform of DR5 lacks a TM domain cysteine, but does form noncovalent dimers in cells, with likely TM domain dimer contacts (3). Among the other closely related and structurally homologous members of the TNFR superfamily, TNFR1 contains a cysteine at the center of the TM domain, but lacks any discernible small residue motifs (e.g., GxxxG). TNFR2 lacks a TM cysteine on the extracellular side, but does have a GxxxG motif positioned similarly to that of DR5. On the other hand, Death Receptor 4, whose functional distinction from DR5 has remained somewhat elusive, lacks both a cysteine and any recognizable small-residue hydrophobic motif.In summary, we have extended recent findings that point to the TM domain of DR5 as an essential structural component in the conformational change associated with activation. Our findings that the DR5-L TM domain occupies distinct open and closed states, separated by a substantial energy barrier, points the way to further studies across the TNF-receptor superfamily.     

The single transmembrane domains of ErbB receptors self-associate in cell membranes.

Members of the epidermal growth factor receptor, or ErbB, family of receptor tyrosine kinases have a single transmembrane (TM) alpha-helix that is usually assumed to play a passive role in ligand-induced dimerization and activation of the receptor. However, recent studies with the epidermal growth factor receptor (ErbB1) and the erythropoietin receptor have indicated that interactions between TM alpha-helices do contribute to stabilization of ligand-independent and/or ligand-induced receptor dimers. In addition, not all of the expected ErbB receptor ligand-induced dimerization events can be recapitulated using isolated extracellular domains, suggesting that other regions of the receptor, such as the TM domain, may contribute to dimerization in vivo. Using an approach for analyzing TM domain interactions in Escherichia coli cell membranes, named TOXCAT, we find that the TM domains of ErbB receptors self-associate strongly in the absence of their extracellular domains, with the rank order ErbB4-TM > ErbB1-TM equivalent to ErbB2-TM > ErbB3-TM. A limited mutational analysis suggests that dimerization of these TM domains involves one or more GXXXG motifs, which occur frequently in the TM domains of receptor tyrosine kinases and are critical for stabilizing the glycophorin A TM domain dimer. We also analyzed the effect of the valine to glutamic acid mutation in ErbB2 that constitutively activates this receptor. Contrary to our expectations, this mutation reduced rather than increased ErbB2-TM dimerization. Our findings suggest a role for TM domain interactions in ErbB receptor function, possibly in stabilizing inactive ligand-independent receptor dimers that have been observed by several groups.     

Receptor tyrosine kinase transmembrane domains: Function,dimer structure and dimerization energetics

The transmembrane (TM) domains of receptor tyrosine kinases (RTKs) play an active role in signaling. They contribute to the stability of full-length receptor dimers and to maintaining a signaling-competent dimeric receptor conformation. In an exciting new development, two structures of RTK TM domains have been solved, a break-through achievement in the field. Here we review these structures, and we discuss recent studies of RTK TM domain dimerization energetics, possible synergies between domains, and the effects of pathogenic RTK TM mutations on structure and dimerization.Key words: transmembrane domain, dimerization thermodynamics, receptor tyrosine kinases, pathogenic mutations, dimer structure     

Sodium dodecyl sulfate-polyacrylamide gel electrophoresis and forster resonance energy transfer suggest weak interactions between fibroblast growth factor receptor 3 (FGFR3) transmembrane domains in the absence of extracellular domains and ligands

Lateral dimerization of membrane proteins has evolved as a means of signal transduction across the plasma membrane for all receptor tyrosine kinases (RTKs). The transmembrane (TM) domains of RTKs are proposed to play an important role in the dimerization process. We have investigated whether the TM domains of one RTK, fibroblast growth factor receptor 3 (FGFR3), dimerize in lipid vesicles in the absence of the extracellular domains and ligands. We have performed sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) with peptides produced via solid-phase peptide synthesis that correspond to the TM domain of FGFR3. We have carried out Forster resonance energy transfer (FRET) measurements using two donor-acceptor pairs, fluorescein/rhodamine and Cy3/Cy5, as a function of peptide concentration and donor-to-acceptor mole ratios. Our results suggest that FGFR3 TM domains form sequence-specific dimers in lipid bilayers. However, the dimerization propensity of FGFR3 TM domain is much weaker than the dimerization propensity of glycophorin A (GpA), the well-characterized "membrane dimer standard". We discuss our findings in the context of cell signaling across the plasma membrane and diseases or disorders that occur due to single amino acid mutations in the TM domain of FGFR3.     

Dimerization properties of the transmembrane domains of Arabidopsis CRINKLY4 receptor-like kinase and homologs

CRINKLY4 (CR4) is a plant serine–threonine receptor kinase. In Zea mays, CR4 functions in the differentiation of the leaf epidermis and the aleurone cell layer and, in Arabidopsis thaliana, the ortholog ACR4 is involved in the development of the integument and seed coat. The Arabidopsis genome also encodes four CR4-related proteins (CRR) whose functions are not known. Based on studies of animal receptor kinase proteins it is likely that the molecular basis of function of CR4 and related proteins is mediated by receptor dimerization. The importance of the transmembrane (TM) domain in the dimerization of several receptor kinases has been demonstrated by the TOXCAT system, a genetic assay that measures helix interactions in a natural membrane environment. In this study, we have used the TOXCAT assay to investigate the potential of the CR4 and CR4-related TM domains to homo-dimerize. Our investigation indicates that the CR4 TM domain and the CRR TM domains have higher propensities for homo-dimerization than the ACR4 TM domain. Interestingly, the dimerization potential of the ACR4 TM domain is significantly weaker even though 13 of 24 amino acids are identical to that of the CR4 TM domain. In order to determine the contributions of specific amino acids to the higher dimerization potential of CR4 compared to ACR4, mutations were made at specific sites in ACR4 TM domain and the strength of the dimer assessed by the TOXCAT assay. One mutation restored the activity to the CR4 level, while other mutations produced either no change or significantly increased the dimerization potential of the ACR4 TM domain. Our results indicate that the TM domains of CR4, ACR4 and the CRR receptor family of proteins have the intrinsic capacity to homo-dimerize, albeit with varying degrees of affinity.     

Oligomerization of the fifth transmembrane domain from the adenosine A2A receptor

The human adenosine A_2A receptor (A_2AR) belongs to one of the largest family of membrane proteins, the G-protein coupled receptors (GPCRs), characterized by seven transmembrane (TM) helices. Little is known about the determinants of their structures, folding, assembly, activation mechanisms, and oligomeric states. Previous studies in our group showed that peptides corresponding to all seven TM domains form stable helical structures in detergent micelles and lipid vesicles. However, the peptides behave differently; TM5 is the only peptide to have a ratio [θ]₂₂₂/[θ]₂₀₈ obtained by circular dichroism (CD) spectroscopy>1. This finding suggested to us that TM5 might self-associate. In the present study, we investigate the unique properties of the TM5 domain. We performed detailed analyses of TM5 peptide behavior in membrane-mimetic environments using CD spectroscopy, fluorescence spectroscopy and Förster resonance energy transfer, and gel electrophoresis. We find that TM5 peptide has the ability to self-associate to form oligomeric structures in various hydrophobic milieus and that these oligomers are highly resistant to temperature and chemical denaturation. We also find that mutation of the full-length A_2AR at position M193, which is located in the fifth TM domain, noticeably alters A_2AR monomer: dimer ratio as observed on SDS-PAGE. Our results suggest that parallel association of TM5 dimers may play a role in the known adenosine A_2A receptor dimerization. This study represents the first evidence of an individual GPCR transmembrane domain self-association.     

Towards a structural understanding of the smallest known oncoprotein: investigation of the bovine papillomavirus E5 protein using solution-state NMR

The homodimeric E5 protein from bovine papillomavirus activates the platelet-derived growth factor β receptor through transmembrane (TM) helix-helix interactions leading to uncontrolled cell growth. Detailed structural information for the E5 dimer is essential if we are to uncover its unique mechanism of action. In vivo mutagenesis has been used to identify residues in the TM domain critical for dimerization, and we previously reported that a truncated synthetic E5 peptide forms dimers via TM domain interactions. Here we extend this work with the first application of high-resolution solution-state NMR to the study of the E5 TM domain in SDS micelles. Using selectively 15N-labelled peptides, we first probe sample homogeneity revealing two predominate species, which we interpret to be monomer and dimer. The equilibrium between the two states is shown to be dependent on detergent concentration, revealed by intensity shifts between two sets of peaks in 15N-(1)H HSQC experiments, highlighting the importance of sample preparation when working with these types of proteins. This information is used to estimate a free energy of association (ΔGx°=-3.05 kcal mol(-1)) for the dimerization of E5 in SDS micelles. In addition, chemical shift changes have been observed that indicate a more pronounced change in chemical environment for those residues expected to be at the dimer interface in vivo versus those that are not. Thus we are able to demonstrate our in vitro dimer is comparable to that defined in vivo, validating the biological significance of our synthetic peptide and providing a solid foundation upon which to base further structural studies. Using detergent concentration to modulate oligomeric state and map interfacial residues by NMR could prove useful in the study of other homo-oligomeric transmembrane proteins.     

Disruption of Rhodopsin Dimerization with Synthetic Peptides Targeting an Interaction Interface

Although homo- and heterodimerizations of G protein-coupled receptors (GPCRs) are well documented, GPCR monomers are able to assemble in different ways, thus causing variations in the interactive interface between receptor monomers among different GPCRs. Moreover, the functional consequences of this phenomenon, which remain to be clarified, could be specific for different GPCRs. Synthetic peptides derived from transmembrane (TM) domains can interact with a full-length GPCR, blocking dimer formation and affecting its function. Here we used peptides corresponding to TM helices of bovine rhodopsin (Rho) to investigate the Rho dimer interface and functional consequences of its disruption. Incubation of Rho with TM1, TM2, TM4, and TM5 peptides in rod outer segment (ROS) membranes shifted the resulting detergent-solubilized protein migration through a gel filtration column toward smaller molecular masses with a reduced propensity for dimer formation in a cross-linking reaction. Binding of these TM peptides to Rho was characterized by both mass spectrometry and a label-free assay from which dissociation constants were calculated. A BRET (bioluminescence resonance energy transfer) assay revealed that the physical interaction between Rho molecules expressed in membranes of living cells was blocked by the same four TM peptides identified in our in vitro experiments. Although disruption of the Rho dimer/oligomer had no effect on the rates of G protein activation, binding of G_t to the activated receptor stabilized the dimer. However, TM peptide-induced disruption of dimer/oligomer decreased receptor stability, suggesting that Rho supramolecular organization could be essential for ROS stabilization and receptor trafficking.     

FGFR3 dimer stabilization due to a single amino acid pathogenic mutation

Mutations in the transmembrane (TM) domains of receptor tyrosine kinases (RTKs) have been implicated in the induction of pathological phenotypes. These mutations are believed to stabilize the RTK dimers, and thus promote unregulated signaling. However, the energetics behind the pathology induction has not been determined. An example of a TM domain pathogenic mutation is the Ala391-->Glu mutation in fibroblast growth factor receptor 3 (FGFR3), linked to Crouzon syndrome with acanthosis nigricans, as well as to bladder cancer. Here, we determine the free energy of dimerization of wild-type and mutant FGFR3 TM domain in lipid bilayers using F?rster resonance energy transfer, and we show that hydrogen bonding between Glu391 and the adjacent helix in the dimer is a feasible mechanism for dimer stabilization. The measured change in the free energy of dimerization due to the Ala391-->Glu pathogenic mutation is -1.3 kcal/mol, consistent with previous reports of hydrogen bond strengths in proteins. This is the first quantitative measurement of mutant RTK stabilization in a membrane environment. We show that this seemingly modest value can lead to a large increase in dimer fraction and thus profoundly affect RTK-mediated signal transduction.     

Packing density of the erythropoietin receptor transmembrane domain correlates with amplification of biological responses

The formation of signal-promoting dimeric or oligomeric receptor complexes at the cell surface is modulated by self-interaction of their transmembrane (TM) domains. To address the importance of TM domain packing density for receptor functionality, we examined a set of asparagine mutants in the TM domain of the erythropoietin receptor (EpoR). We identified EpoR-T242N as a receptor variant that is present at the cell surface similar to wild-type EpoR but lacks visible localization in vesicle-like structures and is impaired in efficient activation of specific signaling cascades. Analysis by a molecular modeling approach indicated an increased interhelical distance for the EpoR-T242N TM dimer. By employing the model, we designed additional mutants with increased or decreased packing volume and confirmed a correlation between packing volume and biological responsiveness. These results propose that the packing density of the TM domain provides a novel layer for fine-tuned regulation of signal transduction and cellular decisions.     

Structural requirements of the extracellular to transmembrane domain junction for erythropoietin receptor function

The erythropoietin receptor (EpoR) is crucial for erythrocyte formation. The x-ray crystal structures of the EpoR extracellular domain lack the juxtamembrane (JM) region and the junction to the transmembrane (TM) domain. Yet the JM-TM regions are important for transmitting the conformational change imposed on the receptor dimer by Epo binding. Cysteine-scanning mutagenesis of the JM-TM regions identified three novel constitutively active mutants, demonstrating close disulfide-bonded juxtapositioning of these residues in the JM (L223C) and N-terminal TM domain (L226C, I227C). Chemical cross-linking defined the interface of the active helical TM dimer and revealed that the JM-TM segment encompassing Leu(226)-Leu(230) is non-helical. Molecular dynamics and NMR studies indicated that the TM-JM junction forms an N-terminal helix cap. This structure is important for EpoR function because replacement of this motif by consecutive leucines rendered the receptor constitutively active.     

Active and inactive orientations of the transmembrane and cytosolic domains of the erythropoietin receptor dimer

Binding of erythropoietin to the erythropoietin receptor (EpoR) extracellular domain orients the transmembrane (TM) and cytosolic regions of the receptor dimer into an unknown activated conformation. By replacing the EpoR extracellular domain with a dimeric coiled coil, we engineered TM EpoR fusion proteins where the helical TM domains were constrained into seven possible relative orientations. We identify one dimeric TM conformation that imparts full activity to the cytosolic domain of the receptor and signals via JAK2, STAT proteins, and MAP kinase, one partially active orientation that preferentially activates MAP kinase, and one conformation corresponding to the inactive receptor. The active and inactive conformations were independently identified by computational searches for low-energy TM dimeric structures. We propose a specific EpoR-activated interface and suggest its use for structural and signaling studies.     

Ligand-induced conformational changes in the Bacillus subtilis chemoreceptor McpB determined by disulfide crosslinking in vivo

Previously, we characterized the organization of the transmembrane (TM) domain of the Bacillus subtilis chemoreceptor McpB using disulfide crosslinking. Cysteine residues were engineered into serial positions along the two helices through the membrane, TM1 and TM2, as well as double mutants in TM1 and TM2, and the extent of crosslinking determined to characterize the organization of the TM domain. In this study, the organization of the TM domain was studied in the presence and absence of ligand to address what ligand-induced structural changes occur. We found that asparagine caused changes in crosslinking rate on all residues along the TM1-TM1' helical interface, whereas the crosslinking rate for almost all residues along the TM2-TM2' interface did not change. These results indicated that helix TM1 rotated counterclockwise and that TM2 did not move in respect to TM2' in the dimer on binding asparagine. Interestingly, intramolecular crosslinking of paired substitutions in 34/280 and 38/273 were unaffected by asparagine, demonstrating that attractant binding to McpB did not induce a "piston-like" vertical displacement of TM2 as seen for Trg and Tar in Escherichia coli. However, these paired substitutions produced oligomeric forms of receptor in response to ligand. This must be due to a shift of the interface between different receptor dimers, within previously suggested trimers of dimers, or even higher order complexes. Furthermore, the extent of disulfide bond formation in the presence of asparagine was unaffected by the presence of the methyl-modification enzymes, CheB and CheR, or the coupling proteins, CheW and CheV, demonstrating that these proteins must have local structural effects on the cytoplasmic domain that is not translated to the entire receptor. Finally, disulfide bond formation was also unaffected by binding proline to McpC. We conclude that ligand-binding induced a conformational change in the TM domain of McpB dimers as an excitation signal that is likely propagated within the cytoplasmic region of receptors and that subsequent adaptational events do not affect this new TM domain conformation.     

Dimerization of the erythropoietin receptor transmembrane domain in micelles

Erythropoietin receptor (EpoR) homodimerization is an initial regulatory step in erythrocyte formation. Receptor dimers form before ligand binding, suggesting that association between receptor proteins is dependent on the receptor itself. EpoR dimerization is an essential step in erythropoiesis, and misregulation of this dimerization has been implicated in several disease states, including multi-lineage leukemias; nevertheless, how EpoR regulates its own dimerization is unclear. In vivo experiments suggest the single-pass transmembrane helix is the strongest candidate for driving ligand-independent association. To address the self-association potential of this transmembrane segment, we studied its interaction energetics in micelles by utilizing a previously successful Staphylococcal nuclease (SN-EpoR TM) fusion protein. This fusion protein strategy allows expression of the EpoR transmembrane domain in Escherichia coli independent of the other EpoR domains. Sedimentation equilibrium analytical ultracentrifugation of the detergent-solubilized SN-EpoR TM demonstrated that the murine EpoR transmembrane domain self-associates to form dimers. Although this interaction is not as stable as the dimerization of the well-studied glycophorin A transmembrane dimer, the murine EpoR transmembrane domain dimer is more stable than the interactions of the colon carcinoma kinase 4 transmembrane domain. The same experiments with the human EpoR transmembrane domain, which differs from the mouse sequence by only three residues, revealed a less favorable interaction than that of the murine sequence and is only slightly more favorable than that expected for non-preferential binding. These results suggest that the mouse and human receptor proteins may differ in the roles they play in signaling.     

Molecular dynamics simulation approach for the prediction of transmembrane helix–helix heterodimers assembly

Computational methods are useful to identify favorable structures of transmembrane (TM) helix oligomers when experimental data are not available or when they cannot help to interpret helix-helix association. We report here a global search method using molecular dynamics (MD) simulations to predict the structures of transmembrane homo and heterodimers. The present approach is based only on sequence information without any experimental data and is first applied to glycophorin A to validate the protocol and to the HER2-HER3 heterodimer receptor. The method successfully reproduces the experimental structures of the TM domain of glycophorin A (GpA(TM)) with a root mean square deviation of 1.5 A. The search protocol identifies three energetically stable models of the TM domain of HER2-HER3 receptor with favorable helix-helix arrangement, including right-handed and left-handed coiled-coils. The predicted TM structures exhibit the GxxxG-like motif at the dimer interface which is presumed to drive receptor oligomerization. We demonstrate that native structures of TM domain can be predicted without quantitative experimental data. This search protocol could help to predict structures of the TM domain of HER heterodimer family.     

The achondroplasia mutation does not alter the dimerization energetics of the fibroblast growth factor receptor 3 transmembrane domain

The Gly380 --> Arg mutation in the TM domain of fibroblast growth factor receptor 3 (FGFR3) of the RTK family is linked to achondroplasia, the most common form of human dwarfism. The molecular mechanism of pathology induction is under debate, and two different mechanisms have been proposed to contribute to pathogenesis: (1) Arg380-mediated FGFR3 dimer stabilization and (2) slow downregulation of the activated mutant receptors. Here we show that the Gly380 --> Arg mutation does not alter the dimerization energetics of the FGFR3 transmembrane domain in detergent micelles or in lipid bilayers. This result indicates that pathogenesis in achondroplasia cannot be explained simply by a higher dimerization propensity of the mutant FGFR3 TM domain, thus highlighting the importance of the observed slow downregulation in phenotype induction.