Transmembrane Recognition of the Semaphorin Co-Receptors Neuropilin 1 and Plexin A1: Coarse-Grained Simulations

The cancer associated class 3 semaphorins require direct binding to neuropilins and association to plexins to trigger cell signaling. Here, we address the role of the transmembrane domains of neuropilin 1 and plexin A1 for the dimerization of the two receptors by characterizing the assembly in lipid bilayers using coarse-grained molecular dynamics simulations. From experimental evidence using a two-hybrid system showing the biochemical association of the two receptors transmembrane domains, we performed molecular simulations in DOPC and POPC demonstrating spontaneously assembly to form homodimers and heterodimers with a very high propensity for right-handed packing of the helices. Inversely, left-handed packing was observed with a very low propensity. This mode of packing was observed uniquely when the plexin A1 transmembrane domain was involved in association. Potential of mean force calculations were used to predict a hierarchy of self-association for the monomers: the two neuropilin 1 transmembrane domains strongly associated, neuropilin 1 and plexin A1 transmembrane domains associated less and the two plexin A1 transmembrane domains weakly but significantly associated. We demonstrated that homodimerization and heterodimerization are driven by GxxxG motifs, and that the sequence context modulates the packing mode of the plexin A1 transmembrane domains. This work presents major advances towards our understanding of membrane signaling platforms assembly through membrane domains and provides exquisite information for the design of antagonist drugs defining a novel class of therapeutic agents.     

Evidence for New Homotypic and Heterotypic Interactions between Transmembrane Helices of Proteins Involved in Receptor Tyrosine Kinase and Neuropilin Signaling

Signaling in eukaryotic cells frequently relies on dynamic interactions of single-pass membrane receptors involving their transmembrane (TM) domains. To search for new such interactions, we have developed a bacterial two-hybrid system to screen for both homotypic and heterotypic interactions between TM helices. We have explored the dimerization of TM domains from 16 proteins involved in both receptor tyrosine kinase and neuropilin signaling. This study has revealed several new interactions. We found that the TM domain of Mucin-4, a putative intramembrane ligand for erbB2, dimerizes not only with erbB2 but also with all four members of the erbB family. In the Neuropilin/Plexin family of receptors, we showed that the TM domains of Neuropilins 1 and 2 dimerize with themselves and also with Plexin-A1, Plexin-B1, and L1CAM, but we were unable to observe interactions with several other TM domains notably those of members of the VEGF receptor family. The potentially important Neuropilin 1/Plexin-A1 interaction was confirmed using a surface plasmon resonance assay. This work shows that TM domain interactions can be highly specific. Exploring further the propensities of TM helix–helix association in cell membrane should have important practical implications related to our understanding of the structure–function of bitopic proteins' assembly and subsequent function, especially in the regulation of signal transduction.     

Helix-helix interactions in membrane proteins: coarse-grained simulations of glycophorin a helix dimerization

Oligomerization of transmembrane (TM) helices is a key stage in the folding of membrane proteins. Glycophorin A (GpA) is a well-documented test system for this process. Coarse-grained molecular dynamics (CG-MD) allows us to simulate the self-assembly of TM helices into dimers, for both wild-type (WT) and mutant GpA sequences. For the WT sequences, dimers formed rapidly and remained stable in all simulations. The resultant dimers exhibited right-handed crossing and the same interhelix contacts as in NMR structures. Simulations of disruptive mutants revealed the dimers were less stable, with values of DeltaDelta G dimerization consistent with experimental data. The dimers of disruptive mutants were distorted relative to the WT and showed left-handed crossing of their helices. CG-MD can therefore be used to explore the interactions of TM helices, an important stage in the folding of membrane proteins. In particular, CG-MD has been shown to be sensitive enough to detect disruptions introduced by mutation. Future refinement of such models via atomistic simulations will enable a multiscale approach to predict the folding of membrane proteins.     

An automatic method for predicting transmembrane protein structures using cryo-EM and evolutionary data

