Triple repetition time saturation transfer (TRiST) 31P spectroscopy for measuring human creatine kinase reaction kinetics

Fast in vivo assay of creatine kinase activity in the human brain by 31 P magnetic resonance fingerprinting

A new and efficient magnetisation transfer 31 P magnetic resonance fingerprinting (MT-31 P-MRF) approach is introduced to measure the creatine kinase metabolic ratek CK between phosphocreatine (PCr) and adenosine triphosphate (ATP) in human brain. The MRF framework is extended to overcome challenges in conventional 31 P measurement methods in the human brain, enabling reduced acquisition time and specific absorption rate (SAR). To address the challenge of creating and matching large multiparametric dictionaries in an MRF scheme, a nested iteration interpolation method (NIIM) is introduced. As the number of parameters to estimate increases, the size of the dictionary grows exponentially. NIIM can reduce the computational load by breaking dictionary matching into subsolutions of linear computational order. MT-31 P-MRF combined with NIIM providesT 1 PCr ,T 1 ATP andk CK estimates in good agreement with those obtained by the exchange kinetics by band inversion transfer (EBIT) method and literature values. Furthermore, the test-retest reproducibility results showed that MT-31 P-MRF achieves a similar or better coefficient of variation (<12%) forT 1 ATP andk CK measurements in 4 min 15 s, than EBIT with 17 min 4 s scan time, enabling a fourfold reduction in scan time. We conclude that MT-31 P-MRF in combination with NIIM is a fast, accurate, and reproducible approach for in vivok CK assays in the human brain, which enables the potential to investigate energy metabolism in a clinical setting.

Cardiac 31P MR spectroscopy: development of the past five decades and future vision-will it be of diagnostic use in clinics?

In the past five decades, the use of the magnetic resonance (MR) technique for cardiovascular diseases has engendered much attention and raised the opportunity that the technique could be useful for clinical applications. MR has two arrows in its quiver: One is magnetic resonance imaging (MRI), and the other is magnetic resonance spectroscopy (MRS). Non-invasively, highly advanced MRI provides unique and profound information about the anatomical changes of the heart. Excellently developed MRS provides irreplaceable and insightful evidence of the real-time biochemistry of cardiac metabolism of underpinning diseases. Compared to MRI, which has already been successfully applied in routine clinical practice, MRS still has a long way to travel to be incorporated into routine diagnostics. Considering the exceptional potential of ³¹P MRS to measure the real-time metabolic changes of energetic molecules qualitatively and quantitatively, how far its powerful technique should be waited before a successful transition from "bench-to-bedside" is enticing. The present review highlights the seminal studies on the chronological development of cardiac ³¹P MRS in the past five decades and the future vision and challenges to incorporating it for routine diagnostics of cardiovascular disease.

Three-dimensional, 2.5-minute, 7T phosphorus magnetic resonance spectroscopic imaging of the human heart using concentric rings

A three-dimensional (3D), density-weighted, concentric rings trajectory (CRT) magnetic resonance spectroscopic imaging (MRSI) sequence is implemented for cardiac phosphorus (³¹ P)-MRS at 7 T. The point-by-point k-space sampling of traditional phase-encoded chemical shift imaging (CSI) sequences severely restricts the minimum scan time at higher spatial resolutions. Our proposed CRT sequence implements a stack of concentric rings, with a variable number of rings and planes spaced to optimise the density of k-space weighting. This creates flexibility in acquisition time, allowing acquisitions substantially faster than traditional phase-encoded CSI sequences, while retaining high signal-to-noise ratio (SNR). We first characterise the SNR and point-spread function of the CRT sequence in phantoms. We then evaluate it at five different acquisition times and spatial resolutions in the hearts of five healthy participants at 7 T. These different sequence durations are compared with existing published 3D acquisition-weighted CSI sequences with matched acquisition times and spatial resolutions. To minimise the effect of noise on the short acquisitions, low-rank denoising of the spatiotemporal data was also performed after acquisition. The proposed sequence measures 3D localised phosphocreatine to adenosine triphosphate (PCr/ATP) ratios of the human myocardium in 2.5 min, 2.6 times faster than the minimum scan time for acquisition-weighted phase-encoded CSI. Alternatively, in the same scan time, a 1.7-times smaller nominal voxel volume can be achieved. Low-rank denoising reduced the variance of measured PCr/ATP ratios by 11% across all protocols. The faster acquisitions permitted by 7-T CRT ³¹ P-MRSI could make cardiac stress protocols or creatine kinase rate measurements (which involve repeated scans) more tolerable for patients without sacrificing spatial resolution.

Magnetic resonance imaging of cardiac metabolism in heart failure: how far have we come?

As one of the highest energy consumer organs in the body, the heart requires tremendous amount of adenosine triphosphate (ATP) to maintain its continuous mechanical work. Fatty acids, glucose, and ketone bodies are the primary fuel source of the heart to generate ATP with perturbations in ATP generation possibly leading to contractile dysfunction. Cardiac metabolic imaging with magnetic resonance imaging (MRI) plays a crucial role in understanding the dynamic metabolic changes occurring in the failing heart, where the cardiac metabolism is deranged. Also, targeting and quantifying metabolic changes in vivo noninvasively is a promising approach to facilitate diagnosis, determine prognosis, and evaluate therapeutic response. Here, we summarize novel MRI techniques used for detailed investigation of cardiac metabolism in heart failure including magnetic resonance spectroscopy (MRS), hyperpolarized MRS, and chemical exchange saturation transfer based on evidence from preclinical and clinical studies and to discuss the potential clinical application in heart failure.

Myocardial ATP depletion detected noninvasively predicts sudden cardiac death risk in heart failure patients

BACKGROUND

Sudden cardiac death (SCD) remains a worldwide public health problem in need of better noninvasive predictive tools. Current guidelines for primary preventive SCD therapies, such as implantable cardioverter defibrillators (ICDs), are based on left ventricular ejection fraction (LVEF), but these guidelines are imprecise: fewer than 5% of ICDs deliver lifesaving therapy per year. Impaired cardiac metabolism and ATP depletion cause arrhythmias in experimental models, but to our knowledge a link between arrhythmias and cardiac energetic abnormalities in people has not been explored, nor has the potential for metabolically predicting clinical SCD risk.

METHODS

We prospectively measured myocardial energy metabolism noninvasively with phosphorus magnetic resonance spectroscopy in patients with no history of significant arrhythmias prior to scheduled ICD implantation for primary prevention in the setting of reduced LVEF (≤35%).

RESULTS

By 2 different analyses, low myocardial ATP significantly predicted the composite of subsequent appropriate ICD firings for life-threatening arrhythmias and cardiac death over approximately 10 years. Life-threatening arrhythmia risk was approximately 3-fold higher in patients with low ATP and independent of established risk factors, including LVEF. In patients with normal ATP, rates of appropriate ICD firings were several-fold lower than reported rates of ICD complications and inappropriate firings.

CONCLUSION

To the best of our knowledge, these are the first data linking in vivo myocardial ATP depletion and subsequent significant arrhythmic events in people, suggesting an energetic component to clinical life-threatening ventricular arrhythmogenesis. The findings support investigation of metabolic strategies that limit ATP loss to treat or prevent life-threatening cardiac arrhythmias and herald noninvasive metabolic imaging as a complementary SCD risk stratification tool.

TRIAL REGISTRATION

ClinicalTrials.gov NCT00181233.

FUNDING

This work was supported by the DW Reynolds Foundation, the NIH (grants HL61912, HL056882, HL103812, HL132181, HL140034), and Russell H. Morgan and Clarence Doodeman endowments at Johns Hopkins.

Mitochondrial Creatine Kinase Attenuates Pathologic Remodeling in Heart Failure

BACKGROUND

Abnormalities in cardiac energy metabolism occur in heart failure (HF) and contribute to contractile dysfunction, but their role, if any, in HF-related pathologic remodeling is much less established. CK (creatine kinase), the primary muscle energy reserve reaction which rapidly provides ATP at the myofibrils and regenerates mitochondrial ADP, is down-regulated in experimental and human HF. We tested the hypotheses that pathologic remodeling in human HF is related to impaired cardiac CK energy metabolism and that rescuing CK attenuates maladaptive hypertrophy in experimental HF.

METHODS

First, in 27 HF patients and 14 healthy subjects, we measured cardiac energetics and left ventricular remodeling using noninvasive magnetic resonance 31P spectroscopy and magnetic resonance imaging, respectively. Second, we tested the impact of metabolic rescue with cardiac-specific overexpression of either Ckmyofib (myofibrillar CK) or Ckmito (mitochondrial CK) on HF-related maladaptive hypertrophy in mice.

RESULTS

In people, pathologic left ventricular hypertrophy and dilatation correlate closely with reduced myocardial ATP levels and rates of ATP synthesis through CK. In mice, transverse aortic constriction-induced left ventricular hypertrophy and dilatation are attenuated by overexpression of CKmito, but not by overexpression of CKmyofib. CKmito overexpression also attenuates hypertrophy after chronic isoproterenol stimulation. CKmito lowers mitochondrial reactive oxygen species, tissue reactive oxygen species levels, and upregulates antioxidants and their promoters. When the CK capacity of CKmito-overexpressing mice is limited by creatine substrate depletion, the protection against pathologic remodeling is lost, suggesting the ADP regenerating capacity of the CKmito reaction rather than CK protein per se is critical in limiting adverse HF remodeling.

