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.     

Analyses of intramolecular disulfide bonds in proteins by polyacrylamide gel electrophoresis following two-step alkylation

A method that makes use of polyacrylamide gel electrophoresis was developed for the analysis of intramolecular disulfide bonds in proteins. Proteins with different numbers of cleaved disulfide bonds are alkylated with iodoacetic acid or iodoacetamide as the first step. The disulfide bonds remaining were reduced by excess dithiothreitol, and the newly generated free sulfhydryl groups were alkylated with the reagent not yet used (iodoacetamide, iodoacetic acid, or vinyl-pyridine) as the second step. This treatment made it possible for lysozyme (Mr, 14,000; 4 disulfides), the N-terminal half-molecule of conalbumin (Mr, 36,000; 6 disulfides), the C-terminal half-molecule of conalbumin (Mr, 40,000; 9 disulfides), and whole conalbumin (Mr, 78,000; 15 disulfides) to be separated by acid-urea polyacrylamide gel electrophoresis into distinct bands depending on the number of disulfide bonds cleaved. The method allowed us to determine the total number of disulfide bonds in native proteins and to assess the cleaved levels of disulfide bonds in partially reduced proteins. Two-step alkylation used in combination with radioautography was especially useful for the analysis of disulfide bonds in proteins synthesized in complex biological systems.     

Denaturation and unfolding of human anaphylatoxin C3a: an unusually low covalent stability of its native disulfide bonds

The complement C3a anaphylatoxin is a major molecular mediator of innate immunity. It is a potent activator of mast cells, basophils and eosinophils and causes smooth muscle contraction. Structurally, C3a is a relatively small protein (77 amino acids) comprising a N-terminal domain connected by 3 native disulfide bonds and a helical C-terminal segment. The structural stability of C3a has been investigated here using three different methods: Disulfide scrambling; Differential CD spectroscopy; and Reductive unfolding. Two uncommon features regarding the stability of C3a and the structure of denatured C3a have been observed in this study. (a) There is an unusual disconnection between the conformational stability of C3a and the covalent stability of its three native disulfide bonds that is not seen with other disulfide proteins. As measured by both methods of disulfide scrambling and differential CD spectroscopy, the native C3a exhibits a global conformational stability that is comparable to numerous proteins with similar size and disulfide content, all with mid-point denaturation of [GdmCl]_1/2 at 3.4-5 M. These proteins include hirudin, tick anticoagulant protein and leech carboxypeptidase inhibitor. However, the native disulfide bonds of C3a is 150-1000 fold less stable than those proteins as evaluated by the method of reductive unfolding. The 3 native disulfide bonds of C3a can be collectively and quantitatively reduced with as low as 1 mM of dithiothreitol within 5 min. The fragility of the native disulfide bonds of C3a has not yet been observed with other native disulfide proteins. (b) Using the method of disulfide scrambling, denatured C3a was shown to consist of diverse isomers adopting varied extent of unfolding. Among them, the most extensively unfolded isomer of denatured C3a is found to assume beads-form disulfide pattern, comprising Cys³⁶-Cys⁴⁹ and two disulfide bonds formed by two pair of consecutive cysteines, Cys²²-Cys²³ and Cys⁵⁶-Cys⁵⁷, a unique disulfide structure of polypeptide that has not been documented previously.     

A method for the controlled cleavage of disulfide bonds in proteins in the absence of denaturants

A simple method was developed for the controlled cleavage of protein disulfide bonds and the simultaneous blockage of the free sulfhydryl groups in the absence of a denaturant. The disulfide bonds of bovine serum albumin were cleaved unsymmetrically at pH 7.0 using 0.1 M sulfite in 0.1 M phosphate buffer and the free sulfhydryl groups formed were sulfonated in an oxidation-reduction cycle using molecular oxygen and 400 microM cupric sulfate as a catalyst. The reaction was affected by cupric ion concentration, sulfite concentration, reaction pH and temperature. The standardized method was successfully used to cleave the disulfide bonds of other proteins pepsin, trypsin, and chymotrypsin. The method is reliable and can be used for achieving progressive cleavage of disulfide bonds in proteins without employing a denaturant.     

Predicting disulfide connectivity from protein sequence using multiple sequence feature vectors and secondary structure

