Energetics, stability, and prediction of transmembrane helices

We show that the peptide backbone of an alpha-helix places a severe thermodynamic constraint on transmembrane (TM) stability. Neglect of this constraint by commonly used hydrophobicity scales underlies the notorious uncertainty of TM helix prediction by sliding-window hydropathy plots of membrane protein (MP) amino acid sequences. We find that an experiment-based whole-residue hydropathy scale (WW scale), which includes the backbone constraint, identifies TM helices of membrane proteins with an accuracy greater than 99 %. Furthermore, it correctly predicts the minimum hydrophobicity required for stable single-helix TM insertion observed in Escherichia coli. In order to improve membrane protein topology prediction further, we introduce the augmented WW (aWW) scale, which accounts for the energetics of salt-bridge formation. An important issue for genomic analysis is the ability of the hydropathy plot method to distinguish membrane from soluble proteins. We find that the method falsely predicts 17 to 43 % of a set of soluble proteins to be MPs, depending upon the hydropathy scale used.     

UMARS: Un-MAppable Reads Solution

Background

Few high-resolution structures of integral membranes proteins are available, as crystallization of such proteins needs yet to overcome too many technical limitations. Nevertheless, prediction of their transmembrane (TM) structure by bioinformatics tools provides interesting insights on the topology of these proteins.

Methods

We describe here how to extract new information from the analysis of hydrophobicity variations or hydrophobic pulses (HPulses) in the sequence of integral membrane proteins using the Hydrophobic Pulse Predictor, a new tool we developed for this purpose. To analyze the primary sequence of 70 integral membrane proteins we defined two levels of analysis: G1-HPulses for sliding windows of n = 2 to 6 and G2-HPulses for sliding windows of n = 12 to 16.

Results

The G2-HPulse analysis of 541 transmembrane helices allowed the definition of the new concept of transmembrane unit (TMU) that groups together transmembrane helices and segments with potential adjacent structures. In addition, the G1-HPulse analysis identified helix irregularities that corresponded to kinks, partial helices or unannotated structural events. These irregularities could represent key dynamic elements that are alternatively activated depending on the channel status as illustrated by the crystal structures of the lactose permease in different conformations.

Conclusions

Our results open a new way in the understanding of transmembrane secondary structures: hydrophobicity through hydrophobic pulses strongly impacts on such embedded structures and is not confined to define the transmembrane status of amino acids.     

Prediction of the Tertiary Structure of the Complement Control Protein Module

Complement control protein (CCP) modules, or short consensus repeats (SCR), exist in a wide variety of complement and adhesion proteins, principally the selectins. We have predicted the three-dimensional structure of a CCP module based upon secondary structural information derived by two-dimensional NMR [Barlow et al. (1991), Biochemistry 30, 997–1004]. Accordingly, the CCP is predicted to contain seven -strands with extensive hydrogen-bonding interactions, and shows a compact, globular structure. Comparison of this model to the X-ray structure of a kringle domain suggests that the CCP unit is more compact than a kringle structure, and that despite their similarities in size and disulfide bond format, the two are not homologous. Although the function of CCP domains is unknown, it is hoped that the structural model presented herein will facilitate further inquiry into how they contribute to so many systems of biological importance.     

Orientational preferences of neighboring helices can drive ER insertion of a marginally hydrophobic transmembrane helix

α-helical integral membrane proteins critically depend on the correct insertion of their transmembrane α helices into the lipid bilayer for proper folding, yet a surprisingly large fraction of the transmembrane α helices in multispanning integral membrane proteins are not sufficiently hydrophobic to insert into the target membrane by themselves. How can such marginally hydrophobic segments nevertheless form transmembrane helices in the folded structure? Here, we show that a transmembrane helix with a strong orientational preference (N(cyt)-C(lum) or N(lum)-C(cyt)) can both increase and decrease the hydrophobicity threshold for membrane insertion of a neighboring, marginally hydrophobic helix. This effect helps explain the "missing hydrophobicity" in polytopic membrane proteins.     

On the nature of the uncoupling effect of fatty acids

The effect of palmitic acid on the electrical potential difference across the inner mitochondrial membrane appears to depend on the medium in which mitochondria are incubated. In medium A (cf. Luvisettoet al. (1987),Biochemistry,26, 7332–7338) decreases much more than in medium B (cf. Rottenberg and Hashimoto (1986),Biochemistry,25, 1747–1755) at concentrations of fatty acid which equally stimulate the rate of respiration in state 4. Valinomycin and NaCl were both present in medium B and absent in medium A. However, in both media the pattern of the P/O ratio as a function of antimycin in the presence of a constant amount of palmitic acid or of FCCP shows similar behaviour. We conclude that in both media palmitic acid increases the membrane conductance to protons, but for unclear reasons the assay fails to measure the decline of in medium B. However, the increase in membrane conductance induced by palmitic acid does not quantitatively account for the stimulation of the rate of respiration.     