CONCLUSIONS

In the failing human heart, pathologic hypertrophy and adverse remodeling are closely related to deficits in ATP levels and in the CK energy reserve reaction. CKmito, sitting at the intersection of cardiac energetics and redox balance, plays a crucial role in attenuating pathologic remodeling in HF. Registration: URL: https://www.clinicaltrials.gov; Unique identifier: NCT00181259.

Myocardial Energy Response to Glyceryl Trinitrate: Physiology Revisited

Objective: Although intravenous nitrates are commonly used in clinical medicine, they have been shown to increase myocardial oxygen consumption and inhibit complex IV of the electron transport chain. As such we sought to measure whether myocardial energetics were impaired during glyceryl trinitrate (GTN) infusion. Methods: 10 healthy volunteers underwent cardiac magnetic resonance imaging to assess cardiac function and 31phosphorus magnetic resonance spectroscopy to measure Phosphocreatine/ATP (PCr/ATP) ratio and creatine kinase forward rate constant (CK kf ) before and during an intravenous infusion of GTN. Results: During GTN infusion, mean arterial pressure (78 ± 7 vs. 65 ± 6 mmHg, p < 0.001), left ventricular (LV) stroke work (7,708 ± 2,782 vs. 6,071 ± 2,660 ml mmHg, p < 0.001), and rate pressure product (7,214 ± 1,051 vs. 6,929 ± 976 mmHg bpm, p = 0.06) all fell. LV ejection fraction increased (61 ± 3 vs. 66 ± 4%, p < 0.001), with cardiac output remaining constant (6.2 ± 1.5 vs. 6.5 ± 1.4 l/min, p = 0.37). Myocardial PCr/ATP fell during GTN infusion (2.17 ± 0.2 vs. 1.99 ± 0.22, p = 0.03) with an increase in both CK kf (0.16 ± 0.07 vs. 0.25 ± 0.1 s-1, p = 0.006) and CK flux (1.8 ± 0.8 vs. 2.6 ± 1.1 μmol/g/s, p = 0.03). Conclusion: During GTN infusion, despite reduced LV stroke work and maintained cardiac output, there was a 44% increase in myocardial ATP delivery through CK. As PCr/ATP fell, this increase in ATP demand coincided with GTN-induced impairment of mitochondrial oxidative phosphorylation. Overall, this suggests that while GTN reduces cardiac work, it does so at the expense of increasing ATP demand beyond the capacity to increase ATP production.

Obesity modifies the energetic phenotype of dilated cardiomyopathy

AIMS

We sought to determine if myocardial energetics could distinguish obesity cardiomyopathy as a distinct entity from dilated cardiomyopathy.

METHODS AND RESULTS

Sixteen normal weight participants with dilated cardiomyopathy (DCMNW), and 27 with DCM and obesity (DCMOB), were compared to 26 normal weight controls (CTLNW). All underwent cardiac magnetic resonance imaging and 31P spectroscopy to assess function and energetics. Nineteen DCMOB underwent repeat assessment after a dietary weight loss intervention. Adenosine triphosphate (ATP) delivery through creatine kinase (CK flux) was 55% lower in DCMNW than in CTLNW (P = 0.004), correlating with left ventricular ejection fraction (LVEF, r = 0.4, P = 0.015). In contrast, despite similar LVEF (DCMOB 41 ± 7%, DCMNW 38 ± 6%, P = 0.14), CK flux was two-fold higher in DCMOB (P < 0.001), due to higher rate through CK [median kf 0.21 (0.14) vs. 0.11 (0.12) s-1, P = 0.002]. During increased workload, the CTLNW heart increased CK flux by 97% (P < 0.001). In contrast, CK flux was unchanged in DCMNW and fell in DCMOB (by >50%, P < 0.001). Intentional weight loss was associated with positive left ventricular remodelling, with reduced left ventricular end-diastolic volume (by 8%, P < 0.001) and a change in LVEF (40 ± 9% vs. 45 ± 10%, P = 0.002). This occurred alongside a fall in ATP delivery rate with weight loss (by 7%, P = 0.049).

CONCLUSIONS

In normal weight, DCM is associated with reduced resting ATP delivery. In obese DCM, ATP demand through CK is greater, suggesting reduced efficiency of energy utilization. Dietary weight loss is associated with significant improvement in myocardial contractility, and a fall in ATP delivery, suggesting improved metabolic efficiency. This highlights distinct energetic pathways in obesity cardiomyopathy, which are both different from dilated cardiomyopathy, and may be reversible with weight loss.

Effect of empagliflozin on ectopic fat stores and myocardial energetics in type 2 diabetes: the EMPACEF study

BACKGROUND

Empagliflozin is a sodium-glucose cotransporter 2 (SGLT2) inhibitor that has demonstrated cardiovascular and renal protection in patients with type 2 diabetes (T2D). We hypothesized that empaglifozin (EMPA) could modulate ectopic fat stores and myocardial energetics in high-fat-high-sucrose (HFHS) diet mice and in type 2 diabetics (T2D).

METHODS

C57BL/6 HFHS mice (n = 24) and T2D subjects (n = 56) were randomly assigned to 12 weeks of treatment with EMPA (30 mg/kg in mice, 10 mg/day in humans) or with placebo. A 4.7 T or 3 T MRI with 1H-MRS evaluation-myocardial fat (primary endpoint) and liver fat content (LFC)-were performed at baseline and at 12 weeks. In humans, standard cardiac MRI was coupled with myocardial energetics (PCr/ATP) measured with 31P-MRS. Subcutaneous (SAT) abdominal, visceral (VAT), epicardial and pancreatic fat were also evaluated. The primary efficacy endpoint was the change in epicardial fat volume between EMPA and placebo from baseline to 12 weeks. Secondary endpoints were the differences in PCr/ATP ratio, myocardial, liver and pancreatic fat content, SAT and VAT between groups at 12 weeks.

RESULTS

In mice fed HFHS, EMPA significantly improved glucose tolerance and increased blood ketone bodies (KB) and β-hydroxybutyrate levels (p < 0.05) compared to placebo. Mice fed HFHS had increased myocardial and liver fat content compared to standard diet mice. EMPA significantly attenuated liver fat content by 55%, (p < 0.001) but had no effect on myocardial fat. In the human study, all the 56 patients had normal LV function with mean LVEF = 63.4 ± 7.9%. Compared to placebo, T2D patients treated with EMPA significantly lost weight (- 2.6 kg [- 1.2; - 3.7]) and improved their HbA1c by 0.88 ± 0.74%. Hematocrit and EPO levels were significantly increased in the EMPA group compared to placebo (p < 0.0001, p = 0.041). EMPA significantly increased glycosuria and plasma KB levels compared to placebo (p < 0.0001, p = 0.012, respectively), and significantly reduced liver fat content (- 27 ± 23 vs. - 2 ± 24%, p = 0.0005) and visceral fat (- 7.8% [- 15.3; - 5.6] vs. - 0.1% [- 1.1;6.5], p = 0.043), but had no effect on myocardial or epicardial fat. At 12 weeks, no significant change was observed in the myocardial PCr/ATP (p = 0.57 between groups).

CONCLUSIONS

EMPA effectively reduced liver fat in mice and humans without changing epicardial, myocardial fat or myocardial energetics, rebutting the thrifty substrate hypothesis for cardiovascular protection of SGLT2 inhibitors. Trial registration NCT, NCT03118336. Registered 18 April 2017, https://clinicaltrials.gov/ct2/show/NCT03118336.

Rapid, B 1 -insensitive, dual-band quasi-adiabatic saturation transfer with optimal control for complete quantification of myocardial ATP flux

PURPOSE

Phosphorus saturation-transfer experiments can quantify metabolic fluxes noninvasively. Typically, the forward flux through the creatine kinase reaction is investigated by observing the decrease in phosphocreatine (PCr) after saturation of γ-ATP. The quantification of total ATP utilization is currently underexplored, as it requires simultaneous saturation of inorganic phosphate ( P i ) and PCr. This is challenging, as currently available saturation pulses reduce the already-low γ-ATP signal present.

METHODS

Using a hybrid optimal-control and Shinnar-Le Roux method, a quasi-adiabatic RF pulse was designed for the dual saturation of PCr and P i to enable determination of total ATP utilization. The pulses were evaluated in Bloch equation simulations, compared with a conventional hard-cosine DANTE saturation sequence, before being applied to perfused rat hearts at 11.7 T.

RESULTS

The quasi-adiabatic pulse was insensitive to a >2.5-fold variation in B 1 , producing equivalent saturation with a 53% reduction in delivered pulse power and a 33-fold reduction in spillover at the minimum effective B 1 . This enabled the complete quantification of the synthesis and degradation fluxes for ATP in 30-45 minutes in the perfused rat heart. While the net synthesis flux (4.24 ± 0.8 mM/s, SEM) was not significantly different from degradation flux (6.88 ± 2 mM/s, P = .06) and both measures are consistent with prior work, nonlinear error analysis highlights uncertainties in the P_i -to-ATP measurement that may explain a trend suggesting a possible imbalance.

CONCLUSIONS

This work demonstrates a novel quasi-adiabatic dual-saturation RF pulse with significantly improved performance that can be used to measure ATP turnover in the heart in vivo.

