Identification of Allosteric Disulfides from Prestress Analysis

Disulfide bonds serve to form physical cross-links between residues in protein structures, thereby stabilizing the protein fold. Apart from this purely structural role, they can also be chemically active, participating in redox reactions, and they may even potentially act as allosteric switches controlling protein functions. Specific types of disulfide bonds have been identified in static protein structures from their distinctive pattern of dihedral bond angles, and the allosteric function of such bonds is purported to be related to the torsional strain they store. Using all-atom molecular-dynamics simulations for ∼700 disulfide bonded proteins, we analyzed the intramolecular mechanical forces in 20 classes of disulfide bonds. We found that two particular classes, the −RHStaple and the −/+RHHook disulfides, are indeed more stressed than other disulfide bonds, but the stress is carried primarily by stretching of the S-S bond and bending of the neighboring bond angles, rather than by dihedral torsion. This stress corresponds to a tension force of magnitude ∼200 pN, which is balanced by repulsive van der Waals interactions between the cysteine Cα atoms. We confirm stretching of the S-S bond to be a general feature of the −RHStaples and the −/+RHHooks by analyzing ∼20,000 static protein structures. Given that forced stretching of S-S bonds is known to accelerate their cleavage, we propose that prestress of allosteric disulfide bonds has the potential to alter the reactivity of a disulfide, thereby allowing us to readily switch between functional states.     

Structure and conformation of the disulfide bond in dimeric lung surfactant peptides SP-B1-25 and SP-B8-25

Raman spectroscopy was used to determine the conformation of the disulfide linkage between cysteine residues in the homodimeric construct of the N-terminal alpha helical domain of surfactant protein B (dSP-B_1-25). The conformation of the disulfide bond between cysteine residues in position 8 of the homodimer of dSP-B_1-25 was compared with that of a truncated homodimer (dSP-B_8-25) of the peptide having a disulfide linkage at the same position in the alpha helix. Temperature-dependent Raman spectra of the S-S stretching region centered at ∼ 500 cm^− 1 indicated a stable, although highly strained disulfide conformation with a χ(CS-SC) dihedral angle of ± 10° for the dSP-B_1-25 dimer. In contrast, the truncated dimer dSP-B_8-25 exhibited a series of disulfide conformations with the χ(CS-SC) dihedral angle taking on values of either ± 30° or 85± 20°. For conformations with χ(CS-SC) close to the ± 90° value, the Raman spectra of the 8-25 truncated dimers exhibited χ(SS-CC) dihedral angles of 90/180° and 20-30°. In the presence of a lipid mixture, both constructs showed a ν(S-S) band at ∼ 488 cm^− 1, corresponding to a χ(CS-SC) dihedral angle of ± 10°. Polarized infrared spectroscopy was also used to determine the orientation of the helix and β-sheet portion of both synthetic peptides. These calculations indicated that the helix was oriented primarily in the plane of the surface, at an angle of ∼ 60-70° to the surface normal, while the β structure had ∼ 40° tilt. This orientation direction did not change in the presence of a lipid mixture or with temperature. These observations suggest that: (i) the conformational flexibility of the disulfide linkage is dependent on the amino acid residues that flank the cysteine disulfide bond, and (ii) in both constructs, the presence of a lipid matrix locks the disulfide bond into a preferred conformation.     

Allosteric disulfide bonds

Disulfide bonds have been generally considered to be either structural or catalytic. Structural bonds stabilize a protein, while catalytic bonds mediate thiol-disulfide interchange reactions in substrate proteins. There is emerging evidence for a third type of disulfide bond that can control protein function by triggering a conformational change when it breaks and/or forms. These bonds can be thought of as allosteric disulfides. To better define the properties of allosteric disulfides, we have analyzed the geometry and dihedral strain of 6874 unique disulfide bonds in 2776 X-ray structures. A total of 20 types of disulfide bonds were identified in the dataset based on the sign of the five chi angles that make up the bond. The known allosteric disulfides were all contained in 1 of the 20 groups, the -RHStaple bonds. This bond group has a high mean potential energy and narrow energy distribution, which is consistent with a functional role. We suggest that the -RHStaple configuration is a hallmark of allosteric disulfides. About 1 in 15 of all structurally determined disulfides is a potential allosteric bond.     

One-Way Allosteric Communication between the Two Disulfide Bonds in Tissue Factor

