Marko Vendelin

Cardiac system bioenergetics

The fundamental principle of cardiac behaviour is described by the Frank‐Starling law relating force of contraction during systole with end‐diastolic volume. While both work and respiration rates increase linearly with imposed load, the basis of mechano‐energetic coupling in heart muscle has remained a long‐standing enigma. Here, we highlight advances made in understanding of complex cellular and molecular mechanisms that orchestrate coupling of mitochondrial oxidative phosphorylation with ATP utilization for muscle contraction. Cardiac system bioenergetics critically depends on an interrelated metabolic infrastructure regulating mitochondrial respiration and energy fluxes throughout cellular compartments. The data reviewed indicate the significance of two interrelated systems regulating mitochondrial respiration and energy fluxes in cells: (1) the creatine kinase, adenylate kinase and glycolytic pathways that communicate flux changes generated by cellular ATPases within structurally organized enzymatic modules and networks; and (2) a secondary system based on mitochondrial participation in cellular calcium cycle, which adjusts substrate oxidation and energy‐transducing processes to meet increasing cellular energy demands. By conveying energetic signals to metabolic sensors, coupled phosphotransfer reactions provide a high‐fidelity regulation of the excitation–contraction cycle. Such integration of energetics with calcium signalling systems provides the basis for ‘metabolic pacing’, synchronizing the cellular electrical and mechanical activities with energy supply processes.

Heterogeneity of ADP diffusion and regulation of respiration in cardiac cells

Heterogeneity of ADP diffusion and regulation of respiration were studied in permeabilized cardiomyocytes and cardiac fibers in situ and in silico. Regular arrangement of mitochondria in cells was altered by short-time treatment with trypsin and visualized by confocal microscopy. Manipulation of matrix volumes by changing K(+) and sucrose concentrations did not affect the affinity for ADP either in isolated heart mitochondria or in skinned fibers. Pyruvate kinase (PK)-phosphoenolpyruvate (PEP) were used to trap ADP generated in Ca,MgATPase reactions. Inhibition of respiration by PK-PEP increased 2-3 times after disorganization of regular mitochondrial arrangement in cells. ADP produced locally in the mitochondrial creatine kinase reaction was not accessible to PK-PEP in intact permeabilized fibers, but some part of it was released from mitochondria after short proteolysis due to increased permeability of outer mitochondrial membrane. In in silico studies we show by mathematical modeling that these results can be explained by heterogeneity of ADP diffusion due to its restrictions at the outer mitochondrial membrane and in close areas, which is changed after proteolysis. Localized restrictions and heterogeneity of ADP diffusion demonstrate the importance of mitochondrial functional complexes with sarcoplasmic reticulum and myofibrillar structures and creatine kinase in regulation of oxidative phosphorylation.

Intracellular diffusion of adenosine phosphates is locally restricted in cardiac muscle

Recent studies have revealed the structural and functional interactions between mitochondria, myofibrils and sarcoplasmic reticulum in cardiac cells. Direct channeling of adenosine phosphates between organelles identified in the experiments indicates that diffusion of adenosine phosphates is limited in cardiac cells due to very specific intracellular structural organization. However, the mode of diffusion restrictions and nature of the intracellular structures in creating the diffusion barriers is still unclear, and, therefore, a subject of active research. The aim of this work is to analyze the possible role of two principally different modes of restriction distribution for adenosine phosphates (a) the uniform diffusion restriction and (b) the localized diffusion limitation in the vicinity of mitochondria, by fitting the experimental data with the mathematical model. The reaction-diffusion model of compartmentalized energy transfer was used to analyze the data obtained from the experiments with the skinned muscle fibers, which described the following processes: mitochondrial respiration rate dependency on exogenous ADP and ATP concentrations; inhibition of endogenous ADP-stimulated respiration by pyruvate kinase (PK) and phosphoenolpyruvate (PEP) system; kinetics of oxygen consumption stabilization after addition of 2 mM MgATP or MgADP; ATPase activity with inhibited mitochondrial respiration; and buildup of MgADP concentration in the medium after addition of MgATP. The analysis revealed that only the second mechanism considered – localization of diffusion restrictions – is able to account for the experimental data. In the case of uniform diffusion restrictions, the model solution was in agreement only with two measurements: the respiration rate as a function of ADP or ATP concentrations and inhibition of respiration by PK + PEP. It was concluded that intracellular diffusion restrictions for adenosine phosphates are not distributed uniformly, but rather are localized in certain compartments of the cardiac cells.

Anisotropic diffusion of fluorescently labeled ATP in rat cardiomyocytes determined by raster image correlation spectroscopy

