Regulation of mitochondrial respiration in heart cells analyzed by reaction-diffusion model of energy transfer

Higher AMPK activation in mouse oxidative compared with glycolytic muscle does not correlate with LKB1 or CaMKKβ expression

AMP-activated protein kinase (AMPK) is an energy-sensing serine/threonine kinase involved in metabolic regulation. It is phosphorylated by the upstream liver kinase B1 (LKB1) or calcium/calmodulin-dependent kinase kinase 2 (CaMKKβ). In cultured cells, AMPK activation correlates with LKB1 activity. The phosphorylation activates AMPK, shifting metabolism toward catabolism and promoting mitogenesis. In muscles, inactivity reduces AMPK activation, shifting the phenotype of oxidative muscles toward a more glycolytic profile. Here, we compared the basal level of AMPK activation in glycolytic and oxidative muscles and analyzed whether this relates to LKB1 or CaMKKβ. Using Western blotting, we assessed AMPK expression and phosphorylation in soleus, gastrocnemius (GAST), extensor digitorum longus (EDL), and heart from C57BL6J mice. We also assessed LKB1 and CaMKKβ expression, and CaMKKβ activity in tissue homogenates. AMPK activation was higher in oxidative (soleus and heart) than in glycolytic muscles (gastrocnemius and EDL). This correlated with AMPK α1-isoform expression, but not LKB1 and CaMKKβ. LKB1 expression was sex dependent and lower in male than female muscles. CaMKKβ expression was very low in skeletal muscles and did not phosphorylate AMPK in muscle lysates. The higher AMPK activation in oxidative muscles is in line with the fact that activated AMPK maintains an oxidative phenotype. However, this could not be explained by LKB1 and CaMKKβ. These results suggest that the regulation of AMPK activation is more complex in muscle than in cultured cells. As AMPK has been proposed as a therapeutic target for several diseases, future research should consider AMPK isoform expression and localization, and energetic compartmentalization.NEW & NOTEWORTHY It is important to understand how AMP-activated kinase, AMPK, is regulated, as it is a potential therapeutic target for several diseases. AMPK is activated by liver kinase B1, LKB1, and calcium/calmodulin-dependent kinase kinase 2, CaMKKβ. In cultured cells, AMPK activation correlates with LKB1 expression. In contrast, we show that AMPK-activation was higher in oxidative than glycolytic muscle, without correlating with LKB1 or CaMKKβ expression. Thus, AMPK regulation is more complex in highly compartmentalized muscle cells.

Compartmentalization in cardiomyocytes modulates creatine kinase and adenylate kinase activities

Intracellular molecules are transported by motor proteins or move by diffusion resulting from random molecular motion. Cardiomyocytes are packed with structures that are crucial for function, but also confine the diffusional spaces, providing cells with a means to control diffusion. They form compartments in which local concentrations are different from the overall, average concentrations. For example, calcium and cyclic AMP are highly compartmentalized, allowing these versatile second messengers to send different signals depending on their location. In energetic compartmentalization, the ratios of AMP and ADP to ATP are different from the average ratios. This is important for the performance of ATPases fuelling cardiac excitation-contraction coupling and mechanical work. A recent study suggested that compartmentalization modulates the activity of creatine kinase and adenylate kinase in situ. This could have implications for energetic signaling through, for example, AMP-activated kinase. It highlights the importance of taking compartmentalization into account in our interpretation of cellular physiology and developing methods to assess local concentrations of AMP and ADP to enhance our understanding of compartmentalization in different cell types.

Natural convection in the cytoplasm: Theoretical predictions of buoyancy-driven flows inside a cell

The existence of temperature gradients within eukaryotic cells has been postulated as a source of natural convection in the cytoplasm, i.e. bulk fluid motion as a result of temperature-difference-induced density gradients. Recent computations have predicted that a temperature differential of ΔT ≈ 1 K between the cell nucleus and the cell membrane could be strong enough to drive significant intracellular material transport. We use numerical computations and theoretical calculations to revisit this problem in order to further understand the impact of temperature gradients on flow generation and advective transport within cells. Surprisingly, our computations yield flows that are an order of magnitude weaker than those obtained previously for the same relative size and position of the nucleus with respect to the cell membrane. To understand this discrepancy, we develop a semi-analytical solution of the convective flow inside a model cell using a bi-spherical coordinate framework, for the case of an axisymmetric cell geometry (i.e. when the displacement of the nucleus from the cell centre is aligned with gravity). We also calculate exact solutions for the flow when the nucleus is located concentrically inside the cell. The results from both theoretical analyses agree with our numerical results, thus providing a robust estimate of the strength of cytoplasmic natural convection and demonstrating that these are much weaker than previously predicted. Finally, we investigate the ability of the aforementioned flows to redistribute solute within a cell. Our calculations reveal that, in all but unrealistic cases, cytoplasmic convection has a negligible contribution toward enhancing the diffusion-dominated mass transfer of cellular material.

Effects of altered cellular ultrastructure on energy metabolism in diabetic cardiomyopathy: an in silico study

Diabetic cardiomyopathy is a leading cause of heart failure in diabetes. At the cellular level, diabetic cardiomyopathy leads to altered mitochondrial energy metabolism and cardiomyocyte ultrastructure. We combined electron microscopy (EM) and computational modelling to understand the impact of diabetes-induced ultrastructural changes on cardiac bioenergetics. We collected transverse micrographs of multiple control and type I diabetic rat cardiomyocytes using EM. Micrographs were converted to finite-element meshes, and bioenergetics was simulated over them using a biophysical model. The simulations also incorporated depressed mitochondrial capacity for oxidative phosphorylation (OXPHOS) and creatine kinase (CK) reactions to simulate diabetes-induced mitochondrial dysfunction. Analysis of micrographs revealed a 14% decline in mitochondrial area fraction in diabetic cardiomyocytes, and an irregular arrangement of mitochondria and myofibrils. Simulations predicted that this irregular arrangement, coupled with the depressed activity of mitochondrial CK enzymes, leads to large spatial variation in adenosine diphosphate (ADP)/adenosine triphosphate (ATP) ratio profile of diabetic cardiomyocytes. However, when spatially averaged, myofibrillar ADP/ATP ratios of a cardiomyocyte do not change with diabetes. Instead, average concentration of inorganic phosphate rises by 40% owing to lower mitochondrial area fraction and dysfunction in OXPHOS. These simulations indicate that a disorganized cellular ultrastructure negatively impacts metabolite transport in diabetic cardiomyopathy. This article is part of the theme issue 'The cardiomyocyte: new revelations on the interplay between architecture and function in growth, health, and disease'.

Effect of crista morphology on mitochondrial ATP output: A computational study

Folding of the mitochondrial inner membrane (IM) into cristae greatly increases the ATP-generating surface area, S _IM, per unit volume but also creates diffusional bottlenecks that could limit reaction rates inside mitochondria. This study explores possible effects of inner membrane folding on mitochondrial ATP output, using a mathematical model for energy metabolism developed by the Jafri group and two- and three-dimensional spatial models for mitochondria, implemented on the Virtual Cell platform. Simulations demonstrate that cristae are micro-compartments functionally distinct from the cytosol. At physiological steady states, standing gradients of ADP form inside cristae that depend on the size and shape of the compartments, and reduce local flux (rate per unit area) of the adenine nucleotide translocase. This causes matrix ADP levels to drop, which in turn reduces the flux of ATP synthase. The adverse effects of membrane folding on reaction fluxes increase with crista length and are greater for lamellar than tubular crista. However, total ATP output per mitochondrion is the product of flux of ATP synthase and S _IM which can be two-fold greater for mitochondria with lamellar than tubular cristae, resulting in greater ATP output for the former. The simulations also demonstrate the crucial role played by intracristal kinases (adenylate kinase, creatine kinase) in maintaining the energy advantage of IM folding.

