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.     

Role of phosphocreatine in energy transport in skeletal muscle of bullfrog studied by 31P-NMR

To evaluate the energy-shuttle hypothesis of the phosphocreatine/creatine kinase system, diffusion rates for ATP, phosphocreatine and flux through the creatine kinase reaction were determined by 31P-NMR in resting bullfrog biceps muscle. The diffusion coefficient of phosphocreatine measured by 31P-pulsed gradient NMR was 1.4-times larger than ATP in the muscle, indicating the advantage of phosphocreatine molecules for the intracellular energy transport. The flux of the creatine kinase reaction measured by 31P-saturation transfer NMR was 3.6 mmol/kg wet wt. per s in the resting muscle. The flux is equal to the turnover rate of ATP, ADP, phosphocreatine and creatine molecules, therefore, the life-times of these substrates and the average distance traversed after the life-times by the diffusing molecules were calculated using the diffusion coefficients obtained by 31P-NMR. The mean square length of one-dimensional diffusion was 22 microns in ATP molecules and the minimum diffusion length was 1.8 microns in ADP molecules. The latter was calculated using free ADP concentration, 30 mumol/kg wet wt., obtained from the equilibrium constant of the creatine kinase reaction and the diffusion coefficient assumed to be the same of ATP in muscle. Similar diffusion lengths of ADP were calculated using the reported values for the flux of the creatine kinase reaction in heart and smooth-muscle. The diffusion lengths of all substrates involved in the creatine kinase reaction were larger than the radii of myofibrils. Therefore, in the muscles with an alternating arrangement of mitochondria and myofibrils, such as heart and certain skeletal muscles, ATP and ADP molecules can move freely between myofibrils and mitochondria without the aid of the creatine kinase reaction; thus, we conclude that the energy-shuttle hypothesis is not obligatory for energy transport between the mitochondria and the myofibrils.     

The creatine kinase reaction: a simple reaction with functional complexity

The classical role of PCr is seen as a reservoir of high-energy phosphates defending cellular ATP levels under anaerobic conditions, high rates of energy transfer or rapid fluctuations in energy requirement. Although the high concentration of PCr in glycolytic fast-twitch fibers supports the role of PCr as a buffer of ATP, the primary importance of the creatine kinase (CK) reaction may in fact be to counteract large increases in ADP, which could otherwise inhibit cellular ATPase-mediated systems. A primary role for CK in the maintenance of ADP homeostasis may explain why, in many conditions, there is an inverse relationship between PCr and muscle contractility but not between ATP and muscle contractility. The high rate of ATP hydrolysis during muscle contraction combined with restricted diffusion of ADP suggests that ADP concentration increases transiently during the contraction phase (ADP spikes) and that these are synchronized with the contraction. The presence of CK, structurally bound in close vicinity to the sites of ATP utilization, will reduce the amplitude and duration of the ADP spikes through PCr-mediated phosphotransfer. When PCr is reduced, the efficiency of CK as an ATP buffer will be reduced and the changes in ADP will become more prominent. The presence of ADP spikes is supported by the finding that other processes known to be activated by ADP (i.e. AMP deamination and glycolysis) are stimulated during exercise but not during anoxia, despite the same low global energy state. Breakdown of PCr is driven by increases in ADP above that depicted by the CK equilibrium and the current method to calculate ADPfree from the CK reaction in a contracting muscle is therefore questionable.     

The creatine kinase system and pleiotropic effects of creatine

