Coupling of cell energetics with membrane metabolic sensing. Integrative signaling through creatine kinase phosphotransfer disrupted by M-CK gene knock-out

KATP channel dependent heart multiome atlas

Plasmalemmal ATP sensitive potassium (K_ATP) channels are recognized metabolic sensors, yet their cellular reach is less well understood. Here, transgenic Kir6.2 null hearts devoid of the K_ATP channel pore underwent multiomics surveillance and systems interrogation versus wildtype counterparts. Despite maintained organ performance, the knockout proteome deviated beyond a discrete loss of constitutive K_ATP channel subunits. Multidimensional nano-flow liquid chromatography tandem mass spectrometry resolved 111 differentially expressed proteins and their expanded network neighborhood, dominated by metabolic process engagement. Independent multimodal chemometric gas and liquid chromatography mass spectrometry unveiled differential expression of over one quarter of measured metabolites discriminating the Kir6.2 deficient heart metabolome. Supervised class analogy ranking and unsupervised enrichment analysis prioritized nicotinamide adenine dinucleotide (NAD⁺), affirmed by extensive overrepresentation of NAD⁺ associated circuitry. The remodeled metabolome and proteome revealed functional convergence and an integrated signature of disease susceptibility. Deciphered cardiac patterns were traceable in the corresponding plasma metabolome, with tissue concordant plasma changes offering surrogate metabolite markers of myocardial latent vulnerability. Thus, Kir6.2 deficit precipitates multiome reorganization, mapping a comprehensive atlas of the K_ATP channel dependent landscape.

Bioinspired facilitation of intrinsically conductive polymers: Mediating intra/extracellular electron transfer and microbial metabolism in denitrification

Intrinsically conductive polymers, polyaniline and polyaniline sulfonate (PASAni) were used to explore their effect on denitrification. Denitrification was accelerated 1.90 times by 2 mM PASAni and the possible mechanisms were mainly attributed to the accelerated electron transfer and the enhanced microbial metabolism activity. Intracellular electron transfer was accelerated by PASAni and the acceleration sites were from NADH to coenzyme Q (CoQ), quinone loop, from Complex II to CoQ and from QH₂ to Cyt. c₁. Extracellular electron transfer was accelerated because PASAni promoted more secretion of redox species and PASAni embedded in extracellular polymeric substance (EPS). Moreover, PASAni itselfprovided more electron transfer pathways as redox species. Microbial metabolism activity was also enhanced by PASAni, which was reflected in the increased nitrate/nitrite reductase activity (236.13/155.43%), electron transfer system activity (112.49%), adenosine triphosphate level (133.41%) and EPS content (189.06%). Besides, the enriched Proteobacteria in PASAni supplement system was also conducive to denitrification. This work provided fundamental information for conductive polymers mediating microbial electron transfer and enhancing contaminants biotransformation.

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.

Role of the C-terminus of SUR in the differential regulation of β-cell and cardiac KATP channels by MgADP and metabolism

Key points

β‐Cell K_ATP channels are partially open in the absence of metabolic substrates, whereas cardiac K_ATP channels are closed.

Using cloned channels heterologously expressed in Xenopus oocytes we measured the effect of MgADP on the MgATP concentration–inhibition curve immediately after patch excision.

MgADP caused a far more striking reduction in ATP inhibition of Kir6.2/SUR1 channels than Kir6.2/SUR2A channels; this effect declined rapidly after patch excision.

Exchanging the final 42 amino acids of SUR was sufficient to switch the Mg‐nucleotide regulation of Kir6.2/SUR1 and Kir6.2/SUR2A channels, and partially switch their sensitivity to metabolic inhibition.

Deletion of the C‐terminal 42 residues of SUR abolished MgADP activation of both Kir6.2/SUR1 and Kir6.2/SUR2A channels.

We conclude that the different metabolic sensitivity of Kir6.2/SUR1 and Kir6.2/SUR2A channels is at least partially due to their different regulation by Mg‐nucleotides, which is determined by the final 42 amino acids.

Abstract

ATP‐sensitive potassium (K_ATP) channels couple the metabolic state of a cell to its electrical activity and play important physiological roles in many tissues. In contrast to β‐cell (Kir6.2/SUR1) channels, which open when extracellular glucose levels fall, cardiac (Kir6.2/SUR2A) channels remain closed. This is due to differences in the SUR subunit rather than cell metabolism. As ATP inhibition and MgADP activation are similar for both types of channels, we investigated channel inhibition by MgATP in the presence of 100 μm MgADP immediately after patch excision [when the channel open probability (P_O) is near maximal]. The results were strikingly different: 100 μm MgADP substantially reduced MgATP inhibition of Kir6.2/SUR1, but had no effect on MgATP inhibition of Kir6.2/SUR2A. Exchanging the final 42 residues of SUR2A with that of SUR1 switched the channel phenotype (and vice versa), and deleting this region abolished Mg‐nucleotide activation. This suggests the C‐terminal 42 residues are important for the ability of MgADP to influence ATP inhibition at Kir6.2. This region was also necessary, but not sufficient, for activation of the K_ATP channel in intact cells by metabolic inhibition (azide). We conclude that the ability of MgADP to impair ATP inhibition at Kir6.2 accounts, in part, for the differential metabolic sensitivities of β‐cell and cardiac K_ATP channels.

β‐Cell K_ATP channels are partially open in the absence of metabolic substrates, whereas cardiac K_ATP channels are closed.

Using cloned channels heterologously expressed in Xenopus oocytes we measured the effect of MgADP on the MgATP concentration–inhibition curve immediately after patch excision.

MgADP caused a far more striking reduction in ATP inhibition of Kir6.2/SUR1 channels than Kir6.2/SUR2A channels; this effect declined rapidly after patch excision.

Exchanging the final 42 amino acids of SUR was sufficient to switch the Mg‐nucleotide regulation of Kir6.2/SUR1 and Kir6.2/SUR2A channels, and partially switch their sensitivity to metabolic inhibition.

Deletion of the C‐terminal 42 residues of SUR abolished MgADP activation of both Kir6.2/SUR1 and Kir6.2/SUR2A channels.

We conclude that the different metabolic sensitivity of Kir6.2/SUR1 and Kir6.2/SUR2A channels is at least partially due to their different regulation by Mg‐nucleotides, which is determined by the final 42 amino acids.

Stress Myocardial Blood Flow Heterogeneity Is a Positron Emission Tomography Biomarker of Ventricular Arrhythmias in Patients With Hypertrophic Cardiomyopathy

Patients with hypertrophic cardiomyopathy (HC) are at increased risk of sudden cardiac death. Abnormalities in myocardial blood flow (MBF) detected by positron emission tomography (PET) are common in HC, but a PET marker that identifies patients at risk of sudden cardiac death is lacking. We hypothesized that disparities in regional myocardial perfusion detected by PET would identify patients with HC at risk of ventricular arrhythmias. To test this hypothesis, we quantified global and regional MBFs by ¹³NH₃-PET at rest and at stress, and developed a heterogeneity index to assess MBF heterogeneity in 133 symptomatic patients with HC. The MBF heterogeneity index was computed by dividing the highest by the lowest regional MBF value, at rest and after vasodilator stress, in each patient. High stress MBF heterogeneity was defined as an index of ≧1.85. Patients with HC were stratified by the presence or the absence of ventricular arrhythmias, defined as sustained ventricular tachycardia (VT) and/or nonsustained VT, during follow-up. We found that global and regional MBFs at rest and stress were similar in patients with HC with or without ventricular arrhythmias. Variability in regional stress MBF was observed in both groups, but the stress MBF heterogeneity index was significantly higher in patients with HC who developed ventricular arrhythmias (1.82 ± 0.77 vs 1.49 ± 0.25, p <0.001). A stress MBF heterogeneity index of ≧1.85 was an independent predictor of both sustained VT (hazard ratio 16.1, 95% confidence interval 3.2 to 80.3) and all-VT (sustained-VT + nonsustained VT: hazard ratio 3.7, 95% confidence interval 1.4 to 9.7). High heterogeneity of stress MBF, reflected by an MBF heterogeneity index of ≥1.85, is a PET biomarker for ventricular arrhythmias in symptomatic patients with HC.

