Integrating mitochondrial energetics, redox and ROS metabolic networks: a two-compartment model

A computational model of the crosstalk between hepatocyte fatty acid metabolism and oxidative stress highlights the key enzymes, metabolites, and detoxification pathways in the context of MASLD

Metabolic dysfunction-associated steatotic liver disease (MASLD; formerly known as NAFLD) is a common liver disease worldwide and carries the risk of progressing to severe liver conditions, such as fibrosis and liver cancer. In the context of MASLD, evaluating fat accumulation in the liver and the subsequent production of oxidative stress is essential to understand the disease propagation. However, clinical studies using human patients to investigate the fat accumulation and the onset of oxidative stress in MASLD face ethical and technical challenges, highlighting the importance of alternative methods. To understand the relationship between fatty acid metabolism, lipid accumulation, oxidative stress generation, and antioxidant mechanisms in hepatocytes, we proposed a new mathematical model. The importance of this model lies in its ability to track the time-dependent changes in oxidative stress and glutathione concentration in response to the input of fatty acids. Furthermore, the model allows for the evaluation of the effects of altering the activity of the key enzymes involved in those mechanisms. Our model is anticipated to provide new insights into MASLD therapy strategies by identifying key pathways and predicting the effects of drug-induced changes in enzyme activity.

Tafazzin deficiency causes substantial remodeling in the lipidome of a mouse model of Barth Syndrome cardiomyopathy

Barth Syndrome (BTHS) is a rare X-linked disease, characterized clinically by cardiomyopathy, skeletal myopathy, neutropenia, and growth retardation. BTHS is caused by mutations in the phospholipid acyltransferase tafazzin (Gene: TAFAZZIN, TAZ). Tafazzin catalyzes the final step in the remodeling of cardiolipin (CL), a glycerophospholipid located in the inner mitochondrial membrane. As the phospholipid composition strongly determines membrane properties, correct biosynthesis of CL and other membrane lipids is essential for mitochondrial function. Mitochondria provide 95% of the energy demand in the heart, particularly due to their role in fatty acid oxidation. Alterations in lipid homeostasis in BTHS have an impact on mitochondrial membrane proteins and thereby contribute to cardiomyopathy. We analyzed a transgenic TAFAZZIN-knockdown (TAZ-KD) BTHS mouse model and determined the distribution of 193 individual lipid species in TAZ-KD and WT hearts at 10 and 50 weeks of age, using electrospray ionization tandem mass spectrometry (ESI-MS/MS). Our results revealed significant lipid composition differences between the TAZ-KD and WT groups, indicating genotype-dependent alterations in most analyzed lipid species. Significant changes in the myocardial lipidome were identified in both young animals without cardiomyopathy and older animals with heart failure. Notable alterations were found in phosphatidylcholine (PC), phosphatidylethanolamine (PE), lysophosphatidylethanolamine (LPE), lysophosphatidylcholine (LPC) and plasmalogen species. PC species with 2-4 double bonds were significantly increased, while polyunsaturated PC species showed a significant decrease in TAZ-KD mice. Furthermore, Linoleic acid (LA, 18:2) containing PC and PE species, as well as arachidonic acid (AA, 20:4) containing PE 38:4 species are increased in TAZ-KD. We found higher levels of AA containing LPE and PE-based plasmalogens (PE P-). Furthermore, we are the first to show significant changes in sphingomyelin (SM) and ceramide (Cer) lipid species Very long-chained SM species are accumulating in TAZ-KD hearts, whereas long-chained Cer and several hexosyl ceramides (HexCer) species accumulate only in 50-week-old TAZ-KD hearts These findings offer potential avenues for the diagnosis and treatment of BTHS, presenting new possibilities for therapeutic approaches.

A Computational Model of Endogenous Hydrogen Peroxide Metabolism in Hepatocytes, Featuring a Critical Role for GSH

This paper presents an ordinary differential equation (ODE) model of endogenous H2O2 production and elimination in hepatocytes that is unique, at the time of writing, in its ability to accurately compute intracellular H2O2 concentration during incidents of oxidative stress and in its usefulness for constructing PBPK/PD models for ROS-generating xenobiotics. Versions of the model are presented for rat hepatocytes in vitro and mouse liver in vivo. A generic method is given for using the model to create PBPK/PD models which predict intracellular H2O2 concentration and oxidative-stress-induced hepatocyte death; these are identifiable from in vitro data sets reporting cell mortality following xenobiotic exposure at various levels. The procedure is demonstrated for the trivalent arsenical dimethylarsinous acid (DMAIII), which is produced in liver as part of the arsenic elimination pathway. This is the first model of H2O2 metabolism in hepatocytes to feature values for the endogenous rates of H2O2 production by mitochondria and other organelles which are inferred from the physiology literature, and to feature a detailed, realistic treatment of GSH metabolism; the latter is achieved by incorporating a minimal version of Reed and coworkers' pioneering model of GSH metabolism in liver. Model simulations indicate that critical GSH depletion is the immediate trigger for intracellular H2O2 rising to concentrations associated with apoptosis (> 1 μM), that this may only occur hours after the xenobiotic concentration peaks ("delay effect"), that when critical GSH depletion does occur, H2O2 concentration rises rapidly in a sequence of two boundary layers, characterized by the kinetics of glutathione peroxidase (first boundary layer) and catalase (second boundary layer), and that intracellular H2O2 concentration > 1 μM implies critical GSH depletion. There has been speculation that ROS levels in the range associated with apoptosis simply indicate, rather than cause, an apoptotic milieu. Model simulations are consistent with this view. In a result of interest to the wider physiology community, the delay effect is shown to provide a GSH-based mechanism by which cells can distinguish transient elevations in H2O2 concentration, of use in intracellular signaling, from persistent ones indicative of either pathology or the presence of toxins, the second state of affairs eventually triggering apoptosis.

AW, Dash RK. Computational Modeling of Substrate-Dependent Mitochondrial Respiration and Bioenergetics in the Heart and Kidney Cortex and Outer Medulla

Integrated computational modeling provides a mechanistic and quantitative framework to characterize alterations in mitochondrial respiration and bioenergetics in response to different metabolic substrates in-silico. These alterations play critical roles in the pathogenesis of diseases affecting metabolically active organs such as heart and kidney. Therefore, the present study aimed to develop and validate thermodynamically constrained integrated computational models of mitochondrial respiration and bioenergetics in the heart and kidney cortex and outer medulla (OM). The models incorporated the kinetics of major biochemical reactions and transport processes as well as regulatory mechanisms in the mitochondria of these tissues. Intrinsic model parameters such as Michaelis-Menten constants were fixed at previously estimated values, while extrinsic model parameters such as maximal reaction and transport velocities were estimated separately for each tissue. This was achieved by fitting the model solutions to our recently published respirometry data measured in isolated rat heart and kidney cortex and OM mitochondria utilizing various NADH- and FADH2-linked metabolic substrates. The models were validated by predicting additional respirometry and bioenergetics data, which were not used for estimating the extrinsic model parameters. The models were able to predict tissue-specific and substrate-dependent mitochondrial emergent metabolic system properties such as redox states, enzyme and transporter fluxes, metabolite concentrations, membrane potential, and respiratory control index under diverse physiological and pathological conditions. The models were also able to quantitatively characterize differential regulations of NADH- and FADH2-linked metabolic pathways, which contribute differently toward regulations of oxidative phosphorylation and ATP synthesis in the heart and kidney cortex and OM mitochondria.

Prospects of potential adipokines as therapeutic agents in obesity-linked atherogenic dyslipidemia and insulin resistance

BACKGROUND

In normal circumstances, AT secretes anti-inflammatory adipokines (AAKs) which regulates lipid metabolism, insulin sensitivity, vascular hemostasis, and angiogenesis. However, during obesity AT dysfunction occurs and leads to microvascular imbalance and secretes several pro-inflammatory adipokines (PAKs), thereby favoring atherogenic dyslipidemia and insulin resistance. Literature suggests decreased levels of circulating AAKs and increased levels of PAKs in obesity-linked disorders. Importantly, AAKs have been reported to play a vital role in obesity-linked metabolic disorders mainly insulin resistance, type-2 diabetes mellitus and coronary heart diseases. Interestingly, AAKs counteract the microvascular imbalance in AT and exert cardioprotection via several signaling pathways such as PI3-AKT/PKB pathway. Although literature reviews have presented a number of investigations detailing specific pathways involved in obesity-linked disorders, literature concerning AT dysfunction and AAKs remains sketchy. In view of the above, in the present contribution an effort has been made to provide an insight on the AT dysfunction and role of AAKs in modulating the obesity and obesity-linked atherogenesis and insulin resistance.