MOTIVATION: Disulfide bonds are primary covalent crosslinks between two cysteine residues in proteins that play critical roles in stabilizing the protein structures and are commonly found in extracy-toplasmatic or secreted proteins. In protein folding prediction, the localization of disulfide bonds can greatly reduce the search in conformational space. Therefore, there is a great need to develop computational methods capable of accurately predicting disulfide connectivity patterns in proteins that could have potentially important applications. RESULTS: We have developed a novel method to predict disulfide connectivity patterns from protein primary sequence, using a support vector regression (SVR) approach based on multiple sequence feature vectors and predicted secondary structure by the PSIPRED program. The results indicate that our method could achieve a prediction accuracy of 74.4% and 77.9%, respectively, when averaged on proteins with two to five disulfide bridges using 4-fold cross-validation, measured on the protein and cysteine pair on a well-defined non-homologous dataset. We assessed the effects of different sequence encoding schemes on the prediction performance of disulfide connectivity. It has been shown that the sequence encoding scheme based on multiple sequence feature vectors coupled with predicted secondary structure can significantly improve the prediction accuracy, thus enabling our method to outperform most of other currently available predictors. Our work provides a complementary approach to the current algorithms that should be useful in computationally assigning disulfide connectivity patterns and helps in the annotation of protein sequences generated by large-scale whole-genome projects. AVAILABILITY: The prediction web server and Supplementary Material are accessible at http://foo.maths.uq.edu.au/~huber/disulfide     

Assignment of disulfide bonds in proteins by partial acid hydrolysis and mass spectrometry

A new method is described for locating disulfide bonds in proteins which cannot be cleaved between half-cystinyl residues by enzymic methods, as is often the case for tightly coiled proteins, or for proteins in which half-cystinyl residues are not separated by residues required for enzymic cleavage. Partial acid hydrolysis of a model protein, hen egg-white lysozyme, produces a mixture of disulfide-containing peptides from which the disulfide connections may be deduced. The usefulness of a combination of HPLC, fast atom bombardment mass spectrometry, and computer-assisted analysis to identify disulfide-containing peptides present in the partial acid hydrolysate of the model protein is demonstrated. Chromatographic fractions of the hydrolysate were analyzed by mass spectrometry before and after chemical reduction of the disulfide bonds to determine the molecular weights of disulfide-containing peptides. Computer-assisted analysis was then used to relate the molecular weights of these peptides to specific segments of the protein from which the disulfide connectivities could be determined. Partial acid hydrolysis of proteins, which is attractive because it proceeds relatively independent of the amino acid sequence and structure, and because disulfide interchange is unlikely to occur in dilute acid, has become practical because disulfide-containing peptides present in complex mixtures can be identified rapidly and definitively by this method.     

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

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

The genomics of disulfide bonding and protein stabilization in thermophiles

Thermophilic organisms flourish in varied high-temperature environmental niches that are deadly to other organisms. Recently, genomic evidence has implicated a critical role for disulfide bonds in the structural stabilization of intracellular proteins from certain of these organisms, contrary to the conventional view that structural disulfide bonds are exclusively extracellular. Here both computational and structural data are presented to explore the occurrence of disulfide bonds as a protein-stabilization method across many thermophilic prokaryotes. Based on computational studies, disulfide-bond richness is found to be widespread, with thermophiles containing the highest levels. Interestingly, only a distinct subset of thermophiles exhibit this property. A computational search for proteins matching this target phylogenetic profile singles out a specific protein, known as protein disulfide oxidoreductase, as a potential key player in thermophilic intracellular disulfide-bond formation. Finally, biochemical support in the form of a new crystal structure of a thermophilic protein with three disulfide bonds is presented together with a survey of known structures from the literature. Together, the results provide insight into biochemical specialization and the diversity of methods employed by organisms to stabilize their proteins in exotic environments. The findings also motivate continued efforts to sequence genomes from divergent organisms.     

Cysteine separations profiles on protein sequences infer disulfide connectivity

