Role of disulfide bonds in modulating internal motions of proteins to tune their function: molecular dynamics simulation of scorpion toxin Lqh III

A series of 1-ns MD simulations were performed on the scorpion toxin Lqh III in native and disulfide bond broken states. The removal of disulfide bonds has caused hydrogen bond network alteration in the five-residue turn, the long loop, the alpha-helix, the loop connecting strands II and III, and the C-terminal region. In addition and more importantly, it has influenced the amplitude of the fluctuations of five-residue turn, loops, and C-terminal region with a minor effect on the fluctuations of the cysteines in the broken bond sites. These findings suggest that disulfide bonds are not the most important factors in rigidifying their own locations, while they have more important effects at a global scale. Furthermore, our results reveal that disulfide bonds have considerable influence on the functionally important essential modes of motions and the correlations between the motions of the binding site residues. Therefore, we can conclude that disulfide bonds have a crucial role in modulating the function via adjusting the dynamics of scorpion toxin molecules. Although this conclusion cannot be generalized to all peptides and proteins, it demonstrates the importance of more investigations on this aspect of disulfide bond efficacy.     

Investigation of the role of disulphide bond in modulating internal motions of BmK AGAP protein by molecular dynamics simulation

Elucidating structural determinants in the functional regions of toxins can provide useful knowledge for designing novel analgesic peptides. A series of 100 ns MD simulations were performed on the scorpion toxin BmK AGAP in native and disulphide bond broken states. The comparison of disulphide bond broken states with the native state showed the α-helix was found to be the key to the analgesic activity. Furthermore, our results revealed disulphide bonds have considerable influence on the functionally important essential modes of motions and the correlations between the motions of the Core domain and the C-terminal region which are involved in the analgesic activity. Therefore, we can conclude that disulphide bonds have a crucial role in modulating the function via adjusting the dynamics of scorpion toxin BmK AGAP molecule.     

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.     

Restoration of structural stability and ligand binding after removal of the conserved disulfide bond in tear lipocalin

Disulfide bonds play diverse structural and functional roles in proteins. In tear lipocalin (TL), the conserved sole disulfide bond regulates stability and ligand binding. Probing protein structure often involves thiol selective labeling for which removal of the disulfide bonds may be necessary. Loss of the disulfide bond may destabilize the protein so strategies to retain the native state are needed. Several approaches were tested to regain the native conformational state in the disulfide-less protein. These included the addition of trimethylamine N-oxide (TMAO) and the substitution of the Cys residues of disulfide bond with residues that can either form a potential salt bridge or others that can create a hydrophobic interaction. TMAO stabilized the protein relaxed by removal of the disulfide bond. In the disulfide-less mutants of TL, 1.0 M TMAO increased the free energy change (ΔG⁰) significantly from 2.1 to 3.8 kcal/mol. Moderate recovery was observed for the ligand binding tested with NBD-cholesterol. Because the disulfide bond of TL is solvent exposed, the substitution of the disulfide bond with a potential salt bridge or hydrophobic interaction did not stabilize the protein. This approach should work for buried disulfide bonds. However, for proteins with solvent exposed disulfide bonds, the use of TMAO may be an excellent strategy to restore the native conformational states in disulfide-less analogs of the proteins.     

Catalysis of protein disulfide bond isomerization in a homogeneous substrate

