Large Tilts in Transmembrane Helices Can Be Induced during Tertiary Structure Formation

While early structural models of helix-bundle integral membrane proteins posited that the transmembrane α-helices [transmembrane helices (TMHs)] were orientated more or less perpendicular to the membrane plane, there is now ample evidence from high-resolution structures that many TMHs have significant tilt angles relative to the membrane. Here, we address the question whether the tilt is an intrinsic property of the TMH in question or if it is imparted on the TMH during folding of the protein. Using a glycosylation mapping technique, we show that four highly tilted helices found in multi-spanning membrane proteins all have much shorter membrane-embedded segments when inserted by themselves into the membrane than seen in the high-resolution structures. This suggests that tilting can be induced by tertiary packing interactions within the protein, subsequent to the initial membrane-insertion step.     

Direct effects of ionizing radiation on integral membrane proteins. Noncovalent energy transfer requires specific interpeptide interactions

The 12 transmembrane alpha helices (TMHs) of human erythrocyte glucose transporter were individually cut by pepsin digestion as membrane-bound 2.5-3.5-kDa peptide fragments. Radiation-induced chemical degradation of these fragments showed an average target size of 34 kDa. This is 10-12 x larger than the average size of an individual TMH, demonstrating that a significant energy transfer occurs among these TMHs in the absence of covalent linkage. Heating this TMH preparation at 100 degrees C for 15 min reduced the target size to 5 kDa or less, suggesting that the noncovalent energy transfer requires specific helix-helix interactions. Purified phospholamban, a small (6-kDa) integral membrane protein containing a single TMH, formed a pentameric assembly in sodium dodecyl sulfate. The chemical degradation target size of this phospholamban pentamer was 5-6 kDa, illustrating that not all integral membrane protein assemblies permit intersubunit energy transfer. These findings together with other published observations suggest strongly that significant noncovalent energy transfer can occur within the tertiary and quaternary structure of membrane proteins and that as yet undefined proper molecular interactions are required for such covalent energy transfer. Our results with pepsin-digested glucose transporter also illustrate the importance of the interhelical interaction as a predominating force in maintaining the tertiary structure of a transmembrane protein.     

Model of the 3-D structure of the GLUT3 glucose transporter and molecular dynamics simulation of glucose transport

A molecular model of the three-dimensional (3-D) structure of the glucose transport protein, GLUT3, has been derived by homology modeling. The model was built on the basis of structural data from the MscL protein, which is a mechanosensitive ion channel, and general insights from aquaporin (a water permeation pore). Structurally conserved regions were defined by amino acid sequence comparisons, optimum interconnecting loops were selected from the protein databank, and amino (N)- and carboxy (C)-terminal ends of the protein were generated as random coil structures. The model was then subjected to energy minimization and molecular dynamics simulations in the presence of bound substrate (D-glucose). In the proposed structure of GLUT3, the 12 transmembrane (TM) helices form a right-hand barrel with a central hydrophilic pore. The pore is shaped like a funnel with dimensions of approximately 5-6 A by 8 A at its narrowest point. A network of polar and aromatic amino acids line the pore region and may facilitate the movement of glucose along the channel. A putative binding site for inhibitory ligands, such as forskolin and cytochalasin B, was identified on an intracellular aspect of the protein. Molecular dynamics studies showed that changes in the tilt and flexibility of key TM helices may modulate the opening of the pore to effect glucose transport. The proposed structure of GLUT3 may prove useful in guiding future experiments aimed at more precisely defining various functional regions of the transporter and may encourage efforts to develop models of other complex membrane proteins.     

A Simple Method for Predicting Transmembrane Proteins Based on Wavelet Transform

The increasing protein sequences from the genome project require theoretical methods to predict transmembrane helical segments (TMHs). So far, several prediction methods have been reported, but there are some deficiencies in prediction accuracy and adaptability in these methods. In this paper, a method based on discrete wavelet transform (DWT) has been developed to predict the number and location of TMHs in membrane proteins. PDB coded as 1KQG is chosen as an example to describe the prediction process by this method. 80 proteins with known 3D structure from Mptopo database are chosen at random as data sets (including 325 TMHs) and 80 sequences are divided into 13 groups according to their function and type. TMHs prediction is carried out for each group of membrane protein sequences and obtain satisfactory result. To verify the feasibility of this method, 80 membrane protein sequences are treated as test sets, 308 TMHs can be predicted and the prediction accuracy is 96.3%. Compared with the main prediction results of seven popular prediction methods, the obtained results indicate that the proposed method in this paper has higher prediction accuracy.     

