A (13)C NMR double-labeling method to quantitate local myocardial O(2) consumption using frozen tissue samples

Understanding the physiology of the ageing individual: computational modelling of changes in metabolism and endurance

Ageing and lifespan are strongly affected by metabolism. The maximal possible uptake of oxygen is not only a good predictor of performance in endurance sports, but also of life expectancy. Figuratively speaking, healthy ageing is a competitive sport. Although the root cause of ageing is damage to macromolecules, it is the balance with repair processes that is decisive. Reduced or intermittent nutrition, hormones and intracellular signalling pathways that regulate metabolism have strong effects on ageing. Homeostatic regulatory processes tend to keep the environment of the cells within relatively narrow bounds. On the other hand, the body is constantly adapting to physical activity and food consumption. Spontaneous fluctuations in heart rate and other processes indicate youth and health. A (homeo)dynamic aspect of homeostasis deteriorates with age. We are now in a position to develop computational models of human metabolism and the dynamics of heart rhythm and oxygen transport that will advance our understanding of ageing. Computational modelling of the connections between dietary restriction, metabolism and protein turnover may increase insight into homeostasis of the proteins in our body. In this way, the computational reconstruction of human physiological processes, the Physiome, can help prevent frailty and age-related disease.

Progressively heterogeneous mismatch of regional oxygen delivery to consumption during graded coronary stenosis in pig left ventricle

In normal hearts, myocardial perfusion is fairly well matched to regional metabolic demand, although both are distributed heterogeneously. Nonuniform regional metabolic vulnerability during coronary stenosis would help to explain nonuniform necrosis during myocardial infarction. In the present study, we investigated whether metabolism-perfusion correlation diminishes during coronary stenosis, indicating increasing mismatch of regional oxygen supply to demand. Thirty anesthetized male pigs were studied: controls without coronary stenosis ( n = 11); group I, left anterior descending (LAD) coronary stenosis leading to coronary perfusion pressure reduction to 70 mmHg ( n = 6); group II, stenosis with perfusion pressure of about 35 mmHg ( n = 6); and group III, stenosis with perfusion pressure of 45 mmHg combined with adenosine infusion ( n = 7). [2-13C]- and [1,2-13C]acetate infusion was used to calculate regional O2 consumption from glutamate NMR spectra measured for multiple tissue samples of about 100 mg dry mass in the LAD region. Blood flow was measured with microspheres in the same regions. In control hearts without stenosis, regional oxygen extraction did not correlate with basal blood flow. Average myocardial O2 delivery and consumption decreased during coronary stenosis, but vasodilation with adenosine counteracted this. Regional oxygen extraction was on average decreased during stenosis, suggesting adaptation of metabolism to lower oxygen supply after half an hour of ischemia. Whereas regional O2 delivery correlated with O2 consumption in controls, this relation was progressively lost with graded coronary hypotension but partially reestablished by adenosine infusion. Therefore, coronary stenosis leads to heterogeneous metabolic stress indicated by decreasing regional O2 supply to demand matching in myocardium during partial coronary obstruction.

Computational estimation of tricarboxylic acid cycle fluxes using noisy NMR data from cardiac biopsies

Background

The aerobic energy metabolism of cardiac muscle cells is of major importance for the contractile function of the heart. Because energy metabolism is very heterogeneously distributed in heart tissue, especially during coronary disease, a method to quantify metabolic fluxes in small tissue samples is desirable. Taking tissue biopsies after infusion of substrates labeled with stable carbon isotopes makes this possible in animal experiments. However, the appreciable noise level in NMR spectra of extracted tissue samples makes computational estimation of metabolic fluxes challenging and a good method to define confidence regions was not yet available.

Results

Here we present a computational analysis method for nuclear magnetic resonance (NMR) measurements of tricarboxylic acid (TCA) cycle metabolites. The method was validated using measurements on extracts of single tissue biopsies taken from porcine heart in vivo. Isotopic enrichment of glutamate was measured by NMR spectroscopy in tissue samples taken at a single time point after the timed infusion of ¹³C labeled substrates for the TCA cycle. The NMR intensities for glutamate were analyzed with a computational model describing carbon transitions in the TCA cycle and carbon exchange with amino acids. The model dynamics depended on five flux parameters, which were optimized to fit the NMR measurements. To determine confidence regions for the estimated fluxes, we used the Metropolis-Hastings algorithm for Markov chain Monte Carlo (MCMC) sampling to generate extensive ensembles of feasible flux combinations that describe the data within measurement precision limits. To validate our method, we compared myocardial oxygen consumption calculated from the TCA cycle flux with in vivo blood gas measurements for 38 hearts under several experimental conditions, e.g. during coronary artery narrowing.