Altered calcium handling in cardiomyocytes from arginine-glycine amidinotransferase-knockout mice is rescued by creatine

The creatine kinase system facilitates energy transfer between mitochondria and the major ATPases in the heart. Creatine-deficient mice, which lack arginine-glycine amidinotransferase (AGAT) to synthesize creatine and homoarginine, exhibit reduced cardiac contractility. We studied how the absence of a functional CK system influences calcium handling in isolated cardiomyocytes from AGAT-knockouts and wild-type littermates as well as in AGAT-knockout mice receiving lifelong creatine supplementation via the food. Using a combination of whole cell patch clamp and fluorescence microscopy, we demonstrate that the L-type calcium channel (LTCC) current amplitude and voltage range of activation were significantly lower in AGAT-knockout compared with wild-type littermates. Additionally, the inactivation of LTCC and the calcium transient decay were significantly slower. According to our modeling results, these changes can be reproduced by reducing three parameters in knockout mice when compared with wild-type: LTCC conductance, the exchange constant of Ca2+ transfer between subspace and cytosol, and SERCA activity. Because tissue expression of LTCC and SERCA protein were not significantly different between genotypes, this suggests the involvement of posttranslational regulatory mechanisms or structural reorganization. The AGAT-knockout phenotype of calcium handling was fully reversed by dietary creatine supplementation throughout life. Our results indicate reduced calcium cycling in cardiomyocytes from AGAT-knockouts and suggest that the creatine kinase system is important for the development of calcium handling in the heart.NEW & NOTEWORTHY Creatine-deficient mice lacking arginine-glycine amidinotransferase exhibit compromised cardiac function. Here, we show that this is at least partially due to an overall slowing of calcium dynamics. Calcium influx into the cytosol via the L-type calcium current (LTCC) is diminished, and the rate of the sarcoendoplasmic reticulum calcium ATPase (SERCA) pumping calcium back into the sarcoplasmic reticulum is slower. The expression of LTCC and SERCA did not change, suggesting that the changes are regulatory.

Getting around the cell: physical transport in the intracellular world

Eukaryotic cells face the challenging task of transporting a variety of particles through the complex intracellular milieu in order to deliver, distribute, and mix the many components that support cell function. In this review, we explore the biological objectives and physical mechanisms of intracellular transport. Our focus is on cytoplasmic and intra-organelle transport at the whole-cell scale. We outline several key biological functions that depend on physically transporting components across the cell, including the delivery of secreted proteins, support of cell growth and repair, propagation of intracellular signals, establishment of organelle contacts, and spatial organization of metabolic gradients. We then review the three primary physical modes of transport in eukaryotic cells: diffusive motion, motor-driven transport, and advection by cytoplasmic flow. For each mechanism, we identify the main factors that determine speed and directionality. We also highlight the efficiency of each transport mode in fulfilling various key objectives of transport, such as particle mixing, directed delivery, and rapid target search. Taken together, the interplay of diffusion, molecular motors, and flows supports the intracellular transport needs that underlie a broad variety of biological phenomena.

Respiration of permeabilized cardiomyocytes from mice: no sex differences, but substrate-dependent changes in the apparent ADP-affinity

Sex differences in cardiac physiology are getting increased attention. This study assessed whether isolated, permeabilized cardiomyocytes from male and female C57BL/6 mice differ in terms of their respiration with multiple substrates and overall intracellular diffusion restriction estimated by the apparent ADP-affinity of respiration. Using respirometry, we recorded 1) the activities of respiratory complexes I, II and IV, 2) the respiration rate with substrates fuelling either complex I, II, or I + II, and 3) the apparent ADP-affinity with substrates fuelling complex I and I + II. The respiration rates were normalized to protein content and citrate synthase (CS) activity. We found no sex differences in CS activity (a marker of mitochondrial content) normalized to protein content or in any of the respiration measurements. This suggests that cardiomyocytes from male and female mice do not differ in terms of mitochondrial respiratory capacity and apparent ADP-affinity. Pyruvate modestly lowered the respiration rate, when added to succinate, glutamate and malate. This may be explained by intramitochondrial compartmentalization caused by the formation of supercomplexes and their association with specific dehydrogenases. To our knowledge, we show for the first time that the apparent ADP-affinity was substrate-dependent. This suggests that substrates may change or regulate intracellular barriers in cardiomyocytes.

Spatial control of neuronal metabolism through glucose-mediated mitochondrial transport regulation

Eukaryotic cells modulate their metabolism by organizing metabolic components in response to varying nutrient availability and energy demands. In rat axons, mitochondria respond to glucose levels by halting active transport in high glucose regions. We employ quantitative modeling to explore physical limits on spatial organization of mitochondria and localized metabolic enhancement through regulated stopping of processive motion. We delineate the role of key parameters, including cellular glucose uptake and consumption rates, that are expected to modulate mitochondrial distribution and metabolic response in spatially varying glucose conditions. Our estimates indicate that physiological brain glucose levels fall within the limited range necessary for metabolic enhancement. Hence mitochondrial localization is shown to be a plausible regulatory mechanism for neuronal metabolic flexibility in the presence of spatially heterogeneous glucose, as may occur in long processes of projection neurons. These findings provide a framework for the control of cellular bioenergetics through organelle trafficking.

Cells are equipped with power factories called mitochondria that turn nutrients into chemical energy to fuel processes in the cell. Hundreds of mitochondria move throughout the cell, shifting their positions in response to energy demands. This happens via molecular motors that pick the mitochondria up and carry them to new locations. Such movements enable the mitochondria to accumulate in parts of the cell with the greatest energy needs.

Mitochondria of nerve cells or neurons have a particular challenging job, as neurons can be very long and different parts within the cells can have different energy needs. It has been shown that mitochondria stop in regions where nutrients such as sugar are most concentrated. So far, it has been unclear whether this regulated stopping helps control energy balance in neurons.

Here, Agrawal et al. used a computational model of rat neurons to find out whether sugar levels are sufficient in guiding mitochondria. The results showed that the mitochondria only accumulated in high-nutrient regions when the sugar concentrations were moderate – not too low and not too high. A specific range of sugar levels was necessary to make this mechanism useful for increasing the efficiency of energy production. Such concentrations match the ones observed in healthy rat brains.

When neurons are unable to meet their energy demands, they stop working and sometimes even die. This is the case in many diseases, including diabetes, dementia, and Alzheimer’s disease. Computer models allow us to explore the complex energy regulation in detail. A better understanding of how neurons regulate their energy production and demand may help us discover how they become faulty in these diseases.

Insights on the impact of mitochondrial organisation on bioenergetics in high-resolution computational models of cardiac cell architecture

Recent electron microscopy data have revealed that cardiac mitochondria are not arranged in crystalline columns but are organised with several mitochondria aggregated into columns of varying sizes spanning the cell cross-section. This raises the question—how does the mitochondrial arrangement affect the metabolite distributions within cardiomyocytes and what is its impact on force dynamics? Here, we address this question by employing finite element modeling of cardiac bioenergetics on computational meshes derived from electron microscope images. Our results indicate that heterogeneous mitochondrial distributions can lead to significant spatial variation across the cell in concentrations of inorganic phosphate, creatine (Cr) and creatine phosphate (PCr). However, our model predicts that sufficient activity of the creatine kinase (CK) system, coupled with rapid diffusion of Cr and PCr, maintains near uniform ATP and ADP ratios across the cell cross sections. This homogenous distribution of ATP and ADP should also evenly distribute force production and twitch duration with contraction. These results suggest that the PCr shuttle and associated enzymatic reactions act to maintain uniform force dynamics in the cell despite the heterogeneous mitochondrial organization. However, our model also predicts that under hypoxia activity of mitochondrial CK enzymes and diffusion of high-energy phosphate compounds may be insufficient to sustain uniform ATP/ADP distribution and hence force generation.

