A computational model of reactive oxygen species and redox balance in cardiac mitochondria

Unlocking mitochondrial dysfunction-associated senescence (MiDAS) with NAD+ - A Boolean model of mitochondrial dynamics and cell cycle control

The steady accumulation of senescent cells with aging creates tissue environments that aid cancer evolution. Aging cell states are highly heterogeneous. 'Deep senescent' cells rely on healthy mitochondria to fuel a strong proinflammatory secretome, including cytokines, growth and transforming signals. Yet, the physiological triggers of senescence such as reactive oxygen species (ROS) can also trigger mitochondrial dysfunction, and sufficient energy deficit to alter their secretome and cause chronic oxidative stress - a state termed Mitochondrial Dysfunction-Associated Senescence (MiDAS). Here, we offer a mechanistic hypothesis for the molecular processes leading to MiDAS, along with testable predictions. To do this we have built a Boolean regulatory network model that qualitatively captures key aspects of mitochondrial dynamics during cell cycle progression (hyper-fusion at the G1/S boundary, fission in mitosis), apoptosis (fission and dysfunction) and glucose starvation (reversible hyper-fusion), as well as MiDAS in response to SIRT3 knockdown or oxidative stress. Our model reaffirms the protective role of NAD+ and external pyruvate. We offer testable predictions about the growth factor- and glucose-dependence of MiDAS and its reversibility at different stages of reactive oxygen species (ROS)-induced senescence. Our model provides mechanistic insights into the distinct stages of DNA-damage induced senescence, the relationship between senescence and epithelial-to-mesenchymal transition in cancer and offers a foundation for building multiscale models of tissue aging.

The Interplay between Mechanoregulation and ROS in Heart Physiology, Disease, and Regeneration

Cardiovascular diseases are currently the most common cause of death in developed countries. Due to lifestyle and environmental factors, this problem is only expected to increase in the future. Reactive oxygen species (ROS) are a key player in the onset of cardiovascular diseases but also have important functions in healthy cardiac tissue. Here, the interplay between ROS generation and cardiac mechanical forces is shown, and the state of the art and a perspective on future directions are discussed. To this end, an overview of what is currently known regarding ROS and mechanosignaling at a subcellular level is first given. There the role of ROS in mechanosignaling as well as the interplay between both factors in specific organelles is emphasized. The consequences at a larger scale across the population of heart cells are then discussed. Subsequently, the roles of ROS in embryogenesis, pathogenesis, and aging are further discussed, exemplifying some aspects of mechanoregulation. Finally, different models that are currently in use are discussed to study the topics above.

Unlocking Mitochondrial Dysfunction-Associated Senescence (MiDAS) with NAD + - a Boolean Model of Mitochondrial Dynamics and Cell Cycle Control

The steady accumulation of senescent cells with aging creates tissue environments that aid cancer evolution. Aging cell states are highly heterogeneous. 'Deep senescent' cells rely on healthy mitochondria to fuel a strong proinflammatory secretome, including cytokines, growth and transforming signals. Yet, the physiological triggers of senescence such as the reactive oxygen species (ROS) can also trigger mitochondrial dysfunction, and sufficient energy deficit to alter their secretome and cause chronic oxidative stress - a state termed Mitochondrial Dysfunction-Associated Senescence (MiDAS). Here, we offer a mechanistic hypothesis for the molecular processes leading to MiDAS, along with testable predictions. To do this we have built a Boolean regulatory network model that qualitatively captures key aspects of mitochondrial dynamics during cell cycle progression (hyper-fusion at the G1/S boundary, fission in mitosis), apoptosis (fission and dysfunction) and glucose starvation (reversible hyper-fusion), as well as MiDAS in response to SIRT3 knockdown or oxidative stress. Our model reaffirms the protective role of NAD + and external pyruvate. We offer testable predictions about the growth factor- and glucose-dependence of MiDAS and its reversibility at different stages of reactive oxygen species (ROS)-induced senescence. Our model provides mechanistic insights into the distinct stages of DNA-damage induced senescence, the relationship between senescence and epithelial-to-mesenchymal transition in cancer and offers a foundation for building multiscale models of tissue aging.

Highlights

Boolean regulatory network model reproduces mitochondrial dynamics during cell cycle progression, apoptosis, and glucose starvation. Model offers a mechanistic explanation for the positive feedback loop that locks in Mitochondrial Dysfunction-Associated Senescence (MiDAS), involving autophagy-resistant, hyperfused, dysfunctional mitochondria. Model reproduces ROS-mediated mitochondrial dysfunction and suggests that MiDAS is part of the early phase of damage-induced senescence. Model predicts that cancer-driving mutations that bypass the G1/S checkpoint generally increase the incidence of MiDAS, except for p53 loss.

Oxidative Stress in Sepsis: A Focus on Cardiac Pathology

This study aims to analyze post-mortem human cardiac specimens, to verify and evaluate the existence or extent of oxidative stress in subjects whose cause of death has been traced to sepsis, through immunohistological oxidative/nitrosative stress markers. Indeed, in the present study, i-NOS, NOX2, and nitrotyrosine markers were higher expressed in the septic death group when compared to the control group, associated with also a significant increase in 8-OHdG, highlighting the pivotal role of oxidative stress in septic etiopathogenesis. In particular, 70% of cardiomyocyte nuclei from septic death specimens showed positivity for 8-OHdG. Furthermore, intense and massive NOX2-positive myocyte immunoreaction was noticed in the septic group, as nitrotyrosine immunostaining intense reaction was found in the cardiac cells. These results demonstrated a correlation between oxidative and nitrosative stress imbalance and the pathophysiology of cardiac dysfunction documented in cases of sepsis. Therefore, subsequent studies will focus on the expression of oxidative stress markers in other organs and tissues, as well as on the involvement of the intracellular pattern of apoptosis, to better clarify the complex pathogenesis of multi-organ failure, leading to support the rationale for including therapies targeting redox abnormalities in the management of septic patients.

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.

