Interresidue interactions in protein classes

The free energy difference between folded and unfolded state is about the same for most proteins and it is not more than the energy of a few noncovalent interactions. In addition to the numerous noncovalent interactions, some proteins contain one or more disulfide bonds, which, as covalent crosslinks, significantly stabilize their tertiary structure. Correlation between the presence of disulfide bond(s), and the number noncovalent interresidue interactions of various kinds is analyzed here. The number of interactions per residue is almost the same for all protein. Also the number of long-range interactions per residue is the same in all proteins. Proteins with S(SINGLE BOND)S bond(s) (extracellular proteins) have more medium-range and fewer short-range interactions than those without S(SINGLE BOND)S bonds. However, the difference is independent of the number of these covalent crosslinks. We concluded that the different distributions of the various kinds of noncovalent interaction reflect the needs of proteins in the different environments, the extracellular and the intracellular ones, rather than the presence of the disulfide bond(s). We also pointed out that the observed differences in the distributions of short- and medium-range interactions are in good agreement with different secondary structure compositions of extracellular and intracellular proteins. Proteins 27:360–366, 1997. © 1997 Wiley-Liss, Inc.     

支持向量机方法预测蛋白质结构中的二硫键

在蛋白质结构预测的研究中,一个重要的问题就是正确预测二硫键的连接,二硫键的准确预测可以减少蛋白质构像的搜索空间,有利于蛋白质3D结构的预测,本文将预测二硫键的连接问题转化成对连接模式的分类问题,并成功地将支持向量机方法引入到预测工作中。通过对半胱氨酸局域序列连接模式的分类预测,可以由蛋白质的一级结构序列预测该蛋白质的二硫键的连接。结果表明蛋白质的二硫键的连接与半胱氨酸局域序列连接模式有重要联系,应用支持向量机方法对蛋白质结构的二硫键预测取得了良好的结果。     

A sequence-profile-based HMM for predicting and discriminating beta barrel membrane proteins

MOTIVATION: Membrane proteins are an abundant and functionally relevant subset of proteins that putatively include from about 15 up to 30% of the proteome of organisms fully sequenced. These estimates are mainly computed on the basis of sequence comparison and membrane protein prediction. It is therefore urgent to develop methods capable of selecting membrane proteins especially in the case of outer membrane proteins, barely taken into consideration when proteome wide analysis is performed. This will also help protein annotation when no homologous sequence is found in the database. Outer membrane proteins solved so far at atomic resolution interact with the external membrane of bacteria with a characteristic beta barrel structure comprising different even numbers of beta strands (beta barrel membrane proteins). In this they differ from the membrane proteins of the cytoplasmic membrane endowed with alpha helix bundles (all alpha membrane proteins) and need specialised predictors. RESULTS: We develop a HMM model, which can predict the topology of beta barrel membrane proteins using, as input, evolutionary information. The model is cyclic with 6 types of states: two for the beta strand transmembrane core, one for the beta strand cap on either side of the membrane, one for the inner loop, one for the outer loop and one for the globular domain state in the middle of each loop. The development of a specific input for HMM based on multiple sequence alignment is novel. The accuracy per residue of the model is 83% when a jack knife procedure is adopted. With a model optimisation method using a dynamic programming algorithm seven topological models out of the twelve proteins included in the testing set are also correctly predicted. When used as a discriminator, the model is rather selective. At a fixed probability value, it retains 84% of a non-redundant set comprising 145 sequences of well-annotated outer membrane proteins. Concomitantly, it correctly rejects 90% of a set of globular proteins including about 1200 chains with low sequence identity (<30%) and 90% of a set of all alpha membrane proteins, including 188 chains.     

Disulfide proteome yields a detailed understanding of redox regulations: a model study of thioredoxin-linked reactions in seed germination

Accumulating evidence suggests that redox regulations play important roles in a broad spectrum of biological processes. Recently, Yano et al. developed a disulfide proteome technique that comprehensively visualizes redox change in proteins. In this paper, using the disulfide proteome, we examined rice bran and identified fragments of embryo-specific protein and dienelactone hydrolase as putative targets of thioredoxin. Also, monitoring of the endogenous and recombinant effects of thioredoxin on rice bran proteins and supporting in vivo observations propose a mechanism of redox regulation in seed germination, in which thioredoxin activates cysteine protease with a concurrent unfolding of its substrate, the embryo-specific protein. Our findings suggest that thioredoxin controls the lifetime of specific proteins effectively by regulating the redox reactions coordinately. The model study demonstrates that the disulfide proteome technique is useful not only for identifying targets of thioredoxin, but also for clarify the detailed mechanism of redox regulation.     