Measuring Myocardial Energetics with Cardiovascular Magnetic Resonance Spectroscopy

The heart has the highest energy demands per gram of any organ in the body and energy metabolism fuels normal contractile function. Metabolic inflexibility and impairment of myocardial energetics occur with several common cardiac diseases, including ischemia and heart failure. This review explores several decades of innovation in cardiac magnetic resonance spectroscopy modalities and their use to noninvasively identify and quantify metabolic derangements in the normal, failing, and diseased heart. The implications of this noninvasive modality for predicting significant clinical outcomes and guiding future investigation and therapies to improve patient care are discussed.

Exercise intolerance and rapid skeletal muscle energetic decline in human age-associated frailty

BACKGROUND

Physical frailty in older individuals is characterized by subjective symptoms of fatigue and exercise intolerance (EI). Objective abnormalities in skeletal muscle (SM) mitochondrial high-energy phosphate (HEP) metabolism contribute to EI in inherited myopathies; however, their presence or link to EI in the frail older adult is unknown.

METHODS

Here, we studied 3 groups of ambulatory, community-dwelling adults with no history of significant coronary disease: frail older (FO) individuals (81 ± 2.7 years, mean ± SEM), nonfrail older (NFO) individuals (79 ± 2.0 years), and healthy middle-aged individuals, who served as controls (CONT, 51 ± 2.1 years). Lower extremity SM HEP levels and mitochondrial function were measured with ³¹P magnetic resonance (MR) techniques during graded multistage plantar flexion exercise (PFE). EI was quantified by a 6-minute walk (6MW) and peak oxygen consumption during cardiopulmonary testing (peak VO₂).

RESULTS

During graded exercise, FO, NFO, and CONT individuals all fatigued at similar SM HEP levels, as measured by ³¹P-MR. However, FO individuals fatigued fastest, with several-fold higher rates of PFE-induced HEP decline that correlated closely with shorter exercise duration in the MR scanner and with 6MW distance and lower peak oxygen consumption on cardiopulmonary testing (P < 0.001 for all). SM mitochondrial oxidative capacity was lower in older individuals and correlated with rapid HEP decline but less closely with EI.

CONCLUSION

Several-fold faster SM energetic decline during exercise occurs in FO individuals and correlates closely with multiple measures of EI. Rapid energetic decline represents an objective, functional measure of SM metabolic changes and a potential new target for mitigating frailty-associated physical limitations.

FUNDING

This work was supported by NIH R21 AG045634, R01 AG063661, R01 HL61912, the Johns Hopkins University Claude D. Pepper Older Americans Independence Center P30AG021334, and the Clarence Doodeman Endowment in Cardiology at Johns Hopkins.

Rapid exercise-induced skeletal muscle high-energy phosphate decline occurs in frail, older individuals and is closely linked to exercise intolerance and fatigue.

Use of cardiac magnetic resonance to detect changes in metabolism in heart failure

The heart has a massive adenosine triphosphate (ATP) requirement, produced from the oxidation of metabolic substrates such as fat and glucose. Magnetic resonance spectroscopy offers a unique opportunity to probe this biochemistry: ³¹Phosphorus spectroscopy can demonstrate the production of ATP and quantify levels of the transport molecule phosphocreatine while 13Carbon spectroscopy can demonstrate the metabolic fates of glucose in real time. These techniques allow the metabolic deficits in heart failure to be interrogated and can be a potential future clinical tool.

Non-invasive investigation of myocardial energetics in cardiac disease using 31P magnetic resonance spectroscopy

Cardiac metabolism and function are intrinsically linked. High-energy phosphates occupy a central and obligate position in cardiac metabolism, coupling oxygen and substrate fuel delivery to the myocardium with external work. This insight underlies the widespread clinical use of ischaemia testing. However, other deficits in high-energy phosphate metabolism (not secondary to supply-demand mismatch of oxygen and substrate fuels) may also be documented, and are of particular interest when found in the context of structural heart disease. This review introduces the scope of deficits in high-energy phosphate metabolism that may be observed in the myocardium, how to assess for them, and how they might be interpreted.

Cardiac Energetics in Patients With Aortic Stenosis and Preserved Versus Reduced Ejection Fraction

Supplemental Digital Content is available in the text.

Why some but not all patients with severe aortic stenosis (SevAS) develop otherwise unexplained reduced systolic function is unclear. We investigate the hypothesis that reduced creatine kinase (CK) capacity and flux is associated with this transition.

Creatine kinase rate constant in the human heart at 7T with 1D-ISIS/2D CSI localization

Purpose

Creatine Kinase (CK) reaction plays an important role in energy metabolism and estimate of its reaction rate constant in heart provides important insight into cardiac energetics. Fast saturation transfer method (T1nom–T1 nominal) to measure CK reaction rate constant (k_f) was previously demonstrated in open chest swine hearts. The goal of this work is to further develop this method for measuring the k_f in human myocardium at 7T. T1nom approach is combined with 1D-ISIS/2D-CSI for in vivo spatial localization and myocardial CK forward rate constant was then measured in 7 volunteers at 7T.

Methods

T1nom method uses two partially relaxed saturation transfer (ST) spectra and correction factor to determine CK rate constant. Correction factor is determined by numerical simulation of Bloch McConnell equations using known spin and experimental parameters. Optimal parameters and error estimate in calculation of CK reaction rate constant were determined by simulations. The technique was validated in calf muscles by direct comparison with saturation transfer measurements. T1nom pulse sequence was incorporated with 1D-image selected in vivo spectroscopy, combined with 2D-chemical shift spectroscopic imaging (1D-ISIS/2D-CSI) for studies in heart. The myocardial CK reaction rate constant was then measured in 7 volunteers.

Results

Skeletal muscle k_f determined by conventional approach and T1nom approach were the same 0.31 ± 0.02 s^-1 and 0.30 ± 0.04 s^-1 demonstrating the validity of the technique. Results are reported as mean ± SD. Myocardial CK reaction rate constant was 0.29 ± 0.05 s^-1, consistent with previously reported studies.

Conclusion

T1nom method enables acquisition of ³¹P saturation transfer MRS under partially relaxed conditions and enables 2D-CSI of k_f in myocardium. This work enables applications for in vivo CSI imaging of energetics in heart and other organs in clinically relevant acquisition time.

Myocardial Energetics in Obesity: Enhanced ATP Delivery Through Creatine Kinase With Blunted Stress Response

BACKGROUND

Obesity is strongly associated with exercise intolerance and the development of heart failure. Whereas myocardial energetics and diastolic function are impaired in obesity, systolic function is usually preserved. This suggests that the rate of ATP delivery is maintained, but this has never been explored in human obesity. We hypothesized that ATP transfer rate through creatine kinase (CK) (k_f^CKrest) would be increased, compensating for depleted energy stores (phosphocreatine/ATP), but potentially limiting greater ATP delivery during increased workload. We hypothesized that these changes would normalize with weight loss.

METHODS

We recruited 80 volunteers (35 controls [body mass index 24±3 kg/m²], 45 obese [body mass index 35±5 kg/m²]) without coexisting cardiovascular disease. Participants underwent body composition analysis, magnetic resonance imaging of abdominal, liver, and myocardial fat content, left ventricular function, and ³¹P magnetic resonance spectroscopy to assess phosphocreatine/ATP and CK kinetics, at rest and during dobutamine stress. Obese volunteers were assigned to a dietary weight loss intervention, before reexamination.

RESULTS

At rest, although myocardial phosphocreatine/ATP was 14% lower in obesity (1.9±0.3 versus 2.2±0.2, P<0.001), k_f^Ckrest was 33% higher (0.23±0.07 s^-1 versus 0.16±0.08 s^-1, P=0.002), yielding no difference in overall resting ATP delivery (obese 2.5±0.9 µmol·g^-1·s^-1 versus control 2.2±1.1 µmol·g^-1·s^-1, P=0.232). In controls, increasing cardiac workload led to an increase in both k_f^CK (+86%, P<0.001) and ATP delivery (+80%, P<0.001). However, in obesity, similar stress led to no significant increase in either k_f^CK (P=0.117) or ATP delivery (P=0.608). This was accompanied by reduced systolic augmentation (absolute increase in left ventricular ejection fraction, obese +16±7% versus control +21±4%, P=0.031). Successful weight loss (-11±5% body weight) was associated with improvement of these energetic changes such that there was no significant difference in comparison with controls.

CONCLUSIONS

In the obese resting heart, the myocardial CK reaction rate is increased, maintaining ATP delivery despite reduced phosphocreatine/ATP. During increased workload, although the nonobese heart increases ATP delivery through CK, the obese heart does not; this is associated with reduced systolic augmentation and exercise tolerance. Weight loss reverses these energetic changes. This highlights myocardial energy delivery through CK as a potential therapeutic target to improve symptoms in obesity-related heart disease, and a fascinating modifiable pathway involved in the progression to heart failure, as well.

Neural-network classification of cardiac disease from 31P cardiovascular magnetic resonance spectroscopy measures of creatine kinase energy metabolism

BACKGROUND

