Philosophical basis and some historical aspects of systems biology: from Hegel to Noble - applications for bioenergetic research

Multi-level advances in databases related to systems pharmacology in traditional Chinese medicine: a 60-year review

The therapeutic effects of traditional Chinese medicine (TCM) involve intricate interactions among multiple components and targets. Currently, computational approaches play a pivotal role in simulating various pharmacological processes of TCM. The application of network analysis in TCM research has provided an effective means to explain the pharmacological mechanisms underlying the actions of herbs or formulas through the lens of biological network analysis. Along with the advances of network analysis, computational science has coalesced around the core chain of TCM research: formula-herb-component-target-phenotype-ZHENG, facilitating the accumulation and organization of the extensive TCM-related data and the establishment of relevant databases. Nonetheless, recent years have witnessed a tendency toward homogeneity in the development and application of these databases. Advancements in computational technologies, including deep learning and foundation model, have propelled the exploration and modeling of intricate systems into a new phase, potentially heralding a new era. This review aims to delves into the progress made in databases related to six key entities: formula, herb, component, target, phenotype, and ZHENG. Systematically discussions on the commonalities and disparities among various database types were presented. In addition, the review raised the issue of research bottleneck in TCM computational pharmacology and envisions the forthcoming directions of computational research within the realm of TCM.

Molecular Mechanisms of Mitochondrial Quality Control in Ischemic Cardiomyopathy

Ischemic cardiomyopathy (ICM) is a special type of coronary heart disease or an advanced stage of the disease, which is related to the pathological mechanism of primary dilated cardiomyopathy. Ischemic cardiomyopathy mainly occurs in the long-term myocardial ischemia, resulting in diffuse myocardial fibrosis. This in turn affects the cardiac ejection function, resulting in a significant impact on myocardial systolic and diastolic function, resulting in a decrease in the cardiac ejection fraction. The pathogenesis of ICM is closely related to coronary heart disease. Mainly due to coronary atherosclerosis caused by coronary stenosis or vascular occlusion, causing vascular inflammatory lesions and thrombosis. As the disease progresses, it leads to long-term myocardial ischemia and eventually ICM. The pathological mechanism is mainly related to the mechanisms of inflammation, myocardial hypertrophy, fibrosis and vascular remodeling. Mitochondria are organelles with a double-membrane structure, so the composition of the mitochondrial outer compartment is basically similar to that of the cytoplasm. When ischemia-reperfusion induces a large influx of calcium into the cell, the concentration of calcium ions in the mitochondrial outer compartment also increases. The subsequent opening of the membrane permeability transition pore in the inner mitochondrial membrane and the resulting calcium overload induces the homeostasis of cardiomyocytes and activates the mitochondrial pathway of apoptosis. Mitochondrial Quality Control (MQC), as an important mechanism for regulating mitochondrial function in cardiomyocytes, affects the morphological structure/function and lifespan of mitochondria. In this review, we discuss the role of MQC (including mitophagy, mitochondrial dynamics, and mitochondrial biosynthesis) in the pathogenesis of ICM and provide important evidence for targeting MQC for ICM.

A mathematical model of GLUT1 modulation in rods and RPE and its differential impact in cell metabolism

We present a mathematical model of key glucose metabolic pathways in two cells of the human retina: the rods and the retinal pigmented epithelium (RPE). Computational simulations of glucose transporter 1 (GLUT1) inhibition in the model accurately reproduce experimental data from conditional knockout mice and reveal that modification of GLUT1 expression levels of both cells differentially impacts their metabolism. We hypothesize that, under glucose scarcity, the RPE's energy producing pathways are altered in order to preserve its functionality, impacting the photoreceptors' outer segment renewal. On the other hand, when glucose is limited in the rods, aerobic glycolysis is preserved, which maintains the lactate contribution to the RPE.

A Convergent Functional Genomics Analysis to Identify Biological Regulators Mediating Effects of Creatine Supplementation

