ADP compartmentation analysis reveals coupling between pyruvate kinase and ATPases in heart muscle

A novel Nav1.5-dependent feedback mechanism driving glycolytic acidification in breast cancer metastasis

Solid tumours have abnormally high intracellular [Na+]. The activity of various Na+ channels may underlie this Na+ accumulation. Voltage-gated Na+ channels (VGSCs) have been shown to be functionally active in cancer cell lines, where they promote invasion. However, the mechanisms involved, and clinical relevance, are incompletely understood. Here, we show that protein expression of the Nav1.5 VGSC subtype strongly correlates with increased metastasis and shortened cancer-specific survival in breast cancer patients. In addition, VGSCs are functionally active in patient-derived breast tumour cells, cell lines, and cancer-associated fibroblasts. Knockdown of Nav1.5 in a mouse model of breast cancer suppresses expression of invasion-regulating genes. Nav1.5 activity increases ATP demand and glycolysis in breast cancer cells, likely by upregulating activity of the Na+/K+ ATPase, thus promoting H+ production and extracellular acidification. The pH of murine xenograft tumours is lower at the periphery than in the core, in regions of higher proliferation and lower apoptosis. In turn, acidic extracellular pH elevates persistent Na+ influx through Nav1.5 into breast cancer cells. Together, these findings show positive feedback between extracellular acidification and the movement of Na+ into cancer cells which can facilitate invasion. These results highlight the clinical significance of Nav1.5 activity as a potentiator of breast cancer metastasis and provide further evidence supporting the use of VGSC inhibitors in cancer treatment.

Rat and mouse cardiomyocytes show subtle differences in creatine kinase expression and compartmentalization

Creatine kinase (CK) and adenylate kinase (AK) are energy transfer systems. Different studies on permeabilized cardiomyocytes suggest that ADP-channelling from mitochondrial CK alone stimulates respiration to its maximum, VO2_max, in rat but not mouse cardiomyocytes. Results are ambiguous on ADP-channelling from AK to mitochondria. This study was undertaken to directly compare the CK and AK systems in rat and mouse hearts. In homogenates, we assessed CK- and AK-activities, and the CK isoform distribution. In permeabilized cardiomyocytes, we assessed mitochondrial respiration stimulated by ADP from CK and AK, VO2_CK and VO2_AK, respectively. The ADP-channelling from CK or AK to mitochondria was assessed by adding PEP and PK to competitively inhibit the respiration rate. We found that rat compared to mouse hearts had a lower aerobic capacity, higher VO2_CK/VO2_max, and different CK-isoform distribution. Although rat hearts had a larger fraction of mitochondrial CK, less ADP was channeled from CK to the mitochondria. This suggests different intracellular compartmentalization in rat and mouse cardiomyocytes. VO2_AK/VO2_max was similar in mouse and rat cardiomyocytes, and AK did not channel ADP to the mitochondria. In the absence of intracellular compartmentalization, the AK- and CK-activities in homogenate should have been similar to the ADP-phosphorylation rates estimated from VO2_AK and VO2_CK in permeabilized cardiomyocytes. Instead, we found that the ADP-phosphorylation rates estimated from permeabilized cardiomyocytes were 2 and 9 times lower than the activities recorded in homogenate for CK and AK, respectively. Our results highlight the importance of energetic compartmentalization in cardiac metabolic regulation and signalling.

IOCBIO Kinetics: An open-source software solution for analysis of data traces

Biological measurements frequently involve measuring parameters as a function of time, space, or frequency. Later, during the analysis phase of the study, the researcher splits the recorded data trace into smaller sections, analyzes each section separately by finding a mean or fitting against a specified function, and uses the analysis results in the study. Here, we present the software that allows to analyze these data traces in a manner that ensures repeatability of the analysis and simplifies the application of FAIR (findability, accessibility, interoperability, and reusability) principles in such studies. At the same time, it simplifies the routine data analysis pipeline and gives access to a fast overview of the analysis results. For that, the software supports reading the raw data, processing the data as specified in the protocol, and storing all intermediate results in the laboratory database. The software can be extended by study- or hardware-specific modules to provide the required data import and analysis facilities. To simplify the development of the data entry web interfaces, that can be used to enter data describing the experiments, we released a web framework with an example implementation of such a site. The software is covered by open-source license and is available through several online channels.

