Computational differentiation of N-terminal signal peptides and transmembrane helices

Signal peptides and transmembrane helices both contain a stretch of hydrophobic amino acids. This common feature makes it difficult for signal peptide and transmembrane helix predictors to correctly assign identity to stretches of hydrophobic residues near the N-terminal methionine of a protein sequence. The inability to reliably distinguish between N-terminal transmembrane helix and signal peptide is an error with serious consequences for the prediction of protein secretory status or transmembrane topology. In this study, we report a new method for differentiating protein N-terminal signal peptides and transmembrane helices. Based on the sequence features extracted from hydrophobic regions (amino acid frequency, hydrophobicity, and the start position), we set up discriminant functions and examined them on non-redundant datasets with jackknife tests. This method can incorporate other signal peptide prediction methods and achieve higher prediction accuracy. For Gram-negative bacterial proteins, 95.7% of N-terminal signal peptides and transmembrane helices can be correctly predicted (coefficient 0.90). Given a sensitivity of 90%, transmembrane helices can be identified from signal peptides with a precision of 99% (coefficient 0.92). For eukaryotic proteins, 94.2% of N-terminal signal peptides and transmembrane helices can be correctly predicted with coefficient 0.83. Given a sensitivity of 90%, transmembrane helices can be identified from signal peptides with a precision of 87% (coefficient 0.85). The method can be used to complement current transmembrane protein prediction and signal peptide prediction methods to improve their prediction accuracies. (C) 2003 Elsevier Inc. All rights reserved.

SVMtm: Support vector machines to predict transmembrane segments

A new method has been developed for prediction of transmembrane helices using support vector machines. Different coding schemes of protein sequences were explored, and their performances were assessed by crossvalidation tests. The best performance method can predict the transmembrane helices with sensitivity of 93.4% and precision of 92.0%. For each predicted transmembrane segment, a score is given to show the strength of transmembrane signal and the prediction reliability. In particular, this method can distinguish transmembrane proteins from soluble proteins with an accuracy of similar to99%. This method can be used to complement current transmembrane helix prediction methods and can be Used for consensus analysis of entire proteomes . The predictor is located at http://genet.imb.uq.edu.au/predictors/ SVMtm. (C) 2004 Wiley Periodicals, Inc.

Predicting the solvent accessibility of transmembrane residues from protein sequence

In this study, we propose a novel method to predict the solvent accessible surface areas of transmembrane residues. For both transmembrane alpha-helix and beta-barrel residues, the correlation coefficients between the predicted and observed accessible surface areas are around 0.65. On the basis of predicted accessible surface areas, residues exposed to the lipid environment or buried inside a protein can be identified by using certain cutoff thresholds. We have extensively examined our approach based on different definitions of accessible surface areas and a variety of sets of control parameters. Given that experimentally determining the structures of membrane proteins is very difficult and membrane proteins are actually abundant in nature, our approach is useful for theoretically modeling membrane protein tertiary structures, particularly for modeling the assembly of transmembrane domains. This approach can be used to annotate the membrane proteins in proteomes to provide extra structural and functional information.

Structural plasticity in MS channels

Gating of the mechanosensitive channel MscS involves cooperative action of glycine and alanine residues along the pore-lining transmembrane helix. Opening of the channel is facilitated by an iris-like rotation and tilt of the pore-lining helices. Site-directed mutagenesis indicates that substantial structural plasticity can be tolerated by MscS without impairing its function.

Prediction of Golgi Type II membrane proteins based on their transmembrane domains

Motivation: A major issue in cell biology today is how distinct intracellular regions of the cell, like the Golgi Apparatus, maintain their unique composition of proteins and lipids. The cell differentially separates Golgi resident proteins from proteins that move through the organelle to other subcellular destinations. We set out to determine if we could distinguish these two types of transmembrane proteins using computational approaches. Results: A new method has been developed to predict Golgi membrane proteins based on their transmembrane domains. To establish the prediction procedure, we took the hydrophobicity values and frequencies of different residues within the transmembrane domains into consideration. A simple linear discriminant function was developed with a small number of parameters derived from a dataset of Type II transmembrane proteins of known localization. This can discriminate between proteins destined for Golgi apparatus or other locations (post-Golgi) with a success rate of 89.3% or 85.2%, respectively on our redundancy-reduced data sets.