The transmembrane (TM) domains of many integral membrane proteins are composed of alpha-helix bundles. Structure determination at high resolution (<4 A) of TM domains is still exceedingly difficult experimentally. Hence, some TM-protein structures have only been solved at intermediate (5-10 A) or low (>10 A) resolutions using, for example, cryo-electron microscopy (cryo-EM). These structures reveal the packing arrangement of the TM domain, but cannot be used to determine the positions of individual amino acids. The observation that typically, the lipid-exposed faces of TM proteins are evolutionarily more variable and less charged than their core provides a simple rule for orienting their constituent helices. Based on this rule, we developed score functions and automated methods for orienting TM helices, for which locations and tilt angles have been determined using, e.g., cryo-EM data. The method was parameterized with the aim of retrieving the native structure of bacteriorhodopsin among near- and far-from-native templates. It was then tested on proteins that differ from bacteriorhodopsin in their sequences, architectures, and functions, such as the acetylcholine receptor and rhodopsin. The predicted structures were within 1.5-3.5 A from the native state in all cases. We conclude that the computational method can be used in conjunction with cryo-EM data to obtain approximate model structures of TM domains of proteins for which a sufficiently heterogeneous set of homologs is available. We also show that in those proteins in which relatively short loops connect neighboring helices, the scoring functions can discriminate between near- and far-from-native conformations even without the constraints imposed on helix locations and tilt angles that are derived from cryo-EM.     

Caveolin-Na/K-ATPase interactions: Role of transmembrane topology in non-genomic steroid signal transduction

Progesterone and its polar metabolite(s) trigger the meiotic divisions in the amphibian oocyte through a non-genomic signaling system at the plasma membrane. Published site-directed mutagenesis studies of ouabain binding and progesterone-ouabain competition studies indicate that progesterone binds to a 23 amino acid extracellular loop of the plasma membrane α-subunit of Na/K-ATPase. Integral membrane proteins such as caveolins are reported to form Na/K-ATPase-peptide complexes essential for signal transduction. We have characterized the progesterone-induced Na/K-ATPase-caveolin (CAV-1)-steroid 5α-reductase interactions initiating the meiotic divisions. Peptide sequence analysis algorithms indicate that CAV-1 contains two plasma membrane spanning helices, separated by as few as 1-2 amino acid residues at the cell surface. The CAV-1 scaffolding domain, reported to interact with CAV-1 binding (CB) motifs in signaling proteins, overlaps transmembrane (TM) helix 1. The α-subunit of Na/K-ATPase (10 TM helices) contains double CB motifs within TM-1 and TM-10. Steroid 5α-reductase (6 TM helices), an initial step in polar steroid formation, contains CB motifs overlapping TM-1 and TM-6. Computer analysis predicts that interaction between antipathic strands may bring CB motifs and scaffolding domains into close proximity, initiating allostearic changes. Progesterone binding to the α-subunit may thus facilitate CB motif:CAV-1 interaction, which in turn induces helix-helix interaction and generates both a signaling cascade and formation of polar steroids.     

A computational analysis of non-genomic plasma membrane progestin binding proteins: Signaling through ion channel-linked cell surface receptors

A number of plasma membrane progestin receptors linked to non-genomic events have been identified. These include: (1) α1-subunit of the Na⁺/K⁺-ATPase (ATP1A1), (2) progestin binding PAQR proteins, (3) membrane progestin receptor alpha (mPRα), (4) progesterone receptor MAPR proteins and (5) the association of nuclear receptor (PRB) with the plasma membrane. This study compares: the pore-lining regions (ion channels), transmembrane (TM) helices, caveolin binding (CB) motifs and leucine-rich repeats (LRRs) of putative progesterone receptors. ATP1A1 contains 10 TM helices (TM-2, 4, 5, 6 and 8 are pores) and 4 CB motifs; whereas PAQR5, PAQR6, PAQR7, PAQRB8 and fish mPRα each contain 8 TM helices (TM-3 is a pore) and 2–4 CB motifs. MAPR proteins contain a single TM helix but lack pore-lining regions and CB motifs. PRB contains one or more TM helices in the steroid binding region, one of which is a pore. ATP1A1, PAQR5/7/8, mPRα, and MAPR-1 contain highly conserved leucine-rich repeats (LRR, common to plant membrane proteins) that are ligand binding sites for ouabain-like steroids associated with LRR kinases. LRR domains are within or overlap TM helices predicted to be ion channels (pore-lining regions), with the variable LRR sequence either at the C-terminus (PAQR and MAPR-1) or within an external loop (ATP1A1). Since ouabain-like steroids are produced by animal cells, our findings suggest that ATP1A1, PAQR5/7/8 and mPRα represent ion channel-linked receptors that respond physiologically to ouabain-like steroids (not progestin) similar to those known to regulate developmental and defense-related processes in plants.     

Site-directed spin labeling of a bacterial chemoreceptor reveals a dynamic, loosely packed transmembrane domain