MAIN BODY

"Obesity-linked insulin resistance", "obesity-linked cardiometabolic disease", "anti-inflammatory adipokines", "pro-inflammatory adipokines", "adipose tissue dysfunction" and "obesity-linked microvascular dysfunction" are the keywords used for searching article. Google scholar, Google, Pubmed and Scopus were used as search engines for the articles.

CONCLUSIONS

This review offers an overview on the pathophysiology of obesity, management of obesity-linked disorders, and areas in need of attention such as novel therapeutic adipokines and their possible future perspectives as therapeutic agents.

Kinetic Mathematical Modeling of Oxidative Phosphorylation in Cardiomyocyte Mitochondria

Oxidative phosphorylation (OXPHOS) is an oxygen-dependent process that consumes catabolized nutrients to produce adenosine triphosphate (ATP) to drive energy-dependent biological processes such as excitation-contraction coupling in cardiomyocytes. In addition to in vivo and in vitro experiments, in silico models are valuable for investigating the underlying mechanisms of OXPHOS and predicting its consequences in both physiological and pathological conditions. Here, we compare several prominent kinetic models of OXPHOS in cardiomyocytes. We examine how their mathematical expressions were derived, how their parameters were obtained, the conditions of their experimental counterparts, and the predictions they generated. We aim to explore the general landscape of energy production mechanisms in cardiomyocytes for future in silico models.

Quantitative modeling of pentose phosphate pathway response to oxidative stress reveals a cooperative regulatory strategy

Living cells use signaling and regulatory mechanisms to adapt to environmental stresses. Adaptation to oxidative stress involves the regulation of many enzymes in both glycolysis and pentose phosphate pathways (PPP), so as to support PPP-driven NADPH recycling for antioxidant defense. The underlying regulatory logic is investigated by developing a kinetic modeling approach fueled with metabolomics and 13C-fluxomics datasets from human fibroblast cells. Bayesian parameter estimation and phenotypic analysis of models highlight complementary roles for several metabolite-enzyme regulations. Specifically, carbon flux rerouting into PPP involves a tight coordination between the upregulation of G6PD activity concomitant to a decreased NADPH/NADP+ ratio and the differential control of downward and upward glycolytic fluxes through the joint inhibition of PGI and GAPD enzymes. Such functional interplay between distinct regulatory feedbacks promotes efficient detoxification and homeostasis response over a broad range of stress level, but can also explain paradoxical pertubation phenotypes for instance reported for 6PGD modulation in mammalian cells.

Digitoxin Affects Metabolism, ROS Production and Proliferation in Pancreatic Cancer Cells Differently Depending on the Cell Phenotype

Digitoxin has repeatedly shown to have negative effects on cancer cell viability; however, the actual mechanism is still unknown. In this study, we investigated the effects of digitoxin (1-100 nM) in four pancreatic cancer cell lines, BxPC-3, CFPAC-1, Panc-1, and AsPC-1. The cell lines differ in their KRAS/BRAF mutational status and primary tumor or metastasis origin. We could detect differences in the basal rates of cell proliferation, glycolysis, and ROS production, giving the cell lines different phenotypes. Digitoxin treatment induced apoptosis in all four cell lines, but to different degrees. Cells derived from primary tumors (Panc-1 and BxPC-3) were highly proliferating with a high proportion of cells in the S/G2 phase, and were more sensitive to digitoxin treatment than the cell lines derived from metastases (CFPAC-1 and AsPC-1), with a high proportion of cells in G0/G1. In addition, the effects of digitoxin on the rate of glycolysis, ROS production, and proliferation were dependent on the basal metabolism and origin of the cells. The KRAS downstream signaling pathways were not altered by digitoxin treatment, thus the effects exerted by digitoxin were probably disconnected from these signaling pathways. We conclude that digitoxin is a promising treatment in highly proliferating pancreatic tumors.

Mechanisms and mathematical modelling of ROS production by the mitochondrial electron transport chain

Reactive oxygen species (ROS) are recognised both as damaging molecules and intracellular signalling entities. In addition to its role in ATP generation, the mitochondrial electron transport chain (ETC) constitutes a relevant source of mitochondrial ROS, in particular during pathological conditions. Mitochondrial ROS homeostasis depends on species- and site-dependent ROS production, their bioreactivity, diffusion, and scavenging. However, our quantitative understanding of mitochondrial ROS homeostasis has thus far been hampered by technical limitations, including lack of truly site- and/or ROS-specific reporter molecules. In this context, the use of computational models is of great value to complement and interpret empirical data, as well as to predict variables that are difficult to assess experimentally. During the last decades, various mechanistic models of ETC-mediated ROS production have been developed. Although these often-complex models have generated novel insights, their parameterisation, analysis, and integration with other computational models is not straightforward. In contrast, phenomenological (sometimes termed "minimal") models use a relatively small set of equations to describe empirical relationship(s) between ROS-related and other parameters, and generally aim to explore system behaviour and generate hypotheses for experimental validation. In this review, we first discuss ETC-linked ROS homeostasis and introduce various detailed mechanistic models. Next, we present how bioenergetic parameters (e.g. NADH/NAD⁺ ratio, mitochondrial membrane potential) relate to site-specific ROS production within the ETC and how these relationships can be used to design minimal models of ROS homeostasis. Finally, we illustrate how minimal models have been applied to explore pathophysiological aspects of ROS.

The Modulatory Role of sti-1 in Methylmercury-Induced Toxicity in Caenorhabditis elegans

Human exposure to the neurotoxin methylmercury (MeHg) poses a significant health risk to the development of the nervous system. The mechanisms of MeHg-induced neurotoxicity are associated with the disruption of cellular homeostasis, and include oxidative stress, loss of calcium homeostasis, and impaired protein quality control. The stress inducible protein 1 (STI-1) is involved in the regulation of protein quality control by acting as a protein cochaperone to maintain optimal protein unfolding and refolding. Here, we utilized the Caenorhabditis elegans (C. elegans) model of MeHg toxicity to characterize the role of the sti-1 gene in MeHg-induced toxicity. We showed that lifespan and developmental milestone timings were significantly altered in sti-1 knockout (KO) animals with MeHg exposure. However, knocking down sti-1 by RNAi did not result in an analogous effect for lifespan, but did still sensitize to delays in developmental milestone progression by acute MeHg, suggesting that insufficiency of sti-1 does not recapitulate all phenotypes of the null mutation. Furthermore, inhibition of ATP levels by MeHg exposure was modulated by sti-1. Considering that the skn-1/gst-4 pathway is highly involved in metal's toxicity, such pathway was also explored in our model. We showed that sti-1 mutant worms exhibited impaired capacity to upregulate the antioxidant genes skn-1/gst-4, highlighting a central role of sti-1 in modulating antioxidant response. Lastly, we showed that loss-of-function mutation in the rrf-3 gene, which encodes a putative RNA-directed RNA polymerase, has significant effect in altering MeHg-induced toxicity by potentiating the animal's detoxification system. Altogether, our novel data show an indispensable role of protein quality control in the defense against MeHg toxicity.

Computational modeling of mitochondrial K+- and H+-driven ATP synthesis

ATP synthase (F1Fo) is a rotary molecular engine that harnesses energy from electrochemical-gradients across the inner mitochondrial membrane for ATP synthesis. Despite the accepted tenet that F1Fo transports exclusively H+, our laboratory has demonstrated that, in addition to H+, F1Fo ATP synthase transports a significant fraction of ΔΨm-driven charge as K+ to synthesize ATP. Herein, we utilize a computational modeling approach as a proof of principle of the feasibility of the core mechanism underlying the enhanced ATP synthesis, and to explore its bioenergetic consequences. A minimal model comprising the 'core' mechanism constituted by ATP synthase, driven by both proton (PMF) and potassium motive force (KMF), respiratory chain, adenine nucleotide translocator, Pi carrier, and K+/H+ exchanger (KHEmito) was able to simulate enhanced ATP synthesis and respiratory fluxes determined experimentally with isolated heart mitochondria. This capacity of F1Fo ATP synthase confers mitochondria with a significant energetic advantage compared to K+ transport through a channel not linked to oxidative phosphorylation (OxPhos). The K+-cycling mechanism requires a KHEmito that exchanges matrix K+ for intermembrane space H+, leaving PMF as the overall driving energy of OxPhos, in full agreement with the standard chemiosmotic mechanism. Experimental data of state 4➔3 energetic transitions, mimicking low to high energy demand, could be reproduced by an integrated computational model of mitochondrial function that incorporates the 'core' mechanism. Model simulations display similar behavior compared to the experimentally observed changes in ΔΨm, mitochondrial K+ uptake, matrix volume, respiration, and ATP synthesis during the energetic transitions at physiological pH and K+ concentration. The model also explores the role played by KHEmito in modulating the energetic performance of mitochondria. The results obtained support the available experimental evidence on ATP synthesis driven by K+ and H+ transport through the F1Fo ATP synthase.