Tissue factor (TF) is a transmembrane glycoprotein that plays distinct roles in the initiation of extrinsic coagulation cascade and thrombosis. TF contains two disulfide bonds, one each in the N-terminal and C-terminal extracellular domains. The C-domain disulfide, Cys186-Cys209, has a ?RHStaple configuration in crystal structures, suggesting that this disulfide carries high pre-stress. The redox state of this disulfide has been proposed to regulate TF encryption/decryption. Ablating the N-domain Cys49-Cys57 disulfide bond was found to increase the redox potential of the Cys186-Cys209 bond, implying an allosteric communication between the domains. Using molecular dynamics simulations, we observed that the Cys186-Cys209 disulfide bond retained the ?RHStaple configuration, whereas the Cys49-Cys57 disulfide bond fluctuated widely. The Cys186-Cys209 bond featured the typical ?RHStaple disulfide properties, such as a longer S-S bond length, larger C-S-S angles, and higher bonded prestress, in comparison to the Cys49-Cys57 bond. Force distribution analysis was used to sense the subtle structural changes upon ablating the disulfide bonds, and allowed us to identify a one-way allosteric communication mechanism from the N-terminal to the C-terminal domain. We propose a force propagation pathway using a shortest-pathway algorithm, which we suggest is a useful method for searching allosteric signal transduction pathways in proteins. As a possible explanation for the pathway being one-way, we identified a pronounced lower degree of conformational fluctuation, or effectively higher stiffness, in the N-terminal domain. Thus, the changes of the rigid domain (N-terminal domain) can induce mechanical force propagation to the soft domain (C-terminal domain), but not vice versa.     

Structure and conformation of the disulfide bond in dimeric lung surfactant peptides SP-B1-25 and SP-B8-25

Raman spectroscopy was used to determine the conformation of the disulfide linkage between cysteine residues in the homodimeric construct of the N-terminal alpha helical domain of surfactant protein B (dSP-B(1-25)). The conformation of the disulfide bond between cysteine residues in position 8 of the homodimer of dSP-B(1-25) was compared with that of a truncated homodimer (dSP-B(8-25)) of the peptide having a disulfide linkage at the same position in the alpha helix. Temperature-dependent Raman spectra of the S-S stretching region centered at approximately 500 cm(-1) indicated a stable, although highly strained disulfide conformation with a chi(CS-SC) dihedral angle of +/-10 degrees for the dSP-B(1-25) dimer. In contrast, the truncated dimer dSP-B(8-25) exhibited a series of disulfide conformations with the chi(CS-SC) dihedral angle taking on values of either +/-30 degrees or 85+/-20 degrees . For conformations with chi(CS-SC) close to the +/-90 degrees value, the Raman spectra of the 8-25 truncated dimers exhibited chi(SS-CC) dihedral angles of 90/180 degrees and 20-30 degrees . In the presence of a lipid mixture, both constructs showed a nu(S-S) band at approximately 488 cm(-1), corresponding to a chi(CS-SC) dihedral angle of +/-10 degrees . Polarized infrared spectroscopy was also used to determine the orientation of the helix and beta-sheet portion of both synthetic peptides. These calculations indicated that the helix was oriented primarily in the plane of the surface, at an angle of approximately 60-70 degrees to the surface normal, while the beta structure had approximately 40 degrees tilt. This orientation direction did not change in the presence of a lipid mixture or with temperature. These observations suggest that: (i) the conformational flexibility of the disulfide linkage is dependent on the amino acid residues that flank the cysteine disulfide bond, and (ii) in both constructs, the presence of a lipid matrix locks the disulfide bond into a preferred conformation.     

The crystallographically determined structures of atypical strained disulfides engineered into subtilisin