Levels of sulfhydryls and disulfides in proteins from Neurospora crassa conidia and mycelia

Proteins extracted with 6 M guanidine at 90 degrees C from conidia (asexual spores) of Neurospora crassa contained ca. 25% more total protein thiol and a fivefold-higher content of disulfide bonds than proteins extracted from mycelia, as determined by labeling with iodo[14C]acetic acid. The total thiol content was 88 mumol/g of protein in conidia and 70 mumol/g of protein in mycelia. The level of protein disulfide was 18.5 mumol/g of protein in conidia and 3.5 mumol/g of protein in mycelia, by the iodo[14C]acetic acid labeling method. Confirmatory results were obtained with 5'5-dithio-bis-2-nitrobenzoic acid titration of protein thiol groups in 1% sodium dodecyl sulfate as well as by amino acid analysis of cysteic acid derivatives. Buffer-extracted proteins from conidia, but not mycelia, were found to contain enriched levels of protein thiols and disulfides per gram of protein as compared with guanidine hydrochloride extracts. It was demonstrated that the high disulfide content of crude conidial extracts was not due to measurable levels of mixed disulfides formed between protein sulfhydryl groups and cysteine. During germination of the conidia, the high disulfide levels of the conidial proteins remained constant. These data suggest that, unlike the disulfides of glutathione, the bulk of conidial protein disulfides were not reduced, excreted, or extensively degraded during germination.     

Labile disulfide bonds are common at the leucocyte cell surface

Redox conditions change in events such as immune and platelet activation, and during viral infection, but the biochemical consequences are not well characterized. There is evidence that some disulfide bonds in membrane proteins are labile while others that are probably structurally important are not exposed at the protein surface. We have developed a proteomic/mass spectrometry method to screen for and identify non-structural, redox-labile disulfide bonds in leucocyte cell-surface proteins. These labile disulfide bonds are common, with several classes of proteins being identified and around 30 membrane proteins regularly identified under different reducing conditions including using enzymes such as thioredoxin. The proteins identified include integrins, receptors, transporters and cell-cell recognition proteins. In many cases, at least one cysteine residue was identified by mass spectrometry as being modified by the reduction process. In some cases, functional changes are predicted (e.g. in integrins and cytokine receptors) but the scale of molecular changes in membrane proteins observed suggests that widespread effects are likely on many different types of proteins including enzymes, adhesion proteins and transporters. The results imply that membrane protein activity is being modulated by a 'redox regulator' mechanism.     

Prediction of disulfide connectivity in proteins.

MOTIVATION: A major problem in protein structure prediction is the correct location of disulfide bridges in cysteine-rich proteins. In protein-folding prediction, the location of disulfide bridges can strongly reduce the search in the conformational space. Therefore the correct prediction of the disulfide connectivity starting from the protein residue sequence may also help in predicting its 3D structure. RESULTS: In this paper we equate the problem of predicting the disulfide connectivity in proteins to a problem of finding the graph matching with the maximum weight. The graph vertices are the residues of cysteine-forming disulfide bridges, and the weight edges are contact potentials. In order to solve this problem we develop and test different residue contact potentials. The best performing one, based on the Edmonds-Gabow algorithm and Monte-Carlo simulated annealing reaches an accuracy significantly higher than that obtained with a general mean force contact potential. Significantly, in the case of proteins with four disulfide bonds in the structure, the accuracy is 17 times higher than that of a random predictor. The method presented here can be used to locate putative disulfide bridges in protein-folding. AVAILABILITY: The program is available upon request from the authors. CONTACT: Casadio@alma.unibo.it; Piero@biocomp.unibo.it.     

Direct analysis of the disulfide content of proteins: methods for monitoring the stability and refolding process of cystine-containing proteins