MemBrain: improving the accuracy of predicting transmembrane helices

Prediction of transmembrane helices (TMH) in alpha helical membrane proteins provides valuable information about the protein topology when the high resolution structures are not available. Many predictors have been developed based on either amino acid hydrophobicity scale or pure statistical approaches. While these predictors perform reasonably well in identifying the number of TMHs in a protein, they are generally inaccurate in predicting the ends of TMHs, or TMHs of unusual length. To improve the accuracy of TMH detection, we developed a machine-learning based predictor, MemBrain, which integrates a number of modern bioinformatics approaches including sequence representation by multiple sequence alignment matrix, the optimized evidence-theoretic K-nearest neighbor prediction algorithm, fusion of multiple prediction window sizes, and classification by dynamic threshold. MemBrain demonstrates an overall improvement of about 20% in prediction accuracy, particularly, in predicting the ends of TMHs and TMHs that are shorter than 15 residues. It also has the capability to detect N-terminal signal peptides. The MemBrain predictor is a useful sequence-based analysis tool for functional and structural characterization of helical membrane proteins; it is freely available at http://chou.med.harvard.edu/bioinf/MemBrain/.     

An Analysis of Reentrant Loops

Reentrant loops are an important structural motif in alpha-helical transmembrane proteins. A reentrant loop is a structural motif that goes only halfway through the membrane and then turns back to the side from which it originates. The question of what causes the reentrant loops to form such a unique topology is still unanswered. In this study, we try to answer this question by analyzing the hydrophobicity distribution on the amino acid sequences of the reentrant loops. Our results show that reentrant loops have very low hydrophobicity around the deepest point buried in the membrane and relative high hydrophobicity close to the membrane surfaces. We speculate that this hydrophobicity distribution is a major force that stabilizes the unique reentrant loop structure. Our results also show that this hydrophobicity distribution results in special patterns on protein sequences, which can be captured using profile hidden Markov models (HMMs). The resulting profile HMMs can detect reentrant loops on protein sequences with high sensitivity and perfect specificity.     

Membrane Insertion of Marginally Hydrophobic Transmembrane Helices Depends on Sequence Context

In mammalian cells, most integral membrane proteins are initially inserted into the endoplasmic reticulum membrane by the so-called Sec61 translocon. However, recent predictions suggest that many transmembrane helices (TMHs) in multispanning membrane proteins are not sufficiently hydrophobic to be recognized as such by the translocon. In this study, we have screened 16 marginally hydrophobic TMHs from membrane proteins of known three-dimensional structure. Indeed, most of these TMHs do not insert efficiently into the endoplasmic reticulum membrane by themselves. To test if loops or TMHs immediately upstream or downstream of a marginally hydrophobic helix might influence the insertion efficiency, insertion of marginally hydrophobic helices was also studied in the presence of their neighboring loops and helices. The results show that flanking loops and nearest-neighbor TMHs are sufficient to ensure the insertion of many marginally hydrophobic helices. However, for at least two of the marginally hydrophobic helices, the local interactions are not enough, indicating that post-insertional rearrangements are involved in the folding of these proteins.     

Mapping the interactions between Escherichia coli TolQ transmembrane segments

The tolQRAB-pal operon is conserved in Gram-negative genomes. The TolQRA proteins of Escherichia coli form an inner membrane complex in which TolQR uses the proton-motive force to regulate TolA conformation and the in vivo interaction of TolA C-terminal region with the outer membrane Pal lipoprotein. The stoichiometry of the TolQ, TolR, and TolA has been estimated and suggests that 4-6 TolQ molecules are associated in the complex, thus involving interactions between the transmembrane helices (TMHs) of TolQ, TolR, and TolA. It has been proposed that an ion channel forms at the interface between two TolQ and one TolR TMHs involving the TolR-Asp(23), TolQ-Thr(145), and TolQ-Thr(178) residues. To define the organization of the three TMHs of TolQ, we constructed epitope-tagged versions of TolQ. Immunodetection of in vivo and in vitro chemically cross-linked TolQ proteins showed that TolQ exists as multimers in the complex. To understand how TolQ multimerizes, we initiated a cysteine-scanning study. Results of single and tandem cysteine substitution suggest a dynamic model of helix interactions in which the hairpin formed by the two last TMHs of TolQ change conformation, whereas the first TMH of TolQ forms intramolecular interactions.     

Conserved evolutionary units in the heme-copper oxidase superfamily revealed by novel homologous protein families

