Modeling non-linear kinetics of hyperpolarized [1-(13)C] pyruvate in the crystalloid-perfused rat heart

Extended Bloch-McConnell equations for mechanistic analysis of hyperpolarized 13C magnetic resonance experiments on enzyme systems

We describe an approach to formulating the kinetic master equations of the time evolution of NMR signals in reacting (bio)chemical systems. Special focus is given to studies that employ signal enhancement (hyperpolarization) methods such as dissolution dynamic nuclear polarization (dDNP) and involving nuclear spin-bearing solutes that undergo reactions mediated by enzymes and membrane transport proteins. We extend the work given in a recent presentation on this topic (Kuchel and Shishmarev, 2020) to now include enzymes with two or more substrates and various enzyme reaction mechanisms as classified by Cleland, with particular reference to non-first-order processes. Using this approach, we can address some pressing questions in the field from a theoretical standpoint. For example, why does binding of a hyperpolarized substrate to an enzyme not cause an appreciable loss of the signal from the substrate or product? Why does the concentration of an unlabelled pool of substrate, for example 12 C lactate, cause an increase in the rate of exchange of the 13 C-labelled pool? To what extent is the equilibrium position of the reaction perturbed during administration of the substrate? The formalism gives a full mechanistic understanding of the time courses derived and is of relevance to ongoing clinical trials using these techniques.

Joint image and field map estimation for multi-echo hyperpolarized 13 C metabolic imaging of the heart

PURPOSE

Image reconstruction of metabolic images from hyperpolarized ¹³ C multi-echo data acquisition is sensitive to susceptibility-induced phase offsets, which are particularly challenging in the heart. A model-based framework for joint estimation of metabolite images and field map from echo shift-encoded data is proposed. Using simulations, it is demonstrated that correction of signal spilling due to incorrect decomposition of metabolites and geometrical distortions over a wide range of off-resonance gradients is possible. In vivo feasibility is illustrated using hyperpolarized [1-¹³ C]pyruvate in the pig heart.

METHODS

The model-based reconstruction for multi-echo, multicoil data was implemented as a nonconvex minimization problem jointly optimizing for metabolic images and B₀ . A comprehensive simulation framework for echo shift-encoded hyperpolarized [1-¹³ C]pyruvate imaging was developed and applied to assess reconstruction performance and distortion correction of the proposed method. In vivo data were obtained in four pigs using hyperpolarized [1-¹³ C]pyruvate on a clinical 3T MR system with a six-channel receiver coil. Dynamic images were acquired during suspended ventilation using cardiac-triggered multi-echo single-shot echo-planar imaging in short-axis orientation.

RESULTS

Simulations revealed that off-resonance gradients up to ±0.26 ppm/pixel can be corrected for with reduced signal spilling and geometrical distortions yielding an accuracy of ≥90% in terms of Dice similarity index. In vivo, improved geometrical consistency (10% Dice improvement) compared to image reconstruction without field map correction and with reference to anatomical data was achieved.

CONCLUSION

Joint image and field map estimation allows addressing off-resonance-induced geometrical distortions and metabolite spilling in hyperpolarized ¹³ C metabolic imaging of the heart.

Correlation between lactate dehydrogenase/pyruvate dehydrogenase activities ratio and tissue pH in the perfused mouse heart: A potential noninvasive indicator of cardiac pH provided by hyperpolarized magnetic resonance

