The glutaredoxin -C-P-Y-C- motif: influence of peripheral residues

Molecular Fates of Organometallic Mercury in Human Brain

Mercury is ubiquitous in the environment, with rising levels due to pollution and climate change being a current global concern. Many mercury compounds are notorious for their toxicity, with the potential of organometallic mercury compounds for devastating effects on the structures and functions of the central nervous system being of particular concern. Chronic exposure of human populations to low levels of methylmercury compounds occurs through consumption of fish and other seafood, although the health consequences, if any, from this exposure remain controversial. We have used high energy resolution fluorescence detected X-ray absorption spectroscopy to determine the speciation of mercury and selenium in human brain tissue. We show that the molecular fate of mercury differs dramatically between individuals who suffered acute organometallic mercury exposure (poisoning) and individuals with chronic low-level exposure from a diet rich in marine fish. For long-term low-level methylmercury exposure from fish consumption, mercury speciation in brain tissue shows methylmercury coordinated to an aliphatic thiolate, resembling the coordination environment observed in marine fish. In marked contrast, for short-term high-level exposure, we observe the presence of biologically less available mercuric selenide deposits, confirmed by X-ray fluorescence imaging, as well as mercury(II)-bis-thiolate complexes, which may be signatures of severe poisoning in humans. These differences between low-level and high-level exposures challenge the relevance of studies involving acute exposure as a proxy for low-level chronic exposure.

Divergence of active site motifs among different classes of Populus glutaredoxins results in substrate switches

Enzymes are essential components of all biological systems. The key characteristics of proteins functioning as enzymes are their substrate specificities and catalytic efficiencies. In plants, most genes encoding enzymes are members of large gene families. Within such families, the contributions of active site motifs to the functional divergence of duplicate genes have not been well elucidated. In this study, we identified 41 glutaredoxin (GRX) genes in the Populus trichocarpa genome. GRXs are ubiquitous enzymes in plants that play important roles in developmental and stress tolerance processes. In poplar, GRX genes were divided into four classes based on clear differences in gene structure and expression pattern, subcellular localization, enzymatic activity, and substrate specificity of the encoded proteins. Using site-directed mutagenesis, this study revealed that the divergence of the active site motif among different classes of GRX proteins resulted in substrate switches and thus provided new insights into the molecular evolution of these important plant enzymes.

A Combined Proteomic and Metabolomic Strategy for Allergens Characterization in Natural and Fermented Brassica napus Bee Pollen

Bee pollen is consumed for its nutritional and pharmacological benefits, but it also contains hazardous allergens which have not been identified. Here, we identified two potential allergens, glutaredoxin and oleosin-B2, in Brassica napus bee pollen using mass spectrometry-based proteomics analyses, and used bioinformatics to predict their antigenic epitopes. Comparison of fermented (by Saccharomyces cerevisiae) and unfermented bee pollen samples indicated that glutaredoxin and oleosin-B2 contents were significantly decreased following fermentation, while the contents of their major constituent oligopeptides and amino acids were significantly increased based on metabolomics analyses. Immunoblot analysis indicated that the IgE-binding affinity with extracted bee pollen proteins was also significantly decreased after fermentation, suggesting a reduction in the allergenicity of fermented bee pollen. Furthermore, fermentation apparently promoted the biosynthesis of L-valine, L-isoleucine, L-tryptophan, and L-phenylalanine, as well as their precursors or intermediates. Thus, fermentation could potentially alleviate allergenicity, while also positively affecting nutritional properties of B. napus bee pollen. Our findings might provide a scientific foundation for improving the safety of bee pollen products to facilitate its wider application.

Activation of 3-Mercaptopyruvate Sulfurtransferase by Glutaredoxin Reducing System

Glutaredoxin (EC 1.15-1.21) is known as an oxidoreductase that protects cysteine residues within proteins against oxidative stress. Glutaredoxin catalyzes an electron transfer reaction that donates an electron to substrate proteins in the reducing system composed of glutaredoxin, glutathione, glutathione reductase, and nicotinamide-adenine dinucleotide phosphate (reduced form). 3-mercaptopyruvate sulfurtransferase (EC 2.8.1.2) is a cysteine enzyme that catalyzes transsulfuration, and glutaredoxin activates 3-mercaptopyruvate sulfurtransferase in the reducing system. Interestingly, even when glutathione or glutathione reductase was absent, 3-mercaptopyruvate sulfurtransferase activity increased, probably because reduced glutaredoxin was partly present and able to activate 3-mercaptopyruvate sulfurtransferase until depletion. A study using mutant Escherichia coli glutaredoxin1 (Cys¹⁴ is the binding site of glutathione and was replaced with a Ser residue) confirmed these results. Some inconsistency was noted, and glutaredoxin with higher redox potential than either 3-mercaptopyruvate sulfurtransferase or glutathione reduced 3-mercaptopyruvate sulfurtransferase. However, electron-transfer enzymatically proceeded from glutaredoxin to 3-mercaptopyruvate sulfurtransferase.