In biological and other types of experiments, we frequently record changes of some parameters in time or space. It is common to analyze the data by splitting the recording into smaller sections and relating it to some changes induced by the researchers. The steps involved in the analysis are: splitting of the data, fitting them to some function, relating the fit result to the change in the environment, and normalization. These steps are frequently done through several software packages, spreedsheets, and manual copy and paste between the programs. The software presented in this work allows to make all these analysis steps in one database in a manner that is easy, can be reproduced by others, and clearly tracks the history of all the analysis steps. In addition, it allows to link the experimental data with the description of the experiment, making it simple to perform tasks such as normalization of the recorded values, relating experimental recordings to the sample or animal, as well as extracting data from the laboratory database for publishing. The software is written to be easily extendable by user-defined modules to fit the analysis pipelines and is expected to improve the data analysis practices in research.

Respiration of permeabilized cardiomyocytes from mice: no sex differences, but substrate-dependent changes in the apparent ADP-affinity

Sex differences in cardiac physiology are getting increased attention. This study assessed whether isolated, permeabilized cardiomyocytes from male and female C57BL/6 mice differ in terms of their respiration with multiple substrates and overall intracellular diffusion restriction estimated by the apparent ADP-affinity of respiration. Using respirometry, we recorded 1) the activities of respiratory complexes I, II and IV, 2) the respiration rate with substrates fuelling either complex I, II, or I + II, and 3) the apparent ADP-affinity with substrates fuelling complex I and I + II. The respiration rates were normalized to protein content and citrate synthase (CS) activity. We found no sex differences in CS activity (a marker of mitochondrial content) normalized to protein content or in any of the respiration measurements. This suggests that cardiomyocytes from male and female mice do not differ in terms of mitochondrial respiratory capacity and apparent ADP-affinity. Pyruvate modestly lowered the respiration rate, when added to succinate, glutamate and malate. This may be explained by intramitochondrial compartmentalization caused by the formation of supercomplexes and their association with specific dehydrogenases. To our knowledge, we show for the first time that the apparent ADP-affinity was substrate-dependent. This suggests that substrates may change or regulate intracellular barriers in cardiomyocytes.

Mitochondrial function remains impaired in the hypertrophied right ventricle of pulmonary hypertensive rats following short duration metoprolol treatment

Pulmonary hypertension (PH) increases the work of the right ventricle (RV) and causes right-sided heart failure. This study examined RV mitochondrial function and ADP transfer in PH animals advancing to right heart failure, and investigated a potential therapy with the specific β1-adrenergic-blocker metoprolol. Adult Wistar rats (317 ± 4 g) were injected either with monocrotaline (MCT, 60 mg kg-1) to induce PH, or with an equivalent volume of saline for controls (CON). At three weeks post-injection the MCT rats began oral metoprolol (10 mg kg-1 day-1-) or placebo treatment until heart failure was observed in the MCT group. Mitochondrial function was then measured using high-resolution respirometry from permeabilised RV fibres. Relative to controls, MCT animals had impaired mitochondrial function but maintained coupling between myofibrillar ATPases and mitochondria, despite an increase in ADP diffusion distances. Cardiomyocytes from the RV of MCT rats were enlarged, primarily due to an increase in myofibrillar protein. The ratio of mitochondria per myofilament area was decreased in both MCT groups (p ≤ 0.05) in comparison to control (CON: 1.03 ± 0.04; MCT: 0.74 ± 0.04; MCT + BB: 0.74 ± 0.03). This not only implicates impaired energy production in PH, but also increases the diffusion distance for metabolites within the MCT cardiomyocytes, adding an additional hindrance to energy supply. Together, these changes may limit energy supply in MCT rat hearts, particularly at high cardiac workloads. Metoprolol treatment did not delay the onset of heart failure symptoms, improve mitochondrial function, or regress RV hypertrophy.

Proteomic analysis of rat serum revealed the effects of chronic sleep deprivation on metabolic, cardiovascular and nervous system

