Computational searching and mutagenesis suggest a structure for the pentameric transmembrane domain of phospholamban

Protein Design: From the Aspect of Water Solubility and Stability

Water solubility and structural stability are key merits for proteins defined by the primary sequence and 3D-conformation. Their manipulation represents important aspects of the protein design field that relies on the accurate placement of amino acids and molecular interactions, guided by underlying physiochemical principles. Emulated designer proteins with well-defined properties both fuel the knowledge-base for more precise computational design models and are used in various biomedical and nanotechnological applications. The continuous developments in protein science, increasing computing power, new algorithms, and characterization techniques provide sophisticated toolkits for solubility design beyond guess work. In this review, we summarize recent advances in the protein design field with respect to water solubility and structural stability. After introducing fundamental design rules, we discuss the transmembrane protein solubilization and de novo transmembrane protein design. Traditional strategies to enhance protein solubility and structural stability are introduced. The designs of stable protein complexes and high-order assemblies are covered. Computational methodologies behind these endeavors, including structure prediction programs, machine learning algorithms, and specialty software dedicated to the evaluation of protein solubility and aggregation, are discussed. The findings and opportunities for Cryo-EM are presented. This review provides an overview of significant progress and prospects in accurate protein design for solubility and stability.

Characterization of interactions within the Igα/Igβ transmembrane domains of the human B-cell receptor provides insights into receptor assembly

The B-cell receptor (BCR), a complex comprised of a membrane-associated immunoglobulin and the Igα/β heterodimer, is one of the most important immune receptors in humans and controls B-cell development, activity, selection, and death. BCR signaling plays key roles in autoimmune diseases and lymphoproliferative disorders, yet, despite the clinical significance of this protein complex, key regions (i.e., the transmembrane domains) have yet to be structurally characterized. The mechanism for BCR signaling also remains unclear and has been variously described by the mutually exclusive cross-linking and dissociation activation models. Common to these models is the significance of local plasma membrane composition, which implies that interactions between BCR transmembrane domains (TMDs) play a role in receptor functionality. Here we used an in vivo assay of TMD oligomerization called GALLEX alongside spectroscopic and computational methods to characterize the structures and interactions of human Igα and Igβ TMDs in detergent micelles and natural membranes. We observed weak self-association of the Igβ TMD and strong self-association of the Igα TMD, which scanning mutagenesis revealed was entirely stabilized by an E–X₁₀–P motif. We also demonstrated strong heterotypic interactions between the Igα and Igβ TMDs both in vitro and in vivo, which scanning mutagenesis and computational models suggest is multiconfigurational but can accommodate distinct interaction sites for self-interactions and heterotypic interactions of the Igα TMD. Taken together, these results demonstrate that the TMDs of the human BCR are sites of strong protein–protein interactions that may direct BCR assembly, endoplasmic reticulum retention, and immune signaling.

Structural basis for assembly of vertical single β-barrel viruses

The vertical double β-barrel major capsid protein (MCP) fold, fingerprint of the PRD1-adeno viral lineage, is widespread in many viruses infecting organisms across the three domains of life. The discovery of PRD1-like viruses with two MCPs challenged the known assembly principles. Here, we present the cryo-electron microscopy (cryo-EM) structures of the archaeal, halophilic, internal membrane-containing Haloarcula californiae icosahedral virus 1 (HCIV-1) and Haloarcula hispanica icosahedral virus 2 (HHIV-2) at 3.7 and 3.8 Å resolution, respectively. Our structures reveal proteins located beneath the morphologically distinct two- and three-tower capsomers and homopentameric membrane proteins at the vertices that orchestrate the positioning of pre-formed vertical single β-barrel MCP heterodimers. The cryo-EM based structures together with the proteomics data provide insights into the assembly mechanism of this type of viruses and into those with membrane-less double β-barrel MCPs.

Here, the authors present the cryo-EM structures of two archaeal, halophilic, internal membrane-containing icosahedral viruses at 3.7 and 3.8 Å resolution, providing insights into the assembly process of these and related PRD1-adeno lineage viruses.

Relative orientation of collagen molecules within a fibril: a homology model for homo sapiens type I collagen

Type I collagen is an essential extracellular protein that plays an important structural role in tissues that require high tensile strength. However, owing to the molecule's size, to date no experimental structural data are available for the Homo sapiens species. Therefore, there is a real need to develop a reliable homology model and a method to study the packing of the collagen molecules within the fibril. Through the use of the homology model and implementation of a novel simulation technique, we have ascertained the orientations of the collagen molecules within a fibril, which is currently below the resolution limit of experimental techniques. The longitudinal orientation of collagen molecules within a fibril has a significant effect on the mechanical and biological properties of the fibril, owing to the different amino acid side chains available at the interface between the molecules.

Combination of Cα-H Hydrogen Bonds and van der Waals Packing Modulates the Stability of GxxxG-Mediated Dimers in Membranes

The GxxxG motif is frequently found at the dimerization interface of a transmembrane structural motif called GAS_right, which is characterized by a short interhelical distance and a right-handed crossing angle between the helices. In GAS_right dimers, such as glycophorin A (GpA), BNIP3, and members of the ErbB family, the backbones of the helices are in contact, and they invariably display networks of 4 to 8 weak hydrogen bonds between Cα-H carbon donors and carbonyl acceptors on opposing helices (Cα-H···O═C hydrogen bonds). These networks of weak hydrogen bonds at the helix-helix interface are presumably stabilizing, but their energetic contribution to dimerization has yet to be determined experimentally. Here, we present a computational and experimental structure-based analysis of GAS_right dimers of different predicted stabilities, which show that a combination of van der Waals packing and Cα-H hydrogen bonding predicts the experimental trend of dimerization propensities. This finding provides experimental support for the hypothesis that the networks of Cα-H hydrogen bonds are major contributors to the free energy of association of GxxxG-mediated dimers. The structural comparison between groups of GAS_right dimers of different stabilities reveals distinct sequence as well as conformational preferences. Stability correlates with shorter interhelical distances, narrower crossing angles, better packing, and the formation of larger networks of Cα-H hydrogen bonds. The identification of these structural rules provides insight on how nature could modulate stability in GAS_right and finely tune dimerization to support biological function.

His499 Regulates Dimerization and Prevents Oncogenic Activation by Asparagine Mutations of the Human Thrombopoietin Receptor