Mitochondrial Ca2+ Homeostasis: Emerging Roles and Clinical Significance in Cardiac Remodeling

Mitochondria are the sites of oxidative metabolism in eukaryotes where the metabolites of sugars, fats, and amino acids are oxidized to harvest energy. Notably, mitochondria store Ca²⁺ and work in synergy with organelles such as the endoplasmic reticulum and extracellular matrix to control the dynamic balance of Ca²⁺ concentration in cells. Mitochondria are the vital organelles in heart tissue. Mitochondrial Ca²⁺ homeostasis is particularly important for maintaining the physiological and pathological mechanisms of the heart. Mitochondrial Ca²⁺ homeostasis plays a key role in the regulation of cardiac energy metabolism, mechanisms of death, oxygen free radical production, and autophagy. The imbalance of mitochondrial Ca²⁺ balance is closely associated with cardiac remodeling. The mitochondrial Ca²⁺ uniporter (mtCU) protein complex is responsible for the uptake and release of mitochondrial Ca²⁺ and regulation of Ca²⁺ homeostasis in mitochondria and consequently, in cells. This review summarizes the mechanisms of mitochondrial Ca²⁺ homeostasis in physiological and pathological cardiac remodeling and the regulatory effects of the mitochondrial calcium regulatory complex on cardiac energy metabolism, cell death, and autophagy, and also provides the theoretical basis for mitochondrial Ca²⁺ as a novel target for the treatment of cardiovascular diseases.

Understanding mitochondria and the utility of optimization as a canonical framework for identifying and modeling mitochondrial pathways

The goal of this paper is to provide an overview of our current understanding of mitochondrial function as a framework to motivate the hypothesis that mitochondrial behavior is governed by optimization principles that are constrained by the laws of the physical and biological sciences. Then, mathematical optimization tools can generally be useful to model some of these processes under reasonable assumptions and limitations. We are specifically interested in optimizations via variational methods, which are briefly summarized. Within such an optimization framework, we suggest that the numerous mechanical instigators of cell and intracellular functioning can be modeled utilizing some of the principles of mechanics that govern engineered systems, as well as by the frequently observed feedback and feedforward mechanisms that coordinate the multitude of processes within cells. These mechanical aspects would need to be coupled to governing biochemical rules. Of course, biological systems are significantly more complex than engineered systems, and require considerably more experimentation to ascertain and characterize parameters and subsequent behavior. That complexity requires well-defined limitations and assumptions for any derived models. Optimality is being motivated as a framework to help us understand how cellular decisions are made, especially those that transition between physiological behaviors and dysfunctions along pathophysiological pathways. We elaborate on our interpretation of optimality and cellular decision making within the body of this paper, as we revisit these ideas in the numerous different contexts of mitochondrial functions.

Computational Approaches and Tools as Applied to the Study of Rhythms and Chaos in Biology

The temporal dynamics in biological systems displays a wide range of behaviors, from periodic oscillations, as in rhythms, bursts, long-range (fractal) correlations, chaotic dynamics up to brown and white noise. Herein, we propose a comprehensive analytical strategy for identifying, representing, and analyzing biological time series, focusing on two strongly linked dynamics: periodic (oscillatory) rhythms and chaos. Understanding the underlying temporal dynamics of a system is of fundamental importance; however, it presents methodological challenges due to intrinsic characteristics, among them the presence of noise or trends, and distinct dynamics at different time scales given by molecular, dcellular, organ, and organism levels of organization. For example, in locomotion circadian and ultradian rhythms coexist with fractal dynamics at faster time scales. We propose and describe the use of a combined approach employing different analytical methodologies to synergize their strengths and mitigate their weaknesses. Specifically, we describe advantages and caveats to consider for applying probability distribution, autocorrelation analysis, phase space reconstruction, Lyapunov exponent estimation as well as different analyses such as harmonic, namely, power spectrum; continuous wavelet transforms; synchrosqueezing transform; and wavelet coherence. Computational harmonic analysis is proposed as an analytical framework for using different types of wavelet analyses. We show that when the correct wavelet analysis is applied, the complexity in the statistical properties, including temporal scales, present in time series of signals, can be unveiled and modeled. Our chapter showcase two specific examples where an in-depth analysis of rhythms and chaos is performed: (1) locomotor and food intake rhythms over a 42-day period of mice subjected to different feeding regimes; and (2) chaotic calcium dynamics in a computational model of mitochondrial function.

Loss of Mitochondrial Ca2+ Uniporter Limits Inotropic Reserve and Provides Trigger and Substrate for Arrhythmias in Barth Syndrome Cardiomyopathy

BACKGROUND

Barth syndrome (BTHS) is caused by mutations of the gene encoding tafazzin, which catalyzes maturation of mitochondrial cardiolipin and often manifests with systolic dysfunction during early infancy. Beyond the first months of life, BTHS cardiomyopathy typically transitions to a phenotype of diastolic dysfunction with preserved ejection fraction, blunted contractile reserve during exercise, and arrhythmic vulnerability. Previous studies traced BTHS cardiomyopathy to mitochondrial formation of reactive oxygen species (ROS). Because mitochondrial function and ROS formation are regulated by excitation-contraction coupling, integrated analysis of mechano-energetic coupling is required to delineate the pathomechanisms of BTHS cardiomyopathy.

METHODS

We analyzed cardiac function and structure in a mouse model with global knockdown of tafazzin (Taz-KD) compared with wild-type littermates. Respiratory chain assembly and function, ROS emission, and Ca²⁺ uptake were determined in isolated mitochondria. Excitation-contraction coupling was integrated with mitochondrial redox state, ROS, and Ca²⁺ uptake in isolated, unloaded or preloaded cardiac myocytes, and cardiac hemodynamics analyzed in vivo.

RESULTS

Taz-KD mice develop heart failure with preserved ejection fraction (>50%) and age-dependent progression of diastolic dysfunction in the absence of fibrosis. Increased myofilament Ca²⁺ affinity and slowed cross-bridge cycling caused diastolic dysfunction, in part, compensated by accelerated diastolic Ca²⁺ decay through preactivated sarcoplasmic reticulum Ca²⁺-ATPase. Taz deficiency provoked heart-specific loss of mitochondrial Ca²⁺ uniporter protein that prevented Ca²⁺-induced activation of the Krebs cycle during β-adrenergic stimulation, oxidizing pyridine nucleotides and triggering arrhythmias in cardiac myocytes. In vivo, Taz-KD mice displayed prolonged QRS duration as a substrate for arrhythmias, and a lack of inotropic response to β-adrenergic stimulation. Cellular arrhythmias and QRS prolongation, but not the defective inotropic reserve, were restored by inhibiting Ca²⁺ export through the mitochondrial Na⁺/Ca²⁺ exchanger. All alterations occurred in the absence of excess mitochondrial ROS in vitro or in vivo.

CONCLUSIONS

Downregulation of mitochondrial Ca²⁺ uniporter, increased myofilament Ca²⁺ affinity, and preactivated sarcoplasmic reticulum Ca²⁺-ATPase provoke mechano-energetic uncoupling that explains diastolic dysfunction and the lack of inotropic reserve in BTHS cardiomyopathy. Furthermore, defective mitochondrial Ca²⁺ uptake provides a trigger and a substrate for ventricular arrhythmias. These insights can guide the ongoing search for a cure of this orphaned disease.

Brain energetics, mitochondria, and traumatic brain injury

We review current thinking about, and draw connections between, brain energetics and metabolism, and between mitochondria and traumatic brain injury. Energy is fundamental to proper brain function. Its creation in a useful form for neurons and glia, and consistently in response to the brain's high energy needs, is critical for physiological pathways. Dysfunction in the mechanisms of energy production is at the center of neurological and neuropsychiatric pathologies. We examine the connections between energetics and mitochondria - the organelle responsible for almost all the energy production in the cell - and how secondary pathologies in traumatic brain injury result from energetic dysfunction. This paper interweaves these topics, a necessity since they are closely coupled, and identifies where there exist a lack of understanding and of data. In addition to summarizing current thinking in these disciplines, our goal is to suggest a framework for the mathematical modeling of mechanisms and pathways based on optimal energetic decisions.