Conclusions

Despite the appreciable NMR noise level, the oxygen consumption in the tissue samples, estimated from the NMR spectra, correlates with blood-gas oxygen uptake measurements for the whole heart. The MCMC method provides confidence regions for the estimated metabolic fluxes in single cardiac biopsies, taking the quantified measurement noise level and the nonlinear dependencies between parameters fully into account.

Measuring non-steady-state metabolic fluxes in starch-converting faecal microbiota in vitro

This paper explores human gut bacterial metabolism of starch using a combined analytical and computational modelling approach for metabolite and flux analysis. Non-steady-state isotopic labelling experiments were performed with human faecal microbiota in a well-established in vitro model of the human colon. After culture stabilisation, [U-13C] starch was added and samples were taken at regular intervals. Metabolite concentrations and 13C isotopomeric distributions were measured amongst other things for acetate, propionate and butyrate by mass spectrometry and NMR. The vast majority of metabolic flux analysis methods based on isotopomer analysis published to date are not applicable to metabolic non-steady-state experiments. We therefore developed a new ordinary differential equation-based representation of a metabolic model of human faecal microbiota to determine eleven metabolic parameters that characterised the metabolic flux distribution in the isotope labelling experiment. The feasibility of the model parameter quantification was demonstrated on noisy in silico data using a downhill simplex optimisation, matching simulated labelling patterns of isotopically labelled metabolites with measured metabolite and isotope labelling data. Using the experimental data, we determined an increasing net label influx from starch during the experiment from 94±1 µmol/l/min to 133±3 µmol/l/min. Only about 12% of the total carbon flux from starch reached propionate. Propionate production mainly proceeded via succinate with a small contribution via acrylate. The remaining flux from starch yielded acetate (35%) and butyrate (53%). Interpretation of 13C NMR multiplet signals further revealed that butyrate, valerate and caproate were mainly synthesised via cross-feeding, using acetate as a co-substrate. This study demonstrates for the first time that the experimental design and the analysis of the results by computational modelling allows the determination of time-resolved effects of nutrition on the flux distribution within human faecal microbiota in metabolic non-steady-state.

Robust modelling, measurement and analysis of human and animal metabolic systems

Modelling human and animal metabolism is impeded by the lack of accurate quantitative parameters and the large number of biochemical reactions. This problem may be tackled by: (i) study of modules of the network independently; (ii) ensemble simulations to explore many plausible parameter combinations; (iii) analysis of 'sloppy' parameter behaviour, revealing interdependent parameter combinations with little influence; (iv) multiscale analysis that combines molecular and whole network data; and (v) measuring metabolic flux (rate of flow) in vivo via stable isotope labelling. For the latter method, carbon transition networks were modelled with systems of ordinary differential equations, but we show that coloured Petri nets provide a more intuitive graphical approach. Analysis of parameter sensitivities shows that only a few parameter combinations have a large effect on predictions. Model analysis of high-energy phosphate transport indicates that membrane permeability, inaccurately known at the organellar level, can be well determined from whole-organ responses. Ensemble simulations that take into account the imprecision of measured molecular parameters contradict the popular hypothesis that high-energy phosphate transport in heart muscle is mostly by phosphocreatine. Combining modular, multiscale, ensemble and sloppy modelling approaches with in vivo flux measurements may prove indispensable for the modelling of the large human metabolic system.

Regional sympathetic denervation affects the relation between canine local myocardial blood flow and oxygen consumption