Ligand binding to the extracellular domain of the thrombopoietin receptor (TpoR) imparts a specific orientation on the transmembrane (TM) and intracellular domains of the receptors that is required for physiologic activation via receptor dimerization. To map the inactive and active dimeric orientations of the TM helices, we performed asparagine (Asn)-scanning mutagenesis of the TM domains of the murine and human TpoR. Substitution of Asn at only one position (S505N) activated the human receptor, whereas Asn substitutions at several positions activated the murine receptor. Second site mutational studies indicate that His(499) near the N terminus of the TM domain is responsible for protecting the human receptor from activation by Asn mutations. Structural studies reveal that the sequence preceding His(499) is helical in the murine receptor but non-helical in peptides corresponding to the TM domain of the inactive human receptor. The activating S505N mutation and the small molecule agonist eltrombopag both induce helix in this region of the TM domain and are associated with dimerization and activation of the human receptor. Thus, His(499) regulates the activation of human TpoR and provides additional protection against activating mutations, such as oncogenic Asn mutations in the TM domain.

Prediction of the structures of helical membrane proteins based on a minimum unfavorable contacts approach

An understanding of structure-function relationships of membrane proteins continues to be a challenging problem, owing to the difficulty in obtaining their structures experimentally. This study suggests a method for modeling membrane protein structures that can be used to generate a reliable initial conformation prior to the use of other approaches for sampling conformations. It involves optimizing the orientation of hydrophilic residues so as to minimize unfavorable contacts with the hydrophobic tails of the lipid bilayer. Starting with the optimized initial conformation for three different proteins modeled based on this method, two independent approaches have been used for sampling the conformational space of the proteins. Both approaches are able to predict structures reasonably close to experimental structures, indicating that the initial structure enables the sampling of conformations that are close to the native structure. Possible improvements in the method for making it broadly applicable to helical membrane proteins are discussed.

A Gly-zipper motif mediates homodimerization of the transmembrane domain of the mitochondrial kinase ADCK3

Interactions between α-helices within the hydrophobic environment of lipid bilayers are integral to the folding and function of transmembrane proteins; however, the major forces that mediate these interactions remain debated, and our ability to predict these interactions is still largely untested. We recently demonstrated that the frequent transmembrane association motif GASright, the GxxxG-containing fold of the glycophorin A dimer, is optimal for the formation of extended networks of Cα-H hydrogen bonds, supporting the hypothesis that these bonds are major contributors to association. We also found that optimization of Cα-H hydrogen bonding and interhelical packing is sufficient to computationally predict the structure of known GASright dimers at near atomic level. Here, we demonstrate that this computational method can be used to characterize the structure of a protein not previously known to dimerize, by predicting and validating the transmembrane dimer of ADCK3, a mitochondrial kinase. ADCK3 is involved in the biosynthesis of the redox active lipid, ubiquinone, and human ADCK3 mutations cause a cerebellar ataxia associated with ubiquinone deficiency, but the biochemical functions of ADCK3 remain largely undefined. Our experimental analyses show that the transmembrane helix of ADCK3 oligomerizes, with an interface based on an extended Gly-zipper motif, as predicted by our models. The data provide strong evidence for the hypothesis that optimization of Cα-H hydrogen bonding is an important factor in the association of transmembrane helices. This work also provides a structural foundation for investigating the role of transmembrane association in regulating the biological activity of ADCK3.

General rules for the arrangements and gating motions of pore-lining helices in homomeric ion channels

The pore-lining helix (PLH) bundles are central to the function of all ion channels, as their conformational rearrangements dictate channel gating. Here we explore all plausible oligomeric arrangements of the PLH bundles within homomeric ion channels by building models using generic restraints. In particular, the distance between two neighbouring PLHs was bounded both below and above in order to avoid steric clash and allow proper packing. The resulting models provide a theoretical representation of the accessible space for oligomeric arrangements. While the represented space is confined, it encompasses nearly all the ion channel PLH bundles for which the structures are currently known. For a multitude of channels, gating models suggested by paths within the confined accessible space are in qualitative agreement with those established in previous structural and computational studies.

Rearrangements of the pore-lining helix (PLH) bundles of ion channels are central to their gating mechanisms. Here, Dai et al. use a modelling approach to define the general rules that govern the arrangements and gating motions of the PLHs in homomeric ion channels.

Mapping the homodimer interface of an optimized, artificial, transmembrane protein activator of the human erythropoietin receptor

Transmembrane proteins constitute a large fraction of cellular proteins, and specific interactions involving membrane-spanning protein segments play an important role in protein oligomerization, folding, and function. We previously isolated an artificial, dimeric, 44-amino acid transmembrane protein that activates the human erythropoietin receptor (hEPOR) in trans. This artificial protein supports limited erythroid differentiation of primary human hematopoietic progenitor cells in vitro, even though it does not resemble erythropoietin, the natural ligand of this receptor. Here, we used a directed-evolution approach to explore the structural basis for the ability of transmembrane proteins to activate the hEPOR. A library that expresses thousands of mutants of the transmembrane activator was screened for variants that were more active than the original isolate at inducing growth factor independence in mouse cells expressing the hEPOR. The most active mutant, EBC5-16, supports erythroid differentiation in human cells with activity approaching that of EPO, as assessed by cell-surface expression of glycophorin A, a late-stage marker of erythroid differentiation. EBC5-16 contains a single isoleucine to serine substitution at position 25, which increases its ability to form dimers. Genetic studies confirmed the importance of dimerization for activity and identified the residues constituting the homodimer interface of EBC5-16. The interface requires a GxxxG dimer packing motif and a small amino acid at position 25 for maximal activity, implying that tight packing of the EBC5-16 dimer is a crucial determinant of activity. These experiments identified an artificial protein that causes robust activation of its target in a natural host cell, demonstrated the importance of dimerization of this protein for engagement of the hEPOR, and provided the framework for future structure-function studies of this novel mechanism of receptor activation.

A conserved tetrameric interaction of cry toxin helix α3 suggests a functional role for toxin oligomerization