Creatine (Cr) and phosphocreatine (PCr) are physiologically essential molecules for life, given they serve as rapid and localized support of energy- and mechanical-dependent processes. This evolutionary advantage is based on the action of creatine kinase (CK) isozymes that connect places of ATP synthesis with sites of ATP consumption (the CK/PCr system). Supplementation with creatine monohydrate (CrM) can enhance this system, resulting in well-known ergogenic effects and potential health or therapeutic benefits. In spite of our vast knowledge about these molecules, no integrative analysis of molecular mechanisms under a systems biology approach has been performed to date; thus, we aimed to perform for the first time a convergent functional genomics analysis to identify biological regulators mediating the effects of Cr supplementation in health and disease. A total of 35 differentially expressed genes were analyzed. We identified top-ranked pathways and biological processes mediating the effects of Cr supplementation. The impact of CrM on miRNAs merits more research. We also cautiously suggest two dose-response functional pathways (kinase- and ubiquitin-driven) for the regulation of the Cr uptake. Our functional enrichment analysis, the knowledge-based pathway reconstruction, and the identification of hub nodes provide meaningful information for future studies. This work contributes to a better understanding of the well-reported benefits of Cr in sports and its potential in health and disease conditions, although further clinical research is needed to validate the proposed mechanisms.

Self-Organization and Information Processing: From Basic Enzymatic Activities to Complex Adaptive Cellular Behavior

One of the main aims of current biology is to understand the origin of the molecular organization that underlies the complex dynamic architecture of cellular life. Here, we present an overview of the main sources of biomolecular order and complexity spanning from the most elementary levels of molecular activity to the emergence of cellular systemic behaviors. First, we have addressed the dissipative self-organization, the principal source of molecular order in the cell. Intensive studies over the last four decades have demonstrated that self-organization is central to understand enzyme activity under cellular conditions, functional coordination between enzymatic reactions, the emergence of dissipative metabolic networks (DMN), and molecular rhythms. The second fundamental source of order is molecular information processing. Studies on effective connectivity based on transfer entropy (TE) have made possible the quantification in bits of biomolecular information flows in DMN. This information processing enables efficient self-regulatory control of metabolism. As a consequence of both main sources of order, systemic functional structures emerge in the cell; in fact, quantitative analyses with DMN have revealed that the basic units of life display a global enzymatic structure that seems to be an essential characteristic of the systemic functional metabolism. This global metabolic structure has been verified experimentally in both prokaryotic and eukaryotic cells. Here, we also discuss how the study of systemic DMN, using Artificial Intelligence and advanced tools of Statistic Mechanics, has shown the emergence of Hopfield-like dynamics characterized by exhibiting associative memory. We have recently confirmed this thesis by testing associative conditioning behavior in individual amoeba cells. In these Pavlovian-like experiments, several hundreds of cells could learn new systemic migratory behaviors and remember them over long periods relative to their cell cycle, forgetting them later. Such associative process seems to correspond to an epigenetic memory. The cellular capacity of learning new adaptive systemic behaviors represents a fundamental evolutionary mechanism for cell adaptation.

Perspectives From Systems Biology to Improve Knowledge of Leishmania Drug Resistance

Neglected Tropical Diseases include a broad range of pathogens, hosts, and vectors, which represent evolving complex systems. Leishmaniasis, caused by different Leishmania species and transmitted to humans by sandflies, are among such diseases. Leishmania and other Trypanosomatidae display some peculiar features, which make them a complex system to study. Leishmaniasis chemotherapy is limited due to high toxicity of available drugs, long-term treatment protocols, and occurrence of drug resistant parasite strains. Systems biology studies the interactions and behavior of complex biological processes and may improve knowledge of Leishmania drug resistance. System-level studies to understand Leishmania biology have been challenging mainly because of its unusual molecular features. Networks integrating the biochemical and biological pathways involved in drug resistance have been reported in literature. Antioxidant defense enzymes have been identified as potential drug targets against leishmaniasis. These and other biomarkers might be studied from the perspective of systems biology and systems parasitology opening new frontiers for drug development and treatment of leishmaniasis and other diseases. Our main goals include: 1) Summarize current advances in Leishmania research focused on chemotherapy and drug resistance. 2) Share our viewpoint on the application of systems biology to Leishmania studies. 3) Provide insights and directions for future investigation.

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.

