Quantifying redox transcription factor dynamics as a tool to investigate redox signalling.

A critical feature of the cellular antioxidant response is the induction of gene expression by redox-sensitive transcription factors. In many cells, activating these transcription factors is a dynamic process involving multiple redox steps, but it is unclear how these dynamics should be measured. Here, we show how the dynamic profile of the Schizosaccharomyces pombe Pap1 transcription factor is quantifiable by three parameters: signal amplitude, signal time and signal duration. In response to increasing hydrogen peroxide concentrations, the Pap1 amplitude decreased while the signal time and duration showed saturable increases. In co-response plots, these parameters showed a complex, non-linear relationship to the mRNA levels of four Pap1-regulated genes. We also demonstrate that hydrogen peroxide and tert-butyl hydroperoxide trigger quantifiably distinct Pap1 activation profiles and transcriptional responses. Based on these findings, we propose that different oxidants and oxidant concentrations modulate the Pap1 dynamic profile, leading to specific transcriptional responses. We further show how the effect of combination and pre-exposure stresses on Pap1 activation dynamics can be quantified using this approach. This method is therefore a valuable addition to the redox signalling toolbox that may illuminate the role of dynamics in determining appropriate responses to oxidative stress.

Computational models as catalysts for investigating redoxin systems.

Thioredoxin, glutaredoxin and peroxiredoxin systems play central roles in redox regulation, signaling and metabolism in cells. In these systems, reducing equivalents from NAD(P)H are transferred by coupled thiol-disulfide exchange reactions to redoxins which then reduce a wide array of targets. However, the characterization of redoxin activity has been unclear, with redoxins regarded as enzymes in some studies and redox metabolites in others. Consequently, redoxin activities have been quantified by enzyme kinetic parameters in vitro, and redox potentials or redox ratios within cells. By analyzing all the reactions within these systems, computational models showed that many kinetic properties attributed to redoxins were due to system-level effects. Models of cellular redoxin networks have also been used to estimate intracellular hydrogen peroxide levels, analyze redox signaling and couple omic and kinetic data to understand the regulation of these networks in disease. Computational modeling has emerged as a powerful complementary tool to traditional redoxin enzyme kinetic and cellular assays that integrates data from a number of sources into a single quantitative framework to accelerate the analysis of redoxin systems.

Modelling the Decamerisation Cycle of PRDX1 and the Inhibition-like Effect on Its Peroxidase Activity.

Peroxiredoxins play central roles in the detoxification of reactive oxygen species and have been modelled across multiple organisms using a variety of kinetic methods. However, the peroxiredoxin dimer-to-decamer transition has been underappreciated in these studies despite the 100-fold difference in activity between these forms. This is due to the lack of available kinetics and a theoretical framework for modelling this process. Using published isothermal titration calorimetry data, we obtained association and dissociation rate constants of 0.050 µM-4·s-1 and 0.055 s-1, respectively, for the dimer-decamer transition of human PRDX1. We developed an approach that greatly reduces the number of reactions and species needed to model the peroxiredoxin decamer oxidation cycle. Using these data, we simulated horse radish peroxidase competition and NADPH-oxidation linked assays and found that the dimer-decamer transition had an inhibition-like effect on peroxidase activity. Further, we incorporated this dimer-decamer topology and kinetics into a published and validated in vivo model of PRDX2 in the erythrocyte and found that it almost perfectly reconciled experimental and simulated responses of PRDX2 oxidation state to hydrogen peroxide insult. By accounting for the dimer-decamer transition of peroxiredoxins, we were able to resolve several discrepancies between experimental data and available kinetic models.

Atypical network topologies enhance the reductive capacity of pathogen thiol antioxidant defense networks.