A series of experimental data points to the existence of profound diffusion restrictions of ADP/ATP in rat cardiomyocytes. This assumption is required to explain the measurements of kinetics of respiration, sarcoplasmic reticulum loading with calcium, and kinetics of ATP-sensitive potassium channels. To be able to analyze and estimate the role of intracellular diffusion restrictions on bioenergetics, the intracellular diffusion coefficients of metabolites have to be determined. The aim of this work was to develop a practical method for determining diffusion coefficients in anisotropic medium and to estimate the overall diffusion coefficients of fluorescently labeled ATP in rat cardiomyocytes. For that, we have extended raster image correlation spectroscopy (RICS) protocols to be able to discriminate the anisotropy in the diffusion coefficient tensor. Using this extended protocol, we estimated diffusion coefficients of ATP labeled with the fluorescent conjugate Alexa Fluor 647 (Alexa-ATP). In the analysis, we assumed that the diffusion tensor can be described by two values: diffusion coefficient along the myofibril and that across it. The average diffusion coefficients found for Alexa-ATP were as follows: 83 ± 14 μm2/s in the longitudinal and 52 ± 16 μm2/s in the transverse directions (n = 8, mean ± SD). Those values are ∼2 (longitudinal) and ∼3.5 (transverse) times smaller than the diffusion coefficient value estimated for the surrounding solution. Such uneven reduction of average diffusion coefficient leads to anisotropic diffusion in rat cardiomyocytes. Although the source for such anisotropy is uncertain, we speculate that it may be induced by the ordered pattern of intracellular structures in rat cardiomyocytes.

The creatine kinase phosphotransfer network: thermodynamic and kinetic considerations, the impact of the mitochondrial outer membrane and modelling approaches.

In this review, we summarize the main structural and functional data on the role of the phosphocreatine (PCr)--creatine kinase (CK) pathway for compartmentalized energy transfer in cardiac cells. Mitochondrial creatine kinase, MtCK, fixed by cardiolipin molecules in the vicinity of the adenine nucleotide translocator, is a key enzyme in this pathway. Direct transfer of ATP and ADP between these proteins has been revealed both in experimental studies on the kinetics of the regulation of mitochondrial respiration and by mathematical modelling as a main mechanism of functional coupling of PCr production to oxidative phosphorylation. In cells in vivo or in permeabilized cells in situ, this coupling is reinforced by limited permeability of the outer membrane of the mitochondria for adenine nucleotides due to the contacts with cytoskeletal proteins. Due to these mechanisms, at least 80% of total energy is exported from mitochondria by PCr molecules. Mathematical modelling of intracellular diffusion and energy transfer shows that the main function of the PCr-CK pathway is to connect different pools (compartments) of ATP and, by this way, to overcome the local restrictions and diffusion limitation of adenine nucleotides due to the high degree of structural organization of cardiac cells.

Application of regularized Richardson–Lucy algorithm for deconvolution of confocal microscopy images

Although confocal microscopes have considerably smaller contribution of out‐of‐focus light than widefield microscopes, the confocal images can still be enhanced mathematically if the optical and data acquisition effects are accounted for. For that, several deconvolution algorithms have been proposed. As a practical solution, maximum‐likelihood algorithms with regularization have been used. However, the choice of regularization parameters is often unknown although it has considerable effect on the result of deconvolution process. The aims of this work were: to find good estimates of deconvolution parameters; and to develop an open source software package that would allow testing different deconvolution algorithms and that would be easy to use in practice. Here, Richardson–Lucy algorithm has been implemented together with the total variation regularization in an open source software package IOCBio Microscope. The influence of total variation regularization on deconvolution process is determined by one parameter. We derived a formula to estimate this regularization parameter automatically from the images as the algorithm progresses. To assess the effectiveness of this algorithm, synthetic images were composed on the basis of confocal images of rat cardiomyocytes. From the analysis of deconvolved results, we have determined under which conditions our estimation of total variation regularization parameter gives good results. The estimated total variation regularization parameter can be monitored during deconvolution process and used as a stopping criterion. An inverse relation between the optimal regularization parameter and the peak signal‐to‐noise ratio of an image is shown. Finally, we demonstrate the use of the developed software by deconvolving images of rat cardiomyocytes with stained mitochondria and sarcolemma obtained by confocal and widefield microscopes.

Cardiac mechanoenergetics replicated by cross-bridge model

AbstractThe aim of this work is to reproduce the experimentally measured linear dependence of cardiac muscle oxygen consumption on stress–strain area using a model, composed of a three-state Huxley-type model for cross-bridge interaction and a phenomenological model of Ca2+-induced activation. By selecting particular cross-bridge cycling rate constants and modifying the cross-bridge activation model, we replicated the linear dependence between oxygen consumption and stress–strain area together with other important mechanical properties of cardiac muscle such as developed stress dependence on the sarcomere length and force-velocity relationship. The model predicts that (1) the amount of the “passenger” cross bridges, i.e., cross bridges that detach without hydrolyzing ATP molecule, is relatively small and (2) ATP consumption rate profile within a beat and the amount of the passenger cross bridges depend on the contraction protocol.

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.

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

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

ADP Compartmentation Analysis Reveals Coupling between Pyruvate Kinase and ATPases in Heart Muscle

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