The pleiotropic effects of creatine (Cr) are based mostly on the functions of the enzyme creatine kinase (CK) and its high-energy product phosphocreatine (PCr). Multidisciplinary studies have established molecular, cellular, organ and somatic functions of the CK/PCr system, in particular for cells and tissues with high and intermittent energy fluctuations. These studies include tissue-specific expression and subcellular localization of CK isoforms, high-resolution molecular structures and structure–function relationships, transgenic CK abrogation and reverse genetic approaches. Three energy-related physiological principles emerge, namely that the CK/PCr systems functions as (a) an immediately available temporal energy buffer, (b) a spatial energy buffer or intracellular energy transport system (the CK/PCr energy shuttle or circuit) and (c) a metabolic regulator. The CK/PCr energy shuttle connects sites of ATP production (glycolysis and mitochondrial oxidative phosphorylation) with subcellular sites of ATP utilization (ATPases). Thus, diffusion limitations of ADP and ATP are overcome by PCr/Cr shuttling, as most clearly seen in polar cells such as spermatozoa, retina photoreceptor cells and sensory hair bundles of the inner ear. The CK/PCr system relies on the close exchange of substrates and products between CK isoforms and ATP-generating or -consuming processes. Mitochondrial CK in the mitochondrial outer compartment, for example, is tightly coupled to ATP export via adenine nucleotide transporter or carrier (ANT) and thus ATP-synthesis and respiratory chain activity, releasing PCr into the cytosol. This coupling also reduces formation of reactive oxygen species (ROS) and inhibits mitochondrial permeability transition, an early event in apoptosis. Cr itself may also act as a direct and/or indirect anti-oxidant, while PCr can interact with and protect cellular membranes. Collectively, these factors may well explain the beneficial effects of Cr supplementation. The stimulating effects of Cr for muscle and bone growth and maintenance, and especially in neuroprotection, are now recognized and the first clinical studies are underway. Novel socio-economically relevant applications of Cr supplementation are emerging, e.g. for senior people, intensive care units and dialysis patients, who are notoriously Cr-depleted. Also, Cr will likely be beneficial for the healthy development of premature infants, who after separation from the placenta depend on external Cr. Cr supplementation of pregnant and lactating women, as well as of babies and infants are likely to be of benefit for child development. Last but not least, Cr harbours a global ecological potential as an additive for animal feed, replacing meat- and fish meal for animal (poultry and swine) and fish aqua farming. This may help to alleviate human starvation and at the same time prevent over-fishing of oceans.     

Role of creatine phosphokinase in cellular function and metabolism.

This paper summarizes the data concerning the role of the creatine phosphokinase system in muscle cells with main attention to the cardiac muscle. Creatine phosphokinase isoenzymes play a key role in the intracellular energy transport from mitochondria to myofibrils and other sites of energy utilization. Due to the existence of the creatine phosphate pathway for energy transport, intracellular creatine phosphate concentration is apparently an important regulatory factor for muscle contraction which influences the contractile force by determining the rate of regeneration of ATP directly available for myosin ATPase, and at the same time controls the activator calcium entry into the myoplasm across the surface membrane of the cells.     

On the role of arginine kinase in insect flight muscle

In the flight muscle of Locusta migratoria L., arginine kinase activity increased 10-fold when 5th instar larvae and adult animals were compared. During the onset of flight, ATP decreased slightly with the amount of phospho-l-arginine remaining constant. Thus, high arginine kinase activity characterizes the adult muscle, giving rise to the speculation that the phospho-l-arginine/l-arginine kinase system does not act only as a buffer system for high-energy phosphate but also as a shuttle mechanism for high-energy phosphate between mitochondria and myofibrils. Judged from electrophoretic mobility, only one isoenzyme exists that is not bound to subcellular structures. Calculations of the diffusive fluxes of ATP, ADP, phosphate, phospho-l-arginine and l-arginine between the sites of ATP-consumption and production, respectively, can be interpreted in such a way, that the low concentration of ADP_free might limit ATP-turnover during flight. Judging from the high arginine kinase activity, the major acceptor for high-energy phosphate at the mitochondria could be l-arginine, while phospho-l-arginine is transphosphorylated to ATP at the myofibrils, thus presumably serving as an energy shuttle.     

Discrimination of cardiac subcellular creatine kinase fluxes by NMR spectroscopy: a new method of analysis.