Infectious diseases are a significant health burden for developing countries, particularly with the rise of multidrug resistance. There is an urgent need to elucidate the factors underlying the persistence of pathogens such as Mycobacterium tuberculosis, Plasmodium falciparum and Trypanosoma brucei. In contrast to host cells, these pathogens traverse multiple and varied redox environments during their infectious cycles, including exposure to high levels of host-derived reactive oxygen species. Pathogen antioxidant defenses such as the peroxiredoxin and thioredoxin systems play critical roles in the redox stress tolerance of these cells. However, many of the kinetic rate constants obtained for the pathogen peroxiredoxins are broadly similar to their mammalian homologs and therefore, their contributions to the redox tolerances within these cells are enigmatic. Using graph theoretical analysis, we show that compared to a canonical Escherichia coli redoxin network, pathogen redoxin networks contain unique network connections (motifs) between their thioredoxins and peroxiredoxins. Analysis of these motifs reveals that they increase the hydroperoxide reduction capacity of these networks and, in response to an oxidative insult, can distribute fluxes into specific thioredoxin-dependent pathways. Our results emphasize that the high oxidative stress tolerance of these pathogens depends on both the kinetic parameters for hydroperoxide reduction and the connectivity within their thioredoxin/peroxiredoxin systems.

Can thiol-based redox systems be utilized as parts for synthetic biology applications?

OBJECTIVES: Synthetic biology has emerged from molecular biology and engineering approaches and aims to develop novel, biologically-inspired systems for industrial and basic research applications ranging from biocomputing to drug production. Surprisingly, redoxin (thioredoxin, glutaredoxin, peroxiredoxin) and other thiol-based redox systems have not been widely utilized in many of these synthetic biology applications. METHODS: We reviewed thiol-based redox systems and the development of synthetic biology applications that have used thiol-dependent parts. RESULTS: The development of circuits to facilitate cytoplasmic disulfide bonding, biocomputing and the treatment of intestinal bowel disease are amongst the applications that have used thiol-based parts. We propose that genetically encoded redox sensors, thiol-based biomaterials and intracellular hydrogen peroxide generators may also be valuable components for synthetic biology applications. DISCUSSION: Thiol-based systems play multiple roles in cellular redox metabolism, antioxidant defense and signaling and could therefore offer a vast and diverse portfolio of components, parts and devices for synthetic biology applications. However, factors limiting the adoption of redoxin systems for synthetic biology applications include the orthogonality of thiol-based components, limitations in the methods to characterize thiol-based systems and an incomplete understanding of the design principles of these systems.

Analyzing biological models and data sets using Jupyter notebooks as an alternate to laboratory-based exercises during COVID-19.

Jupyter notebooks are widely used for data analysis across a large number of scientific disciplines. As a result of the COVID-19 pandemic, I developed a series of computational exercises using the Jupyter notebook to replace the laboratory exercises usually undertaken in my course. My students had no prior coding knowledge and therefore these exercises were structured in a "cookbook" format using the susceptible-infected-resistant model for disease, data from the Lenski long-term evolutionary experiment, and a fission yeast transcriptomic data set. Despite limited internet connectivity and on-line instruction, my students completed these computational exercises and then tested their own hypotheses. Because Jupyter notebooks can be annotated with text and images, student notebooks were submitted for assessment in the form of a structured scientific report. An advantage of this approach was that all the computational analyses presented in these reports could be easily replicated. The notebook and complete instructions used in my course are provided for others who want to adopt this approach.

The thioredoxin redox potential and redox charge are surrogate measures for flux in the thioredoxin system.

The thioredoxin system plays a central role in intracellular redox regulation and its dysregulation is associated with a number of pathologies. However, the connectivity within this system poses a significant challenge for quantification and consequently several disparate measures have been used to characterize the system. For in vitro studies, the thioredoxin system flux has been measured by NADPH oxidation while the thioredoxin redox state has been used to estimate the activity of the system in vivo. The connection between these measures has been obscure although substrate saturation in the thioredoxin system results from the saturation of the thioredoxin redox cycle. We used computational modeling and in vitro kinetic assays to clarify the relationship between flux and the current in vivo measures of the thioredoxin system together with a novel measure, the thioredoxin redox charge (reduced thioredoxin/total thioredoxin). Our results revealed that the thioredoxin redox potential and redox charge closely tracked flux perturbations showing that these indices could be used as surrogate measures of the flux in vivo and, provide a mechanistic explanation for the previously observed correlations between thioredoxin oxidation and certain pathologies. While we found no significant difference in the linear correlations obtained for the thioredoxin redox potential and redox charge with the flux, the redox charge may be preferred because it is bounded between zero and one and can be determined over a wider range of conditions allowing for quantitative flux comparisons between cell types and conditions.