Cardiac metabolic imaging: current imaging modalities and future perspectives

In this review, current imaging techniques and their future perspectives in the field of cardiac metabolic imaging in humans are discussed. This includes a range of noninvasive imaging techniques, allowing a detailed investigation of cardiac metabolism in health and disease. The main imaging modalities discussed are magnetic resonance spectroscopy techniques for determination of metabolite content (triglycerides, glucose, ATP, phosphocreatine, and so on), MRI for myocardial perfusion, and single-photon emission computed tomography and positron emission tomography for quantitation of perfusion and substrate uptake.

Human skeletal muscle fibre contractile properties and proteomic profile: adaptations to 3 weeks of unilateral lower limb suspension and active recovery

KEY POINTS

It is generally assumed that muscle fibres go through atrophy following disuse with a loss of specific force and an increase in unloaded shortening velocity. However, the underlying mechanisms remain to be clarified. Most studies have focused on events taking place during the development of disuse, whereas the subsequent recovery phase, which is equally important, has received little attention. Our findings support the hypotheses that the specific force of muscle fibres decreased following unilateral lower limb suspension (ULLS) and returned to normal after 3 weeks of active recovery as a result of a loss and recovery of myosin and actin content. Furthermore, muscle fibres went through extensive qualitative changes in muscle protein pattern following ULLS, and these were reversed by active recovery. Resistance training was very effective in restoring both muscle mass and qualitative muscle changes, indicating that long-term ULLS did not prevent the positive effect of exercise on human muscle.

ABSTRACT

Following disuse, muscle fibre function goes through adaptations such as a loss of specific force (PO /CSA) and an increase in unloaded shortening velocity, which could be a result of both quantitative changes (i.e. atrophy) and qualitative changes in protein pattern. The underlying mechanisms remain to be clarified. In addition, little is known about the recovery of muscle mass and strength following disuse. In the present study, we report an extensive dataset describing, in detail,the functional and protein content adaptations of skeletal muscle in response to both disuse and re-training. Eight young healthy subjects were subjected to 3 weeks of unilateral lower limb suspension (ULLS), a widely used human model of disuse skeletal muscle atrophy. Needle biopsies samples were taken from the vastus lateralis muscle Pre-ULLS, Post-ULLS and after 3 weeks of recovery during which heavy resistance training was performed. After disuse, cross-sectional area (CSA), PO /CSA and myosin concentration (MC) decreased in both type 1 and 2A skinned muscle fibres. After recovery, CSA and MC returned to levels comparable to those observed before disuse, whereas Po/CSA and unloaded shortening velocity reached a higher level. Myosin heavy chain isoform composition of muscle samples did not differ among the experimental groups. To study the mechanisms underlying such adaptations, a two-dimensional proteomic analysis was performed. ULLS induced a reduction of myofibrillar, metabolic (glycolytic and oxidative) and anti-oxidant defence system protein content. Resistance training was very effective in counteracting ULLS-induced alterations, indicating that long-term ULLS did not prevent the positive effect of exercise on human muscle.

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.

Restrictions in ATP diffusion within sarcomeres can provoke ATP-depleted zones impairing exercise capacity in chronic obstructive pulmonary disease

BACKGROUND

Chronic obstructive pulmonary disease (COPD) is characterized by the inability of patients to sustain a high level of ventilation resulting in perceived exertional discomfort and limited exercise capacity of leg muscles at average intracellular ATP levels sufficient to support contractility.

METHODS

Myosin ATPase activity in biopsy samples from healthy and COPD individuals was implemented as a local nucleotide sensor to determine ATP diffusion coefficients within myofibrils. Ergometric parameters clinically measured during maximal exercise tests in both groups were used to define the rates of myosin ATPase reaction and aerobic ATP re-synthesis. The obtained parameters in combination with AK- and CK-catalyzed reactions were implemented to compute the kinetic and steady-state spatial ATP distributions within control and COPD sarcomeres.

RESULTS

The developed reaction-diffusion model of two-dimensional sarcomeric space identified similar, yet extremely low nucleotide diffusion in normal and COPD myofibrils. The corresponding spatio-temporal ATP distributions, constructed during imposed exercise, predicted in COPD sarcomeres a depletion of ATP in the zones of overlap between actin and myosin filaments along the center axis at average cytosolic ATP levels similar to healthy muscles.

CONCLUSIONS

ATP-depleted zones can induce rigor tension foci impairing muscle contraction and increase a risk for sarcomere damages. Thus, intra-sarcomeric diffusion restrictions at limited aerobic ATP re-synthesis can be an additional risk factor contributing to the muscle contractile deficiency experienced by COPD patients.

GENERAL SIGNIFICANCE

This study demonstrates how restricted substrate mobility within a cellular organelle can provoke an energy imbalance state paradoxically occurring at abounding average metabolic resources.

Decline of Phosphotransfer and Substrate Supply Metabolic Circuits Hinders ATP Cycling in Aging Myocardium

Integration of mitochondria with cytosolic ATP-consuming/ATP-sensing and substrate supply processes is critical for muscle bioenergetics and electrical activity. Whether age-dependent muscle weakness and increased electrical instability depends on perturbations in cellular energetic circuits is unknown. To define energetic remodeling of aged atrial myocardium we tracked dynamics of ATP synthesis-utilization, substrate supply, and phosphotransfer circuits through adenylate kinase (AK), creatine kinase (CK), and glycolytic/glycogenolytic pathways using ¹⁸O stable isotope-based phosphometabolomic technology. Samples of intact atrial myocardium from adult and aged rats were subjected to ¹⁸O-labeling procedure at resting basal state, and analyzed using the ¹⁸O-assisted HPLC-GC/MS technique. Characteristics for aging atria were lower inorganic phosphate Pi[¹⁸O], γ-ATP[¹⁸O], β-ADP[¹⁸O], and creatine phosphate CrP[¹⁸O] ¹⁸O-labeling rates indicating diminished ATP utilization-synthesis and AK and CK phosphotransfer fluxes. Shift in dynamics of glycolytic phosphotransfer was reflected in the diminished G6P[¹⁸O] turnover with relatively constant glycogenolytic flux or G1P[¹⁸O] ¹⁸O-labeling. Labeling of G3P[¹⁸O], an indicator of G3P-shuttle activity and substrate supply to mitochondria, was depressed in aged myocardium. Aged atrial myocardium displayed reduced incorporation of ¹⁸O into second (¹⁸O₂), third (¹⁸O₃), and fourth (¹⁸O₄) positions of Pi[¹⁸O] and a lower Pi[¹⁸O]/γ-ATP[¹⁸ O]-labeling ratio, indicating delayed energetic communication and ATP cycling between mitochondria and cellular ATPases. Adrenergic stress alleviated diminished CK flux, AK catalyzed β-ATP turnover and energetic communication in aging atria. Thus, ¹⁸O-assisted phosphometabolomics uncovered simultaneous phosphotransfer through AK, CK, and glycolytic pathways and G3P substrate shuttle deficits hindering energetic communication and ATP cycling, which may underlie energetic vulnerability of aging atrial myocardium.