A challenge in the understanding of creatine kinase (CK) fluxes reflected by NMR magnetization transfer in the perfused rat heart is the choice of a kinetic model of analysis. The complexity of the energetic pathways, due to the presence of adenosine triphosphate (ATP)-inorganic phosphate (Pi) exchange, of metabolite compartmentation and of subcellular localization of CK isozymes cannot be resolve from the sole information obtained from a single NMR protocol. To analyze multicompartment exchanges, we propose a new strategy based on the simultaneous analysis of four inversion transfer protocols. The time course of ATP and Phosphocreatine (PCr) magnetizations computed from the McConnell equations were adjusted to their experimental value for exchange networks of increasing complexity (up to six metabolite pools). Exchange schemes were selected by the quality of their fit and their consistency with data from other sources: the size of mitochondrial pools and the ATP synthesis flux. The consideration of ATP-Pi exchange and of ATP compartmentation were insufficient to describe the data. The most appropriate exchange scheme in our normoxic heart involved the discrimination of three specific CK activities (cytosolic, mitochondrial, and close to ATPases). At the present level of heart contractility, the energy is transferred from mitochondria to myofibrils mainly by PCr.     

Creatine kinase-deficient hearts exhibit increased susceptibility to ischemia-reperfusion injury and impaired calcium homeostasis

The creatine kinase (CK) system is involved in the rapid transport of high-energy phosphates from the mitochondria to the sites of maximal energy requirements such as myofibrils and sarcolemmal ion pumps. Hearts of mice with a combined knockout of cytosolic M-CK and mitochondrial CK (M/Mito-CK(-/-)) show unchanged basal left ventricular (LV) performance but reduced myocardial high-energy phosphate concentrations. Moreover, skeletal muscle from M/Mito-CK(-/-) mice demonstrates altered Ca2+ homeostasis. Our hypothesis was that in CK-deficient hearts, a cardiac phenotype can be unmasked during acute stress conditions and that susceptibility to ischemia-reperfusion injury is increased because of altered Ca2+ homeostasis. We simultaneously studied LV performance and myocardial Ca2+ metabolism in isolated, perfused hearts of M/Mito-CK(-/-) (n = 6) and wild-type (WT, n = 8) mice during baseline, 20 min of no-flow ischemia, and recovery. Whereas LV performance was not different during baseline conditions, LV contracture during ischemia developed significantly earlier (408 +/- 72 vs. 678 +/- 54 s) and to a greater extent (50 +/- 2 vs. 36 +/- 3 mmHg) in M/Mito-CK(-/-) mice. During reperfusion, recovery of diastolic function was impaired (LV end-diastolic pressure: 22 +/- 3 vs. 10 +/- 2 mmHg), whereas recovery of systolic performance was delayed, in M/Mito-CK(-/-) mice. In parallel, Ca2+ transients were similar during baseline conditions; however, M/Mito-CK(-/-) mice showed a greater increase in diastolic Ca2+ concentration ([Ca2+]) during ischemia (237 +/- 54% vs. 167 +/- 25% of basal [Ca2+]) compared with WT mice. In conclusion, CK-deficient hearts show an increased susceptibility of LV performance and Ca2+ homeostasis to ischemic injury, associated with a blunted postischemic recovery. This demonstrates a key function of an intact CK system for maintenance of Ca2+ homeostasis and LV mechanics under metabolic stress conditions.     

Random walk analysis of restricted metabolite diffusion in skeletal myofibril systems

