Mass spectrometric based mapping of the disulfide bonding patterns of integrin alpha chains

Integrins are one of the major mediators of cellular adherence. Structurally the component alpha and beta chains are characterized by extensive intrachain disulfide bonding. The assignment of these bonds is currently based on homology with the chains of the integrin alphaIIbbeta3. However, recent crystallographic analysis of the soluble alphaVbeta3 construct indicates that the alphaV chain displays bonding patterns different from those predicted for alphaIIb. In an effort to define the disulfide bonding patterns in integrins, we have used mass spectrometric based approaches to map the human alpha3, alpha5, alphaV, and alphaIIb. The results indicate that there are differences in the disulfide patterns of the alpha chains. These do not correlate with the integrin capacity to bind ligands as all integrins used in the present study displayed functional activity. The differences were observed in the bonding patterns linking the heavy (H) and light (L) components of the of the alpha chains. It was also possible to assign the location in alpha5 of an additional disulfide bond involving a pair of cysteines not present in alphaV or alphaIIb. This second bond between the H and L chains of alpha5 has not been previously described. These results indicate that not all integrin species display the same disulfide bonding patterns. They also highlight the need for caution in the use of assignments based on sequence homology.     

Disulfide bonds in homo- and heterodimers of EF-hand subdomains of calbindin D9k: stability, calcium binding, and NMR studies.

The effect of decreased protein flexibility on the stability and calcium binding properties of calbindin D9k has been addressed in studies of a disulfide bridged calbindin D9k mutant, denoted (L39C + P43M + I73C), with substitutions Leu 39-->Cys, Ile 73-->Cys, and Pro 43-->Met. Backbone 1H NMR assignments show that the disulfide bond, which forms spontaneously under air oxidation, is well accommodated. The disulfide is inserted on the opposite end of the protein molecule with respect to the calcium sites, to avoid direct interference with these sites, as confirmed by 113Cd NMR. The effect of the disulfide bond on calcium binding was assessed by titrations in the presence of a chromophoric chelator. A small but significant effect on the cooperativity was found, as well as a very modest reduction in calcium affinity. The disulfide bond increases Tm, the transition midpoint of thermal denaturation, of calcium free calbindin D9k from 85 to 95 degrees C and Cm, the urea concentration of half denaturation, from 5.3 to 8.0 M. Calbindins with one covalent bond linking the two EF-hand subdomains are equally stable regardless if the covalent link is the 43-44 peptide bond or the disulfide bond. Kinetic remixing experiments show that separated CNBr fragments of (L39C + P43M + I73C), each comprising one EF-hand, form disulfide linked homodimers. Each homodimer binds two calcium ions with positive co-operativity, and an average affinity of 10(6) M-1. Disulfide linkage dramatically increases the stability of each homodimer. For the homodimer of the C-terminal fragment Tm increases from 59 +/- 2 without covalent linkage to 91 +/- 2 degrees C with disulfide, and Cm from approximately 1.5 to 7.5 M. The overall topology of this homodimer is derived from 1H NMR assignments and a few key NOEs.     

Determination of the disulfide bond pairings in human tissue factor pathway inhibitor purified fromEscherichia coli

The disulfide bond assignments of human alanyl tissue factor pathway inhibitor purified fromEscherichia coli have been determined. This inhibitor of the extrinsic blood coagulation pathway possesses three Kunitz-type inhibitor domains, each containing three disulfide bonds. The disulfide bond pairings in domains 1 and 3 were determined by amino acid sequencing and mass spectrometry of peptides derived from a thermolysin digest. However, thermolysin digestion did not cleave any peptide bonds within domain 2. The disulfide bond pairings in domain 2 were determined by isolating it from the thermolysin treatment and subsequently cleaving it with pepsin and trypsin into peptides which yielded the three disulfide bond pairings in this domain. These results demonstrate that the disulfide pairings in each of the three domains of human tissue factor pathway inhibitor purified fromEscherichia coli are homologous to each other and also to those in bovine pancreatic trypsin inhibitor.     

Determination of the disulfide bond pairings in human tissue factor pathway inhibitor purified fromEscherichia coli