Crystal (Cry) toxins are widely used for insect control, but their mechanism of toxicity is still uncertain. These toxins can form lytic pores in vitro, and water soluble tetrameric pre-pore intermediates have been reported. Even the precise oligomeric state of the toxin in membranes, trimeric or tetrameric, is still a debated issue. Based on previous reports, we have assumed that interactions between toxin monomers in solution are at least partly mediated by domain I, and we have analyzed in silico the homo-oligomerization tendencies of the domain I α-helices individually. Using many homologous sequences for each α-helix, our strategy allows selection of evolutionarily conserved interactions. These interactions appeared only in helices α3 and α5, but only α3 produced a suitably oriented or α-helical sample in lipid bilayers, forming homotetramers in C14-betaine, and allowing determination of its rotational orientation in lipid bilayers using site-specific infrared dichroism (SSID). The determined orientation in the tetrameric model is in agreement with only one of the evolutionarily conserved models. In addition mutation R99E, which was found to inhibit oligomerization experimentally, greatly destabilized the tetramer in molecular dynamic simulations. In this model, helix 3 is able to form inter-monomer interactions without significant rearrangements of domain I, which is compatible with the available crystal structure of Cry toxins in solution. The model presented here at least partially explains the reported tetrameric oligomerization of Cry toxins in solution and the inhibition of this oligomerization by a synthetic α3 peptide.

Hepatitis C virus RNA replication and virus particle assembly require specific dimerization of the NS4A protein transmembrane domain

Hepatitis C virus (HCV) NS4A is a single-pass transmembrane (TM) protein essential for viral replication and particle assembly. The sequence of the NS4A TM domain is highly conserved, suggesting that it may be important for protein-protein interactions. To test this hypothesis, we measured the potential dimerization of the NS4A TM domain in a well-characterized two-hybrid TM protein interaction system. The NS4A TM domain exhibited a strong homotypic interaction that was comparable in affinity to glycophorin A, a well-studied human blood group antigen that forms TM homodimers. Several mutations predicted to cluster on a common surface of the NS4A TM helix caused significant reductions in dimerization, suggesting that these residues form an interface for NS4A dimerization. Mutations in the NS4A TM domain were further examined in the JFH-1 genotype 2a replicon system; importantly, all mutations that destabilized NS4A dimers also caused defects in RNA replication and/or virus assembly. Computational modeling of NS4A TM interactions suggests a right-handed dimeric interaction of helices with an interface that is consistent with the mutational effects. Furthermore, defects in NS4A oligomerization and virus particle assembly of two mutants were rescued by NS4A A15S, a TM mutation recently identified through forward genetics as a cell culture-adaptive mutation. Together, these data provide the first example of a functionally important TM dimer interface within an HCV nonstructural protein and reveal a fundamental role of the NS4A TM domain in coordinating HCV RNA replication and virus particle assembly.

Binding of MgtR, a Salmonella transmembrane regulatory peptide, to MgtC, a Mycobacterium tuberculosis virulence factor: a structural study

MgtR, a highly hydrophobic peptide expressed in Salmonella enterica serovar Typhimurium, inhibits growth in macrophages through binding to the membrane protein MgtC that has been identified as essential for replication in macrophages. While the Mycobacterium tuberculosis MgtC is highly homologous to its S. Typhi analogue, there does not appear to be an Mtb homologue for MgtR, raising significant pharmacological interest in this system. Here, solid-state NMR and EPR spectroscopy in lipid bilayer preparations were used to demonstrate the formation of a heterodimer between S. Typhi MgtR and the transmembrane helix 4 of Mtb MgtC. Based on the experimental restraints, a structural model of this heterodimer was developed using computational techniques. The result is that MgtR appears to be ideally situated in the membrane to influence the functionality of MgtC.

Inside-out signaling promotes dynamic changes in the carcinoembryonic antigen-related cellular adhesion molecule 1 (CEACAM1) oligomeric state to control its cell adhesion properties

Cell-cell contacts are fundamental to multicellular organisms and are subject to exquisite levels of control. The carcinoembryonic antigen-related cell adhesion molecule 1 (CEACAM1) can engage in both cis-homophilic (parallel) oligomerization and trans-homophilic (anti-parallel) binding. In this study, we establish that the CEACAM1 transmembrane domain has a propensity to form cis-dimers via the transmembrane-embedded (432)GXXXG(436) motif and that this basal state is overcome when activated calmodulin binds to the CEACAM1 cytoplasmic domain. Although mutation of the (432)GXXXG(436) motif reduced CEACAM1 oligomerization, it did not affect surface localization of the receptor or influence CEACAM1-dependent cellular invasion by the pathogenic Neisseria. The mutation did, however, have a striking effect on CEACAM1-dependent cellular aggregation, increasing both the kinetics of cell-cell association and the size of cellular aggregates formed. CEACAM1 association with tyrosine kinase c-Src and tyrosine phosphatases SHP-1 and SHP-2 was not affected by the (432)GXXXG(436) mutation, consistent with their association with the monomeric form of wild type CEACAM1. Collectively, our results establish that a dynamic oligomer-to-monomer shift in surface-expressed CEACAM1 facilitates trans-homophilic binding and downstream effector signaling.

Structural organization of FtsB, a transmembrane protein of the bacterial divisome

We report the first structural analysis of an integral membrane protein of the bacterial divisome. FtsB is a single-pass membrane protein with a periplasmic coiled coil. Its heterologous association with its partner FtsL represents an essential event for the recruitment of the late components to the division site. Using a combination of mutagenesis, computational modeling, and X-ray crystallography, we determined that FtsB self-associates, and we investigated its structural organization. We found that the transmembrane domain of FtsB homo-oligomerizes through an evolutionarily conserved interaction interface where a polar residue (Gln 16) plays a critical role through the formation of an interhelical hydrogen bond. The crystal structure of the periplasmic domain, solved as a fusion with Gp7, shows that 30 juxta-membrane amino acids of FtsB form a canonical coiled coil. The presence of conserved Gly residue in the linker region suggests that flexibility between the transmembrane and coiled coil domains is functionally important. We hypothesize that the transmembrane helices of FtsB form a stable dimeric core for its association with FtsL into a higher-order oligomer and that FtsL is required to stabilize the periplasmic domain of FtsB, leading to the formation of a complex that is competent for binding to FtsQ, and to their consequent recruitment to the divisome. The study provides an experimentally validated structural model and identifies point mutations that disrupt association, thereby establishing important groundwork for the functional characterization of FtsB in vivo.

The small hydrophobic protein of the human respiratory syncytial virus forms pentameric ion channels