Cardiovascular diseases account for more than 30% of all deaths worldwide and many could be ameliorated with early diagnosis. Current cardiac imaging modalities can assess blood flow, heart anatomy and mechanical function. However, for early diagnosis and improved treatment, further functional biomarkers are needed. One such functional biomarker could be the myocardium pH. Although tissue pH is already determinable via MR techniques, and has been since the early 1990s, it remains elusive to use practically. The objective of this study was to explore the possibility to evaluate cardiac pH noninvasively, using in-cell enzymatic rates of hyperpolarized [1-¹³ C]pyruvate metabolism (ie, moles of product produced per unit time) determined directly in real time using magnetic resonance spectroscopy in a perfused mouse heart model. As a gold standard for tissue pH we used ³¹ P spectroscopy and the chemical shift of the inorganic phosphate (Pi) signal. The nonhomogenous pH distribution of the perfused heart was analyzed using a multi-parametric analysis of this signal, thus taking into account the heterogeneous nature of this characteristic. As opposed to the signal ratio of hyperpolarized [¹³ C]bicarbonate to [¹³ CO₂ ], which has shown correlation to pH in other studies, we investigated here the ratio of two intracellular enzymatic rates: lactate dehydrogenase (LDH) and pyruvate dehydrogenase (PDH), by way of determining the production rates of [1-¹³ C]lactate and [¹³ C]bicarbonate, respectively. The enzyme activities determined here are intracellular, while the pH determined using the Pi signal may contain an extracellular component, which could not be ruled out. Nevertheless, we report a strong correlation between the tissue pH and the LDH/PDH activities ratio. This work may pave the way for using the LDH/PDH activities ratio as an indicator of cardiac intracellular pH in vivo, in an MRI examination.

Development of a human heart-sized perfusion system for metabolic imaging studies using hyperpolarized [1-13 C]pyruvate MRI

PURPOSE

Increasing worldwide demand for cardiac transplantation has spurred new developments to increase the donor pool. Normothermic preservation of heart grafts for transplantation is an emerging strategy to improve the utilization of marginal grafts. Hyperpolarized MR using metabolic tracers such as [1-¹³ C]pyruvate, provide a novel means of investigating metabolic status without the use of ionizing radiation. We demonstrate the use of this methodology to examine ex vivo perfused porcine heart grafts.

METHODS

Hearts from three 40-kg Danish domestic pigs were harvested and subsequently perfused in Langendorff mode under normothermic conditions, using an MR-compatible perfusion system adapted to the heart. Proton MRI and hyperpolarized [1-¹³ C]pyruvate were used to investigate and quantify the functional and metabolic status of the grafts.

RESULTS

Hearts were perfused with whole blood for 120 min, using a dynamic contrast-enhanced perfusion experiment to verify successful myocardial perfusion. Hyperpolarized [1-¹³ C]pyruvate MRI was used to assess the metabolic state of the myocardium. Functional assessment was performed using CINE imaging and ventricular pressure data. High lactate and modest alanine levels were observed in the hyperpolarized experiment. The functional assessment produced reduced functional parameters. This suggests an altered functional and metabolic profile compared with corresponding in vivo values.

CONCLUSION

We investigated the metabolic and functional status of machine-perfused porcine hearts. Utilizing hyperpolarized methodology to acquire detailed myocardial metabolic information-in combination with already established MR methods for cardiac investigation-provides a powerful tool to aid the progress of donor heart preservation.

Hyperpolarized product selective saturating-excitations for determination of changes in metabolic reaction rates in real-time

Investigation of hyperpolarized substrate metabolism has been showing utility in real-time determination of in-cell and in vivo enzymatic activities. Intracellular reaction rates may vary during the course of a measurement, even on the very short time scales of visibility on hyperpolarized MR, due to many factors such as the availability of the substrate and co-factors in the intracellular space. Despite this potential variation, the kinetic analysis of hyperpolarized signals typically assumes that the same rate constant (and in many cases, the same rate) applies throughout the course of the reaction as observed via the build-up and decay of the hyperpolarized signals. We demonstrate here an acquisition approach that can null the need for such an assumption and enable the detection of instantaneous changes in the rate of the reaction during an ex vivo hyperpolarized investigation, (i.e. in the course of the decay of one hyperpolarized substrate dose administered to a viable tissue sample ex vivo). This approach utilizes hyperpolarized product selective saturating-excitation pulses. Similar pulses have been previously utilized in vivo for spectroscopic imaging. However, we show here favorable consequences to kinetic rate determinations in the preparations used. We implement this acquisition strategy for studies on perfused tissue slices and develop a theory that explains why this particular approach enables the determination of changes in enzymatic rates that are monitored via the chemical conversions of hyperpolarized substrates. Real-time changes in intracellular reaction rates are demonstrated in perfused brain, liver, and xenograft breast cancer tissue slices and provide another potential differentiation parameter for tissue characterization.

Kinetic Modeling of Hyperpolarized Carbon-13 Pyruvate Metabolism in the Human Brain