Thioredoxin-related protein of 14 kDa as a modulator of redox signalling pathways

Thioredoxin-related protein of 14 kDa (TRP14; also named TXNDC17 for thioredoxin domain-containing protein 17) is a highly conserved and ubiquitously expressed oxidoreductase. It is expressed in parallel with thioredoxin 1 (Trx1, TXN; TXN1), an efficient substrate for the mammalian cytosolic selenoprotein thioredoxin reductase 1 (TrxR1; TXNRD1). However, TRP14, in sharp contrast to Trx1, cannot support the activities of ribonucleotide reductase, peroxiredoxins or methionine sulfoxide reductases, thus is unable to directly support cell proliferation or antioxidant defence through these pathways. However, TRP14 has been shown to efficiently reduce l-cystine, which thereby indirectly supports glutathione synthesis. TRP14 can also suppress NF-κB signalling, is functionally linked to STAT3 signalling, and can directly reactivate oxidized protein-tyrosine phosphatase PTP1B. Furthermore, TRP14 can efficiently reduce persulfidated or nitrosylated cysteine residues in many proteins, thereby having the capacity to modulate signalling through hydrogen sulfide or NO. Additional bioinformatics analyses and observations suggest further roles for TRP14; therefore, further studies of its functions are warranted. Collectively, the results available suggest that TRP14 is a member of the thioredoxin system dedicated to the control of cellular redox signalling pathways. LINKED ARTICLES: This article is part of a themed section on Chemical Biology of Reactive Sulfur Species. To view the other articles in this section visit http://onlinelibrary.wiley.com/doi/10.1111/bph.v176.4/issuetoc.

Proteometabolomic response of Deinococcus radiodurans exposed to UVC and vacuum conditions: Initial studies prior to the Tanpopo space mission

The multiple extremes resistant bacterium Deinococcus radiodurans is able to withstand harsh conditions of simulated outer space environment. The Tanpopo orbital mission performs a long-term space exposure of D. radiodurans aiming to investigate the possibility of interplanetary transfer of life. The revealing of molecular machinery responsible for survivability of D. radiodurans in the outer space environment can improve our understanding of underlying stress response mechanisms. In this paper, we have evaluated the molecular response of D. radiodurans after the exposure to space-related conditions of UVC irradiation and vacuum. Notably, scanning electron microscopy investigations showed that neither morphology nor cellular integrity of irradiated cells was affected, while integrated proteomic and metabolomic analysis revealed numerous molecular alterations in metabolic and stress response pathways. Several molecular key mechanisms of D. radiodurans, including the tricarboxylic acid cycle, the DNA damage response systems, ROS scavenging systems and transcriptional regulators responded in order to cope with the stressful situation caused by UVC irradiation under vacuum conditions. These results reveal the effectiveness of the integrative proteometabolomic approach as a tool in molecular analysis of microbial stress response caused by space-related factors.

The Functions of Metamorphic Metallothioneins in Zinc and Copper Metabolism

Recent discoveries in zinc biology provide a new platform for discussing the primary physiological functions of mammalian metallothioneins (MTs) and their exquisite zinc-dependent regulation. It is now understood that the control of cellular zinc homeostasis includes buffering of Zn²⁺ ions at picomolar concentrations, extensive subcellular re-distribution of Zn²⁺, the loading of exocytotic vesicles with zinc species, and the control of Zn²⁺ ion signalling. In parallel, characteristic features of human MTs became known: their graded affinities for Zn²⁺ and the redox activity of their thiolate coordination environments. Unlike the single species that structural models of mammalian MTs describe with a set of seven divalent or eight to twelve monovalent metal ions, MTs are metamorphic. In vivo, they exist as many species differing in redox state and load with different metal ions. The functions of mammalian MTs should no longer be considered elusive or enigmatic because it is now evident that the reactivity and coordination dynamics of MTs with Zn²⁺ and Cu⁺ match the biological requirements for controlling—binding and delivering—these cellular metal ions, thus completing a 60-year search for their functions. MT represents a unique biological principle for buffering the most competitive essential metal ions Zn²⁺ and Cu⁺. How this knowledge translates to the function of other families of MTs awaits further insights into the specifics of how their properties relate to zinc and copper metabolism in other organisms.

Mechanistic insights on the reduction of glutathione disulfide by protein disulfide isomerase

