To characterize intracellular energy transfer in the heart, two organ-level methods have frequently been employed: xc2xb3xc2xb9P xe2x88x92 NMR inversion and saturation transfer, and dynamic xc2xb9xe2x81xb8O labeling. Creatine kinase (CK) fluxes obtained by following oxygen labeling have been considerably smaller than the fluxes determined by xc2xb3xc2xb9P xe2x88x92 NMR saturation transfer. It has been proposed that dynamic xc2xb9xe2x81xb8O labeling determines net flux through CK shuttle, whereas xc2xb3xc2xb9P xe2x88x92 NMR saturation transfer measures total unidirectional flux. However, to our knowledge, no sensitivity analysis of flux determination by oxygen labeling has been performed, limiting our ability to compare flux distributions predicted by different methods. Here we analyze oxygen labeling in a physiological heart phosphotransfer network with active CK and adenylate kinase (AdK) shuttles and establish which fluxes determine the labeling state. A mathematical model consisting of a system of ordinary differential equations was composed describing xc2xb9xe2x81xb8O enrichment in each phosphoryl group and inorganic phosphate. By varying flux distributions in the model and calculating the labeling, we analyzed labeling sensitivity to different fluxes in the heart. We observed that the labeling state is predominantly sensitive to total unidirectional CK and AdK fluxes and not to net fluxes. We conclude that measuring dynamic incorporation of xc2xb9xe2x81xb8O into the high-energy phosphotransfer network in heart does not permit unambiguous determination of energetic fluxes with a higher magnitude than the ATP synthase rate when the bidirectionality of fluxes is taken into account. Our analysis suggests that the flux distributions obtained using dynamic xc2xb9xe2x81xb8O labeling, after removing the net flux assumption, are comparable with those from xc2xb3xc2xb9P xe2x88x92 NMR inversion and saturation transfer.