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

Metabolic Basis of Creatine in Health and Disease: A Bioinformatics-Assisted Review

Creatine (Cr) is a ubiquitous molecule that is synthesized mainly in the liver, kidneys, and pancreas. Most of the Cr pool is found in tissues with high-energy demands. Cr enters target cells through a specific symporter called Na⁺/Cl⁻-dependent Cr transporter (CRT). Once within cells, creatine kinase (CK) catalyzes the reversible transphosphorylation reaction between [Mg²⁺:ATP⁴⁻]²⁻ and Cr to produce phosphocreatine (PCr) and [Mg²⁺:ADP³⁻]⁻. We aimed to perform a comprehensive and bioinformatics-assisted review of the most recent research findings regarding Cr metabolism. Specifically, several public databases, repositories, and bioinformatics tools were utilized for this endeavor. Topics of biological complexity ranging from structural biology to cellular dynamics were addressed herein. In this sense, we sought to address certain pre-specified questions including: (i) What happens when creatine is transported into cells? (ii) How is the CK/PCr system involved in cellular bioenergetics? (iii) How is the CK/PCr system compartmentalized throughout the cell? (iv) What is the role of creatine amongst different tissues? and (v) What is the basis of creatine transport? Under the cellular allostasis paradigm, the CK/PCr system is physiologically essential for life (cell survival, growth, proliferation, differentiation, and migration/motility) by providing an evolutionary advantage for rapid, local, and temporal support of energy- and mechanical-dependent processes. Thus, we suggest the CK/PCr system acts as a dynamic biosensor based on chemo-mechanical energy transduction, which might explain why dysregulation in Cr metabolism contributes to a wide range of diseases besides the mitigating effect that Cr supplementation may have in some of these disease states.

Intracellular Energy-Transfer Networks and High-Resolution Respirometry: A Convenient Approach for Studying Their Function

Compartmentalization of high-energy phosphate carriers between intracellular micro-compartments is a phenomenon that ensures efficient energy use. To connect these sites, creatine kinase (CK) and adenylate kinase (AK) energy-transfer networks, which are functionally coupled to oxidative phosphorylation (OXPHOS), could serve as important regulators of cellular energy fluxes. Here, we introduce how selective permeabilization of cellular outer membrane and high-resolution respirometry can be used to study functional coupling between CK or AK pathways and OXPHOS in different cells and tissues. Using the protocols presented here the ability of creatine or adenosine monophosphate to stimulate OXPHOS through CK and AK reactions, respectively, is easily observable and quantifiable. Additionally, functional coupling between hexokinase and mitochondria can be investigated by monitoring the effect of glucose on respiration. Taken together, high-resolution respirometry in combination with permeabilization is a convenient approach for investigating energy-transfer networks in small quantities of cells and tissues in health and in pathology.

Dominant and sensitive control of oxidative flux by the ATP-ADP carrier in human skeletal muscle mitochondria: Effect of lysine acetylation

The adenine nucleotide translocase (ANT) of the mitochondrial inner membrane exchanges ADP for ATP. Mitochondria were isolated from human vastus lateralis muscle (n = 9). Carboxyatractyloside titration of O₂ consumption rate (J_o) at clamped [ADP] of 21 μM gave ANT abundance of 0.97 ± 0.14 nmol ANT/mg and a flux control coefficient of 82% ± 6%. Flux control fell to 1% ± 1% at saturating (2 mM) [ADP]. The KmADP for J_o was 32.4 ± 1.8 μM. In terms of the free (-3) ADP anion this KmADP was 12.0 ± 0.7 μM. A novel luciferase-based assay for ATP production gave KmADP of 13.1 ± 1.9 μM in the absence of ATP competition. The free anion KmADP in this case was 2.0 ± 0.3 μM. Targeted proteomic analyses showed significant acetylation of ANT Lysine23 and that ANT1 was the most abundant isoform. Acetylation of Lysine23 correlated positively with KmADP, r = 0.74, P = 0.022. The findings underscore the central role played by ANT in the control of oxidative phosphorylation, particularly at the energy phosphate levels associated with low ATP demand. As predicted by molecular dynamic modeling, ANT Lysine23 acetylation decreased the apparent affinity of ADP for ANT binding.

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.