We explore the enzymatic mechanism of the reduction of glutathione disulfide (GSSG) by the reduced a domain of human protein disulfide isomerase (hPDI) with atomistic resolution. We use classical molecular dynamics and hybrid quantum mechanics/molecular mechanics calculations at the mPW1N/6-311+G(2d,2p):FF99SB//mPW1N/6-31G(d):FF99SB level. The reaction proceeds in two stages: (i) a thiol-disulfide exchange through nucleophilic attack of the Cys53-thiolate to the GSSG-disulfide followed by the deprotonation of Cys56-thiol by Glu47-carboxylate and (ii) a second thiol-disulfide exchange between the Cys56-thiolate and the mixed disulfide intermediate formed in the first step. The Gibbs activation energy for the first stage was 18.7 kcal·mol^-1, and for the second stage, it was 7.2 kcal·mol^-1, in excellent agreement with the experimental barrier (17.6 kcal·mol^-1). Our results also suggest that the catalysis by protein disulfide isomerase (PDI) and thiol-disulfide exchange is mostly enthalpy-driven (entropy changes below 2 kcal·mol^-1 at all stages of the reaction). Hydrogen bonds formed between the backbone of His55 and Cys56 and the Cys56-thiol result in an increase in the Gibbs energy barrier of the first thiol-disulfide exchange. The solvent plays a key role in stabilizing the leaving glutathione thiolate formed. This role is not exclusively electrostatic, because an explicit inclusion of several water molecules at the density-functional theory level is a requisite to form the mixed disulfide intermediate. In the intramolecular oxidation of PDI, a transition state is only observed if hydrogen bond donors are nearby the mixed disulfide intermediate, which emphasizes that the thermochemistry of thiol-disulfide exchange in PDI is influenced by the presence of hydrogen bond donors.

Metal-coupled folding as the driving force for the extreme stability of Rad50 zinc hook dimer assembly

The binding of metal ions at the interface of protein complexes presents a unique and poorly understood mechanism of molecular assembly. A remarkable example is the Rad50 zinc hook domain, which is highly conserved and facilitates the Zn²⁺-mediated homodimerization of Rad50 proteins. Here, we present a detailed analysis of the structural and thermodynamic effects governing the formation and stability (logK₁₂ = 20.74) of this evolutionarily conserved protein assembly. We have dissected the determinants of the stability contributed by the small β-hairpin of the domain surrounding the zinc binding motif and the coiled-coiled regions using peptides of various lengths from 4 to 45 amino acid residues, alanine substitutions and peptide bond-to-ester perturbations. In the studied series of peptides, an >650 000-fold increase of the formation constant of the dimeric complex arises from favorable enthalpy because of the increased acidity of the cysteine thiols in metal-free form and the structural properties of the dimer. The dependence of the enthalpy on the domain fragment length is partially compensated by the entropic penalty of domain folding, indicating enthalpy-entropy compensation. This study facilitates understanding of the metal-mediated protein-protein interactions in which the metal ion is critical for the tight association of protein subunits.

Regulation of neovascularization by S-glutathionylation via the Wnt5a/sFlt-1 pathway

S-glutathionylation occurs when reactive oxygen or nitrogen species react with protein-cysteine thiols. Glutaredoxin-1 (Glrx) is a cytosolic enzyme which enzymatically catalyses the reduction in S-glutathionylation, conferring reversible signalling function to proteins with redox-sensitive thiols. Glrx can regulate vascular hypertrophy and inflammation by regulating the activity of nuclear factor κB (NF-κB) and actin polymerization. Vascular endothelial growth factor (VEGF)-induced endothelial cell (EC) migration is inhibited by Glrx overexpression. In mice overexpressing Glrx, blood flow recovery, exercise function and capillary density were significantly attenuated after hindlimb ischaemia (HLI). Wnt5a and soluble Fms-like tyrosine kinase-1 (sFlt-1) were enhanced in the ischaemic-limb muscle and plasma respectively from Glrx transgenic (TG) mice. A Wnt5a/sFlt-1 pathway had been described in myeloid cells controlling retinal blood vessel development. Interestingly, a Wnt5a/sFlt-1 pathway was found also to play a role in EC to inhibit network formation. S-glutathionylation of NF-κB components inhibits its activation. Up-regulated Glrx stimulated the Wnt5a/sFlt-1 pathway through enhancing NF-κB signalling. These studies show a novel role for Glrx in post-ischaemic neovascularization, which could define a potential target for therapy of impaired angiogenesis in pathological conditions including diabetes.

Protein flexibility and cysteine reactivity: influence of mobility on the H-bond network and effects on pKa prediction