We used site-directed spin labeling and electron paramagnetic resonance spectroscopy to investigate dynamics and helical packing in the four-helix transmembrane domain of the homodimeric bacterial chemoreceptor Trg. We focused on the first transmembrane helix, TM1, particularly on the nine-residue sequence nearest the periplasm, because patterns of disulfide formation between introduced cysteines had identified that segment as the region of closest approach among neighboring transmembrane helices. Along this sequence, mobility and accessibility of the introduced spin label were characteristic of loosely packed or solvent-exposed side chains. This was also the case for eight additional positions around the circumference and along the length of TM1. For the continuous nine-residue sequence near the periplasm, mobility and accessibility varied only modestly as a function of position. We conclude that side chains of TM1 that face the interior of the four-helix domain interact with neighboring helices but dynamic movement results in loose packing. Compared to transmembrane segments of other membrane proteins reconstituted into lipid bilayers and characterized by site-directed spin labeling, TM1 of chemoreceptor Trg is the most dynamic and loosely packed. A dynamic, loosely packed chemoreceptor domain can account for many experimental observations about the transmembrane domains of chemoreceptors.     

Modeling and docking of the three-dimensional structure of the human melanocortin 4 receptor

A three-dimensional structure of the human melanocortin 4 receptor (hMC4R) is constructed in this study using a computer-aided molecular modeling approach. Human melanocortin 4 receptor is a G Protein-Coupled Receptor (GPCR). We structurally aligned transmembrane helices with bovine rhodopsin transmembrane domains, simulated both intracellular and extracellular loop domains on homologous loop regions in other proteins of known 3D structure and modeled the C terminus on the corresponding part of bovine rhodopsin. Then tandem minimization and dynamics calculations were run to refine the crude structure. The simulative model was tested by docking with a triplet peptide (RFF) ligand. It was found that the ligand is located among transmembrane regions TM3, TM4, TM5, and TM6 of hMC4R. In consistence with mutational and biochemical data, binding site is mainly formed as a hydrophobic and negatively charged pocket. The model constructed here might provide a structural framework for making rational predictions in relevant fields.     

Polar residues in membrane domains of proteins: molecular basis for helix-helix association in a mutant CFTR transmembrane segment

Polar side chains constitute over 20% of residues in the transmembrane (TM) helices of membrane proteins, where they may serve as hydrogen bond interaction sites for phenotypic polar mutations that arise in membrane protein-related diseases. To systematically explore the structural consequences of H-bonds between TM helices, we focused on TM4 of the cystic fibrosis conductance regulator (CFTR) and its cystic fibrosis- (CF-) phenotypic mutation, V232D, as a model system. Synthetic peptides corresponding to wild-type (TM4-wt) (residues 219-242: LQASAFCGLGFLIVLALFQAGLGR) and mutant (TM4-V232D) sequences both adopt helical structures in SDS micelles and display dimer bands on SDS-PAGE arising from disulfide bond formation via wild-type residue Cys-225. However, the TM4-V232D peptide additionally forms a ladder of noncovalent oligomers, including tetramers, hexamers, and octamers, mediated by a hydrogen bond network involving Asp-Gln side chain-side chain interactions. Ala-scanning mutagenesis of the TM4 sequence indicated that ladder formation minimally required the simultaneous presence of the Cys-225, Asp-232, and Gln-237 residues. As random hydrophobic sequences containing these three residues at TM4 equivalent positions did not oligomerize, specific van der Waals packing interactions between helix side chains were also shown to play a crucial role. Overall, the results suggest that polar mutations in membrane domains, in conjunction with critically positioned polar partner residues, potentially constitute a source of aberrant helix interactions that could contribute to loss of function when they arise in protein transmembrane domains.     

Single-spanning transmembrane domains in cell growth and cell-cell interactions: More than meets the eye?

As a whole, integral membrane proteins represent about one third of sequenced genomes, and more than 50% of currently available drugs target membrane proteins, often cell surface receptors. Some membrane protein classes, with a defined number of transmembrane (TM) helices, are receiving much attention because of their great functional and pharmacological importance, such as G protein-coupled receptors possessing 7 TM segments. Although they represent roughly half of all membrane proteins, bitopic proteins (with only 1 TM helix) have so far been less well characterized. Though they include many essential families of receptors, such as adhesion molecules and receptor tyrosine kinases, many of which are excellent targets for biopharmaceuticals (peptides, antibodies, et al.). A growing body of evidence suggests a major role for interactions between TM domains of these receptors in signaling, through homo and heteromeric associations, conformational changes, assembly of signaling platforms, etc. Significantly, mutations within single domains are frequent in human disease, such as cancer or developmental disorders. This review attempts to give an overview of current knowledge about these interactions, from structural data to therapeutic perspectives, focusing on bitopic proteins involved in cell signaling.Key words: bitopic membrane proteins, transmembrane domains, transmembrane signaling, helix-helix interactions, receptors     

