Crystal structures of subtilisin BPN' variants containing disulfide bonds and cavities: concerted structural rearrangements induced by mutagenesis

The X-ray structures of four genetically engineered disulfide variants of subtilisin have been analyzed to determine the energetic and structural constraints involved in inserting disulfide bonds into proteins. Each of the engineered disulfides exhibited atypical sets of dihedral angles compared with known structures of natural disulfide bridges in proteins and affected its local structural environment to a different extent. The disulfides located in buried regions, Cys26-Cys232 and Cys29-Cys119, induced larger changes than did Cys24-Cys87 and Cys22-Cys87, which are located on the surface of the molecule. An analysis of the concerted changes in secondary structure units such as alpha-helices and beta-sheets indicated systematic long-range effects. The observed changes in the mutants were largely distributed asymmetrically around the inserted disulfides, reflecting different degrees of inherent flexibility of neighboring secondary structure types. The disulfide substitution in each variant molecule created some invaginations or cavities, causing a reorganization of the surrounding water structure. These changes are described, as well as the changes in side chain positions of groups that border the cavities.     

Three-dimensional structure of the kringle sequence: structure of prothrombin fragment 1

The three-dimensional structure of bovine prothrombin fragment 1 has been solved at 2.8-A resolution. The electron density clearly reveals four disulfide bridges along with more than 80% of the side chains completely in density, which correspond faithfully to the kringle sequence, its preceding 30 residues, and the dodecapeptide carboxy terminal; the polysaccharide and the first 35 residues of the amino terminal of fragment 1 are disordered or about 40% of the structure. The folding of the kringle sequence is based upon close disulfide van der Waals contacts between Cys-87-Cys-127 and Cys-115-Cys-139 (4.1 A between midpoints of the bridges), two antiparallel strands of highly conserved (113-118, 124-129) beta-structure, and the stacking of some conserved aromatic residues, all near the center of the folded structure. Moreover, the overall folding appears to be duplicated as a pair of stacked duplex loops with an antiparallel open loop. The overall shape of the kringle structure approximates an eccentric oblate ellipsoid of dimensions 11 X 28 X 30 A. The residues immediately preceding the kringle are dominated by alpha-helical structure (Phe-41-Cys-48; Leu-56-Glu-63). Residues Phe-41-Trp-42 and Tyr-45, which are conserved in factor IX, factor X, protein C, and protein Z, form another aromatic stacked cluster while the Cys-48-Cys-61 disulfide loop corresponds to the well-known alpha/beta structural unit. The dodecapeptide carboxy-terminal interkringle chain extends along the periphery of the kringle in its plane and forms a beta-structure with the kringle-closing Ser-140-Val-143 tetrapeptide.     

Disulfide structure of the pheromone binding protein from the silkworm moth, Bombyx mori

Disulfide bond formation is the only known posttranslational modification of insect pheromone binding proteins (PBPs). In the PBPs from moths (Lepidoptera), six cysteine residues are highly conserved at positions 19, 50, 54, 97, 108 and 117, but to date nothing is known about their respective linkage or redox status. We used a multiple approach of enzymatic digestion, chemical cleavage, partial reduction with Tris-(2-carboxyethyl)phosphine, followed by digestion with endoproteinase Lys-C to determine the disulfide connectivity in the PBP from Bombyx mori (BmPBP). Identification of the reaction products by on-line liquid chromatography-electrospray ionization mass spectrometry (LC/ESI-MS) and protein sequencing supported the assignment of disulfide bridges at Cys-19-Cys-54, Cys-50-Cys-108 and Cys-97-Cys-117. The disulfide linkages were identical in the protein obtained by periplasmic expression in Escherichia coli and in the native BmPBP.     

Determination of disulfide array and subunit structure of taste-modifying protein, miraculin