The heart's energy demand per gram of tissue is the body's highest and creatine kinase (CK) metabolism, its primary energy reserve, is compromised in common heart diseases. Here, neural-network analysis is used to test whether noninvasive phosphorus (31P) cardiovascular magnetic resonance spectroscopy (CMRS) measurements of cardiac adenosine triphosphate (ATP) energy, phosphocreatine (PCr), the first-order CK reaction rate kf, and the rate of ATP synthesis through CK (CK flux), can predict specific human heart disease and clinical severity.

METHODS

The data comprised the extant 178 complete sets of PCr and ATP concentrations, kf, and CK flux data from human CMRS studies performed on clinical 1.5 and 3 Tesla scanners. Healthy subjects and patients with nonischemic cardiomyopathy, dilated (DCM) or hypertrophic disease, New York Heart Association (NYHA) class I-IV heart failure (HF), or with anterior myocardial infarction are included. Three-layer neural-networks were created to classify disease and to differentiate DCM, hypertrophy and clinical NYHA class in HF patients using leave-one-out training. Network performance was assessed using 'confusion matrices' and 'area-under-the-curve' (AUC) analyses of 'receiver operating curves'. Possible methodological bias and network imbalance were tested by segregating 1.5 and 3 Tesla data, and by data augmentation by random interpolation of nearest neighbors, respectively.

RESULTS

The network differentiated healthy, HF and non-HF cardiac disease with an overall accuracy of 84% and AUC > 90% for each category using the four CK metabolic parameters, alone. HF patients with DCM, hypertrophy, and different NYHA severity were differentiated with ~ 80% overall accuracy independent of CMRS methodology.

CONCLUSIONS

While sample-size was limited in some sub-classes, a neural network classifier applied to noninvasive cardiac 31P CMRS data, could serve as a metabolic biomarker for common disease types and HF severity with clinically-relevant accuracy. Moreover, the network's ability to individually classify disease and HF severity using CK metabolism alone, implies an intimate relationship between CK metabolism and disease, with subtle underlying phenotypic differences that enable their differentiation.

TRIAL REGISTRATION

ClinicalTrials.gov Identifier: NCT00181259.

Localized rest and stress human cardiac creatine kinase reaction kinetics at 3 T

Changes in the kinetics of the creatine kinase (CK) shuttle are sensitive markers of cardiac energetics but are typically measured at rest and in the prone position. This study aims to measure CK kinetics during pharmacological stress at 3 T, with measurement in the supine position. A shorter "stressed saturation transfer" (StreST) extension to the triple repetition time saturation transfer (TRiST) method is proposed. We assess scanning in a supine position and validate the MR measurement against biopsy assay of CK activity. We report normal ranges of stress CK forward rate (kfCK ) for healthy volunteers and obese patients. TRiST measures kfCK in 40 min at 3 T. StreST extends the previously developed TRiST to also make a further kfCK measurement during <20 min of dobutamine stress. We test our TRiST implementation in skeletal muscle and myocardium in both prone and supine positions. We evaluate StreST in the myocardium of six healthy volunteers and 34 obese subjects. We validated MR-measured kfCK against biopsy assays of CK activity. TRiST kfCK values matched literature values in skeletal muscle (kfCK  = 0.25 ± 0.03 s-1 vs 0.27 ± 0.03 s-1 ) and myocardium when measured in the prone position (0.32 ± 0.15 s-1 ), but a significant difference was found for TRiST kfCK measured supine (0.24 ± 0.12 s-1 ). This difference was because of different respiratory- and cardiac-motion-induced B0 changes in the two positions. Using supine TRiST, cardiac kfCK values for normal-weight subjects were 0.15 ± 0.09 s-1 at rest and 0.17 ± 0.15 s-1 during stress. For obese subjects, kfCK was 0.16 ± 0.07 s-1 at rest and 0.17 ± 0.10 s-1 during stress. Rest myocardial kfCK and CK activity from LV biopsies of the same subjects correlated (R = 0.43, p = 0.03). We present an independent implementation of TRiST on the Siemens platform using a commercially available coil. Our extended StreST protocol enables cardiac kfCK to be measured during dobutamine-induced stress in the supine position.

Cardiac work is related to creatine kinase energy supply in human heart failure: a cardiovascular magnetic resonance spectroscopy study

BACKGROUND

It has been hypothesized that the supply of chemical energy may be insufficient to fuel normal mechanical pump function in heart failure (HF). The creatine kinase (CK) reaction serves as the heart's primary energy reserve, and the supply of adenosine triphosphate (ATP flux) it provides is reduced in human HF. However, the relationship between the CK energy supply and the mechanical energy expended has never been quantified in the human heart. This study tests whether reduced CK energy supply is associated with reduced mechanical work in HF patients.

METHODS

Cardiac mechanical work and CK flux in W/kg, and mechanical efficiency were measured noninvasively at rest using cardiac pressure-volume loops, magnetic resonance imaging and phosphorus spectroscopy in 14 healthy subjects and 27 patients with mild-to-moderate HF.

RESULTS

In HF, the resting CK flux (126 ± 46 vs. 179 ± 50 W/kg, p < 0.002), the average (6.8 ± 3.1 vs. 10.1 ± 1.5 W/kg, p  <0.001) and the peak (32 ± 14 vs. 48 ± 8 W/kg, p < 0.001) cardiac mechanical work-rates, as well as the cardiac mechanical efficiency (53% ± 16 vs. 79% ± 3, p < 0.001), were all reduced by a third compared to healthy subjects. In addition, cardiac CK flux correlated with the resting peak and average mechanical power (p < 0.01), and with mechanical efficiency (p = 0.002).

CONCLUSION

These first noninvasive findings showing that cardiac mechanical work and efficiency in mild-to-moderate human HF decrease proportionately with CK ATP energy supply, are consistent with the energy deprivation hypothesis of HF. CK energy supply exceeds mechanical work at rest but lies within a range that may be limiting with moderate activity, and thus presents a promising target for HF treatment.

TRIAL REGISTRATION

ClinicalTrials.gov Identifier: NCT00181259 .

31 P magnetic resonance fingerprinting for rapid quantification of creatine kinase reaction rate in vivo

The purpose of this work was to develop a 31 P spectroscopic magnetic resonance fingerprinting (MRF) method for fast quantification of the chemical exchange rate between phosphocreatine (PCr) and adenosine triphosphate (ATP) via creatine kinase (CK). A 31 P MRF sequence (CK-MRF) was developed to quantify the forward rate constant of ATP synthesis via CK ( kfCK), the T1 relaxation time of PCr ( T1PCr), and the PCr-to-ATP concentration ratio ( MRPCr). The CK-MRF sequence used a balanced steady-state free precession (bSSFP)-type excitation with ramped flip angles and a unique saturation scheme sensitive to the exchange between PCr and γATP. Parameter estimation was accomplished by matching the acquired signals to a dictionary generated using the Bloch-McConnell equation. Simulation studies were performed to examine the susceptibility of the CK-MRF method to several potential error sources. The accuracy of nonlocalized CK-MRF measurements before and after an ischemia-reperfusion (IR) protocol was compared with the magnetization transfer (MT-MRS) method in rat hindlimb at 9.4 T (n = 14). The reproducibility of CK-MRF was also assessed by comparing CK-MRF measurements with both MT-MRS (n = 17) and four angle saturation transfer (FAST) (n = 7). Simulation results showed that CK-MRF quantification of kfCK was robust, with less than 5% error in the presence of model inaccuracies including dictionary resolution, metabolite T2 values, inorganic phosphate metabolism, and B1 miscalibration. Estimation of kfCK by CK-MRF (0.38 ± 0.02 s-1 at baseline and 0.42 ± 0.03 s-1 post-IR) showed strong agreement with MT-MRS (0.39 ± 0.03 s-1 at baseline and 0.44 ± 0.04 s-1 post-IR). kfCK estimation was also similar between CK-MRF and FAST (0.38 ± 0.02 s-1 for CK-MRF and 0.38 ± 0.11 s-1 for FAST). The coefficient of variation from 20 s CK-MRF quantification of kfCK was 42% of that by 150 s MT-MRS acquisition and was 12% of that by 20 s FAST acquisition. This study demonstrates the potential of a 31 P spectroscopic MRF framework for rapid, accurate and reproducible quantification of chemical exchange rate of CK in vivo.

Assessing tissue metabolism by phosphorous-31 magnetic resonance spectroscopy and imaging: a methodology review

Many human diseases are caused by an imbalance between energy production and demand. Magnetic resonance spectroscopy (MRS) and magnetic resonance imaging (MRI) provide the unique opportunity for in vivo assessment of several fundamental events in tissue metabolism without the use of ionizing radiation. Of particular interest, phosphate metabolites that are involved in ATP generation and utilization can be quantified noninvasively by phosphorous-31 (31P) MRS/MRI. Furthermore, 31P magnetization transfer (MT) techniques allow in vivo measurement of metabolic fluxes via creatine kinase (CK) and ATP synthase. However, a major impediment for the clinical applications of 31P-MRS/MRI is the prohibitively long acquisition time and/or the low spatial resolution that are necessary to achieve adequate signal-to-noise ratio. In this review, current 31P-MRS/MRI techniques used in basic science and clinical research are presented. Recent advances in the development of fast 31P-MRS/MRI methods are also discussed.

OXSA: An open-source magnetic resonance spectroscopy analysis toolbox in MATLAB