Mammalian cardiomyocytes contain a high volume of mitochondria, which maintains the continuous and bulk supply of ATP to sustain normal heart function. Previously, cardiac mitochondria were understood to be distributed in a regular, crystalline pattern, which facilitated a steady supply of ATP at different workloads. Using electron microscopy images of cell cross sections, we recently found that they are not regularly distributed inside cardiomyocytes. We created new spatially accurate computational models of cardiac cell bioenergetics and tested whether this heterogeneous distribution of mitochondria causes non-uniform energy supply and contractile force production in the cardiomyocyte. We found that ATP and ADP concentrations remain uniform throughout the cell because of the activity of creatine kinase (CK) enzymes that convert ATP produced in the mitochondria into creatine phosphate. Creatine phosphate rapidly diffuses to the myofibril region where it can be converted back to ATP for the contraction cycle in a timely manner. This mechanism is called the phosphocreatine shuttle (PCr shuttle). The PCr shuttle ensures that different areas of the cell produce the same amount of force regardless of the mitochondrial distribution. However, our model also shows that when the cellular oxygen supply is limited—as can be the case in conditions such as heart failure—the PCr shuttle cannot maintain uniform ATP and ADP concentrations across the cell. This causes a non-uniform acto-myosin force distribution and non-uniform twitch duration across the cell cross section. Our study suggests that mechanisms other than the PCr shuttle may be necessary to maintain uniform supply of ATP in a hypoxic environment.

A computational study of the role of mitochondrial organization on cardiac bioenergetics

All cells in the body have a specific shape and internal organization which is specific to that cell's function. Heart cells are rod-shaped, and contain arrays of contractile protines (myofibrils) and mitochondria (organelles that produce energy) that are aligned along the length of the rod. This arrangement is presumed to allow the cell to generate maximal contractile force for each heartbeat and for energy metabolites to be readily available to generate this force. Heart disease phenotypes, such as diabetic cardiomyopathy and heart failure, exhibit altered organization of mitochondria. However, physiological and computational studies have predominantly investigated the effect of the biochemical changes that accompany the disease alone, such as reduced rates of ATP production by mitochondria. We present a modeling study that examines the effect of mitochondrial organization on energy metabolite distribution during the heartbeat. A 2D micrograph of the cell cross-section was selected from a 3D image stack of structural data of a cardiac cell. The image was used to generate a 2D finite element model, on which mitochondrial oxidative phosphorylation and energy metabolite diffusion was modelled. Results illustrate that mitochondrial density can induce heterogeneity in the distribution of metabolites across the cell area. We discuss the implications of these findings and avenues for future, more indepth studies.

Changes in mitochondrial morphology and organization can enhance energy supply from mitochondrial oxidative phosphorylation in diabetic cardiomyopathy

Diabetic cardiomyopathy is accompanied by metabolic and ultrastructural alterations, but the impact of the structural changes on metabolism itself is yet to be determined. Morphometric analysis of mitochondrial shape and spatial organization within transverse sections of cardiomyocytes from control and streptozotocin-induced type I diabetic Sprague-Dawley rats revealed that mitochondria are 20% smaller in size while their spatial density increases by 53% in diabetic cells relative to control myocytes. Diabetic cells formed larger clusters of mitochondria (60% more mitochondria per cluster) and the effective surface-to-volume ratio of these clusters increased by 22.5%. Using a biophysical computational model we found that this increase can have a moderate compensatory effect by increasing the availability of ATP in the cytosol when ATP synthesis within the mitochondrial matrix is compromised.

Metabolic compartmentation in rainbow trout cardiomyocytes: coupling of hexokinase but not creatine kinase to mitochondrial respiration

Rainbow trout (Oncorhynchus mykiss) cardiomyocytes have a simple morphology with fewer membrane structures such as sarcoplasmic reticulum and t-tubules penetrating the cytosol. Despite this, intracellular ADP diffusion is restricted. Intriguingly, although diffusion is restricted, trout cardiomyocytes seem to lack the coupling between mitochondrial creatine kinase (CK) and respiration. Our aim was to study the distribution of diffusion restrictions in permeabilized trout cardiomyocytes and verify the role of CK. We found a high activity of hexokinase (HK), which led us to reassess the situation in trout cardiomyocytes. We show that diffusion restrictions are more prominent than previously thought. In the presence of a competitive ADP-trapping system, ADP produced by HK, but not CK, was channeled to the mitochondria. In agreement with this, we found no positively charged mitochondrial CK in trout heart homogenate. The results were best fit by a simple mathematical model suggesting that trout cardiomyocytes lack a functional coupling between ATPases and pyruvate kinase. The model simulations show that diffusion is restricted to almost the same extent in the cytosol and by the outer mitochondrial membrane. Furthermore, they confirm that HK, but not CK, is functionally coupled to respiration. In perspective, our results suggest that across a range of species, cardiomyocyte morphology and metabolism go hand in hand with cardiac performance, which is adapted to the circumstances. Mitochondrial CK is coupled to respiration in adult mammalian hearts, which are specialized to high, sustained performance. HK associates with mitochondria in hearts of trout and neonatal mammals, which are more hypoxia-tolerant.

Bioenergetics of the aging heart and skeletal muscles: Modern concepts and controversies

Age-related alterations in the bioenergetics of the heart and oxidative skeletal muscle tissues are of crucial influence on their performance. Until now the prevailing concept of aging was the mitochondrial theory, the increased production of reactive oxygen species, mediated by deficiency in the activity of respiratory chain complexes. However, studies with mitochondria in situ have presented results which, to some extent, disagree with previous ones, indicating that the mitochondrial theory of aging may be overestimated. The studies reporting age-related decline in mitochondrial function were performed using mainly isolated mitochondria. Measurements on this level are not able to take into account the system level properties. The relevant information can be obtained only from appropriate studies using cells or tissue fibers. The functional interactions between the components of Intracellular Energetic Unit (ICEU) regulate the energy production and consumption in oxidative muscle cells. The alterations of these interactions in ICEU should be studied in order to find a more effective protocol to decelerate the age-related changes taking place in the energy metabolism. In this article, an overview is given of the present theories and controversies of causes of age-related alterations in bioenergetics. Also, branches of study, which need more emphasis, are indicated.

Impaired ADP channeling to mitochondria and elevated reactive oxygen species in hypertensive hearts

Systemic hypertension initially promotes a compensatory cardiac hypertrophy, yet it progresses to heart failure (HF), and energetic deficits appear to be central to this failure. However, the transfer of energy between the mitochondria and the myofibrils is not often considered as part of the energetic equation. We compared hearts from old spontaneously hypertensive rats (SHRs) and normotensive Wistar controls. SHR hearts showed a 35% depression in mitochondrial function, yet produced at least double the amount of reactive oxygen species (ROS) in all respiration states in left ventricular (LV) homogenates. To test the connectivity between mitochondria and myofibrils, respiration was further tested in situ with LV permeabilized fibers by addition of multiple substrates and ATP, which requires hydrolysis to mediate oxidative phosphorylation. By trapping ADP using a pyruvate kinase enzyme system, we tested ADP channeling towards mitochondria, and this suppressed respiration and elevated ROS production more in the SHR fibers. The ADP-trapped state was also less relieved on creatine addition, likely reflecting the 30% depression in total CK activity in the SHR heart fibers. Confocal imaging identified a 34% longer distance between the centers of myofibril to mitochondria in the SHR hearts, which increases transverse metabolite diffusion distances (e.g., for ATP, ADP, and creatine phosphate). We propose that impaired connectivity between mitochondria and myofibrils may contribute to elevated ROS production. Impaired energy exchange could be the result of ultrastructural changes that occur with hypertrophy in this model of hypertension.