The small hydrophobic (SH) protein is encoded by the human respiratory syncytial virus. Its absence leads to viral attenuation in the context of whole organisms, and it prevents apoptosis in infected cells. Herein, we have examined the structure of SH protein in detergent micelles and in lipid bilayers, by solution NMR and attenuated total reflection-Fourier transform infrared spectroscopy, respectively. We found that SH protein has a single α-helical transmembrane domain and forms homopentamers in several detergents. In detergent micelles, the transmembrane domain is flanked N-terminally by an α-helix that forms a ring around the lumen of the pore and C-terminally by an extended β-turn. SH protein was found in the plasma membrane of transiently expressing HEK 293 cells, which showed pH-dependent (acid-activated) channel activity. Channel activity was abolished in mutants lacking both native His residues, His(22) and His(51), but not when either His was present. Herein, we propose that the pentameric model of SH protein presented is a physiologically relevant conformation, albeit probably not the only one, in which SH contributes to RSV infection and replication. Viroporins are short (∼100 amino acids) viral membrane proteins that form oligomers of a defined size, act as proton or ion channels, and in general enhance membrane permeability in the host. However, with some exceptions, their precise biological role of their channel activity is not understood. In general, viroporins resemble poorly specialized proteins but are nevertheless critical for viral fitness. In vivo, viruses lacking viroporins usually exhibit an attenuated or weakened phenotype, altered tropism, and diminished pathological effects. We have chosen to study the SH protein, 64 amino acids long, found in the human respiratory syncytial virus because of the effect of RSV on human health and the lack of adequate antivirals. We show that SH protein forms oligomers that behave as ion channels when activated at low pH. This study adds SH protein to a growing group of viroporins that have been structurally characterized. Although the precise biological role of this pentameric channel is still unknown, this report is nevertheless essential to fill some of the many gaps that exist in the understanding of SH protein function.

Probing the helical tilt and dynamic properties of membrane-bound phospholamban in magnetically aligned bicelles using electron paramagnetic resonance spectroscopy

Wild-type phospholamban (WT-PLB), a Ca(2+)-ATPase (SERCA) regulator in the sarcoplasmic reticulum membrane, was studied using TOAC nitroxide spin labeling, magnetically aligned bicelles, and electron paramagnetic resonance (EPR) spectroscopy to ascertain structural and dynamic information. Different structural domains of PLB (transmembrane segment: positions 42 and 45, loop region: position 20, and cytoplasmic domain: position 10) were probed with rigid TOAC spin labels to extract the transmembrane helical tilt and structural dynamic information, which is crucial for understanding the regulatory function of PLB in modulating Ca(2+)-ATPase activity. Aligned experiments indicate that the transmembrane domain of wild-type PLB has a helical tilt of 13°±4° in DMPC/DHPC bicelles. TOAC spin labels placed on the WT-PLB transmembrane domain showed highly restricted motion with more than 100ns rotational correlation time (τ(c)); whereas the loop, and the cytoplasmic regions each consists of two distinct motional dynamics: one fast component in the sub-nanosecond scale and the other component is slower dynamics in the nanosecond range.

A photon-free approach to transmembrane protein structure determination

The structures of membrane proteins are generally solved using samples dissolved in micelles, bicelles, or occasionally phospholipid bilayers using X-ray diffraction or magnetic resonance. Because these are less than perfect mimics of true biological membranes, the structures are often confirmed by evaluating the effects of mutations on the properties of the protein in their native cellular environments. Low-resolution structures are also sometimes generated from the results of site-directed mutagenesis when other structural data are incomplete or not available. Here, we describe a rapid and automated approach to determine structures from data on site-directed mutants for the special case of homo-oligomeric helical bundles. The method uses as input an experimental profile of the effects of mutations on some property of the protein. This profile is then interpreted by assuming that positions that have large effects on structure/function when mutated project toward the center of the oligomeric bundle. Model bundles are generated, and correlation analysis is used to score which structures have inter-subunit C(β) distances between adjoining monomers that best correlate with the experimental profile. These structures are then clustered and refined using energy-based minimization methods. For a set of 10 homo-oligomeric TM protein structures ranging from dimers to pentamers, we show that our method predicts structures to within 1-2 Å backbone RMSD relative to X-ray and NMR structures. This level of agreement approaches the precision of NMR structures solved in different membrane mimetics.

Packing of transmembrane domain 2 of carnitine palmitoyltransferase-1A affects oligomerization and malonyl-CoA sensitivity of the mitochondrial outer membrane protein

The purpose of this study was to investigate the sequence-dependence of oligomerization of transmembrane domain 2 (TM2) of rat carnitine palmitoyltransferase 1A (rCPT1A), to elucidate the role of this domain in the function of the full-length enzyme. Oligomerization of TM2 was studied qualitatively using complementary genetic assays that facilitate measurement of helix-helix interactions in the Escherichia coli inner membrane, and multiple quantitative biophysical methods. The effects of TM2-mutations on oligomerization and malonyl-CoA inhibition of the full-length enzyme (expressed in the yeast Pichia pastoris) were quantified. Changes designed to disrupt close-packing of the GXXXG(A) motifs reduced the oligomeric state of the corresponding TM2 peptides from hexamer to trimer (or lower), a reduction also observed on mutation of the TM2 sequence in the full-length enzyme. Disruption of these GXXXG(A) motifs had a parallel effect on the malonyl-CoA sensitivity of rCPT1A, reducing the IC(50) from 30.3 ± 5.0 to 3.0 ± 0.6 μM. For all measurements, wild-type rCPT1A was used as a control alongside various appropriate (e.g., molecular mass) standards. Our results suggest that sequence-determined, TM2-mediated oligomerization is likely to be involved in the modulation of malonyl-CoA inhibition of CPT1A in response to short- and long-term changes in protein-protein and protein-lipid interactions that occur in vivo.

Thrombopoietin receptor activation: transmembrane helix dimerization, rotation, and allosteric modulation

We report how rotational variations in transmembrane (TM) helix interactions participate in the activity states of the thrombopoietin receptor (TpoR), a type 1 cytokine receptor that controls the production of blood platelets. We also explore the mechanism of small-molecule agonists that do not mimic the natural ligand. We show, by a combination of cysteine cross-linking, alanine-scanning mutagenesis, and computational simulations, that the TpoR TM dimerizes strongly and can adopt 3 different stable, rotationally related conformations, which may correspond to specific states of the full-length receptor (active, inactive, and partially active). Thus, our data suggest that signaling and inactive states of the receptor are related by receptor subunit rotations, rather than a simple monomer-dimer transition. Moreover, results from experiments with and without agonists in vitro and in cells allow us to propose a novel allosteric mechanism of action for a class of small molecules, in which they activate TpoR by binding to the TM region and by exploiting the rotational states of the dimeric receptor. Overall, our results support the emerging view of the participation of mutual rotations of the TM domains in cytokine receptor activation.

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.