Cholesterol and the Safety Factor for Neuromuscular Transmission

A present review is devoted to the analysis of literature data and results of own research. Skeletal muscle neuromuscular junction is specialized to trigger the striated muscle fiber contraction in response to motor neuron activity. The safety factor at the neuromuscular junction strongly depends on a variety of pre- and postsynaptic factors. The review focuses on the crucial role of membrane cholesterol to maintain a high efficiency of neuromuscular transmission. Cholesterol metabolism in the neuromuscular junction, its role in the synaptic vesicle cycle and neurotransmitter release, endplate electrogenesis, as well as contribution of cholesterol to the synaptogenesis, synaptic integrity, and motor disorders are discussed.

Conceptual Foundations of Systems Biology Explaining Complex Cardiac Diseases

Systems biology is an important concept that connects molecular biology and genomics with computing science, mathematics and engineering. An endeavor is made in this paper to associate basic conceptual ideas of systems biology with clinical medicine. Complex cardiac diseases are clinical phenotypes generated by integration of genetic, molecular and environmental factors. Basic concepts of systems biology like network construction, modular thinking, biological constraints (downward biological direction) and emergence (upward biological direction) could be applied to clinical medicine. Especially, in the field of cardiology, these concepts can be used to explain complex clinical cardiac phenotypes like chronic heart failure and coronary artery disease. Cardiac diseases are biological complex entities which like other biological phenomena can be explained by a systems biology approach. The above powerful biological tools of systems biology can explain robustness growth and stability during disease process from modulation to phenotype. The purpose of the present review paper is to implement systems biology strategy and incorporate some conceptual issues raised by this approach into the clinical field of complex cardiac diseases. Cardiac disease process and progression can be addressed by the holistic realistic approach of systems biology in order to define in better terms earlier diagnosis and more effective therapy.

Specialized Functional Diversity and Interactions of the Na,K-ATPase

Na,K-ATPase is a protein ubiquitously expressed in the plasma membrane of all animal cells and vitally essential for their functions. A specialized functional diversity of the Na,K-ATPase isozymes is provided by molecular heterogeneity, distinct subcellular localizations, and functional interactions with molecular environment. Studies over the last decades clearly demonstrated complex and isoform-specific reciprocal functional interactions between the Na,K-ATPase and neighboring proteins and lipids. These interactions are enabled by a spatially restricted ion homeostasis, direct protein-protein/lipid interactions, and protein kinase signaling pathways. In addition to its "classical" function in ion translocation, the Na,K-ATPase is now considered as one of the most important signaling molecules in neuronal, epithelial, skeletal, cardiac and vascular tissues. Accordingly, the Na,K-ATPase forms specialized sub-cellular multimolecular microdomains which act as receptors to circulating endogenous cardiotonic steroids (CTS) triggering a number of signaling pathways. Changes in these endogenous cardiotonic steroid levels and initiated signaling responses have significant adaptive values for tissues and whole organisms under numerous physiological and pathophysiological conditions. This review discusses recent progress in the studies of functional interactions between the Na,K-ATPase and molecular microenvironment, the Na,K-ATPase-dependent signaling pathways and their significance for diversity of cell function.

Elements of the cellular metabolic structure