Mean-square hydrophobic moment for partially helical polypeptides

A mean-square helical hydrophobic moment, 〈h²〉, is defined for polypeptides in analogy to the mean-square dipole moment, 〈μ²〉, for polymer chains. For a freely jointed polymer chain, 〈μ²〉 is given by Σm, where m_i denotes the dipole moment associated with bond i. In the absence of any correlations in the hydrophobic moments of individual amino acid residues in the helix, 〈h²〉 is specified by ΣH, where H_i denotes the hydrophobicity of residue i. The tendency for correlations in orientations of residue hydrophobic moments in helices therefore dictates the size of 〈h²〉/〈H²〉, where 〈H²〉 denotes the average value of ΣH for all helices. The value of 〈h²〉/〈H²〉 will be greater than one in amphiphilic helices. A necessary prerequisite for this diagnostic usage of 〈h²〉/〈H²〉 is that the residue hydrophobic moment be oriented prependicular to the principal axis of the helix. Matrix-generation schemes are formulated that permit rapid evaluation of 〈h²〉 and 〈H²〉. The behavior of 〈h²〉/〈H²〉 is illustrated by calculations performed for model sequential copolypeptides.     

Structural aspects of the cytochromeb 6 f complex; structure of the lumen-side domain of cytochromef

The following findings concerning the structure of the cytochromeb ₆ f complex and its component polypeptides, cytb ₆, subunit IV and cytochromef subunit are discussed:

Targeting of tail‐anchored membrane proteins to subcellular organelles in Toxoplasma gondii

Proper protein localization is essential for critical cellular processes, including vesicle‐mediated transport and protein translocation. Tail‐anchored (TA) proteins are integrated into organellar membranes via the C‐terminus, orienting the N‐terminus towards the cytosol. Localization of TA proteins occurs posttranslationally and is governed by the C‐terminus, which contains the integral transmembrane domain (TMD) and targeting sequence. Targeting of TA proteins is dependent on the hydrophobicity of the TMD as well as the length and composition of flanking amino acid sequences. We previously identified an unusual homologue of elongator protein, Elp3, in the apicomplexan parasite Toxoplasma gondii as a TA protein targeting the outer mitochondrial membrane. We sought to gain further insight into TA proteins and their targeting mechanisms using this early‐branching eukaryote as a model. Our bioinformatics analysis uncovered 59 predicted TA proteins in Toxoplasma, 9 of which were selected for follow‐up analyses based on representative features. We identified novel TA proteins that traffic to specific organelles in Toxoplasma, including the parasite endoplasmic reticulum, mitochondrion, and Golgi apparatus. Domain swap experiments elucidated that targeting of TA proteins to these specific organelles was strongly influenced by the TMD sequence, including charge of the flanking C‐terminal sequence.      

The topology of the C-terminal sections of the NCX1 Na+/Ca2+ exchanger and the NCKX2 Na+/Ca2+-K+ exchanger

Mammalian Na⁺/Ca²⁺ (NCX) and Na⁺/Ca²⁺-K⁺ exchangers (NCKX) are polytopic membrane proteins that play critical roles in calcium homeostasis in many cells. Although hydropathy plots for NCX and NCKX are very similar, reported topological models for NCX1 and NCKX2 differ in the orientation of the three C-terminal transmembrane segments (TMS). NCX1 is thought to have 9 TMS and a re-entrant loop, whereas NCKX2 is thought to have 10 TMS. The current topological model of NCKX2 is very similar to the 10 membrane spanning helices seen in the recently reported crystal structure of NCX_MJ, a distantly related archaebacterial Na⁺/Ca²⁺ exchanger. Here we reinvestigate the orientation of the three C-terminal TMS of NCX1 and NCKX2 using mass-tagging experiments of substituted cysteine residues. Our results suggest that NCX1, NCKX2 and NCX_MJ all share the same 10 TMS topology.     

Sequence-specific 1H and 15N resonance assignments and secondary structure of GDP-bound human c-Ha-Ras protein in solution