Structure of the HERG K+ channel S5P extracellular linker - Role of an amphipathic alpha-helix in C-type inactivation

The HERG K+ channel has very unusual kinetic behavior that includes slow activation but rapid inactivation. These features are critical for normal cardiac repolarization as well as in preventing lethal ventricular arrhythmias. Mutagenesis studies have shown that the extracellular peptide linker joining the fifth transmembrane domain to the pore helix is critical for rapid inactivation of the HERG K+ channel. This peptide linker is also considerably longer in HERG K+ channels, 40 amino acids, than in most other voltage-gated K+ channels. In this study we show that a synthetic 42-residue peptide corresponding to this linker region of the HERG K+ channel does not have defined structural elements in aqueous solution; however, it displays two well defined helical regions when in the presence of SDS micelles. The helices correspond to Trp(585)-Ile(593) and Gly(604)-Tyr(611) of the channel. The Trp(585)-Ile(593) helix has distinct hydrophilic and hydrophobic surfaces. The Gly(604)-Tyr(611) helix corresponds to an N-terminal extension of the pore helix. Electrophysiological studies of HERG currents following application of exogenous S5P peptides show that the amphipathic helix in the S5P linker interacts with the pore region of the channel in a voltage-dependent manner.

Expression and biochemical analysis of the entire HIV-2 gp41 ectodomain: determinants of stability map to N- and C-terminal sequences outside the 6-helix bundle core

The folding of HIV gp41 into a 6-helix bundle drives virus-cell membrane fusion. To examine the structural relationship between the 6-helix bundle core domain and other regions of gp41, we expressed in Escherichia coli, the entire ectodomain of HIV-2(ST) gp41 as a soluble, trimeric maltose-binding protein (MBP)/gp41 chimera. Limiting proteolysis indicated that the Cys-591-Cys-597 disulfide-bonded region is outside a core domain comprising two peptides, Thr-529-Trp-589 and Val-604-Ser-666. A biochemical examination of MBP/gp41 chimeras encompassing these core peptides; indicated that the N-terminal polar segment, 521-528, and C-terminal membrane-proximal segment, 658-666, cooperate in stabilizing the ectodomain. A functional interaction between sequences outside the gp41 core may contribute energy to membrane fusion. (C) 2004 Published by Elsevier B.V. on behalf of the Federation of European Biochemical Societies.

The prediction of disease risk in genomic medicine

In recent years, the phrase 'genomic medicine' has increasingly been used to describe a new development in medicine that holds great promise for human health. This new approach to health care uses the knowledge of an individual's genetic make-up to identify those that are at a higher risk of developing certain diseases and to intervene at an earlier stage to prevent these diseases. Identifying genes that are involved in disease aetiology will provide researchers with tools to develop better treatments and cures. A major role within this field is attributed to 'predictive genomic medicine', which proposes screening healthy individuals to identify those who carry alleles that increase their susceptibility to common diseases, such as cancers and heart disease. Physicians could then intervene even before the disease manifests and advise individuals with a higher genetic risk to change their behaviour - for instance, to exercise or to eat a healthier diet - or offer drugs or other medical treatment to reduce their chances of developing these diseases. These promises have fallen on fertile ground among politicians, health-care providers and the general public, particularly in light of the increasing costs of health care in developed societies. Various countries have established databases on the DNA and health information of whole populations as a first step towards genomic medicine. Biomedical research has also identified a large number of genes that could be used to predict someone's risk of developing a certain disorder. But it would be premature to assume that genomic medicine will soon become reality, as many problems remain to be solved. Our knowledge about most disease genes and their roles is far from sufficient to make reliable predictions about a patient’s risk of actually developing a disease. In addition, genomic medicine will create new political, social, ethical and economic challenges that will have to be addressed in the near future.

Prediction of the chemical composition of lamb carcasses from multi-frequency impedance data