The heme-copper oxidase (HCO) superfamily includes HCOs in aerobic respiratory chains and nitric oxide reductases (NORs) in the denitrification pathway. The HCO/NOR catalytic subunit has a core structure consisting of 12 transmembrane helices (TMHs) arranged in three-fold rotational pseudosymmetry, with six conserved histidines for heme and metal binding. Using sensitive sequence similarity searches, we detected a number of novel HCO/NOR homologs and named them HCO Homology (HCOH) proteins. Several HCOH families possess only four TMHs that exhibit the most pronounced similarity to the last four TMHs (TMHs 9–12) of HCOs/NORs. Encoded by independent genes, four-TMH HCOH proteins represent a single evolutionary unit (EU) that relates to each of the three homologous EUs of HCOs/NORs comprising TMHs 1–4, TMHs 5–8, and TMHs 9–12. Single-EU HCOH proteins could form homotrimers or heterotrimers to maintain the general structure and ligand-binding sites defined by the HCO/NOR catalytic subunit fold. The remaining HCOH families, including NnrS, have 12-TMHs and three EUs. Most three-EU HCOH proteins possess two conserved histidines and could bind a single heme. Limited experimental studies and genomic context analysis suggest that many HCOH proteins could function in the denitrification pathway and in detoxification of reactive molecules such as nitric oxide. HCO/NOR catalytic subunits exhibit remarkable structural similarity to the homotrimers of MAPEG (membrane-associated proteins in eicosanoid and glutathione metabolism) proteins. Gene duplication, fusion, and fission likely play important roles in the evolution of HCOs/NORs and HCOH proteins.     

Cysteine-scanning mutagenesis of transmembrane segment 7 of the GLUT1 glucose transporter

The human erythrocyte facilitative glucose transporter (Glut1) is predicted to contain 12 transmembrane spanning alpha-helices based upon hydropathy plot analysis of the primary sequence. Five of these helices (3, 5, 7, 8, and 11) are capable of forming amphipathic structures. A model of GLUT1 tertiary structure has therefore been proposed in which the hydrophilic faces of several amphipathic helices are arranged to form a central aqueous channel through which glucose traverses the hydrophobic lipid bilayer. In order to test this model, we individually mutated each of the amino acid residues in transmembrane segment 7 to cysteine in an engineered GLUT1 molecule devoid of all native cysteines (C-less). Measurement of 2-deoxyglucose uptake in a Xenopus oocyte expression system revealed that nearly all of these mutants retain measurable transport activity. Over one-half of the cysteine mutants had significantly reduced specific activity relative to the C-less protein. The solvent accessibility and relative orientation of the residues within the helix was investigated by determining the sensitivity of the mutant transporters to inhibition by the sulfhydryl directed reagent p-chloromercuribenzene sulfonate (pCMBS). Cysteine replacement at six positions (Gln(282), Gln(283), Ile(287), Ala(289), Val(290), and Phe(291)), all near the exofacial side of the cell membrane, produced transporters that were inhibited by incubation with extracellular pCMBS. Residues predicted to be near the cytoplasmic side of the cell membrane were minimally affected by pCMBS. These data demonstrate that the exofacial portion of transmembrane segment 7 is accessible to the external solvent and provide evidence for the positioning of this alpha-helix within the glucose permeation pathway.     

Functional architecture of the CFTR chloride channel

Cystic fibrosis is caused by mutations in the cystic fibrosis transmembrane conductance regulator (CFTR), a member of the ATP-binding cassette (ABC) family of membrane transport proteins. CFTR is unique among ABC proteins in that it functions not as an active transporter but as an ATP-gated Cl^? channel. As an ion channel, the function of the CFTR transmembrane channel pore that mediates Cl^? movement has been studied in great detail. On the other hand, only low resolution structural data is available on the transmembrane parts of the protein. The structure of the channel pore has, however, been modeled on the known structure of active transporter ABC proteins. Currently, significant barriers exist to building a unified view of CFTR pore structure and function. Reconciling functional data on the channel with indirect structural data based on other proteins with very different transport functions and substrates has proven problematic. This review summarizes current structural and functional models of the CFTR Cl^? channel pore, including a comprehensive review of previous electrophysiological investigations of channel structure and function. In addition, functional data on the three-dimensional arrangement of pore-lining helices, as well as contemporary hypotheses concerning conformational changes in the pore that occur during channel opening and closing, are discussed. Important similarities and differences between different models of the pore highlight current gaps in our knowledge of CFTR structure and function. In order to fill these gaps, structural and functional models of the membrane-spanning pore need to become better integrated.     