Sulfonylureas suppress the stimulatory action of Mg-nucleotides on Kir6.2/SUR1 but not Kir6.2/SUR2A KATP channels: a mechanistic study

Sulfonylureas suppress the stimulatory effect of Mg-nucleotides on recombinant β-cell (Kir6.2/SUR1) but not cardiac (Kir6.2/SUR2A) K_ATP channels.

Sulfonylureas, which stimulate insulin secretion from pancreatic β-cells, are widely used to treat both type 2 diabetes and neonatal diabetes. These drugs mediate their effects by binding to the sulfonylurea receptor subunit (SUR) of the ATP-sensitive K⁺ (K_ATP) channel and inducing channel closure. The mechanism of channel inhibition is unusually complex. First, sulfonylureas act as partial antagonists of channel activity, and second, their effect is modulated by MgADP. We analyzed the molecular basis of the interactions between the sulfonylurea gliclazide and Mg-nucleotides on β-cell and cardiac types of K_ATP channel (Kir6.2/SUR1 and Kir6.2/SUR2A, respectively) heterologously expressed in Xenopus laevis oocytes. The SUR2A-Y1206S mutation was used to confer gliclazide sensitivity on SUR2A. We found that both MgATP and MgADP increased gliclazide inhibition of Kir6.2/SUR1 channels and reduced inhibition of Kir6.2/SUR2A-Y1206S. The latter effect can be attributed to stabilization of the cardiac channel open state by Mg-nucleotides. Using a Kir6.2 mutation that renders the K_ATP channel insensitive to nucleotide inhibition (Kir6.2-G334D), we showed that gliclazide abolishes the stimulatory effects of MgADP and MgATP on β-cell K_ATP channels. Detailed analysis suggests that the drug both reduces nucleotide binding to SUR1 and impairs the efficacy with which nucleotide binding is translated into pore opening. Mutation of one (or both) of the Walker A lysines in the catalytic site of the nucleotide-binding domains of SUR1 may have a similar effect to gliclazide on MgADP binding and transduction, but it does not appear to impair MgATP binding. Our results have implications for the therapeutic use of sulfonylureas.

Sulfonylurea receptors regulate the channel pore in ATP-sensitive potassium channels via an intersubunit salt bridge

ATP-sensitive potassium channels play key roles in many tissues by coupling metabolic status to membrane potential. In contrast with other potassium channels, the pore-forming Kir6 subunits must co-assemble in hetero-octameric complexes with ATP-binding cassette (ABC) family sulfonylurea receptor (SUR) subunits to facilitate cell surface expression. Binding of nucleotides and drugs to SUR regulates channel gating but how these responses are communicated within the complex has remained elusive to date. We have now identified an electrostatic interaction, forming part of a functional interface between the cytoplasmic nucleotide-binding domain-2 of SUR2 subunits and the distal C-terminus of Kir6 polypeptides that determines channel response to nucleotide, potassium channel opener and antagonist. Mutation of participating residues disrupted physical interaction and regulation of expressed channels, properties that were restored in paired charge-swap mutants. Equivalent interactions were identified in Kir6.1- and Kir6.2-containing channels suggesting a conserved mechanism of allosteric regulation.

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.

Sarcolemmal ATP-sensitive potassium channels modulate skeletal muscle function under low-intensity workloads

ATP-sensitive potassium (K_ATP) channels have the unique ability to adjust membrane excitability and functions in accordance with the metabolic status of the cell. Skeletal muscles are primary sites of activity-related energy consumption and have K_ATP channels expressed in very high density. Previously, we demonstrated that transgenic mice with skeletal muscle–specific disruption of K_ATP channel function consume more energy than wild-type littermates. However, how K_ATP channel activation modulates skeletal muscle resting and action potentials under physiological conditions, particularly low-intensity workloads, and how this can be translated to muscle energy expenditure are yet to be determined. Here, we developed a technique that allows evaluation of skeletal muscle excitability in situ, with minimal disruption of the physiological environment. Isometric twitching of the tibialis anterior muscle at 1 Hz was used as a model of low-intensity physical activity in mice with normal and genetically disrupted K_ATP channel function. This workload was sufficient to induce K_ATP channel opening, resulting in membrane hyperpolarization as well as reduction in action potential overshoot and duration. Loss of K_ATP channel function resulted in increased calcium release and aggravated activity-induced heat production. Thus, this study identifies low-intensity workload as a trigger for opening skeletal muscle K_ATP channels and establishes that this coupling is important for regulation of myocyte function and thermogenesis. These mechanisms may provide a foundation for novel strategies to combat metabolic derangements when energy conservation or dissipation is required.

Metabolic rates of ATP transfer through creatine kinase (CK Flux) predict clinical heart failure events and death

Morbidity and mortality from heart failure (HF) are high, and current risk stratification approaches for predicting HF progression are imperfect. Adenosine triphosphate (ATP) is required for normal cardiac contraction, and abnormalities in creatine kinase (CK) energy metabolism, the primary myocardial energy reserve reaction, have been observed in experimental and clinical HF. However, the prognostic value of abnormalities in ATP production rates through CK in human HF has not been investigated. Fifty-eight HF patients with nonischemic cardiomyopathy underwent ³¹P magnetic resonance spectroscopy (MRS) to quantify cardiac high-energy phosphates and the rate of ATP synthesis through CK (CK flux) and were prospectively followed for a median of 4.7 years. Multiple-event analysis (MEA) was performed for HF-related events including all-cause and cardiac death, HF hospitalization, cardiac transplantation, and ventricular-assist device placement. Among baseline demographic, clinical, and metabolic parameters, MEA identified four independent predictors of HF events: New York Heart Association (NYHA) class, left ventricular ejection fraction (LVEF), African-American race, and CK flux. Reduced myocardial CK flux was a significant predictor of HF outcomes, even after correction for NYHA class, LVEF, and race. For each increase in CK flux of 1 μmol g⁻¹ s⁻¹, risk of HF-related composite outcomes decreased by 32 to 39%. These findings suggest that reduced CK flux may be a potential HF treatment target. Newer imaging strategies, including noninvasive ³¹P MRS that detect altered ATP kinetics, could thus complement risk stratification in HF and add value in conditions involving other tissues with high energy demands, including skeletal muscle and brain.

Molecular mechanism of sulphonylurea block of K(ATP) channels carrying mutations that impair ATP inhibition and cause neonatal diabetes

Sulphonylurea drugs are the therapy of choice for treating neonatal diabetes (ND) caused by mutations in the ATP-sensitive K(+) channel (KATP channel). We investigated the interactions between MgATP, MgADP, and the sulphonylurea gliclazide with KATP channels expressed in Xenopus oocytes. In the absence of MgATP, gliclazide block was similar for wild-type channels and those carrying the Kir6.2 ND mutations R210C, G334D, I296L, and V59M. Gliclazide abolished the stimulatory effect of MgATP on all channels. Conversely, high MgATP concentrations reduced the gliclazide concentration, producing a half-maximal block of G334D and R201C channels and suggesting a mutual antagonism between nucleotide and gliclazide binding. The maximal extent of high-affinity gliclazide block of wild-type channels was increased by MgATP, but this effect was smaller for ND channels; channels that were least sensitive to ATP inhibition showed the smallest increase in sulphonylurea block. Consequently, G334D and I296L channels were not fully blocked, even at physiological MgATP concentrations (1 mmol/L). Glibenclamide block was also reduced in β-cells expressing Kir6.2-V59M channels. These data help to explain why patients with some mutations (e.g., G334D, I296L) are insensitive to sulphonylurea therapy, why higher drug concentrations are needed to treat ND than type 2 diabetes, and why patients with severe ND mutations are less prone to drug-induced hypoglycemia.