In vivo magnetic resonance spectroscopy provides insight into metabolism in the human body. New acquisition protocols are often proposed to improve the quality or efficiency of data collection. Processing pipelines must also be developed to use these data optimally. Current fitting software is either targeted at general spectroscopy fitting, or for specific protocols. We therefore introduce the MATLAB-based OXford Spectroscopy Analysis (OXSA) toolbox to allow researchers to rapidly develop their own customised processing pipelines. The toolbox aims to simplify development by: being easy to install and use; seamlessly importing Siemens Digital Imaging and Communications in Medicine (DICOM) standard data; allowing visualisation of spectroscopy data; offering a robust fitting routine; flexibly specifying prior knowledge when fitting; and allowing batch processing of spectra. This article demonstrates how each of these criteria have been fulfilled, and gives technical details about the implementation in MATLAB. The code is freely available to download from https://github.com/oxsatoolbox/oxsa.

In-vivo31P-MRS of skeletal muscle and liver: A way for non-invasive assessment of their metabolism

In addition to direct assessment of high energy phosphorus containing metabolite content within tissues, phosphorus magnetic resonance spectroscopy (³¹P-MRS) provides options to measure phospholipid metabolites and cellular pH, as well as the kinetics of chemical reactions of energy metabolism in vivo. Even though the great potential of ³¹P-MR was recognized over 30 years ago, modern MR systems, as well as new, dedicated hardware and measurement techniques provide further opportunities for research of human biochemistry. This paper presents a methodological overview of the ³¹P-MR techniques that can be used for basic, physiological, or clinical research of human skeletal muscle and liver in vivo. Practical issues of ³¹P-MRS experiments and examples of potential applications are also provided. As signal localization is essential for liver ³¹P-MRS and is important for dynamic muscle examinations as well, typical localization strategies for ³¹P-MR are also described.

Creatine kinase rate constant in the human heart measured with 3D-localization at 7 tesla

PURPOSE

We present a new Bloch-Siegert four Angle Saturation Transfer (BOAST) method for measuring the creatine kinase (CK) first-order effective rate constant k_f in human myocardium at 7 tesla (T). BOAST combines a variant of the four-angle saturation transfer (FAST) method using amplitude-modulated radiofrequency pulses, phosphorus Bloch-Siegert B1+-mapping to determine the per-voxel flip angles, and nonlinear fitting to Bloch simulations for postprocessing.

METHODS

Optimal flip angles and repetition time parameters were determined from Monte Carlo simulations. BOAST was validated in the calf muscle of two volunteers at 3T and 7T. The myocardial CK forward rate constant was then measured in 10 volunteers at 7T in 82 min (after ¹ H localization).

RESULTS

BOAST kfCK values were 0.281 ± 0.002 s^-1 in the calf and 0.35 ± 0.05 s^-1 in myocardium. These are consistent with literature values from lower fields. Using a literature values for adenosine triphosphate concentration, we computed CK flux values of 4.55 ± 1.52 mmol kg^-1 s^-1 . The sensitive volume for BOAST depends on the B₁ inhomogeneity of the transmit coil.

CONCLUSION

BOAST enables measurement of the CK rate constant in the human heart at 7T, with spatial localization in three dimensions to 5.6 mL voxels, using a 10-cm loop coil. Magn Reson Med 78:20-32, 2017. © 2016 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Test-retest analysis of multiple 31 P magnetization exchange pathways using asymmetric adiabatic inversion

PURPOSE

A ³¹ P-MR inversion transfer (IT) method with a short adiabatic inversion pulse is proposed and its test-retest reliability was evaluated for two spectral fitting strategies.

METHODS

Assessment in a test-retest design (3 Tesla, vastus muscles, 12 healthy volunteers, 14 inversion times, 22 ms asymmetric adiabatic inversion pulse, adiabatic excitation); spectral fitting in Fitting Tool for Interrelated Arrays of Datasets (FitAID) and Java Magnetic Resonance User Interface (jMRUI); least squares solution of the Bloch-McConnell-Solomon matrix formalism including all 14 measured time-points with equal weighting.

RESULTS

The cohort averages of k[PCr→γ-ATP] (phosphocreatine, PCr; adenosine triphosphate, ATP) are 0.246 ± 0.050s^-1 versus 0.254 ± 0.050s^-1 , and k[Pi→γ-ATP] 0.086 ± 0.033s^-1 versus 0.066 ± 0.034s^-1 (average ± standard deviation, jMRUI versus FitAID). Coefficients of variation of the differences between test and retest are lowest (9.5%) for k[PCr→γ-ATP] fitted in FitAID, larger (15.2%) for the fit in jMRUI, and considerably larger for k[Pi→γ-ATP] fitted in FitAID (43.4%) or jMRUI (47.9%). The beginning of the IT effect can be observed with magnetizations above 92% for noninverted lines while inversion of the ATP resonances is better than -72%.

CONCLUSION

The performance of the asymmetric adiabatic pulse allows an accurate observation of IT effects even in the early phase; the least squares fit of the Bloch-McConnell-Solomon matrix formalism is robust; and the type of spectral fitting can influence the results significantly. Magn Reson Med 78:33-39, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Feasibility and repeatability of localized (31) P-MRS four-angle saturation transfer (FAST) of the human gastrocnemius muscle using a surface coil at 7 T

Phosphorus ((31) P) MRS, combined with saturation transfer (ST), provides non-invasive insight into muscle energy metabolism. However, even at 7 T, the standard ST method with T1 (app) measured by inversion recovery takes about 10 min, making it impractical for dynamic examinations. An alternative method, i.e. four-angle saturation transfer (FAST), can shorten the examination time. The aim of this study was to test the feasibility, repeatability, and possible time resolution of the localized FAST technique measurement on an ultra-high-field MR system, to accelerate the measurement of both Pi -to-ATP and PCr-to-ATP reaction rates in the human gastrocnemius muscle and to test the feasibility of using the FAST method for dynamic measurements. We measured the exchange rates and metabolic fluxes in the gastrocnemius muscle of eight healthy subjects at 7 T with the depth-resolved surface coil MRS (DRESS)-localized FAST method. For comparison, a standard ST localized method was also used. The measurement time for the localized FAST experiment was 3.5 min compared with the 10 min for the standard localized ST experiment. In addition, in five healthy volunteers, Pi -to-ATP and PCr-to-ATP metabolic fluxes were measured in the gastrocnemius muscle at rest and during plantar flexion by the DRESS-localized FAST method. The repeatability of PCr-to-ATP and Pi -to-ATP exchange rate constants, determined by the slab-selective localized FAST method at 7 T, is high, as the coefficients of variation remained below 20%, and the results of the exchange rates measured with the FAST method are comparable to those measured with standard ST. During physical activity, the PCr-to-ATP metabolic flux decreased (from FCK  = 8.21 ± 1.15 mM s(-1) to FCK  = 3.86 ± 1.38 mM s(-1) ) and the Pi -to-ATP flux increased (from FATP  = 0.43 ± 0.14 mM s(-1) to FATP  = 0.74 ± 0.13 mM s(-1) ). In conclusion, we could demonstrate that measurements in the gastrocnemius muscle are feasible at rest and are short enough to be used during exercise with the DRESS-localized FAST method at 7 T.

Two repetition time saturation transfer (TwiST) with spill-over correction to measure creatine kinase reaction rates in human hearts

BACKGROUND

Phosphorus saturation transfer (ST) magnetic resonance spectroscopy can measure the rate of ATP generated from phosphocreatine (PCr) via creatine kinase (CK) in the human heart. Recently, the triple-repetition time ST (TRiST) method was introduced to measure the CK pseudo-first-order rate constant kf in three acquisitions. In TRiST, the longitudinal relaxation time of PCr while γ-ATP is saturated, T1`, is measured for each subject, but suffers from low SNR because the PCr signal is reduced due to exchange with saturated γ-ATP, and the short repetition time of one of the acquisitions. Here, a two-repetition time ST (TwiST) method is presented. In TwiST, the acquisition with γ-ATP saturation and short repetition time is dropped. Instead of measuring T1`, an intrinsic relaxation time T1 for PCr, T1 (intrinsic), is assumed. The objective was to validate TwiST measurements of CK kinetics in healthy subjects and patients with heart failure (HF).

METHODS

Bloch equation simulations that included the effect of spillover irradiation on PCr were used to derive formulae for T1 (intrinsic) and kf measured by both TRiST and TwiST methods. Spillover was quantified from an unsaturated PCr measurement used in the current protocol for determining PCr and ATP concentrations. Cardiac TRiST and TwiST data were acquired at 3 T from 12 healthy and 17 HF patients.

RESULTS

Simulations showed that both kf measured by TwiST and T1 (intrinsic) require spill-over corrections. In human heart at 3 T, the spill-over corrected T1 (intrinsic) = 8.4 ± 1.4 s (mean ± SD) independent of study group. TwiST and TRiST kf measurements were the same, but TwiST was 9 min faster. Spill-over corrected TwiST kf was 0.33 ± 0.08 s(-1) vs. 0.20 ± 0.06 s(-1) in healthy vs HF hearts, respectively (p < 0.0001).

CONCLUSION

TwiST was validated against TRiST in the human heart at 3 T, generating the same results 9 min faster. TwiST detected significant reductions in CK kf in HF compared to healthy subjects, consistent with prior 1.5 T studies using different methodology.

Small animal cardiovascular MR imaging and spectroscopy