Summary All the backbone ¹H and ¹⁵N magnetic resonances (except for Pro residues) of the GDP-bound form of a truncated human c-Ha-ras proto-oncogene product (171 amino acid residues, the Ras protein) were assigned by ¹⁵N-edited two-dimensional NMR experiments on selectively ¹⁵N-labeled Ras proteins in combination with three-dimensional NMR experiments on the uniformly ¹⁵N-labeled protein. The sequence-specific assignments were made on the basis of the nuclear Overhauser effect (NOE) connectivities of amide protons with preceding amide and/or Cprotons. In addition to sequential NOEs, vicinal spin coupling constants for amide protons and C protons and deuterium exchange rates of amide protons were used to characterize the secondary structure of the GDP-bound Ras protein; six strands and five helices were identified and the topology of these elements was determined. The secondary structure of the Ras protein in solution was mainly consistent with that in crystal as determined by X-ray analyses. The deuterium exchange rates of amide protons were examined to elucidate the dynamic properties of the secondary structure elements of the Ras protein in solution. In solution, the -sheet structure in the Ras protein is rigid, while the second helix (A66-R73) is much more flexible, and the first and fifth helices (S17-124 and V152-L171) are more rigid than other helices. Secondary structure elements at or near the ends of the effector-region loop were found to be much more flexible in solution than in the crystalline state.     

A fusogenic protein from rat brain microsomal membranes: Partial purification and reconstitution into liposomes

The procedures for purification and reconstitution of rat brain microsomal membrane protein that causes fusion of liposomes at acidic pH are described. A 1,860-fold purification was achieved, starting from the detergent-solubilized microsomal membranes. The fusion process was assayed spectrofluorimetrically by monitoring the formation of terbium-dipicolinic acid complex (Wilschut, J. et al. 1980. Biochemistry 19:6011–6021) evoked by the protein after mixing of two populations of liposomes. The fusogenic activity of the protein inserted into the membrane of Tb³⁺-containing vesicles was found to be strongly dependent on phospholipid composition and was higher in vesicles enriched with exogenous phosphatidylserine, phosphatidylglycerol and phosphatidylethanolamine than in those prepared with an excess of phosphatidylcholine. The vesicles enriched in negatively charged phospholipids were bound to Concanavalin A coupled to Sepharose-4B and could be released from this column only in the presence of a high concentration of -methylmannopyranoside and detergent, indicating a glycoprotein nature of the fusogenic protein. Furthermore, these data show that protein inserted into membrane has its oligosaccharide chains exposed to the environment.Mr. Carlo Ricci is thanked for his skillful technical assistance. This work was supported by a grant from the Ministry of Education, Rome, Italy.     

Hydrophobicity and prediction of the secondary structure of membrane proteins and peptides.

Reliability of the hydropathy method to predict the formation of membrane-spanning alpha-helices by integral membrane proteins and peptides whose structure is known from X-ray crystallography is analysed. It is shown that Kyte-Doolittle hydropathy plots do not predict accurately 22 transmembrane alpha-helices in the reaction centres (RC) of the photosynthetic bacteria Rhodopseudomonas viridis and Rhodobacter sphaeroides (R-26). The accuracy of prediction for these proteins was improved using an optimised Kyte-Doolittle hydrophobicity scale. However, this hydrophobicity scale did not improve the predictions for the alphabeta-peptides of the B800-850 (LH2) complexes of the photosynthetic bacteria Rhodopseudomonas acidophila and Rhodospirillum molischianum, which were excluded from the optimisation procedure. The best and worst predictions of membrane-spanning alpha-helices for the RC proteins and LH2 peptides, respectively, were obtained with a propensity scale (PRC) calculated from the amino acid sequences and X-ray data for the RC proteins. A propensity scale (PLH) obtained using the amino acid sequences and X-ray data for the alphabeta-peptides of the LH2 complexes did not give an acceptable prediction of the transmembrane segments in the LH2 peptides; moreover, it markedly contradicted the PRC scale. Amino acids have been concluded to have no significant preference to localisation in transmembrane segments. Therefore, the predictive ability of the hydropathy methodology appears to be limited: the number of transmembrane segments can be correctly calculated for the best case only, and the lengths and positions of membrane-spanning alpha-helices in a protein amino acid sequence can not be predicted exactly.     

Spasmin and a Putative Spasmin Binding Protein(s) Isolated from Solubilized Spasmonemes1

ABSTRACT The amino acid composition and hydrophobicity scale (hydropathy) of calcium-binding proteins contained in the contractile spasmoneme of Carchesium polypinum was compared with other calcium-binding proteins from eukaryotes. Spasmins which may hind at most 4 calcium ions simultaneously and initiate spasmoneme contraction cooperatively belong to a super family of proteins including; centrin/caltractin and calmodulin. Based on chemical modification of tryptophan and methionine, these residues are involved in contraction but the spasmin proteins contain little or none of these amino acids. Based on this evidence, it is suggested that another, non-calcium binding protein(s) is involved in spasmoneme contraction.     