The disulfide bond assignments of human alanyl tissue factor pathway inhibitor purified fromEscherichia coli have been determined. This inhibitor of the extrinsic blood coagulation pathway possesses three Kunitz-type inhibitor domains, each containing three disulfide bonds. The disulfide bond pairings in domains 1 and 3 were determined by amino acid sequencing and mass spectrometry of peptides derived from a thermolysin digest. However, thermolysin digestion did not cleave any peptide bonds within domain 2. The disulfide bond pairings in domain 2 were determined by isolating it from the thermolysin treatment and subsequently cleaving it with pepsin and trypsin into peptides which yielded the three disulfide bond pairings in this domain. These results demonstrate that the disulfide pairings in each of the three domains of human tissue factor pathway inhibitor purified fromEscherichia coli are homologous to each other and also to those in bovine pancreatic trypsin inhibitor.     

Chemical synthesis of La1 isolated from the venom of the scorpion Liocheles australasiae and determination of its disulfide bonding pattern

La1 is a 73‐residue cysteine‐rich peptide isolated from the scorpion Liocheles australasiae venom. Although La1 is the most abundant peptide in the venom, its biological function remains unknown. Here, we describe a method for efficient chemical synthesis of La1 using the native chemical ligation (NCL) strategy, in which three peptide components of less than 40 residues were sequentially ligated. The peptide thioester necessary for NCL was synthesized using an aromatic N‐acylurea approach with Fmoc‐SPPS. After completion of sequential NCL, disulfide bond formation was carried out using a dialysis method, in which the linear peptide dissolved in an acidic solution was dialyzed against a slightly alkaline buffer to obtain correctly folded La1. Next, we determined the disulfide bonding pattern of La1. Enzymatic and chemical digests of La1 without reduction of disulfide bonds were analyzed by liquid chromatography/mass spectrometry (LC/MS), which revealed two of four disulfide bond linkages. The remaining two linkages were assigned based on MS/MS analysis of a peptide fragment containing two disulfide bonds. Consequently, the disulfide bonding pattern of La1 was found to be similar to that of a von Willebrand factor type C (VWC) domain. To our knowledge, this is the first report of the experimental determination of the disulfide bonding pattern of peptides having a single VWC domain as well as their chemical synthesis. La1 synthesized in this study will be useful for investigation of its biological role in the venom. Copyright © 2015 European Peptide Society and John Wiley & Sons, Ltd.     

Subunit interaction sites between the regulatory and catalytic subunits of cAMP-dependent protein kinase. Identification of a specific interchain disulfide bond

The catalytic (C) subunit and the type II regulatory (RII) subunit of cAMP-dependent protein kinase can be cross-linked by interchain disulfide bonding. This disulfide bond can be catalyzed by cupric phenanthroline and also can be generated by a disulfide interchange using either RII-subunit or C-subunit that has been modified with either 5,5'-dithiobis(2-nitrobenzoic acid) (DTNB) or N-4(azidophenylthio)phthalimide (APTP). When the 2 cysteine residues of the C-subunit are reacted with DTNB prior to incubation with the RII-subunit, interchain disulfide bonding occurs. Similar observations are seen with C-subunit that had been modified with APTP. Interchain disulfide bonds also form when the RII-subunit is modified with DTNB prior to incubation with the C-subunit. The presence of cAMP facilitates this cross-linking while autophosphorylation of the RII-subunit retards the rate at which the interchain disulfide bond forms. Interchain disulfide bonds also form spontaneously when the RII-subunit and the C-subunit are dialyzed at pH 8.0 in the absence of reducing agents. The specific amino acid residues that participate in intersubunit disulfide bonding have been identified as Cys-97 in the RII-subunit and Cys-199 in the C-subunit. Based on the sequence homologies of the RII-subunit with other kinase substrates and on the proximity of Cys-97 to the catalytic site, a model is proposed in which the autophosphorylation site of the RII-subunit occupies the substrate-binding site in the holoenzyme. The model also proposes that this same site may be occupied by the region flanking Cys-199 in the C-subunit when the C-subunit is dissociated.     

Effects of disulfide bonds on folding behavior and mechanism of the beta-sheet protein tendamistat