The use of MR imaging and spectroscopy for studying cardiovascular disease processes in small animals has increased tremendously over the past decade. This is the result of the remarkable advances in MR technologies and the increased availability of genetically modified mice. MR techniques provide a window on the entire timeline of cardiovascular disease development, ranging from subtle early changes in myocardial metabolism that often mark disease onset to severe myocardial dysfunction associated with end-stage heart failure. MR imaging and spectroscopy techniques play an important role in basic cardiovascular research and in cardiovascular disease diagnosis and therapy follow-up. This is due to the broad range of functional, structural and metabolic parameters that can be quantified by MR under in vivo conditions non-invasively. This review describes the spectrum of MR techniques that are employed in small animal cardiovascular disease research and how the technological challenges resulting from the small dimensions of heart and blood vessels as well as high heart and respiratory rates, particularly in mice, are tackled.

On the theoretical limits of detecting cyclic changes in cardiac high-energy phosphates and creatine kinase reaction kinetics using in vivo ³¹P MRS

Adenosine triphosphate (ATP) is absolutely required to fuel normal cyclic contractions of the heart. The creatine kinase (CK) reaction is a major energy reserve reaction that rapidly converts creatine phosphate (PCr) to ATP during the cardiac cycle and at times of stress and ischemia, but is significantly impaired in conditions such as hypertrophy and heart failure. Because the magnitudes of possible in vivo cyclic changes in cardiac high-energy phosphates (HEPs) during the cardiac cycle are not well known from previous work, this study uses mathematical modeling to assess whether, and to what extent, cyclic variations in HEPs and in the rate of ATP synthesis through CK (CK flux) could exist in the human heart, and whether they could be measured with current in vivo (31)P MRS methods. Multi-site exchange models incorporating enzymatic rate equations were used to study the cyclic dynamics of the CK reaction, and Bloch equations were used to simulate (31)P MRS saturation transfer measurements of the CK reaction. The simulations show that short-term buffering of ATP by CK requires temporal variations over the cardiac cycle in the CK reaction velocities modeled by enzymatic rate equations. The maximum variation in HEPs in the normal human heart beating at 60 min(-1) was approximately 0.4 mM and proportional to the velocity of ATP hydrolysis. Such HEP variations are at or below the current limits of detection by in vivo (31)P MRS methods. Bloch equation simulations show that (31)P MRS saturation transfer estimates the time-averaged, pseudo-first-order forward rate constant, k(f,ap)', of the CK reaction, and that periodic short-term fluctuations in kf ' and CK flux are not likely to be detectable in human studies employing current in vivo (31)P MRS methods.

Suppression of skeletal muscle signal using a crusher coil: A human cardiac (31) p-MR spectroscopy study at 7 tesla

Purpose

The translation of sophisticated phosphorus MR spectroscopy (³¹P‐MRS) protocols to 7 Tesla (T) is particularly challenged by the issue of radiofrequency (RF) heating. Legal limits on RF heating make it hard to reliably suppress signals from skeletal muscle that can contaminate human cardiac ³¹P spectra at 7T. We introduce the first surface‐spoiling crusher coil for human cardiac ³¹P‐MRS at 7T.

Methods

A planar crusher coil design was optimized with simulations and its performance was validated in phantoms. Crusher gradient pulses (100 μs) were then applied during human cardiac ³¹P‐MRS at 7T.

Results

In a phantom, residual signals were 50 ± 10% with BISTRO (B₁‐insensitive train to obliterate signal), and 34 ± 8% with the crusher coil. In vivo, residual signals in skeletal muscle were 49 ± 4% using BISTRO, and 24 ± 5% using the crusher coil. Meanwhile, in the interventricular septum, spectral quality and metabolite quantification did not differ significantly between BISTRO (phosphocreatine/adenosine triphosphate [PCr/ATP] = 2.1 ± 0.4) and the crusher coil (PCr/ATP = 1.8 ± 0.4). However, the specific absorption rate (SAR) decreased from 96 ± 1% of the limit (BISTRO) to 16 ± 1% (crusher coil).

Conclusion

A crusher coil is an SAR‐efficient alternative for selectively suppressing skeletal muscle during cardiac ³¹P‐MRS at 7T. A crusher coil allows the use of sequence modules that would have been SAR‐prohibitive, without compromising skeletal muscle suppression. Magn Reson Med 75:962–972, 2016. © 2015 The Authors. Magnetic Resonance in Medicine Published by Wiley Periodicals, Inc. on behalf of International Society of Medicine in Resonance.

Comparison of (31)P saturation and inversion magnetization transfer in human liver and skeletal muscle using a clinical MR system and surface coils

(31)P MRS magnetization transfer ((31)P-MT) experiments allow the estimation of exchange rates of biochemical reactions, such as the creatine kinase equilibrium and adenosine triphosphate (ATP) synthesis. Although various (31)P-MT methods have been successfully used on isolated organs or animals, their application on humans in clinical scanners poses specific challenges. This study compared two major (31)P-MT methods on a clinical MR system using heteronuclear surface coils. Although saturation transfer (ST) is the most commonly used (31)P-MT method, sequences such as inversion transfer (IT) with short pulses might be better suited for the specific hardware and software limitations of a clinical scanner. In addition, small NMR-undetectable metabolite pools can transfer MT to NMR-visible pools during long saturation pulses, which is prevented with short pulses. (31)P-MT sequences were adapted for limited pulse length, for heteronuclear transmit-receive surface coils with inhomogeneous B1 , for the need for volume selection and for the inherently low signal-to-noise ratio (SNR) on a clinical 3-T MR system. The ST and IT sequences were applied to skeletal muscle and liver in 10 healthy volunteers. Monte-Carlo simulations were used to evaluate the behavior of the IT measurements with increasing imperfections. In skeletal muscle of the thigh, ATP synthesis resulted in forward reaction constants (k) of 0.074 ± 0.022 s(-1) (ST) and 0.137 ± 0.042 s(-1) (IT), whereas the creatine kinase reaction yielded 0.459 ± 0.089 s(-1) (IT). In the liver, ATP synthesis resulted in k = 0.267 ± 0.106 s(-1) (ST), whereas the IT experiment yielded no consistent results. ST results were close to literature values; however, the IT results were either much larger than the corresponding ST values and/or were widely scattered. To summarize, ST and IT experiments can both be implemented on a clinical body scanner with heteronuclear transmit-receive surface coils; however, ST results are much more robust against experimental imperfections than the current implementation of IT.

Three-dimensional saturation transfer ³¹P-MRI in muscles of the lower leg at 3.0 T

The creatine kinase (CK) reaction plays a critical role in skeletal muscle function, and can be studied non-invasively using phosphorus (³¹P) saturation transfer (ST) techniques. However, due to the low MR sensitivity of the ³¹P nucleus, most studies on clinically approved magnetic fields (≤3.0 T) have been performed with coarse resolution and limited tissue coverage. However, such methods are not able to detect spatially resolved metabolic heterogeneities, which may be important in diseases of the skeletal muscle. In this study, our aim was to develop and implement a ³¹P-MRI method for mapping the kinetics of the CK reaction, and the unidirectional phosphocreatine (PCr) to adenosine triphosphate (ATP) metabolic fluxes in muscles of the lower leg on a clinical 3.0 T MR scanner. We imaged the lower leg muscles of ten healthy volunteers (total experimental time: 40 min, nominal voxel sizes 0.5 mL), and found statistically significant differences between the kinetics of the CK reaction among muscle groups. Our developed technique may allow in the future the early detection of focal metabolic abnormalities in diseases that affect the function of the skeletal muscle.

Reproducibility of creatine kinase reaction kinetics in human heart: a (31) P time-dependent saturation transfer spectroscopy study

Creatine kinase (CK) is essential for the buffering and rapid regeneration of adenosine triphosphate (ATP) in heart tissue. Herein, we demonstrate a (31) P MRS protocol to quantify CK reaction kinetics in human myocardium at 3 T. Furthermore, we sought to quantify the test-retest reliability of the measured metabolic parameters. The method localizes the (31) P signal from the heart using modified one-dimensional image-selected in vivo spectroscopy (ISIS), and a time-dependent saturation transfer (TDST) approach was used to measure CK reaction parameters. Fifteen healthy volunteers (22 measurements in total) were tested. The CK reaction rate constant (kf ) was 0.32 ± 0.05 s(-1) and the coefficient of variation (CV) was 15.62%. The intrinsic T1 for phosphocreatine (PCr) was 7.36 ± 1.79 s with CV = 24.32%. These values are consistent with those reported previously. The PCr/ATP ratio was equal to 1.94 ± 0.15 with CV = 7.73%, which is within the range of healthy subjects. The reproducibility of the technique was tested in seven subjects and inferred parameters, such as kf and T1 , exhibited good reliability [intraclass correlation coefficient (ICC) of 0.90 and 0.79 for kf and T1 , respectively). The reproducibility data provided in this study will enable the calculation of the power and sample sizes required for clinical and research studies. The technique will allow for the examination of cardiac energy metabolism in clinical and research studies, providing insight into the relationship between energy deficit and functional deficiency in the heart.

Exchange kinetics by inversion transfer: integrated analysis of the phosphorus metabolite kinetic exchanges in resting human skeletal muscle at 7 T

PURPOSE