Creatine kinase adenosine triphosphate and phosphocreatine energy supply in a single kindred of patients with hypertrophic cardiomyopathy

A lethal and extensively characterized familial form of hypertrophic cardiomyopathy (HC) is due to a point mutation (Arg403Gln) in the cardiac β-myosin heavy chain gene. Although this is associated with abnormal energy metabolism and progression to heart failure in an animal model, in vivo cardiac energetics have not been characterized in patients with this mutation. Noninvasive phosphorus saturation transfer magnetic resonance spectroscopy was used to measure the adenosine triphosphate supplied by the creatine kinase (CK) reaction and phosphocreatine, the heart's primary energy reserve, in 9 of 10 patients from a single kindred with HC caused by the Arg403GIn mutation and 17 age-matched healthy controls. Systolic and diastolic function was assessed by echocardiography in all 10 patients with HC. The patients with HC had impairment of diastolic function and mild systolic dysfunction, when assessed using global systolic longitudinal strain. Myocardial phosphocreatine was significantly decreased by 24% in patients (7.1 ± 2.3 μmol/g) compared with the controls (9.4 ± 1.2 μmol/g; p = 0.003). The pseudo-first-order CK rate-constant was 26% lower (0.28 ± 0.15 vs 0.38 ± 0.07 s⁻¹, p = 0.035) and the forward CK flux was 44% lower (2.0 ± 1.4 vs 3.6 ± 0.9 μmol/g/s, p = 0.001) than in the controls. The contractile abnormalities did not correlate with the metabolic indexes. In conclusion, myocardial phosphocreatine and CK-ATP delivery are significantly reduced in patients with HC caused by the Arg403Gln mutation, akin to previous results from mice with the same mutation. A lack of a relation between energetic and contractile abnormalities suggests the former result from the sarcomeric mutation and not a late consequence of mechanical dysfunction.

Unchanged mitochondrial organization and compartmentation of high-energy phosphates in creatine-deficient GAMT-/- mouse hearts

Disruption of the creatine kinase (CK) system in hearts of CK-deficient mice leads to changes in the ultrastructure and regulation of mitochondrial respiration. We expected to see similar changes in creatine-deficient mice, which lack the enzyme guanidinoacetate methyltransferase (GAMT) to produce creatine. The aim of this study was to characterize the changes in cardiomyocyte mitochondrial organization, regulation of respiration, and intracellular compartmentation associated with GAMT deficiency. Three-dimensional mitochondrial organization was assessed by confocal microscopy. On populations of permeabilized cardiomyocytes, we recorded ADP and ATP kinetics of respiration, competition between mitochondria and pyruvate kinase for ADP produced by ATPases, ADP kinetics of endogenous pyruvate kinase, and ATP kinetics of ATPases. These data were analyzed by mathematical models to estimate intracellular compartmentation. Quantitative analysis of morphological and kinetic data as well as derived model fits showed no difference between GAMT-deficient and wild-type mice. We conclude that inactivation of the CK system by GAMT deficiency does not alter mitochondrial organization and intracellular compartmentation in relaxed cardiomyocytes. Thus, our results suggest that the healthy heart is able to preserve cardiac function at a basal level in the absence of CK-facilitated energy transfer without compromising intracellular organization and the regulation of mitochondrial energy homeostasis. This raises questions on the importance of the CK system as a spatial energy buffer in unstressed cardiomyocytes.

AMPK modulates glucose-sensing in insulin-secreting cells by altered phosphotransfer to KATP channels

Glucose-sensing (GS) behaviour in pancreatic β-cells is dependent on ATP-sensitive K(+) channel (KATP) activity, which is controlled by the relative levels of the KATP ligands ATP and ADP, responsible for closing and opening KATP, respectively. However, the mechanism by which β-cells transfer energy status from mitochondria to KATP, and hence to altered electrical excitability and insulin secretion, is presently unclear. Recent work has demonstrated a critical role for AMP-activated protein kinase (AMPK) in GS behaviour of cells. Electrophysiological recordings, coupled with measurements of gene and protein expression were made from rat insulinoma cells to investigate whether AMPK activity regulates this energy transfer process. Using the whole-cell recording configuration with sufficient intracellular ATP to keep KATP closed, raised AMPK activity induced GS electrical behaviour. This effect was prevented by the AMPK inhibitor, compound C and required a phosphotransfer process. Indeed, high levels of intracellular phosphocreatine or the presence of the adenylate kinase (AK) inhibitor AP5A blocked this action of AMPK. Using conditions that maximised AMPK-induced KATP opening, there was a significant increase in AK1, AK2 and UCP2 mRNA expression. Thus we propose that KATP opening in response to lowered glucose concentration requires AMPK activity, perhaps in concert with increased AK and UCP2 to enable mitochondrial-derived ADP signals to be transferred to plasma membrane KATP by phosphotransfer cascades.

Relationship of delayed enhancement by magnetic resonance to myocardial perfusion by positron emission tomography in hypertrophic cardiomyopathy

BACKGROUND

Presence of delayed enhancement (DE) on cardiac magnetic resonance (CMR) is associated with worse clinical outcomes in hypertrophic cardiomyopathy. We investigated the relationship between DE on CMR and myocardial ischemia in hypertrophic cardiomyopathy.

METHODS AND RESULTS

Hypertrophic cardiomyopathy patients (n=47) underwent CMR for assessment of DE and vasodilator stress ammonia positron emission tomography to quantify myocardial blood flow and coronary flow reserve. The summed difference score for regional myocardial perfusion was also assessed. Patients in the DE group (n=35) had greater left ventricular wall thickness (2.09±0.44 versus 1.78±0.34 cm; P=0.03). Stress myocardial blood flow (2.25±0.46 versus 1.78±0.43 mL/min per gram; P=0.01) and coronary flow reserve (2.78±0.32 versus 2.01±0.52; P<0.001) were significantly lower in DE-positive patients. Summed difference score (7.3±6.6 versus 0.9±1.4; P<0.0001) was significantly higher in patients with DE. A coronary flow reserve <2.00 was seen in 18 patients (51%) with DE but in none of the DE-negative patients (P<0.0001). CMR and positron emission tomography showed visually concordant DE and regional myocardial perfusion abnormalities in 31 patients and absence of DE and perfusion defects in 9 patients. Four DE-positive patients demonstrated normal regional myocardial perfusion, and 3 DE-negative patients had (apical) regional myocardial perfusion abnormalities.

CONCLUSIONS

We found a close relationship between DE by CMR and microvascular function in most of the patients studied. However, a small proportion of patients had DE in the absence of perfusion abnormalities, suggesting that microvascular dysfunction and ischemia are not the sole causes of DE in hypertrophic cardiomyopathy patients.

Sensitivity analysis of flux determination in heart by H₂ ¹⁸O -provided labeling using a dynamic Isotopologue model of energy transfer pathways