Sleep is an essential and fundamental physiological process that plays crucial roles in the balance of psychological and physical health. Sleep disorder may lead to adverse health outcomes. The effects of sleep deprivation were extensively studied, but its mechanism is still not fully understood. The present study aimed to identify the alterations of serum proteins associated with chronic sleep deprivation, and to seek for potential biomarkers of sleep disorder mediated diseases. A label-free quantitative proteomics technology was used to survey the global changes of serum proteins between normal rats and chronic sleep deprivation rats. A total of 309 proteins were detected in the serum samples and among them, 117 proteins showed more than 1.8-folds abundance alterations between the two groups. Functional enrichment and network analyses of the differential proteins revealed a close relationship between chronic sleep deprivation and several biological processes including energy metabolism, cardiovascular function and nervous function. And four proteins including pyruvate kinase M1, clusterin, kininogen1 and profilin-1were identified as potential biomarkers for chronic sleep deprivation. The four candidates were validated via parallel reaction monitoring (PRM) based targeted proteomics. In addition, protein expression alteration of the four proteins was confirmed in myocardium and brain of rat model. In summary, the comprehensive proteomic study revealed the biological impacts of chronic sleep deprivation and discovered several potential biomarkers. This study provides further insight into the pathological and molecular mechanisms underlying sleep disorders at protein level.

Metabolic compartmentation in rainbow trout cardiomyocytes: coupling of hexokinase but not creatine kinase to mitochondrial respiration

Rainbow trout (Oncorhynchus mykiss) cardiomyocytes have a simple morphology with fewer membrane structures such as sarcoplasmic reticulum and t-tubules penetrating the cytosol. Despite this, intracellular ADP diffusion is restricted. Intriguingly, although diffusion is restricted, trout cardiomyocytes seem to lack the coupling between mitochondrial creatine kinase (CK) and respiration. Our aim was to study the distribution of diffusion restrictions in permeabilized trout cardiomyocytes and verify the role of CK. We found a high activity of hexokinase (HK), which led us to reassess the situation in trout cardiomyocytes. We show that diffusion restrictions are more prominent than previously thought. In the presence of a competitive ADP-trapping system, ADP produced by HK, but not CK, was channeled to the mitochondria. In agreement with this, we found no positively charged mitochondrial CK in trout heart homogenate. The results were best fit by a simple mathematical model suggesting that trout cardiomyocytes lack a functional coupling between ATPases and pyruvate kinase. The model simulations show that diffusion is restricted to almost the same extent in the cytosol and by the outer mitochondrial membrane. Furthermore, they confirm that HK, but not CK, is functionally coupled to respiration. In perspective, our results suggest that across a range of species, cardiomyocyte morphology and metabolism go hand in hand with cardiac performance, which is adapted to the circumstances. Mitochondrial CK is coupled to respiration in adult mammalian hearts, which are specialized to high, sustained performance. HK associates with mitochondria in hearts of trout and neonatal mammals, which are more hypoxia-tolerant.

Restricted ADP movement in cardiomyocytes: Cytosolic diffusion obstacles are complemented with a small number of open mitochondrial voltage-dependent anion channels

Adequate intracellular energy transfer is crucial for proper cardiac function. In energy starved failing hearts, partial restoration of energy transfer can rescue mechanical performance. There are two types of diffusion obstacles that interfere with energy transfer from mitochondria to ATPases: mitochondrial outer membrane (MOM) with voltage-dependent anion channel (VDAC) permeable to small hydrophilic molecules and cytoplasmatic diffusion barriers grouping ATP-producers and -consumers. So far, there is no method developed to clearly distinguish the contributions of cytoplasmatic barriers and MOM to the overall diffusion restriction. Furthermore, the number of open VDACs in vivo remains unknown. The aim of this work was to establish the partitioning of intracellular diffusion obstacles in cardiomyocytes. We studied the response of mitochondrial oxidative phosphorylation of permeabilized rat cardiomyocytes to changes in extracellular ADP by recording 3D image stacks of NADH autofluorescence. Using cell-specific mathematical models, we determined the permeability of MOM and cytoplasmatic barriers. We found that only ~2% of VDACs are accessible to cytosolic ADP and cytoplasmatic diffusion barriers reduce the apparent diffusion coefficient by 6-10×. In cardiomyocytes, diffusion barriers in the cytoplasm and by the MOM restrict ADP/ATP diffusion to similar extents suggesting a major role of both barriers in energy transfer and other intracellular processes.

The location of energetic compartments affects energetic communication in cardiomyocytes