Thin Filament Proteins and Thin Filament-Linked Regulation of Vertebrate Muscle Contractio

Abstract

The major intrinsic protein (MIP) of the bovine lens fiber cell membrane was the first member of the MIP family of proteins to be sequenced and characterized. It is probably a homotetramer with transmembrane channel activity that plays a role in lens biogenesis or maintenance. The polypeptide chain of each subunit may span the membrane six times, and both the N- and C-termini face the cell cytoplasm. Eighteen sequenced or partially sequenced proteins from bacteria, yeast, plants, and animals have now been shown to be members of the MIP family. These proteins appear to function in (1) metazoan development and neurogenesis (MIP and BIB), (2) water transport across the human erythrocyte membrane (ChIP), (3) communication between host plant cells and symbiotic nitrogen-fixing bacteria (NOD), (4) transport across the tonoplast membrane during plant seed development (α-TIP), (5) water stress-induced resistance to desiccation in plants (Wsi-TIP), (6) suppression of a genetic growth defect on fermentable sugars in yeast (FPS1), and (7) transport of glycerol across bacterial cell membranes (GlpF). One other sequenced member of the MIP family (ORF1 of Lactococcus lactis) has no known physiological function. The biochemical functions of the eukaryotic proteins are not well established.

Computer analyses have revealed that the first and second halves of all MTP family proteins probably arose by a tandem, intragenic, duplication event. Thus, the primary structure of putative transmembrane helices 1 to 3 is similar to that of putative transmembrane helices 4 to 6 even though they are of opposite orientation in the membrane. Among the most conserved residues in these two repeated halves are a membrane-embedded glutamate (E) in helices 1 and 4, an asparagine-proline-alanine (NPA) sequence in the loops between helices 2 and 3 (cytoplasmically localized) and helices 5 and 6 (extracellularly localized), and a glycine within helices 3 and 6. Statistical analyses suggest that the two halves of these proteins have evolved to serve distinct functions: the first half is more important for the generalized or common functions of these proteins, while the second half of these proteins is more differentiated to provide specific or dissimilar functions of the proteins. The apparent origin of MIP family proteins by duplication of a three-spanner precursor protein suggests an evolutionary origin distinct from other transport proteins with six transmembrane spanners. Based on the phylogenetic tree for the 18 sequenced members of the MTP family, we propose that a single, primordial gene arose in prokaryotes shortly before the emergence of eukaryotes, mat this gene was vertically transmitted to the principal eukaryotic kingdoms, and that subsequent gene duplication and divergence events gave rise to kingdom-related subfamilies or clusters of the MIP family.     

Aliphatic 1H and 13C resonance assignments for the 26-10 antibody VL domain derived from heteronuclear multidimensional NMR spectroscopy

Summary Extensive ¹H and ¹³C assignments have been obtained for the aliphatic resonances of a uniformly ¹³C-and ¹⁵N-labeled recombinant V_L domain from the anti-digoxin antibody 26-10. Four-dimensional triple resonance NMR data acquired with the HNCAHA and HN(CO)CAHA pulse sequences [Kay et al. (1992) J. Magn. Reson., 98, 443–450] afforded assignments for the backbone H^N, N, H and C resonances. These data confirm and extend H^N, N and H assignments derived previously from three-dimensional ¹H-¹⁵N NMR studies of uniformly ¹⁵N-labeled V_L domain [Constantine et al. (1992), Biochemistry, 31, 5033–5043]. The identified H and C resonances provided a starting point for assigning the side-chain aliphatic ¹H and ¹³C resonances using three-dimensional HCCH-COSY and HCCH-TOCSY experiments [Clore et al. (1990), Biochemistry, 29, 8172–8184]. The C and C chemical shifts are correlated with the V_L domain secondary structure. The extensive set of side-chain assignments obtained will allow a detailed comparison to be made between the solution structure of the isolated V_L domain and the X-ray structure of the V_L domain within the 26–10 Fab.     

Not all transmembrane helices are born equal: Towards the extension of the sequence homology concept to membrane proteins

Background

Sequence homology considerations widely used to transfer functional annotation to uncharacterized protein sequences require special precautions in the case of non-globular sequence segments including membrane-spanning stretches composed of non-polar residues. Simple, quantitative criteria are desirable for identifying transmembrane helices (TMs) that must be included into or should be excluded from start sequence segments in similarity searches aimed at finding distant homologues.

Results