Calcium and Reactive Oxygen Species Signaling Interplays in Cardiac Physiology and Pathologies

Mitochondria are key players in energy production, critical activity for the smooth functioning of energy-demanding organs such as the muscles, brain, and heart. Therefore, dysregulation or alterations in mitochondrial bioenergetics primarily perturb these organs. Within the cell, mitochondria are the major site of reactive oxygen species (ROS) production through the activity of different enzymes since it is one of the organelles with the major availability of oxygen. ROS can act as signaling molecules in a number of different pathways by modulating calcium (Ca2+) signaling. Interactions among ROS and calcium signaling can be considered bidirectional, with ROS regulating cellular Ca2+ signaling, whereas Ca2+ signaling is essential for ROS production. In particular, we will discuss how alterations in the crosstalk between ROS and Ca2+ can lead to mitochondrial bioenergetics dysfunctions and the consequent damage to tissues at high energy demand, such as the heart. Changes in Ca2+ can induce mitochondrial alterations associated with reduced ATP production and increased production of ROS. These changes in Ca2+ levels and ROS generation completely paralyze cardiac contractility. Thus, ROS can hinder the excitation-contraction coupling, inducing arrhythmias, hypertrophy, apoptosis, or necrosis of cardiac cells. These interplays in the cardiovascular system are the focus of this review.

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.

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.

Computational Models on Pathological Redox Signalling Driven by Pregnancy: A Review

Oxidative stress is associated with a myriad of diseases including pregnancy pathologies with long-term cardiovascular repercussions for both the mother and baby. Aberrant redox signalling coupled with deficient antioxidant defence leads to chronic molecular impairment. Abnormal placentation has been considered the primary source for reactive species; however, placental dysfunction has been deemed secondary to maternal cardiovascular maladaptation in pregnancy. While various therapeutic interventions, aimed at combating deregulated oxidative stress during pregnancy have shown promise in experimental models, they often result as inconclusive or detrimental in clinical trials, warranting the need for further research to identify candidates. The strengths and limitations of current experimental methods in redox research are discussed. Assessment of redox status and oxidative stress in experimental models and in clinical practice remains challenging; the state-of-the-art of computational models in this field is presented in this review, comparing static and dynamic models which provide functional information such as protein-protein interactions, as well as the impact of changes in molecular species on the redox-status of the system, respectively. Enhanced knowledge of redox biology in during pregnancy through computational modelling such as generation of Systems Biology Markup Language model which integrates existing models to a larger network in the context of placenta physiology.

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.

High-Repetition-Rate Transient Absorption Spectroscopy of Respiratory Supercomplexes

Hemeproteins are frequent subjects for ultrafast transient absorption spectroscopy (TAS) because of biological importance, strong UV-vis absorption, high photostability, and interesting transient dynamics that depend on redox, conformation, and ligand binding. TAS on hemeproteins is usually performed on isolated, purified proteins, though their response is likely to be different in their native molecular environment, which involves the formation of protein complexes and supercomplexes. Recently, we reported a transient absorption microscopy (TAM) experiment which elicited a transient response from hemeproteins in intact biological tissue using a visible-wavelength pump (530 nm) and probe (490 nm). Here, we find that adaptive noise canceling plus resonant galvanometer scanning enables a high-repetition-rate fiber laser source to make redox-sensitive measurements of cytochrome c (Cyt-c). We investigate the origins of the visible-wavelength response of biological tissue through TAS of intact mitochondrial respiratory supercomplexes, separated via gel electrophoresis. We find that each of these high-molecular-weight gel bands yields a TAS response characteristic of cytochrome hemes, implying that the TAS response of intact cells and tissue originates from not just Cyt-c but a mixture of respiratory cytochromes. We also find differences in excited-state lifetime between wild-type (WT) and a tafazzin-deficient (TAZ) mouse model of mitochondrial disease.

A computational model of cardiomyocyte metabolism predicts unique reperfusion protocols capable of reducing cell damage during ischemia/reperfusion

If a coronary blood vessel is occluded and the neighboring cardiomyocytes deprived of oxygen, subsequent reperfusion of the ischemic tissue can lead to oxidative damage due to excessive generation of reactive oxygen species. Cardiomyocytes and their mitochondria are the main energy producers and consumers of the heart, and their metabolic changes during ischemia seem to be a key driver of reperfusion injury. Here, we hypothesized that tracking changes in cardiomyocyte metabolism, such as oxygen and ATP concentrations, would help in identifying points of metabolic failure during ischemia and reperfusion. To track some of these changes continuously from the onset of ischemia through reperfusion, we developed a system of differential equations representing the chemical reactions involved in the production and consumption of 67 molecular species. This model was validated and used to identify conditions present during periods of critical transition in ischemia and reperfusion that could lead to oxidative damage. These simulations identified a range of oxygen concentrations that lead to reverse mitochondrial electron transport at complex I of the respiratory chain and a spike in mitochondrial membrane potential, which are key suspects in the generation of reactive oxygen species at the onset of reperfusion. Our model predicts that a short initial reperfusion treatment with reduced oxygen content (5% of physiological levels) could reduce the cellular damage from both of these mechanisms. This model should serve as an open-source platform to test ideas for treatment of the ischemia reperfusion process by following the temporal evolution of molecular concentrations in the cardiomyocyte.

In Sickness and in Health: The Oxygen Reactive Species and the Bone

Oxidative stress plays a central role in physiological and pathological bone conditions. Its role in signalment and control of bone cell population differentiation, activity, and fate is increasingly recognized. The possibilities of its use and manipulation with therapeutic goals are virtually unending. However, how redox balance interplays with the response to mechanical stimuli is yet to be fully understood. The present work summarizes current knowledge on these aspects, in an integrative and broad introductory perspective.

Invited review: Thermal effects on oxidative stress in vertebrate ectotherms