Physical basis behind achondroplasia, the most common form of human dwarfism

Fibroblast growth factor receptor 3 (FGFR3) is a receptor tyrosine kinase that plays an important role in long bone development. The G380R mutation in FGFR3 transmembrane domain is known as the genetic cause for achondroplasia, the most common form of human dwarfism. Despite many studies, there is no consensus about the exact mechanism underlying the pathology. To gain further understanding into the physical basis behind the disorder, here we measure the activation of wild-type and mutant FGFR3 in mammalian cells using Western blots, and we analyze the activation within the frame of a physical-chemical model describing dimerization, ligand binding, and phosphorylation probabilities within the dimers. The data analysis presented here suggests that the mutation does not increase FGFR3 dimerization, as proposed previously. Instead, FGFR3 activity in achondroplasia is increased due to increased probability for phosphorylation of the unliganded mutant dimers. This finding has implications for the design of targeted molecular treatments for achondroplasia.

Helical membrane peptides to modulate cell function

In recent years there has been an abundance of research into the potential of helical peptides to influence cell function. These peptides have been used to achieve a variety of different outcomes from cell repair to cell death, depending upon the peptide sequence and the nature of its interactions with cell membranes and membrane proteins. In this critical review, we summarise several mechanisms by which helical peptides, acting as either transporters, inhibitors, agonists or antibiotics, can have significant effects on cell membranes and can radically affect the internal mechanisms of the cell. The various approaches to peptide design are discussed, including the role of naturally-occurring proteins in the design of these helical peptides and current breakthroughs in the use of non-natural (and therefore more stable) peptide scaffolds. Most importantly, the current successful applications of these peptides, and their potential uses in the field of medicine, are reviewed (131 references).

(15)N Solid-state NMR spectroscopic studies on phospholamban at its phosphorylated form at ser-16 in aligned phospholipid bilayers

Wild-type phospholamban (WT-PLB) is a pentameric transmembrane protein that regulates the cardiac cycle (contraction and relaxation). From a physiological prospective, unphosphorylated WT-PLB inhibits sarcoplasmic reticulum ATPase activity; whereas, its phosphorylated form relieves the inhibition in a mechanism that is not completely understood. In this study, site-specifically (15)N-Ala-11- and (15)N-Leu-7-labeled WT-PLB and the corresponding phosphorylated forms (P-PLB) were incorporated into 1,2-dioleoyl-sn-glycero-3-phosphocholine/2-dioleoyl-sn-glycero-3-phosphoethanolamine (DOPC/DOPE) mechanically oriented lipid bilayers. The aligned (15)N-labeled Ala-11 and Leu-7 WT-PLB samples show (15)N resonance peaks at approximately 71ppm and 75ppm, respectively, while the corresponding phosphorylated forms P-PLB show (15)N peaks at 92ppm and 99ppm, respectively. These (15)N chemical shift changes upon phosphorylation are significant and in agreement with previous reports, which indicate that phosphorylation of WT-PLB at Ser-16 alters the structural properties of the cytoplasmic domain with respect to the lipid bilayers.

Free-energy profiles for ions in the influenza M2-TMD channel

M(2) transmembrane domain channel (M(2)-TMD) permeation properties are studied using molecular dynamics simulations of M(2)-TMD (1NYJ) embedded in a lipid bilayer (DMPC) with 1 mol/kg NaCl or KCl saline solution. This study allows examination of spontaneous cation and anion entry into the selectivity filter. Three titration states of the M(2)-TMD tetramer are modeled for which the four His(37) residues, forming the selectivity filter, are net uncharged, +2 charged, or +3 charged. M(2)-TMD structural properties from our simulations are compared with the properties of other models extracted from NMR and X-ray studies. During 10 ns simulations, chloride ions occasionally occupy the positively-charged selectivity filter region, and from umbrella sampling simulations, Cl(-) has a lower free-energy barrier in the selectivity-filter region than either Na(+) or NH(4) (+), and NH(4) (+) has a lower free-energy barrier than Na(+). For Na(+) and Cl(-), the free-energy barriers are less than 5 kcal/mol, suggesting that the 1NYJ conformation would probably not be exquisitely proton selective. We also point out a rotameric configuration of Trp(41) that could fully occlude the channel.

Structure of the lethal phage pinhole

Perhaps the simplest of biological timing systems, bacteriophage holins accumulate during the phage morphogenesis period and then trigger to permeabilize the cytoplasmic membrane with lethal holes; thus, terminating the infection cycle. Canonical holins form very large holes that allow nonspecific release of fully-folded proteins, but a recently discovered class of holins, the pinholins, make much smaller holes, or pinholes, that serve only to depolarize the membrane. Here, we interrogate the structure of the prototype pinholin by negative-stain transmission electron-microscopy, cysteine-accessibility, and chemical cross-linking, as well as by computational approaches. Together, the results suggest that the pinholin forms symmetric heptameric structures with the hydrophilic surface of one transmembrane domain lining the surface of a central channel approximately 15 A in diameter. The structural model also suggests a rationale for the prehole state of the pinholin, the persistence of which defines the duration of the viral latent period, and for the sensitivity of the holin timing system to the energized state of the membrane.

Structure, dynamics, and ion conductance of the phospholamban pentamer