A large number of studies have demonstrated the existence of metabolic covalent modifications in different molecular structures, which are able to store biochemical information that is not encoded by DNA. Some of these covalent mark patterns can be transmitted across generations (epigenetic changes). Recently, the emergence of Hopfield-like attractor dynamics has been observed in self-organized enzymatic networks, which have the capacity to store functional catalytic patterns that can be correctly recovered by specific input stimuli. Hopfield-like metabolic dynamics are stable and can be maintained as a long-term biochemical memory. In addition, specific molecular information can be transferred from the functional dynamics of the metabolic networks to the enzymatic activity involved in covalent post-translational modulation, so that determined functional memory can be embedded in multiple stable molecular marks. The metabolic dynamics governed by Hopfield-type attractors (functional processes), as well as the enzymatic covalent modifications of specific molecules (structural dynamic processes) seem to represent the two stages of the dynamical memory of cellular metabolism (metabolic memory). Epigenetic processes appear to be the structural manifestation of this cellular metabolic memory. Here, a new framework for molecular information storage in the cell is presented, which is characterized by two functionally and molecularly interrelated systems: a dynamic, flexible and adaptive system (metabolic memory) and an essentially conservative system (genetic memory). The molecular information of both systems seems to coordinate the physiological development of the whole cell.

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.

Systems approaches in integrative cardiac biology: illustrations from cardiac heterocellular signalling studies

Understanding the complexity of cardiac physiology requires system-level studies of multiple cardiac cell types. Frequently, however, the end result of published research lacks the detail of the collaborative and integrative experimental design process, and the underlying conceptual framework. We review the recent progress in systems modelling and omics analysis of the heterocellular heart environment through complementary forward and inverse approaches, illustrating these conceptual and experimental frameworks with case studies from our own research program. The forward approach begins by collecting curated information from the niche cardiac biology literature, and connecting the dots to form mechanistic network models that generate testable system-level predictions. The inverse approach starts from the vast pool of public omics data in recent cardiac biological research, and applies bioinformatics analysis to produce novel candidates for further investigation. We also discuss the possibility of combining these two approaches into a hybrid framework, together with the benefits and challenges. These interdisciplinary research frameworks illustrate the interplay between computational models, omics analysis, and wet lab experiments, which holds the key to making real progress in improving human cardiac wellbeing.

Systems biology: a biologist's viewpoint

The debate over reductionism and antireductionism in biology is very old. Even the systems approach in biology is more than five decades old. However, mainstream biology, particularly experimental biology, has broadly sidestepped those debates and ideas. Post-genome data explosion and development of high-throughput techniques led to resurfacing of those ideas and debates as a new incarnation called Systems Biology. Though experimental biologists have co-opted systems biology and hailed it as a paradigm shift, it is practiced in different shades and understood with divergent meanings. Biology has certain questions linked with organization of multiple components and processes. Often such questions involve multilevel systems. Here in this essay we argue that systems theory provides required framework and abstractions to explore those questions. We argue that systems biology should follow the logical and mathematical approach of systems theory and transmogrification of systems biology to mere collection of higher dimensional data must be avoided. Therefore, the questions that we ask and the priority of those questions should also change. Systems biology should focus on system-level properties and investigate complexity without shying away from it.

Compartmentalization and metabolic channeling for multienzymatic biosynthesis: practical strategies and modeling approaches

: The construction of efficient enzyme complexes for multienzymatic biosynthesis is of increasing interest in order to achieve maximum yield and to minimize the interference due to shortcomings that are typical for straightforward one-pot multienzyme catalysis. These include product or intermediate feedback inhibition, degeneration, and diffusive losses of reaction intermediates, consumption of co-factors, and others. The main mechanisms in nature to tackle these effects in transient or stable protein associations are the formation of metabolic channeling and microcompartments, processes that are desirable also for multienzymatic biosynthesis in vitro. This chapter provides an overview over two main aspects. First, numerous recent strategies for establishing compartmentalized multienzyme associations and constructed synthetic enzyme complexes are reviewed. Second, the computational methods at hand to investigate and optimize such associations systematically, especially with focus on large multienzyme complexes and metabolic channeling, are discussed. Perspectives on future studies of multienzymatic biosynthesis concerning compartmentalization and metabolic channeling are presented.

Defining systems biology: a brief overview of the term and field

Here we provide a broad overview of the definition of the term "systems biology" as well as pinpoint specific events in biological research and beyond that are consistently cited to have contributed and led to the current science of in silico systems biology. Since there have been many reviews and historical accounts describing the term, it would be impossible to include all single references. However, we do attempt to provide a consensus vision of how the field has evolved and consequently the terminology that followed it. We also highlight the development and general acceptance, and use, of standards for model representations as being crucial to the continued success of the in silico systems biology field.

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 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.