Human-induced climate change is occurring rapidly. Ectothermic organisms are particularly vulnerable to these temperature changes due to their reliance on environmental temperature. The extent of ectothermic thermal adaptation and plasticity in the literature is well documented; however, the role of oxidative stress in these processes needs more attention. Oxidative stress occurs when reactive oxygen species, generated mainly through aerobic respiration, overwhelm antioxidant defences and damage crucial biomolecules. The effects of oxidative damage include the alteration of life-history traits and reductions in whole-organism fitness. Here we review the literature addressing experimental temperature effects on oxidative stress in vertebrate ectotherms. Acute and acclimation temperature treatments produce distinctly different results and highlight the role of phylogeny and thermal adaptation in shaping oxidative stress responses. Acute treatments on organisms adapted to stable environments generally produced significant oxidative stress responses, whilst organisms adapted to variable conditions exhibited capacity to cope with temperature changes and mitigate oxidative stress. In acclimation treatments, the temperature treatments higher than optimal temperatures tended to produce significantly less oxidative stress than lower temperatures in reptiles, whilst in some eurythermal fish species, no oxidative stress response was observed. These results highlight the importance of phylogeny and adaptation to past environmental conditions for temperature-dependent oxidative stress responses. We conclude with recommendations on experimental procedures to investigate these phenomena with reference to thermal plasticity, adaptation and biogeographic variation that provide the most significant benefits to adaptable populations. These results have potential conservation ramifications as they may shed light on the physiological effects of temperature alterations in some vertebrate ectotherms.

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.

A simple model of cardiac mitochondrial respiration with experimental validation

Cardiac mitochondria are intracellular organelles that play an important role in energy metabolism and cellular calcium regulation. In particular, they influence the excitation-contraction cycle of the heart cell. A large number of mathematical models have been proposed to better understand the mitochondrial dynamics, but they generally show a high level of complexity, and their parameters are very hard to fit to experimental data. We derived a model based on historical free energy-transduction principles, and results from the literature. We proposed simple expressions that allow to reduce the number of parameters to a minimum with respect to the mitochondrial behavior of interest for us. The resulting model has thirty-two parameters, which are reduced to twenty-three after a global sensitivity analysis of its expressions based on Sobol indices. We calibrated our model to experimental data that consists of measurements of mitochondrial respiration rates controlled by external ADP additions. A sensitivity analysis of the respiration rates showed that only seven parameters can be identified using these observations. We calibrated them using a genetic algorithm, with five experimental data sets. At last, we used the calibration results to verify the ability of the model to accurately predict the values of a sixth dataset. Results show that our model is able to reproduce both respiration rates of mitochondria and transitions between those states, with very low variability of the parameters between each experiment. The same methodology may apply to recover all the parameters of the model, if corresponding experimental data were available.

Cardioprotective effect of anisodamine against ischemia/reperfusion injury through the mitochondrial ATP-sensitive potassium channel

Previous clinical studies have shown that anisodamine could improve no-reflow phenomenon and prevent reperfusion arrhythmias, but whether this protective effect is related to the antagonism of the M-type cholinergic receptor or other potential mechanisms is uncertain. The aim of the present study was to investigate the role of the mitochondrial ATP-sensitive potassium channel (mitoK ATP ) in cardioprotective effect of anisodamine against ischemia/reperfusion injury. Anisodamine and 5- hydroxydecanoic acid were used to explore the relationship between anisodamine and mitoK ATP . Using a Langendorff isolated heart ischemia/reperfusion injury model, hemodynamic parameters and reperfusion ventricular arrhythmia were evaluated; in addition, changes in myocardial infarct size, cTnI from coronary effluent and myocardial ultrastructure, as well as ATP, MDA and SOD in myocardial tissues, were detected. In the hypoxia/reoxygenation injury model of neonatal rat cardiomyocyte, cTnI release in the culture medium and levels of ATP, MDA and SOD in cardiomyocytes and mitochondrial membrane potential, were analyzed. Overall, anisodamine could significantly improve the hemodynamic indexes of isolated rat heart injured by ischemia/reperfusion, reduce the occurrence of ventricular reperfusion arrhythmia and myocardial infarction area, and improve the ultrastructural damage of myocardium and mitochondria. The in vitro results demonstrated that anisodamine could improve mitochondrial energy metabolism, reduce oxidative stress and stabilize mitochondrial membrane potential. The cardioprotective effects were significantly inhibited by 5-hydroxydecanoic acid. In conclusion, this study suggests that the opening of mitoK ATP could play an important role in the protective effect of anisodamine against myocardial ischemia/reperfusion injury.

Mitochondrial reactive oxygen species in physiology and disease

Mitochondrial reactive oxygen species (mROS) are routinely produced at several sites within the organelle. The balance in their formation and elimination is maintained by a complex and robust antioxidant system. mROS may act as second messengers and regulate a number of physiological processes, such as insulin signaling, cell differentiation and proliferation, wound healing, etc. Nevertheless, when a sudden or sustained increase in ROS formation is not efficiently neutralized by the endogenous antioxidant defense system, the detrimental impact of high mROS levels on cell function and viability eventually results in disease development. In this review, we will focus on the dual role of mROS in pathophysiology, emphasizing the physiological role exerted by a regulated mROS production/elimination, and discussing the detrimental effects evoked by an imbalance in mitochondrial redox state. Furthermore, we will touch upon the interplay between mROS and Ca²⁺ homeostasis.

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.

Network integration and modelling of dynamic drug responses at multi-omics levels