The heart relies on accurate regulation of mitochondrial energy supply to match energy demand. The main regulators are Ca²⁺ and feedback of ADP and P_i. Regulation via feedback has intrigued for decades. First, the heart exhibits a remarkable metabolic stability. Second, diffusion of ADP and other molecules is restricted specifically in heart and red muscle, where a fast feedback is needed the most. To explain the regulation by feedback, compartmentalization must be taken into account. Experiments and theoretical approaches suggest that cardiomyocyte energetic compartmentalization is elaborate with barriers obstructing diffusion in the cytosol and at the level of the mitochondrial outer membrane (MOM). A recent study suggests the barriers are organized in a lattice with dimensions in agreement with those of intracellular structures. Here, we discuss the possible location of these barriers. The more plausible scenario includes a barrier at the level of MOM. Much research has focused on how the permeability of MOM itself is regulated, and the importance of the creatine kinase system to facilitate energetic communication. We hypothesize that at least part of the diffusion restriction at the MOM level is not by MOM itself, but due to the close physical association between the sarcoplasmic reticulum (SR) and mitochondria. This will explain why animals with a disabled creatine kinase system exhibit rather mild phenotype modifications. Mitochondria are hubs of energetics, but also ROS production and signaling. The close association between SR and mitochondria may form a diffusion barrier to ADP added outside a permeabilized cardiomyocyte. But in vivo, it is the structural basis for the mitochondrial-SR coupling that is crucial for the regulation of mitochondrial Ca²⁺-transients to regulate energetics, and for avoiding Ca²⁺-overload and irreversible opening of the mitochondrial permeability transition pore.

Tight coupling of Na+/K+-ATPase with glycolysis demonstrated in permeabilized rat cardiomyocytes

The effective integrated organization of processes in cardiac cells is achieved, in part, by the functional compartmentation of energy transfer processes. Earlier, using permeabilized cardiomyocytes, we demonstrated the existence of tight coupling between some of cardiomyocyte ATPases and glycolysis in rat. In this work, we studied contribution of two membrane ATPases and whether they are coupled to glycolysis--sarcoplasmic reticulum Ca2+ ATPase (SERCA) and plasmalemma Na+/K+-ATPase (NKA). While SERCA activity was minor in this preparation in the absence of calcium, major role of NKA was revealed accounting to ∼30% of the total ATPase activity which demonstrates that permeabilized cell preparation can be used to study this pump. To elucidate the contribution of NKA in the pool of ATPases, a series of kinetic measurements was performed in cells where NKA had been inhibited by 2 mM ouabain. In these cells, we recorded: ADP- and ATP-kinetics of respiration, competition for ADP between mitochondria and pyruvate kinase (PK), ADP-kinetics of endogenous PK, and ATP-kinetics of total ATPases. The experimental data was analyzed using a series of mathematical models with varying compartmentation levels. The results show that NKA is tightly coupled to glycolysis with undetectable flux of ATP between mitochondria and NKA. Such tight coupling of NKA to PK is in line with its increased importance in the pathological states of the heart when the substrate preference shifts to glucose.

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

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

Real-time determination of sarcomere length of a single cardiomyocyte during contraction

Sarcomere length of a cardiomyocyte is an important control parameter for physiology studies on a single cell level; for instance, its accurate determination in real time is essential for performing single cardiomyocyte contraction experiments. The aim of this work is to develop an efficient and accurate method for estimating a mean sarcomere length of a contracting cardiomyocyte using microscopy images as an input. The novelty in developed method lies in 1) using unbiased measure of similarities to eliminate systematic errors from conventional autocorrelation function (ACF)-based methods when applied to region of interest of an image, 2) using a semianalytical, seminumerical approach for evaluating the similarity measure to take into account spatial dependence of neighboring image pixels, and 3) using a detrend algorithm to extract the sarcomere striation pattern content from the microscopy images. The developed sarcomere length estimation procedure has superior computational efficiency and estimation accuracy compared with the conventional ACF and spectral analysis-based methods using fast Fourier transform. As shown by analyzing synthetic images with the known periodicity, the estimates obtained by the developed method are more accurate at the subpixel level than ones obtained using ACF analysis. When applied in practice on rat cardiomyocytes, our method was found to be robust to the choice of the region of interest that may 1) include projections of carbon fibers and nucleus, 2) have uneven background, and 3) be slightly disoriented with respect to average direction of sarcomere striation pattern. The developed method is implemented in open-source software.

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

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

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

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.

Molecular dynamics simulations of creatine kinase and adenine nucleotide translocase in mitochondrial membrane patch