To develop an inversion pulse-based, chemical exchange saturation transfer-like method for detection of (31) P magnetization exchanges among all nuclear magnetic resonance visible metabolites suitable for providing an integrated kinetic analysis of phosphorus exchange reactions in vivo.

METHODS

The exchange kinetics by inversion transfer (EKIT) sequence includes application of a frequency-selective inversion pulse arrayed over the range of relevant (31) P frequencies, followed by a constant delay and a hard readout pulse. A series of EKIT spectra, each given by a plot of Z-magnetization for each metabolite of interest versus frequency of the inversion pulse, can be generated from this single data set.

RESULTS

EKIT spectra reflect chemical exchange due to known biochemical reactions, cross-relaxation effects, and relayed magnetization transfers due to both processes. The rate constants derived from EKIT data collected on resting human skeletal muscle were: ATP synthesis via ATP synthase (0.050 ± 0.016 s(-1) ), ATP synthesis via creatine kinase (0.264 ± 0.023 s(-1) ), and cross-relaxation between neighboring spin pairs within ATP (0.164 ± 0.022 s(-1) ).

CONCLUSION

EKIT provides a simple, alternative method to detect chemical exchange, cross relaxation, and relayed magnetization transfer effects in human skeletal muscle at 7 T.

One-dimensional image-selected in vivo spectroscopy localized phosphorus saturation transfer at 7T

PURPOSE

To evaluate the feasibility of a one-dimensional image-selected in vivo spectroscopy (1D-ISIS) saturation transfer (ST) sequence at 7T for localized in vivo measurements of energy metabolism in different tissues in clinically reasonable examination times.

METHODS

The performance of a gradient offset independent adiabacity-based 1D-ISIS localization was tested on phantom and the localized ST sequence was compared with the nonlocalized version in vivo. We performed localized measurements of basal metabolism of human liver and different muscle groups of the calf. Localized ST experiments took 15-25 minutes.

RESULTS

The selectivity of the 1D-ISIS sequence was 81.63% and the outer volume suppression was 97.57%. The ST parameters acquired with the 1D-ISIS sequence and with the nonlocalized acquisition in the muscle were not statistically different. The forward rate constants for phosphocreatine (PCr)-adenosine triphosphate (ATP) and inorganic phosphate (Pi)-ATP exchange reactions were measured in the soleus (kCK  = 0.30 ± 0.06 s(-1) and kATP  = 0.11 ± 0.02 s(-1) , respectively) and in the medial gastrocnemius (kCK  = 0.27 ± 0.06 s(-1) and kATP  = 0.09 ± 0.03s(-1) , respectively) in 15 minutes per muscle group. The corresponding fluxes were FCK  = 6.26 ± 1.28 μmol/g/s, FATP  = 0.22 ± 0.05 μmol/g/s and FCK  = 6.29 ± 1.66 μmol/g/s, FATP  = 0.21 ± 0.07 μmol/g/s, for soleus and gastrocnemius, respectively. The hepatic ATP synthesis measurement was feasible in 24 minutes.

CONCLUSION

The fast assessment of PCr-ATP and Pi-ATP exchange rates at 7T makes the 1D-ISIS ST sequence a promising tool for examining local resting-state metabolism in clinically acceptable measurement times.

In vivo evidence for cerebral bioenergetic abnormalities in schizophrenia measured using 31P magnetization transfer spectroscopy

IMPORTANCE

Abnormalities in neural activity and cerebral bioenergetics have been observed in schizophrenia (SZ). Further defining energy metabolism anomalies would provide crucial information about molecular mechanisms underlying SZ and may be valuable for developing novel treatment strategies.

OBJECTIVE

To investigate cerebral bioenergetics in SZ via measurement of creatine kinase activity using in vivo 31P magnetization transfer spectroscopy.

DESIGN, SETTING, AND PARTICIPANTS

Cross-sectional case-control study in the setting of clinical services and a brain imaging center of an academic psychiatric hospital. Twenty-six participants with chronic SZ (including a subgroup diagnosed as having schizoaffective disorder) and 26 age-matched and sex-matched healthy control subjects (25 usable magnetic resonance spectroscopy data sets from the latter).

INTERVENTION

31P magnetization transfer spectroscopy.

MAIN OUTCOMES AND MEASURES

The primary outcome measure was the forward rate constant (k(f)) of the creatine kinase enzyme in the frontal lobe. We also collected independent measures of brain intracellular pH and steady-state metabolite ratios of high-energy phosphate-containing compounds (phosphocreatine and adenosine triphosphate [ATP]), inorganic phosphate, and the 2 membrane phospholipids phosphodiester and phosphomonoester.

RESULTS

A substantial (22%) and statistically significant (P = .003) reduction in creatine kinase kf was observed in SZ. In addition, intracellular pH was significantly reduced (7.00 in the SZ group vs 7.03 in the control group, P = .007) in this condition. The phosphocreatine to ATP ratio, inorganic phosphate to ATP ratio, and phosphomonoester to ATP ratio were not substantially altered in SZ, but a significant (P = .02) reduction was found in the phosphodiester to ATP ratio. The abnormalities were similar between SZ and schizoaffective disorder.

CONCLUSIONS AND RELEVANCE

Using a novel 31P magnetization transfer magnetic resonance spectroscopy approach, we provide direct and compelling evidence for a specific bioenergetic abnormality in SZ. Reduced kf of the creatine kinase enzyme is consistent with an abnormality in storage and use of brain energy. The intracellular pH reduction suggests a relative increase in the contribution of glycolysis to ATP synthesis, providing convergent evidence for bioenergetic abnormalities in SZ. The similar phosphocreatine to ATP ratios in SZ and healthy controls suggest that the underlying bioenergetics abnormality is not associated with change in this metabolite ratio.

Metabolic rates of ATP transfer through creatine kinase (CK Flux) predict clinical heart failure events and death

Morbidity and mortality from heart failure (HF) are high, and current risk stratification approaches for predicting HF progression are imperfect. Adenosine triphosphate (ATP) is required for normal cardiac contraction, and abnormalities in creatine kinase (CK) energy metabolism, the primary myocardial energy reserve reaction, have been observed in experimental and clinical HF. However, the prognostic value of abnormalities in ATP production rates through CK in human HF has not been investigated. Fifty-eight HF patients with nonischemic cardiomyopathy underwent ³¹P magnetic resonance spectroscopy (MRS) to quantify cardiac high-energy phosphates and the rate of ATP synthesis through CK (CK flux) and were prospectively followed for a median of 4.7 years. Multiple-event analysis (MEA) was performed for HF-related events including all-cause and cardiac death, HF hospitalization, cardiac transplantation, and ventricular-assist device placement. Among baseline demographic, clinical, and metabolic parameters, MEA identified four independent predictors of HF events: New York Heart Association (NYHA) class, left ventricular ejection fraction (LVEF), African-American race, and CK flux. Reduced myocardial CK flux was a significant predictor of HF outcomes, even after correction for NYHA class, LVEF, and race. For each increase in CK flux of 1 μmol g⁻¹ s⁻¹, risk of HF-related composite outcomes decreased by 32 to 39%. These findings suggest that reduced CK flux may be a potential HF treatment target. Newer imaging strategies, including noninvasive ³¹P MRS that detect altered ATP kinetics, could thus complement risk stratification in HF and add value in conditions involving other tissues with high energy demands, including skeletal muscle and brain.

Quantification of human high-energy phosphate metabolite concentrations at 3 T with partial volume and sensitivity corrections

Practical noninvasive methods for the measurement of absolute metabolite concentrations are key to the assessment of the depletion of myocardial metabolite pools which occurs with several cardiac diseases, including infarction and heart failure. Localized MRS offers unique noninvasive access to many metabolites, but is often confounded by nonuniform sensitivity and partial volume effects in the large, poorly defined voxels commonly used for the detection of low-concentration metabolites with surface coils. These problems are exacerbated at higher magnetic field strengths by greater radiofrequency (RF) field inhomogeneity and differences in RF penetration with heteronuclear concentration referencing. An example is the (31)P measurement of cardiac adenosine triphosphate (ATP) and phosphocreatine (PCr) concentrations, which, although central to cardiac energetics, have not been measured at field strengths above 1.5 T. Here, practical acquisition and analysis protocols are presented for the quantification of [PCr] and [ATP] with one-dimensionally resolved surface coil spectra and concentration referencing at 3 T. The effects of nonuniform sensitivity and partial tissue volumes are addressed at 3 T by the application of MRI-based three-dimensional sensitivity weighting and tissue segmentation. The method is validated in phantoms of different sizes and concentrations, and used to measure [PCr] and [ATP] in healthy subjects. In calf muscle (n = 8), [PCr] = 24.7 ± 3.4 and [ATP] = 5.7 ± 1.3 µmol/g wet weight, whereas, in heart (n = 18), [PCr] = 10.4 ± 1.5 and [ATP] = 6.0 ± 1.1 µmol/g wet weight (all mean ± SD), consistent with previous reports at lower fields. The method enables, for the first time, the efficient, semi-automated quantification of high-energy phosphate metabolites in humans at 3 T with nonuniform excitation and detection.

Human cardiac 31P magnetic resonance spectroscopy at 7 Tesla

Purpose