The purpose of this work was the development of a basal mathematical model for the diffusion of low-molecular metabolites in a skeletal muscle cell. A three-dimension diffusion of low-molecular particles was simulated by a Monte-Carlo method (random walks of diffusing molecules). The model takes into account the following structural elements: (i) a regular lattice of actin and myosin filaments inside a myofibril; (ii) the membranes of sarcoplasmic reticulum and mitochondria surrounding the myofibrils; (iii) a set of myofibrils inside a skeletal muscle cell. We simulated diffusion of particles in the bulk of intracellular water phase and their reflections from the rigid surfaces of intracellular structures. The model allowed to calculate the apparent coefficients of particle diffusion in the axial and radial directions, Dparallel(app) and Dperpendicular(app), respectively. In accordance with experimental data from literature, the coefficient Dparallel(app) was independent of time. The coefficient of radial diffusion Dperpendicular(app) decreased with time to steady state values similar to that determined by the NMR diffusion spectroscopy methods. The interactions of diffusing particles with thin and thick filaments of myofibrils could explain the decrease in the Dperpendicular(app) value by a factor of 20%. The collisions of particles with myofilaments began to reveal themselves as a gradual decrease in the Dperpendicular(app) value at early stages of diffusion (t1/2 approximately equal to 0.05 microsec). The contribution of particle reflections from the membranes of sarcoplasmic reticulum and mitochondria to the retardation of the radial diffusion was about of 20-30%, depending on porosity of a membranous shield around the myofibril. For conventional sizes of a membranous shield (diameter 2 microm), the interactions of particles with the shield caused a decrease in the Dperpendicular(app) value with a half-time t1/2 approximately equal to 0.5 msec. This time is essentially lower by a factor about of 100 than that found in published NMR measurements. When we considered diffusion of particles inside a cell compartment confined to impermeable membranous shield, the reflection of particles from this shield led the drastic decrease in the radial diffusion coefficient (Dperpendicular(app) --> porportional to when t --> porportional to). This pattern of the Dperpendicular(app)(t) time-course might be expected in the NMR measurements on skeletal muscle tissue where a sarcolemma represents an impermeable shield for ATP and PCr molecules.     

Creatine supplementation in health and disease. Effects of chronic creatine ingestion in vivo: Down-regulation of the expression of creatine transporter isoforms in skeletal muscle

Interest in creatine (Cr) as a nutritional supplement and ergogenic aid for athletes has surged over recent years. After cellular uptake, Cr is phosphorylated to phosphocreatine (PCr) by the creatine kinase (CK) reaction using ATP. At subcellular sites with high energy requirements, e.g. at the myofibrillar apparatus during muscle contraction, CK catalyzes the transphosphorylation of PCr to ADP to regenerate ATP, thus preventing a depletion of ATP levels. PCr is thus available as an immediate energy source, serving not only as an energy buffer but also as an energy transport vehicle. Ingestion of creatine increases intramuscular Cr, as well as PCr concentrations, and leads to exercise enhancement, especially in sprint performance. Additional benefits of Cr supplementation have also been noticed for high-intensity long-endurance tasks, e.g. shortening of recovery periods after physical exercise.The present article summarizes recent findings on the influence of Cr supplementation on energy metabolism, and introduces the Cr transporter protein (CreaT), responsible for uptake of Cr into cells, as one of the key-players for the multi-faceted regulation of cellular Cr homeostasis. Furthermore, it is suggested that patients with disturbances in Cr metabolism or with different neuro-muscular diseases may benefit from Cr supplementation as an adjuvant therapy to relieve or delay the onset of symptoms. Although it is still unclear how Cr biosynthesis and transport are regulated in health and disease, so far there are no reports of harmful side effects of Cr loading in humans. However, in this study, we report that chronic Cr supplementation in rats down-regulates in vivo the expression of the CreaT. In addition, we describe the presence of CreaT isoforms most likely generated by alternative splicing.     

Functional characterization of the creatine phosphokinase reactions in heart mitochondria and myofibrils]

The kinetic properties of MM-isozyme of creatine phosphokinase (CPK) bound to heart myofibrils have been determined experimentally. It has been shown that CPK isozymes bound to the heart myofibrils and mitochondria are electrophoretically different, but have very similar kinetic properties. For both isozymes the ATP formation reaction is preferable. However, in heart mitochondria the kinetic properties of CPK are compensated for by a tight functional coupling with ATP-ADP translocase. Due to this coupling the ATP formed in the course of oxidative phosphorylation can be used completely for creatine phosphate production in mitochondria. On the other hand, the kinetic properties of myofibrillar CPK isozyme are such that they provide for the effective utilization of creatine phosphate produced in mitochondria for rephosphorylation of AKP formed in the myofibrils during contraction. It is concluded that in the heart cells energy can be transferred from the mitochondria to the myofibrils by creatine phosphate molecules.     