Interaction between mitochondrial creatine kinase (MtCK) and adenine nucleotide translocase (ANT) can play an important role in determining energy transfer pathways in the cell. Although the functional coupling between MtCK and ANT has been demonstrated, the precise mechanism of the coupling is not clear. To study the details of the coupling, we turned to molecular dynamics simulations. We introduce a new coarse-grained molecular dynamics model of a patch of the mitochondrial inner membrane containing a transmembrane ANT and an MtCK above the membrane. The membrane model consists of three major types of lipids (phosphatidylcholine, phosphatidylethanolamine, and cardiolipin) in a roughly 2:1:1 molar ratio. A thermodynamics-based coarse-grained force field, termed MARTINI, has been used together with the GROMACS molecular dynamics package for all simulated systems in this work. Several physical properties of the system are reproduced by the model and are in agreement with known data. This includes membrane thickness, dimension of the proteins, and diffusion constants. We have studied the binding of MtCK to the membrane and demonstrated the effect of cardiolipin on the stabilization of the binding. In addition, our simulations predict which part of the MtCK protein sequence interacts with the membrane. Taken together, the model has been verified by dynamical and structural data and can be used as the basis for further studies.

Permeabilized rat cardiomyocyte response demonstrates intracellular origin of diffusion obstacles

Intracellular diffusion restrictions for ADP and other molecules have been predicted earlier based on experiments on permeabilized fibers or cardiomyocytes. However, it is possible that the effective diffusion distance is larger than the cell dimensions due to clumping of cells and incomplete separation of cells in fiber preparations. The aim of this work was to check whether diffusion restrictions exist inside rat cardiomyocytes or are caused by large effective diffusion distance. For that, we determined the response of oxidative phosphorylation (OxPhos) to exogenous ADP and ATP stimulation in permeabilized rat cardiomyocytes using fluorescence microscopy. The state of OxPhos was monitored via NADH and flavoprotein autofluorescence. By varying the ADP or ATP concentration in flow chamber, we determined that OxPhos has a low affinity in cardiomyocytes. The experiments were repeated in a fluorometer on cardiomyocyte suspensions leading to similar autofluorescence changes induced by ADP as recorded under the microscope. ATP stimulated OxPhos more in a fluorometer than under the microscope, which was attributed to accumulation of ADP in fluorometer chamber. By calculating the flow profile around the cell in the microscope chamber and comparing model solutions to measured data, we demonstrate that intracellular structures impose significant diffusion obstacles in rat cardiomyocytes.

Analyzing the functional properties of the creatine kinase system with multiscale 'sloppy' modeling

In this study the function of the two isoforms of creatine kinase (CK; EC 2.7.3.2) in myocardium is investigated. The 'phosphocreatine shuttle' hypothesis states that mitochondrial and cytosolic CK plays a pivotal role in the transport of high-energy phosphate (HEP) groups from mitochondria to myofibrils in contracting muscle. Temporal buffering of changes in ATP and ADP is another potential role of CK. With a mathematical model, we analyzed energy transport and damping of high peaks of ATP hydrolysis during the cardiac cycle. The analysis was based on multiscale data measured at the level of isolated enzymes, isolated mitochondria and on dynamic response times of oxidative phosphorylation measured at the whole heart level. Using 'sloppy modeling' ensemble simulations, we derived confidence intervals for predictions of the contributions by phosphocreatine (PCr) and ATP to the transfer of HEP from mitochondria to sites of ATP hydrolysis. Our calculations indicate that only 15±8% (mean±SD) of transcytosolic energy transport is carried by PCr, contradicting the PCr shuttle hypothesis. We also predicted temporal buffering capabilities of the CK isoforms protecting against high peaks of ATP hydrolysis (3750 µM*s(-1)) in myofibrils. CK inhibition by 98% in silico leads to an increase in amplitude of mitochondrial ATP synthesis pulsation from 215±23 to 566±31 µM*s(-1), while amplitudes of oscillations in cytosolic ADP concentration double from 77±11 to 146±1 µM. Our findings indicate that CK acts as a large bandwidth high-capacity temporal energy buffer maintaining cellular ATP homeostasis and reducing oscillations in mitochondrial metabolism. However, the contribution of CK to the transport of high-energy phosphate groups appears limited. Mitochondrial CK activity lowers cytosolic inorganic phosphate levels while cytosolic CK has the opposite effect.

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.

Symbolic flux analysis for genome-scale metabolic networks

Background

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

Results

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

Conclusions

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

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 (31)P 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.