Uncovering cellular responses from heterogeneous genomic data is crucial for molecular medicine in particular for drug safety. This can be realized by integrating the molecular activities in networks of interacting proteins. As proof-of-concept we challenge network modeling with time-resolved proteome, transcriptome and methylome measurements in iPSC-derived human 3D cardiac microtissues to elucidate adverse mechanisms of anthracycline cardiotoxicity measured with four different drugs (doxorubicin, epirubicin, idarubicin and daunorubicin). Dynamic molecular analysis at in vivo drug exposure levels reveal a network of 175 disease-associated proteins and identify common modules of anthracycline cardiotoxicity in vitro, related to mitochondrial and sarcomere function as well as remodeling of extracellular matrix. These in vitro-identified modules are transferable and are evaluated with biopsies of cardiomyopathy patients. This to our knowledge most comprehensive study on anthracycline cardiotoxicity demonstrates a reproducible workflow for molecular medicine and serves as a template for detecting adverse drug responses from complex omics data.

Using a network propagation approach with integrated multi-omic data, Selevsek et al. develop a reproducible workflow for identifying drug toxicity effects in cellular systems. This is demonstrated with the analysis of anthracycline cardiotoxicity in cardiac microtissues under the effect of multiple drugs.

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.

Mitochondrial CaMKII causes adverse metabolic reprogramming and dilated cardiomyopathy

Despite the clear association between myocardial injury, heart failure and depressed myocardial energetics, little is known about upstream signals responsible for remodeling myocardial metabolism after pathological stress. Here, we report increased mitochondrial calmodulin kinase II (CaMKII) activation and left ventricular dilation in mice one week after myocardial infarction (MI) surgery. By contrast, mice with genetic mitochondrial CaMKII inhibition are protected from left ventricular dilation and dysfunction after MI. Mice with myocardial and mitochondrial CaMKII overexpression (mtCaMKII) have severe dilated cardiomyopathy and decreased ATP that causes elevated cytoplasmic resting (diastolic) Ca2+ concentration and reduced mechanical performance. We map a metabolic pathway that rescues disease phenotypes in mtCaMKII mice, providing insights into physiological and pathological metabolic consequences of CaMKII signaling in mitochondria. Our findings suggest myocardial dilation, a disease phenotype lacking specific therapies, can be prevented by targeted replacement of mitochondrial creatine kinase or mitochondrial-targeted CaMKII inhibition.

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.

Cardiac-specific Conditional Knockout of the 18-kDa Mitochondrial Translocator Protein Protects from Pressure Overload Induced Heart Failure

Heart failure (HF) is characterized by abnormal mitochondrial calcium (Ca²⁺) handling, energy failure and impaired mitophagy resulting in contractile dysfunction and myocyte death. We have previously shown that the 18-kDa mitochondrial translocator protein of the outer mitochondrial membrane (TSPO) can modulate mitochondrial Ca²⁺ uptake. Experiments were designed to test the role of the TSPO in a murine pressure-overload model of HF induced by transverse aortic constriction (TAC). Conditional, cardiac-specific TSPO knockout (KO) mice were generated using the Cre-loxP system. TSPO-KO and wild-type (WT) mice underwent TAC for 8 weeks. TAC-induced HF significantly increased TSPO expression in WT mice, associated with a marked reduction in systolic function, mitochondrial Ca²⁺ uptake, complex I activity and energetics. In contrast, TSPO-KO mice undergoing TAC had preserved ejection fraction, and exhibited fewer clinical signs of HF and fibrosis. Mitochondrial Ca²⁺ uptake and energetics were restored in TSPO KO mice, associated with decreased ROS, improved complex I activity and preserved mitophagy. Thus, HF increases TSPO expression, while preventing this increase limits the progression of HF, preserves ATP production and decreases oxidative stress, thereby preventing metabolic failure. These findings suggest that pharmacological interventions directed at TSPO may provide novel therapeutics to prevent or treat HF.

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.

Analysis of a Functional Dimer Model of Ubiquinol Cytochrome c Oxidoreductase

Ubiquinol cytochrome c oxidoreductase (bc₁ complex) serves as an important electron junction in many respiratory systems. It funnels electrons coming from NADH and ubiquinol to cytochrome c, but it is also capable of producing significant amounts of the free radical superoxide. In situ and in other experimental systems, the enzyme exists as a dimer. But until recently, it was believed to operate as a functional monomer. Here we show that a functional dimer model is capable of explaining both kinetic and superoxide production rate data. The model consists of six electronic states characterized by the number of electrons deposited on the complex. It is fully reversible and strictly adheres to the thermodynamics governing the reactions. A total of nine independent data sets were used to parameterize the model. To explain the data with a consistent set of parameters, it was necessary to incorporate intramonomer Coulombic effects between hemes b_L and b_H and intermonomer Coulombic effects between b_L hemes. The fitted repulsion energies fall within the theoretical range of electrostatic calculations. In addition, model analysis demonstrates that the Q pool is mostly oxidized under normal physiological operation but can switch to a more reduced state when reverse electron transport conditions are in place.

Computational Modeling in Liver Surgery

The need for extended liver resection is increasing due to the growing incidence of liver tumors in aging societies. Individualized surgical planning is the key for identifying the optimal resection strategy and to minimize the risk of postoperative liver failure and tumor recurrence. Current computational tools provide virtual planning of liver resection by taking into account the spatial relationship between the tumor and the hepatic vascular trees, as well as the size of the future liver remnant. However, size and function of the liver are not necessarily equivalent. Hence, determining the future liver volume might misestimate the future liver function, especially in cases of hepatic comorbidities such as hepatic steatosis. A systems medicine approach could be applied, including biological, medical, and surgical aspects, by integrating all available anatomical and functional information of the individual patient. Such an approach holds promise for better prediction of postoperative liver function and hence improved risk assessment. This review provides an overview of mathematical models related to the liver and its function and explores their potential relevance for computational liver surgery. We first summarize key facts of hepatic anatomy, physiology, and pathology relevant for hepatic surgery, followed by a description of the computational tools currently used in liver surgical planning. Then we present selected state-of-the-art computational liver models potentially useful to support liver surgery. Finally, we discuss the main challenges that will need to be addressed when developing advanced computational planning tools in the context of liver surgery.