Kinetic modeling of the in vivo pyruvate-to-lactate conversion is crucial to investigating aberrant cancer metabolism that demonstrates Warburg effect modifications. Non-invasive detection of alterations to metabolic flux might offer prognostic value and improve the monitoring of response to treatment. In this clinical research project, hyperpolarized [1-¹³C] pyruvate was intravenously injected in a total of 10 brain tumor patients to measure its rate of conversion to lactate ( k_PL ) and bicarbonate ( k_PB ) via echo-planar imaging. Our aim was to investigate new methods to provide k_PL and k_PB maps with whole-brain coverage. The approach was data-driven and addressed two main issues: selecting the optimal model for fitting our data and determining an appropriate goodness-of-fit metric. The statistical analysis suggested that an input-less model had the best agreement with the data. It was also found that selecting voxels based on post-fitting error criteria provided improved precision and wider spatial coverage compared to using signal-to-noise cutoffs alone.

Intra-cardiac Side-Firing Light Catheter for Monitoring Cellular Metabolism using Transmural Absorbance Spectroscopy of Perfused Mammalian Hearts

Absorbance spectroscopy of cardiac muscle provides non-destructive assessment of cytosolic and mitochondrial oxygenation via myoglobin and cytochrome absorbance respectively. In addition, numerous aspects of the mitochondrial metabolic status such as membrane potential and substrate entry can also be estimated. To perform cardiac wall transmission optical spectroscopy, a commercially available side-firing optical fiber catheter is placed in the left ventricle of the isolated perfused heart as a light source. Light passing through the heart wall is collected with an external optical fiber to perform optical spectroscopy of the heart in near real- time. The transmission approach avoids numerous surface scattering interference occurring in widely used reflection approaches. Changes in transmural absorbance spectra were deconvolved using a library of chromophore reference spectra, providing quantitative measures of all the known cardiac chromophores simultaneously. This spectral deconvolution approach eliminated intrinsic errors that may result from using common dual wavelength methods applied to overlapping absorbance spectra, as well as provided a quantitative evaluation of the goodness of fit. A custom program was designed for data acquisition and analysis, which permitted the investigator to monitor the metabolic state of the preparation during the experiment. These relatively simple additions to the standard heart perfusion system provide a unique insight into the metabolic state of the heart wall in addition to conventional measures of contraction, perfusion, and substrate/oxygen extraction.

Effects of deuteration on transamination and oxidation of hyperpolarized 13C-Pyruvate in the isolated heart

This study was designed to determine the effects of deuteration in pyruvate on exchange reactions in alanine aminotransferase (ALT), lactate dehydrogenase (LDH) and flux through pyruvate dehydrogenase (PDH). Although deuteration of a ¹³C enriched substrate is commonly used to increase the lifetime of a probe for hyperpolarization experiments, the potential impact of kinetic isotope effects on such substitutions has not been studied in detail. Metabolism of deuterated pyruvate was investigated in isolated rat hearts. Hearts were perfused with a 1:1 mixture of [U-¹³C₃]pyruvate and [2-¹³C₁]pyruvate or a 1:1 mixture of [U-¹³C₃]pyruvate plus [2-¹³C₁, U-²H₃]pyruvate for 30 min before being freeze clamped. Another set of hearts received [2-¹³C₁, U-²H₃]pyruvate and was freeze-clamped at 3 min or 6 min. Tissue extracts were analyzed by ¹H and ¹³C{¹H} NMR spectroscopy. The chemical shift isotope effect of ²H was monitored in the ¹³C NMR spectra of the C2 resonance of lactate and alanine plus the C5 of glutamate. There was little kinetic isotope effect of ²H in pyruvate on flux through PDH, LDH or ALT as detected by the distribution of ¹³C, but the distribution of ²H differed markedly between alanine and lactate. At steady-state, alanine was a mixture of deuterated species, while lactate was largely perdeuterated. Consistent with results at steady-state, hearts freeze-clamped at 3 min or 6 min showed rapid removal of deuterium in alanine but not in lactate. Metabolism of hyperpolarized [1-¹³C₁]pyruvate was compared to [1-¹³C₁,U-²H₃]pyruvate in isolated hearts. Consistent with the results from tissue extracts, there was little effect of deuteration on the kinetics of appearance of lactate, alanine or bicarbonate, but there was a small, time-dependent upfield chemical shift in the HP[1-¹³C₁]alanine signal reflecting exchange of methyl deuterons with water protons. Together, these results demonstrate that (1) the kinetics of pyruvate metabolism in hearts detected by ¹³C NMR are not affected by replacement of the pyruvate methyl protons with deuterons and (2) that the loss of deuterium from the methyl position occurs rapidly during the conversion of pyruvate to alanine. The majority of the deuterium atoms are lost on the time-scale of a hyperpolarization experiment.

