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.