Abstract
The goal of this study was to establish a quantitative method for measuring fatty acid (FA) metabolism with partial volume (PV) and spill-over (SP) corrections using dynamic [11C]palmitate positron emission tomographic (PET) images of mouse heart in vivo. Twenty-minute dynamic [11C]palmitate PET scans of four 18- to 20-week-old male C57BL/6 mice under isoflurane anesthesia were performed using a Focus F-120 PET scanner. A model-corrected blood input function, by which the input function with SP and PV corrections and the metabolic rate constants (k1–k5) are simultaneously estimated from the dynamic [11C]palmitate PET images of mouse hearts in a four-compartment tracer kinetic model, was used to determine rates of myocardial fatty acid oxidation (MFAO), myocardial FA esterification, myocardial FA use, and myocardial FA uptake. The MFAO thus measured in C57BL/6 mice was 375.03 ± 43.83 nmol/min/g. This compares well to the MFAO measured in perfused working C57BL/6 mouse hearts ex vivo of about 350 nmol/g/min and 400 nmol/min/g. FA metabolism was measured for the first time in mouse heart in vivo using dynamic [11C]palmitate PET in a four-compartment tracer kinetic model. MFAO obtained with this model was validated by results previously obtained with mouse hearts ex vivo.
HEART FAILURE (HF) affects nearly 6 million people in the United States. Coronary heart disease, high blood pressure, and diabetes are common causes of HF. The estimated cost for HF management in the United States is about 34 billion dollars per year. 1 Current medical therapies are ineffective, with a 5-year survival rate of only 25%. 2 It is well established that cardiac metabolism is inextricably linked to cardiac function, and modulation of substrate energy metabolism may be an attractive strategy for treatment of HF. 3 5 However, the molecular mechanisms connecting cardiac function and metabolism are not well understood, and our current understanding of metabolic alterations in the heart is largely based on ex vivo studies in animal models. 6 9
Transgenic and knockout mouse models offer opportunities to address this issue in vivo and in parallel at the molecular level. However, a severe limitation to this approach is that there are no suitable methods that allow noninvasive assessment of cardiac energy metabolism and correlation to cardiac function in an intact mouse.
Imaging myocardial metabolism in mouse heart is challenging due to the limited intrinsic resolution of small animal scanners and the small size of the mouse heart, resulting in image blur. Recently, we optimized a model-corrected blood input function (MCBIF) from time-resolved gated 2-[18F]fluoro-2-deoxy-D-glucose (FDG) positron emission tomographic (PET) images of mouse heart to measure myocardial glucose metabolism in vivo.10,11 Measuring fatty acid (FA) metabolism using [11C]palmitate in mouse heart is even more challenging than measuring glucose metabolism due to the long range of positrons emitted by 11C. Poor quality of images results in severe partial volume (PV) averaging, and spill-over (SP) contamination is further confounded by the rapid heart rate. Previous studies using [11C] palmitate PET have only been performed in dog and rat hearts. 12 15 In this study, we developed a method to measure FA metabolism from dynamic [11C]palmitate PET images of mouse heart using a four-compartment tracer kinetic model, for the first time, in vivo. Our computed rate of myocardial fatty acid oxidation (MFAO) in mouse heart in vivo was in agreement with the measured palmitate oxidation in the working mouse heart ex vivo. 16
Materials and Methods
Animal Model
Four adult C57BL/6N male mice (18–20 weeks of age) were imaged using the Focus-F120 PET scanner (Siemens Medical Solutions USA, Inc., Malvern, PA) under isoflurane anesthesia. 17 The mice were kept on a normal chow diet (7912 Teklad LM-485 from Harlan Laboratories, Indianapolis, IN) and had free access to food and water before the study. All experiments were performed in compliance with the Guide for the Care and Use of Laboratory Animals, published by the National Institutes of Health, and were conducted under protocols approved by the Institutional Animal Care and Use Committee at the University of Virginia.
[11C]Palmitate Labeling
[11C]Palmitate radiolabeling was performed in the Radiochemistry Core Laboratory at the University of Virginia Molecular Imaging Center using the Modular-Lab Pharm-Tracer system (Eckert & Ziegler, E&Z, Hopkinton, MA). We modified the E&Z cassette and software designed for the production of [11C]acetate to make [11C]palmitate following the method of Mock and colleagues. 18 Briefly, [11C]CO2 from the cyclotron was trapped by carbospheres at −20oC to remove excess helium carrier gas. The carbospheres were heated to 60oC, at which point, helium was passed over the carbospheres to force the [11C]CO2 through the precursor solutions, which consisted of 0.143 M pentadecylmagnesium bromide in diethyl ether (NOVEL Chemical Solutions, Crete, NE). The reaction vessel was then sealed and heated to 30°C for 3 minutes to complete the reaction. The reaction was quenched with 1 normal HCl diethyl ether. The quenched solution was passed through a dry Alumina-N Sep-Pak Plus cartridge that trapped the polar [11C]palmitic acid. Nonpolar chemicals were removed from the Sep-Pak by rinsing with three 1 mL diethyl ether washes. The Sep-Pak was rinsed with two 1 mL portions of 0.5 M NaH2PO4 solution to disrupt the carboxylic acid binding to the alumina surface, followed by two 1 mL water washes to remove soluble Mg2+ salts from the Sep-Pak. The [11C]palmitic acid was eluted from the Sep-Pak with three portions of absolute ethanol (2 mL total volume), and the solution was then passed through a sterile 0.2 mm Millex-FG filter (EMD Millipore, Darmstadt, Germany) into a sterile intermediate product vial. The non-decay-corrected radiochemical yield of the [11C]palmitic acid product averaged ~ 30% (250 mCi of [11C]CO2 yielded 75 mCi of [11C]palmitate). The [11C]palmitate ethanol solution was formulated with bovine serum albumin. The synthesis was completed within about 30 minutes.
PET Imaging
Using the [11C]palmitate label described above, we imaged FA metabolism in mouse heart in vivo as follows. Briefly, a mouse with electrocardiographic surface electrodes and a respiratory pillow attached to its limbs and chest, respectively, and under 1 to 1.5% isoflurane anesthesia was positioned in the bore of the PET scanner. 19 A 10-minute 57Co transmission scan was performed for attenuation correction. A 20-minute dynamic PET acquisition was then initiated prior to the administration of 300 to 350 μCi in about 200 to 300 μL of [11C]palmitate via the tail vein over 30 seconds. The list-mode dynamic PET data were reconstructed using ordered subset expectation maximization maximum a posteriori (OSEM-MAP) algorithm 20 with attenuation correction into the following dynamic frames (frames, time [s]: 12, 5; 8, 30; 4, 150; 1, 300). The reconstructed images were composed of 95 transverse slices with a thickness of 0.79 mm and an in-plane voxel resolution of 0.4 × 0.4 mm (128 × 128 pixels) corresponding to a zoom factor of 2.13. Regions of interest (ROI) in the region corresponding to the left ventricular blood pool (LVBP) and the myocardium were drawn in the last time frame of the dynamic image data in the transverse plane and time-activity curves (TACs) generated for the LVBP and the myocardium for the 20-minute scan. The blood and myocardium TACs were used in a four-compartment tracer kinetic model (described below), written in MATLAB (The Mathworks, Natick, MA) using nonlinear regression with SP and PV corrections, to compute myocardial FA metabolism, including FA oxidation, esterification, use, and uptake in mouse heart in vivo.
Four-Compartment Tracer Kinetic Model
Following Bergmann and colleagues, 15 the differential equations for [11C]palmitate kinetics can be written as follows:
where q1–q4 are the concentrations in the four compartments as shown in Figure 1, kn represents the turnover rate constants between the compartments, F is myocardial blood flow, V is the vascular volume, and Ca(t) is the arterial concentration of tracer over time.

