A Reaction-Diffusion Model of ROS-Induced ROS Release in a Mitochondrial Network

Loss of mitochondrial function is a fundamental determinant of cell injury and death. In heart cells under metabolic stress, we have previously described how the abrupt collapse or oscillation of the mitochondrial energy state is synchronized across the mitochondrial network by local interactions dependent upon reactive oxygen species (ROS). Here, we develop a mathematical model of ROS-induced ROS release (RIRR) based on reaction-diffusion (RD-RIRR) in one- and two-dimensional mitochondrial networks. The nodes of the RD-RIRR network are comprised of models of individual mitochondria that include a mechanism of ROS-dependent oscillation based on the interplay between ROS production, transport, and scavenging; and incorporating the tricarboxylic acid (TCA) cycle, oxidative phosphorylation, and Ca2+ handling. Local mitochondrial interaction is mediated by superoxide (O2.−) diffusion and the O2.−-dependent activation of an inner membrane anion channel (IMAC). In a 2D network composed of 500 mitochondria, model simulations reveal ΔΨm depolarization waves similar to those observed when isolated guinea pig cardiomyocytes are subjected to a localized laser-flash or antioxidant depletion. The sensitivity of the propagation rate of the depolarization wave to O2.− diffusion, production, and scavenging in the reaction-diffusion model is similar to that observed experimentally. In addition, we present novel experimental evidence, obtained in permeabilized cardiomyocytes, confirming that ΔΨm depolarization is mediated specifically by O2.−. The present work demonstrates that the observed emergent macroscopic properties of the mitochondrial network can be reproduced in a reaction-diffusion model of RIRR. Moreover, the findings have uncovered a novel aspect of the synchronization mechanism, which is that clusters of mitochondria that are oscillating can entrain mitochondria that would otherwise display stable dynamics. The work identifies the fundamental mechanisms leading from the failure of individual organelles to the whole cell, thus it has important implications for understanding cell death during the progression of heart disease.Author Summary: Cardiac cell injury and death is a key component of cardiac diseases such as heart failure or myocardial infarction, thus it is important to understand the earliest steps leading up to irreversible cell damage. Mitochondria are the organelles responsible for generating the energy required to keep the cell running, yet they are particularly vulnerable to damage by toxic byproducts of metabolism, which include reactive oxygen species (ROS). ROS wreak havoc on cellular functions by attacking proteins, lipids and DNA, so the cardiac cell has evolved sophisticated defenses to remove them. The work we present in this paper using a computer model of the mitochondrial network describes how ROS generated inside the cell can spread from one mitochondrion to the next in a positive feedback process known as ROS-induced ROS release. Understanding this process will help us to find ways to intervene in this catastrophic mechanism to prevent loss of cell function and the associated cardiac arrhythmias and contractile problems leading to sudden cardiac death.