The taste-modifying protein, miraculin (Theerasilp, S. et al. (1989) J. Biol. Chem. 264, 6655-6659) has seven cysteine residues in a molecule composed of 191 amino acid residues. The formation of three intrachain disulfide bridges at Cys-47-Cys-92, Cys-148-Cys-159 and Cys-152-Cys-155 and one interchain disulfide bridge at Cys-138 was determined by amino acid sequencing and composition analysis of cystine-containing peptides isolated by HPLC. The presence of an interchain disulfide bridge was also supported by the fact that the cystine peptide containing Cys-138 showed a negative color test for the free sulfhydryl group and a positive test after reduction with dithiothreitol. The molecular mass of non-reduced miraculin (43 kDa) in sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) was nearly twice the calculated molecular mass based on the amino acid sequence and the carbohydrate content of reduced miraculin (25 kDa). The molecular mass of native miraculin determined by low-angle laser light scattering was 90 kDa. Application of a crude extract of miraculin to a Sephadex G-75 column indicated that the taste-modifying activity appears at 52 kDa. It was concluded that native miraculin in pure form is a tetramer of the 25 kDa-peptide and native miraculin in crude state or denatured, non-reduced miraculin in pure form is a dimer of the peptide. Both tetramer miraculin and native dimer miraculin in crude state had the taste-modifying activity.     

The molecular structure of a dimer composed of the variable portions of the Bence-Jones protein REI refined at 2.0-A resolution.

The structure of the variable portions of a K-type Bence-Jones protein REI forming a dimer has been determined by X-ray diffraction to a resolution of 2.0 A. The structure has been refined using a constrained crystallographic refinement procedure. The final R value is 0.24 for 15000 significantly measured reflections; the estimated standard deviation of atomic positions is 0.09 A. A more objective assessment of the error in the atomic positions is possible by comparing the two independently refined monomers. The mean deviation of main-chain atoms of the two chains in internal segments in 0.22 A, of main-chain dihedral angles 6.3 degrees for these segments. The unrefined molecular structure of the VREI dimer has been published (Epp, O., Colman, P., Fehlhammer, H., Bode, W., Schiffer, M., Huber, R., and Palm, W. (1974), Eur. J. Biochem. 45, 513). Now a detailed analysis is presented in terms of hydrogen bonds and conformational angles. Secondary structural elements (antiparallel beta structure, reverse turns) are defined. A more precise atomic arrangement of the amino acid residues forming the contact region and the hapten binding site is given as well as the localization of solvent molecules. Two cis-prolines (Pro-8 and Pro-95) were detected. The intrachain disulfide bridge (Cys-23-Cys-88) occurs statistically in two alternative conformations. The structure suggests reasons for strong conservation of several amino acid residues. The knowledge of the refined molecular structure enables crystal structure analyses of related molecules to be made by Patterson search techniques. The calculated phases based on the refined structure are much improved compared to isomorphous phases. Therefore the effects of hapten binding on the molecular structure can be analyzed by the difference Fourier technique with more reliability. Hapten binding studies have been started.     

Role of disulfide bridges in the folding, structure and biological activity of omega-conotoxin GVIA.

Omega-Conotoxin GVIA (GVIA), an N-type calcium channel blocker from the cone shell Conus geographus, is a 27 residue polypeptide cross-linked by three disulfide bonds. Here, we report the synthesis, structural analysis by (1)H NMR and bioassay of analogues of GVIA with disulfide bridge deletions and N- and C-terminal truncations. Two analogues that retain the crucial Lys-2 and Tyr-13 residues in loops constrained by two native disulfide bridges were synthesised using orthogonal protection of cysteine residues. In the first analogue, the Cys-15-Cys-26 disulfide bridge was deleted (by replacing the appropriate Cys residues with Ser), while in the second, this disulfide bridge and the eight C-terminal residues were deleted. No activity was detected for either analogue in a rat vas deferens assay, which measures N-type calcium channel activity in sympathetic nerve, and NMR studies showed that this was due to a gross loss of secondary and tertiary structure. Five inactive analogues that were synthesised without orthogonal protection of Cys residues as part of a previous study (Flinn et al. (1995) J. Pept. Sci. 1, 379-384) were also investigated. Three had single disulfide deletions (via Ser substitutions) and two had N- or C-terminal deletions in addition to the disulfide deletion. Peptide mapping and NMR analyses demonstrated that at least four of these analogues had non-native disulfide pairings, which presumably accounts for their lack of activity. The NMR studies also showed that all five analogues had substantially altered tertiary structures, although the backbone chemical shifts and nuclear Overhauser enhancements (NOEs) implied that native-like turn structures persisted in some of these analogues despite the non-native disulfide pairings. This work demonstrates the importance of the disulfides in omega-conotoxin folding and shows that the Cys-15-Cys-26 disulfide is essential for activity in GVIA. The NMR analyses also emphasise that backbone chemical shifts and short- and medium-range NOEs are dictated largely by local secondary structure elements and are not necessarily reliable monitors of the tertiary fold.     