Tendamistat, a small disulfide-bonded beta-sheet protein, and its three single/double-disulfide mutants are investigated by using a modified Gō-like model, aiming to understand the folding mechanism of disulfide-bonded protein as well as the effects of removal of disulfide bond on the folding process. Our simulations show that tendamistat and its two single-disulfide mutants are all two-state folders, consistent with the experimental observations. It is found that the disulfide bonds as well as three hydrogen bonds between the N-terminal loop-0 and strand-6 are of significant importance for the folding of tendamistat. Without these interactions, their two-state behaviors become unstable and the predictions of the model are inconsistent with experiments. In addition, the effect of disulfide bonds on the folding process are studied by comparing the wild-type tendamistat and its two mutants; it is found that the removal of either of the C11-C27 or C45-C73 disulfide bond leads to a large decrease in the thermodynamical stability and loss of structure in the unfolded state, and the effect of the former is stronger than that of the later. These simulation results are in good agreement with experiments and, thus, validate our model. Based on the same model, the detailed folding pathways of the wild-type tendamistat and two mutants are studied, and the effect of disulfide bonds on the folding kinetics are discussed. The obtained results provide a detailed folding picture of these proteins and complement experimental findings. Finally, the folding nuclei predicted to be existent in this protein tendamistat as well as its mutants are firstly identified in this work. The positions of the nucleus are consistent with those argued in experimental studies. Therefore, a nucleation/growth folding mechanism that can explain the two-state folding manner is clearly characterized. Moreover, the effect by the removal of each disulfide bond on the folding thermodynamics and dynamics can also be well interpreted from their influence on the folding nucleus. The implementation of this work indicates that the modified Gō-like model really describes the folding behavior of protein tendamistat and could be used to study the folding of other disulfide-bonded proteins.     

Predicting disulfide connectivity patterns

Disulfide bonds play an important role in stabilizing protein structure and regulating protein function. Therefore, the ability to infer disulfide connectivity from protein sequences will be valuable in structural modeling and functional analysis. However, to predict disulfide connectivity directly from sequences presents a challenge to computational biologists due to the nonlocal nature of disulfide bonds, i.e., the close spatial proximity of the cysteine pair that forms the disulfide bond does not necessarily imply the short sequence separation of the cysteine residues. Recently, Chen and Hwang (Proteins 2005;61:507-512) treated this problem as a multiple class classification by defining each distinct disulfide pattern as a class. They used multiple support vector machines based on a variety of sequence features to predict the disulfide patterns. Their results compare favorably with those in the literature for a benchmark dataset sharing less than 30% sequence identity. However, since the number of disulfide patterns grows rapidly when the number of disulfide bonds increases, their method performs unsatisfactorily for the cases of large number of disulfide bonds. In this work, we propose a novel method to represent disulfide connectivity in terms of cysteine pairs, instead of disulfide patterns. Since the number of bonding states of the cysteine pairs is independent of that of disulfide bonds, the problem of class explosion is avoided. The bonding states of the cysteine pairs are predicted using the support vector machines together with the genetic algorithm optimization for feature selection. The complete disulfide patterns are then determined from the connectivity matrices that are constructed from the predicted bonding states of the cysteine pairs. Our approach outperforms the current approaches in the literature.     

红细胞生长刺激蛋白二硫键定位连接方式的测定

目的:建立红细胞生长刺激蛋白(ESP)的二硫键连接方式的测定方法。方法:先将红细胞生长刺激蛋白的糖链用糖苷酶切除,再分别对ESP和还原烷基化后ESP用胰蛋白酶进行酶切,然后用MALDI-TOF测得该蛋白质的肽质量指纹图谱,通过比较还原烷基化前后各肽质量指纹图谱,找到差异肽段的分子离子峰[M+H]+,通过比对该蛋白理论酶切肽的[M+H]+,确定二硫键的连接方式。结果:ESP有4个半胱氨酸,通过比较还原烷基化前后的肽质量指纹图谱定位二硫键的位置为Cys7-Cys161和Cys29-Cys33,与理论上的二硫键连接相符。结论:建立了酶切结合质谱法测定蛋白质二硫键定位的方法,为今后生物技术产品的二硫键连接方式的质量控制提供了有效的方法。     

Assignment of disulfide bonds in proteins by fast atom bombardment mass spectrometry