The influence of sex, stroke and distance on the lactate characteristics in high performance swimming

Background

In order to achieve world-class performances, regular performance diagnostics is required as an essential prerequisite for guiding high performance sport. In high performance swimming, the lactate performance diagnostic is an important instrument in testing the sport specific endurance capacity. Although the role of lactate as a signaling molecule, fuel and a gluconeogenic substrate is accepted, lactate parameters are discussed concerning stability, explanatory power and interpretability.

Methods

We calculated the individual anaerobic threshold (IAT) of Bunc using the swimming-specific lactate threshold test by Pansold.

Results

The cross-sectional analysis (ANOVA) of n = 398 high performance swimmers showed significant effects for sex, stroke and distance on the IAT, the percentage of personal best time on the IAT (% of PB on IAT) and maximal lactate values (max. bLA). For the freestyle events the IAT decreased, % of PB on IAT and max. bLA increased from 100 to 400 m significantly in men and women. Women showed significantly higher % of PB on IAT with descriptive lower IAT in 7 of 8 analyzed events. Men showed significantly higher max. bLA in 5 of 8 events. In the second step, the analysis of 1902 data sets of these 398 athletes with a multi-level analysis (MLA) showed also significant effects for sex, swimming distance and stroke. For initial status and development over time, the effect sizes for the variables distance and sex were medium to large, whereas for stroke there were no or small effect sizes.

Discussion

These significant results suggest that lactate tests in swimming specifically have to consider the lactate affecting factors sex and distance under consideration of the time period between measurements. Anthropometrical factors and the physiology of women are possible explanations for the relative better performance for lower lactate concentrations compared to men.

Reverse engineering biomolecular systems using -omic data: challenges, progress and opportunities

Recent advances in high-throughput biotechnologies have led to the rapid growing research interest in reverse engineering of biomolecular systems (REBMS). 'Data-driven' approaches, i.e. data mining, can be used to extract patterns from large volumes of biochemical data at molecular-level resolution while 'design-driven' approaches, i.e. systems modeling, can be used to simulate emergent system properties. Consequently, both data- and design-driven approaches applied to -omic data may lead to novel insights in reverse engineering biological systems that could not be expected before using low-throughput platforms. However, there exist several challenges in this fast growing field of reverse engineering biomolecular systems: (i) to integrate heterogeneous biochemical data for data mining, (ii) to combine top-down and bottom-up approaches for systems modeling and (iii) to validate system models experimentally. In addition to reviewing progress made by the community and opportunities encountered in addressing these challenges, we explore the emerging field of synthetic biology, which is an exciting approach to validate and analyze theoretical system models directly through experimental synthesis, i.e. analysis-by-synthesis. The ultimate goal is to address the present and future challenges in reverse engineering biomolecular systems (REBMS) using integrated workflow of data mining, systems modeling and synthetic biology.

From cycling between coupled reactions to the cross-bridge cycle: mechanical power output as an integral part of energy metabolism

ATP delivery and its usage are achieved by cycling of respective intermediates through interconnected coupled reactions. At steady state, cycling between coupled reactions always occurs at zero resistance of the whole cycle without dissipation of free energy. The cross-bridge cycle can also be described by a system of coupled reactions: one energising reaction, which energises myosin heads by coupled ATP splitting, and one de-energising reaction, which transduces free energy from myosin heads to coupled actin movement. The whole cycle of myosin heads via cross-bridge formation and dissociation proceeds at zero resistance. Dissipation of free energy from coupled reactions occurs whenever the input potential overcomes the counteracting output potential. In addition, dissipation is produced by uncoupling. This is brought about by a load dependent shortening of the cross-bridge stroke to zero, which allows isometric force generation without mechanical power output. The occurrence of maximal efficiency is caused by uncoupling. Under coupled conditions, Hill’s equation (velocity as a function of load) is fulfilled. In addition, force and shortening velocity both depend on [Ca²⁺]. Muscular fatigue is triggered when ATP consumption overcomes ATP delivery. As a result, the substrate of the cycle, [MgATP²⁻], is reduced. This leads to a switch off of cycling and ATP consumption, so that a recovery of [ATP] is possible. In this way a potentially harmful, persistent low energy state of the cell can be avoided.

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.