Thanks to its chemical plasticity, cysteine (Cys) is a very versatile player in proteins. A major determinant of Cys reactivity is pKa: the ability to predict it is deemed critical in redox bioinformatics. I considered different computational methods for pKa predictions and ultimately applied one (propka, ppka1) to various datasets; for all residues I assessed the effect of (1) hydrogen bonding, electrostatics and solvation on predictions and (2) protein mobility on pKa variability. Particularly for Cys, exposure and H-bond contributions heavily dictated propka predictions. The prominence of H-bond contributions was previously reported: this may explain the effectiveness of ppka1 (with Cys, tested in a benchmark). However ppka1 was also very sensitive to protein mobility; I assessed the effects of mobility on particularly large (compared to previous studies) datasets of structural ensembles; I found that exposed Cys presented the highest pKa variability, ascribable to correspondingly high H-bond fluctuations associated with protein flexibility. The benefit of including protein dynamics in pKa predictions was previously proposed, but empirical methods were never tested in this sense; instead, giving their outstanding speed, they could lend particularly well to this purpose. I devised a strategy combining short range molecular dynamics with ppka1; the protocol aimed to mitigate high ppka1 variability by including a "statistical view" of fast conformational changes. Tested in a benchmark, the strategy lead to improved performances. These results provide new insights on Cys bioinformatics (pKa prediction protocols) and Cys biology (effect of mobility on exposed Cys properties).

Robustness analysis on interspecies interaction network for iron and glucose competition between Candida albicans and zebrafish during infection

Candida albicans has emerged as an important model organism for the study of infectious disease. Using high-throughput simultaneously quantified time-course transcriptomics, this study constructed host-pathogen interspecies interaction networks between C. albicans and zebrafish during the adhesion, invasion, and damage stages. Given that iron and glucose have been identified as crucial resources required during the infection process between C. albicans and zebrafish, we focused on the construction of the interspecies networks associated with them. Furthermore, a randomization technique was proposed to identify differentially regulated proteins that are statistically eminent for the three infection stages. The behaviors of the highly connected or differentially regulated proteins identified from the resulting networks were further investigated.

"Robustness" is an important system property that measures the ability of the system tolerating the intrinsic perturbations in a dynamic network. This characteristic provides a systematic and quantitative view to elucidate the dynamics of iron and glucose competition in terms of the interspecies interaction networks. Here, we further estimated the robustness of our constructed interspecies interaction networks for the three infection stages.

The constructed networks and robustness analysis provided significant insight into dynamic interactions related to iron and glucose competition during infection and enabled us to quantify the system's intrinsic perturbation tolerance ability during iron and glucose competition throughout the three infection stages. Moreover, the networks also assist in elucidating the offensive and defensive mechanisms of C. albicans and zebrafish during their competition for iron and glucose. Our proposed method can be easily extended to identify other such networks involved in the competition for essential resources during infection.

Oxidant sensing by reversible disulfide bond formation

Maintenance of the cellular redox balance is crucial for cell survival. An increase in reactive oxygen, nitrogen, or chlorine species can lead to oxidative stress conditions, potentially damaging DNA, lipids, and proteins. Proteins are very sensitive to oxidative modifications, particularly methionine and cysteine residues. The reversibility of some of these oxidative protein modifications makes them ideally suited to take on regulatory roles in protein function. This is especially true for disulfide bond formation, which has the potential to mediate extensive yet fully reversible structural and functional changes, rapidly adjusting the protein's activity to the prevailing oxidant levels.

Top-down proteomics reveals a unique protein S-thiolation switch in Salmonella Typhimurium in response to infection-like conditions

Characterization of the mature protein complement in cells is crucial for a better understanding of cellular processes on a systems-wide scale. Toward this end, we used single-dimension ultra-high-pressure liquid chromatography mass spectrometry to investigate the comprehensive "intact" proteome of the Gram-negative bacterial pathogen Salmonella Typhimurium. Top-down proteomics analysis revealed 563 unique proteins including 1,665 proteoforms generated by posttranslational modifications (PTMs), representing the largest microbial top-down dataset reported to date. We confirmed many previously recognized aspects of Salmonella biology and bacterial PTMs, and our analysis also revealed several additional biological insights. Of particular interest was differential utilization of the protein S-thiolation forms S-glutathionylation and S-cysteinylation in response to infection-like conditions versus basal conditions. This finding of a S-glutathionylation-to-S-cysteinylation switch in a condition-specific manner was corroborated by bottom-up proteomics data and further by changes in corresponding biosynthetic pathways under infection-like conditions and during actual infection of host cells. This differential utilization highlights underlying metabolic mechanisms that modulate changes in cellular signaling, and represents a report of S-cysteinylation in Gram-negative bacteria. Additionally, the functional relevance of these PTMs was supported by protein structure and gene deletion analyses. The demonstrated utility of our simple proteome-wide intact protein level measurement strategy for gaining biological insight should promote broader adoption and applications of top-down proteomics approaches.

pH-sensitive residues in the p19 RNA silencing suppressor protein from carnation Italian ringspot virus affect siRNA binding stability