Myocardial blood flow and oxygen consumption are heterogeneously distributed. Perfusion and myocardial oxygen consumption are closely correlated in the normal heart. It is unknown how this metabolism-perfusion relation is influenced by sympathetic denervation. We investigated this question in seven chloralose-anaesthetized dogs, 3-4 weeks after regional sympathetic denervation of the left circumflex coronary artery area of supply of the left ventricle. Measurements were made of local myocardial blood flow (MBF, in ml min(-1) (g dry wt)(-1)), measured with microspheres, and myocardial oxygen consumption ( , in mumol min(-1) (g dry wt)(-1)) in the same location, calculated from the (13)C spectrum of tissue extracts after intracoronary infusion of 3-(13)C-lactate. Since both innervated and denervated regions are subject to the same arterial pressure, lower blood flow indicates higher resistance. Mean MBF was 5.56 ml min(-1) (g dry wt)(-1) (heterogeneity of 3.47 ml min(-1) (g dry wt)(-1)) innervated, 7.48 ml min(-1) (g dry wt)(-1) (heterogeneity of 3.62 ml min(-1) (g dry wt)(-1)) denervated (n.s.). Significant linear relations were found between MBF and M Vo2 of individual samples within the innervated and denervated regions. The slopes of these relations were not significantly different, but the adjusted mean was significantly higher in the denervated regions (+1.92 ml min(-1) (g dry wt)(-1), an increase of 38% of the mean MBF at the pooled mean M Vo2, P = 0.028, ANCOVA). The ratio MBF/M Vo2(in ml micromol(-1)) was significantly higher, being 0.296 +/- 0.167 ml micromol(-1) in the denervated region compared with the innervated region, 0.216 +/- 0.126 ml micromol(-1), P = 0.0182, Mann-Whitney U test. These results indicate that sympathetic tone under chloralose anaesthesia imposes a moderate vasoconstrictive effect in the myocardium that is not detected by comparison of the mean blood flow or resistance.

Basal and hyperaemic myocardial blood flow in regionally denervated canine hearts: an in vivo study with positron emission tomography

PURPOSE

Positron emission tomography (PET) studies in patients with diabetic autonomic neuropathy (DAN) have demonstrated the impact of this disease on cardiac sympathetic innervation and myocardial blood flow (MBF). To investigate the effects of selective partial sympathetic denervation of the left ventricle (LV) on baseline and hyperaemic MBF, we measured myocardial presynaptic catecholamine re-uptake (uptake-1), beta-adrenoceptor (beta-AR) density and MBF non-invasively by means of PET in a canine model of regional sympathetic denervation.

METHODS

In 11 anaesthetised dogs, the sympathetic nerves of the free wall and septum of the LV were removed by means of dissection and phenol painting. Three weeks later, the animals were studied with PET. MBF was measured at baseline and following i.v. adenosine (140 microg kg(-1) min(-1)) and dobutamine (20 microg kg(-1) min(-1)) using(15)O-labelled water. Sympathetic denervation was confirmed by an 80+/-12% decrease in the volume of distribution (V(d)) of [(11)C]hydroxyephedrine (HED) compared with innervated regions. Myocardial beta-AR density was measured using [(11)C]CGP12177.

RESULTS

Innervated and denervated regions showed no differences in MBF at baseline and during adenosine or dobutamine. [(11)C]HED V(d)was inversely correlated with MBF in both regions at baseline, and the correlation was lost during hyperaemia in denervated regions. However, for any given value of MBF, [(11)C]HED V(d)was significantly lower in the denervated regions. beta-AR density was comparable in denervated and innervated regions (17.9+/-4.2 vs 18.4+/-3.3 pmol g(-1); p=NS).

CONCLUSION

In this experimental model, selective, regional sympathetic denervation of the LV, which results in a profound reduction in [(11)C]HED V(d), did not affect baseline or hyperaemic MBF. In addition, we demonstrated that, under baseline conditions, there was a significant inverse correlation between [(11)C]HED V(d)and MBF in both denervated and innervated regions.

Myocardial O2 consumption in porcine left ventricle is heterogeneously distributed in parallel to heterogeneous O2 delivery

Myocardial blood flow is unevenly distributed, but the cause of this heterogeneity is unknown. Heterogeneous blood flow may reflect heterogeneity of oxygen demand. The aim of the present study was to assess the relation between oxygen consumption and blood flow in small tissue regions in porcine left ventricle. In seven male, anesthetized, open-chest pigs, local oxygen consumption was quantitated by computational model analysis of the incorporation of 13C in glutamate via the tricarboxylic acid cycle during timed infusion of [13C]acetate into the left anterior descending coronary artery. Blood flow was measured with radioactive microspheres before and during acetate infusion. High-resolution nuclear magnetic resonance 13C spectra were obtained from extracts of tissue samples (159 mg mean dry wt) taken at the end of the acetate infusion. Mean regional myocardial blood flow was stable [5.0 ± 1.6 (SD) and 5.0 ± 1.4 ml·min−1·g dry wt−1 before and after 30 min of acetate infusion, respectively]. Mean left ventricular oxygen consumption measured with the NMR method was 18.6 ± 7.7 μmol·min−1·g dry wt−1 and correlated well ( r = 0.85, P = 0.02, n = 7) with oxygen consumption calculated from blood flow, hemoglobin, and blood gas measurements (mean 22.8 ± 4.7 μmol·min−1·g dry wt−1). Local blood flow and oxygen consumption were significantly correlated ( r = 0.63 for pooled normalized data, P < 0.0001, n = 60). We calculate that, in the heart at normal workload, the variance of left ventricular oxygen delivery at submilliliter resolution is explained for 43% by heterogeneity in oxygen demand.

