A combination of compositional index and genetic algorithm for predicting transmembrane helical segments

Transmembrane helix (TMH) topology prediction is becoming a focal problem in bioinformatics because the structure of TM proteins is difficult to determine using experimental methods. Therefore, methods that can computationally predict the topology of helical membrane proteins are highly desirable. In this paper we introduce TMHindex, a method for detecting TMH segments using only the amino acid sequence information. Each amino acid in a protein sequence is represented by a Compositional Index, which is deduced from a combination of the difference in amino acid occurrences in TMH and non-TMH segments in training protein sequences and the amino acid composition information. Furthermore, a genetic algorithm was employed to find the optimal threshold value for the separation of TMH segments from non-TMH segments. The method successfully predicted 376 out of the 378 TMH segments in a dataset consisting of 70 test protein sequences. The sensitivity and specificity for classifying each amino acid in every protein sequence in the dataset was 0.901 and 0.865, respectively. To assess the generality of TMHindex, we also tested the approach on another standard 73-protein 3D helix dataset. TMHindex correctly predicted 91.8% of proteins based on TM segments. The level of the accuracy achieved using TMHindex in comparison to other recent approaches for predicting the topology of TM proteins is a strong argument in favor of our proposed method. Availability: The datasets, software together with supplementary materials are available at: http://faculty.uaeu.ac.ae/nzaki/TMHindex.htm.     

A novel method for predicting transmembrane segments in proteins based on a statistical analysis of the SwissProt database: the PRED-TMR algorithm.

We present a novel method that predicts transmembrane domains in proteins using solely information contained in the sequence itself. The PRED-TMR algorithm described, refines a standard hydrophobicity analysis with a detection of potential termini ('edges', starts and ends) of transmembrane regions. This allows one both to discard highly hydrophobic regions not delimited by clear start and end configurations and to confirm putative transmembrane segments not distinguishable by their hydrophobic composition. The accuracy obtained on a test set of 101 non-homologous transmembrane proteins with reliable topologies compares well with that of other popular existing methods. Only a slight decrease in prediction accuracy was observed when the algorithm was applied to all transmembrane proteins of the SwissProt database (release 35). A WWW server running the PRED-TMR algorithm is available at http://o2.db.uoa. gr/PRED-TMR/     

The position of the Gly-xxx-Gly motif in transmembrane segments modulates dimer affinity.

Although the intrinsic low solubility of membrane proteins presents challenges to their high-resolution structure determination, insight into the amino acid sequence features and forces that stabilize their folds has been provided through study of sequence-dependent helix-helix interactions between single transmembrane (TM) helices. While the stability of helix-helix partnerships mediated by the Gly-xxx-Gly (GG4) motif is known to be generally modulated by distal interfacial residues, it has not been established whether the position of this motif, with respect to the ends of a given TM segment, affects dimer affinity. Here we examine the relationship between motif position and affinity in the homodimers of 2 single-spanning membrane protein TM sequences: glycophorin A (GpA) and bacteriophage M13 coat protein (MCP). Using the TOXCAT assay for dimer affinity on a series of GpA and MCP TM segments that have been modified with either 4 Leu residues at each end or with 8 Leu residues at the N-terminal end, we show that in each protein, centrally located GG4 motifs are capable of stronger helix-helix interactions than those proximal to TM helix ends, even when surrounding interfacial residues are maintained. The relative importance of GG4 motifs in stabilizing helix-helix interactions therefore must be considered not only in its specific residue context but also in terms of the location of the interactive surface relative to the N and C termini of alpha-helical TM segments.     

Ion channel formation by synthetic transmembrane segments of the inhibitory glycine receptor--a model study.

The inhibitory glycine receptor (GlyR) of rat spinal cord contains an intrinsic transmembrane channel mediating agonist-gated anion flux. Here, synthetic peptides modelled after the predicted transmembrane domains M2 and M4 of its ligand-binding subunit were incorporated into lipid vesicle membranes and black lipid bilayers to analyze their channel forming capabilities. Both types of peptides prohibited the establishment of, or dissipated, preexisting transmembrane potentials in the vesicle system. Incorporation of peptide M2 into the black lipid bilayer elicited randomly gated single channel events with various conductance states and life-times. Peptide M4 increased the conductance of the bilayer without producing single channels. Exchange of the terminal arginine residues of peptide M2 by glutamate resulted in a significant shift towards cation selectivity of the respective channels as compared to peptide M2. In conclusion, the peptide channels observed differed significantly from native GlyR in both conductivity and ion-selectivity indicating that individual synthetic transmembrane segments are not sufficient to mimic a channel protein composed of subunits with multiple transmembrane segments.     