Tombusviruses, such as Carnation Italian ringspot virus (CIRV), encode a protein homodimer called p19 that is capable of suppressing RNA silencing in their infected hosts by binding to and sequestering short-interfering RNA (siRNA) away from the RNA silencing pathway. P19 binding stability has been shown to be sensitive to changes in pH but the specific amino acid residues involved have remained unclear. Using constant pH molecular dynamics simulations, we have identified key pH-dependent residues that affect CIRV p19-siRNA binding stability at various pH ranges based on calculated changes in the free energy contribution from each titratable residue. At high pH, the deprotonation of Lys60, Lys67, Lys71, and Cys134 has the largest effect on the binding stability. Similarly, deprotonation of several acidic residues (Asp9, Glu12, Asp20, Glu35, and/or Glu41) at low pH results in a decrease in binding stability. At neutral pH, residues Glu17 and His132 provide a small increase in the binding stability and we find that the optimal pH range for siRNA binding is between 7.0 and 10.0. Overall, our findings further inform recent experiments and are in excellent agreement with data on the pH-dependent binding profile.

Electrostatics of cysteine residues in proteins: parameterization and validation of a simple model

One of the most popular and simple models for the calculation of pK(a) s from a protein structure is the semi-macroscopic electrostatic model MEAD. This model requires empirical parameters for each residue to calculate pK(a) s. Analysis of current, widely used empirical parameters for cysteine residues showed that they did not reproduce expected cysteine pK(a) s; thus, we set out to identify parameters consistent with the CHARMM27 force field that capture both the behavior of typical cysteines in proteins and the behavior of cysteines which have perturbed pK(a) s. The new parameters were validated in three ways: (1) calculation across a large set of typical cysteines in proteins (where the calculations are expected to reproduce expected ensemble behavior); (2) calculation across a set of perturbed cysteines in proteins (where the calculations are expected to reproduce the shifted ensemble behavior); and (3) comparison to experimentally determined pK(a) values (where the calculation should reproduce the pK(a) within experimental error). Both the general behavior of cysteines in proteins and the perturbed pK(a) in some proteins can be predicted reasonably well using the newly determined empirical parameters within the MEAD model for protein electrostatics. This study provides the first general analysis of the electrostatics of cysteines in proteins, with specific attention paid to capturing both the behavior of typical cysteines in a protein and the behavior of cysteines whose pK(a) should be shifted, and validation of force field parameters for cysteine residues.

Thioredoxin 1-mediated post-translational modifications: reduction, transnitrosylation, denitrosylation, and related proteomics methodologies

Despite the significance of redox post-translational modifications (PTMs) in regulating diverse signal transduction pathways, the enzymatic systems that catalyze reversible and specific oxidative or reductive modifications have yet to be firmly established. Thioredoxin 1 (Trx1) is a conserved antioxidant protein that is well known for its disulfide reductase activity. Interestingly, Trx1 is also able to transnitrosylate or denitrosylate (defined as processes to transfer or remove a nitric oxide entity to/from substrates) specific proteins. An intricate redox regulatory mechanism has recently been uncovered that accounts for the ability of Trx1 to catalyze these different redox PTMs. In this review, we will summarize the available evidence in support of Trx1 as a specific disulfide reductase, and denitrosylation and transnitrosylation agent, as well as the biological significance of the diverse array of Trx1-regulated pathways and processes under different physiological contexts. The dramatic progress in redox proteomics techniques has enabled the identification of an increasing number of proteins, including peroxiredoxin 1, whose disulfide bond formation and nitrosylation status are regulated by Trx1. This review will also summarize the advancements of redox proteomics techniques for the identification of the protein targets of Trx1-mediated PTMs. Collectively, these studies have shed light on the mechanisms that regulate Trx1-mediated reduction, transnitrosylation, and denitrosylation of specific target proteins, solidifying the role of Trx1 as a master regulator of redox signal transduction.

Multiple pH regime molecular dynamics simulation for pK calculations

Ionisation equilibria in proteins are influenced by conformational flexibility, which can in principle be accounted for by molecular dynamics simulation. One problem in this method is the bias arising from the fixed protonation state during the simulation. Its effect is mostly exhibited when the ionisation behaviour of the titratable groups is extrapolated to pH regions where the predetermined protonation state of the protein may not be statistically relevant, leading to conformational sampling that is not representative of the true state. In this work we consider a simple approach which can essentially reduce this problem. Three molecular dynamics structure sets are generated, each with a different protonation state of the protein molecule expected to be relevant at three pH regions, and pK calculations from the three sets are combined to predict pK over the entire pH range of interest. This multiple pH molecular dynamics approach was tested on the GCN4 leucine zipper, a protein for which a full data set of experimental data is available. The pK values were predicted with a mean deviation from the experimental data of 0.29 pH units, and with a precision of 0.13 pH units, evaluated on the basis of equivalent sites in the dimeric GCN4 leucine zipper.