Multi-frequency bioimpedance analysis (MFBIA) was used to determine the impedance, reactance and resistance of 103 lamb carcasses (17.1-34.2 kg) immediately after slaughter and evisceration. Carcasses were halved, frozen and one half subsequently homogenized and analysed for water, crude protein and fat content. Three measures of carcass length were obtained. Diagonal length between the electrodes (right side biceps femoris to left side of neck) explained a greater proportion of the variance in water mass than did estimates of spinal length and was selected for use in the index L-2/Z to predict the mass of chemical components in the carcass. Use of impedance (Z) measured at the characteristic frequency (Z(c)) instead of 50 kHz (Z(50)) did not improve the power of the model to predict the mass of water, protein or fat in the carcass. While L-2/Z(50) explained a significant proportion of variation in the masses of body water (r(2) 0.64), protein (r(2) 0.34) and fat (r(2) 0.35), its inclusion in multi-variate indices offered small or no increases in predictive capacity when hot carcass weight (HCW) and a measure of rib fat-depth (GR) were present in the model. Optimized equations were able to account for 65-90 % of the variance observed in the weight of chemical components in the carcass. It is concluded that single frequency impedance data do not provide better prediction of carcass composition than can be obtained from measures of HCW and GR. Indices of intracellular water mass derived from impedance at zero frequency and the characteristic frequency explained a similar proportion of the variance in carcass protein mass as did the index L-2/Z(50).

Prediction of MHC class II-binding peptides using an evolutionary algorithm and artificial neural network

Motivation: Prediction methods for identifying binding peptides could minimize the number of peptides required to be synthesized and assayed, and thereby facilitate the identification of potential T-cell epitopes. We developed a bioinformatic method for the prediction of peptide binding to MHC class II molecules. Results: Experimental binding data and expert knowledge of anchor positions and binding motifs were combined with an evolutionary algorithm (EA) and an artificial neural network (ANN): binding data extraction --> peptide alignment --> ANN training and classification. This method, termed PERUN, was implemented for the prediction of peptides that bind to HLA-DR4(B1*0401). The respective positive predictive values of PERUN predictions of high-, moderate-, low- and zero-affinity binder-a were assessed as 0.8, 0.7, 0.5 and 0.8 by cross-validation, and 1.0, 0.8, 0.3 and 0.7 by experimental binding. This illustrates the synergy between experimentation and computer modeling, and its application to the identification of potential immunotheraaeutic peptides.

Crystallization of a trimeric human T cell leukemia virus type 1 gp21 ectodomain fragment as a chimera with maltose-binding protein

We present a novel protein crystallization strategy, applied to the crystallization of human T cell leukemia virus type 1 (HTLV-1) transmembrane protein gp21 lacking the fusion peptide and the transmembrane domain, as a chimera with the Escherichia coli maltose binding protein (MBP). Crystals could not be obtained with a MBP/gp21 fusion protein in which fusion partners were separated by a flexible linker, but were obtained after connecting the MBP C-terminal alpha-helix to the predicted N-terminal alpha-helical sequence of gp21 via three alanine residues. The gp21 sequences conferred a trimeric structure to the soluble fusion proteins as assessed by sedimentation equilibrium and X-ray diffraction, consistent with the trimeric structures of other retroviral transmembrane proteins. The envelope protein precursor, gp62, is likewise trimeric when expressed in mammalian cells. Our results suggest that MBP may have a general application for the crystallization of proteins containing N-terminal alpha-helical sequences.

Solution structure of amyloid beta-peptide(1-40) in a water-micelle environment. Is the membrane-spanning domain where we think it is?

The three-dimensional solution structure of the 40 residue amyloid beta-peptide, A beta(1-40), has been determined using NMR spectroscopy at pH 5.1, in aqueous sodium dodecyl sulfate (SDS) micelles, In this environment, which simulates to some extent a water-membrane medium, the peptide is unstructured between residues 1 and 14 which are mainly polar and likely solvated by water. However, the rest of the protein adopts an alpha-helical conformation between residues 15 and 36 with a kink or hinge at 25-27. This largely hydrophobic region is likely solvated by SDS. Based on the derived structures, evidence is provided in support of a possible new location for the transmembrane domain of A beta within the amyloid precursor protein (APP). Studies between pH 4.2 and 7.9 reveal a pH-dependent helix-coil conformational switch. At the lower pH values, where the carboxylate residues are protonated, the helix is uncharged, intact, and lipid-soluble. As the pH increases above 6.0, part of the helical region (15-24) becomes less structured, particularly near residues E22 and D23 where deprotonation appears to facilitate unwinding of the helix. This pH-dependent unfolding to a random coil conformation precedes any tendency of this peptide to aggregate to a beta-sheet as the pH increases. The structural biology described herein for A beta(1-40) suggests that (i) the C-terminal two-thirds of the peptide is an alpha-helix in membrane-like environments, (ii) deprotonation of two acidic amino acids in the helix promotes a helix-coil conformational transition that precedes aggregation, (iii) a mobile hinge exists in the helical region of A beta(1-40) and this may be relevant to its membrane-inserting properties and conformational rearrangements, and (iv) the location of the transmembrane domain of amyloid precursor proteins may be different from that accepted in the Literature. These results may provide new insight to the structural properties of amyloid beta-peptides of relevance to Alzheimer's disease.