Endogenous nitric oxide formation in cardiac myocytes does not control respiration during β-adrenergic stimulation

KEY POINTS

In the heart, endothelial nitric oxide (NO) controls oxygen consumption in the working heart through paracrine mechanisms. While cardiac myocytes contain several isoforms of NO synthases, it is unclear whether these can control respiration in an intracrine fashion. A long-standing controversy is whether a NOS exists within mitochondria. By combining fluorescence technologies with electrical field stimulation or the patch-clamp technique in beating cardiac myocytes, we identified a neuronal NO synthase (nNOS) as the most relevant source of intracellular NO during β-adrenergic stimulation, while no evidence for a mitochondria-located NOS was obtained. The amounts of NO produced by non-mitochondrial nNOS were insufficient to regulate respiration during β-adrenergic stimulation, arguing against intracrine control of respiration by NO within cardiac myocytes.

ABSTRACT

Endothelial nitric oxide (NO) controls cardiac oxygen (O₂ ) consumption in a paracrine way by slowing respiration at the mitochondrial electron transport chain. While NO synthases (NOSs) are also expressed in cardiac myocytes, it is unclear whether they control respiration in an intracrine way. Furthermore, the existence of a mitochondrial NOS is controversial. Here, by combining fluorescence imaging with electrical field stimulation, the patch-clamp method and knock-out technology, we determined the sources and consequences of intracellular NO formation during workload transitions in isolated murine and guinea pig cardiac myocytes and mitochondria. Using 4-amino-5-methylamino-2',7'-difluorofluorescein diacetate (DAF) as a fluorescent NO-sensor that locates to the cytosol and mitochondria, we observed that NO increased by ∼12% within 3 min of β-adrenergic stimulation in beating cardiac myocytes. This NO stems from neuronal NOS (nNOS), but not endothelial (eNOS). After patch clamp-mediated dialysis of cytosolic DAF, the remaining NO signals (mostly mitochondrial) were blocked by nNOS deletion, but not by inhibiting the mitochondrial Ca²⁺ uniporter with Ru360. While in isolated mitochondria exogenous NO inhibited respiration and reduced the NAD(P)H redox state, pyridine nucleotide redox states were unaffected by pharmacological or genetic disruption of endogenous nNOS or eNOS during workload transitions in cardiac myoctyes. We conclude that under physiological conditions, nNOS is the most relevant source for NO in cardiac myocytes, but this nNOS is not located in mitochondria and does not control respiration. Therefore, cardiac O₂ consumption is controlled by endothelial NO in a paracrine, but not intracrine, fashion.

Modeling oxygen requirements in ischemic cardiomyocytes

Heart disease remains the leading cause of death globally. Although reperfusion following myocardial ischemia can prevent death by restoring nutrient flow, ischemia/reperfusion injury can cause significant heart damage. The mechanisms that drive ischemia/reperfusion injury are not well understood; currently, few methods can predict the state of the cardiac muscle cell and its metabolic conditions during ischemia. Here, we explored the energetic sustainability of cardiomyocytes, using a model for cellular metabolism to predict the levels of ATP following hypoxia. We modeled glycolytic metabolism with a system of coupled ordinary differential equations describing the individual metabolic reactions within the cardiomyocyte over time. Reduced oxygen levels and ATP consumption rates were simulated to characterize metabolite responses to ischemia. By tracking biochemical species within the cell, our model enables prediction of the cell's condition up to the moment of reperfusion. The simulations revealed a distinct transition between energetically sustainable and unsustainable ATP concentrations for various energetic demands. Our model illustrates how even low oxygen concentrations allow the cell to perform essential functions. We found that the oxygen level required for a sustainable level of ATP increases roughly linearly with the ATP consumption rate. An extracellular O₂ concentration of ∼0.007 mm could supply basic energy needs in non-beating cardiomyocytes, suggesting that increased collateral circulation may provide an important source of oxygen to sustain the cardiomyocyte during extended ischemia. Our model provides a time-dependent framework for studying various intervention strategies to change the outcome of reperfusion.

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.

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.

A Biophysical Systems Approach to Identifying the Pathways of Acute and Chronic Doxorubicin Mitochondrial Cardiotoxicity

The clinical use of the anthracycline doxorubicin is limited by its cardiotoxicity which is associated with mitochondrial dysfunction. Redox cycling, mitochondrial DNA damage and electron transport chain inhibition have been identified as potential mechanisms of toxicity. However, the relative roles of each of these proposed mechanisms are still not fully understood. The purpose of this study is to identify which of these pathways independently or in combination are responsible for doxorubicin toxicity. A state of the art mathematical model of the mitochondria including the citric acid cycle, electron transport chain and ROS production and scavenging systems was extended by incorporating a novel representation for mitochondrial DNA damage and repair. In silico experiments were performed to quantify the contributions of each of the toxicity mechanisms to mitochondrial dysfunction during the acute and chronic stages of toxicity. Simulations predict that redox cycling has a minor role in doxorubicin cardiotoxicity. Electron transport chain inhibition is the main pathway for acute toxicity for supratherapeutic doses, being lethal at mitochondrial concentrations higher than 200μM. Direct mitochondrial DNA damage is the principal pathway of chronic cardiotoxicity for therapeutic doses, leading to a progressive and irreversible long term mitochondrial dysfunction.

Evaluating biomarkers to model cancer risk post cosmic ray exposure