Identifying Site-Specific Superoxide and Hydrogen Peroxide Production Rates From the Mitochondrial Electron Transport System Using a Computational Strategy

Mitochondrial reactive oxygen species (ROS) play important roles in cellular signaling; however, certain pathological conditions such as ischemia/reperfusion (I/R) injury disrupt ROS homeostasis and contribute to cell death. A major impediment to developing therapeutic measures against oxidative stress-induced cellular damage is the lack of a quantitative framework to identify the specific sources and regulatory mechanisms of mitochondrial ROS production. We developed a thermodynamically consistent, mass-and-charge balanced, kinetic model of mitochondrial ROS homeostasis focused on redox sites of electron transport chain complexes I, II, and III. The model was calibrated and corroborated using comprehensive data sets relevant to ROS homeostasis. The model predicts that complex I ROS production dominates other sources under conditions favoring a high membrane potential with elevated nicotinamide adenine dinucleotide (NADH) and ubiquinol (QH₂) levels. In general, complex I contributes to significant levels of ROS production under pathological conditions, while complexes II and III are responsible for basal levels of ROS production, especially when QH₂ levels are elevated. The model also reveals that hydrogen peroxide production by complex I underlies the non-linear relationship between ROS emission and O₂ at low O₂ concentrations. Lastly, the model highlights the need to quantify scavenging system activity under different conditions to establish a complete picture of mitochondrial ROS homeostasis. In summary, we describe the individual contributions of the electron transport system complex redox sites to total ROS emission in mitochondria respiring under various combinations of NADH- and Q-linked respiratory fuels under varying workloads.

Mitochondrial remodelling-a vicious cycle in diabetic complications

Diabetes mellitus (DM) is a chronic, metabolic condition characterized by excessive blood glucose that causes perturbations in physiological functioning of almost all the organs of human body. This devastating metabolic disease has its implications in cognitive decline, heart damage, renal, retinal and neuronal complications that severely affects quality of life and associated with decreased life expectancy. Mitochondria possess adaptive mechanisms to meet the cellular energy demand and combat cellular stress. In recent years mitochondrial homeostasis has been point of focus where several mechanisms regulating mitochondrial health and function are evaluated. Mitochondrial dynamics plays crucial role in maintaining healthy mitochondria in cell under physiological as well as stress condition. Mitochondrial dynamics and corresponding regulating mechanisms have been implicated in progression of metabolic disorders including diabetes and its complications. In current review we have discussed about role of mitochondrial dynamics under physiological and pathological conditions. Also, modulation of mitochondrial fission and fusion in diabetic complications are described. The available literature supports mitochondrial remodelling as reliable target for diabetic complications.

Oxidative Stress Evaluation in Ischemia Reperfusion Models: Characteristics, Limits and Perspectives

Ischemia reperfusion injury is a complex process consisting of a seemingly chaotic but actually organized and compartmentalized shutdown of cell function, of which oxidative stress is a key component. Studying oxidative stress, which results in an imbalance between reactive oxygen species (ROS) production and antioxidant defense activity, is a multi-faceted issue, particularly considering the double function of ROS, assuming roles as physiological intracellular signals and as mediators of cellular component damage. Herein, we propose a comprehensive overview of the tools available to explore oxidative stress, particularly in the study of ischemia reperfusion. Applying chemistry as well as biology, we present the different models currently developed to study oxidative stress, spanning the vitro and the silico, discussing the advantages and the drawbacks of each set-up, including the issues relating to the use of in vitro hypoxia as a surrogate for ischemia. Having identified the limitations of historical models, we shall study new paradigms, including the use of stem cell-derived organoids, as a bridge between the in vitro and the in vivo comprising 3D intercellular interactions in vivo and versatile pathway investigations in vitro. We shall conclude this review by distancing ourselves from "wet" biology and reviewing the in silico, computer-based, mathematical modeling, and numerical simulation options: (a) molecular modeling with quantum chemistry and molecular dynamic algorithms, which facilitates the study of molecule-to-molecule interactions, and the integration of a compound in a dynamic environment (the plasma membrane...); (b) integrative systemic models, which can include many facets of complex mechanisms such as oxidative stress or ischemia reperfusion and help to formulate integrated predictions and to enhance understanding of dynamic interaction between pathways.

Mitochondrial Ca2+, redox environment and ROS emission in heart failure: Two sides of the same coin?

Heart failure (HF) is a progressive, debilitating condition characterized, in part, by altered ionic equilibria, increased ROS production and impaired cellular energy metabolism, contributing to variable profiles of systolic and diastolic dysfunction with significant functional limitations and risk of premature death. We summarize current knowledge concerning changes of intracellular Na+ and Ca2+ control mechanisms during the disease progression and their consequences on mitochondrial Ca2+ homeostasis and the shift in redox balance. Absent existing biological data, our computational modeling studies advance a new 'in silico' analysis to reconcile existing opposing views, based on different experimental HF models, regarding variations in mitochondrial Ca2+ concentration that participate in triggering and perpetuating oxidative stress in the failing heart and their impact on cardiac energetics. In agreement with our hypothesis and the literature, model simulations demonstrate the possibility that the heart's redox status together with cytoplasmic Na+ concentrations act as regulators of mitochondrial Ca2+ levels in HF and of the bioenergetics response that will ultimately drive ATP supply and oxidative stress. The resulting model predictions propose future directions to study the evolution of HF as well as other types of heart disease, and to develop novel testable mechanistic hypotheses that may lead to improved therapeutics.

Introduction to In Silico Modeling to Study ROS Dynamics

This book chapter is drafted for biologists with experimental experiences in ROS biology but being newcomers in the field of modeling. We start with a general introduction about computational modeling in biology and an overview of software tools suitable for beginners. This chapter encompasses an introduction to computational models with special focus on simulation of ROS dynamics. A step-by-step tutorial follows providing guidance for all relevant model development processes. This course of action gives a comprehensible way to understand the benefits of computational models and to gain the necessary knowledge to build own small equation-based models. Small models can be created without any special programming expertise or in-depth technical and mathematical knowledge. Afterward in the final section, a short overview of pitfalls, challenges, and limitations is provided, combined with suggestions for further reading to improve and expand modeling skills of biologists.

Hydrogen peroxide diffusion and scavenging shapes mitochondrial network instability and failure by sensitizing ROS-induced ROS release

The mitochondrial network of cardiac cells is finely tuned for ATP delivery to sites of energy demand; however, emergent phenomena, such as mitochondrial transmembrane potential oscillations or propagating waves of depolarization have been observed under metabolic stress. While regenerative signaling by reactive oxygen species (ROS)-induced ROS release (RIRR) has been suggested as a potential trigger, it is unknown how it could lead to widespread responses. Here, we present a novel computational model of RIRR transmission that explains the mechanisms of this phenomenon. The results reveal that superoxide mediates neighbor-neighbor activation of energy-dissipating ion channels, while hydrogen peroxide distributes oxidative stress to sensitize the network to mitochondrial criticality. The findings demonstrate the feasibility of RIRR as a synchronizing factor across the dimensions of the adult heart cell and illustrate how a cascade of failures at the organellar level can scale to impact cell and organ level functions of the heart.

Kinetic modeling of H2O2 dynamics in the mitochondria of HeLa cells

Hydrogen peroxide (H2O2) promotes a range of phenotypes depending on its intracellular concentration and dosing kinetics, including cell death. While this qualitative relationship has been well established, the quantitative and mechanistic aspects of H2O2 signaling are still being elucidated. Mitochondria, a putative source of intracellular H2O2, have recently been demonstrated to be particularly vulnerable to localized H2O2 perturbations, eliciting a dramatic cell death response in comparison to similar cytosolic perturbations. We sought to improve our dynamic and mechanistic understanding of the mitochondrial H2O2 reaction network in HeLa cells by creating a kinetic model of this system and using it to explore basal and perturbed conditions. The model uses the most current quantitative proteomic and kinetic data available to predict reaction rates and steady-state concentrations of H2O2 and its reaction partners within individual mitochondria. Time scales ranging from milliseconds to one hour were simulated. We predict that basal, steady-state mitochondrial H2O2 will be in the low nM range (2-4 nM) and will be inversely dependent on the total pool of peroxiredoxin-3 (Prx3). Neglecting efflux of H2O2 to the cytosol, the mitochondrial reaction network is expected to control perturbations well up to H2O2 generation rates ~50 μM/s (0.25 nmol/mg-protein/s), above which point the Prx3 system would be expected to collapse. Comparison of these results with redox Western blots of Prx3 and Prx2 oxidation states demonstrated reasonable trend agreement at short times (≤ 15 min) for a range of experimentally perturbed H2O2 generation rates. At longer times, substantial efflux of H2O2 from the mitochondria to the cytosol was evidenced by peroxiredoxin-2 (Prx2) oxidation, and Prx3 collapse was not observed. A refined model using Monte Carlo parameter sampling was used to explore rates of H2O2 efflux that could reconcile model predictions of Prx3 oxidation states with the experimental observations.