Study of the thiol/disulfide redox systems of the anaerobe Desulfovibrio vulgaris points out pyruvate:ferredoxin oxidoreductase as a new target for thioredoxin 1

Sulfate reducers have developed a multifaceted adaptative strategy to survive against oxidative stresses. Along with this oxidative stress response, we recently characterized an elegant reversible disulfide bond-dependent protective mechanism in the pyruvate:ferredoxin oxidoreductase (PFOR) of various Desulfovibrio species. Here, we searched for thiol redox systems involved in this mechanism. Using thiol fluorescent labeling, we show that glutathione is not the major thiol/disulfide balance-controlling compound in four different Desulfovibrio species and that no other plentiful low molecular weight thiol can be detected. Enzymatic analyses of two thioredoxins (Trxs) and three thioredoxin reductases allow us to propose the existence of two independent Trx systems in Desulfovibrio vulgaris Hildenborough (DvH). The TR1/Trx1 system corresponds to the typical bacterial Trx system. We measured a TR1 apparent K(m) value for Trx1 of 8.9 μM. Moreover, our results showed that activity of TR1 was NADPH-dependent. The second system named TR3/Trx3 corresponds to an unconventional Trx system as TR3 used preferentially NADH (K(m) for NADPH, 743 μM; K(m) for NADH, 5.6 μM), and Trx3 was unable to reduce insulin. The K(m) value of TR3 for Trx3 was 1.12 μM. In vitro experiments demonstrated that the TR1/Trx1 system was the only one able to reactivate the oxygen-protected form of Desulfovibrio africanus PFOR. Moreover, ex vivo pulldown assays using the mutant Trx1(C33S) as bait allowed us to capture PFOR from the DvH extract. Altogether, these data demonstrate that PFOR is a new target for Trx1, which is probably involved in the protective switch mechanism of the enzyme.

A residue outside the active site CXXC motif regulates the catalytic efficiency of Glutaredoxin 3

The glutaredoxin (Grx) family of oxidoreductases has a conserved residue at position 8 that varies between Arginine in Grx1 and Lysine in Grx3. It has been proposed that this Arg/Lys change is the main cause for the 35 mV difference in redox potential between the two enzymes. To gain insights into the catalytic machinery of Grx3 and directly evaluate the role of residue 8 in the catalysis of thiol-disulfide exchange by this enzyme, we synthesized the "wild type" enzyme (sGrx3), and four analogues substituting the lysine at position 8 with arginine, ornithine (Orn), citrulline (Cit) and norvaline (Nva). The redox potential and equilibration kinetics with thioredoxin (Trx1) were determined for each enzyme by fluorescence intensity. While minor effects on redox potential were observed, we found that residue 8 had a more marked effect on the catalytic efficiency of this enzyme. Surprisingly, truncation of the functional group resulted in a more efficient enzyme, Lys8Nva, exhibiting rate constants that are an order of magnitude higher than sGrx3 for both forward and reverse reactions. These observations pose the question why would a residue that reduces the rate of enzyme turnover be evolutionarily conserved? The significant changes in the kinetic parameters suggest that this position plays an important role in the thiol-disulfide exchange reaction by affecting the nucleophilic thiolate through electrostatic or hydrogen bonding interactions. Since the reduced Grx has an exposed thiol that could easily be alkylated, either Arg or Lys could act as a gatekeeper that deters unwanted electrophiles from attacking the active site thiolate.

How thioredoxin dissociates its mixed disulfide

The dissociation mechanism of the thioredoxin (Trx) mixed disulfide complexes is unknown and has been debated for more than twenty years. Specifically, opposing arguments for the activation of the nucleophilic cysteine as a thiolate during the dissociation of the complex have been put forward. As a key model, the complex between Trx and its endogenous substrate, arsenate reductase (ArsC), was used. In this structure, a Cys29^Trx-Cys89^ArsC intermediate disulfide is formed by the nucleophilic attack of Cys29^Trx on the exposed Cys82^ArsC-Cys89^ArsC in oxidized ArsC. With theoretical reactivity analysis, molecular dynamics simulations, and biochemical complex formation experiments with Cys-mutants, Trx mixed disulfide dissociation was studied. We observed that the conformational changes around the intermediate disulfide bring Cys32^Trx in contact with Cys29^Trx. Cys32^Trx is activated for its nucleophilic attack by hydrogen bonds, and Cys32^Trx is found to be more reactive than Cys82^ArsC. Additionally, Cys32^Trx directs its nucleophilic attack on the more susceptible Cys29^Trx and not on Cys89^ArsC. This multidisciplinary approach provides fresh insights into a universal thiol/disulfide exchange reaction mechanism that results in reduced substrate and oxidized Trx.