To characterize intracellular energy transfer in the heart, two organ-level methods have frequently been employed: inversion and saturation transfer, and dynamic labeling. Creatine kinase (CK) fluxes obtained by following oxygen labeling have been considerably smaller than the fluxes determined by saturation transfer. It has been proposed that dynamic labeling determines net flux through CK shuttle, whereas saturation transfer measures total unidirectional flux. However, to our knowledge, no sensitivity analysis of flux determination by oxygen labeling has been performed, limiting our ability to compare flux distributions predicted by different methods. Here we analyze oxygen labeling in a physiological heart phosphotransfer network with active CK and adenylate kinase (AdK) shuttles and establish which fluxes determine the labeling state. A mathematical model consisting of a system of ordinary differential equations was composed describing enrichment in each phosphoryl group and inorganic phosphate. By varying flux distributions in the model and calculating the labeling, we analyzed labeling sensitivity to different fluxes in the heart. We observed that the labeling state is predominantly sensitive to total unidirectional CK and AdK fluxes and not to net fluxes. We conclude that measuring dynamic incorporation of into the high-energy phosphotransfer network in heart does not permit unambiguous determination of energetic fluxes with a higher magnitude than the ATP synthase rate when the bidirectionality of fluxes is taken into account. Our analysis suggests that the flux distributions obtained using dynamic labeling, after removing the net flux assumption, are comparable with those from inversion and saturation transfer.

In heart, the movement of energy metabolites between force-producing myosin, other ATPases, and mitochondria is vital for its function and closely related to heart pathologies. In addition to diffusion, transport of ATP, ADP, Pi, and phosphocreatine occurs along parallel pathways such as the adenylate kinase and creatine kinase shuttles. Two organ-level methods have been developed to study the relative flux through these pathways. However, their results differ. It was recently demonstrated that studies often suffer from the exclusion of compartmentation from their metabolic models. One study overcame this limitation by using compartmental models and statistical methods on multiple experiments. Here, we analyzed the sensitivity of the other method - dynamic labeling of phosphoryl groups and inorganic phosphate. For that, we composed a mathematical model tracking enrichment of the metabolites and evaluated sensitivity of labeling to different flux distribution scenarios. Our study shows that the dynamic method provides a measure of total flux, and not net flux as presumed previously, making the fluxes predicted from both methods consistent. Importantly, conclusions derived on the basis of labeling analysis, particularly those regarding the net flux through the shuttles in control and pathological cases, need to be reevaluated.

Regulation of cardiac ATP-sensitive potassium channel surface expression by calcium/calmodulin-dependent protein kinase II

Cardiac ATP-sensitive potassium (K(ATP)) channels are key sensors and effectors of the metabolic status of cardiomyocytes. Alteration in their expression impacts their effectiveness in maintaining cellular energy homeostasis and resistance to injury. We sought to determine how activation of calcium/calmodulin-dependent protein kinase II (CaMKII), a central regulator of calcium signaling, translates into reduced membrane expression and current capacity of cardiac K(ATP) channels. We used real-time monitoring of K(ATP) channel current density, immunohistochemistry, and biotinylation studies in isolated hearts and cardiomyocytes from wild-type and transgenic mice as well as HEK cells expressing wild-type and mutant K(ATP) channel subunits to track the dynamics of K(ATP) channel surface expression. Results showed that activation of CaMKII triggered dynamin-dependent internalization of K(ATP) channels. This process required phosphorylation of threonine at 180 and 224 and an intact (330)YSKF(333) endocytosis motif of the K(ATP) channel Kir6.2 pore-forming subunit. A molecular model of the μ2 subunit of the endocytosis adaptor protein, AP2, complexed with Kir6.2 predicted that μ2 docks by interaction with (330)YSKF(333) and Thr-180 on one and Thr-224 on the adjacent Kir6.2 subunit. Phosphorylation of Thr-180 and Thr-224 would favor interactions with the corresponding arginine- and lysine-rich loops on μ2. We concluded that calcium-dependent activation of CaMKII results in phosphorylation of Kir6.2, which promotes endocytosis of cardiac K(ATP) channel subunits. This mechanism couples the surface expression of cardiac K(ATP) channels with calcium signaling and reveals new targets to improve cardiac energy efficiency and stress resistance.

Analysis of molecular movement reveals latticelike obstructions to diffusion in heart muscle cells

Intracellular diffusion in muscle cells is known to be restricted. Although characteristics and localization of these restrictions is yet to be elucidated, it has been established that ischemia-reperfusion injury reduces the overall diffusion restriction. Here we apply an extended version of raster image correlation spectroscopy to determine directional anisotropy and coefficients of diffusion in rat cardiomyocytes. Our experimental results indicate that diffusion of a smaller molecule (1127 MW fluorescently labeled ATTO633-ATP) is restricted more than that of a larger one (10,000 MW Alexa647-dextran), when comparing diffusion in cardiomyocytes to that in solution. We attempt to provide a resolution to this counterintuitive result by applying a quantitative stochastic model of diffusion. Modeling results suggest the presence of periodic intracellular barriers situated ∼1 μm apart having very low permeabilities and a small effect of molecular crowding in volumes between the barriers. Such intracellular structuring could restrict diffusion of molecules of energy metabolism, reactive oxygen species, and apoptotic signals, enacting a significant role in normally functioning cardiomyocytes as well as in pathological conditions of the heart.

Compartmentation of membrane processes and nucleotide dynamics in diffusion-restricted cardiac cell microenvironment

Orchestrated excitation-contraction coupling in heart muscle requires adequate spatial arrangement of systems responsible for ion movement and metabolite turnover. Co-localization of regulatory and transporting proteins into macromolecular complexes within an environment of microanatomical cell components raises intracellular diffusion barriers that hamper the mobility of metabolites and signaling molecules. Compared to substrate diffusion in the cytosol, diffusional restrictions underneath the sarcolemma are much larger and could impede ion and nucleotide movement by a factor of 10(3)-10(5). Diffusion barriers thus seclude metabolites within the submembrane space enabling rapid and vectorial effector targeting, yet hinder energy supply from the bulk cytosolic space implicating the necessity for a shunting transfer mechanism. Here, we address principles of membrane protein compartmentation, phosphotransfer enzyme-facilitated interdomain energy transfer, and nucleotide signal dynamics at the subsarcolemma-cytosol interface. This article is part of a Special Issue entitled "Local Signaling in Myocytes".

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.

Dynamic phosphometabolomic profiling of human tissues and transgenic models by 18O-assisted ³¹P NMR and mass spectrometry

Next-generation screening of disease-related metabolomic phenotypes requires monitoring of both metabolite levels and turnover rates. Stable isotope (18)O-assisted (31)P nuclear magnetic resonance (NMR) and mass spectrometry uniquely allows simultaneous measurement of phosphometabolite levels and turnover rates in tissue and blood samples. The (18)O labeling procedure is based on the incorporation of one (18)O into P(i) from [(18)O]H(2)O with each act of ATP hydrolysis and the distribution of (18)O-labeled phosphoryls among phosphate-carrying molecules. This enables simultaneous recording of ATP synthesis and utilization, phosphotransfer fluxes through adenylate kinase, creatine kinase, and glycolytic pathways, as well as mitochondrial substrate shuttle, urea and Krebs cycle activity, glycogen turnover, and intracellular energetic communication. Application of expanded (18)O-labeling procedures has revealed significant differences in the dynamics of G-6-P[(18)O] (glycolysis), G-3-P[(18)O] (substrate shuttle), and G-1-P[(18)O] (glycogenolysis) between human and rat atrial myocardium. In human atria, the turnover of G-3-P[(18)O], which defects are associated with the sudden death syndrome, was significantly higher indicating a greater importance of substrate shuttling to mitochondria. Phosphometabolomic profiling of transgenic hearts deficient in adenylate kinase (AK1-/-), which altered levels and mutations are associated to human diseases, revealed a stress-induced shift in metabolomic profile with increased CrP[(18)O] and decreased G-1-P[(18)O] metabolic dynamics. The metabolomic profile of creatine kinase M-CK/ScCKmit-/--deficient hearts is characterized by a higher G-6-[(18)O]P turnover rate, G-6-P levels, glycolytic capacity, γ/β-phosphoryl of GTP[(18)O] turnover, as well as β-[(18)O]ATP and β-[(18)O]ADP turnover, indicating altered glycolytic, guanine nucleotide, and adenylate kinase metabolic flux. Thus, (18)O-assisted gas chromatography-mass spectrometry and (31)P NMR provide a suitable platform for dynamic phosphometabolomic profiling of the cellular energetic system enabling prediction and diagnosis of metabolic diseases states.