We found that there are two types of TMs in membrane-associated proteins. On the one hand, there are so-called simple TMs with elevated hydrophobicity, low sequence complexity and extraordinary enrichment in long aliphatic residues. They merely serve as membrane-anchoring device. In contrast, so-called complex TMs have lower hydrophobicity, higher sequence complexity and some functional residues. These TMs have additional roles besides membrane anchoring such as intra-membrane complex formation, ligand binding or a catalytic role. Simple and complex TMs can occur both in single- and multi-membrane-spanning proteins essentially in any type of topology. Whereas simple TMs have the potential to confuse searches for sequence homologues and to generate unrelated hits with seemingly convincing statistical significance, complex TMs contain essential evolutionary information.

Conclusion

For extending the homology concept onto membrane proteins, we provide a necessary quantitative criterion to distinguish simple TMs (and a sufficient criterion for complex TMs) in query sequences prior to their usage in homology searches based on assessment of hydrophobicity and sequence complexity of the TM sequence segments.

Reviewers

This article was reviewed by Shamil Sunyaev, L. Aravind and Arcady Mushegian.     

Modeling 2hJiso(N, N) in nucleic acid base pairs: Ab initio characterization of the 2hJ(N, N) tensor in the methyleneimine dimer as a function of hydrogen bond geometry

The observation of ^2h J _iso(N, N) coupling has prompted considerable interest in this phenomenon from experimentalists and theoreticians due to the potential these couplings hold for the determination of secondary and tertiary structure in biologically important molecules. Here, we present an ab initio (MCSCF) study of the complete ^2h J(N, N) tensor for a model methyleneimine dimer system as a function of (i) the N-N separation, r _NN, and (ii) the hydrogen bond angle, . This simple system models the ^2h J(N, N) tensor of nucleic acid base pairs. Results indicate that although the Fermi-contact mechanism dominates ^2h J _iso(N, N), the coupling tensor is anisotropic due to contributions from the Fermi-contact spin-dipolar cross term. The variation in ^2h J _iso(N, N) as a function of r _NN is fit to an exponential decay. The influence of on the coupling constant is less pronounced but must be considered if experimental coupling constants are to be used for quantitative structure determination. Our results for this simple model system demonstrate that ^2h J _iso(N, N) is a valuable probe of hydrogen bonding in nucleic acid base pairs.     

Characterization of the immunochemical reactivity of fibrinogen fragments by competitive radioimmunoassay: An improved method of analysis

Published results on the immunochemical reactivities of fibrinogen and fibrinogen fragments with fibrinogen-elicited antibodies that had been fractionated on the basis of preferential interaction with A [Nagy, J. A., Meinwald, Y. C., and Scheraga, H. A. (1982),Biochemistry 21, 1794–1806] and B [Nagy, J. A., Meinwald, Y. C., and Scheraga, H. A. (1985)Biochemistry 24, 882–887] peptides of this bivalent antigen have been reinterpreted. First, the multivalent counterpart of the Scatchard analysis has been used to determine the intrinsic association constant for the interaction of antibody with [¹²⁵I]fibrinogen, the radiolabeled ligand used in subsequent competitive binding studies. Second, the corresponding affinity constant for native fibrinogen has been evaluated from the relevant competitive radioimmunoassays by means of a quantitative analysis that takes into account the bivalency of both the radiolabeled and native fibrinogen molecules. Finally, affinity constants for the interactions of various fibrinogen fragments with antibody are also obtained by the procedure, and their magnitudes rationalized in terms of the equilibrium coexistence of unreactive (disordered) and native (functional) states of the fibrinogen peptides.     

The carboxy-terminal region of human interleukin-5 is essential for maintenance of tertiary structure but not for dimerization

The C-terminal region of interleukin-5 has previously been suggested to be important for biological activity [Mackenzieet al., (1991),Mol. Immunol. 28, 155–158; Kodamaet al. (1991),Biochem. Biophys. Res. Commun. 178, 514–519]. We have investigated this region by making a series of truncation mutants. The proteins were expressed inEscherichia coli, purified from inclusion bodies, and were able to refold with the disulfide homodimeric topology typical of interleukin-5. Analysis of the truncated carboxy-terminal proteins in an interleukin-5-dependent proliferation assay on TF-1 cells showed a rapid loss of activity as the C-terminal was shortened by more than two amino acids. This loss of biological activity correlated with a drop in binding affinity to both the chain of the receptor and the high-affinity complex consisting of the and subunits. Analysis of the proteins by¹H-NMR showed that the truncated mutants have higher exchange rates with solvent, indicating a less rigid structure. The carboxy-terminal region is therefore necessary to maintain the stability of the four-helix bundle and to orient correctly the important residues of the fourth helix. Inspection of the structure determined by X-ray crystallography shows that Trp-110 acts as the major residue in anchoring the fourth helix.