Thioredoxins are found in all types of cells and control several essential functions of life, including promotion of cell growth, inhibition of apoptosis, and modulation of inflammation. Thioredoxin has two ‘free’ cysteines in its active site, which are used to break disulfide bonds in oxidized substrate proteins. In the first step, an intermediate thioredoxin-protein complex is formed, which is broken in the second step, releasing the substrate protein in its reduced state. In other words, the disulfide is being transferred from the substrate protein to thioredoxin, or the electrons coming from thioredoxin are shuttled to the protein substrate. The exact reaction mechanism, i.e., the detailed succession of steps in which the reaction takes place, of how this mixed disulfide is broken is not known and has been debated over the last twenty years. With a multidisciplinary approach, combining computational and experimental work, we provide fresh insights into how conformational changes activate the catalytic cysteines with which this universal reduction mechanism is completed with the correct regioselectivity. This work illustrates the strengths of computational approaches in probing phenomena which are otherwise very difficult to investigate experimentally.

CHARMM: the biomolecular simulation program

CHARMM (Chemistry at HARvard Molecular Mechanics) is a highly versatile and widely used molecular simulation program. It has been developed over the last three decades with a primary focus on molecules of biological interest, including proteins, peptides, lipids, nucleic acids, carbohydrates, and small molecule ligands, as they occur in solution, crystals, and membrane environments. For the study of such systems, the program provides a large suite of computational tools that include numerous conformational and path sampling methods, free energy estimators, molecular minimization, dynamics, and analysis techniques, and model-building capabilities. The CHARMM program is applicable to problems involving a much broader class of many-particle systems. Calculations with CHARMM can be performed using a number of different energy functions and models, from mixed quantum mechanical-molecular mechanical force fields, to all-atom classical potential energy functions with explicit solvent and various boundary conditions, to implicit solvent and membrane models. The program has been ported to numerous platforms in both serial and parallel architectures. This article provides an overview of the program as it exists today with an emphasis on developments since the publication of the original CHARMM article in 1983.

De novo design and evolution of artificial disulfide isomerase enzymes analogous to the bacterial DsbC

The Escherichia coli disulfide isomerase, DsbC is a V-shaped homodimer with each monomer comprising a dimerization region that forms part of a putative peptide-binding pocket and a thioredoxin catalytic domain. Disulfide isomerases from prokaryotes and eukaryotes exhibit little sequence homology but display very similar structural organization with two thioredoxin domains facing each other on top of the dimerization/peptide-binding region. To aid the understanding of the mechanistic significance of thioredoxin domain dimerization and of the peptide-binding cleft of DsbC, we constructed a series of protein chimeras comprising unrelated protein dimerization domains fused to thioredoxin superfamily enzymes. Chimeras consisting of the dimerization domain and the alpha-helical linker of the bacterial proline cis/trans isomerase FkpA and the periplasmic oxidase DsbA gave rise to enzymes that catalyzed the folding of multidisulfide substrate proteins in vivo with comparable efficiency to E. coli DsbC. In addition, expression of FkpA-DsbAs conferred modest resistance to CuCl2, a phenotype that depends on disulfide bond isomerization. Selection for resistance to elevated CuCl2 concentrations led to the isolation of FkpA-DsbA mutants containing a single amino acid substitution that changed the active site of the DsbA domain from CPHC into CPYC, increasing the similarity to the DsbC active site (CGYC). Unlike DsbC, which is resistant to oxidation by DsbB-DsbA and does not normally catalyze disulfide bond formation under physiological conditions, the FkpA-DsbA chimeras functioned both as oxidases and isomerases. The engineering of these efficient artificial isomerases delineates the key features of catalysis of disulfide bond isomerization and enhances our understanding of its evolution.

Probing the chemistry of thioredoxin catalysis with force

Thioredoxins are enzymes that catalyse disulphide bond reduction in all living organisms. Although catalysis is thought to proceed through a substitution nucleophilic bimolecular (S(N)2) reaction, the role of the enzyme in modulating this chemical reaction is unknown. Here, using single-molecule force-clamp spectroscopy, we investigate the catalytic mechanism of Escherichia coli thioredoxin (Trx). We applied mechanical force in the range of 25-600 pN to a disulphide bond substrate and monitored the reduction of these bonds by individual enzymes. We detected two alternative forms of the catalytic reaction, the first requiring a reorientation of the substrate disulphide bond, causing a shortening of the substrate polypeptide by 0.79 +/- 0.09 A (+/- s.e.m.), and the second elongating the substrate disulphide bond by 0.17 +/- 0.02 A (+/- s.e.m.). These results support the view that the Trx active site regulates the geometry of the participating sulphur atoms with sub-ångström precision to achieve efficient catalysis. Our results indicate that substrate conformational changes may be important in the regulation of Trx activity under conditions of oxidative stress and mechanical injury, such as those experienced in cardiovascular disease. Furthermore, single-molecule atomic force microscopy techniques, as shown here, can probe dynamic rearrangements within an enzyme's active site during catalysis that cannot be resolved with any other current structural biological technique.