Assembly of transmembrane helices of simple polytopic membrane proteins from sequence conservation patterns

The transmembrane (TM) domains of most membrane proteins consist of helix bundles. The seemingly simple task of TM helix bundle assembly has turned out to be extremely difficult. This is true even for simple TM helix bundle proteins, i.e., those that have the simple form of compact TM helix bundles. Herein, we present a computational method that is capable of generating native-like structural models for simple TM helix bundle proteins having modest numbers of TM helices based on sequence conservation patterns. Thus, the only requirement for our method is the presence of more than 30 homologous sequences for an accurate extraction of sequence conservation patterns. The prediction method first computes a number of representative well-packed conformations for each pair of contacting TM helices, and then a library of tertiary folds is generated by overlaying overlapping TM helices of the representative conformations. This library is scored using sequence conservation patterns, and a subsequent clustering analysis yields five final models. Assuming that neighboring TM helices in the sequence contact each other (but not that TM helices A and G contact each other), the method produced structural models of Calpha atom root-mean-square deviation (CA RMSD) of 3-5 A from corresponding crystal structures for bacteriorhodopsin, halorhodopsin, sensory rhodopsin II, and rhodopsin. In blind predictions, this type of contact knowledge is not available. Mimicking this, predictions were made for the rotor of the V-type Na(+)-adenosine triphosphatase without such knowledge. The CA RMSD between the best model and its crystal structure is only 3.4 A, and its contact accuracy reaches 55%. Furthermore, the model correctly identifies the binding pocket for sodium ion. These results demonstrate that the method can be readily applied to ab initio structure prediction of simple TM helix bundle proteins having modest numbers of TM helices.     

Biophysical approaches to membrane protein structure determination.

Recently, there have been several technical advances in the use of solution and solid-state NMR spectroscopy to determine the structures of membrane proteins. The structures of several isolated transmembrane (TM) helices and pairs of TM helices have been solved by solution NMR methods. Similarly, the complete folds of two TM beta-barrel proteins with molecular weights of 16 and 19 kDa have been determined by solution NMR in detergent micelles. Solution NMR has also provided a first glimpse at the dynamics of an integral membrane protein. Structures of individual TM helices have also been determined by solid-state NMR. A combination of NMR with site-directed spin-label electron paramagnetic resonance or Fourier transform IR spectroscopy allows one to assemble quite detailed protein structures in the membrane.     

Protein-protein interactions in the membrane: sequence, structural, and biological motifs

Single-span transmembrane (TM) helices have structural and functional roles well beyond serving as mere anchors to tether water-soluble domains in the vicinity of the membrane. They frequently direct the assembly of protein complexes and mediate signal transduction in ways analogous to small modular domains in water-soluble proteins. This review highlights different sequence and structural motifs that direct TM assembly and discusses their roles in diverse biological processes. We believe that TM interactions are potential therapeutic targets, as evidenced by natural proteins that modulate other TM interactions and recent developments in the design of TM-targeting peptides.     

Topological analysis of DctQ, the small integral membrane protein of the C4-dicarboxylate TRAP transporter of Rhodobacter capsulatus

Tripartite ATP-independent periplasmic ('TRAP') transporters are a novel group of bacterial and archaeal secondary solute uptake systems which possess a periplasmic binding protein, but which are unrelated to ATP-binding cassette (ABC) systems. In addition to the binding protein, TRAP transporters contain two integral membrane proteins or domains, one of which is 40-50 kDa with 12 predicted transmembrane (TM) helices, thought to be the solute import protein, while the other is 20-30 kDa and of unknown function. Using a series of plasmid-encoded beta-lactamase fusions, we have determined the topology of DctQ, the smaller integral membrane protein from the high-affinity C4-dicarboxylate transporter of Rhodobacter capsulatus, which to date is the most extensively characterised TRAP transporter. DctQ was predicted by several topology prediction programmes to have four TM helices with the N- and C-termini located in the cytoplasm. The levels of ampicillin resistance conferred by the fusions when expressed in Escherichia coli were found to correlate with this predicted topology. The data have provided a topological model which can be used to test hypotheses concerning the function of the different regions of DctQ and which can be applied to other members of the DctQ family.     

Peptide mimics of the M13 coat protein transmembrane segment. Retention of helix-helix interaction motifs