Block diagram of the four-compartment model. Compartment 1 represents the vascular space; 2, the interstitial and intracellular spaces; 3, neutral lipids, amino acids, and other slow turnover pools; and 4, mitochondrial β-oxidation. kn represents the forward and backward rate constants between compartments, F is myocardial blood flow, V is fractional vascular volume, and Ca(t) is arterial tracer concentration over time.
The total tracer concentration in myocardium can be defined by the sum of the tracer in each compartment:
Performing Laplace and inverse Laplace transforms (see Supplementary Material, online version only), we get the final form as follows:
where M is composed of the rate constants (k1–k5) and F/V.
Dual Output Model
Ideally, when an ROI is drawn within the cavity of the left ventricle, the image-derived input function (IDIF) would equal the whole-blood TAC Ca(t). However, due to SP and PV effects, the model equation for an image-derived TAC from the blood pool can be written as fractions of the tissue concentration in the blood compartment and partial recovery of radio activity concentration from the blood as follows:
Similarly, the myocardium tissue of the model equation is
where rb and rm are the recovery coefficients (accounting for PV effect) for the myocardium and blood pool, respectively. Smb and Sbm are the SP coefficients from the blood pool to the myocardium and vice versa, respectively; tei and tbi are the beginning and end times, respectively, for frame in a dynamic PET scan.
The model equation for the blood input function can be written as follows:
where each term determines the amplitude, shape, and washout of the tracer over time; r indicates the time lag between the initiation of imaging and injection of the tracer. The model equations can be fitted to the blood (PETIDIF) and myocardial tissue (PETmyo) TACs obtained from OSEM- MAP dynamic PET images with attenuation correction by substituting equations 6 and 9 into equations 7 and 8, as indicated below:
Interpolation and Simultaneous Estimation
The shorter reconstructed frames at the early time points limited the data resolution, resulting in fluctuations in the blood and tissue TACs at these points. A cubic spline interpolation was used to approximate the peak of curve and improve the accuracy of fitting. The interpolation was used near the peak (PETIDIF,max, PETmyo,max) of blood (PETIDIF,i) and tissue (PETmyo,i) concentrations obtained from the dynamic PET data. The interpolated sequences were written as follows:
All numeric analyses were done using MATLAB R2013b. Optimization was performed by minimizing the objective function (equation 10) using the MATLAB function “fmincon,” which is based on an interior-reflective Newton method. The initial guesses and bounds for all the parameters used in the optimization routine are shown in Table 1.
The following constraints and previous knowledge were used to set the bounds: (a) rates of palmitate oxidation are greater than esterification, that is, k5 > k3 under control conditions 15 ; (b) in the steady state, the rate of backflow from compartment 3 (neutral lipids and amino acids) to compartment 2 (interstitial and cytosolic) is close to zero, that is, k4 = 0 15 ; and (c) the bounds for the recovery coefficients were based on structural measurements using magnetic resonance imaging (MRI) of control mouse heart. 21 Since the spillover factors, Sbm and Smb, are changing with time (multiplied by time-dependent image-derived blood and tissue curves; see equations 7 and 8), there is no good way to get a previous estimate of these factors. Hence, we set the bounds for these factors from 0 to 1. We also observed that the computed results have little influence on the changes in the bounds and initial guess values of the seven-parameter blood function (see equation 9).
The minimization of equation 10 results in simultaneous estimation of the parameters of blood input function Ca(t), compartment model parameters k1–k5, and the SP and PV coefficients (Smb, rb, rm, Sbm). The optimized rate constants and SP and PV coefficients are listed in Table 2. A myocardial blood flow (F) value of 4.8 mL/min/g 22 and vascular blood volume (V) of 0.1 mL/g 15 was used for the calculations.
Free Fatty Acid in Mouse Blood
For the measurement of free fatty acid (FFA) levels in the mouse blood, samples (~ 20 μL) from the tail vein were collected during the dynamic PET scan at 10 and 15 minutes following tracer administration using nonheparinized capillary tubes. The blood was allowed to clot for 30 minutes; serum was then separated from the clot by centrifugation (15 minutes at 1,600g at room temperature in a microcentrifuge) and stored in a −20oC freezer until analyzed. FFA levels were measured using HR series NEFA-HR2 reagents (Wako Life Sciences Inc., Richmond, VA).
Results
In this study, using dynamic [11C]palmitate PET of mouse heart in vivo, we measured myocardial FA metabolism including rates of oxidation, esterification, use, and uptake.
Biodistribution of [99mTc]-Met Peptide in MKN-45 Xenografts after 30, 60, 120, and 180 Minutes
Each value represents (% injected dose per gram [%ID/g]) ± SD of [99mTc]-Met peptide (n = 4 [60-minute group] or 5 [all other time points]).
Initial Guess Values and Bounds for Parameters Used in the Optimization
Figure 2A shows an example of a [11C]palmitate PET image of mouse heart in vivo. ROI are shown in the blood and the myocardium in the last time frame. The ROI drawn in the last time frame was used to obtain TACs for the blood and the myocardium from the dynamic PET data. In Figure 2B, the model fits (see equations 710) to the blood and myocardium TACs are shown without interpolation resulting in poor fit to the blood TAC around the peak region. Cubic spline interpolation with a time interval of 0.5 seconds around the peak resulted in better fits to both the blood and myocardium TACs (Figure 2C). Model fits resulted in simultaneous estimation of blood input function Ca(t), compartment model parameters k1–k5, and the SP and PV coefficients (Smb, rb, rm, Sbm). The MCBIF, Ca(t), was compared to the image-derived blood input function (PET blood) in Figure 3. The recovery coefficients for the blood pool and the myocardium are denoted by rb and rm, respectively. Since they are a function of structure, the bounds used in the optimization routine were based on structural measurements using MRI of control mouse heart. 21 As indicated in Figure 3, the model correction results in an ~ 20% increase in the MCBIF at the early time points compared to image-derived input function (PET blood), indicating improved radioactivity recovery. Model correction also eliminated SP contamination from the myocardium into the blood pool at the late time points.