A 52-residue membrane protein, phospholamban (PLN) is an inhibitor of an adenosine-5'-triphosphate-driven calcium pump, the Ca2+-ATPase. Although the inhibition of Ca2+-ATPase involves PLN monomers, in a lipid bilayer membrane, PLN monomers form stable pentamers of unknown biological function. The recent NMR structure of a PLN pentamer depicts cytoplasmic helices extending normal to the bilayer in what is known as the bellflower conformation. The structure shows transmembrane helices forming a hydrophobic pore 4 A in diameter, which is reminiscent of earlier reports of possible ion conductance through PLN pentamers. However, recent FRET measurements suggested an alternative structure for the PLN pentamer, known as the pinwheel model, which features a narrower transmembrane pore and cytoplasmic helices that lie against the bilayer. Here, we report on structural dynamics and conductance properties of the PLN pentamers from all-atom (AA) and coarse-grained (CG) molecular dynamics simulations. Our AA simulations of the bellflower model demonstrate that in a lipid bilayer membrane or a detergent micelle, the cytoplasmic helices undergo large structural fluctuations, whereas the transmembrane pore shrinks and becomes asymmetric. Similar asymmetry of the transmembrane region was observed in the AA simulations of the pinwheel model; the cytoplasmic helices remained in contact with the bilayer. Using the CG approach, structural dynamics of both models were investigated on a microsecond timescale. The cytoplasmic helices of the CG bellflower model were observed to fall against the bilayer, whereas in the CG pinwheel model the conformation of the cytoplasmic helices remained stable. Using steered molecular dynamics simulations, we investigated the feasibility of ion conductance through the pore of the bellflower model. The resulting approximate potentials of mean force indicate that the PLN pentamer is unlikely to function as an ion channel.

Model generation of viral channel forming 2B protein bundles from polio and coxsackie viruses

2B is a 99 amino acid membrane protein encoded by enteroviruses such as polio and coxsackie viruses with two transmembrane domains. The protein is found to make membranes of infected cells permeable. Using a computational approach which positions the models and assesses stability by molecular dynamics (MD) simulations a putative tetrameric bundle model of 2B is generated. The bundles show a pore lining motif of three lysines followed by a serine. The bundle is discussed in terms of different possible orientations of the helices in the membrane and the consequences this has on the in vivo activity of 2B.

The strong dimerization of the transmembrane domain of the fibroblast growth factor receptor (FGFR) is modulated by C-terminal juxtamembrane residues

The fibroblast growth factor receptor 3 (FGFR3) is a member of the FGFR subfamily of the receptor tyrosine kinases (RTKs) involved in signaling across the plasma membrane. Generally, ligand binding leads to receptor dimerization and activation. Dimerization involves the transmembrane (TM) domain, where mutations can lead to constitutive activation in certain cancer types and also in skeletal malformations. Thus, it has been postulated that FGFR homodimerization must be inherently weak to allow regulation, a feature reminiscent of alpha and beta integrin TM interactions. However, we show herein that in FGFR3-TM, four C-terminal residues, CRLR, have a profound destabilizing effect in an otherwise strongly dimerizing TM peptide. In the absence of these four residues, the dimerizing propensity of FGFR3-TM is comparable to glycophorin, as shown using various detergents. In addition, the expected enhanced dimerization induced by the mutation associated to the Crouzon syndrome A391E, was observed only when these four C-terminal residues were present. In the absence of these four residues, A391E was dimer-destabilizing. Finally, using site specific infrared dichroism and convergence with evolutionary conservation data, we have determined the backbone model of the FGFR3-TM homodimer in model lipid bilayers. This model is consistent with, and correlates with the effects of, most known pathological mutations found in FGFR-TM.

Artificial transmembrane oncoproteins smaller than the bovine papillomavirus E5 protein redefine sequence requirements for activation of the platelet-derived growth factor beta receptor

The bovine papillomavirus E5 protein (BPV E5) is a 44-amino-acid homodimeric transmembrane protein that binds directly to the transmembrane domain of the platelet-derived growth factor (PDGF) beta receptor and induces ligand-independent receptor activation. Three specific features of BPV E5 are considered important for its ability to activate the PDGF beta receptor and transform mouse fibroblasts: a pair of C-terminal cysteines, a transmembrane glutamine, and a juxtamembrane aspartic acid. By using a new genetic technique to screen libraries expressing artificial transmembrane proteins for activators of the PDGF beta receptor, we isolated much smaller proteins, from 32 to 36 residues, that lack all three of these features yet still dimerize noncovalently, specifically activate the PDGF beta receptor via its transmembrane domain, and transform cells efficiently. The primary amino acid sequence of BPV E5 is virtually unrecognizable in some of these proteins, which share as few as seven consecutive amino acids with the viral protein. Thus, small artificial proteins that bear little resemblance to a viral oncoprotein can nevertheless productively interact with the same cellular target. We speculate that similar cellular proteins may exist but have been overlooked due to their small size and hydrophobicity.

Interaction and conformational dynamics of membrane-spanning protein helices

Within 1 or 2 decades, the reputation of membrane-spanning alpha-helices has changed dramatically. Once mostly regarded as dull membrane anchors, transmembrane domains are now recognized as major instigators of protein-protein interaction. These interactions may be of exquisite specificity in mediating assembly of stable membrane protein complexes from cognate subunits. Further, they can be reversible and regulatable by external factors to allow for dynamic changes of protein conformation in biological function. Finally, these helices are increasingly regarded as dynamic domains. These domains can move relative to each other in different functional protein conformations. In addition, small-scale backbone fluctuations may affect their function and their impact on surrounding lipid shells. Elucidating the ways by which these intricate structural features are encoded by the amino acid sequences will be a fascinating subject of research for years to come.

The assembly motif of a bacterial small multidrug resistance protein

Multidrug transporters such as the small multidrug resistance (SMR) family of bacterial integral membrane proteins are capable of conferring clinically significant resistance to a variety of common therapeutics. As antiporter proteins of approximately 100 amino acids, SMRs must self-assemble into homo-oligomeric structures for efflux of drug molecules. Oligomerization centered at transmembrane helix four (TM4) has been implicated in SMR assembly, but the full complement of residues required to mediate its self-interaction remains to be characterized. Here, we use Hsmr, the 110-residue SMR family member of the archaebacterium Halobacterium salinarum, to determine the TM4 residue motif required to mediate drug resistance and SMR self-association. Twelve single point mutants that scan the central portion of the TM4 helix (residues 85-104) were constructed and were tested for their ability to confer resistance to the cytotoxic compound ethidium bromide. Six residues were found to be individually essential for drug resistance activity (Gly(90), Leu(91), Leu(93), Ile(94), Gly(97), and Val(98)), defining a minimum activity motif of (90)GLXLIXXGV(98) within TM4. When the propensity of these mutants to dimerize on SDS-PAGE was examined, replacements of all but Ile resulted in approximately 2-fold reduction of dimerization versus the wild-type antiporter. Our work defines a minimum activity motif of (90)GLXLIXXGV(98) within TM4 and suggests that this sequence mediates TM4-based SMR dimerization along a single helix surface, stabilized by a small residue heptad repeat sequence. These TM4-TM4 interactions likely constitute the highest affinity locus for disruption of SMR function by directly targeting its self-assembly mechanism.