Perspective: peptides as mimics of transmembrane segments in proteins.

Peptide-based approaches to protein structure within membranes have proven enormously valuable. When one focusses on the detailed manner through which membrane proteins actually traverse the cell bilayer, a simple observation emerges: helical peptide segments of 20 amino acids each constitute the only tangible connection between the inside and outside of the cell. Thus, a major step towards understanding the key relationships between biological function and membrane protein structure can be taken through characterization, by composition, sequence, chain length, hydrophobicity and conformation, of hydrophobic peptides designed as mimics of transmembrane segments.     

Insertion of an arginine residue into the transmembrane segments corrects protein misfolding

Deletion of Phe-508 (DeltaF508) in cystic fibrosis transmembrane conductance regulator causes cystic fibrosis because of misfolding of the protein. P-glycoprotein (P-gp) containing the equivalent mutation (DeltaY490) is also misfolded but can be rescued with drug substrates. Whether rescue is due to direct binding of drug substrate to the transmembrane (TM) segments or to indirect effects on cellular protein folding pathways is still controversial. P-gp-drug substrate interactions likely involve hydrogen bonds. If the mechanism of drug rescue involves changes to TM packing then we should be able to identify suppressor mutations in the TM segments that can mimic the drug rescue effects. We predicted that an arginine residue in the TM segments predicted to line the drug-binding pocket of P-gp (I306(TM5) or F343(TM6)) might suppress DeltaY490 P-gp protein misfolding because it has the highest propensity to form hydrogen bonds. We show that R306(TM5) or R343(TM6) increased the relative amount of mature DeltaY490 P-gp by 6-fold. Most other changes to Ile-306 or Phe-343 did not enhance maturation of DeltaY490 P-gp. The I306R mutant also promoted maturation of misprocessed mutants that had mutations in the second nucleotide-binding domain (L1260A), the cytoplasmic loops (G251V, F804A), the linker region (P709A), or in TM segments (G300V, G722A). These results show that arginine residues in the TM domains can mimic the drug rescue effects and are effective suppressor mutations for processing mutations located throughout the molecule.     

Membrane orientation of carboxyl-terminal half P-glycoprotein: topogenesis of transmembrane segments.

Multiple topological orientations of the carboxyl-terminal half of P-glycoprotein have been observed. One orientation is consistent with the hydropathy-predicted model and contains six transmembrane (TM)-spanning regions. In another orientation, the cytoplasmic-predicted loop between TM8 and TM9 is extracellular and glycosylated. In support of this "alternative" topology, TM8 was previously established to function as a signal-anchor sequence to insert with its amino-terminal end in the cytoplasm and the carboxyl-terminal end in the extracytoplasmic space. However, it is unclear how downstream TM segments fold in the membrane when TM8 functions as a signal-anchor sequence. Here, we created several chimeric Pgp molecules to examine the membrane insertion of TM segments 9 and 10 using a cell-free system. We found that TM9 functions as a stop-transfer sequence when following the signal-anchor sequence, TM8. However, the stop-transfer activity of TM9 depends on the presence of TM10. In the absence of TM10, TM9 partially translocated across the membrane into the endoplasmic reticulum lumen. In contrast, TM9 efficiently stopped the translocation event of the nascent chain in the presence of TM10. Our results suggest that the membrane insertion of TM8 and TM9 establishes the extracellular loop between TM8 and TM9. Formation of this loop apparently involves the interactions between Pgp TM segments, which facilitate proper folding of the Pgp carboxyl-terminal half.     

The DeltaF508 mutation disrupts packing of the transmembrane segments of the cystic fibrosis transmembrane conductance regulator