The reducing activity of glutaredoxin 3 toward cytoplasmic substrate proteins is restricted by methionine 43

The reducing proteins glutaredoxin 3 (Grx3) and glutaredoxin 1 (Grx1) are structurally similar but exhibit different specificities toward substrates. Grx1 efficiently reduces ribonucleotide reductase and PAPS reductase, while Grx3 reduces these enzymes inefficiently or not at all. We previously described a selection for Grx3 mutants with increased activity toward substrates of Grx1 in vivo. Remarkably, we repeatedly isolated mutants with changes in only one of the amino acids of Grx3, methionine 43, converting it to either valine, leucine, or isoleucine. In this paper we present additional genetic studies and a biochemical characterization of Grx3-Met43Val, the most efficient mutant. We show that Grx3-Met43Val is able to reduce ribonucleotide reductae and PAPS reductase much more efficiently than the wild-type protein in vitro. The altered protein has an increased Vmax over that of Grx3, nearly the same Vmax as Grx1 while the Km remains high. Molecular dynamics simulations suggest that the Met43Val substitution results in changes in properties of the N-terminal cysteine of the active site leading to a considerably lower pKa. Furthermore, Grx3-Met43Val shows an 11 mV lower redox potential than the wild-type Grx3. These findings provide biochemical and structural explanations for the increased reductive efficiency of the mutant Grx3.

The origami of thioredoxin-like folds

Origami is the Japanese art of folding a piece of paper into complex shapes and forms. Much like origami of paper, Nature has used conserved protein folds to engineer proteins for a particular task. An example of a protein family, which has been used by Nature numerous times, is the thioredoxin superfamily. Proteins in the thioredoxin superfamily are all structured with a beta-sheet core surrounded with alpha-helices, and most contain a canonical CXXC motif. The remarkable feature of these proteins is that the link between them is the fold; however, their reactivity is different for each member due to small variations in this general fold as well as their active site. This review attempts to unravel the minute differences within this protein family, and it also demonstrates the ingenuity of Nature to use a conserved fold to generate a diverse collection of proteins to perform a number of different biochemical tasks.

Oxidized and synchrotron cleaved structures of the disulfide redox center in the N-terminal domain of Salmonella typhimurium AhpF

The flavoprotein component (AhpF) of Salmonella typhimurium alkyl hydroperoxide reductase contains an N-terminal domain (NTD) with two contiguous thioredoxin folds but only one redox-active disulfide (within the sequence -Cys129-His-Asn-Cys132-). This active site is responsible for mediating the transfer of electrons from the thioredoxin reductase-like segment of AhpF to AhpC, the peroxiredoxin component of the two-protein peroxidase system. The previously reported crystal structure of AhpF possessed a reduced NTD active site, although fully oxidized protein was used for crystallization. To further investigate this active site, we crystallized an isolated recombinant NTD (rNTD); using diffraction data sets collected first at our in-house X-ray source and subsequently at a synchrotron, we showed that the active site disulfide bond (Cys129-Cys132) is oxidized in the native crystals but becomes reduced during synchrotron data collection. The NTD disulfide bond is apparently particularly sensitive to radiation cleavage compared with other protein disulfides. The two data sets provide the first view of an oxidized (disulfide) form of NTD and show that the changes in conformation upon reduction of the disulfide are localized and small. Furthermore, we report the apparent pKa of the active site thiol to be approximately 5.1, a relatively low pKa given its redox potential (approximately 265 mV) compared with most members of the thioredoxin family.

Prediction of pKa and redox properties in the thioredoxin superfamily

Electrostatic interactions play important roles in diverse biological phenomena controlling the function of many proteins. Polar molecules can be studied with the FDPB method solving the Poisson-Boltzmann equation on a finite difference grid. A method for the prediction of pK(a)s and redox potentials in the thioredoxin superfamily is introduced. The results are compared with experimental pK(a) data where available, and predictions are made for members lacking such data. Studying CxxC motif variation in the context of different background structures permits analysis of contributions to cysteine DeltapK(a)s. The motif itself and the overall framework regulate pK(a) variation. The reported method includes generation of multiple side-chain rotamers for the CxxC motif and is an effective predictive tool for functional pK(a) variation across the superfamily. Redox potential follows the trend in cysteine pK(a) variation, but some residual discrepancy indicates that a pH-independent factor plays a role in determining redox potentials for at least some members of the superfamily. A possible molecular basis for this feature is discussed.