HIGH TORSIONAL ENERGY DISULFIDES: RELATIONSHIP BETWEEN CROSS-STRAND DISULFIDES AND RIGHT-HANDED STAPLES

Intramolecular Repulsion Visible Through the Electrostatic Potential in Disulfides: Analysis via Varying Iso-density Envelopes

For a series of diethyl disulfide conformations, the nearly touching contours of the electrostatic potential plotted on iso-density molecular surfaces allow the assessment of intramolecular repulsion. The electrostatic potential is plotted on varying iso-density envelopes to find the nearly touching contours for which (a) the surface electrostatic potential does not show overlap between atoms or functional groups and (b) the typical features are visible (σ-hole, lone pair, hydrogen VS,max). When these nearly touching contours X are closer to the nuclei, the more electron density is excluded from the iso-density envelopes and the smaller are the volumes corresponding to these envelopes. Both the contours X and the corresponding volumes are found to correlate with relative conformational energy, reflecting the degree of intramolecular repulsion present in the various diethyl disulfides. Quantitative estimates of intramolecular repulsion can be made based on relationships between the nearly touching contour X vs relative energy and volume (of the nearly touching contour X) vs relative energy, obtained for series of diethyl disulfide conformers. These relations were used to analyze intramolecular repulsion in a set of disulfides broader than diethyl disulfide conformers. We have shown that the approach of varying electronic density contours can be used in the study of repulsive intramolecular interactions, hereby extending earlier work involving attractive intermolecular interactions.

Towards solvent regulated self-activation of N-terminal disulfide bond oxidoreductase-D

N-terminal disulfide bond oxidoreductase-D (nDsbD), an essential redox enzyme in Gram-negative bacteria, consists of a single disulfide bond (Cys₁₀₃-Cys₁₀₉) in its active site. The enzymatic functions are believed to be regulated by an electron transfer mediated redox switching of the disulfide bond, which is vital in controlling bacterial virulence factors. In light of the disulfide bond's inclination towards nucleophilic cleavage, it is also plausible that an internal nucleophile could second the existing electron transfer mechanism in nDsbD. Using QM/MM MD metadynamics simulations, we explore different possibilities of generating an internal nucleophile near the nDsbD active site, which could serve as a fail-over mechanism in cleaving the disulfide bond. The simulations show the formation of the internal nucleophile Tyr₄₂O^- (F ≈ 9 kcal mol^-1) and its stabilization through the solvent medium. The static gas-phase calculations show that Tyr₄₂O^- could be a potential nucleophile for cleaving the S-S bond. Most strikingly, it is also seen that Tyr₄₂O^- and Asp₆₈OH communicate with each other through a proton-hole like water wire (F ≈ 12 kcal mol^-1), thus modulating the nucleophile formation. Accordingly, we propose the role of a solvent in regulating the internal nucleophilic reactions and the subsequent self-activation of nDsbD. We believe that this could be deterministic while designing enzyme-targeted inhibitor compounds.

A redox-active switch in fructosamine-3-kinases expands the regulatory repertoire of the protein kinase superfamily

Aberrant regulation of metabolic kinases by altered redox homeostasis substantially contributes to aging and various diseases, such as diabetes. We found that the catalytic activity of a conserved family of fructosamine-3-kinases (FN3Ks), which are evolutionarily related to eukaryotic protein kinases, is regulated by redox-sensitive cysteine residues in the kinase domain. The crystal structure of the FN3K homolog from Arabidopsis thaliana revealed that it forms an unexpected strand-exchange dimer in which the ATP-binding P-loop and adjoining β strands are swapped between two chains in the dimer. This dimeric configuration is characterized by strained interchain disulfide bonds that stabilize the P-loop in an extended conformation. Mutational analysis and solution studies confirmed that the strained disulfides function as redox "switches" to reversibly regulate the activity and dimerization of FN3K. Human FN3K, which contains an equivalent P-loop Cys, was also redox sensitive, whereas ancestral bacterial FN3K homologs, which lack a P-loop Cys, were not. Furthermore, CRISPR-mediated knockout of FN3K in human liver cancer cells altered the abundance of redox metabolites, including an increase in glutathione. We propose that redox regulation evolved in FN3K homologs in response to changing cellular redox conditions. Our findings provide insights into the origin and evolution of redox regulation in the protein kinase superfamily and may open new avenues for targeting human FN3K in diabetic complications.

Prediction of disulfide dihedral angles using chemical shifts