Restricted ADP movement in cardiomyocytes: Cytosolic diffusion obstacles are complemented with a small number of open mitochondrial voltage-dependent anion channels

Adequate intracellular energy transfer is crucial for proper cardiac function. In energy starved failing hearts, partial restoration of energy transfer can rescue mechanical performance. There are two types of diffusion obstacles that interfere with energy transfer from mitochondria to ATPases: mitochondrial outer membrane (MOM) with voltage-dependent anion channel (VDAC) permeable to small hydrophilic molecules and cytoplasmatic diffusion barriers grouping ATP-producers and -consumers. So far, there is no method developed to clearly distinguish the contributions of cytoplasmatic barriers and MOM to the overall diffusion restriction. Furthermore, the number of open VDACs in vivo remains unknown. The aim of this work was to establish the partitioning of intracellular diffusion obstacles in cardiomyocytes. We studied the response of mitochondrial oxidative phosphorylation of permeabilized rat cardiomyocytes to changes in extracellular ADP by recording 3D image stacks of NADH autofluorescence. Using cell-specific mathematical models, we determined the permeability of MOM and cytoplasmatic barriers. We found that only ~2% of VDACs are accessible to cytosolic ADP and cytoplasmatic diffusion barriers reduce the apparent diffusion coefficient by 6-10×. In cardiomyocytes, diffusion barriers in the cytoplasm and by the MOM restrict ADP/ATP diffusion to similar extents suggesting a major role of both barriers in energy transfer and other intracellular processes.

The impact of cardiac ischemia/reperfusion on the mitochondria-cytoskeleton interactions

Cardiac ischemia-reperfusion (IR) injury compromises mitochondrial oxidative phosphorylation (OxPhos) and compartmentalized intracellular energy transfer via the phosphocreatine/creatine kinase (CK) network. The restriction of ATP/ADP diffusion at the level of the mitochondrial outer membrane (MOM) is an essential element of compartmentalized energy transfer. In adult cardiomyocytes, the MOM permeability to ADP is regulated by the interaction of voltage-dependent anion channel with cytoskeletal proteins, particularly with β tubulin II. The IR-injury alters the expression and the intracellular arrangement of cytoskeletal proteins. The objective of the present study was to investigate the impact of IR on the intracellular arrangement of β tubulin II and its effect on the regulation of mitochondrial respiration. Perfused rat hearts were subjected to total ischemia (for 20min (I20) and 45min (I45)) or to ischemia followed by 30min of reperfusion (I20R and I45R groups). High resolution respirometry and fluorescent confocal microscopy were used to study respiration, β tubulin II and mitochondrial arrangements in cardiac fibers. The results of these experiments evidence a heterogeneous response of mitochondria to IR-induced damage. Moreover, the intracellular rearrangement of β tubulin II, which in the control group colocalized with mitochondria, was associated with increased apparent affinity of OxPhos for ADP, decreased regulation of respiration by creatine without altering mitochondrial CK activity and the ratio between octameric to dimeric isoenzymes. The results of this study allow us to highlight changes of mitochondrial interactions with cytoskeleton as one of the possible mechanisms underlying cardiac IR injury.

Distinct functional roles of cardiac mitochondrial subpopulations revealed by a 3D simulation model

Experimental characterization of two cardiac mitochondrial subpopulations, namely, subsarcolemmal mitochondria (SSM) and interfibrillar mitochondria (IFM), has been hampered by technical difficulties, and an alternative approach is eagerly awaited. We previously developed a three-dimensional computational cardiomyocyte model that integrates electrophysiology, metabolism, and mechanics with subcellular structure. In this study, we further developed our model to include intracellular oxygen diffusion, and determined whether mitochondrial localization or intrinsic properties cause functional variations. For this purpose, we created two models: one with equal SSM and IFM properties and one with IFM having higher activity levels. Using these two models to compare the SSM and IFM responses of [Ca(2+)], tricarboxylic acid cycle activity, [NADH], and mitochondrial inner membrane potential to abrupt changes in pacing frequency (0.25-2 Hz), we found that the reported functional differences between these subpopulations appear to be mostly related to local [Ca(2+)] heterogeneity, and variations in intrinsic properties only serve to augment these differences. We also examined the effect of hypoxia on mitochondrial function. Under normoxic conditions, intracellular oxygen is much higher throughout the cell than the half-saturation concentration for oxidative phosphorylation. However, under limited oxygen supply, oxygen is mostly exhausted in SSM, leaving the core region in an anoxic condition. Reflecting this heterogeneous oxygen environment, the inner membrane potential continues to decrease in IFM, whereas it is maintained to nearly normal levels in SSM, thereby ensuring ATP supply to this region. Our simulation results provide clues to understanding the origin of functional variations in two cardiac mitochondrial subpopulations and their differential roles in maintaining cardiomyocyte function as a whole.

Regulation of cardiac cellular bioenergetics: mechanisms and consequences

The regulation of cardiac cellular bioenergetics is critical for maintaining normal cell function, yet the nature of this regulation is not fully understood. Different mechanisms have been proposed to explain how mitochondrial ATP production is regulated to match changing cellular energy demand while metabolite concentrations are maintained. We have developed an integrated mathematical model of cardiac cellular bioenergetics, electrophysiology, and mechanics to test whether stimulation of the dehydrogenase flux by Ca²⁺ or Pi, or stimulation of complex III by Pi can increase the rate of mitochondrial ATP production above that determined by substrate availability (ADP and Pi). Using the model, we show that, under physiological conditions the rate of mitochondrial ATP production can match varying demand through substrate availability alone; that ATP production rate is not limited by the supply of reducing equivalents in the form of NADH, as a result of Ca²⁺ or Pi activation of the dehydrogenases; and that ATP production rate is sensitive to feedback activation of complex III by Pi. We then investigate the mechanistic implications on cytosolic ion homeostasis and force production by simulating the concentrations of cytosolic Ca²⁺, Na⁺ and K⁺, and activity of the key ATPases, SERCA pump, Na⁺/K⁺ pump and actin-myosin ATPase, in response to increasing cellular energy demand. We find that feedback regulation of mitochondrial complex III by Pi improves the coupling between energy demand and mitochondrial ATP production and stabilizes cytosolic ADP and Pi concentrations. This subsequently leads to stabilized cytosolic ionic concentrations and consequentially reduced energetic cost from cellular ATPases.

In silico analysis of exercise intolerance in myalgic encephalomyelitis/chronic fatigue syndrome

Post-exertional malaise is commonly observed in patients with myalgic encephalomyelitis/chronic fatigue syndrome, but its mechanism is not yet well understood. A reduced capacity for mitochondrial ATP synthesis is associated with the pathogenesis of CFS and is suspected to be a major contribution to exercise intolerance in CFS patients. To demonstrate the connection between a reduced mitochondrial capacity and exercise intolerance, we present a model which simulates metabolite dynamics in skeletal muscles during exercise and recovery. CFS simulations exhibit critically low levels of ATP, where an increased rate of cell death would be expected. To stabilize the energy supply at low ATP concentrations the total adenine nucleotide pool is reduced substantially causing a prolonged recovery time even without consideration of other factors, such as immunological dysregulations and oxidative stress. Repeated exercises worsen this situation considerably. Furthermore, CFS simulations exhibited an increased acidosis and lactate accumulation consistent with experimental observations.

'Idealized' state 4 and state 3 in mitochondria vs. rest and work in skeletal muscle