The geometries of two disulfide bridges genetically engineered into subtilisin have been characterized by x-ray crystallography to determine the structural and energetic constraints involved in introducing disulfide bonds into proteins. Both disulfide bridges (Cys-24-Cys-87 and Cys-22-Cys-87) exhibit atypical sets of dihedral angles compared to those for other reported disulfide structures in proteins. The geometric trends for naturally occurring disulfides in protein crystal structures are examined. Comparison of the disulfide-containing mutant protein structures with the wild-type structure shows that, in both cases, disulfide incorporation is accommodated by relatively minor changes in local main-chain conformation. The Cys-22-Cys-87 disulfide has two high energy dihedral angles (X2 = 121 degrees, X2' = 143 degrees). Both disulfides produce short non-bonded contacts with the main-chain.     

Characterization of a Reduced Form of Plasma Plasminogen as the Precursor for Angiostatin Formation

Plasma plasminogen is the precursor of the tumor angiogenesis inhibitor, angiostatin. Generation of angiostatin in blood involves activation of plasminogen to the serine protease plasmin and facilitated cleavage of two disulfide bonds and up to three peptide bonds in the kringle 5 domain of the protein. The mechanism of reduction of the two allosteric disulfides has been explored in this study. Using thiol-alkylating agents, mass spectrometry, and an assay for angiostatin formation, we show that the Cys⁴⁶²-Cys⁵⁴¹ disulfide bond is already cleaved in a fraction of plasma plasminogen and that this reduced plasminogen is the precursor for angiostatin formation. From the crystal structure of plasminogen, we propose that plasmin ligands such as phosphoglycerate kinase induce a conformational change in reduced kringle 5 that leads to attack by the Cys⁵⁴¹ thiolate anion on the Cys⁵³⁶ sulfur atom of the Cys⁵¹²-Cys⁵³⁶ disulfide bond, resulting in reduction of the bond by thiol/disulfide exchange. Cleavage of the Cys⁵¹²-Cys⁵³⁶ allosteric disulfide allows further conformational change and exposure of the peptide backbone to proteolysis and angiostatin release. The Cys⁴⁶²-Cys⁵⁴¹ and Cys⁵¹²-Cys⁵³⁶ disulfides have −/+RHHook and −LHHook configurations, respectively, which are two of the 20 different measures of the geometry of a disulfide bond. Analysis of the structures of the known allosteric disulfide bonds identified six other bonds that have these configurations, and they share some functional similarities with the plasminogen disulfides. This suggests that the −/+RHHook and −LHHook disulfides, along with the −RHStaple bond, are potential allosteric configurations.     

Search for allosteric disulfide bonds in NMR structures

Background  

Allosteric disulfide bonds regulate protein function when they break and/or form. They typically have a -RHStaple configuration, which is defined by the sign of the five chi angles that make up the disulfide bond.     

Real-time Monitoring of Intermediates Reveals the Reaction Pathway in the Thiol-Disulfide Exchange between Disulfide Bond Formation Protein A (DsbA) and B (DsbB) on a Membrane-immobilized Quartz Crystal Microbalance (QCM) System

Disulfide bond formation protein B (DsbB_S-S,S-S) is an inner membrane protein in Escherichia coli that has two disulfide bonds (S-S, S-S) that play a role in oxidization of a pair of cysteine residues (SH, SH) in disulfide bond formation protein A (DsbA_SH,SH). The oxidized DsbA_S-S, with one disulfide bond (S-S), can oxidize proteins with SH groups for maturation of a folding preprotein. Here, we have described the transient kinetics of the oxidation reaction between DsbA_SH,SH and DsbB_S-S,S-S. We immobilized DsbB_S-S,S-S embedded in lipid bilayers on the surface of a 27-MHz quartz crystal microbalance (QCM) device to detect both formation and degradation of the reaction intermediate (DsbA-DsbB), formed via intermolecular disulfide bonds, as a mass change in real time. The obtained kinetic parameters (intermediate formation, reverse, and oxidation rate constants (k_f, k_r, and k_cat, respectively) indicated that the two pairs of cysteine residues in DsbB_S-S,S-S were more important for the stability of the DsbA-DsbB intermediate than ubiquinone, an electron acceptor for DsbB_S-S,S-S. Our data suggested that the reaction pathway of almost all DsbA_SH,SH oxidation processes would proceed through this stable intermediate, avoiding the requirement for ubiquinone.     

Reduced representation model of protein structure prediction: statistical potential and genetic algorithms.

A reduced representation model, which has been described in previous reports, was used to predict the folded structures of proteins from their primary sequences and random starting conformations. The molecular structure of each protein has been reduced to its backbone atoms (with ideal fixed bond lengths and valence angles) and each side chain approximated by a single virtual united-atom. The coordinate variables were the backbone dihedral angles phi and psi. A statistical potential function, which included local and nonlocal interactions and was computed from known protein structures, was used in the structure minimization. A novel approach, employing the concepts of genetic algorithms, has been developed to simultaneously optimize a population of conformations. With the information of primary sequence and the radius of gyration of the crystal structure only, and starting from randomly generated initial conformations, I have been able to fold melittin, a protein of 26 residues, with high computational convergence. The computed structures have a root mean square error of 1.66 A (distance matrix error = 0.99 A) on average to the crystal structure. Similar results for avian pancreatic polypeptide inhibitor, a protein of 36 residues, are obtained. Application of the method to apamin, an 18-residue polypeptide with two disulfide bonds, shows that it folds apamin to native-like conformations with the correct disulfide bonds formed.     

Redox Regulation of Methionine Aminopeptidase 2 Activity

Protein translation is initiated with methionine in eukaryotes, and the majority of proteins have their N-terminal methionine removed by methionine aminopeptidases (MetAP1 and MetAP2) prior to action. Methionine removal can be important for protein function, localization, or stability. No mechanism of regulation of MetAP activity has been identified. MetAP2, but not MetAP1, contains a single Cys²²⁸-Cys⁴⁴⁸ disulfide bond that has an −RHStaple configuration and links two β-loop structures, which are hallmarks of allosteric disulfide bonds. From analysis of crystal structures and using mass spectrometry and activity assays, we found that the disulfide bond exists in oxidized and reduced states in the recombinant enzyme. The disulfide has a standard redox potential of −261 mV and is efficiently reduced by the protein reductant, thioredoxin, with a rate constant of 16,180 m⁻¹ s⁻¹. The MetAP2 disulfide bond also exists in oxidized and reduced states in glioblastoma tumor cells, and stressing the cells by oxygen or glucose deprivation results in more oxidized enzyme. The Cys²²⁸-Cys⁴⁴⁸ disulfide is at the rim of the active site and is only three residues distant from the catalytic His²³¹, which suggested that cleavage of the bond would influence substrate hydrolysis. Indeed, oxidized and reduced isoforms have different catalytic efficiencies for hydrolysis of MetAP2 peptide substrates. These findings indicate that MetAP2 is post-translationally regulated by an allosteric disulfide bond, which controls substrate specificity and catalytic efficiency.     

Role of the propeptide in controlling conformation and assembly state of hepatitis B virus e-antigen

Hepatitis B virus “e-antigen” (HBeAg) is thought to be a soluble dimeric protein that is associated with chronic infection. It shares 149 residues with the viral capsid protein “core-antigen” (HBcAg), but has an additional 10-residue, hydrophobic, cysteine-containing amino-terminal propeptide whose presence correlates with physical, serological, and immunological differences between the two proteins. In HBcAg dimers, the subunits pair by forming a four-helix bundle stabilized by an intermolecular disulfide bond. The structure of HBeAg is probably similar but, instead, has two intramolecular disulfide bonds involving the propeptide. To compare the proteins directly and thereby clarify the role of the propeptide, we identified mutations and solution conditions that render both proteins as either soluble dimers or assembled capsids. Thermally induced unfolding monitored by circular dichroism, and electrophoresis of oxidized and reduced dimers, showed that the propeptide has a destabilizing effect and that the intramolecular disulfide bond forms preferentially and blocks the formation of the intermolecular disulfide bond that otherwise stabilizes the dimer. The HBeAg capsids are less regular than the HBcAg capsids; nevertheless, cryo-electron microscopy reconstructions confirm that they are constructed of dimers resembling those of HBcAg capsids. In them, a portion of the propeptide is visible near the dimer interface, suggesting that it intercalates there, consistent with the known formation of a disulfide bond between C(− 7) in the propeptide and C61 in the dimer interface. However, this intercalation distorts the dimer into an assembly-reluctant conformation.     

Hierarchical formation of disulfide bonds in the immunoglobulin Fc fragment is assisted by protein-disulfide isomerase

Antibodies provide an excellent system to study the folding and assembly of all beta-sheet proteins and to elucidate the hierarchy of intra/inter chain disulfide bonds formation during the folding process of multimeric and multidomain proteins. Here, the folding process of the Fc fragment of the heavy chain of the antibody MAK33 was investigated. The Fc fragment consists of the C(H)3 and C(H)2 domains of the immunoglobulin heavy chain, both containing a single S-S bond. The folding process was investigated both in the absence and presence of the folding catalyst protein-disulfide isomerase (PDI), monitoring the evolution of intermediates by electrospray mass spectrometry. Moreover, the disulfide bonds present at different times in the folding mixture were identified by mass mapping to determine the hierarchy of disulfide bond formation. The analysis of the uncatalyzed folding showed that the species containing one intramolecular disulfide predominated throughout the entire process, whereas the fully oxidized Fc fragment never accumulated in significant amounts. This result suggests the presence of a kinetic trap during the Fc folding, preventing the one-disulfide-containing species (1S2H) to reach the fully oxidized protein (2S). The assignment of disulfide bonds revealed that 1S2H is a homogeneous species characterized by the presence of a single disulfide bond (Cys-130-Cys-188) belonging to the C(H)3 domain. When the folding experiments were carried out in the presence of PDI, the completely oxidized species accumulated and predominated at later stages of the process. This species contained the two native S-S bonds of the Fc protein. Our results indicate that the two domains of the Fc fragment fold independently, with a precise hierarchy of disulfide formation in which the disulfide bond, especially, of the C(H)2 domain requires catalysis by PDI.     

Disulfide chromophore and its optical activity

The compounds I-IV derived from α-D-cyclodextrin moiety by bridging and/or interconnecting with various patterns of disulfide bonds were chosen as models for the spectroscopic study of conformation of the disulfide bridge. The energy gap between the disulfide and cyclodextrin's electronic transitions allows us to investigate absorption and electronic circular dichroism spectra without disturbing spectral overlaps with amides or aromatic amino acids in peptides or proteins. Raman optical activity (ROA) spectra were measured and the bands due to S-S and C-S stretching motion identified. Comparison with the quantum mechanical calculations of simple models indicates that sense of disulfide twist follows sign of the measured S-S ROA band.     

Impact of an easily reducible disulfide bond on the oxidative folding rate of multi-disulfide-containing proteins.

The burial of native disulfide bonds, formed within stable structure in the regeneration of multi-disulfide-containing proteins from their fully reduced states, is a key step in the folding process, as the burial greatly accelerates the oxidative folding rate of the protein by sequestering the native disulfide bonds from thiol-disulfide exchange reactions. Nevertheless, several proteins retain solvent-exposed disulfide bonds in their native structures. Here, we have examined the impact of an easily reducible native disulfide bond on the oxidative folding rate of a protein. Our studies reveal that the susceptibility of the (40-95) disulfide bond of Y92G bovine pancreatic ribonuclease A (RNase A) to reduction results in a reduced rate of oxidative regeneration, compared with wild-type RNase A. In the native state of RNase A, Tyr 92 lies atop its (40-95) disulfide bond, effectively shielding this bond from the reducing agent, thereby promoting protein oxidative regeneration. Our work sheds light on the unique contribution of a local structural element in promoting the oxidative folding of a multi-disulfide-containing protein.     

The timing of hamster sperm nuclear decondensation and male pronucleus formation is related to sperm nuclear disulfide bond content

The relationship between the timing of both sperm nuclear decondensation and male pronucleus formation in the oocyte and the relative level of disulfide bonds within the sperm nucleus was evaluated. Since reduction of sperm nuclear disulfide (S-S) bonds is a prerequisite for sperm nuclear decondensation in vitro and in vivo, we hypothesized that sperm nuclei with relatively few S-S bonds would require less time to decondense in the oocyte than sperm nuclei with higher numbers of S-S bonds, and that male pronucleus formation would occur more rapidly as well. Four types of hamster sperm nuclei, in which the extent of S-S bonding differed, were microinjected into hamster oocytes, and the time course of sperm nuclear decondensation and male pronucleus formation was charted. Cauda epididymal sperm nuclei, which are rich in S-S bonds, required 45-60 min to decondense. In contrast, nuclei containing few S-S bonds (namely sonication-resistant spermatid nuclei and cauda epididymal sperm nuclei treated in vitro with the S-S bond-reducing agent dithiothreitol) decondensed within 5-10 min of microinjection. Caput epididymal sperm nuclei, with intermediate S-S bond content, decondensed in 10-20 min. Regardless of when decondensation occurred, formation of the male pronucleus never preceded that of the female pronucleus, which occurred 1.25-1.5 h after microinjection. However, sperm nuclei with few S-S bonds were more likely than S-S rich nuclei to transform into male pronuclei in synchrony with the formation of the female pronucleus. We conclude that the timing sperm nuclear decondensation and pronucleus formation depends in part upon the S-S bond content of the sperm nucleus.     

Analysis of side-chain conformational distributions in neutrophil peptide-5 NMR structures

The side-chain conformations have been analyzed in the antimicrobial peptide, Neutrophil Peptide-5 (NP-5), whose structure was independently generated from nmr-derived distance constraints using a distance geometry algorithm. The side-chain and peptide dihedral angle distributions in the nmr structures were compared with those constructed from a data base of high-resolution protein crystal structures. The side-chain conformational preferences for NP-5 in solution are significantly different from those observed in the crystal structure data base. These results indicate that the side-chain conformations are quite disordered for many of the residues of NP-5. The absence of a correlation between the width of the conformational distribution and surface accessibility suggests that the disorder may be due to limitations in the structural information extracted from the nmr data rather than to molecular motion. However, it is also observed that the degree of conformational disorder is only weakly correlated with the number of nuclear Overhauser enhancements to a given side chain. Possible reasons for this are discussed. Molecular mechanics refinement of these structures did not significantly change the side-chain populations. Anomolously wide distributions are observed for rotations about the peptide bonds and the disulfide bonds in the NP-5 distance geometry structures, which are improved by the refinement. The very high degree of order observed for the central dihedral angle of the disulfide bond in the high-resolution crystal data base suggests that the rotation about this bond in proteins is determined by the local potential.     

Carbohydrates Facilitate Correct Disulfide Bond Formation and Folding of Rotavirus VP7

It is well established that glycosylation is essential for assembly of enveloped viruses, but no information is yet available as to the function of carbohydrates on the nonenveloped but glycosylated rotavirus. We show that tunicamycin and, more pronouncedly, a combination of tunicamycin and brefeldin A treatment caused misfolding of the luminal VP7 protein, leading to interdisulfide bond aggregation. While formation of VP7 aggregates could be prevented under reducing conditions, they reoccurred in less than 30 min after a shift to an oxidizing milieu. Furthermore, while glycosylated VP7 interacted during maturation with protein disulfide isomerase, nonglycosylated VP7 did not, suggesting that glycosylation is a prerequisite for protein disulfide isomerase interaction. While native NSP4, which does not possess S-S bonds, was not dependent on N-linked glycosylation or on protein disulfide isomerase assistance for maturation, nonglycosylated NSP4 was surprisingly found to interact with protein disulfide isomerase, further suggesting that protein disulfide isomerase can act both as an enzyme and as a chaperone. In conclusion, our data suggest that the major function of carbohydrates on VP7 is to facilitate correct disulfide bond formation and protein folding.     

Enzyme reduction of disulfide bonds by thioredoxin. The reactivity of disulfide bonds in human choriogonadotropin and its subunits.

The NADPH-dependent enzymic reduction of disulfide bonds in human choriogonadotropin and its two subunits, alpha and beta, was examined with thioredoxin and thioredoxin reductase from Escherichia coli. With 12 muM thioredoxin and 0.1 muM thioredoxin reductase at pH 7 all disulfide bonds in the alpha subunit could be reduced in 15 min. The reduction of disulfide bonds was recorded by a simple spectrophotometric assay at 340 nm, which allowed quantitation of the reduction rate and the number of disulfide bonds reduced. Partial reduction of the alpha subunit with thioredoxin followed by S-carboxymethylation with iodol[2-3H]acetic acid and analysis of tryptic peptides indicated that all S-S bonds in the alpha subunit were surface oriented and equally reactive. The usefulness of thioredoxin reduction of disulfide bonds as a chemical probe of protein structure was shown by the much slower reaction of disulfide bonds in the intact hormone as compared to its two biologically inactive subunits.     

Additional disulfide bonds in insulin: Prediction,recombinant expression,receptor binding affinity,and stability

The structure of insulin, a glucose homeostasis-controlling hormone, is highly conserved in all vertebrates and stabilized by three disulfide bonds. Recently, we designed a novel insulin analogue containing a fourth disulfide bond located between positions A10-B4. The N-terminus of insulin''s B-chain is flexible and can adapt multiple conformations. We examined how well disulfide bond predictions algorithms could identify disulfide bonds in this region of insulin. In order to identify stable insulin analogues with additional disulfide bonds, which could be expressed, the C_β cut-off distance had to be increased in many instances and single X-ray structures as well as structures from MD simulations had to be used. The analogues that were identified by the algorithm without extensive adjustments of the prediction parameters were more thermally stable as assessed by DSC and CD and expressed in higher yields in comparison to analogues with additional disulfide bonds that were more difficult to predict. In contrast, addition of the fourth disulfide bond rendered all analogues resistant to fibrillation under stress conditions and all stable analogues bound to the insulin receptor with picomolar affinities. Thus activity and fibrillation propensity did not correlate with the results from the prediction algorithm.