Computationally modeling mammalian succinate dehydrogenase kinetics identifies the origins and primary determinants of ROS production

Succinate dehydrogenase (SDH) is an inner mitochondrial membrane protein complex that links the Krebs cycle to the electron transport system. It can produce significant amounts of superoxide ([Formula: see text]) and hydrogen peroxide (H₂O₂); however, the precise mechanisms are unknown. This fact hinders the development of next-generation antioxidant therapies targeting mitochondria. To help address this problem, we developed a computational model to analyze and identify the kinetic mechanism of [Formula: see text] and H₂O₂ production by SDH. Our model includes the major redox centers in the complex, namely FAD, three iron-sulfur clusters, and a transiently bound semiquinone. Oxidation state transitions involve a one- or two-electron redox reaction, each being thermodynamically constrained. Model parameters were simultaneously fit to many data sets using a variety of succinate oxidation and free radical production data. In the absence of respiratory chain inhibitors, model analysis revealed the 3Fe-4S iron-sulfur cluster as the primary [Formula: see text] source. However, when the quinone reductase site is inhibited or the quinone pool is highly reduced, [Formula: see text] is generated primarily by the FAD. In addition, H₂O₂ production is only significant when the enzyme is fully reduced, and fumarate is absent. Our simulations also reveal that the redox state of the quinone pool is the primary determinant of free radical production by SDH. In this study, we showed the importance of analyzing enzyme kinetics and associated side reactions in a consistent, quantitative, and biophysically detailed manner using a diverse set of experimental data to interpret and explain experimental observations from a unified perspective.

Integrating Ultra-Weak Photon Emission Analysis in Mitochondrial Research

Once regarded solely as the energy source of the cell, nowadays mitochondria are recognized to perform multiple essential functions in addition to energy production. Since the discovery of pathogenic mitochondrial DNA defects in the 1980s, research advances have revealed an increasing number of common human diseases, which share an underlying pathogenesis involving mitochondrial dysfunction. A major factor in this dysfunction is reactive oxygen species (ROS), which influence the mitochondrial-nuclear crosstalk and the link with the epigenome, an influence that provides explanations for pathogenic mechanisms. Regarding these mechanisms, we should take into account that mitochondria produce the majority of ultra-weak photon emission (UPE), an aspect that is often ignored - this type of emission may serve as assay for ROS, thus providing new opportunities for a non-invasive diagnosis of mitochondrial dysfunction. In this article, we overviewed three relevant areas of mitochondria-related research over the period 1960-2020: (a) respiration and energy production, (b) respiration-related production of free radicals and other ROS species, and (c) ultra-weak photon emission in relation to ROS and stress. First, we have outlined how these research areas initially developed independently of each other - following that, our review aims to show their stepwise integration during later stages of development. It is suggested that a further stimulation of research on UPE may have the potential to enhance the progress of modern mitochondrial research and its integration in medicine.

Compromised DNA repair is responsible for diabetes-associated fibrosis

Diabetes-associated organ fibrosis, marked by elevated cellular senescence, is a growing health concern. Intriguingly, the mechanism underlying this association remained unknown. Moreover, insulin alone can neither reverse organ fibrosis nor the associated secretory phenotype, favoring the exciting notion that thus far unknown mechanisms must be operative. Here, we show that experimental type 1 and type 2 diabetes impairs DNA repair, leading to senescence, inflammatory phenotypes, and ultimately fibrosis. Carbohydrates were found to trigger this cascade by decreasing the NAD⁺ /NADH ratio and NHEJ-repair in vitro and in diabetes mouse models. Restoring DNA repair by nuclear over-expression of phosphomimetic RAGE reduces DNA damage, inflammation, and fibrosis, thereby restoring organ function. Our study provides a novel conceptual framework for understanding diabetic fibrosis on the basis of persistent DNA damage signaling and points to unprecedented approaches to restore DNA repair capacity for resolution of fibrosis in patients with diabetes.

Metabolic remodelling of glucose, fatty acid and redox pathways in the heart of type 2 diabetic mice

KEY POINTS

Hearts from type 2 diabetic animals display perturbations in excitation-contraction coupling, impairing myocyte contractility and delaying relaxation, along with altered substrate consumption patterns. Under high glucose and β-adrenergic stimulation conditions, palmitate can, at least in part, offset left ventricle (LV) dysfunction in hearts from diabetic mice, improving contractility and relaxation while restoring coronary perfusion pressure. Fluxome calculations of central catabolism in diabetic hearts show that, in the presence of palmitate, there is a metabolic remodelling involving tricarboxylic acid cycle, polyol and pentose phosphate pathways, leading to improved redox balance in cytoplasmic and mitochondrial compartments. Under high glucose and increased energy demand, the metabolic/fluxomic redirection leading to restored redox balance imparted by palmitate helps explain maintained LV function and may contribute to designing novel therapeutic approaches to prevent cardiac dysfunction in diabetic patients.

ABSTRACT

Type-2 diabetes (T2DM) leads to reduced myocardial performance, and eventually heart failure. Excessive accumulation of lipids and glucose is central to T2DM cardiomyopathy. Previous data showed that palmitate (Palm) or glutathione preserved heart mitochondrial energy/redox balance under excess glucose, rescuing β-adrenergic-stimulated cardiac excitation-contraction coupling. However, the mechanisms underlying the accompanying improved contractile performance have been largely ignored. Herein we explore in intact heart under substrate excess the metabolic remodelling associated with cardiac function in diabetic db/db mice subjected to stress given by β-adrenergic stimulation with isoproterenol and high glucose compared to their non-diabetic controls (+/+, WT) under euglycaemic conditions. When perfused with Palm, T2DM hearts exhibited improved contractility/relaxation compared to WT, accompanied by extensive metabolic remodelling as demonstrated by metabolomics-fluxomics combined with bioinformatics and computational modelling. The T2DM heart metabolome showed significant differences in the abundance of metabolites in pathways related to glucose, lipids and redox metabolism. Using a validated computational model of heart's central catabolism, comprising glucose and fatty acid (FA) oxidation in cytoplasmic and mitochondrial compartments, we estimated that fluxes through glucose degradation pathways are ∼2-fold lower in heart from T2DM vs. WT under all conditions studied. Palm addition elicits improvement of the redox status via enhanced β-oxidation and decreased glucose uptake, leading to flux-redirection away from redox-consuming pathways (e.g. polyol) while maintaining the flux through redox-generating pathways together with glucose-FA 'shared fuelling' of oxidative phosphorylation. Thus, available FAs such as Palm may help improve function via enhanced redox balance in T2DM hearts during peaks of hyperglycaemia and increased workload.

Protein S-glutathionylation reactions as a global inhibitor of cell metabolism for the desensitization of hydrogen peroxide signals

The pathogenesis of many human diseases has been attributed to the over production of reactive oxygen species (ROS), particularly superoxide (O₂^●-) and hydrogen peroxide (H₂O₂), by-products of metabolism that are generated by the premature reaction of electrons with molecular oxygen (O₂) before they reach complex IV of the respiratory chain. To date, there are 32 known ROS generators in mammalian cells, 16 of which reside inside mitochondria. Importantly, although these ROS are deleterious at high levels, controlled and temporary bursts in H₂O₂ production is beneficial to mammalian cells. Mammalian cells use sophisticated systems to take advantage of the second messaging properties of H₂O₂. This includes controlling its availability using antioxidant systems and negative feedback loops that inhibit the genesis of ROS at sites of production. At its core, ROS production depends on fuel metabolism. Therefore, desensitizing H₂O₂ signals would also require the temporary inhibition of fuel combustion and fluxes through metabolic pathways that promote ROS production. Additionally, this would also demand the diversion of fuels and nutrients into pathways that generate NADPH and other molecules required to maintain cellular redox buffering capacity. Therefore, fuel selection and metabolic flux plays an integral role in dictating the strength and duration of cellular redox signals. In the present review I provide an updated view on the function of protein S-glutathionylation, a ubiquitous redox sensitive modification involving the formation of a disulfide between the small molecular antioxidant glutathione and a cysteine residue, in the regulation of cellular metabolism on a global scale. To date, these concepts have mostly been reviewed at the level of mitochondrial bioenergetics in the contexts of health and disease. Careful examination of the literature revealed that glutathionylation is a temporary inhibitor of most metabolic pathways including glycolysis, the Krebs cycle, oxidative phosphorylation, amino acid metabolism, and fatty acid combustion, resulting in the diversion of fuels towards NADPH-producing pathways and the inhibition of ROS production. Armed with this information, I propose that protein S-glutathionylation reactions desensitize H₂O₂ signals emanating from catabolic pathways using a three-pronged regulatory mechanism; 1) inhibition of metabolic flux through pathways that promote ROS production, 2) diversion of metabolites towards pathways that support antioxidant defenses, and 3) direct inhibition of ROS-generating enzymes.