Protein disulfide isomerase (PDI) catalyzes the rearrangement of nonnative disulfide bonds in the endoplasmic reticulum of eukaryotic cells, a process that often limits the rate at which polypeptide chains fold into a native protein conformation. The mechanism of the reaction catalyzed by PDI is unclear. In assays involving protein substrates, the reaction appears to involve the complete reduction of some or all of its nonnative disulfide bonds followed by oxidation of the resulting dithiols. The substrates in these assays are, however, heterogeneous, which complicates mechanistic analyses. Here, we report the first analysis of disulfide bond isomerization in a homogeneous substrate. Our substrate is based on tachyplesin I, a 17-mer peptide that folds into a beta hairpin stabilized by two disulfide bonds. We describe the chemical synthesis of a variant of tachyplesin I in which its two disulfide bonds are in a nonnative state and side chains near its N and C terminus contain a fluorescence donor (tryptophan) and acceptor (N(epsilon)-dansyllysine). Fluorescence resonance energy transfer from 280 to 465 nm increases by 28-fold upon isomerization of the disulfide bonds into their native state (which has a lower E(o') = -0.313 V than does PDI). We use this continuous assay to analyze catalysis by wild-type human PDI and a variant in which the C-terminal cysteine residue within each Cys-Gly-His-Cys active site is replaced with alanine. We find that wild-type PDI catalyzes the isomerization of the substrate with kcat/K(M) = 1.7 x 10(5) M(-1) s(-1), which is the largest value yet reported for catalysis of disulfide bond isomerization. The variant, which is a poor catalyst of disulfide bond reduction and dithiol oxidation, retains virtually all of the activity of wild-type PDI in catalysis of disulfide bond isomerization. Thus, the C-terminal cysteine residues play an insignificant role in the isomerization of the disulfide bonds in nonnative tachyplesin I. We conclude that catalysis of disulfide bond isomerization by PDI does not necessarily involve a cycle of substrate reduction/oxidation.     

Nonnative interactions between cysteines direct productive assembly of P22 tailspike protein

Nonnative disulfide bond formation can play a critical role in the assembly of disulfide bonded proteins. During the folding and assembly of the P22 tailspike protein, nonnative disulfide bonds form both in vivo and in vitro. However, the mechanism and identity of cysteine disulfide pairs remains elusive, particularly for P22 tailspike, which contains no disulfide bonds in its native, functional form. Understanding the interactions between cysteine residues is important for developing a mechanistic model for the role of nonnative cysteines in P22 tailspike assembly. Prior in vivo studies have suggested that cysteines 496, 613, and 635 are the most likely site for sulfhydryl reactivity. Here we demonstrate that these three cysteines are critical for efficient assembly of tailspike trimers, and that interactions between cysteine pairs lead to productive assembly of native tailspike.     

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.     

Prediction of the disulfide-bonding state of cysteine in proteins

The bonding states of cysteine play important functional and structural roles in proteins. In particular, disulfide bond formation is one of the most important factors influencing the three-dimensional fold of proteins. Proteins of known structure were used to teach computer-simulated neural networks rules for predicting the disulfide-bonding state of a cysteine given only its flanking amino acid sequence. Resulting networks make accurate predictions on sequences different from those used in training, suggesting that local sequence greatly influences cysteines in disulfide bond formation. The average prediction rate after seven independent network experiments is 81.4% for disulfide-bonded and 80.0% for non-disulfide-bonded scenarios. Predictive accuracy is related to the strength of network output activities. Network weights reveal interesting position-dependent amino acid preferences and provide a physical basis for understanding the correlation between the flanking sequence and a cysteine's disulfide-bonding state. Network predictions may be used to increase or decrease the stability of existing disulfide bonds or to aid the search for potential sites to introduce new disulfide bonds.     

The role of short-range Cys171-Cys178 disulfide bond in maintaining cutinase active site integrity: A molecular dynamics simulation

Understanding structural determinants in enzyme active site integrity can provide a good knowledge to design efficient novel catalytic machineries. Fusarium solani pisi cutinase with classic triad Ser-His-Asp is a promising enzyme to scrutinize these structural determinants. We performed two MD simulations: one, with the native structure, and the other with the broken Cys171-Cys178 disulfide bond. This disulfide bond stabilizes a turn in active site on which catalytic Asp175 is located. Functionally important H-bonds and atomic fluctuations in catalytic pocket have been changed. We proposed that this disulfide bond within active site can be considered as an important determinant of cutinase active site structural integrity.     

Atomistic details of effect of disulfide bond reduction on active site of Phytase B from Aspergillus niger: A MD Study

The molecular integrity of the active site of phytases from fungi is critical for maintaining phytase function as efficient catalytic machines. In this study, the molecular dynamics (MD) of two monomers of phytase B from Aspergillus niger, the disulfide intact monomer (NAP) and a monomer with broken disulfide bonds (RAP), were simulated to explore the conformational basis of the loss of catalytic activity when disulfide bonds are broken. The simulations indicated that the overall secondary and tertiary structures of the two monomers were nearly identical but differed in some crucial secondary–structural elements in the vicinity of the disulfide bonds and catalytic site. Disulfide bonds stabilize the β-sheet that contains residue Arg66 of the active site and destabilize the α-helix that contains the catalytic residue Asp319. This stabilization and destabilization lead to changes in the shape of the active–site pocket. Functionally important hydrogen bonds and atomic fluctuations in the catalytic pocket change during the RAP simulation. None of the disulfide bonds are in or near the catalytic pocket but are most likely essential for maintaining the native conformation of the catalytic site.

Abbreviations

PhyB - 2.5 pH acid phophatese from Aspergillus niger, NAP - disulphide intact monomer of Phytase B, RAP - disulphide reduced monomer of Phytase B, Rg - radius of gyration, RMSD - root mean square deviation, MD - molecular dynamics.     

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.     

The contribution of cross-links to protein stability: A normal mode analysis of the configurational entropy of the native state

The vibrational entropy of native BPTI, with three disulfide bonds, was determined by use of normal mode calculations and compared with that of folded variants having either one less disulfide bond or lacking a peptide bond at the trypsin-reactive site. Favorable contributions to the free energy of 2.5–5.1 kcal/mol at 300 K were calculated for the reduction of disulfide bonds in the folded state, whereas no favorable contribution was found for the hydrolysis of the peptide bond cleaved by trypsin. This is on the order of the effect of disulfides in the unfolded state. The implications of these results for the stabilization of a folded protein by the introduction of crosslinks are discussed. © 1993 Wiley-Liss, Inc.     

Expression and characterization of recombinant kurtoxin, an inhibitor of T-type voltage-gated calcium channels

Kurtoxin, a 63-amino acid peptide stabilized by four disulfide bonds, is the first reported peptide inhibitor of T-type voltage-gated calcium channels. Although T-type calcium channels have been implicated in a number of disease states, including epilepsy, chronic pain, hypertension and cancer, the lack of selective inhibitors has slowed progress in understanding their precise roles. Kurtoxin is a potentially valuable tool with which to study T-type calcium channels. However, because of the limited availability of the native protein, little is known about the structure and molecular mechanism of kurtoxin. Here we report the expression of kurtoxin in Escherichia coli and the structural and functional characterization of the recombinant protein. The disulfide bond pairings and secondary structure of recombinant kurtoxin were characterized through enzymatic cleavage, mass analysis and CD spectroscopy. Recombinant kurtoxin almost completely inhibited the T-type calcium channel in a manner identical to the native toxin. The availability of recombinant kurtoxin that is identical to the native toxin should help in the study of T-type calcium channels and enable development of new strategies for producing even more-selective T-type calcium channel inhibitors and for investigating the molecular basis of the toxin-channel interactions.     

Disulfide bond effects on protein stability: designed variants of Cucurbita maxima trypsin inhibitor-V

Attempts to increase protein stability by insertion of novel disulfide bonds have not always been successful. According to the two current models, cross-links enhance stability mainly through denatured state effects. We have investigated the effects of removal and addition of disulfide cross-links, protein flexibility in the vicinity of a cross-link, and disulfide loop size on the stability of Cucurbita maxima trypsin inhibitor-V (CMTI-V; 7 kD) by differential scanning calorimetry. CMTI-V offers the advantage of a large, flexible, and solvent-exposed loop not involved in extensive intra-molecular interactions. We have uncovered a negative correlation between retention time in hydrophobic column chromatography, a measure of protein hydrophobicity, and melting temperature (T(m)), an indicator of native state stabilization, for CMTI-V and its variants. In conjunction with the complete set of thermodynamic parameters of denaturation, this has led to the following deductions: (1) In the less stable, disulfide-removed C3S/C48S (Delta Delta G(d)(50 degrees C) = -4 kcal/mole; Delta T(m) = -22 degrees C), the native state is destabilized more than the denatured state; this also applies to the less-stable CMTI-V* (Delta Delta G(d)(50 degrees C) = -3 kcal/mole; Delta T(m) = -11 degrees C), in which the disulfide-containing loop is opened by specific hydrolysis of the Lys(44)-Asp(45) peptide bond; (2) In the less stable, disulfide-inserted E38C/W54C (Delta Delta G(d)(50 degrees C) = -1 kcal/mole; Delta T(m) = +2 degrees C), the denatured state is more stabilized than the native state; and (3) In the more stable, disulfide-engineered V42C/R52C (Delta Delta G(d)(50 degrees C) = +1 kcal/mole; Delta T(m) = +17 degrees C), the native state is more stabilized than the denatured state. These results show that a cross-link stabilizes both native and denatured states, and differential stabilization of the two states causes either loss or gain in protein stability. Removal of hydrogen bonds in the same flexible region of CMTI-V resulted in less destabilization despite larger changes in the enthalpy and entropy of denaturation. The effect of a cross-link on the denatured state of CMTI-V was estimated directly by means of a four-state thermodynamic cycle consisting of native and denatured states of CMTI-V and CMTI-V*. Overall, the results show that an enthalpy-entropy compensation accompanies disulfide bond effects and protein stabilization is profoundly modulated by altered hydrophobicity of both native and denatured states, altered flexibility near the cross-link, and residual structure in the denatured state.     

Local interactions and the role of the 6-120 disulfide bond in alpha-lactalbumin: implications for formation of the molten globule state

Molten globule states are partially folded states of proteins which are compact and contain a high degree of secondary structure but which lack many of the fixed tertiary interactions associated with the native state. A set of peptides has been prepared in order to probe the role of local interactions in the vicinity of the Cys(6)-Cys(120) disulfide bond in stabilizing the molten globule state of human alpha-lactalbumin. Peptides derived from the N-terminal and C-terminal regions of human alpha-lactalbumin have been analyzed using nuclear magnetic resonance, circular dichroism, fluorescence spectroscopy and sedimentation equilibrium experiments. A peptide corresponding to the first helical region in the native protein, residues 1-13, is only slightly helical in isolation. Extending the peptide to include residues 14-18 results in a modest increase in helicity. A peptide derived from the C-terminal 12 residues, residues 112-123, is predominantly unstructured. Crosslinking the N- and C-terminal peptides by the native disulfide bond results in almost no increase in structure and there is no evidence for any significant cooperative structure formation over the range of pH 2.2-11.7. These results demonstrate that there is very little enhancement of local structure due to the formation of the Cys(6)-Cys(120) disulfide bond. This is in striking contrast to peptides derived from the region of the Cys(28)-Cys(111) disulfide.     

Effects of cysteine to serine substitutions in the two inter-chain disulfide bonds of insulin

Using site-directed mutagenesis we deleted the two inter-chain disulfide bonds of insulin, separately or both, by substitution of the cysteine residues with serine. Deletion of A20-B19 or both of the two inter-chain disulfide bonds resulted in the complete loss of secretion of the mutant single-chain porcine insulin precursor (PIP) from Saccharomyces cerevisiae cells. Removal of the A7-B7 disulfide bond resulted in a large reduction of secretion, but we could obtain the mutant for analysis of its biological and some physico-chemical properties. The A7-B7 disulfide bond deleted insulin mutant retained only 0.1% receptor-binding activity compared with porcine insulin, and its in vivo biological potency measured by mouse convulsion assay was also very low. We also studied some physico-chemical properties of the mutant using circular dichroism, native polyacrylamide gel electrophoresis and reversed-phase HPLC, which revealed some structural changes of the mutant peptides compared to native insulin. The present study shows that the two inter-chain disulfide bonds are important for efficient in vivo folding/secretion of PIP from yeast, especially the A20-B19 disulfide bond, and that the A7-B7 disulfide bond is crucial for maintaining the native conformation and biological activity of insulin.     

Role of disulfide bridges in the activity and stability of a cold-active alpha-amylase

The cold-adapted alpha-amylase from Pseudoalteromonas haloplanktis unfolds reversibly and cooperatively according to a two-state mechanism at 30 degrees C and unfolds reversibly and sequentially with two transitions at temperatures below 12 degrees C. To examine the role of the four disulfide bridges in activity and conformational stability of the enzyme, the eight cysteine residues were reduced with beta-mercaptoethanol or chemically modified using iodoacetamide or iodoacetic acid. Matrix-assisted laser desorption-time of flight mass spectrometry analysis confirmed that all of the cysteines were modified. The iodoacetamide-modified enzyme reversibly folded/unfolded and retained approximately one-third of its activity. Removal of all disulfide bonds resulted in stabilization of the least stable region of the enzyme (including the active site), with a concomitant decrease in activity (increase in activation enthalpy). Disulfide bond removal had a greater impact on enzyme activity than on stability (particularly the active-site region). The functional role of the disulfide bridges appears to be to prevent the active site from developing ionic interactions. Overall, the study demonstrated that none of the four disulfide bonds are important in stabilizing the native structure of enzyme, and instead, they appear to promote a localized destabilization to preserve activity.     

The effects of disulfide bonds on the denatured state of barnase

The effects of engineered disulfide bonds on protein stability are poorly understood because they can influence the structure, dynamics, and energetics of both the native and denatured states. To explore the effects of two engineered disulfide bonds on the stability of barnase, we have conducted a combined molecular dynamics and NMR study of the denatured state of the two mutants. As expected, the disulfide bonds constrain the denatured state. However, specific extended beta-sheet structure can also be detected in one of the mutant proteins. This mutant is also more stable than would be predicted. Our study suggests a possible cause of the very high stability conferred by this disulfide bond: the wild-type denatured ensemble is stabilized by a nonnative hydrophobic cluster, which is constrained from occurring in the mutant due to the formation of secondary structure.     

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.     

Contribution of the 6-120 disulfide bond of alpha-lactalbumin to the stabilities of its native and molten globule states.

The unfolding and refolding of a derivative of alpha-lactalbumin, in which the disulfide bond between Cys6 and Cys120 is selectively reduced and S-carboxymethylated, are investigated by equilibrium and kinetic circular dichroism measurements. The native conformation of this derivative is known to be essentially identical to that of intact alpha-lactalbumin. The equilibrium unfolding of the derivative involves a stable intermediate, which is also similar to the molten globule state of the disulfide intact protein. The results of stopped-flow circular dichroism experiments show that the same intermediate is formed rapidly as a transient intermediate in kinetic refolding. The conformational stabilities for the native and intermediate states have been estimated and compared with the stabilities for the corresponding states of intact alpha-lactalbumin. The stabilization of the native state by the disulfide has been interpreted in terms of a decrease in chain entropy in the unfolded state and elimination of the strain imposed on the disulfide bond in the native state. The molten globule state is also stabilized by the disulfide bond, although the degree of stabilization of the molten globule state is smaller than of the native state. The results suggest that, in the molten globule state, some ordered structures are present within the loop moiety formed by the 6-120 disulfide.