A rapid and sensitive method for assignment of disulfide bonds using fast atom bombardment mass spectrometry is described for hen egg white lysozyme and bovine ribonuclease A. The protein is initially digested to a mixture of peptides using chemical and enzymatic methods under conditions which minimize disulfide bond reduction and exchange. The digested sample is analyzed directly by fast atom bombardment mass spectrometry before and after chemical reduction of cystine residues. An important feature of the method is that it is not necessary to completely resolve the peptides in the digest chromatographically prior to analysis. The disulfide-containing peptides are also characterized directly by prolonged exposure of the sample to the high energy xenon atom beam which results in the reduction of cystine residues. Intra- as well as interchain disulfide bond assignments are made on the basis of the mass difference between the molecular ions (MH+) of the oxidized and reduced peptides. Confirmation of the mass assignments may be obtained from the mass spectra of the digests after one cycle of manual Edman degradation. Although the quantity of protein required to unambiguously assign all of the disulfide linkages will depend on the ease with which the appropriate peptide fragments can be formed, results from these studies indicate that approximately 1 nmol of protein is usually sufficient.     

Chemical shift assignments of a reduced N-terminal truncation mutant of the disulfide bond isomerase TrbB from plasmid F,TrbBΔ29

TrbB from the conjugative plasmid F is a 181-residue disulfide bond isomerase that plays a role in the correct folding and maintenance of disulfide bonds within F plasmid encoded proteins in the bacterial periplasm. As a member of the thioredoxin-like superfamily, TrbB has a predicted thioredoxin-like fold that contains a C–X–X–C active site required for performing specific redox chemistries on protein substrates. Here we report the sequence-specific assignments of the reduced form of the N-terminally truncated TrbB construct, TrbBΔ29.     

Analysis of the disulfide bond arrangement of the HIV-1 envelope protein CON-S gp140 ΔCFI shows variability in the V1 and V2 regions

Disulfide bonding of cysteines is one of the most important protein modifications, and it plays a key role in establishing/maintaining protein structures in biologically active forms. Therefore, the determination of disulfide bond arrangement is one important aspect to understanding the chemical structure of a protein and defining its functional domains. Herein, aiming to understand how the HIV-1 envelope protein's structure influences its immunogenicity, we used an MS-based approach, liquid chromatography electrospray ionization Fourier transform ion cyclotron resonance (LC/ESI-FTICR) mass spectrometry, to determine the disulfide linkages on an oligomeric form of the group M consensus HIV-1 envelope protein (Env), CON-S gp140 ΔCFI. This protein has marked improvement in its immunogenicity compared to monomeric gp120 and wild-type forms of gp140 Envs. Our results demonstrate that the disulfide connectivity in the N-terminal region of CON-S gp140 ΔCFI is different from the disulfide bonding previously reported in the monomeric form of gp120 HIV-1 Env. Additionally, heterogeneity of the disulfide bonding was detected in this region. These data suggest that the V1/V2 region does not have a single, conserved disulfide bonding pattern and that variability could impact immunogenicity of expressed Envs.     

Secondary structure extensions in Pyrococcus furiosus ferredoxin destabilize the disulfide bond relative to that in other hyperthermostable ferredoxins. Global consequences for the disulfide orientational heterogeneity.

The single cubane cluster ferredoxin (Fd) from the hyperthermophilic archaeon Pyrococcus furiosus (Pf) possesses several unique properties when compared even to Fds from other hyperthermophilic archaea or bacteria. These include an equilibrium molecular heterogeneity, a six- to seven-residue increase in size, an Asp rather than the Cys as one cluster ligand, and a readily reducible disulfide bond. NMR assignments and determination of both secondary structure and tertiary contacts remote from the paramagnetic oxidized cluster of Pf 3Fe Fd with an intact disulfide bond reported previously (Teng Q., Zhou, Z. H., Smith, E. T., Busse, S. C., Howard, J. B. Adams, M. W. W., and La Mar, G. (1994) Biochemistry 33, 6316-6328) are extended here to the 4Fe oxidized cluster WT (1H and 15N) and D14C (1H only) Fds with an intact disulfide bond and to the 4Fe oxidized WT Fd (1H and 15N) with a cleaved disulfide bond. All forms are shown to possess a long (13-member) alpha-helix, two beta-sheets (one double-, one triple-stranded), and three turns outside the cluster vicinity, each with tertiary contacts among themselves as found in other Fds. While the same secondary structural elements, with similar tertiary contacts, are found in other hyperthermostable Fds, Pf Fd has two elements, the long helix and the triple-stranded beta-sheet, that exhibit extensions and form multiple tertiary contacts. All Pf Fd forms with an intact disulfide bond exhibit a dynamic equilibrium heterogeneity which is shown to modulate a hydrogen-bonding network in the hydrophobic core that radiates from the Cys21-Cys48 disulfide bond and encompasses residues Lys36, Val24, Cys21, and Cys17 and the majority of the long helix. The heterogeneity is attributed to population of the alternate S and R chiralities of the disulfide bond, each destabilized by steric interactions with the extended alpha-helix. Comparison of the chemical shifts and their temperature gradients reveals that the molecular structure of the protein with the less stable R disulfide resembles that of the Fd with a cleaved disulfide bond. Both cluster architecture (3Fe vs 4Fe) and ligand mutation (Cys for Asp14) leave the disulfide orientational heterogeneity largely unperturbed. It is concluded that the six- to seven-residue extension that results in a longer helix and larger beta-sheet in Pf Fd, relative to other hyperthermostable Fds, more likely serves to destabilize the disulfide bond, and hence make it more readily reducible, than to significantly increase protein thermostability.     