Assays of the metabolic viability of single giant mitochondria. Experiments with intact and impaled mitochondria

Single giant mitochondria isolated from mice fed cuprizone were assayed for their metabolic viability. Two tests were devised. One test optically detected the accumulation of calcium phosphate within the mitochondria under massive loading conditions (including the presence of succinate and ATP). The accumulation corresponds to a test of energy coupling from either electron transport or the hydrolysis of ATP since it is blocked by either antimycin A or oligomycin. The other assay tested for the production of ATP from ADP and Pi, using myofibrils. Myofibrils prepared from glycerinated rabbit psoas muscle contract only in the presence of ATP and not in the presence of ADP. Myofibrillar contraction is unaffected by the presence of antimycin A or oligomycin. However, myofibrils in the presence of mitochondria that are phosphorylating ADP to ATP do contract. This contraction is blocked by antimycin A and/or oligomycin. Hence, the ATP which causes myofibrillar contraction is produced by oxidative phosphorylation. At low mitochondrial concentration, only the myofibrils in close proximity with mitochondria contract in the presence of ADP. Therefore the assay can be used to test the viability of individual mitochondria. Individual giant mitochondria were found to be viable, using both of these assays. Comparable results were obtained in mitochondria impaled with microelectrodes. The potentials and resistances were unaffected by concomitant calcium phosphate accumulation or oxidative phosphorylation.     

Kinetic, thermodynamic, and developmental consequences of deleting creatine kinase isoenzymes from the heart. Reaction kinetics of the creatine kinase isoenzymes in the intact heart

Creatine kinase (CK) exists as a family of isoenzymes in excitable tissue. We studied isolated perfused hearts from mice lacking genes for either the main muscle isoform of CK (M-CK) or both M-CK and the main mitochondrial isoform (Mt-CK) to determine 1) the biological significance of CK isoenzyme shifts, 2) the necessity of maintaining a high CK reaction rate, and 3) the role of CK isoenzymes in establishing the thermodynamics of ATP hydrolysis. (31)P NMR was used to measure [ATP], [PCr], [P(i)], [ADP], pH, as well as the unidirectional reaction rate of PCr--> [gamma-P]ATP. Developmental changes in the main fetal isoform of CK (BB-CK) were unaffected by loss of other CK isoenzymes. In hearts lacking both M- and Mt-CK, the rate of ATP synthesis from PCr was only 9% of the rate of ATP synthesis from oxidative phosphorylation demonstrating a lack of any high energy phosphate shuttle. We also found that the intrinsic activities of the BB-CK and the MM-CK isoenzymes were equivalent. Finally, combined loss of M- and Mt-CK (but not loss of only M-CK) prevented the amount of free energy released from ATP hydrolysis from increasing when pyruvate was provided as a substrate for oxidative phosphorylation.     

Protection of tubular epithelial cells during renal injury via post-transcriptional control of BMP7

The role of mitochondria in alterations that take place in the muscle cell during healthy aging is a matter of debate during recent years. Most of the studies in bioenergetics have a focus on the model of isolated mitochondria, while changes in the crosstalk between working myofibrils and mitochondria in senescent cardiomyocytes have been less studied. The aim of our research was to investigate the modifications in the highly regulated ATP production and energy transfer systems in heart cells in old rat cardiomyocytes. The results of our work demonstrated alterations in the diffusion restrictions of energy metabolites, manifested by changes in the apparent Michaelis–Menten constant of mitochondria to exogenous ADP. The creatine kinase (CK) phosphotransfer pathway efficiency declines significantly in senescence. The ability of creatine to stimulate OXPHOS as well as to increase the affinity of mitochondria for ADP is falling and the most critical decline is already in the 1-year group (middle-age model in rats). Also, a moderate decrease in the adenylate kinase phosphotransfer system was detected. The importance of glycolysis increases in senescence, while the hexokinase activity does not change during healthy aging. The main result of our study is that the decline in the heart muscle performance is not caused by the changes in the respiratory chain complexes activity but mainly by the decrease in the energy transfer efficiency, especially by the CK pathway.     