[11C]Palmitate PET images of mouse heart in vivo. (A) Example of a PET image at the last time frame of the dynamic data set. The image also exhibits regions of interest in the left ventricular blood pool and the myocardium. (B) Model fits to the blood and the myocardial time-activity curves without interpolation. (C) Cubic spline interpolation resulted in improved fits, especially at the peak region for the time-activity curve obtained from the left ventricular blood pool.

The model-corrected blood input function (MCIF), Ca(t), computed by simultaneous estimation compared to the image-derived blood input function (IDIF: PET blood) obtained from the dynamic [11C]palmitate PET images. MCIF estimation improved radioactivity recovery at the early time points and eliminated SP contamination at the late time points from the myocardium to the left ventricular blood pool.
Using the formulas described in detail in the Supplementary Material (online version only), under steady-state condition, the rates of MFAO, esterification (MFAE), and hence use (MFAU) were computed. Figure 4 shows the computed average MFAO (375.03 ± 43.83 nmol/g/min), MFAE (10.97 ± 7.11 nmol/g/min), and MFAU (385.99 ± 49.51 nmol/g/min) obtained by optimization of the dynamic [11C]palmitate PET data from four 18- to 20-week-old male C57BL/6N mice. The above computed MFAO in vivo agrees with the range of MFAO measured ex vivo in perfused working 12- and 15-week-old strain-matched mouse heart. 16