Interplay between Disulfide Bonding and N-Glycosylation Defines SLC4 Na+-coupled Transporter Extracellular Topography

The extracellular loop 3 (EL-3) of SLC4 Na⁺-coupled transporters contains 4 highly conserved cysteines and multiple N-glycosylation consensus sites. In the electrogenic Na⁺-HCO₃⁻ cotransporter NBCe1-A, EL-3 is the largest extracellular loop and is predicted to consist of 82 amino acids. To determine the structural-functional importance of the conserved cysteines and the N-glycosylation sites in NBCe1-A EL-3, we analyzed the potential interplay between EL-3 disulfide bonding and N-glycosylation and their roles in EL-3 topological folding. Our results demonstrate that the 4 highly conserved cysteines form two intramolecular disulfide bonds, Cys⁵⁸³-Cys⁵⁸⁵ and Cys⁶¹⁷-Cys⁶⁴², respectively, that constrain EL-3 in a folded conformation. The formation of the second disulfide bond is spontaneous and unaffected by the N-glycosylation state of EL-3 or the first disulfide bond, whereas formation of the first disulfide bond relies on the presence of the second disulfide bond and is affected by N-glycosylation. Importantly, EL-3 from each monomer is adjacently located at the NBCe1-A dimeric interface. When the two disulfide bonds are missing, EL-3 adopts an extended conformation highly accessible to protease digestion. This unique adjacent parallel location of two symmetrically folded EL-3 loops from each monomer resembles a domain-like structure that is potentially important for NBCe1-A function in vivo. Moreover, the formation of this unique structure is critically dependent on the finely tuned interplay between disulfide bonding and N-glycosylation in the membrane processed NBCe1-A dimer.     

Purification of correctly oxidized MHC class I heavy-chain molecules under denaturing conditions: a novel strategy exploiting disulfide assisted protein folding

The aim of this study has been to develop a strategy for purifying correctly oxidized denatured major histocompability complex class I (MHC-I) heavy-chain molecules, which on dilution, fold efficiently and become functional. Expression of heavy-chain molecules in bacteria results in the formation of insoluble cellular inclusion bodies, which must be solubilized under denaturing conditions. Their subsequent purification and refolding is complicated by the fact that (1). correct folding can only take place in combined presence of beta(2)-microglobulin and a binding peptide; and (2). optimal in vitro conditions for disulfide bond formation ( approximately pH 8) and peptide binding ( approximately pH 6.6) are far from complementary. Here we present a two-step strategy, which relies on uncoupling the events of disulfide bond formation and peptide binding. In the first phase, heavy-chain molecules with correct disulfide bonding are formed under non-reducing denaturing conditions and separated from scrambled disulfide bond forms by hydrophobic interaction chromatography. In the second step, rapid refolding of the oxidized heavy chains is afforded by disulfide bond-assisted folding in the presence of beta(2)-microglobulin and a specific peptide. Under conditions optimized for peptide binding, refolding and simultaneous peptide binding of the correctly oxidized heavy chain was much more efficient than that of the fully reduced molecule.     

Identification of disulfide bonds among the nine core 2 N-acetylglucosaminyltransferase-M cysteines conserved in the mucin beta6-N-acetylglucosaminyltransferase family