Measuring and evaluating the role of ATP-sensitive K+ channels in cardiac muscle

Since ion channels move electrical charge during their activity, they have traditionally been studied using electrophysiological approaches. This was sometimes combined with mathematical models, for example with the description of the ionic mechanisms underlying the initiation and propagation of action potentials in the squid giant axon by Hodgkin and Huxley. The methods for studying ion channels also have strong roots in protein chemistry (limited proteolysis, the use of antibodies, etc.). The advent of the molecular cloning and the identification of genes coding for specific ion channel subunits in the late 1980s introduced a multitude of new techniques with which to study ion channels and the field has been rapidly expanding ever since (e.g. antibody development against specific peptide sequences, mutagenesis, the use of gene targeting in animal models, determination of their protein structures) and new methods are still in development. This review focuses on techniques commonly employed to examine ion channel function in an electrophysiological laboratory. The focus is on the K(ATP) channel, but many of the techniques described are also used to study other ion channels.

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.

Dystrophin is required for the normal function of the cardio-protective K(ATP) channel in cardiomyocytes

Duchenne and Becker muscular dystrophy patients often develop a cardiomyopathy for which the pathogenesis is still unknown. We have employed the murine animal model of Duchenne muscular dystrophy (mdx), which develops a cardiomyopathy that includes some characteristics of the human disease, to study the molecular basis of this pathology. Here we show that the mdx mouse heart has defects consistent with alteration in compounds that regulate energy homeostasis including a marked decrease in creatine-phosphate (PC). In addition, the mdx heart is more susceptible to anoxia than controls. Since the cardio-protective ATP sensitive potassium channel (K_ATP) complex and PC have been shown to interact we investigated whether deficits in PC levels correlate with other molecular events including K_ATP ion channel complex presence, its functionality and interaction with dystrophin. We found that this channel complex is present in the dystrophic cardiac cell membrane but its ability to sense a drop in the intracellular ATP concentration and consequently open is compromised by the absence of dystrophin. We further demonstrate that the creatine kinase muscle isoform (CKm) is displaced from the plasma membrane of the mdx cardiac cells. Considering that CKm is a determinant of K_ATP channel complex function we hypothesize that dystrophin acts as a scaffolding protein organizing the K_ATP channel complex and the enzymes necessary for its correct functioning. Therefore, the lack of proper functioning of the cardio-protective K_ATP system in the mdx cardiomyocytes may be part of the mechanism contributing to development of cardiac disease in dystrophic patients.

Compartmentation of glycogen metabolism revealed from 13C isotopologue distributions

Background

Stable isotope tracers are used to assess metabolic flux profiles in living cells. The existing methods of measurement average out the isotopic isomer distribution in metabolites throughout the cell, whereas the knowledge of compartmental organization of analyzed pathways is crucial for the evaluation of true fluxes. That is why we accepted a challenge to create a software tool that allows deciphering the compartmentation of metabolites based on the analysis of average isotopic isomer distribution.

Results

The software Isodyn, which simulates the dynamics of isotopic isomer distribution in central metabolic pathways, was supplemented by algorithms facilitating the transition between various analyzed metabolic schemes, and by the tools for model discrimination. It simulated ¹³C isotope distributions in glucose, lactate, glutamate and glycogen, measured by mass spectrometry after incubation of hepatocytes in the presence of only labeled glucose or glucose and lactate together (with label either in glucose or lactate). The simulations assumed either a single intracellular hexose phosphate pool, or also channeling of hexose phosphates resulting in a different isotopic composition of glycogen. Model discrimination test was applied to check the consistency of both models with experimental data. Metabolic flux profiles, evaluated with the accepted model that assumes channeling, revealed the range of changes in metabolic fluxes in liver cells.

Conclusions

The analysis of compartmentation of metabolic networks based on the measured ¹³C distribution was included in Isodyn as a routine procedure. The advantage of this implementation is that, being a part of evaluation of metabolic fluxes, it does not require additional experiments to study metabolic compartmentation. The analysis of experimental data revealed that the distribution of measured ¹³C-labeled glucose metabolites is inconsistent with the idea of perfect mixing of hexose phosphates in cytosol. In contrast, the observed distribution indicates the presence of a separate pool of hexose phosphates that is channeled towards glycogen synthesis.

Rearrangement of energetic and substrate utilization networks compensate for chronic myocardial creatine kinase deficiency

Plasticity of the cellular bioenergetic system is fundamental to every organ function, stress adaptation and disease tolerance. Here, remodelling of phosphotransfer and substrate utilization networks in response to chronic creatine kinase (CK) deficiency, a hallmark of cardiovascular disease, has been revealed in transgenic mouse models lacking either cytosolic M-CK (M-CK(-/-)) or both M-CK and sarcomeric mitochondrial CK (M-CK/ScCKmit(-/-)) isoforms. The dynamic metabolomic signatures of these adaptations have also been defined. Tracking perturbations in metabolic dynamics with (18)O and (13)C isotopes and (31)P NMR and mass spectrometry demonstrate that hearts lacking M-CK have lower phosphocreatine (PCr) turnover but increased glucose-6-phosphate (G-6-P) turnover, glucose utilization and inorganic phosphate compartmentation with normal ATP γ-phosphoryl dynamics. Hearts lacking both M-CK and sarcomeric mitochondrial CK have diminished PCr turnover, total phosphotransfer capacity and intracellular energetic communication but increased dynamics of β-phosphoryls of ADP/ATP, G-6-P and γ-/β-phosphoryls of GTP, indicating redistribution of flux through adenylate kinase (AK), glycolytic and guanine nucleotide phosphotransfer circuits. Higher glycolytic and mitochondrial capacities and increased glucose tolerance contributed to metabolic resilience of M-CK/ScCKmit(-/-) mice. Multivariate analysis revealed unique metabolomic signatures for M-CK(-/-) and M-CK/ScCKmit(-/-) hearts suggesting that rearrangements in phosphotransfer and substrate utilization networks provide compensation for genetic CK deficiency. This new information highlights the significance of integrated CK-, AK-, guanine nucleotide- and glycolytic enzyme-catalysed phosphotransfer networks in supporting the adaptivity and robustness of the cellular energetic system.

Symbolic flux analysis for genome-scale metabolic networks

Background