Phosphorus magnetic resonance spectroscopy (³¹P-MRS) affords unique insight into cardiac energetics but has a low intrinsic signal-to-noise ratio (SNR) in humans. Theory predicts an increased ³¹P-MRS SNR at 7T, offering exciting possibilities to better investigate cardiac metabolism. We therefore compare the performance of human cardiac ³¹P-MRS at 7T to 3T, and measure T₁s for ³¹P metabolites at 7T.

Methods

Matched ³¹P-MRS data were acquired at 3T and 7T, on nine normal volunteers. A novel Look-Locker CSI acquisition and fitting approach was used to measure T₁s on six normal volunteers.

Results

T₁s in the heart at 7T were: phosphocreatine (PCr) 3.05 ± 0.41s, γ-ATP 1.82 ± 0.09s, α-ATP 1.39 ± 0.09s, β-ATP 1.02 ± 0.17s and 2,3-DPG (2,3-diphosphoglycerate) 3.05 ± 0.41s (N = 6). In the field comparison (N = 9), PCr SNR increased 2.8× at 7T relative to 3T, the Cramer-Ráo uncertainty (CRLB) in PCr concentration decreased 2.4×, the mean CRLB in PCr/ATP decreased 2.7× and the PCr/ATP SD decreased 2×.

Conclusion

Cardiac ³¹P-MRS at 7T has higher SNR and the spectra can be quantified more precisely than at 3T. Cardiac ³¹P T₁s are shorter at 7T than at 3T. We predict that 7T will become the field strength of choice for cardiac ³¹P-MRS. Magn Reson Med 72:304–315, 2014. © 2013 The Authors. Magnetic Resonance in Medicine Published by Wiley Periodicals, Inc. on behalf of International Society of Medicine in Resonance. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution, and reproduction in any medium, provided the original work is properly cited.

Three-dimensional mapping of the creatine kinase enzyme reaction rate in muscles of the lower leg

Phosphorus ((31) P) magnetization transfer (MT) techniques enable the non-invasive measurement of metabolic turnover rates of important enzyme-catalyzed reactions, such as the creatine kinase reaction (CK), a major transducing reaction involving adenosine triphosphate and phosphocreatine. Alteration in the kinetics of the CK reaction rate appears to play a central role in many disease states. In this study, we developed and implemented at ultra-high field (7T) a novel three-dimensional (31) P-MT imaging sequence that maps the kinetics of CK in the entire volume of the lower leg at relatively high resolution (0.52 mL voxel size), and within acquisition times that can be tolerated by patients (below 60 min). We tested the sequence on five healthy and two clinically diagnosed type 2 diabetic subjects. Overall, we obtained measurements that are in close agreement with measurements reported previously using spectroscopic methods. Importantly, our spatially resolved method allowed us to measure local CK reaction rate constants and metabolic fluxes in individual muscles in a non-invasive manner. Furthermore, it allowed us to detect variations of the CK rates of different muscles, which would not have been possible using unlocalized MRS methods. The results of this work suggest that 3D mapping of the CK reaction rates and metabolic fluxes can be achieved in the skeletal muscle in vivo at relatively high spatial resolution and with acquisition times well tolerated by patients. The ability to measure bioenergetics simultaneously in large areas of muscles will bring new insights into possible heterogeneous patterns of muscle metabolism associated with several diseases and serve as a valuable tool for monitoring the efficacy of interventions.

Magnetic resonance imaging and spectroscopy of the murine cardiovascular system

Magnetic resonance imaging (MRI) has emerged as a powerful and reliable tool to noninvasively study the cardiovascular system in clinical practice. Because transgenic mouse models have assumed a critical role in cardiovascular research, technological advances in MRI have been extended to mice over the last decade. These have provided critical insights into cardiac and vascular morphology, function, and physiology/pathophysiology in many murine models of heart disease. Furthermore, magnetic resonance spectroscopy (MRS) has allowed the nondestructive study of myocardial metabolism in both isolated hearts and in intact mice. This article reviews the current techniques and important pathophysiological insights from the application of MRI/MRS technology to murine models of cardiovascular disease.

Creatine kinase and ATP synthase reaction rates in human frontal lobe measured by ³¹P magnetization transfer spectroscopy at 4T

The human frontal lobe is critical for cognitive function in the healthy brain. Many psychiatric disorders including schizophrenia and bipolar disorder are associated with apparent mitochondrial dysfunction and bioenergetic abnormalities in the frontal lobe. Therefore, measuring cerebral bioenergetics associated with creatine kinase and adenosine triphosphate (ATP) synthase reactions could provide crucial information regarding the underlying molecular mechanisms associated with psychiatric disorders. In this study, the unidirectional forward chemical exchange metabolic fluxes of creatine kinase and ATP synthase reactions as well as reverse chemical exchange metabolic flux associated with ATP hydrolysis were determined at 4T by (31)P magnetization transfer. The current experiments indicate that the kinetic network of PCr↔ATP↔Pi can be measured reliably in the human frontal lobe at 4T, which will enable detailed in vivo characterization of bioenergetic abnormalities in a variety of neuropsychiatric disorders.

Creatine kinase-mediated improvement of function in failing mouse hearts provides causal evidence the failing heart is energy starved

ATP is required for normal cardiac contractile function, and it has long been hypothesized that reduced energy delivery contributes to the contractile dysfunction of heart failure (HF). Despite experimental and clinical HF data showing reduced metabolism through cardiac creatine kinase (CK), the major myocardial energy reserve and temporal ATP buffer, a causal relationship between reduced ATP-CK metabolism and contractile dysfunction in HF has never been demonstrated. Here, we generated mice conditionally overexpressing the myofibrillar isoform of CK (CK-M) to test the hypothesis that augmenting impaired CK-related energy metabolism improves contractile function in HF. CK-M overexpression significantly increased ATP flux through CK ex vivo and in vivo but did not alter contractile function in normal mice. It also led to significantly increased contractile function at baseline and during adrenergic stimulation and increased survival after thoracic aortic constriction (TAC) surgery-induced HF. Withdrawal of CK-M overexpression after TAC resulted in a significant decline in contractile function as compared with animals in which CK-M overexpression was maintained. These observations provide direct evidence that the failing heart is "energy starved" as it relates to CK. In addition, these data identify CK as a promising therapeutic target for preventing and treating HF and possibly diseases involving energy-dependent dysfunction in other organs with temporally varying energy demands.

High-energy phosphate transfer in human muscle: diffusion of phosphocreatine

The creatine kinase (CK) reaction is central to muscle energetics, buffering ATP levels during periods of intense activity via consumption of phosphocreatine (PCr). PCr is believed to serve as a spatial shuttle of high-energy phosphate between sites of energy production in the mitochondria and sites of energy utilization in the myofibrils via diffusion. Knowledge of the diffusion coefficient of PCr (D(PCr)) is thus critical for modeling and understanding energy transport in the myocyte, but D(PCr) has not been measured in humans. Using localized phosphorus magnetic resonance spectroscopy, we measured D(PCr) in the calf muscle of 11 adults as a function of direction and diffusion time. The results show that the diffusion of PCr is anisotropic, with significantly higher diffusion along the muscle fibers, and that the diffusion of PCr is restricted to a ∼28-μm pathlength assuming a cylindrical model, with an unbounded diffusion coefficient of ∼0.69 × 10(-3) mm(2)/s. This distance is comparable in size to the myofiber radius. On the basis of prior measures of CK reaction kinetics in human muscle, the expected diffusion distance of PCr during its half-life in the CK reaction is ∼66 μm. This distance is much greater than the average distances between mitochondria and myofibrils. Thus these first measurements of PCr diffusion in human muscle in vivo support the view that PCr diffusion is not a factor limiting high-energy phosphate transport between the mitochondria and the myofibrils in healthy resting myocytes.

Impaired ATP kinetics in failing in vivo mouse heart

BACKGROUND

The hypothesis that the failing heart may be energy-starved is supported in part by observations of reduced rates of adenosine 5'-triphosphate (ATP) synthesis through the creatine kinase (CK) reaction, the primary myocardial energy reservoir, in patients with heart failure (HF). Although murine models have been used to probe HF pathophysiology, it has not been possible to noninvasively measure the rate of ATP synthesis through CK in the in vivo mouse heart. The purpose of this work was to exploit noninvasive spatially localized magnetic resonance spectroscopy techniques to measure ATP flux through CK in in vivo mouse hearts and determine the extent of any reductions in murine HF.

METHODS AND RESULTS

The Triple Repetition Time Saturation Transfer (TRiST) magnetic resonance spectroscopy method of measuring ATP kinetics was first validated in skeletal muscle, rendering similar results to conventional saturation transfer magnetic resonance spectroscopy. In normal mouse hearts, the in vivo CK pseudo-first-order-rate constant, k(F), was 0.32±0.03 s(-1) (mean±SD) and the rate of ATP synthesis through CK was 3.16±0.47 μmol/g/s. Thoracic aortic constriction reduced k(F) by 31% (0.23±0.03 s(-1), P<0.0001) and ATP synthesis through CK by 51% (1.54±0.25 μmol/g/s, P<0.0001), values analogous to those in failing human hearts.

CONCLUSIONS

Despite the small size and high murine heart rate, the ATP synthesis rate through CK is similar in vivo in murine and human hearts and comparably reduced in HF. Because murine thoracic aortic constriction shares fundamental energetic similarities with human HF, this model and new magnetic resonance spectroscopy approach promise a powerful means to noninvasively probe altered energetics in HF.