The metabolic core and catalytic switches are fundamental elements in the self-regulation of the systemic metabolic structure of cells

BACKGROUND

Experimental observations and numerical studies with dissipative metabolic networks have shown that cellular enzymatic activity self-organizes spontaneously leading to the emergence of a metabolic core formed by a set of enzymatic reactions which are always active under all environmental conditions, while the rest of catalytic processes are only intermittently active. The reactions of the metabolic core are essential for biomass formation and to assure optimal metabolic performance. The on-off catalytic reactions and the metabolic core are essential elements of a Systemic Metabolic Structure which seems to be a key feature common to all cellular organisms.

METHODOLOGY/PRINCIPAL FINDINGS

In order to investigate the functional importance of the metabolic core we have studied different catalytic patterns of a dissipative metabolic network under different external conditions. The emerging biochemical data have been analysed using information-based dynamic tools, such as Pearson's correlation and Transfer Entropy (which measures effective functionality). Our results show that a functional structure of effective connectivity emerges which is dynamical and characterized by significant variations of bio-molecular information flows.

CONCLUSIONS/SIGNIFICANCE

We have quantified essential aspects of the metabolic core functionality. The always active enzymatic reactions form a hub--with a high degree of effective connectivity--exhibiting a wide range of functional information values being able to act either as a source or as a sink of bio-molecular causal interactions. Likewise, we have found that the metabolic core is an essential part of an emergent functional structure characterized by catalytic modules and metabolic switches which allow critical transitions in enzymatic activity. Both, the metabolic core and the catalytic switches in which also intermittently-active enzymes are involved seem to be fundamental elements in the self-regulation of the Systemic Metabolic Structure.

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''.

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."

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).

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.

Quantitative analysis of cellular metabolic dissipative, self-organized structures

One of the most important goals of the postgenomic era is understanding the metabolic dynamic processes and the functional structures generated by them. Extensive studies during the last three decades have shown that the dissipative self-organization of the functional enzymatic associations, the catalytic reactions produced during the metabolite channeling, the microcompartmentalization of these metabolic processes and the emergence of dissipative networks are the fundamental elements of the dynamical organization of cell metabolism. Here we present an overview of how mathematical models can be used to address the properties of dissipative metabolic structures at different organizational levels, both for individual enzymatic associations and for enzymatic networks. Recent analyses performed with dissipative metabolic networks have shown that unicellular organisms display a singular global enzymatic structure common to all living cellular organisms, which seems to be an intrinsic property of the functional metabolism as a whole. Mathematical models firmly based on experiments and their corresponding computational approaches are needed to fully grasp the molecular mechanisms of metabolic dynamical processes. They are necessary to enable the quantitative and qualitative analysis of the cellular catalytic reactions and also to help comprehend the conditions under which the structural dynamical phenomena and biological rhythms arise. Understanding the molecular mechanisms responsible for the metabolic dissipative structures is crucial for unraveling the dynamics of cellular life.

Mitochondrial respiration and membrane potential are regulated by the allosteric ATP-inhibition of cytochrome c oxidase