TFEC is a macrophage-restricted member of the microphthalmia-TFE subfamily of basic helix-loop-helix leucine zipper transcription factors

The murine homologue of the TFEC was cloned as part of an analysis of the expression of the microphthalmia-TFE (MiT) subfamily of transcription factors in macrophages. TFEC, which most likely acts as a transcriptional repressor in heterodimers with other MiT family members, was identified in cells of the mononuclear phagocyte lineage, coexpressed,vith all other known MiT subfamily members (Mitf, TFE3, TFEB), Northern blot analysis of several different cell lineages indicated that the expression of murine TFEC (mTFEC) was restricted to macrophages. A 600-bp fragment of the TATA-less putative proximal promoter of TFEC shares features with many known macrophage-specific promoters and preferentially directs luciferase expression in the RAW264.7 macrophage cell line in transient transfection assays. Five of six putative Ets motifs identified in the TFEC promoter bind the macrophage-restricted transcription factor PU,I under in vitro conditions and in transfected 3T3 fibroblasts; the minimal luciferase activity of the TFEC promoter could be induced by coexpression of PU.1 or the related transcription factor Ets-2. The functional importance of the tissue-restricted expression of TFEC and a possible role in macrophage-specific gene regulation require further investigation, but are likely to be linked to the role of the other MiT family members in this lineage.

Self-efficacy for refusal mediated by outcome expectancies in the prediction of alcohol-dependence amongst young adults

The present study examined the relative importance of outcome expectancies and self-efficacy [1] in the prediction of alcohol dependence [2] and alcohol consumption in a sample of young adult drinkers drawn from a milieu previously reported as supportive of risky drinking. In predicting alcohol dependence, outcome expectancies were found to mediate self-efficacy and the same pattern was found for both males and females. This suggests that male and female drinkers may become more similar as they progress along the drinking continuum from risky drinking to dependent drinking. However, in women, in comparison to men, a greater array of expectancies and self-efficacy scales were found to predict heavy drinking, as measured by quantity and frequency. These results suggest that heavy drinking women are particularly at risk of developing drinking related complications and that preventative education needs to take into account gender differences.

Crystal structure of human T cell leukemia virus type 1 gp21 ectodomain crystallized as a maltose-binding protein chimera reveals structural evolution of retroviral transmembrane proteins

Retroviral entry into cells depends on envelope glycoproteins, whereby receptor binding to the surface-exposed subunit triggers membrane fusion by the transmembrane protein (TM) subunit. We determined the crystal structure at 2.5-Angstrom resolution of the ectodomain of gp21, the TM from human T cell leukemia virus type 1. The gp21 fragment was crystallized as a maltose-binding protein chimera, and the maltose-binding protein domain was used to solve the initial phases by the method of molecular replacement. The structure of gp21 comprises an N-terminal trimeric coiled coil, an adjacent disulfide-bonded loop that stabilizes a chain reversal, and a C-terminal sequence structurally distinct from HIV type 1/simian immunodeficiency virus gp41 that packs against the coil in an extended antiparallel fashion. Comparison of the gp21 structure with the structures of other retroviral TMs contrasts the conserved nature of the coiled coil-forming region and adjacent disulfide-bonded loop with the variable nature of the C-terminal ectodomain segment. The structure points to these features having evolved to enable the dual roles of retroviral TMs: conserved fusion function and an ability to anchor diverse surface-exposed subunit structures to the virion envelope and infected cell surface. The structure of gp21 implies that the N-terminal fusion peptide is in close proximity to the C-terminal transmembrane domain and likely represents a postfusion conformation.