So far, the cystine (disulfide) content of proteins has been routinely determined by indirect methods, i.e., cystines are first converted to more stable half-cystines prior to analysis. We present a study which demonstrates that the cystine content can be directly and accurately analyzed at low picomole to femtomole levels. The method involves (a) the employment of a vacuum during sample hydrolysis which permits quantitative recovery of cystine, and (b) the dabsyl chloride precolumn derivatization method which yields stable DABS-cystine that is subsequently analyzed by HPLC. Direct analysis of the cystine residue is important in numerous areas of protein research. It allows, by a single analysis, simultaneous determination of the sulfhydryl (cysteine, as S-carboxymethyl derivative) and disulfide (cystine) contents. The technique can be used to follow the refolding process of fully reduced, cystine-containing proteins. Direct analysis of the cystine content also serves to monitor the extent of inactivation of cystine-containing proteins caused by alkaline pH and heat treatments. In this mode of protein inactivation, cystines are selectively destroyed and converted to lanthionine and lysinoalanine. Both the decrease of cystine and the recoveries of lanthionine and lysinoalanine can be simultaneously evaluated by the proposed method. Examples of these applications are presented here.     

High throughput screening of disulfide‐containing proteins in a complex mixture

The formation of disulfide bonds between cysteine residues is crucial for the stabilization of native protein structures and, thus, determination of disulfide linkages is an important facet of protein structural characterization. Nonetheless, the identification of disulfide bond linkages remains a significant analytical challenge, particularly in large proteins with complex disulfide patterns. Herein, we have developed a new LC/MS strategy for rapid screening of disulfides in an intact protein mixture after a straightforward reduction step with tris(2‐carboxyethyl)phosphine. LC/MS analysis of reduced and nonreduced protein mixtures quickly revealed disulfide‐containing proteins owing to a 2 Da mass increase per disulfide reduction and, subsequently, the total number of disulfide bonds in the intact proteins could be determined. We have demonstrated the effectiveness of this method in a protein mixture composed of both disulfide‐containing and disulfide‐free proteins. Our method is simple (no need for proteolytic digestion, alkylation, or the removal of reducing agents prior to MS analysis), high throughput (fast on‐line LC/MS analysis), and reliable (no S–S scrambling), underscoring its potential as a rapid disulfide screening method for proteomics applications.     

The challenge of the proteome dynamic range and its implications for in‐depth proteomics

The dynamic range of the cellular proteome approaches seven orders of magnitude—from one copy per cell to ten million copies per cell. Since a proteome's abundance distribution represents a nearly symmetric bell‐shape curve on the logarithmic copy number scale, detection of half of the expressed cellular proteome, i.e. approximately 5000 proteins, should be a relatively straightforward task with modern mass spectrometric instrumentation that exhibits four orders of magnitude of the dynamic range, while deeper proteome analysis should be progressively more difficult. Indeed, metaanalysis of 15 recent papers that claim detection of >5000 protein groups reveals that the half‐proteome analyses currently requires ≈5 h of chromatographic separation, while deeper analyses yield on average ≤20 new proteins per hour of chromatographic gradient. Therefore, a typical proteomics experiment consists of a “high‐content” part, with the detection rate of approximately 1000 proteins/h, and a “low‐content” tail with much lower rate of discovery and respectively, lower cost efficiency. This result calls for disruptive innovation in deep proteomics analysis.     

Structural disulfide bonds in the Bacillus thuringiensis subsp. israelensis protein crystal.

We examined disulfide bonds in mosquito larvicidal crystals produced by Bacillus thuringiensis subsp. israelensis. Intact crystals contained 2.01 X 10(-8) mol of free sulfhydryls and 3.24 X 10(-8) mol of disulfides per mg of protein. Reduced samples of alkali-solubilized crystals resolved into several proteins, the most prominent having apparent molecular sizes of 28, 70, 135, and 140 kilodaltons (kDa). Nonreduced samples contained two new proteins of 52 and 26 kDa. When reduced, both the 52- and 26-kDa proteins were converted to 28-kDa proteins. Furthermore, both bands reacted with antiserum prepared against reduced 28-kDa protein. Approximately 50% of the crystal proteins could be solubilized without disulfide cleavage. These proteins were 70 kDa or smaller. Solubilization of the 135- and 140-kDa proteins required disulfide cleavage. Incubation of crystals at pH 12.0 for 2 h cleaved 40% of the disulfide bonds and solubilized 83% of the crystal protein. Alkali-stable disulfides were present in both the soluble and insoluble portions. The insoluble pellet contained 12 to 14 disulfides per 100 kDa of protein and was devoid of sulfhydryl groups. Alkali-solubilized proteins contained both intrachain and interchain disulfide bonds. Despite their structural significance, it is unlikely that disulfide bonds are involved in the formation or release of the larvicidal toxin.     