Hyperoxidation of Peroxiredoxins: Gain or Loss of Function?

SIGNIFICANCE: In 2003, structural studies revealed that eukaryotic 2-Cys peroxiredoxins (Prx) have evolved to be sensitive to inactivation of their thioredoxin peroxidase activity by hyperoxidation (sulfinylation) of their peroxide-reacting catalytic cysteine. This was accompanied by the unexpected discovery, that the sulfinylation of this cysteine was reversible in vivo and the identification of a new enzyme, sulfiredoxin, that had apparently co-evolved specifically to reduce hyperoxidized 2-Cys Prx, restoring their peroxidase activity. Together, these findings have provided the impetus for multiple studies investigating the purpose of this reversible, Prx hyperoxidation. Recent Advances: It has been suggested that inhibition of the thioredoxin peroxidase activity by hyperoxidation can both promote and inhibit peroxide signal transduction, depending on the context. Prx hyperoxidation has also been proposed to protect cells against reactive oxygen species (ROS)-induced damage, by preserving reduced thioredoxin and/or by increasing non-peroxidase chaperone or signaling activities of Prx. CRITICAL ISSUES: Here, we will review the evidence in support of each of these proposed functions, in view of the in vivo contexts in which Prx hyperoxidation occurs, and the role of sulfiredoxin. Thus, we will attempt to explain the basis for seemingly contradictory roles for Prx hyperoxidation in redox signaling. FUTURE DIRECTIONS: We provide a rationale, based on modeling and experimental studies, for why Prx hyperoxidation should be considered a suitable, early biomarker for damaging levels of ROS. We discuss the implications that this has for the role of Prx in aging and the detection of hyperoxidized Prx as a conserved feature of circadian rhythms. Antioxid. Redox Signal. 28, 574-590.

Quantitative measures for redox signaling.

Redox signaling is now recognized as an important regulatory mechanism for a number of cellular processes including the antioxidant response, phosphokinase signal transduction and redox metabolism. While there has been considerable progress in identifying the cellular machinery involved in redox signaling, quantitative measures of redox signals have been lacking, limiting efforts aimed at understanding and comparing redox signaling under normoxic and pathogenic conditions. Here we have outlined some of the accepted principles for redox signaling, including the description of hydrogen peroxide as a signaling molecule and the role of kinetics in conferring specificity to these signaling events. Based on these principles, we then develop a working definition for redox signaling and review a number of quantitative methods that have been employed to describe signaling in other systems. Using computational modeling and published data, we show how time- and concentration- dependent analyses, in particular, could be used to quantitatively describe redox signaling and therefore provide important insights into the functional organization of redox networks. Finally, we consider some of the key challenges with implementing these methods.

The Development of Computational Biology in South Africa: Successes Achieved and Lessons Learnt.

Bioinformatics is now a critical skill in many research and commercial environments as biological data are increasing in both size and complexity. South African researchers recognized this need in the mid-1990s and responded by working with the government as well as international bodies to develop initiatives to build bioinformatics capacity in the country. Significant injections of support from these bodies provided a springboard for the establishment of computational biology units at multiple universities throughout the country, which took on teaching, basic research and support roles. Several challenges were encountered, for example with unreliability of funding, lack of skills, and lack of infrastructure. However, the bioinformatics community worked together to overcome these, and South Africa is now arguably the leading country in bioinformatics on the African continent. Here we discuss how the discipline developed in the country, highlighting the challenges, successes, and lessons learnt.

The thioredoxin system and not the Michaelis-Menten equation should be fitted to substrate saturation datasets from the thioredoxin insulin assay.