Bovine core 2 beta1,6-N-acetylglucosaminyltransferase-M (bC2GnT-M) catalyzes the formation of all mucin beta1,6-N-acetylglucosaminides, including core 2, core 4, and blood group I structures. These structures expand the complexity of mucin carbohydrate structure and thus the functional potential of mucins. The four known mucin beta1,6-N-acetylglucosaminyltransferases contain nine conserved cysteines. We determined the disulfide bond assignments of these cysteines in [(35)S]cysteine-labeled bC2GnT-M isolated from the serum-free conditioned medium of Chinese hamster ovary cells stably transfected with a pSecTag plasmid. This plasmid contains bC2GnT-M cDNA devoid of the 5'-sequence coding the cytoplasmic tail and transmembrane domain. The C18 reversed phase high performance liquid chromatographic profile of the tryptic peptides of reduced-alkylated (35)S-labeled C2GnT-M was established using microsequencing. Each cystine pair was identified by rechromatography of the C8 high performance liquid chromatographic radiolabeled tryptic peptides of alkylated bC2GnT-M on C18 column. Among the conserved cysteines in bC2GnT-M, the second (Cys(113)) was a free thiol, whereas the other eight cysteines formed four disulfide bridges, which included the first (Cys(73)) and sixth (Cys(230)), third (Cys(164)) and seventh (Cys(384)), fourth (Cys(185)) and fifth (Cys(212)), and eighth (Cys(393)) and ninth (Cys(425)) cysteine residues. This pattern of disulfide bond formation differs from that of mouse C2GnT-L, which may contribute to the difference in substrate specificity between these two enzymes. Molecular modeling using disulfide bond assignments and the fold recognition/threading method to search the Protein Data Bank found a match with aspartate aminotransferase structure. This structure is different from the two major protein folds proposed for glycosyltransferases.     

Identification of an atypical membrane protein involved in the formation of protein disulfide bonds in oxygenic photosynthetic organisms

The evolution of oxygenic photosynthesis in cyanobacteria nearly three billion years ago provided abundant reducing power and facilitated the elaboration of numerous oxygen-dependent reactions in our biosphere. Cyanobacteria contain an internal thylakoid membrane system, the site of photosynthesis, and a typical Gram-negative envelope membrane system. Like other organisms, the extracytoplasmic space in cyanobacteria houses numerous cysteine-containing proteins. However, the existence of a biochemical system for disulfide bond formation in cyanobacteria remains to be determined. Extracytoplasmic disulfide bond formation in non-photosynthetic organisms is catalyzed by coordinated interaction between two proteins, a disulfide carrier and a disulfide generator. Here we describe a novel gene, SyndsbAB, required for disulfide bond formation in the extracytoplasmic space of cyanobacteria. The SynDsbAB orthologs are present in most cyanobacteria and chloroplasts of higher plants with fully sequenced genomes. The SynDsbAB protein contains two distinct catalytic domains that display significant similarity to proteins involved in disulfide bond formation in Escherichia coli and eukaryotes. Importantly, SyndsbAB complements E. coli strains defective in disulfide bond formation. In addition, the activity of E. coli alkaline phosphatase localized to the periplasm of Synechocystis 6803 is dependent on the function of SynDsbAB. Deletion of SyndsbAB in Synechocystis 6803 causes significant growth impairment under photoautotrophic conditions and results in hyper-sensitivity to dithiothreitol, a reductant, whereas diamide, an oxidant had no effect on the growth of the mutant strains. We conclude that SynDsbAB is a critical protein for disulfide bond formation in oxygenic photosynthetic organisms and required for their optimal photoautotrophic growth.     

Connexon integrity is maintained by non-covalent bonds: intramolecular disulfide bonds link the extracellular domains in rat connexin-43.

Rat heart connexin-43 (RCx43) has been isolated using a modified procedure that is rapid and can be used with fresh or frozen hearts. When RCx43 is isolated in the presence of the alkylating reagent iodoacetamide, no intermolecular disulfide bonds are found. However, the alkylated RCx43 does have at least one intramolecular disulfide bond. By using site directed antibodies and proteolytic cleavage the location of the intramolecular disulfide bonding is shown between the two extracellular loops of RCx43.     

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.