Pressure effects on the structure and function of human thioredoxin

Thioredoxin is one of the major proteins that catalyze disulfide reduction and defines the thioredoxin superfamily bearing the CXXC structural motif. Human thioredoxin contains only 1 Trp residue proximal to the active site (WCGPC). We are interested in thioredoxin structure-function relationships, in particular, active site hydration and flexibility. Hence, in this study, we used hydrostatic pressure as a perturbation and monitored the conformational changes around the active site of thioredoxin by analyzing Trp fluorescence. The structure of thioredoxin was drastically altered by increasing pressure and did not completely refold after pressure release. The conformation in the active site vicinity was modified at low pressure (less than 100 MPa) and the Trp residue was completely exposed to aqueous medium at pressures above 350 MPa. Upon pressure release, thioredoxin showed no activity, although it folded 80% of the alpha-helical content relative to the native state. According to these results, pressure denaturation induces critical damage for the activity of thioredoxin, indicating extreme fragility of the active site with respect to pressure. This result is in contrast to the pressure effect on protein disulfide isomerase (PDI) which is organized by four thioredoxin-like domains including two WCGHC motifs.     

PDI improves secretion of redox-inactive beta-glucosidase

Although manipulation of the endoplasmic reticulum (ER) folding environment in the yeast Saccharomyces cerevisiae has been shown to increase the secretory productivity of recombinant proteins, the cellular interactions and processes of native enzymes and chaperones such as protein disulfide isomerase (PDI) are still unclear. Previously, we reported that overexpression of the ER chaperone PDI enabled up to a 3-fold increase in secretion levels of the Pyrococcus furiosus beta-glucosidase in the yeast S. cerevisiae. This result was surprising since beta-glucosidase contains only one cysteine per monomer and no disulfide bonds. Two possible mechanisms were proposed: PDI either forms a transient disulfide bond with the lone cysteine residue of the nascent beta-glucosidase during the folding and assembly process or acts as a chaperone to aid in proper folding. To discern between the two mechanisms, the single cysteine residue was mutated to serine, and the secretion of the two protein variants was determined. The serine mutant still showed increased secretion in vivo when PDI levels were elevated. When the folding bottleneck is removed by increasing expression temperatures to 37 degrees C rather than 30 degrees C, PDI no longer has an improvement on secretion. These results suggest that, unexpectedly, PDI acts in a chaperone-like capacity or possibly cooperates with the cell's folding or degradation mechanisms regardless of whether the protein is redox-active.     

An assessment of the impacts of molecular oxygen on the evolution of proteomes

Oxygen is not only one of life's essential elements but also a source of protein damage, mutagenesis, and ageing. Many proteome adaptations have been proposed to tackle such stresses and we assessed them using comparative genomics in a phylogenetic context. First, we find that aerobiosis is a trait with important phylogenetic inertia but that oxygen content in proteins is not. Instead, oxygen content is close to the expected values given the nucleotide composition. Accordingly, we find no evidence of oxygen being a scarce resource for protein synthesis even among anaerobes. Second, we searched for counterselection of amino acids more prone to oxidation among aerobes. Only cysteine follows the expected trend, whereas tryptophan follows the inverse one. When analyzing composition in the context of protein structures and residue accessibility, we find that all oxidable residues are avoided at the surface of proteins. Yet, there is no difference between aerobes and anaerobes in this respect, and the effect might be explained by the hydrophobicity of these residues. Third, we revisited the hypothesis that atmospheric enrichment in molecular oxygen led to the development of the communication capabilities of eukaryotes. With a larger data set and adequate controls, we confirm the trend of longer oxygen-rich outer domains in transmembrane proteins of eukaryotes. Yet, we find no significant association between oxygen concentration in the environment and this trait within prokaryotes, suggesting that this difference is clade specific and independent of oxygen availability. We find that genes involved in cellular responses to oxygen are much more frequent among aerobes, and we suggest that they erase most expected differences in terms of proteome composition between organisms facing high and low oxygen concentrations.     

Heterogeneity of the covalent structure of the blue copper protein umecyanin from horseradish roots