Robust predictive models are essential to manage the risk of radiation-induced carcinogenesis. Chronic exposure to cosmic rays in the context of the complex deep space environment may place astronauts at high cancer risk. To estimate this risk, it is critical to understand how radiation-induced cellular stress impacts cell fate decisions and how this in turn alters the risk of carcinogenesis. Exposure to the heavy ion component of cosmic rays triggers a multitude of cellular changes, depending on the rate of exposure, the type of damage incurred and individual susceptibility. Heterogeneity in dose, dose rate, radiation quality, energy and particle flux contribute to the complexity of risk assessment. To unravel the impact of each of these factors, it is critical to identify sensitive biomarkers that can serve as inputs for robust modeling of individual risk of cancer or other long-term health consequences of exposure. Limitations in sensitivity of biomarkers to dose and dose rate, and the complexity of longitudinal monitoring, are some of the factors that increase uncertainties in the output from risk prediction models. Here, we critically evaluate candidate early and late biomarkers of radiation exposure and discuss their usefulness in predicting cell fate decisions. Some of the biomarkers we have reviewed include complex clustered DNA damage, persistent DNA repair foci, reactive oxygen species, chromosome aberrations and inflammation. Other biomarkers discussed, often assayed for at longer points post exposure, include mutations, chromosome aberrations, reactive oxygen species and telomere length changes. We discuss the relationship of biomarkers to different potential cell fates, including proliferation, apoptosis, senescence, and loss of stemness, which can propagate genomic instability and alter tissue composition and the underlying mRNA signatures that contribute to cell fate decisions. Our goal is to highlight factors that are important in choosing biomarkers and to evaluate the potential for biomarkers to inform models of post exposure cancer risk. Because cellular stress response pathways to space radiation and environmental carcinogens share common nodes, biomarker-driven risk models may be broadly applicable for estimating risks for other carcinogens.

An Analysis of the Truncated Bid- and ROS-dependent Spatial Propagation of Mitochondrial Permeabilization Waves during Apoptosis

Apoptosis is a form of programmed cell death that is essential for the efficient elimination of surplus, damaged, and transformed cells during metazoan embryonic development and adult tissue homeostasis. Situated at the interface of apoptosis initiation and execution, mitochondrial outer membrane permeabilization (MOMP) represents one of the most fundamental processes during apoptosis signal transduction. It was shown that MOMP can spatiotemporally propagate through cells, in particular in response to extrinsic apoptosis induction. Based on apparently contradictory experimental evidence, two distinct molecular mechanisms have been proposed to underlie the propagation of MOMP signals, namely a reaction-diffusion mechanism governed by anisotropies in the production of the MOMP-inducer truncated Bid (tBid), or a process that drives the spatial propagation of MOMP by sequential bursts of reactive oxygen species. We therefore generated mathematical models for both scenarios and performed in silico simulations of spatiotemporal MOMP signaling to identify which one of the two mechanisms is capable of qualitatively and quantitatively reproducing the existing data. We found that the explanatory power of each model was limited in that only a subset of experimental findings could be replicated. However, the integration of both models into a combined mathematical description of spatiotemporal tBid and reactive oxygen species signaling accurately reproduced all available experimental data and furthermore, provided robustness to spatial MOMP propagation when mitochondria are spatially separated. Our study therefore provides a theoretical framework that is sufficient to describe and mechanistically explain the spatiotemporal propagation of one of the most fundamental processes during apoptotic cell death.

Modeling calcium regulation of contraction, energetics, signaling, and transcription in the cardiac myocyte

Calcium (Ca(2+)) plays many important regulatory roles in cardiac muscle cells. In the initial phase of the action potential, influx of Ca(2+) through sarcolemmal voltage-gated L-type Ca(2+) channels (LCCs) acts as a feed-forward signal that triggers a large release of Ca(2+) from the junctional sarcoplasmic reticulum (SR). This Ca(2+) drives heart muscle contraction and pumping of blood in a process known as excitation-contraction coupling (ECC). Triggered and released Ca(2+) also feed back to inactivate LCCs, attenuating the triggered Ca(2+) signal once release has been achieved. The process of ECC consumes large amounts of ATP. It is now clear that in a process known as excitation-energetics coupling, Ca(2+) signals exert beat-to-beat regulation of mitochondrial ATP production that closely couples energy production with demand. This occurs through transport of Ca(2+) into mitochondria, where it regulates enzymes of the tricarboxylic acid cycle. In excitation-signaling coupling, Ca(2+) activates a number of signaling pathways in a feed-forward manner. Through effects on their target proteins, these interconnected pathways regulate Ca(2+) signals in complex ways to control electrical excitability and contractility of heart muscle. In a process known as excitation-transcription coupling, Ca(2+) acting primarily through signal transduction pathways also regulates the process of gene transcription. Because of these diverse and complex roles, experimentally based mechanistic computational models are proving to be very useful for understanding Ca(2+) signaling in the cardiac myocyte.

Redox mechanisms of cardiomyocyte mitochondrial protection

Oxidative and nitrosative stress are primary contributors to the loss of myocardial tissue in insults ranging from ischemia/reperfusion injury from coronary artery disease and heart transplantation to sepsis-induced myocardial dysfunction and drug-induced myocardial damage. This cell damage caused by oxidative and nitrosative stress leads to mitochondrial protein, DNA, and lipid modifications, which inhibits energy production and contractile function, potentially leading to cell necrosis and/or apoptosis. However, cardiomyocytes have evolved an elegant set of redox-sensitive mechanisms that respond to and contain oxidative and nitrosative damage. These responses include the rapid induction of antioxidant enzymes, mitochondrial DNA repair mechanisms, selective mitochondrial autophagy (mitophagy), and mitochondrial biogenesis. Coordinated cytoplasmic to nuclear cell-signaling and mitochondrial transcriptional responses to the presence of elevated cytoplasmic oxidant production, e.g., H₂O₂, allows nuclear translocation of the Nfe2l2 transcription factor and up-regulation of downstream cytoprotective genes such as heme oxygenase-1 which generates physiologic signals, such as CO that up-regulates Nfe212 gene transcription. Simultaneously, a number of other DNA binding transcription factors are expressed and/or activated under redox control, such as Nuclear Respiratory Factor-1 (NRF-1), and lead to the induction of genes involved in both intracellular and mitochondria-specific repair mechanisms. The same insults, particularly those related to vascular stress and inflammation also produce elevated levels of nitric oxide, which also has mitochondrial protein thiol-protective functions and induces mitochondrial biogenesis through cyclic GMP-dependent and perhaps other pathways. This brief review provides an overview of these pathways and interconnected cardiac repair mechanisms.