Modelling mitochondrial ROS production by the respiratory chain

ROS (superoxide and oxygen peroxide in this paper) play a dual role as signalling molecules and strong oxidizing agents leading to oxidative stress. Their production mainly occurs in mitochondria although they may have other locations (such as NADPH oxidase in particular cell types). Mitochondrial ROS production depends in an interweaving way upon many factors such as the membrane potential, the cell type and the respiratory substrates. Moreover, it is experimentally difficult to quantitatively assess the contribution of each potential site in the respiratory chain. To overcome these difficulties, mathematical models have been developed with different degrees of complexity in order to analyse different physiological questions ranging from a simple reproduction/simulation of experimental results to a detailed model of the possible mechanisms leading to ROS production. Here, we analyse experimental results concerning ROS production including results still under discussion. We then critically review the three models of ROS production in the whole respiratory chain available in the literature and propose some direction for future modelling work.

ROS and pentose phosphate pathway: mathematical modelling of the metabolic regulation in response to xenobiotic-induced oxidative stress and the proposed Impact of the gluconate shunt

Elevated intracellular levels of reactive oxygen species (ROS), e.g. resulting from exposure to xenobiotics, can cause severe damages. Antioxidant defence mechanisms, which involve regulation of enzyme activities, protect cells to a certain extent. Nevertheless, continuous or increased exposure can overwhelm this system resulting in an adverse cellular state. To simulate exposure scenarios and to investigate the transition to an adverse cellular state, a mathematical model for the dynamics of ROS in response to xenobiotic-induced oxidative stress has been developed. It is based on exposure experiments of human urothelial cells (RT4) to the nitrated polycyclic aromatic hydrocarbon 3-nitrobenzanthrone (3-NBA), a component of diesel engine exhaust, and takes into account the following metabolic pathways of the antioxidant defence system: glutathione redox cycle scavenging directly ROS, the pentose phosphate pathway and the gluconate shunt as NADPH supplier and the beginning of glycolysis. In addition, ROS generation due to the bioactivation of 3-NBA has been implemented. The regulation of enzyme activities plays an important role in the presented mathematical model. The in silico model consists of ordinary differential equations on the basis of enzyme kinetics and mass action for the metabolism of 3-NBA. Parameters are either estimated from performed in vitro experiments via least-squares fitting or obtained from the literature. The results underline the importance of the pentose phosphate pathway to cope with oxidative stress and suggest an important role of the gluconate shunt during low-dose exposure.

Redox-Driven Signaling: 2-Oxo Acid Dehydrogenase Complexes as Sensors and Transmitters of Metabolic Imbalance

SIGNIFICANCE

This article develops a holistic view on production of reactive oxygen species (ROS) by 2-oxo acid dehydrogenase complexes. Recent Advances: Catalytic and structural properties of the complexes and their components evolved to minimize damaging effects of side reactions, including ROS generation, simultaneously exploiting the reactions for homeostatic signaling.

CRITICAL ISSUES

Side reactions of the complexes, characterized in vitro, are analyzed in view of protein interactions and conditions in vivo. Quantitative data support prevalence of the forward 2-oxo acid oxidation over the backward NADH oxidation in feeding physiologically significant ROS production by the complexes. Special focus on interactions between the active sites within 2-oxo acid dehydrogenase complexes highlights the central relevance of the complex-bound thiyl radicals in regulation of and signaling by complex-generated ROS. The thiyl radicals arise when dihydrolipoyl residues of the complexes regenerate FADH₂ from the flavin semiquinone coproduced with superoxide anion radical in 1e^- oxidation of FADH₂ by molecular oxygen.

FUTURE DIRECTIONS

Interaction of 2-oxo acid dehydrogenase complexes with thioredoxins (TRXs), peroxiredoxins, and glutaredoxins mediates scavenging of the thiyl radicals and ROS generated by the complexes, underlying signaling of disproportional availability of 2-oxo acids, CoA, and NAD⁺ in key metabolic branch points through thiol/disulfide exchange and medically important hypoxia-inducible factor, mammalian target of rapamycin (mTOR), poly (ADP-ribose) polymerase, and sirtuins. High reactivity of the coproduced ROS and thiyl radicals to iron/sulfur clusters and nitric oxide, peroxynitrite reductase activity of peroxiredoxins and transnitrosylating function of thioredoxin, implicate the side reactions of 2-oxo acid dehydrogenase complexes in nitric oxide-dependent signaling and damage.

Control and Regulation of Substrate Selection in Cytoplasmic and Mitochondrial Catabolic Networks. A Systems Biology Analysis

Appropriate substrate selection between fats and glucose is associated with the success of interventions that maintain health such as exercise or caloric restriction, or with the severity of diseases such as diabetes or other metabolic disorders. Although the interaction and mutual inhibition between glucose and fatty-acids (FAs) catabolism has been studied for decades, a quantitative and integrated understanding of the control and regulation of substrate selection through central catabolic pathways is lacking. We addressed this gap here using a computational model representing cardiomyocyte catabolism encompassing glucose (Glc) utilization, pyruvate transport into mitochondria and oxidation in the tricarboxylic acid (TCA) cycle, β-oxidation of palmitate (Palm), oxidative phosphorylation, ion transport, pH regulation, and ROS generation and scavenging in cytoplasmic and mitochondrial compartments. The model is described by 82 differential equations and 119 enzymatic, electron transport and substrate transport reactions accounting for regulatory mechanisms and key players, namely pyruvate dehydrogenase (PDH) and its modulation by multiple effectors. We applied metabolic control analysis to the network operating with various Glc to Palm ratios. The flux and metabolites’ concentration control were visualized through heat maps providing major insights into main control and regulatory nodes throughout the catabolic network. Metabolic pathways located in different compartments were found to reciprocally control each other. For example, glucose uptake and the ATP demand exert control on most processes in catabolism while TCA cycle activities and membrane-associated energy transduction reactions exerted control on mitochondrial processes namely β-oxidation. PFK and PDH, two highly regulated enzymes, exhibit opposite behavior from a control perspective. While PFK activity was a main rate-controlling step affecting the whole network, PDH played the role of a major regulator showing high sensitivity (elasticity) to substrate availability and key activators/inhibitors, a trait expected from a flexible substrate selector strategically located in the metabolic network. PDH regulated the rate of Glc and Palm consumption, consistent with its high sensitivity toward AcCoA, CoA, and NADH. Overall, these results indicate that the control of catabolism is highly distributed across the metabolic network suggesting that fuel selection between FAs and Glc goes well beyond the mechanisms traditionally postulated to explain the glucose-fatty-acid cycle.

Modelling the role of redox-related mechanisms in musculoskeletal ageing

The decline in the musculoskeletal system with age is driven at the cellular level by random molecular damage. Cells possess mechanisms to repair or remove damage and many of the pathways involved in this response are regulated by redox signals. However, with ageing there is an increase in oxidative stress which can lead to chronic inflammation and disruption of redox signalling pathways. The complexity of the processes involved has led to the use of computational modelling to help increase our understanding of the system, test hypotheses and make testable predictions. This paper will give a brief background of the biological systems that have been modelled, an introduction to computational modelling, a review of models that involve redox-related mechanisms that are applicable to musculoskeletal ageing, and finally a discussion of the future potential for modelling in this field.

A mechanistic model for thiol redox dynamics in the organogenesis stage rat conceptus