Creatine kinase binding and possible role in chemically skinned guinea-pig taenia coli.

The activity and role of creatine kinase (CK) associated with contractile proteins of smooth muscle have been investigated using skinned guinea-pig taenia coli fibers. Total CK activity was 163 +/- 22 IU/g (ww) and agarose electrophoresis showed BB, MB, and MM isoforms (BB-CK being the predominant isoenzyme). After skinning for 1 h with Triton X-100, BB-CK was specifically associated with the myofibrils, representing 22% of the preskinned CK activity. When relaxed fibers were exposed to pCa 9 in the presence of 250 microM ADP, 0 ATP and 12 mM PCr, tension was not significantly different from resting tension, but changing to pCa 4.5 caused the fibers to generate 59.1 +/- 5.2 percent of maximal tension. When a high-tension rigor state was achieved (250 microM ADP, 0 ATP, 0 PCr, and pCa 9), the addition of 12 mM PCr effected significant relaxation. These observations implicate an endogenous form of BB-CK, associated with the myofilaments and capable of producing enough ATP for submaximal tension generation and significant relaxation from rigor conditions. It was also shown that ADP is bound to the myofibrils and available for rephosphorylation by BB-CK. These results suggest co-localization of ATPase, MLCK and CK on the contractile proteins of the taenia coli. This enzymic association may play a role in the compartmentation of adenine nucleotides in smooth muscle.     

31P NMR detection of subcellular creatine kinase fluxes in the perfused rat heart: contractility modifies energy transfer pathways

The subcellular fluxes of exchange of ATP and phosphocreatine (PCr) between mitochondria, cytosol, and ATPases were assessed by (31)P NMR spectroscopy to investigate the pathways of energy transfer in a steady state beating heart. Using a combined analysis of four protocols of inversion magnetization transfer associated with biochemical data, three different creatine kinase (CK) activities were resolved in the rat heart perfused in isovolumic control conditions: (i) a cytosolic CK functioning at equilibrium (forward, F(f) = PCr --> ATP, and reverse flux, F(r) = ATP --> PCr = 3.3 mm.s(-1)), (ii) a CK localized in the vicinity of ATPases (MM-CK bound isoform) favoring ATP synthesis (F(f) = 1.7 x F(r)), and (iii) a mitochondrial CK displaced toward PCr synthesis (F(f) = 0.3 and F(r) = 2.6 mm.s(-1)). This study thus provides the first experimental evidence that the energy is carried from mitochondria to ATPases by PCr (i.e. CK shuttle) in the whole heart. In contrast, a single CK functioning at equilibrium was sufficient to describe the data when ATP synthesis was partly inhibited by cyanide (0.15 mm). In this case, ATP was directly transferred from mitochondria to cytosol suggesting that cardiac activity modified energy transfer pathways. Bioenergetic implications of the localization and activity of enzymes within myocardial cells are discussed.     

Compartmentation of creatine kinases during perinatal development of mammalian heart