The positions of the disulfide bonds and the glycosylation site in a lectin of the acorn barnacle Megabalanus rosa

The positions of the interchain and intrachain disulfide bonds and the glycosylation site in a lectin of the acorn barnacle Megabalanus rosa were determined. The lectin (Mr 140,000) is composed of the same subunit (Mr 22,000) which is cross-linked by disulfide bonds to form a dimer. Intact lectin yielded two fragments, CB1 and CB2, by cleavage with cyanogen bromide. One intrachain and two interchain disulfide bonds were identified as Cys-53-Cys-61, Cys-14-Cys-50' and Cys-50-Cys-14', respectively, by enzymatic digestion and Edman degradation of CB1. Two intrachain disulfide bonds were determined as Cys-78-Cys-168 and Cys-144-Cys-160 by enzymatic digestion of CB2. The two intrachain disulfide bonds are well conserved through all invertebrate lectins and calcium-dependent animal lectins. S-Carboxamidomethylated lectin was digested with Staphylococcus aureus V8 proteinase and separated by reversed-phase HPLC. Glycopeptides were detected by the 4-N,N-dimethylamino-4'-azobenzene sulfonyl hyrazide method. Sequence analyses of the glycopeptides showed that a carbohydrate chain attached to Asn-39.     

Disulfide structures of highly bridged peptides: a new strategy for analysis.

A new approach is described for analyzing disulfide linkage patterns in peptides containing tightly clustered cystines. Such peptides are very difficult to analyze with traditional strategies, which require that the peptide chain be split between close or adjacent Cys residues. The water-soluble tris-(2-carboxyethyl)-phosphine (TCEP) reduced disulfides at pH 3, and partially reduced peptides were purified by high performance liquid chromatography with minimal thiol-disulfide exchange. Alkylation of free thiols, followed by sequencer analysis, provided explicit assignment of disulfides that had been reduced. Thiol-disulfide exchange occurred during alkylation of some peptides, but correct deductions were still possible. Alkylation competed best with exchange when peptide solution was added with rapid mixing to 2.2 M iodoacetamide. Variants were developed in which up to three alkylating agents were used to label different pairs of thiols, allowing a full assignment in one sequencer analysis. Model peptides used included insulin (three bridges, intra- and interchain disulfides; -Cys.Cys- pair), endothelin and apamin (two disulfides; -Cys.x.Cys- pair), conotoxin GI and isomers (two disulfides; -Cys.Cys- pair), and bacterial enterotoxin (three bridges within 13 residues; two -Cys.Cys- pairs). With insulin, all intermediates in the reduction pathway were identified; with conotoxin GI, analysis was carried out successfully for all three disulfide isomers. In addition to these known structures, the method has been applied successfully to the analysis of several previously unsolved structures of similar complexity. Rates of reduction of disulfide bonds varied widely, but most peptides did not show a strongly preferred route for reduction.     

Refined crystal structure of troponin C from turkey skeletal muscle at 2.0 A resolution