Quantitative myocardial first-pass cardiovascular magnetic resonance perfusion imaging using hyperpolarized [1-13C] pyruvate

BACKGROUND

The feasibility of absolute myocardial blood flow quantification and suitability of hyperpolarized [1-13C] pyruvate as contrast agent for first-pass cardiovascular magnetic resonance (CMR) perfusion measurements are investigated with simulations and demonstrated in vivo in a swine model.

METHODS

A versatile simulation framework for hyperpolarized CMR subject to physical, physiological and technical constraints was developed and applied to investigate experimental conditions for accurate perfusion CMR with hyperpolarized [1-13C] pyruvate. Absolute and semi-quantitative perfusion indices were analyzed with respect to experimental parameter variations and different signal-to-noise ratio (SNR) levels. Absolute myocardial blood flow quantification was implemented with an iterative deconvolution approach based on Fermi functions. To demonstrate in vivo feasibility, velocity-selective excitation with an echo-planar imaging readout was used to acquire dynamic myocardial stress perfusion images in four healthy swine. Arterial input functions were extracted from an additional image slice with conventional excitation that was acquired within the same heartbeat.

RESULTS

Simulations suggest that obtainable SNR and B0 inhomogeneity in vivo are sufficient for the determination of absolute and semi-quantitative perfusion with ≤25% error. It is shown that for expected metabolic conversion rates, metabolic conversion of pyruvate can be neglected over the short duration of acquisition in first-pass perfusion CMR. In vivo measurements suggest that absolute myocardial blood flow quantification using hyperpolarized [1-13C] pyruvate is feasible with an intra-myocardial variability comparable to semi-quantitative perfusion indices.

CONCLUSION

The feasibility of quantitative hyperpolarized first-pass perfusion CMR using [1-13C] pyruvate has been investigated in simulations and demonstrated in swine. Using an approved and metabolically active compound is envisioned to increase the value of hyperpolarized perfusion CMR in patients.

Investigation of analysis methods for hyperpolarized 13C-pyruvate metabolic MRI in prostate cancer patients

MRI using hyperpolarized (HP) carbon-13 pyruvate is being investigated in clinical trials to provide non-invasive measurements of metabolism for cancer and cardiac imaging. In this project, we applied HP [1-¹³ C]pyruvate dynamic MRI in prostate cancer to measure the conversion from pyruvate to lactate, which is expected to increase in aggressive cancers. The goal of this work was to develop and test analysis methods for improved quantification of this metabolic conversion. In this work, we compared specialized kinetic modeling methods to estimate the pyruvate-to-lactate conversion rate, k_PL , as well as the lactate-to-pyruvate area-under-curve (AUC) ratio. The kinetic modeling included an "inputless" method requiring no assumptions regarding the input function, as well as a method incorporating bolus characteristics in the fitting. These were first evaluated with simulated data designed to match human prostate data, where we examined the expected sensitivity of metabolism quantification to variations in k_PL , signal-to-noise ratio (SNR), bolus characteristics, relaxation rates, and B₁ variability. They were then applied to 17 prostate cancer patient datasets. The simulations indicated that the inputless method with fixed relaxation rates provided high expected accuracy with no sensitivity to bolus characteristics. The AUC ratio showed an undesired strong sensitivity to bolus variations. Fitting the input function as well did not improve accuracy over the inputless method. In vivo results showed qualitatively accurate k_PL maps with inputless fitting. The AUC ratio was sensitive to bolus delivery variations. Fitting with the input function showed high variability in parameter maps. Overall, we found the inputless k_PL fitting method to be a simple, robust approach for quantification of metabolic conversion following HP [1-¹³ C]pyruvate injection in human prostate cancer studies. This study also provided initial ranges of HP [1-¹³ C]pyruvate parameters (SNR, k_PL , bolus characteristics) in the human prostate.