Low-O(2) affinity erythrocytes improve performance of ischemic myocardium

O(2) transport and O(2) diffusion interact in providing O(2) to tissue, but the extent to which diffusion may be critical in the heart is unclear. If O(2) diffusion limits mitochondrial oxygenation, a change in blood O(2) affinity at constant total O(2) transport should alter cardiac O(2) consumption (VO(2)) and function. To test this hypothesis, we perfused isolated isovolumically working rabbit hearts with erythrocytes at physiological blood-gas values and P(50) (PO(2) required to half-saturate hemoglobin) values at pH of 7.4 of 17 +/- 1 Torr (2,3-bisphosphoglycerate depletion) and 33 +/- 5 Torr (inositol hexaphosphate incorporation). When perfused at 40 and 20% of normal coronary flow, mean VO(2) decreased from the control value by 37 and 46% (P < 0.001), and function, expressed as cardiac work, decreased by 38 and 52%, respectively (P < 0.001). Perfusion at higher P(50) during low-flow ischemia improved VO(2) by 20% (P < 0.001) and function by 36% (P < 0.02). There was also modest improvement at basal flow (P < 0.02 and P < 0.002, respectively). The improvement in VO(2) and function due to the P(50) increase demonstrates the importance of O(2) diffusion in this cardiac ischemia model.

The mechanical and metabolic basis of myocardial blood flow heterogeneity

Precise measurements of regional myocardial blood flow heterogeneity had to be developed before one could seek causation for the heterogeneity. Deposition techniques (particles or molecular microspheres) are the most precise, but imaging techniques have begun to provide high enough resolution to allow in vivo studies. Assigning causation has been difficult. There is no apparent association with the regional concentrations of energy-related enzymes or substrates, but these are measures of status, not of metabolism. There is statistical correlation between flow and regional substrate uptake and utilization. Attribution of regional flow variation to vascular anatomy or to vasomotor control appears not to be causative on a long-term basis. The closest relationships appear to be with mechanical function, but one cannot say for sure whether this is related to ATP hydrolysis at the crossbridge or associated metabolic reactions such as calcium uptake by the sarcoplasmic reticulum.

New double-tracer digital radiography for analysis of spatial and temporal myocardial flow heterogeneity

A new high-resolution digital radiographic technique based on the deposition of (125)I- and (3)H-labeled desmethylimipramine (IDMI and HDMI, respectively) was developed for the assessment of spatial and temporal myocardial flow heterogeneity at a microvascular level. The density distributions of two tracers, or relative flow distributions, were determined by subtraction digital radiography using two imaging plates of different sensitivity. The regions resolved are comparable in size to vascular regulatory units (400 x 400 microm(2)). This method was applied to the measurement of within-layer myocardial flow distributions in Langendorff-perfused rabbit hearts. The validity of this method was confirmed by the strong correlation between regional densities of two tracers injected simultaneously (r = 0.89 +/- 0.03, n = 8). The temporal flow stability was evaluated by a 90-s continuous IDMI injection and subsequent bolus HDMI injection (n = 8). Regional densities of the two tracers were fairly correlated (r = 0.86 +/- 0.03), indicating that the spatial pattern of flow distribution was stable even at a microvascular level over a 90-s period. The effect of microsphere embolization on the flow distribution was also investigated by the sequential injections of IDMI, 15-microm microspheres, and HDMI at 20-s intervals (n = 8). Microembolization increased the coefficient of variation of tracer density from 19 to 25% (P < 0.05), whereas the regional densities of two tracers were still correlated substantially, as in the case of no embolization (r = 0.84 +/- 0.06). Thus the microsphere embolization enhanced flow heterogeneity with increasing flow differences between control high-flow and control low-flow regions but rather maintained the pattern of flow distribution. In conclusion, double-tracer digital radiography will be a promising method for the spatial and temporal myocardial flow analysis at microvascular levels.