The most common mutation in cystic fibrosis (deletion of Phe-508 in the first nucleotide binding domain (DeltaF508)) in the cystic fibrosis transmembrane conductance regulator (CFTR) causes retention of the mutant protein in the endoplasmic reticulum. We previously showed that the DeltaF508 mutation causes the CFTR protein to be retained in the endoplasmic reticulum in an inactive and structurally altered state. Proper packing of the transmembrane (TM) segments is critical for function because the TM segments form the chloride channel. Here we tested whether the DeltaF508 mutation altered packing of the TM segments by disulfide cross-linking analysis between TM6 and TM12 in wild-type and DeltaF508 CFTRs. These TM segments were selected because TM6 appears to line the chloride channel, and cross-linking between these TM segments has been observed in the CFTR sister protein, the multidrug resistance P-glycoprotein. We first mapped potential contact points in wild-type CFTR by cysteine mutagenesis and thiol cross-linking analysis. Disulfide cross-linking was detected in CFTR mutants M348C(TM6)/T1142C(TM12), T351C(TM6)/T1142C(TM12), and W356C(TM6)/W1145C(TM12) in a wild-type background. The disulfide cross-linking occurs intramolecularly and was reducible by dithiothreitol. Introduction of the DeltaF508 mutation into these cysteine mutants, however, abolished cross-linking. The results suggest that the DeltaF508 mutation alters interactions between the TM domains. Therefore, a potential target to correct folding defects in the DeltaF508 mutant of CFTR is to identify compounds that promote correct folding of the TM domains.     

Tertiary interactions between the fifth and sixth transmembrane segments of rhodopsin.

We have used cysteine scanning mutagenesis and disulfide cross-linking in a split rhodopsin construct to investigate the secondary structure and tertiary contacts of the fifth (TM5) and sixth (TM6) transmembrane segments of rhodopsin. Using a simple increase in pH to promote disulfide bond formation, three cross-links between residues on the extracellular side of TM5 (at positions 198, 200, and 204) and TM6 (at position 276) have been identified and characterized. The helical pattern of cross-linking observed indicates that the fifth transmembrane helix extends through residue 200 but does not include residue 198. Rhodopsin mutants containing these disulfides demonstrate nativelike absorption spectra and light-dependent activation of transducin, suggesting that large movements on the extracellular side of TM5 with respect to TM6 are not required for receptor activation.     

Peptides as transmembrane segments: decrypting the determinants for helix-helix interactions in membrane proteins

Although the structural analysis of membrane proteins is advancing, an understanding of the basic principles that underlie their folding and assembly remains limited because of the high insolubility intrinsic to these molecules and concomitant challenges in obtaining crystals. Fortunately, from an experimental standpoint, membrane protein folding can be approximated as the rigid-body docking of pre-formed alpha-helical transmembrane segments one with another to form the final functional protein structure. Peptides derived from the sequences of native alpha-helical transmembrane segments and those that mimic their properties are therefore valuable in the experimental evaluation of protein folding within the membrane. Here we present an overview of the progress made in our laboratory and elsewhere in using peptide models toward defining the sequence requirements and forces stabilizing membrane protein folds.     

Isoelectric focusing management: an investigation for salt interference and an algorithm for optimization

Two-dimensional electrophoresis (2-DE) has evolved into a robust separation technique in proteomic research. However, one of the major challenges in 2-DE experiments, the reproducibility of the first dimensional electrophoresis (IEF), has remained unsolved. It is well-known that the quality of IEF experiments is significantly affected by the salt interference. Nevertheless, the interference mechanisms of salts in IEF have never been systematically investigated. In this study, we comprehensively investigated the interference effects in IEF due to various kinds of simple and buffer salts in protein samples. Two interference schemes were proposed accordingly to elucidate the interference mechanisms of salts in IEF. Furthermore, to increase the reproducibility of IEF, we proposed that conductivity measurement is a feasible method to assess the salt content of 2-DE samples and developed an algorithm to predict the optimal total volt-hours (Vh) required for protein focusing in IEF. The developed algorithm had been evaluated under various IEF conditions for a variety of 2-DE samples and proven to be a reliable guide. In sum, information disclosed in this study should be of use for increasing the reproducibility and thus the applicability of 2-DE in current proteomics.     

Getting greasy: how transmembrane polypeptide segments integrate into the lipid bilayer

