Complex oscillatory redox dynamics with signaling potential at the edge between normal and pathological mitochondrial function

The time-keeping properties bestowed by oscillatory behavior on functional rhythms represent an evolutionarily conserved trait in living systems. Mitochondrial networks function as timekeepers maximizing energetic output while tuning reactive oxygen species (ROS) within physiological levels compatible with signaling. In this work, we explore the potential for timekeeping functions dependent on mitochondrial dynamics with the validated two-compartment mitochondrial energetic-redox (ME-R) computational model, that takes into account (a) four main redox couples [NADH, NADPH, GSH, Trx(SH)2], (b) scavenging systems (glutathione, thioredoxin, SOD, catalase) distributed in matrix and extra-matrix compartments, and (c) transport of ROS species between them. Herein, we describe that the ME-R model can exhibit highly complex oscillatory dynamics in energetic/redox variables and ROS species, consisting of at least five frequencies with modulated amplitudes and period according to power spectral analysis. By stability analysis we describe that the extent of steady state—as against complex oscillatory behavior—was dependent upon the abundance of Mn and Cu, Zn SODs, and their interplay with ROS production in the respiratory chain. Large parametric regions corresponding to oscillatory dynamics of increasingly complex waveforms were obtained at low Cu, Zn SOD concentration as a function of Mn SOD. This oscillatory domain was greatly reduced at higher levels of Cu, Zn SOD. Interestingly, the realm of complex oscillations was located at the edge between normal and pathological mitochondrial energetic behavior, and was characterized by oxidative stress. We conclude that complex oscillatory dynamics could represent a frequency- and amplitude-modulated H2O2 signaling mechanism that arises under intense oxidative stress. By modulating SOD, cells could have evolved an adaptive compromise between relative constancy and the flexibility required under stressful redox/energetic conditions.Fil: Kembro, Jackelyn Melissa. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Córdoba. Instituto de Investigaciones Biológicas y Tecnológicas. Universidad Nacional de Córdoba. Facultad de Ciencias Exactas, Físicas y Naturales. Instituto de Investigaciones Biológicas y Tecnológicas; ArgentinaFil: Cortassa, Sonia del Carmen. University Johns Hopkins; Estados Unidos. Consejo Nacional de Investigaciones Científicas y Técnicas; ArgentinaFil: Aon, Miguel A.. University Johns Hopkins; Estados Unido

Why Homeodynamics, Not Homeostasis?

Ideas of homeostasis derive from the concept of the organism as an open system. These ideas can be traced back to Heraclitus. Hopkins, Bernard, Hill, Cannon, Weiner and von Bertalanffy developed further the mechanistic basis of turnover of biological components, and Schoenheimer and Rittenberg were pioneers of experimental approaches to the problems of measuring pool sizes and dynamic fluxes. From the second half of the twentieth century, a biophysical theory mainly founded on self-organisation and Dynamic Systems Theory allowed us to approach the quantitative and qualitative analysis of the organised complexity that characterises living systems. This combination of theoretical framework and more refined experimental techniques revealed that feedback control of steady states is a mode of operation that, although providing stability, is only one of many modes and may be the exception rather than the rule. The concept of homeodynamics that we introduce here offers a radically new and all-embracing concept that departs from the classical homeostatic idea that emphasises the stability of the internal milieu toward perturbation. Indeed, biological systems are homeody- namic because of their ability to dynamically self-organise at bifurcation points of their behaviour where they lose stability. Consequently, they exhibit diverse behaviour; in addition to monotonic stationary states, living systems display complex behaviour with all its emergent characteristics, i.e., bistable switches, thresholds, waves, gradients, mutual entrainment, and periodic as well as chaotic behaviour, as evidenced in cellular phenomena such as dynamic (supra)molecular organisation and flux coordination. These processes may proceed on different spatial scales, as well as across time scales, from the very rapid processes within and between molecules in membranes to the slow time scales of evolutionary change. It is dynamic organisation under homeodynamic conditions that make possible the organised complexity of life

Control and Regulation of Integrated Mitochondrial Function in Metabolic and Transport Networks

The pattern of flux and concentration control coefficients in an integrated mitochondrial energetics model is examined by applying a generalized matrix method of control analysis to calculate control coefficients, as well as response coefficients The computational model of Cortassa et al. encompasses oxidative phosphorylation, the TCA cycle, and Ca2+ dynamics. Control of ATP synthesis, TCA cycle, and ANT fluxes were found to be distributed among various mitochondrial processes. Control is shared by processes associated with ATP/ADP production and transport, as well as by Ca2+ dynamics. The calculation also analyzed the control of the concentrations of key regulatory ions and metabolites (Ca2+, NADH, ADP). The approach we have used demonstrates how properties of integrated systems may be understood through applications of computational modeling and control analysis

Aldose Reductase Inhibition or Activation of Transketolase Offset Adverse Metabolic Remodeling Improving Function in Type 2 Diabetes Myocytes Exposed to Hyperglycemia

n/

A Computational Model of Reactive Oxygen Species and Redox Balance in Cardiac Mitochondria