Sequence-specific noncovalent helix-helix interactions between transmembrane (TM) segments in proteins are investigated by incorporating selected TM sequences into synthetic peptides using the construct CKKK-TM-KKK. The peptides are of suitable hydrophobicity for spontaneous membrane insertion, whereas formation of an N-terminal S-S bond can bring pairs of TM helices into proximity and promote their parallel orientation. Using the propensity of the protein to undergo thermally induced alpha-helix --> beta-sheet transitions as a parameter for helix stability, we compared the wild type and mutant (V29A and V31A) bacteriophage M13 coat proteins with their corresponding TM peptide constructs (M13 residues 24-42). Our results demonstrated that the relevant helix-helix tertiary contacts found in the intact proteins persist in the peptide mimics. Molecular dynamics simulations support the tight "two in-two out" dimerization motif for V31A consistent with mutagenesis data. The overall results reinforce the notion of TM segments as autonomous folding domains and suggest that the generic peptide construct provides a viable reductionist system for membrane protein structural and computational analysis.     

Dynamic motions of CD39 transmembrane domains regulate and are regulated by the enzymatic active site

The two transmembrane domains flanking the active site of CD39 regulate its activity, but little is known about the structural and dynamic features underlying their importance. Here we use a disulfide crosslinking strategy to examine transmembrane helix interactions and dynamics and to correlate these features with activity and substrate binding. We find strong intrasubunit TM1-TM2 interactions, as well as TM1-TM1' and TM2-TM2' interactions between dimer subunits, near the extracellular side of the membrane but only weak interactions near the cytoplasmic end. The specific helix faces that constitute each interface are highly flexible, indicating a significant degree of rotational mobility within the packed structure. Analysis of activity after locking the helices in various orientations via disulfide bonds suggests that not only the arrangement but also the ability of the helices to move relative to each other is crucial for enzyme function. Helix mobility is in turn modulated by substrate binding. These results suggest that rather than playing a static structural role to support an optimal active site conformation, the transmembrane domains undergo dynamic motions that underlie their functional relationship with the active site.     

Structural understanding of the transmembrane domains of inositol triphosphate receptors and ryanodine receptors towards calcium channeling.

Inositol 1,4,5-triphosphate receptors (Insp(3)Rs) and ryanodine receptors (ryRs) act as cationic channels transporting calcium ions from the endoplasmic reticulum to cytosol by forming tetramers and are proteins localized to the endoplasmic reticulum (ER). Despite the absence of classical calcium-binding motifs, calcium channeling occurs at the transmembrane domain. We have investigated putative calcium binding motifs in these sequences. Prediction methods indicate the presence of six transmembrane helices in the C-terminal domain, one of the three domains conserved between Insp(3)R and ryR receptors. The recently identified crystal structure of the K(+) channel, which also forms tetramers, revealed that two transmembrane helices, an additional pore helix and a selectivity filter are responsible for selective K(+) ion channeling. The last three TM helices of Insp(3)R and ryR are particularly well conserved and we found analogous pore helix and selectivity filter motif in these sequences. We obtained a three-dimensional structural model for the transmembrane tetramer by extrapolating the distant structural similarity to the K(+) channels.     

Identification of proline residues in the core cytoplasmic and transmembrane regions of multidrug resistance protein 1 (MRP1/ABCC1) important for transport function, substrate specificity, and nucleotide interactions

Multidrug resistance protein 1 (MRP1/ABCC1) is an ATP-binding cassette transporter that confers resistance to drugs and mediates the transport of organic anions. MRP1 has a core structure of two membrane spanning domains (MSDs) each followed by a nucleotide binding domain. This core structure is preceded by a third MSD with five transmembrane (TM) helices, whereas MSD2 and MSD3 each contain six TM helices. We investigated the consequences of Ala substitution of 18 Pro residues in both the non-membrane and TM regions of MSD2 and MSD3 on MRP1 expression and organic anion transport function. All MRP1-Pro mutants except P1113A were expressed in human embryonic kidney cells at levels comparable with wild-type MRP1. In addition, five mutants containing substitutions of Pro residues in or proximal to the TM helices of MSD2 (TM6-Pro(343), TM8-Pro(448), TM10-Pro(557), and TM11-Pro(595)) and MSD3 (TM14-Pro(1088)) exhibited significantly reduced transport of five organic anion substrates. In contrast, mutation of Pro(1150) in the cytoplasmic loop (CL7) linking TM15 to TM16 caused a substantial increase in 17beta-estradiol-17-beta-(D-glucuronide) and methotrexate transport, whereas transport of other organic anions was reduced or unchanged. Significant substrate-specific changes in the ATP dependence of transport and binding by the P1150A mutant were also observed. Our findings demonstrate the importance of TM6, TM8, TM10, TM11, and TM14 in MRP1 transport function and suggest that CL7 may play a differential role in coupling the activity of the nucleotide binding domains to the translocation of different substrates across the membrane.