The conformation of the cytoplasmic helix 8 of the CB1 cannabinoid receptor using NMR and circular dichroism

The cytoplasmic helix domain (fourth cytoplasmic loop, helix 8) of numerous GPCRs such as rhodopsin and the beta-adrenergic receptor exhibits unique structural and functional characteristics. Computational models also predict the existence of such a structural motif within the CB1 cannabinoid receptor, another member of the G-protein coupled receptor superfamily. To gain insights into the conformational properties of this GPCR component, a peptide corresponding to helix 8 of the CB1 receptor with a small contiguous segment from transmembrane helix 7 (TM7) was chemically synthesized and its secondary structure determined by circular dichroism (CD) and solution NMR spectroscopy. Our studies in DPC and SDS micelles revealed significant alpha-helical structure while in an aqueous medium, the peptide exhibited a random coil configuration. The relative orientation of helix 8 within the CB1 receptor was obtained from intermolecular 31P-1H and 1H-1H NOE measurements. Our results suggest that in the presence of an amphipathic membrane environment, helix 8 assumes an alpha helical structure with an orientation parallel to the phospholipid membrane surface and perpendicular to TM7. In this model, positively charged side chains interact with the lipid headgroups while the other polar side chains face the aqueous region. The above observations may be relevant to the activation/deactivation of the CB1 receptor.     

Structural basis of FFAT motif-mediated ER targeting

The FFAT motif is a targeting signal responsible for localizing a number of proteins to the cytosolic surface of the endoplasmic reticulum (ER) and to the nuclear membrane. FFAT motifs bind to members of the highly conserved VAP protein family, which are tethered to the cytoplasmic face of the ER by a C-terminal transmembrane domain. We have solved crystal structures of the rat VAP-A MSP homology domain alone and in complex with an FFAT motif. The co-crystal structure was used to design a VAP mutant that disrupts rat and yeast VAP-FFAT interactions in vitro. The FFAT binding-defective mutant also blocked function of the VAP homolog Scs2p in yeast. Finally, overexpression of the FFAT binding-defective VAP in COS7 cells dramatically altered ER morphology. Our data establish the structural basis of FFAT-mediated ER targeting and suggest that FFAT-targeted proteins play an important role in determining ER morphology.     

Structural role of glycine in amyloid fibrils formed from transmembrane alpha-helices

Amyloid fibrils associated with diseases such as Alzheimer's are often derived from the transmembrane helices of membrane proteins. It is known that the fibrils have a cross-beta-sheet structure where main chain hydrogen bonding occurs between beta-strands in the direction of the fibril axis. However, the structural basis for how the membrane-spanning helix is converted into a beta-sheet or how protofibrils associate into fibrils is not known. Here, we use a model peptide corresponding to a portion of the single transmembrane helix of glycophorin A to investigate the structural role of glycine in amyloid-like fibrils formed from transmembrane helices. Glycophorin A contains a GxxxG motif that is found in many transmembrane sequences including that of the amyloid precursor protein and prion protein. We propose that glycine, which mediates helix interactions in membrane proteins, also provides key packing motifs when it occurs in beta-sheets. We show that glycines in the glycophorin A transmembrane helix promote extended beta-strand formation when the helix partitions into aqueous environments and stabilize the packing of beta-sheets in the formation of amyloid-like fibrils. We demonstrate that fibrillization can be disrupted with a new class of inhibitors that target the molecular grooves created by glycine.     

Twelve Transmembrane Helices Form the Functional Core of Mammalian MATE1 (Multidrug and Toxin Extruder 1) Protein

The x-ray structure of the prototypic MATE family member, NorM from Vibrio cholerae, reveals a protein fold composed of 12 transmembrane helices (TMHs), confirming hydropathy analyses of the majority of (prokaryotic and plant) MATE transporters. However, the mammalian MATEs are generally predicted to have a 13(th) TMH and an extracellular C terminus. Here we affirm this prediction, showing that the C termini of epitope-tagged, full-length human, rabbit, and mouse MATE1 were accessible to antibodies from the extracellular face of the membrane. Truncation of these proteins at or near the predicted junction between the 13(th) TMH and the long cytoplasmic loop that precedes it resulted in proteins that (i) trafficked to the membrane and (ii) interacted with antibodies only after permeabilization of the plasma membrane. CHO cells expressing rbMate1 truncated at residue Gly-545 supported levels of pH-sensitive transport similar to that of cells expressing the full-length protein. Although the high transport rate of the Gly-545 truncation mutant was associated with higher levels of membrane expression (than full-length MATE1), suggesting the 13(th) TMH may influence substrate translocation, the selectivity profile of the mutant indicated that TMH13 has little impact on ligand binding. We conclude that the functional core of MATE1 consists of 12 (not 13) TMHs. Therefore, we used the x-ray structure of NorM to develop a homology model of the first 12 TMHs of MATE1. The model proved to be stable in molecular dynamic simulations and agreed with topology evident from preliminary cysteine scanning of intracellular versus extracellular loops.     