INTRODUCTION: The thioredoxin system, consisting of thioredoxin reductase, thioredoxin and NADPH, is present in most living organisms and reduces a large array of target protein disulfides. OBJECTIVE: The insulin reduction assay is commonly used to characterise thioredoxin activity in vitro, but it is not clear whether substrate saturation datasets from this assay should be fitted and modeled with the Michaelis-Menten equation (thioredoxin enzyme model), or fitted to the thioredoxin system with insulin reduction described by mass-action kinetics (redox couple model). METHODS: We utilized computational modeling and in vitro assays to determine which of these approaches yield consistent and accurate kinetic parameter sets for insulin reduction. RESULTS: Using computational modeling, we found that fitting to the redox couple model, rather than to the thioredoxin enzyme model, resulted in consistent parameter sets over a range of thioredoxin reductase concentrations. Furthermore, we established that substrate saturation in this assay was due to the progressive redistribution of the thioredoxin moiety into its oxidised form. We then confirmed these results in vitro using the yeast thioredoxin system. DISCUSSION: This study shows how consistent parameter sets for thioredoxin activity can be obtained regardless of the thioredoxin reductase concentration used in the insulin reduction assay, and validates computational systems biology modeling studies that have described the thioredoxin system with the redox couple modeling approach.

The glutaredoxin mono- and di-thiol mechanisms for deglutathionylation are functionally equivalent: implications for redox systems biology.

Glutathionylation plays a central role in cellular redox regulation and anti-oxidative defence. Grx (Glutaredoxins) are primarily responsible for reversing glutathionylation and their activity therefore affects a range of cellular processes, making them prime candidates for computational systems biology studies. However, two distinct kinetic mechanisms involving either one (monothiol) or both (dithiol) active-site cysteines have been proposed for their deglutathionylation activity and initial studies predicted that computational models based on either of these mechanisms will have different structural and kinetic properties. Further, a number of other discrepancies including the relative activity of active-site mutants and contrasting reciprocal plot kinetics have also been reported for these redoxins. Using kinetic modelling, we show that the dithiol and monothiol mechanisms are identical and, we were also able to explain much of the discrepant data found within the literature on Grx activity and kinetics. Moreover, our results have revealed how an apparently futile side-reaction in the monothiol mechanism may play a significant role in regulating Grx activity in vivo.

From top-down to bottom-up: computational modeling approaches for cellular redoxin networks.

SIGNIFICANCE: Thioredoxin, glutaredoxin, and peroxiredoxin systems play critical roles in a large number of redox-sensitive cellular processes. These systems are linked to each other by coupled redox cycles and common reaction intermediates into a larger network. Given the scale and connectivity of this network, computational approaches are required to analyze its dynamics and organization. RECENT ADVANCES: Theoretical advances, as well as new redox proteomic methods, have led to the development of both top-down and bottom-up systems biology approaches to analyze the these systems and the network as a whole. Top-down approaches have been based on modifications to the Nernst equation or on graph theoretical approaches, while bottom-up approaches have been based on kinetic or stoichiometric modeling techniques. CRITICAL ISSUES: This review will consider the rationale behind these approaches and focus on their advantages and limitations. Further, the review will discuss modeling standards to ensure model accuracy and availability. FUTURE DIRECTIONS: Top-down and bottom-up approaches have distinct strengths and limitations in describing cellular redoxin networks. The availability of methods to overcome these limitations, together with the adoption of common modeling standards, is expected to increase the pace of model-led discovery within the redox biology field.

The logic of kinetic regulation in the thioredoxin system.

BACKGROUND: The thioredoxin system consisting of NADP(H), thioredoxin reductase and thioredoxin provides reducing equivalents to a large and diverse array of cellular processes. Despite a great deal of information on the kinetics of individual thioredoxin-dependent reactions, the kinetic regulation of this system as an integrated whole is not known. We address this by using kinetic modeling to identify and describe kinetic behavioral motifs found within the system. RESULTS: Analysis of a realistic computational model of the Escherichia coli thioredoxin system revealed several modes of kinetic regulation in the system. In keeping with published findings, the model showed that thioredoxin-dependent reactions were adaptable (i.e. changes to the thioredoxin system affected the kinetic profiles of these reactions). Further and in contrast to other systems-level descriptions, analysis of the model showed that apparently unrelated thioredoxin oxidation reactions can affect each other via their combined effects on the thioredoxin redox cycle. However, the scale of these effects depended on the kinetics of the individual thioredoxin oxidation reactions with some reactions more sensitive to changes in the thioredoxin cycle and others, such as the Tpx-dependent reduction of hydrogen peroxide, less sensitive to these changes. The coupling of the thioredoxin and Tpx redox cycles also allowed for ultrasensitive changes in the thioredoxin concentration in response to changes in the thioredoxin reductase concentration. We were able to describe the kinetic mechanisms underlying these behaviors precisely with analytical solutions and core models. CONCLUSIONS: Using kinetic modeling we have revealed the logic that underlies the functional organization and kinetic behavior of the thioredoxin system. The thioredoxin redox cycle and associated reactions allows for a system that is adaptable, interconnected and able to display differential sensitivities to changes in this redox cycle. This work provides a theoretical, systems-biological basis for an experimental analysis of the thioredoxin system and its associated reactions.