The covalent structure of umecyanin has been determined by a combination of classical Edman degradation sequence analysis and plasma desorption, laser desorption, and electrospray ionization mass spectrometry. The preparation appeared to contain two isoforms having either a valine (75%) or an isoleucine (25%) residue at position 48. The polypeptide chain of 115 amino acids is strongly heterogeneous at its C-terminal end as a result of proteolytic cleavages at several places within the last 10 residues. The major fraction of the umecyanin preparation is only 106 residues long. The C-terminal tail 107–115 contains mainly alanine and glycine residues and a single hydroxyproline residue. In the native protein there is a disulfide bridge between Cys 91 and Cys 57, but in the apoprotein there is a disulfide shift that involves Cys 91 and one of the four copper binding residues (Cys 85). The three other ligand binding residues are His 44, His 90, and Gin 95. This tetrad of amino acids is the same as occurs in other type 1 copper proteins from plants such as cucumber peeling cupredoxin and lacquer tree stellacyanin. The umecyanin isoforms are glycoproteins with a glycan core having the same carbohydrate composition as that of horseradish peroxidase, a fact that is convincingly supported thanks to the high accuracy of the electrospray mass spectrometric technique. We suggest that the glycan may play a role in the association of the protein to the cellular membrane, but the precise functional role of umecyanin remains to be determined. We also discuss the evolutionary position of umecyanin in relation to the type 1 copper proteins in general.     

Helix-helix packing and interfacial pairwise interactions of residues in membrane proteins

Helix-helix packing plays a critical role in maintaining the tertiary structures of helical membrane proteins. By examining the overall distribution of voids and pockets in the transmembrane (TM) regions of helical membrane proteins, we found that bacteriorhodopsin and halorhodopsin are the most tightly packed, whereas mechanosensitive channel is the least tightly packed. Large residues F, W, and H have the highest propensity to be in a TM void or a pocket, whereas small residues such as S, G, A, and T are least likely to be found in a void or a pocket. The coordination number for non-bonded interactions for each of the residue types is found to correlate with the size of the residue. To assess specific interhelical interactions between residues, we have developed a new computational method to characterize nearest neighboring atoms that are in physical contact. Using an atom-based probabilistic model, we estimate the membrane helical interfacial pairwise (MHIP) propensity. We found that there are many residue pairs that have high propensity for interhelical interactions, but disulfide bonds are rarely found in the TM regions. The high propensity pairs include residue pairs between an aromatic residue and a basic residue (W-R, W-H, and Y-K). In addition, many residue pairs have high propensity to form interhelical polar-polar atomic contacts, for example, residue pairs between two ionizable residues, between one ionizable residue and one N or Q. Soluble proteins do not share this pattern of diverse polar-polar interhelical interaction. Exploratory analysis by clustering of the MHIP values suggests that residues similar in side-chain branchness, cyclic structures, and size tend to have correlated behavior in participating interhelical interactions. A chi-square test rejects the null hypothesis that membrane protein and soluble protein have the same distribution of interhelical pairwise propensity. This observation may help us to understand the folding mechanism of membrane proteins.     

Thioredoxins and glutaredoxins as facilitators of protein folding

Thiol-disulfide oxidoreductase systems of bacterial cytoplasm and eukaryotic cytosol favor reducing conditions and protein thiol groups, while bacterial periplasm and eukaryotic endoplasmatic reticulum provide oxidizing conditions and a machinery for disulfide bond formation in the secretory pathway. Oxidoreductases of the thioredoxin fold superfamily catalyze steps in oxidative protein folding via protein-protein interactions and covalent catalysis to act as chaperones and isomerases of disulfides to generate a native fold. The active site dithiol/disulfide of thioredoxin fold proteins is CXXC where variations of the residues inside the disulfide ring are known to increase the redox potential like in protein disulfide isomerases. In the catalytic mechanism thioredoxin fold proteins bind to target proteins through conserved backbone-backbone hydrogen bonds and induce conformational changes of the target disulfide followed by nucleophilic attack by the N-terminally located low pK(a) Cys residue. This generates a mixed disulfide covalent bond which subsequently is resolved by attack from the C-terminally located Cys residue. This review will focus on two members of the thioredoxin superfamily of proteins known to be crucial for maintaining a reduced intracellular redox state, thioredoxin and glutaredoxin, and their potential functions as facilitators and regulators of protein folding and chaperone activity.