Fitting the blood-to-plasma or plasma-to-blood ratios
In the analysis of PET data, conversion of blood to plasma data is often needed. For the conversion, we need to know the dynamic relation between blood and plasma concentrations, and therefore the concentration in blood and plasma must be measured from blood samples taken during the PET scan. The plasma must be separated from blood, and the radioactivity and mass of both plasma and blood samples (or the precipitant left in the tube after most of plasma is removed) must be measured. Radioactivity measurement is hampered by the fast decay of radioactivity, especially with C-11 labelled radiopharmaceuticals. Fitting of a mathematical function to the measured ratios reduces variation and enables interpolation and extrapolation of the ratio curve, for the conversion of blood TAC to plasma.
The plasma-to-blood and blood-to-plasma ratio curves can be very different for different radiotracers, requiring development of a model of its own for each radiotracer, or empirically selecting a suitable function.
Erythrocytes (red blood cells, RBCs) are the main cellular component of blood, and the concentration difference between plasma and blood is dependent on partitioning of radiopharmaceutical and its labelled metabolites between water spaces in plasma and RBCs, and the rate of transport across RBC membrane. Binding to plasma proteins has major impact on the distribution, as well as binding to RBC proteins (mainly haemoglobin). RBC-to-plasma ratio (rRBC/P) as function of time determines the blood-to-plasma ratio, but concentrations can be directly measured only from plasma and blood samples, not from RBCs. RBC-to-plasma ratio can be converted to plasma-to-blood ratio with equation
, where HCT is the haematocrit. Plasma-to-blood ratio is preferred to blood-to-plasma ratio, because divider is not zero if either plasma or RBCs contain any radioactivity. Most radiotracers are administered intravenously into the plasma space, and in that case RBC-to-plasma ratio at zero time is zero, and plasma-to-blood ratio is only dependent on the haematocrit:
This is also the maximum value that the plasma-to-blood curve can get at any time point. In the opposite case, if practically all radioactivity is found in RBCs ([15O]CO, [15O]O2), rRBC/P is very high, and then rP/B approaches zero. HCT cannot be much higher than 0.5, and thus rP/B must be between 0 and ∼2.
After intravenous administration, radiotracer immediately starts to form equilibrium with the constituents of the blood, and at least partial equilibrium is achieved before bolus reaches the blood sampling site. While radiolabelled metabolites usually appear in the blood at sampling site only after one full circulation, requiring a delay time term in functions and models for fraction of metabolites in plasma, the delay time term is not needed in functions for RBC-to-plasma ratio.
Time-dependent plasma-to-blood ratio can be described using compartmental models. For instance, Lee et al (2008) proposed a simple model for [18F]FDG rat studies. Asselin et al (2002) developed a unified compartment model for 6-[18F]fluoro-L-meta-tyrosine and its metabolites in blood and partitioning between plasma and blood cells.
Measurement of rapid concentration changes in blood and plasma soon after radiotracer administration is difficult and prone to sampling and timing errors. Unfortunately the modelling of radiotracer kinetics between plasma and blood cells has reached little attention, and in most studies only simplistic empiric functions, at best, have been applied to the data.
The rRBC/P curve could be described with variable functions. Simulated rRBC/P curves based on surge function with recirculation
Program fit_pbr can fit plasma-to-blood ratio curve with this function and HCT.
Exponential or two-exponential function with constant term has sometimes been used for fitting plasma-to-blood or blood-to-plasma ratio curves (Lee et al., 2008; Wu et al., 2016; Mu et al., 2020; Galovic et al., 2021). The two-exponential function for plasma-to-blood ratio
(Mu et al., 2020) can fit similar data as the surge function with recirculation (Eq 1), but parameter limits need to be set carefully to prevent non-physiological fit results, such as initially negative concentration in blood cells, or plasma-to-blood ratio exceeding its maximum 1/(1-HCT).
Usually, at least the terms p5 and p6 can be set to 0. Examples of simple rRBC/P curves that can be fitted with rational function are shown in Figure 2 (left side), and plasma-to-blood ratios calculated from the same curves are shown on the right side of the figure.
More complicated data can be fitted with this function by freeing the parameter p5. However, rational functions can easily result into non-physiological results and the fits cannot be used to extrapolate the ratio.
Due to its many parameters and usually few samples, some of the parameters may need to be constrained to population means. Constraining parameter A1 to zero leaves function
, which too may be useful for some radiotracers.
, where p1 < 0, and t is normalised time (sample time divided by study length, thus t < 1). Function value approaches 1 at very small t, which means that concentration in blood cells at administration time is assumed to be the same as in plasma.
Somewhat similar function was used to fit blood-to-plasma ratio curves by Tarkia et al (2012):
With this function at t=0 the blood-to-plasma ratio is p4, which thus can be fixed to 1-HCT, if HCT is measured, and we assume that the concentration in blood cells is initially zero. For RBC-to-plasma ratio the function can be written as
Function parameters are saved into specific fit file format, which are ASCII text files.
- Blood to plasma conversion
- Plotting curves
- Blood sampling
- Processing input data
- Fitting PET input curves
- Fitting the fractions of unchanged radiopharmaceutical in plasma
Asselin M-C, Wahl LM, Cunningham VJ, Amano S, Nahmias C. In vivo metabolism and partitioning of 6-[18F]fluoro-L-meta-tyrosine in whole blood: a unified compartment model. Phys Med Biol. 2002; 47: 1961-1977. doi: 10.1088/0031-9155/47/11/309.
Lee J-S, Su K-H, Lin J-C, Chuang Y-T, Chueh H-S, Liu R-S, Wang S-J, Chen J-C. A novel blood-cell-two-compartment model for transferring a whole blood time activity curve to plasma in rodents. Comput Methods Programs Biomed. 2008; 92(3): 299-304. doi: 10.1016/j.cmpb.2008.02.006.
Updated at: 2023-05-09
Created at: 2016-04-06
Written by: Vesa Oikonen