Fatty acid metabolic parameters in mouse heart in vivo. Myocardial fatty acid oxidation (MFAO), myocardial fatty acid esterification (MFAE), and myocardial fatty acid use (MFAU) computed from dynamic [11C]palmitate PET images of mouse heart in vivo.
Myocardial fatty acid uptake rate (MFAUp) was also computed by dividing MFAU with measured FFA (849 nM/mL) in mouse blood. The computed MFAUp rate was 0.45 ± 0.06 mL/min/g. When compared to the rate of myocardial FDG uptake, Ki (0.15 ± 0.03 mL/min/g) measured in age- and strain-matched mice in a previously published study from our laboratory, 10 we observed a threefold higher uptake rate of palmitate (MFAUp) compared to FDG uptake rate, Ki.
Discussion
The commonly used tracers for measuring myocardial FA metabolism are 14(R,S)-[18F]fluoro-6-thia-heptadecanoic acid (FTHA) 23 and [11C]palmitate.12,24 There is a distinct advantage in using [11C]palmitate due to the short half-life of the labeled tracer allowing serial measurements of FA and glucose metabolism in the same animal and monitoring the effect oftherapy in vivo. Imaging the mouse heart using a small animal PET scanner is a challenge due to the limited intrinsic resolution, 19 caused primarily by the finite positron range, which results in image blur. This leads to incomplete radioactivity recovery (PV averaging) and hence SP of radioactivity from the myocardium to the blood pool and vice versa. Cardiac motion results in further image blur, confounding the effects of SP contamination in the image-derived blood input function. 25 In a recent study, we optimized an MCBIF from time-resolved gated FDG PET images of mouse heart in vivo. 10 The MCBIF was then used to evaluate glucose metabolism in the pressure overload mouse heart in vivo. 11 We observed that under transverse aortic constriction–induced pressure overload stress, the mouse heart exhibited a sixfold increase in FDG uptake at day 1 postconstriction with preserved contractile function, indicating an adaptive metabolic response. At day 7 postconstriction, a further increase in myocardial FDG uptake was found, together with a decrease in cardiac function, indicating maladaptation. 11 However, how changes in FA metabolism relate to cardiac dysfunction in the stressed mouse heart is not known.
Clinical data using [11C]palmitate PET in hypertensive left ventricular hypertrophy (LVH) suggested a decrease in myocardial FA oxidation and use as possible causes for a decrease in myocardial efficiency. 26 Recent 18F-FTHA PET imaging studies in vivo in the spontaneously hypertensive rat (SHR) heart over a period of 20 months observed a significant increase in both glucose and fatty acid use, which appeared to precede mechanical changes in LVH progression. 27 Studies in a Dahl salt-sensitive rat model ex vivo, however, showed the decrease in FA metabolism to occur late in the progression from LVH to heart failure, whereas glucose uptake was increased at an earlier stage. 28 These studies indicate that there is no clear consensus on the abnormalities in FA metabolism in the stressed heart. Improved tools for quantitatively imaging FA metabolism in a compartment model are needed to complement the conventional graphical Patlak analysis. 27
In this study, we optimized and implemented a four-compartment tracer kinetic model with SP and PV corrections to quantify FA metabolism from dynamic [11C]palmitate images of mouse heart in vivo. We measured rates of MFAO, MFAE, MFAU, and MFAUp for the first time in mouse heart in vivo. This method can be easily implemented to evaluate the role of FA metabolism in the pressure overload mouse heart in vivo.
Our study, however, is not without limitations. The blood input curves obtained from dynamic PET images were not corrected for the conversion of palmitate to [11C]CO2. But studies in dog hearts have shown that this correction performed by arterial blood sampling only modestly affects the input curve, especially in the tail region. 15 The factors that may be affected as a result are k3 (rate of transfer to triglycerides) and the fractional esterification rates because they are represented by the tail of the TAC. Our data in mouse heart, however, suggest that the tail of the input curve may be minimally affected by the correction due to the conversion of palmitate to [11C]CO2. Thus, this correction may not be necessary when imaging a mouse heart using [11C]palmitate, making the process of obtaining the FA metabolic parameters completely noninvasive.
Conclusions
We have measured FA metabolism, including myocardial FA oxidation, esterification, use, and uptake, for the first time from dynamic [11C]palmitate PET images of mouse heart in vivo. The results for MFAO measured in vivo compared well to the measured range of palmitate oxidation from perfused working mouse heart ex vivo. By comparing rates of FA uptake measured in this study with rates of FDG uptake determined in an earlier study, we corroborated the theory that the normal mouse heart preferentially uses FA over glucose.
Footnotes
Acknowledgments
We thank Gina Wimer and Jeremy Gatesman for performing the tail-vein injections during the course of the study.
Financial disclosure of authors: This work was supported in part by grants R21 HL-102627 (to B.K.K.) and R01 DK-81471 (to S.R.K.) from the National Institutes of Health.
Financial disclosure of reviewers: None reported.