High efficiency of energy flux controls within mitochondrial interactosome in cardiac intracellular energetic units

The aim of our study was to analyze a distribution of metabolic flux controls of all mitochondrial complexes of ATP-Synthasome and mitochondrial creatine kinase (MtCK) in situ in permeabilized cardiac cells. For this we used their specific inhibitors to measure flux control coefficients (C(vi)(JATP)) in two different systems: A) direct stimulation of respiration by ADP and B) activation of respiration by coupled MtCK reaction in the presence of MgATP and creatine. In isolated mitochondria the C(vi)(JATP) were for system A: Complex I - 0.19, Complex III - 0.06, Complex IV 0.18, adenine nucleotide translocase (ANT) - 0.11, ATP synthase - 0.01, Pi carrier - 0.20, and the sum of C(vi)(JATP) was 0.75. In the presence of 10mM creatine (system B) the C(vi)(JATP) were 0.38 for ANT and 0.80 for MtCK. In the permeabilized cardiomyocytes inhibitors had to be added in much higher final concentration, and the following values of C(vi)(JATP) were determined for condition A and B, respectively: Complex I - 0.20 and 0.64, Complex III - 0.41 and 0.40, Complex IV - 0.40 and 0.49, ANT - 0.20 and 0.92, ATP synthase - 0.065 and 0.38, Pi carrier - 0.06 and 0.06, MtCK 0.95. The sum of C(vi)(JATP) was 1.33 and 3.84, respectively. Thus, C(vi)(JATP) were specifically increased under conditions B only for steps involved in ADP turnover and for Complex I in permeabilized cardiomyocytes within Mitochondrial Interactosome, a supercomplex consisting of MtCK, ATP-Synthasome, voltage dependent anion channel associated with tubulin βII which restricts permeability of the mitochondrial outer membrane.

Studies of the role of tubulin beta II isotype in regulation of mitochondrial respiration in intracellular energetic units in cardiac cells

The aim of this study was to investigate the possible role of tubulin βII, a cytoskeletal protein, in regulation of mitochondrial oxidative phosphorylation and energy fluxes in heart cells. This isotype of tubulin is closely associated with mitochondria and co-expressed with mitochondrial creatine kinase (MtCK). It can be rapidly removed by mild proteolytic treatment of permeabilized cardiomyocytes in the absence of stimulatory effect of cytochrome c, that demonstrating the intactness of the outer mitochondrial membrane. Contrary to isolated mitochondria, in permeabilized cardiomyocytes (in situ mitochondria) the addition of pyruvate kinase (PK) and phosphoenolpyruvate (PEP) in the presence of creatine had no effect on the rate of respiration controlled by activated MtCK, showing limited permeability of voltage-dependent anion channel (VDAC) in mitochondrial outer membrane (MOM) for ADP regenerated by MtCK. Under normal conditions, this effect can be considered as one of the most sensitive tests of the intactness of cardiomyocytes and controlled permeability of MOM for adenine nucleotides. However, proteolytic treatment of permeabilized cardiomyocytes with trypsin, by removing mitochondrial βII tubulin, induces high sensitivity of MtCK-regulated respiration to PK-PEP, significantly changes its kinetics and the affinity to exogenous ADP. MtCK coupled to ATP synthasome and to VDAC controlled by tubulin βII provides functional compartmentation of ATP in mitochondria and energy channeling into cytoplasm via phosphotransfer network. Therefore, direct transfer of mitochondrially produced ATP to sites of its utilization is largely avoided under physiological conditions, but may occur in pathology when mitochondria are damaged. This article is part of a Special Issue entitled ''Local Signaling in Myocytes''.

Adaptive control model reveals systematic feedback and key molecules in metabolic pathway regulation