Self-association of transmembrane domain 2 (TM2), but not TM1, in carnitine palmitoyltransferase 1A: role of GXXXG(A) motifs

Carnitine palmitoyltransferase 1 (CPT1) controls the rate of entry of long-chain fatty acids into the mitochondrial matrix for beta-oxidation and has been reported to exist as an oligomer. We have investigated the in vivo oligomerization of full-length rat CPT1A (rCPT1A) along with those of the N-terminal truncation/deletion mutants Delta(1-82), Delta(1-18), and Delta(19-30) expressed in yeast mitochondria. The data indicate that in liver mitochondria in vivo CPT1A exists as a hexamer but that during preparation and storage of mitochondria the order of oligomerization is rapidly reduced to the trimer, such that a mixture of hexamer and trimer is observed in isolated mitochondria in vitro. Mutants bearing deletions of different segments of the N terminus (including the more N-terminal of the two transmembrane domains) have the same pattern of oligomerization when expressed in yeast mitochondria. The self-association of the individual rCPT1A transmembrane (TM) domains (TM1, TM2) was also studied using the TOXCAT assay (which measures TM self-association in the Escherichia coli inner membrane). There was minimal self-association of the sequence corresponding to TM1 but significant self-association of TM2 in TOXCAT. Chemical cross-linking and analytical ultracentrifugation of a TM2-derived synthetic peptide showed oligomerization with a similar trimer/hexamer equilibrium to that observed for native rCPT1A in isolated mitochondria. Therefore, there was a correlation between the oligomerization behavior of TM2 peptide and that of the full-length protein. In silico molecular modeling of rCPT1A TM2 highlighted the favorable orientation of GXXXG and GXXXA motifs in the formation of the TM2 hexamer.

Structural flexibility of the pentameric SARS coronavirus envelope protein ion channel

Coronaviruses contain a small envelope membrane protein with cation-selective ion channel activity mediated by its transmembrane domain (ETM). In a computational study, we proposed that ion channel activity can be explained by either of two similar ETM homopentameric transmembrane α-helical bundles, related by a ∼50° rotation of the helices. Later, we tested this prediction, using site-specific infrared dichroism of a lysine-flanked isotopically labeled ETM peptide from the virus responsible for the severe acute respiratory syndrome, SARS, reconstituted in lipid bilayers. However, the data were consistent with the presence of a kink at the center of the ETM α-helix, and it did not fit completely either computational model. Herein, we have used native ETM, without flanking lysines, and show that the helix orientation is now consistent with one of the predicted models. ETM only produced one oligomeric form, pentamers, in the lipid-mimic detergent dodecylphosphocholine and in perfluorooctanoic acid. We thus report the correct backbone model for the pentameric α-helical bundle of ETM. The disruptive effects caused by terminal lysines probably highlight the conformational flexibility required during ion channel function.

A new approach to the prediction of transmembrane structures

About 20%–30% of genome products have been predicted as membrane proteins, which have significant biological functions. The prediction of the amount and position for the transmembrane protein helical segments (TMHs) is the hot spot in bioinformatics. In this paper, a new approach, maximum spectrum of continuous wavelet transform (MSCWT), is proposed to predict TMHs. The predictions for eight SARS-CoV membrane proteins indicate that MSCWT has the same capacity with software TMpred. Moreover, the test on a dataset of 131 structure-known proteins with 548 TMHs shows that the prediction accuracy of MSCWT for TMHs is 91.6% and that for membrane protein is 89.3%.

Transmembrane helices that form two opposite homodimeric interactions: an asparagine scan study of alphaM and beta2 integrins

Integrins are alpha/beta heterodimers, but recent in vitro and in vivo experiments also suggest an ability to associate through their transmembrane domains to form homomeric interactions. While the results of some in vitro experiments are consistent with an interaction mediated by a GxxxG-like motif, homo-oligomers observed after in vivo cross-linking are consistent with an almost opposite helix-helix interface. We have shown recently that both models of interaction are compatible with evolutionary conservation data, and we predicted that the alpha-helices in both models would have a similar rotational orientation. Herein, we have tested our prediction using in vitro asparagine scan of five consecutive residues along the GxxxG-like motif of the transmembrane domain of alpha and beta integrins, alphaM and beta2. We show that Asn-mediated dimerization occurs twice for every turn of the helix, consistent with two almost opposite forms of interaction as suggested previously for alphaIIb and beta3 transmembrane domains. The orientational parameters helix tilt and rotational orientation of each of these two Asn-stabilized dimers were measured by site-specific infrared dichroism (SSID) in model lipid bilayers and were found to be consistent with our predicted computational models. Our results highlight an intrinsic tendency for integrin transmembrane alpha-helices to form two opposite types of homomeric interaction in addition to their heteromeric interactions and suggest that integrins may form complex and specific networks at the transmembrane domain during function.

Structure and ion channel activity of the human respiratory syncytial virus (hRSV) small hydrophobic protein transmembrane domain

The small hydrophobic (SH) protein from the human respiratory syncytial virus (hRSV) is a glycoprotein of approximately 64 amino acids with one putative alpha-helical transmembrane domain. Although SH protein is important for viral infectivity, its exact role during viral infection is not clear. Herein, we have studied the secondary structure, orientation, and oligomerization of the transmembrane domain of SH (SH-TM) in the presence of lipid bilayers. Only one oligomer, a pentamer, was observed in PFO-PAGE. Using polarized attenuated total reflection-Fourier transform infrared (PATR-FTIR) spectroscopy, we show that the SH-TM is alpha-helical. The rotational orientation of SH-TM was determined by site-specific infrared dichroism (SSID) at two consecutive isotopically labeled residues. This orientation is consistent with that of an evolutionary conserved pentameric model obtained from a global search protocol using 13 homologous sequences of RSV. Conductance studies of SH-TM indicate ion channel activity, which is cation selective, and inactive below the predicted pK(a) of histidine. Thus, our results provide experimental evidence that the transmembrane domain of SH protein forms pentameric alpha-helical bundles that form cation-selective ion channels in planar lipid bilayers. We provide a model for this pore, which should be useful in mutagenesis studies to elucidate its role during the virus cycle.