With the advent of genomic technology, the size of metabolic networks that are subject to analysis is growing. A common task when analyzing metabolic networks is to find all possible steady state regimes. There are several technical issues that have to be addressed when analyzing large metabolic networks including accumulation of numerical errors and presentation of the solution to the researcher. One way to resolve those technical issues is to analyze the network using symbolic methods. The aim of this paper is to develop a routine that symbolically finds the steady state solutions of large metabolic networks.

Results

A symbolic Gauss-Jordan elimination routine was developed for analyzing large metabolic networks. This routine was tested by finding the steady state solutions for a number of curated stoichiometric matrices with the largest having about 4000 reactions. The routine was able to find the solution with a computational time similar to the time used by a numerical singular value decomposition routine. As an advantage of symbolic solution, a set of independent fluxes can be suggested by the researcher leading to the formation of a desired flux basis describing the steady state solution of the network. These independent fluxes can be constrained using experimental data. We demonstrate the application of constraints by calculating a flux distribution for the central metabolic and amino acid biosynthesis pathways of yeast.

Conclusions

We were able to find symbolic solutions for the steady state flux distribution of large metabolic networks. The ability to choose a flux basis was found to be useful in the constraint process and provides a strong argument for using symbolic Gauss-Jordan elimination in place of singular value decomposition.

KATP channel Kir6.2 E23K variant overrepresented in human heart failure is associated with impaired exercise stress response

ATP-sensitive K+ (K(ATP)) channels maintain cardiac homeostasis under stress, as revealed by murine gene knockout models of the KCNJ11-encoded Kir6.2 pore. However, the translational significance of K(ATP) channels in human cardiac physiology remains largely unknown. Here, the frequency of the minor K23 allele of the common functional Kir6.2 E23K polymorphism was found overrepresented in 115 subjects with congestive heart failure compared to 2,031 community-based controls (69 vs. 56%, P < 0.001). Moreover, the KK genotype, present in 18% of heart failure patients, was associated with abnormal cardiopulmonary exercise stress testing. In spite of similar baseline heart rates at rest among genotypic subgroups (EE: 72.2 ± 2.3, EK: 75.0 ± 1.8 and KK:77.1 ± 3.0 bpm), subjects with the KK genotype had a significantly reduced heart rate increase at matched workload (EE: 32.8 ± 2.7%, EK: 28.8 ± 2.1%, KK: 21.7 ± 2.6%, P < 0.05), at 75% of maximum oxygen consumption (EE: 53.9 ± 3.9%, EK: 49.9 ± 3.1%, KK: 36.8 ± 5.3%, P < 0.05), and at peak V(O2) (EE: 82.8 ± 6.0%, EK: 80.5 ± 4.7%, KK: 59.7 ± 8.1%, P < 0.05). Molecular modeling of the tetrameric Kir6.2 pore structure revealed the E23 residue within the functionally relevant intracellular slide helix region. Substitution of the wild-type E residue with an oppositely charged, bulkier K residue would potentially result in a significant structural rearrangement and disrupted interactions with neighboring Kir6.2 subunits, providing a basis for altered high-fidelity K(ATP) channel gating, particularly in the homozygous state. Blunted heart rate response during exercise is a risk factor for mortality in patients with heart failure, establishing the clinical relevance of Kir6.2 E23K as a biomarker for impaired stress performance and underscoring the essential role of K(ATP) channels in human cardiac physiology.

Singular features of trypanosomatids' phosphotransferases involved in cell energy management

Trypanosomatids are responsible for economically important veterinary affections and severe human diseases. In Africa, Trypanosoma brucei causes sleeping sickness or African trypanosomiasis, while in America, Trypanosoma cruzi is the etiological agent of Chagas disease. These parasites have complex life cycles which involve a wide variety of environments with very different compositions, physicochemical properties, and availability of metabolites. As the environment changes there is a need to maintain the nucleoside homeostasis, requiring a quick and regulated response. Most of the enzymes required for energy management are phosphotransferases. These enzymes present a nitrogenous group or a phosphate as acceptors, and the most clear examples are arginine kinase, nucleoside diphosphate kinase, and adenylate kinase. Trypanosoma and Leishmania have the largest number of phosphotransferase isoforms ever found in a single cell; some of them are absent in mammals, suggesting that these enzymes are required in many cellular compartments associated to different biological processes. The presence of such number of phosphotransferases support the hypothesis of the existence of an intracellular enzymatic phosphotransfer network that communicates the spatially separated intracellular ATP consumption and production processes. All these unique features make phosphotransferases a promising start point for rational drug design for the treatment of human trypanosomiasis.

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.

Exercise-induced expression of cardiac ATP-sensitive potassium channels promotes action potential shortening and energy conservation

Physical activity is one of the most important determinants of cardiac function. The ability of the heart to increase delivery of oxygen and metabolic fuels relies on an array of adaptive responses necessary to match bodily demand while avoiding exhaustion of cardiac resources. The ATP-sensitive potassium (K(ATP)) channel has the unique ability to adjust cardiac membrane excitability in accordance with ATP and ADP levels, and up-regulation of its expression that occurs in response to exercise could represent a critical element of this adaption. However, the mechanism by which K(ATP) channel expression changes result in a beneficial effect on cardiac excitability and function remains to be established. Here, we demonstrate that an exercise-induced rise in K(ATP) channel expression enhanced the rate and magnitude of action potential shortening in response to heart rate acceleration. This adaptation in membrane excitability promoted significant reduction in cardiac energy consumption under escalating workloads. Genetic disruption of normal K(ATP) channel pore function abolished the exercise-related changes in action potential duration adjustment and caused increased cardiac energy consumption. Thus, an expression-driven enhancement in the K(ATP) channel-dependent membrane response to alterations in cardiac workload represents a previously unrecognized mechanism for adaptation to physical activity and a potential target for cardioprotection.

K(ATP) channel-dependent metaboproteome decoded: systems approaches to heart failure prediction, diagnosis, and therapy

Systems biology provides an integrative platform by which to account for the biological complexity related to cardiac health and disease. In this way, consequences of ATP-sensitive K(+) (K(ATP)) channel deficiency for heart failure prediction, diagnosis, and therapy were resolved recently at a proteomic level. Under stress-free conditions, knockout of the Kir6.2 K(ATP) channel pore induced metabolic proteome remodelling, revealing overrepresentation of markers of cardiovascular disease. Imposed stress precipitated structural and functional defects in Kir6.2-knockout hearts, decreasing survival and validating prediction of disease susceptibility. In the setting of hypertension, a leading risk for heart failure development, proteomic analysis diagnosed the metabolism-centric impact of K(ATP) channel deficiency in disease. Bioinformatic interrogation of K(ATP) channel-dependent proteome prioritized heart-specific adverse effects, exposing cardiomyopathic traits of aggravated contractility, fibrosis, and ventricular hypertrophy. In dilated cardiomyopathy induced by Kir6.2-knockout pressure overload, proteomic remodelling was exacerbated, underlying a multifaceted molecular pathology that indicates the necessity for a broad-based strategy to achieve repair. Embryonic stem cell intervention in cardiomyopathic K(ATP) channel knockout hearts elicited a distinct proteome signature that forecast amelioration of adverse cardiac outcomes. Functional/structural measurements validated improved contractile performance, reduced ventricular size, and decreased cardiac damage in the treated cohort, while systems assessment unmasked cardiovascular development as a prioritized biological function in stem cell-reconstructed hearts. Thus, proteomic deconvolution of K(ATP) channel-deficient hearts provides definitive evidence for the channel's homeostatic contribution to the cardiac metaboproteome and establishes the utility of systems-oriented approaches to predict disease susceptibility, diagnose consequences of heart failure progression, and monitor therapy outcome.