Precise control of the glutathione/glutathione disulfide (GSH/GSSG) redox balance is vital for the developing embryo, but regulatory mechanisms are poorly understood. We developed a novel, mechanistic mass-balance model for GSH metabolism in the organogenesis stage (gestational day 10.0-11.13) rat conceptus predicting the dynamics of 8 unique metabolites in 3 conceptal compartments: the visceral yolk sac (VYS), the extra-embryonic fluid (EEF) and the embryo proper (EMB). Our results show that thiol concentrations in all compartments are well predicted by the model. Protein synthesis is predicted to be a major efflux pathway for all amino acid precursors of GSH synthesis and an essential model element. Our model provides quantitative insights in the transport fluxes and enzymatic fluxes needed to maintain thiol redox balances under normal physiological conditions. This is crucial to further elucidate the mechanisms through which chemical exposure can perturb redox homeostasis, causing oxidative stress, and potentially birth defects.

Mitochondrial chaotic dynamics: Redox-energetic behavior at the edge of stability

Mitochondria serve multiple key cellular functions, including energy generation, redox balance, and regulation of apoptotic cell death, thus making a major impact on healthy and diseased states. Increasingly recognized is that biological network stability/instability can play critical roles in determining health and disease. We report for the first-time mitochondrial chaotic dynamics, characterizing the conditions leading from stability to chaos in this organelle. Using an experimentally validated computational model of mitochondrial function, we show that complex oscillatory dynamics in key metabolic variables, arising at the “edge” between fully functional and pathological behavior, sets the stage for chaos. Under these conditions, a mild, regular sinusoidal redox forcing perturbation triggers chaotic dynamics with main signature traits such as sensitivity to initial conditions, positive Lyapunov exponents, and strange attractors. At the “edge” mitochondrial chaos is exquisitely sensitive to the antioxidant capacity of matrix Mn superoxide dismutase as well as to the amplitude and frequency of the redox perturbation. These results have potential implications both for mitochondrial signaling determining health maintenance, and pathological transformation, including abnormal cardiac rhythms.

Systems Biology Approaches to Redox Metabolism in Stress and Disease States

SIGNIFICANCE

All cellular metabolic processes are tied to the cellular redox environment. Therefore, maintaining redox homeostasis is critically important for normal cell function. Indeed, redox stress contributes to the pathobiology of many human diseases. The cellular redox response system is composed of numerous interconnected components, including free radicals, redox couples, protein thiols, enzymes, metabolites, and transcription factors. Moreover, interactions between and among these factors are regulated in time and space. Owing to their complexity, systems biology approaches to the characterization of the cellular redox response system may provide insights into novel homeostatic mechanisms and methods of therapeutic reprogramming. Recent Advances: The emergence and development of systems biology has brought forth a set of innovative technologies that provide new avenues for studying redox metabolism. This article will review these systems biology approaches and their potential application to the study of redox metabolism in stress and disease states.

CRITICAL ISSUES

Clarifying the scope of biological intermediaries affected by dysregulated redox metabolism requires methods that are suitable for analyzing big datasets as classical methods that do not account for multiple interactions are unlikely to portray the totality of perturbed metabolic systems.

FUTURE DIRECTIONS

Given the diverse redox microenvironments within cells, it will be important to improve the spatial resolution of omic approaches. Futures studies on the integration of multiple systems-based methods and heterogeneous omics data for redox metabolism are required to accelerate the development of the field of redox systems biology. Antioxid. Redox Signal. 29, 953-972.

Quantitative systems models illuminate arrhythmia mechanisms in heart failure: Role of the Na+ -Ca2+ -Ca2+ /calmodulin-dependent protein kinase II-reactive oxygen species feedback

Quantitative systems modeling aims to integrate knowledge in different research areas with models describing biological mechanisms and dynamics to gain a better understanding of complex clinical syndromes. Heart failure (HF) is a chronic complex cardiac disease that results from structural or functional disorders impairing the ability of the ventricle to fill with or eject blood. Highly interactive and dynamic changes in mechanical, structural, neurohumoral, metabolic, and electrophysiological properties collectively predispose the failing heart to cardiac arrhythmias, which are responsible for about a half of HF deaths. Multiscale cardiac modeling and simulation integrate structural and functional data from HF experimental models and patients to improve our mechanistic understanding of this complex arrhythmia syndrome. In particular, they allow investigating how disease-induced remodeling alters the coupling of electrophysiology, Ca²⁺ and Na⁺ handling, contraction, and energetics that lead to rhythm derangements. The Ca²⁺ /calmodulin-dependent protein kinase II, which expression and activity are enhanced in HF, emerges as a critical hub that modulates the feedbacks between these various subsystems and promotes arrhythmogenesis. This article is categorized under: Physiology > Mammalian Physiology in Health and Disease Models of Systems Properties and Processes > Mechanistic Models Models of Systems Properties and Processes > Cellular Models Models of Systems Properties and Processes > Organ, Tissue, and Physiological Models.

Nicotinamide Improves Aspects of Healthspan, but Not Lifespan, in Mice

The role in longevity and healthspan of nicotinamide (NAM), the physiological precursor of NAD⁺, is elusive. Here, we report that chronic NAM supplementation improves healthspan measures in mice without extending lifespan. Untargeted metabolite profiling of the liver and metabolic flux analysis of liver-derived cells revealed NAM-mediated improvement in glucose homeostasis in mice on a high-fat diet (HFD) that was associated with reduced hepatic steatosis and inflammation concomitant with increased glycogen deposition and flux through the pentose phosphate and glycolytic pathways. Targeted NAD metabolome analysis in liver revealed depressed expression of NAM salvage in NAM-treated mice, an effect counteracted by higher expression of de novo NAD biosynthetic enzymes. Although neither hepatic NAD⁺ nor NADP⁺ was boosted by NAM, acetylation of some SIRT1 targets was enhanced by NAM supplementation in a diet- and NAM dose-dependent manner. Collectively, our results show health improvement in NAM-supplemented HFD-fed mice in the absence of survival effects.

Mitochondrial health, the epigenome and healthspan

Food nutrients and metabolic supply-demand dynamics constitute environmental factors that interact with our genome influencing health and disease states. These gene-environment interactions converge at the metabolic-epigenome-genome axis to regulate gene expression and phenotypic outcomes. Mounting evidence indicates that nutrients and lifestyle strongly influence genome-metabolic functional interactions determining disease via altered epigenetic regulation. The mitochondrial network is a central player of the metabolic-epigenome-genome axis, regulating the level of key metabolites [NAD(+), AcCoA (acetyl CoA), ATP] acting as substrates/cofactors for acetyl transferases, kinases (e.g. protein kinase A) and deacetylases (e.g. sirtuins, SIRTs). The chromatin, an assembly of DNA and nucleoproteins, regulates the transcriptional process, acting at the epigenomic interface between metabolism and the genome. Within this framework, we review existing evidence showing that preservation of mitochondrial network function is directly involved in decreasing the rate of damage accumulation thus slowing aging and improving healthspan.

Mitochondrial respiration and ROS emission during β-oxidation in the heart: An experimental-computational study

Lipids are main fuels for cellular energy and mitochondria their major oxidation site. Yet unknown is to what extent the fuel role of lipids is influenced by their uncoupling effects, and how this affects mitochondrial energetics, redox balance and the emission of reactive oxygen species (ROS). Employing a combined experimental-computational approach, we comparatively analyze β-oxidation of palmitoyl CoA (PCoA) in isolated heart mitochondria from Sham and streptozotocin (STZ)-induced type 1 diabetic (T1DM) guinea pigs (GPs). Parallel high throughput measurements of the rates of oxygen consumption (VO2) and hydrogen peroxide (H2O2) emission as a function of PCoA concentration, in the presence of L-carnitine and malate, were performed. We found that PCoA concentration < 200 nmol/mg mito protein resulted in low H2O2 emission flux, increasing thereafter in Sham and T1DM GPs under both states 4 and 3 respiration with diabetic mitochondria releasing higher amounts of ROS. Respiratory uncoupling and ROS excess occurred at PCoA > 600 nmol/mg mito prot, in both control and diabetic animals. Also, for the first time, we show that an integrated two compartment mitochondrial model of β-oxidation of long-chain fatty acids and main energy-redox processes is able to simulate the relationship between VO2 and H2O2 emission as a function of lipid concentration. Model and experimental results indicate that PCoA oxidation and its concentration-dependent uncoupling effect, together with a partial lipid-dependent decrease in the rate of superoxide generation, modulate H2O2 emission as a function of VO2. Results indicate that keeping low levels of intracellular lipid is crucial for mitochondria and cells to maintain ROS within physiological levels compatible with signaling and reliable energy supply.

Triosephosphates as Intermediates of the Crabtree Effect

