Abstract:Introduction Magnetization saturation transfer (ST) is the most commonly used method in P MRS to non-invasively measure the kinetics of phosphoryl-exchange reactions and thus especially suitable for in vivo enzyme kinetic studies such as creatine kinase (CK) and ATPase. However, the conventional ST approach is greatly limited by the lengthy data acquisition time for achieving accurate quantification of rate constant and flux. Aside from the multi-point progressive saturation experiment that retrieves apparent T1 by curve fitting, the steady-state saturation transfer experiment alone is time-consuming enough because it usually requires a very long pre-saturation delay followed by a long saturation time in order to achieve the true steady-state magnetization. This would greatly compromise its in vivo applications in which scan number is large (i.e. studies requiring 3D spatial localization) or minimal scan time is preferred (i.e. patient study). Here, we propose a novel strategy called “d1opt” for improving the steady-state saturation transfer experiment, which, by employing a numerically optimized pre-saturation delay (d1opt), results in substantial reductions in saturation time as well as total repetition time. Theory The d1opt strategy is elucidated with an example of the forward rate constant (kf) measurement in CK reactions: PCr↔ATPγ. The ST pulse sequence usually employed consists of three segments: a pre-saturation delay (d1), frequency-selective saturation on ATPγ (d2) and 90 readout pulse followed by FID acquisition. The total repetition time (TR) approximately equals d1+d2 if acquisition time is short and negligible. The change of PCr and ATPγ magnetizations during this pulse sequence can be described by the well-known modified Block equations shown to the right, where Equations (1) and (2) apply to the d1 duration whereas Equations (3) and (4) apply during d2. The steady-state solution to Equation (3) leads to a PCr magnetization (Mss,PCr) that is independent of its initial values (MPCr(d1)) which in turn is determined by d1. Therefore, choosing an appropriate pre-saturation delay (d1opt) would allows us to achieve an MPCr(d1) very close to Mss,PCr in a short time and thus results in shortest saturation time (d2opt) as well as minimal repetition time. The calculation of d1opt and d2opt is based on numerical simulation of Equations (1)~(4) and depends upon an arbitrarily estimated kf range and accuracy level using known intrinsic T1 values. Experiment The experiments were performed at a 40-cm bore 4.7 T magnet with 4 healthy farm pigs. In vivo P MR spectra from the anterior left ventricle (LV) wall were obtained with an open-chest surgery and a surface coil sutured onto it as described previously [1]. A BISTRO pulse train was used to saturate the ATPγ followed by an adiabatic half passage pulse (AHP) for readout [2,3]. Results The d1opt and d2opt for an estimated kf range of 0.2~0.8 s (valid for most in vivo and in situ myocardial CK kinetics studies based on previous reports) and accuracy level of >95% were calculated to be 2 and 3.3 s, respectively. This is significantly superior to non-optimized approaches (both long and short d1, shown in Fig 1) in regards to saturation time and total repetition time. The experimental validation of d1opt strategy on pig heart is shown in Fig 2. There’s no detectable difference in MssPCr measurement between the d1opt and the conventional strategies whereas great reductions in saturation time (56%) and total repetition time (73%) have been granted for d1opt approach. A further quantitative validation was performed using progressive saturation experiments with different d1 values. The normalized PCr signal intensities together with simulated lines are shown in Fig 3. The excellent consistence between experiment and simulation results again confirms the validity of d1opt strategy. The kf value calculated by using the new strategy yields a value of 0.42±0.01 s which is similar to our previous reports of 0.49±0.07 s in pig hearts using the conventional method [1]. Discussions The estimated kf range is based on reports of previous in vivo or in situ studies. In general, the narrower the estimated kf range is, the shorter a saturation time as well as repetition time can be granted, whereas accompanied with greater risks and vise versa. Considering usually a very small standard deviation in kf measurement, a much narrower kf range can be confidently estimated after initial few experiments and thus further shorten the d2 and TR. In addition, detailed simulation reveals that the d1opt strategy is quite robust for variations of parameters used for calculation such as PCr/ATPγ and intrinsic T1. Another unique advantage of d1opt strategy over other strategies is the minimal saturation time. This will significantly reduce safety concerns for human studies and also benefit in reducing RF spillover effects. Using the BISTRO saturation pulse, the spillover effect was undetectable with d1opt strategy in our experiment such that control spectra usually required in conventional approach can be eliminated. Finally, d1opt strategy is readily applied to more complicated reaction models such as three-site exchange system (PCr↔ATPγ↔Pi) usually utilized in bioenergetic studies of rodent heart or brain. Reference: [1] Ye Y, et al., Circulation, 2001. [2] de Graaf RA, et al., J. Magn. Reson., 1996. [3] Lei H, et al., PNAS, 2003. Acknowledgements: This works is supported in part by NIH grants HL 50470, HL 67828, NS41262 and P41 RR08079.

ATP production rate via creatine kinase or ATP synthase in vivo: a novel superfast magnetization saturation transfer method.

Novel strategy for measuring creatine kinase reaction rate in the in vivo heart.

Novel 31 P Saturation Transfer Strategy for Measurement of Chemical Exchange Reaction Kinetics in vivo

Rapid and Simultaneous Measurement of Phosphorus Metabolite Pool Size Ratio and Reaction Kinetics of Enzymes in Vivo.

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

A technique for in vivo mapping of myocardial creatine kinase metabolism

Synergistic effect on cardiac energetics by targeting the creatine kinase system: in vivo application of high-resolution 31P-CMRS in the mouse

In vivo quantification of creatine kinase kinetics in mouse brain using 31P-MRS at 7 Tesla

Real-time in vivo imaging of extracellular ATP in the brain with a hybrid-type fluorescent sensor

Probing actin‐activated ATP turnover kinetics of human cardiac myosin II by single molecule fluorescence

Real-Time Visualization of Cytosolic and Mitochondrial ATP Dynamics in Response to Metabolic Stress in Cultured Cells

Rapid 3D-imaging of phosphocreatine recovery kinetics in the human lower leg muscles with compressed sensing.

Phosphate Metabolite Concentrations and ATP Hydrolysis Potential in Normal and Ischaemic Hearts.

Demonstration of fast and equilibrium human muscle creatine CEST imaging at 3 T

A Dual-Control Strategy by Phosphate Ions and Local Microviscosity for Tracking Adenosine Triphosphate Metabolism in Mitochondria and Cellular Activity Dynamically

Mitochondrial oxidative phosphorylation capacity in skeletal muscle measured by ultrafast Z‐spectroscopy (UFZ) MRI at 3T

A microfluidic chip-based capillary zone electrophoresis-mass spectrometry method for measuring adenosine 5′-Triphosphate and its similar nucleotide analogues

Validation of a recording spectrophotometric method for measurement of membrane-associated Mg- and NaK-ATPase activity

Determination of bioactivity of cardiac muscle sarcolemma ATPase and its application

Ex vivo and in vivo fluorescence detection and imaging of adenosine triphosphate

Rapid and Specific Imaging of Extracellular Signaling Molecule Adenosine Triphosphate with a Self-Phosphorylating DNAzyme