Computational modeling analysis of mitochondrial superoxide production under varying substrate conditions and upon inhibition of different segments of the electron transport chain

A computational mechanistic model of superoxide (O2•-) formation in the mitochondrial electron transport chain (ETC) was developed to facilitate the quantitative analysis of factors controlling mitochondrial O2•- production and assist in the interpretation of experimental studies. The model takes into account all individual electron transfer reactions in Complexes I and III. The model accounts for multiple, often seemingly contradictory observations on the effects of ΔΨ and ΔpH, and for the effects of multiple substrate and inhibitor conditions, including differential effects of Complex III inhibitors antimycin A, myxothiazol and stigmatellin. Simulation results confirm that, in addition to O2•- formation in Complex III and at the flavin site of Complex I, the quinone binding site of Complex I is an additional superoxide generating site that accounts for experimental observations on O2•- production during reverse electron transfer. However, our simulation results predict that, when cytochrome c oxidase is inhibited during oxidation of succinate, ROS production at this site is eliminated and almost all superoxide in Complex I is generated by reduced FMN, even when the redox pressure for reverse electron transfer from succinate is strong. In addition, the model indicates that conflicting literature data on the kinetics of electron transfer in Complex III involving the iron-sulfur protein-cytochrome bL complex can be resolved in favor of a dissociation of the protein only after electron transfer to cytochrome bH. The model predictions can be helpful in understanding factors driving mitochondrial superoxide formation in intact cells and tissues.

Quantitative analysis of NRF2 pathway reveals key elements of the regulatory circuits underlying antioxidant response and proliferation of ovarian cancer cells

Cells are constantly exposed to Reactive Oxygen Species (ROS) produced both endogenously to meet physiological requirements and from exogenous sources. While endogenous ROS are considered as important signalling molecules, high uncontrollable ROS are detrimental. It is unclear how cells can achieve a balance between maintaining physiological redox homeostasis and robustly activate the antioxidant system to remove exogenous ROS. We have utilised a Systems Biology approach to understand how this robust adaptive system fulfils homeostatic requirements of maintaining steady-state ROS and growth rate, while undergoing rapid readjustment under challenged conditions. Using a panel of human ovarian and normal cell lines, we experimentally quantified and established interrelationships between key elements of ROS homeostasis. The basal levels of NRF2 and KEAP1 were cell line specific and maintained in tight correlation with their growth rates and ROS. Furthermore, perturbation of this balance triggered cell specific kinetics of NRF2 nuclear-cytoplasmic relocalisation and sequestration of exogenous ROS. Our experimental data were employed to parameterise a mathematical model of the NRF2 pathway that elucidated key response mechanisms of redox regulation and showed that the dynamics of NRF2-H2O2 regulation defines a relationship between half-life, total and nuclear NRF2 level and endogenous H2O2 that is cell line specific.

The "Goldilocks Zone" from a redox perspective-Adaptive vs. deleterious responses to oxidative stress in striated muscle

Consequences of oxidative stress may be beneficial or detrimental in physiological systems. An organ system's position on the “hormetic curve” is governed by the source and temporality of reactive oxygen species (ROS) production, proximity of ROS to moieties most susceptible to damage, and the capacity of the endogenous cellular ROS scavenging mechanisms. Most importantly, the resilience of the tissue (the capacity to recover from damage) is a decisive factor, and this is reflected in the disparate response to ROS in cardiac and skeletal muscle. In myocytes, a high oxidative capacity invariably results in a significant ROS burden which in homeostasis, is rapidly neutralized by the robust antioxidant network. The up-regulation of key pathways in the antioxidant network is a central component of the hormetic response to ROS. Despite such adaptations, persistent oxidative stress over an extended time-frame (e.g., months to years) inevitably leads to cumulative damages, maladaptation and ultimately the pathogenesis of chronic diseases. Indeed, persistent oxidative stress in heart and skeletal muscle has been repeatedly demonstrated to have causal roles in the etiology of heart disease and insulin resistance, respectively. Deciphering the mechanisms that underlie the divergence between adaptive and maladaptive responses to oxidative stress remains an active area of research for basic scientists and clinicians alike, as this would undoubtedly lead to novel therapeutic approaches. Here, we provide an overview of major types of ROS in striated muscle and the divergent adaptations that occur in response to them. Emphasis is placed on highlighting newly uncovered areas of research on this topic, with particular focus on the mitochondria, and the diverging roles that ROS play in muscle health (e.g., exercise or preconditioning) and disease (e.g., cardiomyopathy, ischemia, metabolic syndrome).

Inhibiting Na+/K+ ATPase can impair mitochondrial energetics and induce abnormal Ca2+ cycling and automaticity in guinea pig cardiomyocytes

Cardiac glycosides have been used for the treatment of heart failure because of their capabilities of inhibiting Na⁺/K⁺ ATPase (NKA), which raises [Na⁺]_i and attenuates Ca²⁺ extrusion via the Na⁺/Ca²⁺ exchanger (NCX), causing [Ca²⁺]_i elevation. The resulting [Ca²⁺]_i accumulation further enhances Ca²⁺-induced Ca²⁺ release, generating the positive inotropic effect. However, cardiac glycosides have some toxic and side effects such as arrhythmogenesis, confining their extensive clinical applications. The mechanisms underlying the proarrhythmic effect of glycosides are not fully understood. Here we investigated the mechanisms by which glycosides could cause cardiac arrhythmias via impairing mitochondrial energetics using an integrative computational cardiomyocyte model. In the simulations, the effect of glycosides was mimicked by blocking NKA activity. Results showed that inhibiting NKA not only impaired mitochondrial Ca²⁺ retention (thus suppressed reactive oxygen species (ROS) scavenging) but also enhanced oxidative phosphorylation (thus increased ROS production) during the transition of increasing workload, causing oxidative stress. Moreover, concurrent blocking of mitochondrial Na⁺/Ca²⁺ exchanger, but not enhancing of Ca²⁺ uniporter, alleviated the adverse effects of NKA inhibition. Intriguingly, NKA inhibition elicited Ca²⁺ transient and action potential alternans under more stressed conditions such as severe ATP depletion, augmenting its proarrhythmic effect. This computational study provides new insights into the mechanisms underlying cardiac glycoside-induced arrhythmogenesis. The findings suggest that targeting both ion handling and mitochondria could be a very promising strategy to develop new glycoside-based therapies in the treatment of heart failure.