Many integral membrane proteins use the same translocation machinery for membrane insertion as secretory proteins use to get across the membrane. This requires that transmembrane segments can be discriminated from other parts of the protein during membrane translocation, and further requires that the transmembrane segments can be moved laterally out of the translocation channel into the surrounding lipid. The molecular basis for this remarkable intramembraneous sorting event is a major focus of current studies of membrane protein biogenesis.     

Branch-and-terminate: a combinatorial optimization algorithm for protein design.

BACKGROUND: Several deterministic and stochastic combinatorial optimization algorithms have been applied to computational protein design and homology modeling. As structural targets increase in size, however, it has become necessary to find more powerful methods to address the increased combinatorial complexity. RESULTS: We present a new deterministic combinatorial search algorithm called 'Branch-and-Terminate' (B&T), which is derived from the Branch-and-Bound search method. The B&T approach is based on the construction of an efficient but very restrictive bounding expression, which is used for the search of a combinatorial tree representing the protein system. The bounding expression is used both to determine the optimal organization of the tree and to perform a highly effective pruning procedure named 'termination'. For some calculations, the B&T method rivals the current deterministic standard, dead-end elimination (DEE), sometimes finding the solution up to 21 times faster. A more significant feature of the B&T algorithm is that it can provide an efficient way to complete the optimization of problems that have been partially reduced by a DEE algorithm. CONCLUSIONS: The B&T algorithm is an effective optimization algorithm when used alone. Moreover, it can increase the problem size limit of amino acid sidechain placement calculations, such as protein design, by completing DEE optimizations that reach a point at which the DEE criteria become inefficient. Together the two algorithms make it possible to find solutions to problems that are intractable by either algorithm alone.     

N-tail translocation in a eukaryotic polytopic membrane protein: synergy between neighboring transmembrane segments.

We have used the natural N-glycosylation site in the N-tail of cig30, a eukaryotic polytopic membrane protein, as a marker for N-tail translocation across the microsomal membrane. Analysis of C-terminally truncated cig30 constructs reveals that the first transmembrane segment is sufficient for translocation of the wild-type N-tail; in contrast, in a mutant with four arginines introduced into the N-tail the second transmembrane segment is also required for efficient N-tail translocation. Our observations imply a non-sequential assembly mechanism in which the ultimate location of the N-tail relative to the membrane may depend on more than one transmembrane segment.     

kPROT: a knowledge-based scale for the propensity of residue orientation in transmembrane segments. Application to membrane protein structure prediction