AbstractElevated levels of reactive oxygen species (ROS) play a critical role in cardiac myocyte signaling in both healthy and diseased cells. Mitochondria represent the predominant cellular source of ROS, specifically the activity of complexes I and III. The model presented here explores the modulation of electron transport chain ROS production for state 3 and state 4 respiration and the role of substrates and respiratory inhibitors. Model simulations show that ROS production from complex III increases exponentially with membrane potential (ΔΨm) when in state 4. Complex I ROS release in the model can occur in the presence of NADH and succinate (reverse electron flow), leading to a highly reduced ubiquinone pool, displaying the highest ROS production flux in state 4. In the presence of ample ROS scavenging, total ROS production is moderate in state 3 and increases substantially under state 4 conditions. The ROS production model was extended by combining it with a minimal model of ROS scavenging. When the mitochondrial redox status was oxidized by increasing the proton permeability of the inner mitochondrial membrane, simulations with the combined model show that ROS levels initially decline as production drops off with decreasing ΔΨm and then increase as scavenging capacity is exhausted. Hence, this mechanistic model of ROS production demonstrates how ROS levels are controlled by mitochondrial redox balance

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

The Scale-Free Dynamics of Eukaryotic Cells

Temporal organization of biological processes requires massively parallel processing on a synchronized time-base. We analyzed time-series data obtained from the bioenergetic oscillatory outputs of Saccharomyces cerevisiae and isolated cardiomyocytes utilizing Relative Dispersional (RDA) and Power Spectral (PSA) analyses. These analyses revealed broad frequency distributions and evidence for long-term memory in the observed dynamics. Moreover RDA and PSA showed that the bioenergetic dynamics in both systems show fractal scaling over at least 3 orders of magnitude, and that this scaling obeys an inverse power law. Therefore we conclude that in S. cerevisiae and cardiomyocytes the dynamics are scale-free in vivo. Applying RDA and PSA to data generated from an in silico model of mitochondrial function indicated that in yeast and cardiomyocytes the underlying mechanisms regulating the scale-free behavior are similar. We validated this finding in vivo using single cells, and attenuating the activity of the mitochondrial inner membrane anion channel with 4-chlorodiazepam to show that the oscillation of NAD(P)H and reactive oxygen species (ROS) can be abated in these two evolutionarily distant species. Taken together these data strongly support our hypothesis that the generation of ROS, coupled to redox cycling, driven by cytoplasmic and mitochondrial processes, are at the core of the observed rhythmicity and scale-free dynamics. We argue that the operation of scale-free bioenergetic dynamics plays a fundamental role to integrate cellular function, while providing a framework for robust, yet flexible, responses to the environment

Temporal metabolic partitioning of the yeast and protist cellular networks:the cell is a global scale-invariant (fractal or self-similar) multioscillator

Britton Chance, electronics expert when a teenager, became an enthusiastic student of biological oscillations, passing on this enthusiasm to many students and colleagues, including one of us (DL). This historical essay traces BC’s influence through the accumulated work of DL to DL’s many collaborators. The overall temporal organization of mass-energy, information, and signaling networks in yeast in self-synchronized continuous cultures represents, until now, the most characterized example of in vivo elucidation of time structure. Continuous online monitoring of dissolved gases by direct measurement (membrane-inlet mass spectrometry, together with NAD(P)H and flavin fluorescence) gives strain-specific dynamic information from timescales of minutes to hours as does two-photon imaging. The predominantly oscillatory behavior of network components becomes evident, with spontaneously synchronized cellular respiration cycles between discrete periods of increased oxygen consumption (oxidative phase) and decreased oxygen consumption (reductive phase). This temperature-compensated ultradian clock provides coordination, linking temporally partitioned functions by direct feedback loops between the energetic and redox state of the cell and its growing ultrastructure. Multioscillatory outputs in dissolved gases with 13 h, 40 min, and 4 min periods gave statistical self-similarity in power spectral and relative dispersional analyses: i.e., complex nonlinear (chaotic) behavior and a functional scale-free (fractal) network operating simultaneously over several timescales

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

Age-dependent impact of two exercise training regimens on genomic and metabolic remodeling in skeletal muscle and liver of male mice

Skeletal muscle adapts to different exercise training modalities with age; however, the impact of both variables at the systemic and tissue levels is not fully understood. Here, adult and old C57BL/6 male mice were assigned to one of three groups: sedentary, daily high-intensity intermittent training (HIIT), or moderate intensity continuous training (MICT) for 4 weeks, compatible with the older group’s exercise capacity. Improvements in body composition, fasting blood glucose, and muscle strength were mostly observed in the MICT old group, while effects of HIIT training in adult and old animals was less clear. Skeletal muscle exhibited structural and functional adaptations to exercise training, as revealed by electron microscopy, OXPHOS assays, respirometry, and muscle protein biomarkers. Transcriptomics analysis of gastrocnemius muscle combined with liver and serum metabolomics unveiled an age-dependent metabolic remodeling in response to exercise training. These results support a tailored exercise prescription approach aimed at improving health and ameliorating age-associated loss of muscle strength and function in the elderly.This work was supported by funding from the Intramural Research Program of the National Institute on Aging/NIH. Work in JMV laboratory was supported by the Spanish Ministerio de Economía y Competitividad (MINECO) grant BFU2015-64630-R, Ministerio de Ciencia, Innovación y Universidades (MICIU) grant RTI2018-100695-B-I00, Spanish Junta de Andalucía grants P18-RT-4264, 1263735-R and BIO-276, the FEDER Funding Program from the European Union, and Universidad de Córdoba. MCR was supported by a FPU fellowship from the Spanish Ministerio de Educación, Cultura y Deporte (reference FPU14/06308). SRL held a FPI predoctoral contract funded by MINECO (reference BES-2016-078229).Peer reviewe