Cystine residues result from the formation of disulfide bonds between pairs of cysteine residues. This cross linking of the backbone is essential for the structure and activity of peptides and proteins. The conformation of a cystine side chain can be described using five dihedral angles, χ1, χ2, χ3, χ2', and χ1', with cystines favouring certain combinations of these angles. 2D NMR spectroscopy is ideally suited for structure determination of disulfide-rich peptides, because of their small size and constrained nature. However, only limited information of the cystine side chain conformation can be determined by NMR spectroscopy, leading to ambiguity in the deduced 3D structures. Resolving accurate structures is important as disulfide-rich peptides have proven to be promising drug candidates in a number of fields, either as bioactive leads or scaffolds. Using a database of NMR chemical shifts combined with crystallographic structures, we have developed a method called DISH that uses support vector machines to predict the dihedral angles of cysteine side chains. It is able to successfully predict χ2 angles with 91% accuracy, and has improved performance over existing prediction methods for χ1 angles, with 87% accuracy. For 81% of cysteine residues, DISH successfully predicted both the χ1 and χ2 angles. By revisiting published solution structures of peptides determined using NMR spectroscopy, we assessed the impact of additional cystine dihedral restraints on the quality of 3D models. DISH improved the resolution and accuracy, highlighting the potential for improving the understanding of structure-activity relationships and rational development of peptide drugs.

Cross-strand disulfides in the non-hydrogen bonding site of antiparallel β-sheet (aCSDns): poised for biological switching

aCSDns are forbidden disulfides with protein redox-activity. Within the aCSDn structural motif, a cognate substrate of Trx-like enzymes, the disulfide bonds are strained and metastable, facilitating their role as redox-regulated protein switches.

Post-translational control of protein function by disulfide bond cleavage

Protein action in nature is largely controlled by the level of expression and by post-translational modifications. Post-translational modifications result in a proteome that is at least two orders of magnitude more diverse than the genome. There are three basic types of post-translational modifications: covalent modification of an amino acid side chain, hydrolytic cleavage or isomerization of a peptide bond, and reductive cleavage of a disulfide bond. This review addresses the modification of disulfide bonds. Protein disulfide bonds perform either a structural or a functional role, and there are two types of functional disulfide: the catalytic and allosteric bonds. The allosteric disulfide bonds control the function of the mature protein in which they reside by triggering a change when they are cleaved. The change can be in ligand binding, substrate hydrolysis, proteolysis, or oligomer formation. The allosteric disulfides are cleaved by oxidoreductases or by thiol/disulfide exchange, and the configurations of the disulfides and the secondary structures that they link share some recurring features. How these bonds are being identified using bioinformatics and experimental screens and what the future holds for this field of research are also discussed.

Conformational analysis and design of cross-strand disulfides in antiparallel β-sheets

Cross-strand disulfides bridge two cysteines in a registered pair of antiparallel β-strands. A nonredundant data set comprising 5025 polypeptides containing 2311 disulfides was used to study cross-strand disulfides. Seventy-six cross-strand disulfides were found of which 75 and 1 occurred at non-hydrogen-bonded (NHB) and hydrogen-bonded (HB) registered pairs, respectively. Conformational analysis and modeling studies demonstrated that disulfide formation at HB pairs necessarily requires an extremely rare and positive χ¹ value for at least one of the cysteine residues. Disulfides at HB positions also have more unfavorable steric repulsion with the main chain. Thirteen pairs of disulfides were introduced in NHB and HB pairs in four model proteins: leucine binding protein (LBP), leucine, isoleucine, valine binding protein (LIVBP), maltose binding protein (MBP), and Top7. All mutants LIVBP T247C V331C showed disulfide formation either on purification, or on treatment with oxidants. Protein stability in both oxidized and reduced states of all mutants was measured. Relative to wild type, LBP and MBP mutants were destabilized with respect to chemical denaturation, although the sole exposed NHB LBP mutant showed an increase of 3.1°C in T(m). All Top7 mutants were characterized for stability through guanidinium thiocyanate chemical denaturation. Both exposed and two of the three buried NHB mutants were appreciably stabilized. All four HB Top7 mutants were destabilized (ΔΔG⁰ = -3.3 to -6.7 kcal/mol). The data demonstrate that introduction of cross-strand disulfides at exposed NHB pairs is a robust method of improving protein stability. All four exposed Top7 disulfide mutants showed mild redox activity.

The PAPS-independent aryl sulfotransferase and the alternative disulfide bond formation system in pathogenic bacteria

Sulfurylation of biomolecules (often termed sulfonation or sulfation) has been described in many organisms in all kingdoms of life. To date, most studies on sulfotransferases, the enzymes catalyzing sulfurylation, have focused on 3'-phosphate-5'-phosphosulfate (PAPS)-dependent enzymes, which transfer the sulfuryl group from this activated anhydride to hydroxyl groups of acceptor molecules. By contrast, the PAPS-independent aryl sulfotransferases (ASSTs) from bacteria, which catalyze sulfotransfer from phenolic sulfate esters to another phenol in the bacterial periplasm, were not well characterized until recently, although they were first described in 1986 in a search for nonhepatic sulfurylation processes. Recent studies revealed that this unusual class of sulfotransferases differs profoundly in both molecular structure and catalytic mechanism from PAPS-dependent sulfotransferases, and that ASSTs from certain bacterial pathogens are upregulated during infection. In this review, we summarize the literature on the roles of sulfurylation in prokaryotes and analyze the occurrence of ASSTs and their dependence on Dsb proteins catalyzing oxidative folding in the periplasm. Furthermore, we discuss structural differences and similarities between aryl sulfotransferases and PAPS-dependent sulfotransferases.