High-Resolution Modeling of Transmembrane Helical Protein Structures from Distant Homologues

Eukaryotic transmembrane helical (TMH) proteins perform a wide diversity of critical cellular functions, but remain structurally largely uncharacterized and their high-resolution structure prediction is currently hindered by the lack of close structural homologues. To address this problem, we present a novel and generic method for accurately modeling large TMH protein structures from distant homologues exhibiting distinct loop and TMH conformations. Models of the adenosine A2AR and chemokine CXCR4 receptors were first ranked in GPCR-DOCK blind prediction contests in the receptor structure accuracy category. In a benchmark of 50 TMH protein homolog pairs of diverse topology (from 5 to 12 TMHs), size (from 183 to 420 residues) and sequence identity (from 15% to 70%), the method improves most starting templates, and achieves near-atomic accuracy prediction of membrane-embedded regions. Unlike starting templates, the models are of suitable quality for computer-based protein engineering: redesigned models and redesigned X-ray structures exhibit very similar native interactions. The method should prove useful for the atom-level modeling and design of a large fraction of structurally uncharacterized TMH proteins from a wide range of structural homologues.     

Membrane-protein integration and the role of the translocation channel

Most eukaryotic membrane proteins are integrated into the lipid bilayer during their synthesis at the endoplasmic reticulum (ER). Their integration occurs with the help of a protein-conducting channel formed by the heterotrimeric Sec61 membrane-protein complex. The crystal structure of an archaeal homolog of the complex suggests mechanisms that enable the channel to open across the membrane and to release laterally hydrophobic transmembrane segments of nascent membrane proteins into lipid. Many aspects of membrane-protein integration remain controversial and poorly understood, but new structural data provide testable hypotheses. We propose a model of how the channel recognizes transmembrane segments, orients them properly with respect to the plane of the membrane and releases them into lipid. We also discuss how the channel would prevent small molecules from crossing the lipid bilayer while it is integrating proteins.     

A simple method for modeling transmembrane helix oligomers

We describe an effective procedure for modeling the structures of simple transmembrane helix homo-oligomers. The method differs from many previous approaches in that the only structural constraint we use to help select the correct model is the oligomerization state of the protein. The method involves the following steps: (1) perform 100-250 independent Monte Carlo energy minimizations of helix pairs to produce a large collection of well-packed structures; (2) filter the minimized structures to find those that are consistent with the expected symmetry of the oligomer; (3) cluster the structures that pass the symmetry filter; and (4) select a representative of the most populous cluster as the final prediction. We applied the method to the transmembrane helices of five proteins and compare our results to the available experimental data. Our predictions of glycophorin A, neu, the M2 channel and phospholamban resulted in a single model for each protein that agreed with the experimental results. In the case of erbB-2, however, we obtained three structurally distinct clusters of approximately equal sizes, so it was not possible to identify a clearly favored structure. This may reflect a real heterogeneity of packing modes for erbB-2, which is known to interact with different receptor subunits. Our method should be useful for obtaining structural models of transmembrane domains, improving our understanding of structure/function relationships for particular membrane proteins.     

SARS冠状病毒M蛋白的生物信息学研究

针对GenBank上发布的来自不同国家地区的39条SARSCoV推测M蛋白,采用生物信息学软件分析其核酸和氨基酸序列,获得其分子生物学特征,确定突变位点,预测功能结构区、Motif及抗原决定簇,比较基因突变对这些功能结构的影响.结果表明:在39个病毒株M蛋白的666 bp中,共有18个病毒株在7个位点上发生了25次变异.在M蛋白序列上预测获得3个跨膜螺旋序列和一个可能的信号肽序列.氨基酸序列的变异主要发生在其跨膜和胞外区域,胞内区域相对较少.预测发现12个Motif和7个抗原决定簇.提示突变对M蛋白的结构功能区的影响不大,也未造成M蛋白的Motif的数量和构成发生改变.对抗原决定簇的影响也主要体现在序列成分构成的改变上,在设计疫苗时,应考虑由其导致的抗原特性改变.