Robust behavior in metabolic pathways resembles stabilized performance in systems under autonomous control. This suggests we can apply control theory to study existing regulation in these cellular networks. Here, we use model-reference adaptive control (MRAC) to investigate the dynamics of de novo sphingolipid synthesis regulation in a combined theoretical and experimental case study. The effects of serine palmitoyltransferase over-expression on this pathway are studied in vitro using human embryonic kidney cells. We report two key results from comparing numerical simulations with observed data. First, MRAC simulations of pathway dynamics are comparable to simulations from a standard model using mass action kinetics. The root-sum-square (RSS) between data and simulations in both cases differ by less than 5%. Second, MRAC simulations suggest systematic pathway regulation in terms of adaptive feedback from individual molecules. In response to increased metabolite levels available for de novo sphingolipid synthesis, feedback from molecules along the main artery of the pathway is regulated more frequently and with greater amplitude than from other molecules along the branches. These biological insights are consistent with current knowledge while being new that they may guide future research in sphingolipid biology. In summary, we report a novel approach to study regulation in cellular networks by applying control theory in the context of robust metabolic pathways. We do this to uncover potential insight into the dynamics of regulation and the reverse engineering of cellular networks for systems biology. This new modeling approach and the implementation routines designed for this case study may be extended to other systems. Supplementary Material is available at www.liebertonline.com/cmb .

Intracellular Energetic Units regulate metabolism in cardiac cells

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

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

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

Mitochondria-cytoskeleton interaction: distribution of β-tubulins in cardiomyocytes and HL-1 cells

Mitochondria-cytoskeleton interactions were analyzed in adult rat cardiomyocytes and in cancerous non-beating HL-1 cells of cardiac phenotype. We show that in adult cardiomyocytes βII-tubulin is associated with mitochondrial outer membrane (MOM). βI-tubulin demonstrates diffused intracellular distribution, βIII-tubulin is colocalized with Z-lines and βIV-tubulin forms microtubular network. HL-1 cells are characterized by the absence of βII-tubulin, by the presence of bundles of filamentous βIV-tubulin and diffusely distributed βI- and βIII-tubulins. Mitochondrial isoform of creatine kinase (MtCK), highly expressed in cardiomyocytes, is absent in HL-1 cells. Our results show that high apparent K(m) for exogenous ADP in regulation of respiration and high expression of MtCK both correlate with the expression of βII-tubulin. The absence of βII-tubulin isotype in isolated mitochondria and in HL-1 cells results in increased apparent affinity of oxidative phosphorylation for exogenous ADP. This observation is consistent with the assumption that the binding of βII-tubulin to mitochondria limits ADP/ATP diffusion through voltage-dependent anion channel of MOM and thus shifts energy transfer via the phosphocreatine pathway. On the other hand, absence of both βII-tubulin and MtCK in HL-1 cells can be associated with their more glycolysis-dependent energy metabolism which is typical for cancer cells (Warburg effect).

Energetics of glucose metabolism: a phenomenological approach to metabolic network modeling

A new formalism to describe metabolic fluxes as well as membrane transport processes was developed. The new flux equations are comparable to other phenomenological laws. Michaelis-Menten like expressions, as well as flux equations of nonequilibrium thermodynamics, can be regarded as special cases of these new equations. For metabolic network modeling, variable conductances and driving forces are required to enable pathway control and to allow a rapid response to perturbations. When applied to oxidative phosphorylation, results of simulations show that whole oxidative phosphorylation cannot be described as a two-flux-system according to nonequilibrium thermodynamics, although all coupled reactions per se fulfill the equations of this theory. Simulations show that activation of ATP-coupled load reactions plus glucose oxidation is brought about by an increase of only two different conductances: a [Ca²⁺] dependent increase of cytosolic load conductances, and an increase of phosphofructokinase conductance by [AMP], which in turn becomes increased through [ADP] generation by those load reactions. In ventricular myocytes, this feedback mechanism is sufficient to increase cellular power output and O₂ consumption several fold, without any appreciable impairment of energetic parameters. Glucose oxidation proceeds near maximal power output, since transformed input and output conductances are nearly equal, yielding an efficiency of about 0.5. This conductance matching is fulfilled also by glucose oxidation of β-cells. But, as a price for the metabolic mechanism of glucose recognition, β-cells have only a limited capability to increase their power output.