A computer model of oxidative phosphorylation (OXPHOS) in skeletal muscle is used to compare state 3, intermediate state and state 4 in mitochondria with rest and work in skeletal muscle. 'Idealized' state 4 and 3 in relation to various 'experimental' states 4 and 3 are defined. Theoretical simulations show, in accordance with experimental data, that oxygen consumption (V'O2), ADP and Pi are higher, while ATP/ADP and Δp are lower in rest than in state 4, because of the presence of basal ATP consuming reactions in the former. It is postulated that moderate and intensive work in skeletal muscle is very different from state 3 in isolated mitochondria. V'O2, ATP/ADP, Δp and the control of ATP usage over V'O2 are much higher, while ADP and Pi are much lower in the former. The slope of the phenomenological V'O2-ADP relationship is much steeper during the rest-work transition than during the state 4-state 3 transition. The work state in intact muscle is much more similar to intermediate state than to state 3 in isolated mitochondria in terms of ADP, ATP/ADP, Δp and metabolic control pattern, but not in terms of V'O2. The huge differences between intact muscle and isolated mitochondria are proposed to be caused by the presence of the each-step activation (ESA) mechanism of the regulation of OXPHOS in intact skeletal muscle. Generally, the present study suggests that isolated mitochondria (at least in the absence of Ca2+) cannot serve as a good model of OXPHOS regulation in intact skeletal muscle.

Modular organization of cardiac energy metabolism: energy conversion, transfer and feedback regulation

To meet high cellular demands, the energy metabolism of cardiac muscles is organized by precise and coordinated functioning of intracellular energetic units (ICEUs). ICEUs represent structural and functional modules integrating multiple fluxes at sites of ATP generation in mitochondria and ATP utilization by myofibrillar, sarcoplasmic reticulum and sarcolemma ion-pump ATPases. The role of ICEUs is to enhance the efficiency of vectorial intracellular energy transfer and fine tuning of oxidative ATP synthesis maintaining stable metabolite levels to adjust to intracellular energy needs through the dynamic system of compartmentalized phosphoryl transfer networks. One of the key elements in regulation of energy flux distribution and feedback communication is the selective permeability of mitochondrial outer membrane (MOM) which represents a bottleneck in adenine nucleotide and other energy metabolite transfer and microcompartmentalization. Based on the experimental and theoretical (mathematical modelling) arguments, we describe regulation of mitochondrial ATP synthesis within ICEUs allowing heart workload to be linearly correlated with oxygen consumption ensuring conditions of metabolic stability, signal communication and synchronization. Particular attention was paid to the structure-function relationship in the development of ICEU, and the role of mitochondria interaction with cytoskeletal proteins, like tubulin, in the regulation of MOM permeability in response to energy metabolic signals providing regulation of mitochondrial respiration. Emphasis was given to the importance of creatine metabolism for the cardiac energy homoeostasis.

Molecular and subcellular-scale modeling of nucleotide diffusion in the cardiac myofilament lattice

Contractile function of cardiac cells is driven by the sliding displacement of myofilaments powered by the cycling myosin crossbridges. Critical to this process is the availability of ATP, which myosin hydrolyzes during the cross-bridge cycle. The diffusion of adenine nucleotides through the myofilament lattice has been shown to be anisotropic, with slower radial diffusion perpendicular to the filament axis relative to parallel, and is attributed to the periodic hexagonal arrangement of the thin (actin) and thick (myosin) filaments. We investigated whether atomistic-resolution details of myofilament proteins can refine coarse-grain estimates of diffusional anisotropy for adenine nucleotides in the cardiac myofibril, using homogenization theory and atomistic thin filament models from the Protein Data Bank. Our results demonstrate considerable anisotropy in ATP and ADP diffusion constants that is consistent with experimental measurements and dependent on lattice spacing and myofilament overlap. A reaction-diffusion model of the half-sarcomere further suggests that diffusional anisotropy may lead to modest adenine nucleotide gradients in the myoplasm under physiological conditions.

ADP protects cardiac mitochondria under severe oxidative stress

ADP is not only a key substrate for ATP generation, but also a potent inhibitor of mitochondrial permeability transition pore (mPTP). In this study, we assessed how oxidative stress affects the potency of ADP as an mPTP inhibitor and whether its reduction of reactive oxygen species (ROS) production might be involved. We determined quantitatively the effects of ADP on mitochondrial Ca(2+) retention capacity (CRC) until the induction of mPTP in normal and stressed isolated cardiac mitochondria. We used two models of chronic oxidative stress (old and diabetic mice) and two models of acute oxidative stress (ischemia reperfusion (IR) and tert-butyl hydroperoxide (t-BH)). In control mitochondria, the CRC was 344 ± 32 nmol/mg protein. 500 μmol/L ADP increased CRC to 774 ± 65 nmol/mg protein. This effect of ADP seemed to relate to its concentration as 50 μmol/L had a significantly smaller effect. Also, oligomycin, which inhibits the conversion of ADP to ATP by F0F1ATPase, significantly increased the effect of 50 μmol/L ADP. Chronic oxidative stress did not affect CRC or the effect of 500 μmol/L ADP. After IR or t-BH exposure, CRC was drastically reduced to 1 ± 0.2 and 32 ± 4 nmol/mg protein, respectively. Surprisingly, ADP increased the CRC to 447 ± 105 and 514 ± 103 nmol/mg protein in IR and t-BH, respectively. Thus, it increased CRC by the same amount as in control. In control mitochondria, ADP decreased both substrate and Ca(2+)-induced increase of ROS. However, in t-BH mitochondria the effect of ADP on ROS was relatively small. We conclude that ADP potently restores CRC capacity in severely stressed mitochondria. This effect is most likely not related to a reduction in ROS production. As the effect of ADP relates to its concentration, increased ADP as occurs in the pathophysiological situation may protect mitochondrial integrity and function.

Regulation of oxidative phosphorylation during work transitions results from its kinetic properties

The regulation of oxidative phosphorylation (OXPHOS) during work transitions in skeletal muscle and heart is still not well understood. Different computer models of this process have been developed that are characterized by various kinetic properties. In the present research-polemic theoretical study it is argued that models belonging to one group (Model A), which predict that among OXPHOS complexes complex III keeps almost all of the metabolic control over oxygen consumption (Vo2) and involve a strong complex III activation by inorganic phosphate (Pi), lead to the conclusion that an increase in Pi is the main mechanism responsible for OXPHOS activation (feedback-activation mechanism). Models belonging to another group (Model B), which were developed to take into account an approximately uniform distribution of metabolic control over Vo2 among particular OXPHOS complexes (complex I, complex III, complex IV, ATP synthase, ATP/ADP carrier, phosphate carrier) encountered in experimental studies in isolated mitochondria, predict that all OXPHOS complexes are directly activated in parallel with ATP usage and NADH supply by some external cytosolic factor/mechanism during rest-to-work or low-to-high work transitions in skeletal muscle and heart ("each-step-activation" mechanism). Model B demonstrates that different intensities of each-step activation can account for the very different (slopes of) phenomenological Vo2-ADP relationships observed in various skeletal muscles and heart. Thus they are able to explain the differences in the regulation of OXPHOS during work transitions between skeletal muscle (where moderate changes in ADP take place) and intact heart in vivo (where ADP is essentially constant).

On-off asymmetries in oxygen consumption kinetics of single Xenopus laevis skeletal muscle fibres suggest higher-order control

The mechanisms controlling skeletal muscle oxygen consumption (V(o)₂) during exercise are not well understood. We determined whether first-order control could explain V(o)₂kinetics at contractions onset (V(o)₂(on)) and cessation (V(o)₂off)) in single skeletal muscle fibres differing in oxdidative capacity, and across stimulation intensities up to V(o)₂(max). Xenopus laevis fibres (n = 21) were suspended in a sealed chamber with a fast response P(o)₂ electrode to measure V(o)₂ every second before, during and after stimulated isometric contractions. A first-order model did not well characterize on-transient V(o)₂ kinetics. Including a time delay (TD) in the model provided a significantly improved characterization than a first-order fit without TD (F-ratio; P < 0.05), and revealed separate 'activation' and 'exponential' phases in 15/21 fibres contracting at V(o)₂(max) (mean ± SD TD: 14 ± 3s). On-transient kinetics (τV(o)₂(on)) was weakly and linearly related to V(o)₂(max) (R² = 0.271, P = 0.015). Off-transient kinetics, however, were first-order, and τV(o)₂(off) was greater in low-oxidative (V(o)₂max < 0.05 nmol mm⁻³s⁻¹ than high-oxidative fibres (V(o)₂(max > 0.10 nmol mm ⁻³ s⁻¹; 170 ± 70 vs. 29 ± 6 s, P < 0.001). 1/ τV(o)₂(off) was proportional to V(o)₂(max) (R² = 0.727, P < 0.001), unlike in the on-transient. The calculated oxygen deficit was larger (P < 0.05) than the post-contraction volume of consumed oxygen at all intensities except V(o)₂(max). These data show a clear dissociation between the kinetic control of V(o)₂at the onset and cessation of contractions and across stimulation intensities. More complex models are therefore required to understand the activation of mitochondrial respiration in skeletal muscle at the start of exercise.

Computational modeling of mitochondrial energy transduction

Mitochondria are the power plant of the heart, burning fat and sugars to supply the muscle with the adenosine triphosphate (ATP) free energy that drives contraction and relaxation during each heart beat. This function was first captured in a mathematical model in 1967. Today, interest in such a model has been rekindled by ongoing in silico integrative physiology efforts such as the Cardiac Physiome project. Here, the status of the field of computational modeling of mitochondrial ATP synthetic function is reviewed.

Molecular dynamics simulations of creatine kinase and adenine nucleotide translocase in mitochondrial membrane patch

Interaction between mitochondrial creatine kinase (MtCK) and adenine nucleotide translocase (ANT) can play an important role in determining energy transfer pathways in the cell. Although the functional coupling between MtCK and ANT has been demonstrated, the precise mechanism of the coupling is not clear. To study the details of the coupling, we turned to molecular dynamics simulations. We introduce a new coarse-grained molecular dynamics model of a patch of the mitochondrial inner membrane containing a transmembrane ANT and an MtCK above the membrane. The membrane model consists of three major types of lipids (phosphatidylcholine, phosphatidylethanolamine, and cardiolipin) in a roughly 2:1:1 molar ratio. A thermodynamics-based coarse-grained force field, termed MARTINI, has been used together with the GROMACS molecular dynamics package for all simulated systems in this work. Several physical properties of the system are reproduced by the model and are in agreement with known data. This includes membrane thickness, dimension of the proteins, and diffusion constants. We have studied the binding of MtCK to the membrane and demonstrated the effect of cardiolipin on the stabilization of the binding. In addition, our simulations predict which part of the MtCK protein sequence interacts with the membrane. Taken together, the model has been verified by dynamical and structural data and can be used as the basis for further studies.

Regulation of respiration in muscle cells in vivo by VDAC through interaction with the cytoskeleton and MtCK within Mitochondrial Interactosome

This review describes the recent experimental data on the importance of the VDAC-cytoskeleton interactions in determining the mechanisms of energy and metabolite transfer between mitochondria and cytoplasm in cardiac cells. In the intermembrane space mitochondrial creatine kinase connects VDAC with adenine nucleotide translocase and ATP synthase complex, on the cytoplasmic side VDAC is linked to cytoskeletal proteins. Applying immunofluorescent imaging and Western blot analysis we have shown that β2-tubulin coexpressed with mitochondria is highly important for cardiac muscle cells mitochondrial metabolism. Since it has been shown by Rostovtseva et al. that αβ-heterodimer of tubulin binds to VDAC and decreases its permeability, we suppose that the β-tubulin subunit is bound on the cytoplasmic side and α-tubulin C-terminal tail is inserted into VDAC. Other cytoskeletal proteins, such as plectin and desmin may be involved in this process. The result of VDAC-cytoskeletal interactions is selective restriction of the channel permeability for adenine nucleotides but not for creatine or phosphocreatine that favors energy transfer via the phosphocreatine pathway. In some types of cancer cells these interactions are altered favoring the hexokinase binding and thus explaining the Warburg effect of increased glycolytic lactate production in these cells. This article is part of a Special Issue entitled: VDAC structure, function, and regulation of mitochondrial metabolism.

Molecular system bioenergics of the heart: experimental studies of metabolic compartmentation and energy fluxes versus computer modeling

In this review we analyze the recent important and remarkable advancements in studies of compartmentation of adenine nucleotides in muscle cells due to their binding to macromolecular complexes and cellular structures, which results in non-equilibrium steady state of the creatine kinase reaction. We discuss the problems of measuring the energy fluxes between different cellular compartments and their simulation by using different computer models. Energy flux determinations by ¹⁸O transfer method have shown that in heart about 80% of energy is carried out of mitochondrial intermembrane space into cytoplasm by phosphocreatine fluxes generated by mitochondrial creatine kinase from adenosine triphosphate (ATP), produced by ATP Synthasome. We have applied the mathematical model of compartmentalized energy transfer for analysis of experimental data on the dependence of oxygen consumption rate on heart workload in isolated working heart reported by Williamson et al. The analysis of these data show that even at the maximal workloads and respiration rates, equal to 174 μmol O₂ per min per g dry weight, phosphocreatine flux, and not ATP, carries about 80–85% percent of energy needed out of mitochondria into the cytosol. We analyze also the reasons of failures of several computer models published in the literature to correctly describe the experimental data.

Permeabilized rat cardiomyocyte response demonstrates intracellular origin of diffusion obstacles

Intracellular diffusion restrictions for ADP and other molecules have been predicted earlier based on experiments on permeabilized fibers or cardiomyocytes. However, it is possible that the effective diffusion distance is larger than the cell dimensions due to clumping of cells and incomplete separation of cells in fiber preparations. The aim of this work was to check whether diffusion restrictions exist inside rat cardiomyocytes or are caused by large effective diffusion distance. For that, we determined the response of oxidative phosphorylation (OxPhos) to exogenous ADP and ATP stimulation in permeabilized rat cardiomyocytes using fluorescence microscopy. The state of OxPhos was monitored via NADH and flavoprotein autofluorescence. By varying the ADP or ATP concentration in flow chamber, we determined that OxPhos has a low affinity in cardiomyocytes. The experiments were repeated in a fluorometer on cardiomyocyte suspensions leading to similar autofluorescence changes induced by ADP as recorded under the microscope. ATP stimulated OxPhos more in a fluorometer than under the microscope, which was attributed to accumulation of ADP in fluorometer chamber. By calculating the flow profile around the cell in the microscope chamber and comparing model solutions to measured data, we demonstrate that intracellular structures impose significant diffusion obstacles in rat cardiomyocytes.

Analyzing the functional properties of the creatine kinase system with multiscale 'sloppy' modeling

In this study the function of the two isoforms of creatine kinase (CK; EC 2.7.3.2) in myocardium is investigated. The 'phosphocreatine shuttle' hypothesis states that mitochondrial and cytosolic CK plays a pivotal role in the transport of high-energy phosphate (HEP) groups from mitochondria to myofibrils in contracting muscle. Temporal buffering of changes in ATP and ADP is another potential role of CK. With a mathematical model, we analyzed energy transport and damping of high peaks of ATP hydrolysis during the cardiac cycle. The analysis was based on multiscale data measured at the level of isolated enzymes, isolated mitochondria and on dynamic response times of oxidative phosphorylation measured at the whole heart level. Using 'sloppy modeling' ensemble simulations, we derived confidence intervals for predictions of the contributions by phosphocreatine (PCr) and ATP to the transfer of HEP from mitochondria to sites of ATP hydrolysis. Our calculations indicate that only 15±8% (mean±SD) of transcytosolic energy transport is carried by PCr, contradicting the PCr shuttle hypothesis. We also predicted temporal buffering capabilities of the CK isoforms protecting against high peaks of ATP hydrolysis (3750 µM*s(-1)) in myofibrils. CK inhibition by 98% in silico leads to an increase in amplitude of mitochondrial ATP synthesis pulsation from 215±23 to 566±31 µM*s(-1), while amplitudes of oscillations in cytosolic ADP concentration double from 77±11 to 146±1 µM. Our findings indicate that CK acts as a large bandwidth high-capacity temporal energy buffer maintaining cellular ATP homeostasis and reducing oscillations in mitochondrial metabolism. However, the contribution of CK to the transport of high-energy phosphate groups appears limited. Mitochondrial CK activity lowers cytosolic inorganic phosphate levels while cytosolic CK has the opposite effect.

Intracellular Energetic Units regulate metabolism in cardiac cells

This review describes developments in historical perspective as well as recent results of investigations of cellular mechanisms of regulation of energy fluxes and mitochondrial respiration by cardiac work - the metabolic aspect of the Frank-Starling law of the heart. A Systems Biology solution to this problem needs the integration of physiological and biochemical mechanisms that take into account intracellular interactions of mitochondria with other cellular systems, in particular with cytoskeleton components. Recent data show that different tubulin isotypes are involved in the regular arrangement exhibited by mitochondria and ATP-consuming systems into Intracellular Energetic Units (ICEUs). Beta II tubulin association with the mitochondrial outer membrane, when co-expressed with mitochondrial creatine kinase (MtCK) specifically limits the permeability of voltage-dependent anion channel for adenine nucleotides. In the MtCK reaction this interaction changes the regulatory kinetics of respiration through a decrease in the affinity for adenine nucleotides and an increase in the affinity for creatine. Metabolic Control Analysis of the coupled MtCK-ATP Synthasome in permeabilized cardiomyocytes showed a significant increase in flux control by steps involved in ADP recycling. Mathematical modeling of compartmentalized energy transfer represented by ICEUs shows that cyclic changes in local ADP, Pi, phosphocreatine and creatine concentrations during contraction cycle represent effective metabolic feedback signals when amplified in the coupled non-equilibrium MtCK-ATP Synthasome reactions in mitochondria. This mechanism explains the regulation of respiration on beat to beat basis during workload changes under conditions of metabolic stability. This article is part of a Special Issue entitled "Local Signaling in Myocytes."

Integrative modeling of the cardiac ventricular myocyte

Cardiac electrophysiology is a discipline with a rich 50-year history of experimental research coupled with integrative modeling which has enabled us to achieve a quantitative understanding of the relationships between molecular function and the integrated behavior of the cardiac myocyte in health and disease. In this paper, we review the development of integrative computational models of the cardiac myocyte. We begin with a historical overview of key cardiac cell models that helped shape the field. We then narrow our focus to models of the cardiac ventricular myocyte and describe these models in the context of their subcellular functional systems including dynamic models of voltage-gated ion channels, mitochondrial energy production, ATP-dependent and electrogenic membrane transporters, intracellular Ca dynamics, mechanical contraction, and regulatory signal transduction pathways. We describe key advances and limitations of the models as well as point to new directions for future modeling research. WIREs Syst Biol Med 2011 3 392-413 DOI: 10.1002/wsbm.122

Systems bioenergetics of creatine kinase networks: physiological roles of creatine and phosphocreatine in regulation of cardiac cell function

Physiological role of creatine (Cr) became first evident in the experiments of Belitzer and Tsybakova in 1939, who showed that oxygen consumption in a well-washed skeletal muscle homogenate increases strongly in the presence of creatine and with this results in phosphocreatine (PCr) production with PCr/O(2) ratio of about 5-6. This was the beginning of quantitative analysis in bioenergetics. It was also observed in many physiological experiments that the contractile force changes in parallel with the alteration in the PCr content. On the other hand, it was shown that when heart function is governed by Frank-Starling law, work performance and oxygen consumption rate increase in parallel without any changes in PCr and ATP tissue contents (metabolic homeostasis). Studies of cellular mechanisms of all these important phenomena helped in shaping new approach to bioenergetics, Molecular System Bioenergetics, a part of Systems Biology. This approach takes into consideration intracellular interactions that lead to novel mechanisms of regulation of energy fluxes. In particular, interactions between mitochondria and cytoskeleton resulting in selective restriction of permeability of outer mitochondrial membrane anion channel (VDAC) for adenine nucleotides and thus their recycling in mitochondria coupled to effective synthesis of PCr by mitochondrial creatine kinase, MtCK. Therefore, Cr concentration and the PCr/Cr ratio became important kinetic parameters in the regulation of respiration and energy fluxes in muscle cells. Decrease in the intracellular contents of Cr and PCr results in a hypodynamic state of muscle and muscle pathology. Many experimental studies have revealed that PCr may play two important roles in the regulation of muscle energetics: first by maintaining local ATP pools via compartmentalized creatine kinase reactions, and secondly by stabilizing cellular membranes due to electrostatic interactions with phospholipids. The second mechanism decreases the production of lysophosphoglycerides in hypoxic heart, protects the cardiac cells sarcolemma against ischemic damage, decreases the frequency of arrhythmias and increases the post-ischemic recovery of contractile function. PCr is used as a pharmacological product Neoton in cardiac surgery as one of the components of cardioplegic solutions for protection of the heart against intraoperational injury and injected intravenously in acute myocardial ischemic conditions for improving the hemodynamic response and clinical conditions of patients with heart failure.

Virtual mitochondrion: towards an integrated model of oxidative phosphorylation complexes and beyond

The modelling of OXPHOS (oxidative phosphorylation) in order to integrate all kinetic and thermodynamic aspects of chemiosmotic theory has a long history. We briefly review this history and show how new ways of modelling are required to integrate a local model of the individual respiratory complexes into a global model of OXPHOS and, beyond that, into a reliable overall model of central metabolism.

Modulation of energy transfer pathways between mitochondria and myofibrils by changes in performance of perfused heart

In the heart, the energy supplied by mitochondria to myofibrils is continuously and finely tuned to the contraction requirement over a wide range of cardiac loads. This process is mediated both by the creatine kinase (CK) shuttle and by direct ATP transfer. The aim of this study was to identify the contribution of energy transfer pathways at different cardiac performance levels. For this, five protocols of (31)P NMR inversion and saturation transfer experiments were performed at different performance levels on Langendorff perfused rat hearts. The cardiac performance was changed either through variation of external calcium in the presence or absence of isoprenaline or through variation of LV balloon inflation. The recordings were analyzed by mathematical models composed on the basis of different energy transfer pathway configurations. According to our results, the total CK unidirectional flux was relatively stable when the cardiac performance was changed by increasing the calcium concentration or variation of LV balloon volume. The stability of total CK unidirectional flux is lost at extreme energy demand levels leading to a rise in inorganic phosphate, a drop of ATP and phosphocreatine, a drop of total CK unidirectional flux, and to a bypass of CK shuttle by direct ATP transfer. Our results provide experimental evidence for the existence of two pathways of energy transfer, direct ATP transfer, and PCr transfer through the CK shuttle, whose contribution may vary depending on the metabolic status of the heart.

ADP compartmentation analysis reveals coupling between pyruvate kinase and ATPases in heart muscle

Cardiomyocytes have intracellular diffusion restrictions, which spatially compartmentalize ADP and ATP. However, the models that predict diffusion restrictions have used data sets generated in rat heart permeabilized fibers, where diffusion distances may be heterogeneous. This is avoided by using isolated, permeabilized cardiomyocytes. The aim of this work was to analyze the intracellular diffusion of ATP and ADP in rat permeabilized cardiomyocytes. To do this, we measured respiration rate, ATPase rate, and ADP concentration in the surrounding solution. The data were analyzed using mathematical models that reflect different levels of cell compartmentalization. In agreement with previous studies, we found significant diffusion restriction by the mitochondrial outer membrane and confirmed a functional coupling between mitochondria and a fraction of ATPases in the cell. In addition, our experimental data show that considerable activity of endogenous pyruvate kinase (PK) remains in the cardiomyocytes after permeabilization. A fraction of ATPases were inactive without ATP feedback by this endogenous PK. When analyzing the data, we were able to reproduce the measurements only with the mathematical models that include a tight coupling between the fraction of endogenous PK and ATPases. To our knowledge, this is the first time such a strong coupling of PK to ATPases has been demonstrated in permeabilized cardiomyocytes.

Complement regulator CD59 protects against angiotensin II-induced abdominal aortic aneurysms in mice

BACKGROUND

Complement system, an innate immunity, has been well documented to play a critical role in many inflammatory diseases. However, the role of complement in the pathogenesis of abdominal aortic aneurysm, which is considered an immune and inflammatory disease, remains obscure.

METHODS AND RESULTS

Here, we evaluated the pathogenic roles of complement membrane attack complex and CD59, a key regulator that inhibits the membrane attack complex, in the development of abdominal aortic aneurysm. We demonstrated that in the angiotensin II-induced abdominal aortic aneurysm model, deficiency of the membrane attack complex regulator CD59 in ApoE-null mice (mCd59ab(-/-)/ApoE(-/-)) accelerated the disease development, whereas transgenic overexpression of human CD59 (hCD59(ICAM-2+/-)/ApoE(-/-)) in this model attenuated the progression of abdominal aortic aneurysm. The severity of aneurysm among these 3 groups positively correlates with C9 deposition, and/or the activities of MMP2 and MMP9, and/or the levels of phosphorylated c-Jun, c-Fos, IKK-alpha/beta, and p65. Furthermore, we demonstrated that the membrane attack complex directly induced gene expression of matrix metalloproteinase-2 and -9 in vitro, which required activation of the activator protein-1 and nuclear factor-kappaB signaling pathways.

CONCLUSIONS

Together, these results defined the protective role of CD59 and shed light on the important pathogenic role of the membrane attack complex in abdominal aortic aneurysm.

Application of the principles of systems biology and Wiener's cybernetics for analysis of regulation of energy fluxes in muscle cells in vivo

The mechanisms of regulation of respiration and energy fluxes in the cells are analyzed based on the concepts of systems biology, non-equilibrium steady state kinetics and applications of Wiener’s cybernetic principles of feedback regulation. Under physiological conditions cardiac function is governed by the Frank-Starling law and the main metabolic characteristic of cardiac muscle cells is metabolic homeostasis, when both workload and respiration rate can be changed manifold at constant intracellular level of phosphocreatine and ATP in the cells. This is not observed in skeletal muscles. Controversies in theoretical explanations of these observations are analyzed. Experimental studies of permeabilized fibers from human skeletal muscle vastus lateralis and adult rat cardiomyocytes showed that the respiration rate is always an apparent hyperbolic but not a sigmoid function of ADP concentration. It is our conclusion that realistic explanations of regulation of energy fluxes in muscle cells require systemic approaches including application of the feedback theory of Wiener’s cybernetics in combination with detailed experimental research. Such an analysis reveals the importance of limited permeability of mitochondrial outer membrane for ADP due to interactions of mitochondria with cytoskeleton resulting in quasi-linear dependence of respiration rate on amplitude of cyclic changes in cytoplasmic ADP concentrations. The system of compartmentalized creatine kinase (CK) isoenzymes functionally coupled to ANT and ATPases, and mitochondrial-cytoskeletal interactions separate energy fluxes (mass and energy transfer) from signalling (information transfer) within dissipative metabolic structures – intracellular energetic units (ICEU). Due to the non-equilibrium state of CK reactions, intracellular ATP utilization and mitochondrial ATP regeneration are interconnected by the PCr flux from mitochondria. The feedback regulation of respiration occurring via cyclic fluctuations of cytosolic ADP, Pi and Cr/PCr ensures metabolic stability necessary for normal function of cardiac cells.

Matching ATP supply and demand in mammalian heart: in vivo, in vitro, and in silico perspectives

Although the heart rapidly adapts cardiac output to match the body's circulatory demands, the regulatory mechanisms ensuring that sufficient ATP is available to perform the required cardiac work are not completely understood. Two mechanisms have been suggested to serve as key regulators: (1) ADP and Pi concentrations--ATP utilization/hydrolysis in the cytosol increases ADP and Pi fluxes to mitochondria and hence the amount of available substrates for ATP production increases; and (2) Ca2+ concentration--ATP utilization/hydrolysis is coupled to changes in free cytosolic calcium and mitochondrial calcium, the latter controlling Ca2+-dependent activation of mitochondrial enzymes taking part in ATP production. Here we discuss the evolving perspectives of each of the putative regulatory mechanisms and the precise molecular targets (dehydrogenase enzymes, ATP synthase) based on existing experimental and theoretical evidence. The data synthesis can generate novel hypotheses and experimental designs to solve the ongoing enigma of energy supply-demand matching in the heart.

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 Ca(2+) 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 DeltaPsi(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 O(2.-) 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 DeltaPsi(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.

Structure-function relationships in feedback regulation of energy fluxes in vivo in health and disease: mitochondrial interactosome

The aim of this review is to analyze the results of experimental research of mechanisms of regulation of mitochondrial respiration in cardiac and skeletal muscle cells in vivo obtained by using the permeabilized cell technique. Such an analysis in the framework of Molecular Systems Bioenergetics shows that the mechanisms of regulation of energy fluxes depend on the structural organization of the cells and interaction of mitochondria with cytoskeletal elements. Two types of cells of cardiac phenotype with very different structures were analyzed: adult cardiomyocytes and continuously dividing cancerous HL-1 cells. In cardiomyocytes mitochondria are arranged very regularly, and show rapid configuration changes of inner membrane but no fusion or fission, diffusion of ADP and ATP is restricted mostly at the level of mitochondrial outer membrane due to an interaction of heterodimeric tubulin with voltage dependent anion channel, VDAC. VDAC with associated tubulin forms a supercomplex, Mitochondrial Interactosome, with mitochondrial creatine kinase, MtCK, which is structurally and functionally coupled to ATP synthasome. Due to selectively limited permeability of VDAC for adenine nucleotides, mitochondrial respiration rate depends almost linearly upon the changes of cytoplasmic ADP concentration in their physiological range. Functional coupling of MtCK with ATP synthasome amplifies this signal by recycling adenine nucleotides in mitochondria coupled to effective phosphocreatine synthesis. In cancerous HL-1 cells this complex is significantly modified: tubulin is replaced by hexokinase and MtCK is lacking, resulting in direct utilization of mitochondrial ATP for glycolytic lactate production and in this way contributing in the mechanism of the Warburg effect. Systemic analysis of changes in the integrated system of energy metabolism is also helpful for better understanding of pathogenesis of many other diseases.