MOTIVATION: Disulfide bonds play an important role in protein folding. A precise prediction of disulfide connectivity can strongly reduce the conformational search space and increase the accuracy in protein structure prediction. Conventional disulfide connectivity predictions use sequence information, and prediction accuracy is limited. Here, by using an alternative scheme with global information for disulfide connectivity prediction, higher performance is obtained with respect to other approaches. RESULT: Cysteine separation profiles have been used to predict the disulfide connectivity of proteins. The separations among oxidized cysteine residues on a protein sequence have been encoded into vectors named cysteine separation profiles (CSPs). Through comparisons of their CSPs, the disulfide connectivity of a test protein is inferred from a non-redundant template set. For non-redundant proteins in SwissProt 39 (SP39) sharing less than 30% sequence identity, the prediction accuracy of a fourfold cross-validation is 49%. The prediction accuracy of disulfide connectivity for proteins in SwissProt 43 (SP43) is even higher (53%). The relationship between the similarity of CSPs and the prediction accuracy is also discussed. The method proposed in this work is relatively simple and can generate higher accuracies compared to conventional methods. It may be also combined with other algorithms for further improvements in protein structure prediction. AVAILABILITY: The program and datasets are available from the authors upon request. CONTACT: cykao@csie.ntu.edu.tw.     

Online mass spectrometric analysis of proteins/peptides following electrolytic cleavage of disulfide bonds

The disulfide bond bridge is an important post-translational modification for proteins. This study presents a structural analysis of biologically active peptides and proteins containing disulfide bonds using electrochemistry (EC) online combined with desorption electrospray ionization mass spectrometry (DESI-MS), in which the sample undergoes electrolytic disulfide cleavage in an electrochemical flow cell followed by MS detection. Using this EC/DESI-MS method, the disulfide-containing peptides can be quickly identified from enzymatic digestion mixtures, simply based on the abrupt decrease in their relative ion abundances after electrolysis. Peptide mass mapping and tandem MS analysis of the ions of the resulting free peptide chains can possibly establish the disulfide linkage pattern and sequence the precursor peptides. In this regard, the method provides much more chemical information than previous analogous electrochemical analyses. In addition, derivatization of thiols by selective selenamide reagents is useful for easy recognition of reduced peptide ions and the number of their free thiols. Furthermore, electrolytic reduction of proteins (e.g., α-lactalbumin) leads to increased charges on the detected protein ions, revealing the role of disulfide bonds on maintaining protein conformation. This electrochemical mass spectrometric method is fast (completed in few minutes) and does not need chemical reductants, potentially having valuable applications in proteomics research.     

Sensitive quantitative analysis of disulfide bonds in polypeptides and proteins

A sensitive quantitative method has been developed to determine the number of disulfide bonds in peptides and proteins. The disulfide bonds of several peptides and proteins were cleaved quantitatively by excess sodium sulfite at pH 9.5 and room temperature. Guanidine thiocyanate (2 M) was added to the protein solutions in order to denature them and thereby make the disulfide bonds accessible. The reaction with sulfite leads to a thiosulfonate and a free sulfhydryl group; the concentration of the latter was determined by reaction with disodium 2-nitro-5-thiosulfobenzoate (NTSB) in the presence of excess sodium sulfite. The synthesis, purification, and characterization of NTSB are described. The assay is rapid, requiring 3-5 min for oligopeptides and 20 min for proteins, and is as sensitive and quantitative as the sulfhydryl group assay employing 5,5'-dithiobis(2-nitrobenzoic acid) (Ellman's reagent). It can be used for the analysis of as little as 10(8) mol of disulfide bonds, with an error of +/- 3%.     

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.     

Disulfide crosslinks to probe the structure and flexibility of a designed four-helix bundle protein.

The introduction of disulfide crosslinks is a generally useful method by which to identify regions of a protein that are close together in space. Here we describe the use of disulfide crosslinks to investigate the structure and flexibility of a family of designed 4-helix bundle proteins. The results of these analyses lend support to our working model of the proteins' structure and suggest that the proteins have limited main-chain flexibility.     

Predicting disulfide bond connectivity in proteins by correlated mutations analysis