This paper describes the problems of measuring the allosteric ATP-inhibition of cytochrome c oxidase (CcO) in isolated mitochondria. Only by using the ATP-regenerating system phosphoenolpyruvate and pyruvate kinase full ATP-inhibition of CcO could be demonstrated by kinetic measurements. The mechanism was proposed to keep the mitochondrial membrane potential (DeltaPsi(m)) in living cells and tissues at low values (100-140 mV), when the matrix ATP/ADP ratios are high. In contrast, high DeltaPsi(m) values (180-220 mV) are generally measured in isolated mitochondria. By using a tetraphenyl phosphonium electrode we observed in isolated rat liver mitochondria with glutamate plus malate as substrates a reversible decrease of DeltaPsi(m) from 233 to 123 mV after addition of phosphoenolpyruvate and pyruvate kinase. The decrease of DeltaPsi(m) is explained by reversal of the gluconeogenetic enzymes pyruvate carboxylase and phosphoenolpyruvate carboxykinase yielding ATP and GTP, thus increasing the matrix ATP/ADP ratio. With rat heart mitochondria, which lack these enzymes, no decrease of DeltaPsi(m) was found. From the data we conclude that high matrix ATP/ADP ratios keep DeltaPsi(m) at low values by the allosteric ATP-inhibition of CcO, thus preventing the generation of reactive oxygen species which could generate degenerative diseases. It is proposed that respiration in living eukaryotic organisms is normally controlled by the DeltaPsi(m)-independent "allosteric ATP-inhibition of CcO." Only when the allosteric ATP-inhibition is switched off under stress, respiration is regulated by "respiratory control," based on DeltaPsi(m) according to the Mitchell Theory.

Biophysics and systems biology

Biophysics at the systems level, as distinct from molecular biophysics, acquired its most famous paradigm in the work of Hodgkin and Huxley, who integrated their equations for the nerve impulse in 1952. Their approach has since been extended to other organs of the body, notably including the heart. The modern field of computational biology has expanded rapidly during the first decade of the twenty-first century and, through its contribution to what is now called systems biology, it is set to revise many of the fundamental principles of biology, including the relations between genotypes and phenotypes. Evolutionary theory, in particular, will require re-assessment. To succeed in this, computational and systems biology will need to develop the theoretical framework required to deal with multilevel interactions. While computational power is necessary, and is forthcoming, it is not sufficient. We will also require mathematical insight, perhaps of a nature we have not yet identified. This article is therefore also a challenge to mathematicians to develop such insights.

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.

Global self-regulation of the cellular metabolic structure

BACKGROUND

Different studies have shown that cellular enzymatic activities are able to self-organize spontaneously, forming a metabolic core of reactive processes that remain active under different growth conditions while the rest of the molecular catalytic reactions exhibit structural plasticity. This global cellular metabolic structure appears to be an intrinsic characteristic common to all cellular organisms. Recent work performed with dissipative metabolic networks has shown that the fundamental element for the spontaneous emergence of this global self-organized enzymatic structure could be the number of catalytic elements in the metabolic networks.

METHODOLOGY/PRINCIPAL FINDINGS

In order to investigate the factors that may affect the catalytic dynamics under a global metabolic structure characterized by the presence of metabolic cores we have studied different transitions in catalytic patterns belonging to a dissipative metabolic network. The data were analyzed using non-linear dynamics tools: power spectra, reconstructed attractors, long-term correlations, maximum Lyapunov exponent and Approximate Entropy; and we have found the emergence of self-regulation phenomena during the transitions in the metabolic activities.

CONCLUSIONS/SIGNIFICANCE

The analysis has also shown that the chaotic numerical series analyzed correspond to the fractional Brownian motion and they exhibit long-term correlations and low Approximate Entropy indicating a high level of predictability and information during the self-regulation of the metabolic transitions. The results illustrate some aspects of the mechanisms behind the emergence of the metabolic self-regulation processes, which may constitute an important property of the global structure of the cellular metabolism.

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

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

Origins of systems biology in William Harvey's masterpiece on the movement of the heart and the blood in animals

In this article we continue our exploration of the historical roots of systems biology by considering the work of William Harvey. Central arguments in his work on the movement of the heart and the circulation of the blood can be shown to presage the concepts and methods of integrative systems biology. These include: (a) the analysis of the level of biological organization at which a function (e.g. cardiac rhythm) can be said to occur; (b) the use of quantitative mathematical modelling to generate testable hypotheses and deduce a fundamental physiological principle (the circulation of the blood) and (c) the iterative submission of his predictions to an experimental test. This article is the result of a tri-lingual study: as Harvey's masterpiece was published in Latin in 1628, we have checked the original edition and compared it with and between the English and French translations, some of which are given as notes to inform the reader of differences in interpretation.