Modeling of integral membrane proteins and the prediction of their functional sites requires the identification of transmembrane (TM) segments and the determination of their angular orientations. Hydrophobicity scales predict accurately the location of TM helices, but are less accurate in computing angular disposition. Estimating lipid-exposure propensities of the residues from statistics of solved membrane protein structures has the disadvantage of relying on relatively few proteins. As an alternative, we propose here a scale of knowledge-based Propensities for Residue Orientation in Transmembrane segments (kPROT), derived from the analysis of more than 5000 non-redundant protein sequences. We assume that residues that tend to be exposed to the membrane are more frequent in TM segments of single-span proteins, while residues that prefer to be buried in the transmembrane bundle interior are present mainly in multi-span TMs. The kPROT value for each residue is thus defined as the logarithm of the ratio of its proportions in single and multiple TM spans. The scale is refined further by defining it for three discrete sections of the TM segment; namely, extracellular, central, and intracellular. The capacity of the kPROT scale to predict angular helical orientation was compared to that of alternative methods in a benchmark test, using a diversity of multi-span alpha-helical transmembrane proteins with a solved 3D structure. kPROT yielded an average angular error of 41 degrees, significantly lower than that of alternative scales (62 degrees -68 degrees ). The new scale thus provides a useful general tool for modeling and prediction of functional residues in membrane proteins. A WWW server (http://bioinfo.weizmann.ac.il/kPROT) is available for automatic helix orientation prediction with kPROT.     

Evidence for similar function of transmembrane segments in receptor and membrane-anchored proteins

Transmembrane (TM) regions of receptor proteins should have unique structural and/or chemical characteristics if these regions contain residues functional in TM signal transduction. However, in a survey of the membrane-occurring residues in 37 integral membrane proteins, we found that amino acid compositions of TM regions of receptor proteins (n = 11) could not be distinguished statistically from corresponding regions of membrane-anchored proteins (e.g., recognition molecules) with a functional external domain attached to a single hydrophobic membrane-spanning anchor segment (n = 16). TM regions in both categories of proteins differed from the compositions of TM regions in membrane-transport proteins (n = 10). The analysis implies that TM regions in receptor proteins may function mainly to anchor (and position) receptors in their cellular membranes, and therefore residues in receptors that participate in signal transduction need not be restricted to these regions. In addition to mechanisms involving receptor aggregation, ligand-activated conformational perturbation of a receptor external aqueous domain, resulting in membrane penetration of hydrophobic segment(s) of this domain to produce intramembranous contact with its cytoplasmic domain, is hypothesized as a further possible mode of signal transduction.     

Tertiary interactions between transmembrane segments 3 and 5 near the cytoplasmic side of rhodopsin.

Previous studies [Yu, H., Kono, M., and Oprian, D. D. (1999) Biochemistry 38, xxxx-xxxx] using split receptors and disulfide cross-linking have shown that native cysteines 140 and 222 on the cytoplasmic side of transmembrane segments (TM) 3 and 5 of rhodopsin, respectively, can cross-link to each other upon treatment with the oxidant Cu(phen)3(2+). In this paper we show that although the 140-222 cross-link does not affect the spectral properties of rhodopsin, it completely and reversibly inactivates the ability of the receptor to activate transducin. Following on this lead we further investigate the cytoplasmic region of TM3 and TM5 and identify three additional pairs of residues that when changed to Cys are capable of forming disulfide cross-links in the protein: 140/225, 136/222, and 136/225. These disulfides are able to form without addition of the Cu(phen)3(2+) oxidant. Similar to the 140-222 cross-link, none of the additional disulfides affect the spectral properties of rhodopsin. Also like the 140-222 bond, the 136-222 disulfide completely and reversibly inactivates the light-dependent activation of transducin by the receptor. In contrast, the 140-225 and 136-225 disulfides have no effect on the ability of rhodopsin to activate transducin. The pattern of cross-linking observed in Cys and disulfide scans of the protein is consistent with helical secondary structure in TM3 from 130 to 142 and in TM5 from 218 to 225.     

Cooperative combinatorial optimization: evolutionary computation case study

This paper presents a formalization of the notion of cooperation and competition of multiple systems that work toward a common optimization goal of the population using evolutionary computation techniques. It is proved that evolutionary algorithms are more expressive than conventional recursive algorithms, such as Turing machines. Three classes of evolutionary computations are introduced and studied: bounded finite, unbounded finite, and infinite computations. Universal evolutionary algorithms are constructed. Such properties of evolutionary algorithms as completeness, optimality, and search decidability are examined. A natural extension of evolutionary Turing machine (ETM) model is proposed to properly reflect phenomena of cooperation and competition in the whole population.     

Substrate-induced conformational changes in the transmembrane segments of human P-glycoprotein. Direct evidence for the substrate-induced fit mechanism for drug binding

The human multidrug resistance P-glycoprotein (P-gp, ABCB1) is quite promiscuous in that it can transport a broad range of structurally diverse compounds out of the cell. We hypothesized that the transmembrane (TM) segments that constitute the drug-binding site are quite mobile such that drug binding occurs through a "substrate-induced fit" mechanism. Here, we used cysteine-scanning mutagenesis and oxidative cross-linking to test for substrate-induced changes in the TM segments. Pairs of cysteines were introduced into a Cys-less P-gp and the mutants treated with oxidant (copper phenanthroline) in the presence or absence of various drug substrates. We show that cyclosporin A promoted cross-linking between residues P350C(TM6)/G939C(TM11), while colchicine and demecolcine promoted cross-linking between residues P350C(TM6)/V991C(TM12). Progesterone promoted cross-linking between residues P350C(TM6)/A935C(TM11), P350C(TM6)/G939C(TM11), as well as between residues P350C(TM6)/V991C(TM12). Other substrates such as vinblastine, verapamil, cis-(Z)-flupenthixol or trans-(E)-flupenthixol did not induce cross-linking at these sites. These results provide direct evidence that the packing of the TM segments in the drug-binding site is changed when P-gp binds to a particular substrate. The induced-fit mechanism explains how P-gp can accommodate a broad range of compounds.