Maturation of the cardiac cell is characterized by increasing diversity of isozymic expression of creatine kinases. Expression of the M-CK isozyme always precedes that of mitochondrial isozyme (mi-CK), however the expression of an isoform does not inform about its localization or cellular function. The functional role of isozymes binding to sites of energy utilization and production characteristic of the adult myocardium can be evidenced by the functional coupling of M-CK to myofibrillar ATPase and mito-CK to translocase in Triton X-100 and saponin skinned fibers. Functional activity of M-CK and mito-CK were investigated during perinatal development. Both functional activities appear during late fetal life in species mature at birth like guinea pig, and in the first postnatal weeks in immature species like rat or rabbit. Thus, the functional activity of bound CK isozymes is not associated with birthper se but with the general process of cell maturation. Localization of CK in the cytosol appears optimal for the transfer of glycolytic production of ATP to sites of utilization in an immature heart. During cell maturation, the increasing contribution of oxidative phosphorylation to ATP production, the apparition and binding of mi-CK to mitochondria, the binding of M-CK to myofibrils, turn the cell in a compartmentalized system of energy production. This provides the cellular basis for energy transfer by the PCr-Cr-CK system between sites of ATP production and utilization. Compartmentation of both Ca handling and energy turnover leads to a highly structured cell organization and could be essential for the efficiency of heart function.     

Mathematical model of compartmentalized energy transfer: Its use for analysis and interpretation of 31P-NMR studies of isolated heart of creatine kinase deficient mice

A mathematical model of the compartmentalized energy transfer in cardiac cells is described and used for interpretation of novel experimental data obtained by using phosphorus NMR for determination of the energy fluxes in the isolated hearts of transgenic mice with knocked out creatine kinase isoenzymes. These experiments were designed to study the meaning and importance of compartmentation of creatine kinase isoenzymes in the cells in vivo. The model was constructed to describe quantitatively the processes of energy production, transfer, utilization, and feedback between these processes. It describes the production of ATP in mitochondrial matrix space by ATP synthase, use of this ATP for phosphocreatine production in the mitochondrial creatine kinase reaction coupled to the adenine nucleotide translocation, diffusional exchange of metabolites in the cytoplasmic space, and use of phosphocreatine for resynthesis of ATP in the myoplasmic creatine kinase reaction. It accounts also for the recently discovered phenomenon of restricted diffusion of adenine nucleotides through mitochondrial outer membrane porin pores (VDAC). Practically all parameters of the model were determined experimentally. The analysis of energy fluxes between different cellular compartments shows that in all cellular compartments of working heart cells the creatine kinase reaction is far from equilibrium in the systolic phase of the contraction cycle and approaches equilibrium only in cytoplasm and only in the end-diastolic phase of the contraction cycle.Experimental determination of the relationship between energy fluxes by a ³¹P-NMR saturation transfer method and workload in isolated and perfused heart of transgenic mice deficient in MM isoenzyme of the creatine kinase, MM -/- showed that in the hearts from wild mice, containing all creatine kinase isoenzymes, the energy fluxes determined increased 3-4 times with elevation of the workload. By contrast, in the hearts in which only the mitochondrial creatine kinase was active, the energy fluxes became practically independent of the workload in spite of the preservation of 26% of normal creatine kinase activity. These results cannot be explained on the basis of the conventional near-equilibrium theory of creatine kinase in the cells, which excludes any difference between creatine kinase isoenzymes. However, these apparently paradoxical experimental results are quantitatively described by a mathematical model of the compartmentalized energy transfer based on the steady state kinetics of coupled creatine kinase reactions, compartmentation of creatine kinase isoenzymes in the cells, and the kinetics of ATP production and utilization reactions. The use of this model shows that: (1) in the wild type heart cells a major part of energy is transported out of mitochondria via phosphocreatine, which is used for complete regeneration of ATP locally in the myofibrils - this is the quantitative estimate for PCr pathway; (2) however, in the absence of MM-creatine kinase in the myofibrils in transgenic mice the contraction results in a very rapid rise of ADP in cytoplasmic space, that reverses the mitochondrial creatine kinase reaction in the direction of ATP production. In this way, because of increasing concentrations of cytoplasmic ADP, mitochondrial creatine kinase is switched off functionally due to the absence of its counterpart in PCr pathway, MM-creatine kinase. This may explain why the creatine kinase flux becomes practically independent from the workload in the hearts of transgenic mouse without MM-CK. Thus, the analysis of the results of studies of hearts of creatine kinase-deficient transgenic mice, based on the use of a mathematical model of compartmentalized energy transfer, show that in the PCr pathway of intracellular energy transport two isoenzymes of creatine kinase always function in a coordinated manner out of equilibrium, in the steady state, and disturbances in functioning of one of them inevitably result in the disturbances of the other component of the PCr pathway. In the latter case, energy is transferred from mitochondria to myofibrils by alternative metabolic pathways, probably involving adenylate kinase or other systems.     