MOTIVATION: Prediction of disulfide bond connectivity facilitates structural and functional annotation of proteins. Previous studies suggest that cysteines of a disulfide bond mutate in a correlated manner. RESULTS: We developed a method that analyzes correlated mutation patterns in multiple sequence alignments in order to predict disulfide bond connectivity. Proteins with known experimental structures and varying numbers of disulfide bonds, and that spanned various evolutionary distances, were aligned. We observed frequent variation of disulfide bond connectivity within members of the same protein families, and it was also observed that in 99% of the cases, cysteine pairs forming non-conserved disulfide bonds mutated in concert. Our data support the notion that substitution of a cysteine in a disulfide bond prompts the substitution of its cysteine partner and that oxidized cysteines appear in pairs. The method we developed predicts disulfide bond connectivity patterns with accuracies of 73, 69 and 61% for proteins with two, three and four disulfide bonds, respectively.     

Multi-template approach to modeling engineered disulfide bonds

The key issue for disulfide bond engineering is to select the most appropriate location in the protein. By surveying the structure of experimentally engineered disulfide bonds, we found about half of them that have geometry incompatible with any native disulfide bond geometry. To improve the current prediction methods that tend to apply either ideal geometrical or energetical criteria to single three-dimensional structures, we have combined a novel computational protocol with the usage of multiple protein structures to take into account protein backbone flexibility. The multiple structures can be selected from either independently determined crystal structures for identical proteins, models of nuclear magnetic resonance experiments, or crystal structures of homology-related proteins. We have validated our approach by comparing the predictions with known disulfide bonds. The accuracy of prediction for native disulfide bonds reaches 99.6%. In a more stringent test on the reported engineered disulfide bonds, we have obtained a success rate of 93%. Our protocol also determines the oxido-reduction state of a predicted disulfide bond and the corresponding mutational cost. From the energy ranking, the user can easily choose top predicted sites for mutagenesis experiments. Our method provides information about local stability of the engineered disulfide bond surroundings.     

A reporter platform for the monitoring of in vivo conformational changes in AcrB

AcrB is an inner membrane protein in Escherichia coli that is a component of a triplex AcrA-AcrB-TolC (AcrAB-TolC) multidrug efflux pump. The AcrAB-TolC complex and its homologues are highly conserved among Gram-negative bacteria and are major players in conferring multidrug resistance (MDR) in many pathogens. In this study we developed a disulfide trapping method that may reveal AcrB conformational changes under the native condition in the cell membrane. Specifically, we created seven disulfide bond pairs in the periplasmic domain of AcrB, which can be used as probes to determine local conformational changes. We have developed a rigorous protocol to measure the extent of disulfide bond formation in membrane proteins under the native condition. The rigorousness of the method was verified through examining the extent of disulfide bond formation in Cys pairs separated by different distances. The blocking-reducing-labeling scheme combined with fluorescence labeling made the current method convenient, reliable, and quantitative.     

Conformational stability of secretory leucocyte protease inhibitor: a protein with no hydrophobic core and very little secondary structure

Conformational stability of proteins (including disulfide containing proteins) has been routinely characterized by spectroscopic techniques. Proteins which lack adequate signal of circular dichroism may require unconventional technique. Secretory Leucocyte Protease Inhibitor (SLPI) is a 107 amino acids protein with a high density of disulfide pairing (eight). The native SLPI has no hydrophobic core and contains very little hydrogen bonded secondary structure [Gruetter, M., Fendrich, G., Huber, R., and Bode, W. (1988) The 2.5 A X-ray crystal structure of the acid stable proteinase inhibitor from human mucous secretions analyzed in its complex with bovine alpha-chymotrypsin. The EMBO J. 7, 345-352.]. In this study, conformational stability of SLPI has been investigated by the method of disulfide scrambling, which permits quantification of the native and denatured (scrambled) proteins by HPLC. Due to high heterogeneity of denatured SLPI, the native and scrambled SLPI are extensively overlapped on HPLC. This impediment was further overcome by the development of a novel method which distinguishes the native and scrambled isomers of SLPI by exploiting the relative stability of their disulfide bonds. The study reveals mid-point denaturation of SLPI at 1.36 M of GdmSCN, 4.0 M of GdmCl and >8 M urea. Based on the GdmCl denaturation curve, the unfolding free energy (DeltaG(H20)) of SLPI was estimated to be 4.56 kcal/mol. The results of our studies suggest an alternative strategy for analyzing conformational stability of disulfide proteins that are not suitable to the conventional spectroscopic techniques.