Modeling CaMKII-mediated regulation of L-type Ca(2+) channels and ryanodine receptors in the heart

Excitation-contraction coupling (ECC) in the cardiac myocyte is mediated by a number of highly integrated mechanisms of intracellular Ca²⁺ transport. Voltage- and Ca²⁺-dependent L-type Ca²⁺ channels (LCCs) allow for Ca²⁺ entry into the myocyte, which then binds to nearby ryanodine receptors (RyRs) and triggers Ca²⁺ release from the sarcoplasmic reticulum in a process known as Ca²⁺-induced Ca²⁺ release. The highly coordinated Ca²⁺-mediated interaction between LCCs and RyRs is further regulated by the cardiac isoform of the Ca²⁺/calmodulin-dependent protein kinase (CaMKII). Because CaMKII targets and modulates the function of many ECC proteins, elucidation of its role in ECC and integrative cellular function is challenging and much insight has been gained through the use of detailed computational models. Multiscale models that can both reconstruct the detailed nature of local signaling events within the cardiac dyad and predict their functional consequences at the level of the whole cell have played an important role in advancing our understanding of CaMKII function in ECC. Here, we review experimentally based models of CaMKII function with a focus on LCC and RyR regulation, and the mechanistic insights that have been gained through their application.

Mitochondrial reactive oxygen species production and elimination

Reactive oxygen species (ROS) play an important role in cardiovascular diseases, and one important source for ROS are mitochondria. Emission of ROS from mitochondria is the net result of ROS production at the electron transport chain (ETC) and their elimination by antioxidative enzymes. Both of these processes are highly dependent on the mitochondrial redox state, which is dynamically altered under different physiological and pathological conditions. The concept of "redox-optimized ROS balance" integrates these aspects and implies that oxidative stress occurs when the optimal equilibrium of an intermediate redox state is disturbed towards either strong oxidation or reduction. Furthermore, mitochondria integrate ROS signals from other cellular sources, presumably through a process termed "ROS-induced ROS release" that involves mitochondrial ion channels. Here, we attempt to integrate these recent advances in our understanding of the control of mitochondrial ROS emission and develop a concept of how in heart failure, defects in ion handling can lead to mitochondrial oxidative stress. This article is part of a Special Issue entitled "Redox Signalling in the Cardiovascular System".

The interplay between autophagy and mitochondrial dysfunction in oxidative stress-induced cardiac aging and pathology

Aging is accompanied by a progressive increase in the incidence and prevalence of cardiovascular disease (CVD). Prolonged exposure to cardiovascular risk factors, together with intrinsic age-dependent declines in cardiac functionality, increases the vulnerability of the heart to both endogenous and exogenous stressors, ultimately enhancing the susceptibility to developing CVD in late life. Both increased levels of oxidative damage and the accumulation of dysfunctional mitochondria have been observed in a wide range of cardiac diseases, which may therefore represent a common ground upon which many aspects of CVD develop. In this review, we summarize the current knowledge on the mechanisms whereby oxidative stress arising from mitochondrial dysfunction is involved in the process of cardiac aging and in the pathogenesis of CVD highly prevalent in late life (e.g., heart failure and ischemic heart disease). Special emphasis is placed on recent evidence about the role played by alterations in cellular quality control systems, in particular autophagy/mitophagy and mitochondrial dynamics (fusion and fission), and their interconnections in the context of age-related CVD. Cardioprotective interventions acting through the modulation of mitochondrial autophagy (calorie restriction, calorie restriction mimetics, and the gasotransmitter hydrogen sulfide) are also presented. This article is part of a Special Issue entitled "Protein Quality Control, the Ubiquitin Proteasome System, and Autophagy".

An integrated mitochondrial ROS production and scavenging model: implications for heart failure

It has been observed experimentally that cells from failing hearts exhibit elevated levels of reactive oxygen species (ROS) upon increases in energetic workload. One proposed mechanism for this behavior is mitochondrial Ca(2+) mismanagement that leads to depletion of ROS scavengers. Here, we present a computational model to test this hypothesis. Previously published models of ROS production and scavenging were combined and reparameterized to describe ROS regulation in the cellular environment. Extramitochondrial Ca(2+) pulses were applied to simulate frequency-dependent changes in cytosolic Ca(2+). Model results show that decreased mitochondrial Ca(2+)uptake due to mitochondrial Ca(2+) uniporter inhibition (simulating Ru360) or elevated cytosolic Na(+), as in heart failure, leads to a decreased supply of NADH and NADPH upon increasing cellular workload. Oxidation of NADPH leads to oxidation of glutathione (GSH) and increased mitochondrial ROS levels, validating the Ca(2+) mismanagement hypothesis. The model goes on to predict that the ratio of steady-state [H2O2]m during 3Hz pacing to [H2O2]m at rest is highly sensitive to the size of the GSH pool. The largest relative increase in [H2O2]m in response to pacing is shown to occur when the total GSH and GSSG is close to 1 mM, whereas pool sizes below 0.9 mM result in high resting H2O2 levels, a quantitative prediction only possible with a computational model.