Structural and functional adaptations of striated muscles to CK deficiency

In adult mammalian muscle cells, energy consuming processes are mainly localized to the sarcolemma, sarcoplasmic reticulum (SR) and myofibrillar compartments, while energy production occurs within mitochondria or glycolytic complexes. Due to the restricted diffusion of adenine nucleotides near the active sites of ATPases involved in contractile activity and calcium homeostasis, there are multiple local systems that can locally rephosphorylate ADP and provide ATP. The creatine kinase (CK) system, with specific isoenzymes localized within each compartment, efficiently controls local adenylate pools and links energy production and utilization. However, mice lacking one or both of the MM-CK and mi-CK isoforms (CK-/-) are viable and develop almost normal cardiac and skeletal muscle function under the conditions of moderate workload, suggesting adaptations or other mechanisms that may ensure efficient energy transfer. While fixed CK is essentially important, other systems could also be involved as well, such as bound glycolytic enzymes or adenylate kinase. We have shown that, additionally, a direct functional interplay exists between mitochondria and sarcoplasmic reticulum, or between mitochondria and myofilaments in muscle cells, that catalyzes direct energy and signal transfer between organelles. In cardiac cells of CK-/- mice, marked cytoarchitectural modifications were observed, and direct adenine nucleotide channeling between mitochondria and organelles was very effective to rescue SR and myofilament functions. In fast skeletal muscles, increased oxidative capacity also indicates compensatory mechanisms. In mutant mice, mitochondrial capacity increases and a direct energy channeling occurs between mitochondria on one hand and ATP consuming sites on the other. However, these systems appear to be insufficient to fully compensate for the lack of CK at high workload. It can be concluded that local rephosphorylation of ADP is a crucial regulatory point in highly differentiated and organized muscle cells to ensure contractile diversity and efficiency and that the CK system is important to control energy fluxes and energy homeostasis.     

Influence of phosphagen concentration on phosphocreatine breakdown kinetics. Data from human gastrocnemius muscle

At the onset of a square-wave exercise of moderate intensity, in the absence of any detectable lactate production, the hydrolysis of phosphocreatine (PCr) fills the gap between energy requirement and energy yield by oxidative pathways, thus representing a readily available source of energy for the muscle. We verified experimentally the relationships between high-energy phosphates and/or their changes and the time constant of PCr concentration ([PCr]) kinetics in humans (tau(PCr)). High-energy phosphate concentration (by (31)P-NMR spectroscopy) in the calf muscles were measured during three repetitions of the rest-to-work transition of moderate aerobic square-wave exercise on nine healthy volunteers, while resting [PCr] was estimated from the appropriate spectroscopy data. PCr concentration decreased significantly (22 +/- 6%) from rest to steady-state exercise, without differences among the three repetitions. Absolute resting [PCr] and tau(PCr) were consistent with literature values, amounting to 27.5 +/- 2.2 mM and 23.9 +/- 2.9 s, respectively. No significant relationships were detected between individual tau(PCr) and mechanical power, fraction or absolute amount of PCr hydrolyzed, or change in ADP concentration. On the contrary, individual tau(PCr) (s) was linearly related to absolute resting [PCr] (mM), the relationship being described by: tau(PCr) = 0.656 + 0.841.[PCr] (n = 9, R = 0.708, P < 0.05). These data support the view that in humans PCr concentration sets the time course of the oxidative metabolism in skeletal muscle at the start of exercise, being one of the main controllers of oxidative phosphorylation.