Activation of the K(ATP) channel by Mg-nucleotide interaction with SUR1

The mechanism of adenosine triphosphate (ATP)-sensitive potassium (K_ATP) channel activation by Mg-nucleotides was studied using a mutation (G334D) in the Kir6.2 subunit of the channel that renders K_ATP channels insensitive to nucleotide inhibition and has no apparent effect on their gating. K_ATP channels carrying this mutation (Kir6.2-G334D/SUR1 channels) were activated by MgATP and MgADP with an EC₅₀ of 112 and 8 µM, respectively. This activation was largely suppressed by mutation of the Walker A lysines in the nucleotide-binding domains of SUR1: the remaining small (∼10%), slowly developing component of MgATP activation was fully inhibited by the lipid kinase inhibitor LY294002. The EC₅₀ for activation of Kir6.2-G334D/SUR1 currents by MgADP was lower than that for MgATP, and the time course of activation was faster. The poorly hydrolyzable analogue MgATPγS also activated Kir6.2-G334D/SUR1. AMPPCP both failed to activate Kir6.2-G334D/SUR1 and to prevent its activation by MgATP. Maximal stimulatory concentrations of MgATP (10 mM) and MgADP (1 mM) exerted identical effects on the single-channel kinetics: they dramatically elevated the open probability (P_O > 0.8), increased the mean open time and the mean burst duration, reduced the frequency and number of interburst closed states, and eliminated the short burst states. By comparing our results with those obtained for wild-type K_ATP channels, we conclude that the MgADP sensitivity of the wild-type K_ATP channel can be described quantitatively by a combination of inhibition at Kir6.2 (measured for wild-type channels in the absence of Mg²⁺) and activation via SUR1 (determined for Kir6.2-G334D/SUR1 channels). However, this is not the case for the effects of MgATP.

Targeted disruption of K(ATP) channels aggravates cardiac toxicity in cocaine abuse

Cocaine is the most frequently used illicit drug among individuals seeking emergency-room care, with fatal outcome most often attributable to the cardiovascular manifestations of drug abuse. While the symptomatic presentations of cocaine toxicity are increasingly understood, the molecular determinants that define outcome remain largely unknown. Here, we report that the susceptibility to cocaine-induced cardiotoxicity is genetically regulated. Targeted deletion of the KCNJ11-encoded Kir6.2 pore-forming subunit of sarcolemmal K(ATP) channels resulted in amplified vulnerability to the toxic effects of chronic cocaine abuse. Under the hyperadrenergic stress, imposed by daily 3-week-long intraperitoneal administration of 30 mg/kg cocaine in Kir6.2-knockout mice, failure to maintain cardiac homeostasis translated into decreased exercise tolerance revealed by poor treadmill stress performance, and dilated hypokinetic left hearts with aggravated cellular hypertrophy and pathognomonic characteristics of chronic cocaine-induced cardiac toxicity. This study therefore reveals a previously unrecognized role of Kir6.2-encoded K(ATP) channels in determining cardiovascular outcome in chronic cocaine abuse, identifying a novel molecular determinant of cocaine cardiotoxicity.

K(ATP) channels process nucleotide signals in muscle thermogenic response

Uniquely gated by intracellular adenine nucleotides, sarcolemmal ATP-sensitive K(+) (K(ATP)) channels have been typically assigned to protective cellular responses under severe energy insults. More recently, K(ATP) channels have been instituted in the continuous control of muscle energy expenditure under non-stressed, physiological states. These advances raised the question of how K(ATP) channels can process trends in cellular energetics within a milieu where each metabolic system is set to buffer nucleotide pools. Unveiling the mechanistic basis of the K(ATP) channel-driven thermogenic response in muscles thus invites the concepts of intracellular compartmentalization of energy and proteins, along with nucleotide signaling over diffusion barriers. Furthermore, it requires gaining insight into the properties of reversibility of intrinsic ATPase activity associated with K(ATP) channel complexes. Notwithstanding the operational paradigm, the homeostatic role of sarcolemmal K(ATP) channels can be now broadened to a wider range of environmental cues affecting metabolic well-being. In this way, under conditions of energy deficit such as ischemic insult or adrenergic stress, the operation of K(ATP) channel complexes would result in protective energy saving, safeguarding muscle performance and integrity. Under energy surplus, downregulation of K(ATP) channel function may find potential implications in conditions of energy imbalance linked to obesity, cold intolerance and associated metabolic disorders.

[Cloning and tissue expression analysis of creatine kinase (M-CK) cDNA from the mandarin fish, Siniperca chuatsi]

The creatine kinase (CK) cDNA from the mandarin fish Siniperca chuatsi was cloned by RT-PCR and rapid amplification of cDNA ends (RACE) methods. The structural characteristics and phylogeny of this gene were analyzed. Sequence analysis revealed a 1586 bp cDNA sequence containing 92 bp 5'-untranslated region, 348 bp 3'-untranslated region and 1146 bp open reading frame (ORF), which encoded 381 amino acids. Conserved sequence blocks of vertebrate CKs and diagnostic boxes for the muscle CK (M-CK) isozyme were identified in S. chuatsi CK. Siniperca chuatsi CK showed a higher similarity with vertebrates M-CK isozyme than other CK isozymes (Brain CK, Mitochondrial CKs) and grouped with M-CK isozyme in CK phylogeny, which strongly supported that S. chuatsi CK belongs to M-CK isozyme type. Semi-quantitative RT-PCR analysis demonstrated that the M-CK transcript expression varied among the different tissues and was detected at a high level in skin, ovary, kidney, stomach, muscle and heart, but lower in eye, brain and liver.

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.

ATP-sensitive K+ channel knockout induces cardiac proteome remodeling predictive of heart disease susceptibility

Forecasting disease susceptibility requires detection of maladaptive signatures prior to onset of overt symptoms. A case-in-point are cardiac ATP-sensitive K+ (K(ATP)) channelopathies, for which the substrate underlying disease vulnerability remains to be identified. Resolving molecular pathobiology, even for single genetic defects, mandates a systems platform to reliably diagnose disease predisposition. High-throughput proteomic analysis was here integrated with network biology to decode consequences of Kir6.2 K(ATP) channel pore deletion. Differential two-dimensional gel electrophoresis reproducibly resolved >800 protein species from hearts of asymptomatic wild-type and Kir6.2-knockout counterparts. K(ATP) channel ablation remodeled the cardiac proteome, significantly altering 71 protein spots, from which 102 unique identities were assigned following hybrid linear ion trap quadrupole-Orbitrap tandem mass spectrometry. Ontological annotation stratified the K(ATP) channel-dependent protein cohort into a predominant bioenergetic module (63 resolved identities), with additional focused sets representing signaling molecules (6), oxidoreductases (8), chaperones (6), and proteins involved in catabolism (6), cytostructure (8), and transcription and translation (5). Protein interaction mapping, in conjunction with expression level changes, localized a K(ATP) channel-associated subproteome within a nonstochastic scale-free network. Global assessment of the K(ATP) channel deficient environment verified the primary impact on metabolic pathways and revealed overrepresentation of markers associated with cardiovascular disease. Experimental imposition of graded stress precipitated exaggerated structural and functional myocardial defects in the Kir6.2-knockout, decreasing survivorship and validating the forecast of disease susceptibility. Proteomic cartography thus provides an integral view of molecular remodeling in the heart induced by K(ATP) channel deletion, establishing a systems approach that predicts outcome at a presymptomatic stage.