Enzymes or redox couples? The kinetics of thioredoxin and glutaredoxin reactions in a systems biology context.

Systems biology approaches, such as kinetic modelling, could provide valuable insights into how thioredoxins, glutaredoxins and peroxiredoxins (here collectively called redoxins), and the systems that reduce these molecules are regulated. However, it is not clear whether redoxins should be described as redox couples (with redox potentials) or as enzymes (with Michaelis-Menten parameters) in such approaches. We show that in complete redoxin systems, redoxin substrate saturation and other purported enzymatic behaviours result from limitations in the redoxin redox cycles in these systems. Michaelis-Menten parameters are therefore inappropriate descriptors of redoxin activity; data from redoxin kinetic experiments should rather be interpreted in terms of the complete system of reactions under study. These findings were confirmed by fitting kinetic models of the thioredoxin and glutaredoxin systems to in vitro datasets. This systems approach clarifies the inconsistencies with the descriptions of redoxins and emphasizes the roles of redoxin systems in redox regulation.

Cathepsin B stability, but not activity, is affected in cysteine:cystine redox buffers.

In order to test the hypothesis that the lysosomal cysteine protease cathepsin B may be redox regulated in vivo, cathepsin B activity and stability were measured in cysteine- and/or cystine-containing buffers. Cathepsin B activity in cysteine-containing buffers was similar at pH 6.0 and pH 7.0, over all thiol concentrations tested. In contrast, the stability of the enzyme was greater at pH 6.0 than at pH 7.0. This suggests that the enzyme's operational pH in vivo may be < pH 7.0. The activity of the enzyme was depressed in glutathione-containing buffers. When assessed in cysteine:cystine redox buffers (pH 6.0-7.0) cathepsin B was active over a broad redox potential range, suggesting that cathepsin B activity may not be redox regulated. However, at pH 7.0, the stability of cathepsin B decreased with increasing reduction potential and ambient cystine concentration. This suggests that the stability of the enzyme at neutral pH is dependent on redox potential, and on the presence of oxidising agents.

Endolysosomal proteolysis and its regulation.

The endolysosomal system comprises a unique environment for proteolysis, which is regulated in a manner that apparently does not involve protease inhibitors. The system comprises a series of membrane-bound intracellular compartments, within which endocytosed material and redundant cellular components are hydrolysed. Endocytosed material tends to flow vectorially through the system, proceeding through the early endosome, the endosome carrier vesicle, the late endosome and the lysosome. Phagocytosis and autophagy provide alternative entry points into the system. Late endosomes, lysosome/late endosome hybrid organelles, phagosomes and autophagosomes are the principal sites for proteolysis. In each case, hydrolytic competence is due to components of the endolysosomal system, i.e. proteases, lysosome-associated membrane proteins, H(+)-ATPases and possibly cysteine transporters. The view is emerging that lysosomes are organelles for the storage of hydrolases, perhaps in an inactivated form. Once a substrate has entered a proteolytically competent environment, the rate-limiting proteolytic steps are probably effected by cysteine endoproteinases. As these are affected by pH and possibly redox potential, they may be regulated by the organelle luminal environment. Regulation is probably also affected, among other factors, by organelle fusion reactions, whereby the meeting of enzyme and substrate may be controlled. Such systems would permit simultaneous regulation of a number of unrelated hydrolases.