The crystal structure of troponin C from turkey skeletal muscle has been refined at 2.0 A resolution (1 A = 0.1 nm). The resulting crystallographic R factor (R = sigma[[Fo[-[Fc[[/sigma[Fo[, where [Fo[ and [Fc[ are the observed and calculated structure factor amplitudes) is 0.155 for the 8054 reflections with intensities I greater than or equal to 2 sigma(I) within the 10 A to 2.0 A resolution range. With 66% of the residues in helical conformation, troponin C provides a good sample for helix analysis. The mean alpha-helix dihedral angles (phi, psi = -62 degrees, -42 degrees) agree with values observed for helical regions in other proteins. The helices are all curved and/or kinked. In particular, the 31 amino acid long inter-domain helix is smoothly curved, with a rather large radius of curvature of 137 A. Helix packing is different in the Ca2+-free domain (N-terminal) and the Ca2+-bound domain (C-terminal). The inter-helix angles for the two helix-loop-helix motifs in the regulatory domain are 133 degrees and 151 degrees, whereas the value for the two motifs in the C-terminal domain is 110 degrees, as observed in the EF-hands of parvalbumin. These differences affect the packing of the respective hydrophobic cores of each domain, in particular the disposition of aromatic rings. Pairwise arrangement of Ca2+-binding loops is common to both states, but the conformation is markedly different. Conversion of one to the other can be achieved by small cumulative changes of main-chain dihedral angles. The integrity of loop structure is maintained by numerous electrostatic interactions. Both salt bridges and carboxyl-carboxylate interactions are observed in TnC. There are more intramolecular (9) than intermolecular (1) salt bridges. Carboxyl-carboxylate interactions occur because the pH of the crystals is 5.0 and there is a multitude of aspartate and glutamate residues. One is intramolecular and four are intermolecular. Polar side-chain interactions occur more commonly with main-chain carbonyls and amides than with other polar side-chains. These interactions are mostly short range, and are similar to those observed in other proteins with one exception: negatively charged side-chains interact more frequently with main-chain carbonyl oxygen atoms. However, out of 19 such interactions, 10 involve oxygen atoms of the Ca2+ ligands. These unfavorable interactions are compensated by the favorable interactions with the Ca2+ ions and with main-chain amides. They are a trivial consequence of the tight fold of the Ca2+-binding loops.     

The position of the disulfide bonds in human plasma α2HS-glycoprotein and the repeating double disulfide bonds in the domain structure

The positions of the inter- and intra-chain disulfide bonds of human plasma α₂HS-glycoprotein were determined. α₂HS-glycoprotein was digested with acid proteinase and then with thermolysin. The disulfide bonds containing peptides were separated by reversed-phase HPLC and detected by SBD-F (7-fluorobenzo-2-oxa-1,3-diasole-4-sulfonic acid ammonium salt) method. One inter-disulfide bond containing peptide and five intra-disulfide bond containing peptides (A-chain) were purified and identified as Cys-18 (B-chain)-Cys-14 (A-chain), Cys-71-Cys-82, Cys-96-Cys-114, Cys-128-Cys-131, Cys-190-Cys-201 and Cys-212-Cys-229, respectively. The location of the intra-disulfide bonds revealed that the A-chain of α₂HS-glycoprotein is composed of three domains. Two domains were shown to possess intramolecular homology judging from the total chain length of the domains, size of the loops formed by the SS bonds, the location of two disulfide loops near the C-terminal end of domains A and B, the distance between two SS bonds of each domain, the amino acid sequence homology between these two domains (22.6%), number of amino acid residues between the second SS loops and the end of domains A and B, and the positions of the ordered structures.     

Stereochemical modeling of disulfide bridges. Criteria for introduction into proteins by site-directed mutagenesis

A computer modeling procedure for assessing the stereochemical suitability of pairs of residues in proteins as potential sites for introduction of cystine disulfide crosslinks has been developed. Residue pairs with C alpha-C alpha distances of less than or equal to 6.5 A and C beta-C beta distances of less than or equal to 4.5 A are chosen for geometrical fixation of S atoms using the program MODIP. The stereochemistry of the modeled disulfides is evaluated using limits for the structural parameters of the various torsion angles and S-S bond length in the disulfide bridge. The ability of the procedure to correctly model disulfides has been checked with examples of cystine peptides of known crystal structures and 103 disulfide bridges from 25 available protein crystal structures determined at less than or equal to 2 A resolution. An analysis of results on three proteins with engineered disulfides, T4 lysozyme, dihydrofolate reductase and subtilisin, is presented. Two positions for the introduction of 'stereochemically optimal' disulfides are identified in subtilisin.     

Function and solution structure of huwentoxin-IV,a potent neuronal tetrodotoxin (TTX)-sensitive sodium channel antagonist from Chinese bird spider Selenocosmia huwena

We have isolated a highly potent neurotoxin from the venom of the Chinese bird spider, Selenocosmia huwena. This 4.1-kDa toxin, which has been named huwentoxin-IV, contains 35 residues with three disulfide bridges: Cys-2-Cys-17, Cys-9-Cys-24, and Cys-16-Cys-31, assigned by a chemical strategy including partial reduction of the toxin and sequence analysis of the modified intermediates. It specifically inhibits the neuronal tetrodotoxin-sensitive (TTX-S) voltage-gated sodium channel with the IC(50) value of 30 nm in adult rat dorsal root ganglion neurons, while having no significant effect on the tetrodotoxin-resistant (TTX-R) voltage-gated sodium channel. This toxin seems to be a site I toxin affecting the sodium channel through a mechanism quite similar to that of TTX: it suppresses the peak sodium current without altering the activation or inactivation kinetics. The three-dimensional structure of huwentoxin-IV has been determined by two-dimensional (1)H NMR combined with distant geometry and simulated annealing calculation by using 527 nuclear Overhauser effect constraints and 14 dihedral constraints. The resulting structure is composed of a double-stranded antiparallel beta-sheet (Leu-22-Ser-25 and Trp-30-Tyr-33) and four turns (Glu-4-Lys-7, Pro-11-Asp-14, Lys-18-Lys-21 and Arg-26-Arg-29) and belongs to the inhibitor cystine knot structural family. After comparison with other toxins purified from the same species, we are convinced that the positively charged residues of loop IV (residues 25-29), especially residue Arg-26, must be crucial to its binding to the neuronal tetrodotoxin-sensitive voltage-gated sodium channel.     

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.     

Identification of in vivo disulfide conformation of TRPA1 ion channel

TRPA1 (transient receptor potential ankyrin 1) is an ion channel expressed in the termini of sensory neurons and is activated in response to a broad array of noxious exogenous and endogenous thiol-reactive compounds, making it a crucial player in chemical nociception. A number of conserved cysteine residues on the N-terminal domain of the channel have been identified as critical for sensing these electrophilic pungent chemicals, and our recent EM structure with modeled domains predicts that these cysteines form a ligand-binding pocket, allowing for the possibility of disulfide bonding between the cysteine residues. Here, we present a comprehensive mass spectrometry investigation of the in vivo disulfide bonding conformation and in vitro reactivity of 30 of the 31 cysteine residues in the TRPA1 ion channel. Four disulfide bonds were detected in the in vivo TRPA1 structure: Cys-666-Cys-622, Cys-666-Cys-463, Cys-622-Cys-609, and Cys-666-Cys-193. All of the cysteines detected were reactive to N-methylmaleimide (NMM) in vitro, with varying degrees of labeling efficiency. Comparison of the ratio of the labeling efficiency at 300 μM versus 2 mM NMM identified a number of cysteine residues that were outliers from the mean labeling ratio, suggesting that protein conformation changes rendered these cysteines either more or less protected from labeling at the higher NMM concentrations. These results indicate that the activation mechanism of TRPA1 may involve N-terminal conformation changes and disulfide bonding between critical cysteine residues.     

Relevance of cysteine residues for biosynthesis and antigenicity of human hepatitis B virus e protein.

The mature form of the secretory core protein (HBe protein) of human hepatitis B virus contains four cysteines which are located at amino acid positions -7, 48, 61, and 107 relative to the HBc start methionine. In addition, there is a cysteine, Cys-183, located in the C-terminal domain of the HBe precursor, which is cleaved during HBe maturation. Here, the significance of these cysteines for biosynthesis and antigenicity of the HBe protein was examined. The cysteines at positions -7 and 61 were found to be crucial for HBe biosynthesis. As has already been described, if the Cys at position -7 is mutated, disulfide-linked HBe homodimers which have both HBe antigenicity and HBc antigenicity are expressed. Here we show that these dimers are due to Cys-61-Cys-61 disulfide bridges which are formed only if the Cys at position -7 is not present. In the wild-type protein, this dimerization appears to be inhibited by formation of intramolecular disulfide bridges between the Cys at -7 and one of the internal cysteines. Moreover, Cys-61 is important for HBe biosynthesis in general since mutation of this amino acid results in production of HBe proteins which are either only poorly secreted or possess a different antigenicity.     

Different roles of the two disulfide bonds of the cysteine proteinase inhibitor, chicken cystatin, for the conformation of the active protein.

The Cys-71-Cys-81 disulfide bond of the cysteine proteinase inhibitor, chicken cystatin, was specifically reduced by thioredoxin or low concentrations of dithiothreitol. This cleavage, followed by S-carbamoylmethylation, induced a conformational change of the protein, as evidenced by changes in isoelectric point and circular dichroism spectra and by an increased susceptibility to digestion by nontarget proteinases. The proteinase binding ability and the immunological properties of the inhibitor, however, were not detectably altered, indicating that the conformational change was limited to the region around the disrupted bond. In contrast, reduction of both disulfide bonds of cystatin by higher concentrations of dithiothreitol and subsequent alkylation led to the slow conversion of the inhibitor into two forms lacking proteinase binding ability, indicative of more extensive conformational changes. Together, these results suggest that the less accessible Cys-95-Cys-115 disulfide bond of chicken cystatin, but not the more accessible Cys-71-Cys-81 bond, is of importance for maintaining the conformation of the inhibitor required for binding of target proteinases.     

Analysis of the relationship between side-chain conformation and secondary structure in globular proteins

The relationship between the preferred side-chain dihedral angles and the secondary structure of a residue was examined. The structures of 61 proteins solved to a resolution of 2.0 A (1 A = 0.1 nm) or better were analysed using a relational database to store the information. The strongest feature observed was that the chi 1 distribution for most side-chains in an alpha-helix showed an absence of the g- conformation and a shift towards the t conformation when compared to the non-alpha/beta structures. The exceptions to this tendency were for short polar side-chains that form hydrogen bonds with the main-chain which prefer g+. Shifts in the chi 1 preferences for residues in the beta-sheet were observed. Other side-chain dihedral angles (chi 2, chi 3, chi 4) were found to be influenced by the main-chain. This paper presents more accurate distributions for the side-chain dihedral angles which were obtained from the increased number of proteins determined to high resolution. The means and standard deviations for chi 1 and chi 2 angles are presented for all residues according to the secondary structure of the main-chain. The means and standard deviations are given for the most popular conformations for side-chains in which chi 3 and chi 4 rotations affect the position of C atoms.     

Carrier subunit of plasma membrane transporter is required for oxidative folding of its helper subunit

We study the amino acid transport system b(0,+) as a model for folding, assembly, and early traffic of membrane protein complexes. System b(0,+) is made of two disulfide-linked membrane subunits: the carrier, b(0,+) amino acid transporter (b(0,+)AT), a polytopic protein, and the helper, related to b(0,+) amino acid transporter (rBAT), a type II glycoprotein. rBAT ectodomain mutants display folding/trafficking defects that lead to type I cystinuria. Here we show that, in the presence of b(0,+)AT, three disulfides were formed in the rBAT ectodomain. Disulfides Cys-242-Cys-273 and Cys-571-Cys-666 were essential for biogenesis. Cys-673-Cys-685 was dispensable, but the single mutants C673S, and C685S showed compromised stability and trafficking. Cys-242-Cys-273 likely was the first disulfide to form, and unpaired Cys-242 or Cys-273 disrupted oxidative folding. Strikingly, unassembled rBAT was found as an ensemble of different redox species, mainly monomeric. The ensemble did not change upon inhibition of rBAT degradation. Overall, these results indicated a b(0,+)AT-dependent oxidative folding of the rBAT ectodomain, with the initial and probably cotranslational formation of Cys-242-Cys-273, followed by the oxidation of Cys-571-Cys-666 and Cys-673-Cys-685, that was completed posttranslationally.     

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.     

Primary structure of a hemorrhagic metalloproteinase, HT-2, isolated from the venom of Crotalus ruber ruber

Crotalidae and Viperidae snake venoms contains several kinds of metalloproteinases which cause localized hemorrhage by direct action on blood vessel walls. We report here the entire amino acid sequence and the disulfide bridge locations of HT-2, one of the hemorrhagic toxins isolated from the venom of Crotalus ruber ruber (red rattlesnake). The non-reduced protein was first cleaved at methionine residues to provide a set of 8 fragments, which covered the entire sequence of HT-2. The disulfide bridge locations of HT-2 were also determined by using these primary fragments. The unambiguous sequence for the whole protein was then established by conventional methods using lysyl endopeptidase and thermolysin digests. HT-2 consisted of 202 amino acid residues with two disulfide bridges, which were assigned to Cys-117-Cys-197 and Cys-157-Cys-164. HT-2 had a typical zinc-chelating sequence His-Glu-X-X-His (residues 142-146) found in thermolysin, and its overall sequence showed, respectively, 50, 52, and 53% identities to those of HR2a, H2-proteinase, and the metalloproteinase domain of HR1B. However, the disulfide bridge locations of HT-2 were different from those in the other metalloproteinases. The primary structure of HT-2 was more closely related to that of Ht-d from Crotalus atrox recently determined (81% identity). From the structural comparison with five metalloproteinases so far elucidated, six conservative amino acid residues, which may possibly be related to the induction of the hemorrhagic activity, were suggested to be present in these toxins.