Investigation of analysis methods for hyperpolarized 13C-pyruvate metabolic MRI in prostate cancer patients

MRI using hyperpolarized (HP) carbon-13 pyruvate is being investigated in clinical trials to provide non-invasive measurements of metabolism for cancer and cardiac imaging. In this project, we applied HP [1-13 C]pyruvate dynamic MRI in prostate cancer to measure the conversion from pyruvate to lactate, which is expected to increase in aggressive cancers. The goal of this work was to develop and test analysis methods for improved quantification of this metabolic conversion. In this work, we compared specialized kinetic modeling methods to estimate the pyruvate-to-lactate conversion rate, kPL , as well as the lactate-to-pyruvate area-under-curve (AUC) ratio. The kinetic modeling included an "inputless" method requiring no assumptions regarding the input function, as well as a method incorporating bolus characteristics in the fitting. These were first evaluated with simulated data designed to match human prostate data, where we examined the expected sensitivity of metabolism quantification to variations in kPL , signal-to-noise ratio (SNR), bolus characteristics, relaxation rates, and B1 variability. They were then applied to 17 prostate cancer patient datasets. The simulations indicated that the inputless method with fixed relaxation rates provided high expected accuracy with no sensitivity to bolus characteristics. The AUC ratio showed an undesired strong sensitivity to bolus variations. Fitting the input function as well did not improve accuracy over the inputless method. In vivo results showed qualitatively accurate kPL maps with inputless fitting. The AUC ratio was sensitive to bolus delivery variations. Fitting with the input function showed high variability in parameter maps. Overall, we found the inputless kPL fitting method to be a simple, robust approach for quantification of metabolic conversion following HP [1-13 C]pyruvate injection in human prostate cancer studies. This study also provided initial ranges of HP [1-13 C]pyruvate parameters (SNR, kPL , bolus characteristics) in the human prostate.

Assessing the influence of isoflurane anesthesia on cardiac metabolism using hyperpolarized [1-13 C]pyruvate

Isoflurane is a frequently used anesthetic in small-animal dissolution dynamic nuclear polarization-magnetic resonance imaging (DNP-MRI) studies. Although the literature suggests interactions with mitochondrial metabolism, the influence of the compound on cardiac metabolism has not been assessed systematically to date. In the present study, the impact of low versus high isoflurane concentration was examined in a crossover experiment in healthy rats. The results revealed that cardiac metabolism is modulated by isoflurane concentration, showing increased [1-¹³ C]lactate and reduced [¹³ C]bicarbonate production during high isoflurane relative to low isoflurane dose [average differences: +16% [1-¹³ C]lactate/total myocardial carbon, -22% [¹³ C]bicarbonate/total myocardial carbon; +51% [1-¹³ C]lactate/[¹³ C]bicarbonate]. These findings emphasize that reproducible anesthesia is important when studying cardiac metabolism. As the depth of anesthesia is difficult to control in an experimental animal setting, careful study design is required to exclude confounding factors.

Dioxygen and Metabolism; Dangerous Liaisons in Cardiac Function and Disease

The heart must consume a significant amount of energy to sustain its contractile activity. Although the fuel demands are huge, the stock remains very low. Thus, in order to supply its daily needs, the heart must have amazing adaptive abilities, which are dependent on dioxygen availability. However, in myriad cardiovascular diseases, “fuel” depletion and hypoxia are common features, leading cardiomyocytes to favor low-dioxygen-consuming glycolysis rather than oxidation of fatty acids. This metabolic switch makes it challenging to distinguish causes from consequences in cardiac pathologies. Finally, despite the progress achieved in the past few decades, medical treatments have not improved substantially, either. In such a situation, it seems clear that much remains to be learned about cardiac diseases. Therefore, in this review, we will discuss how reconciling dioxygen availability and cardiac metabolic adaptations may contribute to develop full and innovative strategies from bench to bedside.