An increase in glucose concentration in the medium rapidly decreases respiration rate in many cell types, including tumor cells. The molecular mechanism of this phenomenon, the Crabtree effect, is still unclear. It was shown earlier that adding the intermediate product of glycolysis fructose-1,6-bisphosphate to isolated mitochondria suppresses their respiration. To study possible roles of glycolytic intermediates in the Crabtree effect, we used a model organism, the yeast Saccharomyces cerevisiae. To have the option to rapidly increase intracellular concentrations of certain glycolytic intermediates, we used mutant cells with glycolysis blocked at different stages. We studied fast effects of glucose addition on the respiration rate in such cells. We found that addition of glucose affected cells with deleted phosphoglycerate mutase (strain gpm1-delta) more strongly than ones with inactivated aldolase or phosphofructokinase. In the case of preincubation of gpm1-delta cells with 2-deoxyglucose, which blocks glycolysis at the stage of 2-deoxyglucosephosphate formation, the effect of glucose addition was absent. This suggests that triosephosphates are intermediates of the Crabtree effect. Apart from this, the incubation of gpm1-delta cells in galactose-containing medium appeared to cause a large increase in their size. It was previously shown that galactose addition did not have any short-term effect on respiration rate of gpm1-delta cells and, at the same time, strongly suppressed their growth rate. Apparently, the influence of increasing triosephosphate concentration on yeast physiology is not limited to the activation of the Crabtree effect.

Mitochondrial energetics and calcium coupling in the heart

Contraction and relaxation of the heart consume large amounts of energy that need to be replenished by oxidative phosphorylation in mitochondria, and matching energy supply to demand involves the complimentary control of respiration through ADP and Ca²⁺ . In heart failure, an imbalance between ADP and Ca²⁺ leads to oxidation of mitochondrial pyridine nucleotides, where NADH oxidation may limit ATP production and contractile function, while NADPH oxidation can induce oxidative stress with consecutive maladaptive remodelling. Understanding the complex mechanisms that disturb this finely tuned equilibrium may aid the development of drugs that could ameliorate the progression of heart failure beyond the classical neuroendocrine inhibition.

Modelling the molecular mechanisms of aging

The aging process is driven at the cellular level by random molecular damage that slowly accumulates with age. Although cells possess mechanisms to repair or remove damage, they are not 100% efficient and their efficiency declines with age. There are many molecular mechanisms involved and exogenous factors such as stress also contribute to the aging process. The complexity of the aging process has stimulated the use of computational modelling in order to increase our understanding of the system, test hypotheses and make testable predictions. As many different mechanisms are involved, a wide range of models have been developed. This paper gives an overview of the types of models that have been developed, the range of tools used, modelling standards and discusses many specific examples of models that have been grouped according to the main mechanisms that they address. We conclude by discussing the opportunities and challenges for future modelling in this field.

Traditional and novel tools to probe the mitochondrial metabolism in health and disease

Mitochondrial metabolism links energy production to other essential cellular processes such as signaling, cellular differentiation, and apoptosis. In addition to producing adenosine triphosphate (ATP) as an energy source, mitochondria are responsible for the synthesis of a myriad of important metabolites and cofactors such as tetrahydrofolate, α-ketoacids, steroids, aminolevulinic acid, biotin, lipoic acid, acetyl-CoA, iron-sulfur clusters, heme, and ubiquinone. Furthermore, mitochondria and their metabolism have been implicated in aging and several human diseases, including inherited mitochondrial disorders, cardiac dysfunction, heart failure, neurodegenerative diseases, diabetes, and cancer. Therefore, there is great interest in understanding mitochondrial metabolism and the complex relationship it has with other cellular processes. A large number of studies on mitochondrial metabolism have been conducted in the last 50 years, taking a broad range of approaches. In this review, we summarize and discuss the most commonly used tools that have been used to study different aspects of the metabolism of mitochondria: ranging from dyes that monitor changes in the mitochondrial membrane potential and pharmacological tools to study respiration or ATP synthesis, to more modern tools such as genetically encoded biosensors and trans-omic approaches enabled by recent advances in mass spectrometry, computation, and other technologies. These tools have allowed the large number of studies that have shaped our current understanding of mitochondrial metabolism. WIREs Syst Biol Med 2017, 9:e1373. doi: 10.1002/wsbm.1373 For further resources related to this article, please visit the WIREs website.

Functional Implications of Cardiac Mitochondria Clustering

The spatio-temporal organization of mitochondria in cardiac myocytes facilitates myocyte-wide, cluster-bound, mitochondrial inner membrane potential oscillatory depolarizations, commonly triggered by metabolic or oxidative stressors. Local intermitochondrial coupling can be mediated by reactive oxygen species (ROS) that activate inner membrane pores to initiate a ROS-induced-ROS-release process that produces synchronized limit cycle oscillations of mitochondrial clusters within the whole mitochondrial network. The network's dynamic organization, structure and function can be assessed by quantifying dynamic local coupling constants and dynamic functional clustering coefficients, both providing information about the network's response to external stimuli. In addition to its special organization, the mitochondrial network of cardiac myocytes exhibits substrate-sensitive coupling constants and clustering coefficients. The myocyte's ability to form functional clusters of synchronously oscillating mitochondria is sensitive to conditions such as substrate availability (e.g., glucose, pyruvate, β-hydroxybutyrate), antioxidant status, respiratory chain activity, or history of oxidative challenge (e.g., ischemia-reperfusion). This underscores the relevance of quantitative methods to characterize the network's functional status as a way to assess the myocyte's resilience to pathological stressors.

The capillary bed offers the largest hemodynamic resistance to the cortical blood supply

The cortical angioarchitecture is a key factor in controlling cerebral blood flow and oxygen metabolism. Difficulties in imaging the complex microanatomy of the cortex have so far restricted insight about blood flow distribution in the microcirculation. A new methodology combining advanced microscopy data with large scale hemodynamic simulations enabled us to quantify the effect of the angioarchitecture on the cerebral microcirculation. High-resolution images of the mouse primary somatosensory cortex were input into with a comprehensive computational model of cerebral perfusion and oxygen supply ranging from the pial vessels to individual brain cells. Simulations of blood flow, hematocrit and oxygen tension show that the wide variation of hemodynamic states in the tortuous, randomly organized capillary bed is responsible for relatively uniform cortical tissue perfusion and oxygenation. Computational analysis of microcirculatory blood flow and pressure drops further indicates that the capillary bed, including capillaries adjacent to feeding arterioles (d < 10 µm), are the largest contributors to hydraulic resistance.

Computational Models of Reactive Oxygen Species as Metabolic Byproducts and Signal-Transduction Modulators

Reactive oxygen species (ROS) are widely involved in intracellular signaling and human pathologies, but their precise roles have been difficult to enumerate and integrate holistically. The context- and dose-dependent intracellular effects of ROS can lead to contradictory experimental results and confounded interpretations. For example, lower levels of ROS promote cell signaling and proliferation, whereas abundant ROS cause overwhelming damage to biomolecules and cellular apoptosis or senescence. These complexities raise the question of whether the many facets of ROS biology can be joined under a common mechanistic framework using computational modeling. Here, we take inventory of some current models for ROS production or ROS regulation of signaling pathways. Several models captured non-intuitive observations or made predictions that were later verified by experiment. There remains a need for systems-level analyses that jointly incorporate ROS production, handling, and modulation of multiple signal-transduction cascades.

Investigating mitochondrial redox state using NADH and NADPH autofluorescence

The redox states of the NAD and NADP pyridine nucleotide pools play critical roles in defining the activity of energy producing pathways, in driving oxidative stress and in maintaining antioxidant defences. Broadly speaking, NAD is primarily engaged in regulating energy-producing catabolic processes, whilst NADP may be involved in both antioxidant defence and free radical generation. Defects in the balance of these pathways are associated with numerous diseases, from diabetes and neurodegenerative disease to heart disease and cancer. As such, a method to assess the abundance and redox state of these separate pools in living tissues would provide invaluable insight into the underlying pathophysiology. Experimentally, the intrinsic fluorescence of the reduced forms of both redox cofactors, NADH and NADPH, has been used for this purpose since the mid-twentieth century. In this review, we outline the modern implementation of these techniques for studying mitochondrial redox state in complex tissue preparations. As the fluorescence spectra of NADH and NADPH are indistinguishable, interpreting the signals resulting from their combined fluorescence, often labelled NAD(P)H, can be complex. We therefore discuss recent studies using fluorescence lifetime imaging microscopy (FLIM) which offer the potential to discriminate between the two separate pools. This technique provides increased metabolic information from cellular autofluorescence in biomedical investigations, offering biochemical insights into the changes in time-resolved NAD(P)H fluorescence signals observed in diseased tissues.