Transmembrane protein models based on high-throughput molecular dynamics simulations with experimental constraints

Elucidating the structure of transmembrane proteins domains with high-resolution methods is a difficult and sometimes impossible task. Here, we explain the method of combining a limited amount of experimental data with automated high-throughput molecular dynamics (MD) simulations of alpha-helical transmembrane bundles in an explicit lipid bilayer/water environment. The procedure uses a systematic conformational search of the helix rotation with experimentally constrained MDs simulations. The experimentally determined helix tilt and rotational angle of a labeled residue with site-specific infrared dichroism allows us to select a unique high-resolution model from a number of possible energy minima encountered in the systematic conformational search.

A trimeric building block model for Cry toxins in vitro ion channel formation

The crystal (Cry) insecticidal toxins, or delta-endotoxins, are lethal to a wide variety of insect larvae, and are therefore very important in insect control. Toxicity has been explained by formation of transmembrane oligomeric pores or ion channels and, more recently, by the ability of the monomeric toxin to subvert cellular signaling pathways. The structure, topology, and precise role of the putative pore in toxicity are not known. However, in vitro biophysical studies suggest that helices alpha4 and alpha5 in domain I insert into the lipid bilayer as an alpha-helical hairpin. Mutagenesis studies have assigned an important role to alpha5 in maintaining oligomerization, and to alpha4 in channel formation. To detect the possible homo-oligomerizing tendencies of these two helices, we have used the evolutionary conservation data contained in sixteen Cry homologs in order to filter non-native interactions found during a global conformational search. No conserved homo-oligomer was found for alpha4, but a right handed trimeric alpha5 model was present in the simulations of all Cry sequences. We propose a model for Cry toxin oligomerization based on sequence analysis and available mutagenesis data.

Side chain and backbone dynamics of phospholamban in phospholipid bilayers utilizing 2H and 15N solid-state NMR spectroscopy

2H and 15N solid-state NMR spectroscopic techniques were used to investigate both the side chain and backbone dynamics of wild-type phospholamban (WT-PLB) and its phosphorylated form (P-PLB) incorporated into 1-palmitoyl-2-oleoyl-sn-glycerophosphocholine (POPC) phospholipid bilayers. 2H NMR spectra of site-specific CD3-labeled WT-PLB (at Leu51, Ala24, and Ala15) in POPC bilayers were similar under frozen conditions (-25 degrees C). However, significant differences in the line shapes of the 2H NMR spectra were observed in the liquid crystalline phase at and above 0 degrees C. The 2H NMR spectra indicate that Leu51, located toward the lower end of the transmembrane (TM) helix, shows restricted side chain motion, implying that it is embedded inside the POPC lipid bilayer. Additionally, the line shape of the 2H NMR spectrum of CD3-Ala24 reveals more side chain dynamics, indicating that this residue (located in the upper end of the TM helix) has additional backbone and internal side chain motions. 2H NMR spectra of both WT-PLB and P-PLB with CD3-Ala15 exhibit strong isotropic spectral line shapes. The dynamic isotropic nature of the 2H peak can be attributed to side chain and backbone motions to residues located in an aqueous environment outside the membrane. Also, the spectra of 15N-labeled amide WT-PLB at Leu51 and Leu42 residues showed only a single powder pattern component indicating that these two 15N-labeled residues located in the TM helix are motionally restricted at 25 degrees C. Conversely, 15N-labeled amide WT-PLB at Ala11 located in the cytoplasmic domain showed both powder and isotropic components at 25 degrees C. Upon phosphorylation, the mobile component contribution increases at Ala11. The 2H and 15N NMR data indicate significant backbone motion for the cytoplasmic domain of WT-PLB when compared to the transmembrane section.

The structural topology of wild-type phospholamban in oriented lipid bilayers using 15N solid-state NMR spectroscopy

For the first time, 15N solid-state NMR experiments were conducted on wild-type phospholamban (WT-PLB) embedded inside mechanically oriented phospholipid bilayers to investigate the topology of its cytoplasmic and transmembrane domains. 15N solid-state NMR spectra of site-specific 15N-labeled WT-PLB indicate that the transmembrane domain has a tilt angle of 13 degrees+/-6 degrees with respect to the POPC (1-palmitoyl-2-oleoyl-sn-glycero-phosphocholine) bilayer normal and that the cytoplasmic domain of WT-PLB lies on the surface of the phospholipid bilayers. Comparable results were obtained from site-specific 15N-labeled WT-PLB embedded inside DOPC/DOPE (1,2-dioleoyl-sn-glycero-3-phosphocholine/1,2-dioleoyl-sn-glycero-3-phosphoethanolamine) mechanically oriented phospholipids' bilayers. The new NMR data support a pinwheel geometry of WT-PLB, but disagree with a bellflower structure in micelles, and indicate that the orientation of the cytoplasmic domain of the WT-PLB is similar to that reported for the monomeric AFA-PLB mutant.

Spectroscopic validation of the pentameric structure of phospholamban

Phospholamban (PLN) regulates calcium translocation within cardiac myocytes by shifting sarco(endo)plasmic reticulum Ca(2+)-ATPase (SERCA) affinity for calcium. Although the monomeric form of PLN (6 kDa) is the principal inhibitory species, recent evidence suggests that the PLN pentamer (30 kDa) also is able to bind SERCA. To date, several membrane architectures of the pentamer have been proposed, with different topological orientations for the cytoplasmic domain: (i) extended from the bilayer normal by 50-60 degrees; (ii) continuous alpha-helix tilted 28 degrees relative to the bilayer normal; (iii) pinwheel geometry, with the cytoplasmic helix perpendicular to the bilayer normal and in contact with the surface of the bilayer; and (iv) bellflower structure, in which the cytoplasmic domain helix makes approximately 20 degrees angle with respect to the membrane bilayer normal. Using a variety of cell membrane mimicking systems (i.e., lipid vesicles, oriented lipid bilayers, and detergent micelles) and a combination of multidimensional solution/solid-state NMR and EPR spectroscopies, we tested the different structural models. We conclude that the pinwheel topology is the predominant conformation of pentameric PLN, with the cytoplasmic domain interacting with the membrane surface. We propose that the interaction with the bilayer precedes SERCA binding and may mediate the interactions with other proteins such as protein kinase A and protein phosphatase 1.

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.