Disulfide conformation and design at helix N-termini

To understand structural and thermodynamic features of disulfides within an α‐helix, a non‐redundant dataset comprising of 5025 polypeptide chains containing 2311 disulfides was examined. Thirty‐five examples were found of intrahelical disulfides involving a CXXC motif between the N‐Cap and third helical positions. GLY and PRO were the most common amino acids at positions 1 and 2, respectively. The N‐Cap residue for disulfide bonded CXXC motifs had average (ϕ,ψ) values of (−112 ± 25.2°, 106 ± 25.4°). To further explore conformational requirements for intrahelical disulfides, CYS pairs were introduced at positions N‐Cap‐3; 1,4; 7,10 in two helices of an Escherichia coli thioredoxin mutant lacking its active site disulfide (nSS Trx). In both helices, disulfides formed spontaneously during purification only at positions N‐Cap‐3. Mutant stabilities were characterized by chemical denaturation studies (in both oxidized and reduced states) and differential scanning calorimetry (oxidized state only). All oxidized as well as reduced mutants were destabilized relative to nSS Trx. All mutants were redox active, but showed decreased activity relative to wild‐type thioredoxin. Such engineered disulfides can be used to probe helix start sites in proteins of unknown structure and to introduce redox activity into proteins. Conversely, a protein with CYS residues at positions N‐Cap and 3 of an α‐helix is likely to have redox activity. Proteins 2010. © 2009 Wiley‐Liss, Inc.

Estimating Relative Disulfide Energies: An Accurate Ab Initio Potential Energy Surface

Disulfide torsional energy, a good predictor of disulfide redox potential in proteins, may be estimated by interpolation on a potential energy surface (PES) describing the twisting of diethyl disulfide through its three central dihedral angles. Here we update PES calculations at the M05-2X level of theory with the 6-31G(d) basis set. Although the surface shows no qualitative differences from an earlier MP2(full) PES, energy differences greater than 1 kJ mol–1 were seen for conformations with χ2 between –60° and 30°, or with χ3 below 60° or above 130°. This is particularly significant for highly strained disulfides that are likely to be spontaneously reduced by mechanical means. In benchmarking against the high-level G3X method, M05-2X showed significantly reduced root mean squared deviation compared with MP2(full) (1.0 versus 2.0 kJ mol–1 respectively). Results are incorporated into a web application that calculates relative torsional energies from disulfide dihedral angles (http://www.sbinf.org/applications/pes.html).

Disulfides as redox switches: from molecular mechanisms to functional significance

The molecular mechanisms underlying thiol-based redox control are poorly defined. Disulfide bonds between Cys residues are commonly thought to confer extra rigidity and stability to their resident protein, forming a type of proteinaceous spot weld. Redox biologists have been redefining the role of disulfides over the last 30-40 years. Disulfides are now known to form in the cytosol under conditions of oxidative stress. Isomerization of extracellular disulfides is also emerging as an important regulator of protein function. The current paradigm is that the disulfide proteome consists of two subproteomes: a structural group and a redox-sensitive group. The redox-sensitive group is less stable and often associated with regions of stress in protein structures. Some characterized redox-active disulfides are the helical CXXC motif, often associated with thioredoxin-fold proteins; and forbidden disulfides, a group of metastable disulfides that disobey elucidated rules of protein stereochemistry. Here we discuss the role of redox-active disulfides as switches in proteins.

Evaluating the stability of disulfide bridges in proteins: a torsional potential energy surface for diethyl disulfide

Disulfide bonds formed by the oxidation of cysteine residues in proteins are the major form of intra- and inter-molecular covalent linkages in the polypeptide chain. To better understand the conformational energetics of this linkage, we have used the MP2(full)/6-31G(d) method to generate a full potential energy surface (PES) for the torsion of the model compound diethyl disulfide (DEDS) around its three critical dihedral angles (χ₂, χ₃, χ₂'). The use of ten degree increments for each of the parameters resulted in a continuous, fine-grained surface. This allowed us to accurately predict the relative stabilities of disulfide bonds in high resolution structures from the Protein Data Bank. The MP2(full) surface showed significant qualitative differences from the PES calculated using the Amber force field. In particular, a different ordering was seen for the relative energies of the local minima. Thus, Amber energies are not reliable for comparison of the relative stabilities of disulfide bonds. Surprisingly, the surface did not show a minimum associated with χ₂ ∼ - 60°, χ₃ ∼ 90, χ₂' ∼ - 60°. This is due to steric interference between Hα atoms. Despite this, significant populations of disulfides were found to adopt this conformation. In most cases this conformation is associated with an unusual secondary structure motif, the cross-strand disulfide. The relative instability of cross-strand disulfides is of great interest, as they have the potential to act as functional switches in redox processes.