Kinetic models for estimating occupancy from single-scan PET displacement studies

Abstract

The traditional design of PET target engagement studies is based on a baseline scan and one or more scans after drug administration. We here evaluate an alternative design in which the drug is administered during an on-going scan (i.e., a displacement study). This approach results both in lower radiation exposure and lower costs. Existing kinetic models assume steady state. This condition is not present during a drug displacement and consequently, our aim here was to develop kinetic models for analysing PET displacement data. We modified existing compartment models to accommodate a time-variant increase in occupancy following the pharmacological in-scan intervention. Since this implies the use of differential equations that cannot be solved analytically, we developed instead one approximate and one numerical solution. Through simulations, we show that if the occupancy is relatively high, it can be estimated without bias and with good accuracy. The models were applied to PET data from six pigs where [¹¹C]UCB-J was displaced by intravenous brivaracetam. The dose-occupancy relationship estimated from these scans showed good agreement with occupancies calculated with Lassen plot applied to baseline-block scans of two pigs. In summary, the proposed models provide a framework to determine target occupancy from a single displacement scan.

Keywords

Displacement drug occupancy kinetic modelling PET synaptic density

Introduction

The discovery and development of drugs for the treatment of brain disorders is a challenging process requiring extensive resources, long timelines and significant investments.¹ Over the past decades, imaging with positron emission tomography (PET) has become a valuable tool in CNS drug development. PET imaging with appropriate radioligands makes it possible to determine at an early stage whether a candidate drug penetrates the blood-brain barrier and binds to the target of interest in vivo. This helps to ensure that only suitable candidates will be advanced to subsequent trial phases, saving substantial resources.²

In the traditional experimental set-up for determination of target occupancy, two (or more) PET scans are acquired in each subject using a radiotracer that binds to the same target as the drug. Usually, one scan is acquired at baseline (i.e., without drug), and subsequent scan(s) after administration of the drug. The difference between outcome measures from baseline and follow-up scans is then used to determine occupancy, i.e., the fraction of targets occupied by the drug.³

The analysis of data from a PET occupancy study typically follows a three-step approach. First, a mathematical model (either based on arterial blood samples or a reference region) is used to quantify radiotracer uptake for each scan.⁴ Second, outcome measures from the different scans are combined to estimate the occupancy at the time of the post-drug scan.^3,5 Last, all subjects are pooled in a occupancy plot (sometimes referred to as the E_max model) where the administered doses or plasma concentrations of the drug are related to the measured occupancies. This final step provides information of the drug’s affinity to the target, defined as the half maximum inhibitory concentration, IC₅₀.

This established methodology leaves room for improvement. First, relying on multiple PET measurements may introduce unwanted variance in the data. It is often difficult to design the experiment in such a way that the intervention is the only difference between the scans, as other (e.g., time-related) factors may affect radiotracer uptake.^6–9 Second, the test person is exposed to multiple doses of ionizing radiation. Third, PET is a relatively expensive research tool, so a method that requires two scans for each data point places an unnecessary burden on the research budget.

An alternative approach to determine drug occupancy is to administer the drug during an on-going PET scan and based on a single injection of radioligand. When the drug is administered, competitive binding causes displacement of the radiotracer. Deriving occupancy from such a study would result in reduced costs and lower radiation exposure, as each subject would need to undergo only one scan.

Unfortunately, standard pharmacokinetic models, routinely used to analyse PET data, cannot be applied to data obtained from a displacement scan, as these models rely on the assumption of steady-state throughout a scan.¹⁰ This assumption implies that all model parameters remain constant over time, i.e., the modelled system is assumed to be time-invariant. This assumption is violated when a competing drug is administered during an on-going scan. To quantify displacement studies, a new class of pharmacokinetic models needs to be developed to incorporate the perturbation of the steady state.

The idea to model competitive binding was already pioneered in the 1990 s, when models to quantify release of endogenous dopamine were developed.^11–14 Of the models relying on bolus injection, the only established model today is the neurotransmitter PET (ntPET), which has been used to evaluate dopamine release under various conditions.^15–18 The ntPET, and variants thereof,^19–23 rely on reference tissue models, i.e., quantification is performed using a reference region rather than an arterial input function. For many radiotracers, true reference regions do not exist, due to ubiquitous expression of the target, although some degree of specific binding in the reference region may be tolerable in clinical studies. In pharmacological intervention studies, however, specific binding in the reference region is particularly problematic, as blocking in both target and reference region can result in a complicated bias in the occupancy estimates.

Here, we present a pharmacokinetic model that describes radiotracer kinetics during a displacement scan based on an arterial input function rather than a reference region. The model is based on one-tissue compartment model (1TCM) kinetics. The corresponding two-tissue compartment model (2TCM) is presented in section A of the supplementary material. Because the competing drug will perturb the system’s steady state, analytical solutions to the model equations do not exist. Instead, we present two alternative approaches: an approximate analytical solution that is derived by introducing assumptions on the accumulation of the competing drug in brain, and one numerical solution. The performance of the model and solutions are evaluated using simulations, and applied to pig [¹¹C]UCB-J PET scans24 with brivaracetam displacement, where both the radiotracer and drug were administered by bolus injection.

Materials and methods

Theory

The occupancy model

When radioligand displacement is induced by the introduction of a competing cold ligand (drug) it is assumed that the change is caused by reduction in specific radioligand binding only. This can be modelled by defining an occupancy function, $\partial (t),$ with $0 \leq \partial (t) \leq 1,$ that acts on the concentration of available binding sites, B_avail. We typically have little knowledge about the drug concentration time profile in brain tissue in vivo. To set $\partial (t)$ , we defined a set of conditions to be fulfilled. The occupancy model should be,

monotone non-decreasing

continuous and differentiable in all time points (i.e., a smooth growth)

0 at the time of drug administration

To fulfil the conditions above, we modified a model originally developed in agricultural sciences to predict crop growth rates.²⁵ Our model for the occupancy function $\partial (t)$ is

\partial (t) = \partial^{\max} (1 + \frac{t_{e} - t}{t_{e} - t_{m}}) {(\frac{t {- t}_{b}}{t_{e} - t_{b}})}^{\frac{t_{e} - t_{b}}{t_{e} - t_{m}}}, t_{b} \leq t \leq t_{e}

(1)

where

\partial^{\max}

is the maximal occupancy reached (

0 \leq \partial^{\max} \leq 1

t_{e}

is the end time of the growth (i.e., the time at which

\partial^{\max}

is reached),

t_{b}

is the begin-time of the growth (i.e., the time of the effect of intervention), and

t_{m}

is the time during which

\partial' (t)

reaches a maximum (i.e.,

t_{m}

will control the steepness of

\partial (t)

). The function

\partial (t)

will show a sigmoidal growth within the interval

t_{b} \leq t \leq t_{e}

, is exactly 0 at

{t = t}_{b}

, and can allow asymmetric or symmetric growth curves depending on the choices of

t_{e}

and

t_{m}

. To reduce the number of fitted parameters, we have constrained t_m to be the midpoint between t_b and t_e (

t_{m} = t_{b} + ((t_{e} - t_{b}) / 2

), in both simulations and solutions. Examples of

\partial (t)

are shown in Figure 1.

Figure 1.

The left panel shows the occupancy function, ∂(t) as defined in equation (1), for 3 different choices of t_e. The middle panels show schematic diagrams of the intervention models for 2TCM and 1TCM, respectively. The right panels show time-activity curves originating from each compartmental model using the choices for ∂(t) depicted in the left panel.

Displacement model

The 2TCM is the most common pharmacokinetic model used in quantification of brain PET data. In the 2TCM, the rate of exchange between compartments is determined by the constants K₁, k₂, k₃ and k₄ (Figure 1). The rate constant k₃ is linearly dependent on the concentration of available targets ( $B_{avail}$ ), $k_{3} = f_{N D} k_{o n} B_{avail}$ .^3,4 We assume that a reduction of available targets will have a negligible impact on both the association rate constant $k_{o n}$ and the fraction of free tracer in the non-displaceable compartment, $f_{N D}$ . It follows that a time dependent reduction of available targets, i.e., $(1 - \partial (t)) B_{avail}$ will affect k₃ equally, i.e., $f_{N D} k_{o n} \cdot (1 - \partial (t)) B_{avail} = {(1 - \partial (t)) k}_{3}$ . With this, the 2TCM can be modified to accommodate an increase in occupancy, starting at some time $t_{b}$ after radiotracer injection (details and equations provided in section A of the supplementary material).

The pharmacokinetics of some radiotracers are reasonably well approximated by a 1TCM. In the 1TCM, the compartments corresponding to specific and non-displaceable uptake are collapsed into a single compartment, where rate constants K₁ and k₂ describe the transfer rate of radiotracer to and from that compartment. To modify the 1TCM to accommodate displacement, we adapted the framework of the simplified reference tissue model, where the compartments for specific and non-displaceable binding are presumed to latently reside within the model configuration (Figure 1).²⁶ Setting equal the distribution volumes for the 1- and 2 tissue compartment configurations, a relationship between the apparent efflux rate constant, k_2a, and k₂–k₄ from the latently present 2TCM configuration can be derived, $k_{2 a} = k_{2} / (1 + \frac{k_{3}}{k_{4}})$ . This modification enables the introduction of an occupancy parameter to act on k₃ in the 1TCM configuration. A schematic for the model is shown in Figure 1, and the corresponding differential equation becomes

\frac{d C_{T} (t)}{d t} = {K_{1} C}_{p} (t) - \frac{k_{2}}{1 + (1 - \partial (t)) {B P}_{N D}} C_{T} (t)

(2)

where

C_{p} (t)

is the metabolite corrected arterial plasma input function,

{B P}_{N D} = k_{3} / k_{4}

, and

\partial (t)

as defined by

\partial (t) = \{\begin{array}{l} 0, t < t_{b} \\ \partial^{\max} (1 + \frac{t_{e} - t}{t_{e} - t_{m}}) {(\frac{t - t_{b}}{t_{e} - t_{b}})}^{\frac{t_{e} - t_{b}}{t_{e} - t_{m}}}, t_{b} \leq t < t_{e} \\ \partial^{\max}, t \geq t_{e} \end{array}

(3)

As indicated in equation (3), the competing drug is assumed not to wash out during the course of the scan (see Discussion).

The single step solution

In contrast to the standard kinetic models, equations (2) describe a time-variant system, and common tools for finding analytical solutions are not strictly defined. To find solutions for the differential equations (2), we placed some restrictions on the occupancy function, $\partial (t)$ .

For a drug acting rapidly on the target, we assume that $\partial (t)$ takes the form of a step function, i.e.,

\partial (t) = \{\begin{array}{l} 0, t \leq t_{s} \\ \partial^{\max}, t > t_{s} \end{array}

(4)

With this simplification, we can partition the PET time-activity curve (TAC) into two segments (i.e., before and after the time at which the drug is assumed to act on the system, t_s) and apply the 1TCM separately to each segment. We assume that the rate constants are the same for the two segments, but the differential equation for the segment after the step ( $t > t_{s}$ ) will come with non-zero initial values. The initial values for the $t > t_{s}$ segment are set to the values at the endpoint of the $t \leq t_{s}$ segment.

Let $C_{T}^{0} (t)$ and $C_{T}^{1} (t)$ denote the tissue concentrations before and after administration of the competing drug. The equations describing $C_{T}^{0} (t)$ are the standard differential equations for the 1TCM (see equation (2), with $\partial (t)$ = 0). Setting $τ = t - t_{s},$ the differential equation for $t > t_{s}$ becomes

\{\begin{array}{l} \frac{d C_{T}^{1} (τ)}{d τ} = {K_{1} C}_{p} (τ + t_{s}) - \frac{k_{2}}{1 + (1 - \partial^{\max}) {B P}_{N D}} C_{T}^{1} (τ), \\ C_{T}^{1} (0) = C_{T}^{0} (t_{s}) \end{array}

τ > 0

(5)

The solution to these differential equations becomes

C_{T} (t) = \{\begin{array}{l} K_{1} C_{p} (t) \otimes e^{- \frac{k_{2}}{1 + {B P}_{N D}} t}, t \leq t_{s} \\ K_{1} C_{p} (t) \otimes e^{- \frac{k_{2}}{1 + (1 - \partial^{\max}) {B P}_{N D}} (t - t_{s})} t > t_{s} \\ + C_{T}^{0} (t_{s}) {\cdot e}^{- \frac{k_{2}}{1 + (1 - \partial^{\max}) {B P}_{N D}} ({t - t}_{s})} \end{array}

(6)

Fitting the single step solution thus means estimating a total of 4 model parameters: K₁, k₂ and BP_ND and $\partial^{\max}$ . In our implementation, we also included the fractional blood volume (v_B), and the time at which the occupancy step occurs ( $t_{s}$ ), as free parameters. The reason to include $t_{s}$ as a free parameter is to reduce errors caused by setting $\partial (t)$ to be a step function: by allowing the model to perform the step later than the time of intervention we submit that a better description of the data can be obtained. The equations for the single step solution of the 2TCM is provided in the supplementary material (section A.2).

The numerical solution

The Euler Forward method was used to obtain a numerical solution to equation (2).

Starting at known initial conditions, each next point on the model curves is calculated by,

C (t_{n}) = C (t_{n - 1}) + h C' (t_{n - 1})

where

h = t_{n} - t_{n - 1}

(i.e., linearity is assumed in the small time interval

t_{n - 1}

< t < t_n). Insertion of the model equations in (2) gives

C_{T} (t_{n}) = h \cdot K_{1} C_{p} (t_{n - 1}) + (1 - \frac{h \cdot k_{2}}{1 + {(1 - \partial (t_{n - 1})) B P}_{N D}}) C_{T} (t_{n - 1})

The occupancy function shown in equation (3) (Figure 1) was used for $\partial (t_{n})$ , and h was set to 0.5 seconds. The corresponding numerical solution for the 2TCM is presented in the supplementary material (section A.3).

Fitting of multiple regions simultaneously

As seen in equation (2), the free model parameters $k_{2}$ , $\partial^{\max},$ and ${B P}_{N D}$ appear only as a product, and as such, each of these parameters are free to take any value as long as the term $\frac{k_{2}}{(1 + \partial (t)) {B P}_{N D}}$ results in an adequate fit. Consequently, unless some constraints are placed in the optimization, none of these parameters can be properly identified (see supplementary material, section B and supplementary figures 1 and 2).

In PET kinetic modelling, it is common to assume that the non-displaceable distribution volume (V_ND) is the same across all included brain regions. In dose-occupancy studies, it is also common to assume that the fractional occupancy is the same across the brain. In fact, these two assumptions form the basis for many quantification strategies, including all reference tissue modelling as well as the Lassen Plot.²⁷ We therefore constructed our optimizer so that $V_{N D}$ and $\partial^{\max}$ are shared across all included brain regions, whereas the other model parameters, i.e., $K_{1}, B P_{N D},$ and $v_{B}$ , are free to vary across the brain. Thus, while $V_{N D}$ , which is equal to the ratio of $K_{1}$ to $k_{2}$ , is shared across regions, $K_{1}$ and $k_{2}$ themselves are allowed to vary between regions. For the single step approach, the estimated time at which the model performs the jump, t_s, was also treated as a global parameter. Similarly for the numerical solution, the estimated time at which the growth of the occupancy ends, t_e, was estimated globally. For both methods, we used a nested approach to fit the models to the data. In an outer layer, the global parameters (V_ND, ∂^max, t_e/t_s) were estimated with non-linear least squares. For each iteration of the outer layer, the remaining model parameters (K₁, V_S, v_B) were fitted for each ROI separately. More details about models and implementation can be found in the supplementary material (section C).

Simulations

Generation of noise-free time activity curves

In simulations, we attempted to mimic the behaviour of [¹¹C]UCB-J.²⁴ For C_p(t), we used an arterial input function measured from a pig scan (baseline scan of experiment 1 in Table 1), where the measured activities after the peak were fitted with a tri-exponential function. The tissue rate constants were taken from Table 1 in (Finnema et al., 2018),²⁸ and set to result in V_ND = 4, and V_T ranging between 14.2 and 22.4. We simulated TACs for seven regions: putamen, temporal cortex, occipital cortex, frontal cortex, thalamus, cerebellum, and hippocampus.

Table 1.

Overview of pig experiments.

Experiment	Type	Dose [mg/kg]
Experiment 1	Baseline scan
	Intervention scan	0.1
	Block scan
Experiment 2	Baseline scan
	Intervention scan	2
	Block scan
Experiment 3	Intervention scan	5
Experiment 4	Intervention scan	1
Experiment 5	Intervention scan	0.75
Experiment 6	Intervention scan	0.2

To evaluate the performance of our methods, we simulated displacement TACs with two different types of drugs: one fast and one slower. In both cases, the time of injection of the drug was set at 60 minutes after radiotracer administration. For the fast-acting drug, the time t_e at which the drug reached it maximal occupancy was set to 65 minutes, for the slow-acting drug, t_e was set to 90 minutes. Scan durations were set to 150 minutes. For each of the drugs, we simulated displacement scans at three different occupancies (∂^max): 25%, 50% and 75%. For each of the six different cases (combination of drug and occupancy) we simulated 1000 unique scans, with noise added as explained in the section Generation of noise below. To generate the TACs we used the Euler Forward method and the differential equation for the 1TC displacement model (equations (2) and (3)). Some example TACs with different t_e are shown in Figure 1.

Generation of noise

To create realistic noise, a previously proposed noise-model was used that allows for time-dependent variance.^29,30

C_{noisy} (t_{k}) = C_{true} (t_{k}) \cdot (1 + α \sqrt{\frac{e^{λ t_{k}^{mid}}}{C_{true} (t_{k}) Δ t_{k}}} \cdot G (0, 1))

where

λ

is the decay constant for the isotope (in this case ¹¹C),

t_{k}^{mid}

and

Δ t_{k}

are the mid-time and duration of frame k, respectively, and

G (0, 1)

is a number sampled from a Gaussian distribution centered on 0 with a SD of 1. The scaling factor

α

was set to 5 in order to create noise on par with that of real experiments. Figures of example TACs at different noise levels can be found in the supplementary material (supplementary figures 3 and 4).

PET experiments

Experimental procedure

[¹¹C]UCB-J was largely synthesized as in Nabulsi et al., 2016,²⁴ with some modifications.³¹ All animal experiments conformed to the European Commission’s Directive 2010/63/EU and the ARRIVE guidelines. The Danish Council of Animal Ethics had approved all procedures (Journal no. 2016-15-0201-01149). Six female domestic pigs (crossbreed of Yorkshire × Duroc × Landrace, mean weight 23.3 [range: 18–27] kg) were fully anesthetized and scanned using [¹¹C]UCB-J administered as a bolus injection (injected dose: 441 [range: 344–528] MBq; injected mass: 0.69 [range: 0.08–2.77] µg). Two of the pigs (experiments 1 and 2) underwent three scans, i.e., baseline (121 min), displacement (150 min) and blocking (121 min) scans, on the same day. Brivaracetam (Briviact®, 10 mg/mL, UCB Pharma, Belgium) was administered i.v. during the displacement scan and served as a blocking agent in the third scan, which was started approximately 120 min after the brivaracetam intervention. The remaining four pigs only underwent displacement scans. In all displacement scans, brivaracetam was administered i.v. over 20 seconds, 60 min after radioligand injection. A list of all experiments is shown in Table 1.

PET data processing

PET scans were acquired on a High Resolution Research Tomograph (HRRT; Siemens, USA) and reconstructed using OP-3D-OSEM, including modelling of the point-spread function, with 16 subsets, 10 iterations and all standard corrections.³² Data was binned into the following time frames: 6 × 10, 6 × 20, 4 × 30, 9 × 60, 8 × 120, 4 × 180, 2 × 240, 1 × 300, 1 × 420, 1 × 600, 1 × 900 and 1 × 1680 s for the 121 min scans, and 6 × 10, 6 × 20, 6 × 30, 6 × 60, 4 × 120, 14 × 300, 8 × 150 s and 8 × 300 s for the 150 min scans. Attenuation correction was performed using the MAP-TR μ-map.³³ Definition of brain regions of interests (ROIs) was performed using a dedicated pig brain template.³⁴ The seven regions from the simulation experiment were also used here: putamen, temporal cortex, occipital cortex, frontal cortex, thalamus, cerebellum and hippocampus.

Blood and plasma analyses

Radioactivity in arterial whole blood was measured continuously for the first 30 min of each scan using an Allogg ABSS autosampler (Allogg Technology, Sweden). Arterial blood was manually drawn at 3, 8, 10, 15, 30, 45, 59, 61, 65, 75, 90, 105, 120 and 150 min for measuring radioactivity in whole blood and plasma using a gamma counter (Cobra, 5003, Packard Instruments, Meriden, USA) that was cross-calibrated against the HRRT. Radio-HPLC was used to measure radioligand parent fractions.³⁵ A more detailed account of the blood and plasma analyses can be found in the supplementary material (Section D).

In the baseline scan of experiment 2, (see Table 1), the parent fraction could not be estimated due to a technical failure. The parent fractions from the displacement and blocking scans conducted in the same animal were, however, very similar (absolute difference averaged across time was 6.0 ± 4.1%, and difference in AUC was 2%). Therefore, for the baseline scan in experiment 2, the mean parent fraction from the corresponding displacement and blocking scans was used.

The concentration of brivaracetam in arterial plasma was analysed using UPLC-MS/MS (Filadelfia Epilepsy Hospital, Denmark). During the displacement scans seven blood samples (at approximately 1, 5, 15, 30, 45, 60 and 90 min after brivaracetam injection) were collected for this purpose, and during the block scans five samples (at approximately 3, 15, 45, 75 and 90 min after scan start) were collected.

Results

Simulation results

Figure 2 and Table 2 summarise results from the simulation experiment. They show that occupancy estimation improves both with increasing drug speed and with increasing dose. The performance of the two methods (numerical solution and single-step approximation) were comparable throughout, especially for the higher occupancies, where the histograms are almost identical. While the occupancy estimates are approximately normally distributed for the higher occupancies, the distributions of estimates are slightly skewed for both methods at 25% occupancy. At this lower occupancy, there is also a bigger difference between the distributions, with the single-step solution, unlike the numerical solution, showing a slight tendency to overestimate ∂^max (for the fast drug, median ∂^max estimates were 25.2% with the numerical solution and 30.2% with the single step solution). Corresponding results for V_ND and V_S are found in supplementary figures 8 and 9, respectively.

Figure 2.

Results from simulating fast- (t_e = 5 min) and slow-acting (t_e = 30 min) drugs, displacing [¹¹C]UCB-J binding. Each panel shows histograms of occupancy (∂^max) estimates from the numerical solution (green) and single-step approximation (red). In each panel, the dashed black line corresponds to the true value for ∂^max.

Table 2.

Summary of simulation study results.

True values		V_ND		∂^max [%]
∂^max [%]	t_e–t_b [min]	Numerical solution	Single step solution	Numerical solution	Single step solution
25	5	4.05 (4.91)	6.16 (3.50)	25.2 (8.9)	30.2 (9.1)
	30	3.73 (6.31)	5.35 (2.41)	24.7 (9.9)	31.3 (8.7)
	120	2.38 (9.62)	4.93 (9.69)	37.9 (86.1)	13.4 (11.5)
50	5	3.94 (1.64)	3.94 (1.65)	49.9 (6.0)	50.0 (6.0)
	30	4.05 (1.98)	3.97 (1.94)	50.2 (7.6)	49.7 (7.3)
	120	3.40 (7.12)	3.77 (4.54)	59.5 (65.9)	25.8 (12.5)
75	5	3.99 (0.63)	3.97 (0.63)	75.0 (3.4)	74.9 (3.5)
	30	4.03 (0.81)	3.89 (0.79)	75.2 (4.4)	74.1 (4.2)
	120	3.55 (4.90)	3.77 (4.54)	74.3 (44.9)	39.9 (12.0)

All parameter estimates are given as median with interquartile ranges (difference between 75th and 25th percentile) in parenthesis. Occupancy ∂^max estimates are given in percentage. t_e–t_b is the difference between the time of the intervention and the time when ∂^max is reached for the simulated drug and is given in minutes. For each condition, 1000 different noise instances were simulated. In all simulations, the true value of V_ND was 4.

Displacement models applied to real data

Model fits

The 1TC displacement model was consistently able to describe the measured TACs using both the single-step approximation and the numerical solution. Figure 3 shows model fits to temporal cortex TACs, with both methods, for the largest and the lowest dose experiments. For the numerical solution, fits to all TACs, with residuals, for the same two scans can be found in supplementary figures 13 and 14 and normalized residuals for all six pig scans can be found in supplementary figure 15. The average ± SD total distribution volume in temporal cortex, across all displacement scans was 20.2 ± 3.7 mL/cm³ for both methods. The occupancies ranged from 41% to 86% (see Table 3).

Figure 3.

Displacement model fits (solid lines) to [¹¹C]UCB-J temporal cortex TACs (dots) from two pig scans in which brivaracetam was administered i.v. 60 minutes after radiotracer injection. The dashed lines show model curves in the absence of displacement. These curves were generated from the estimated model parameters, with occupancy set to zero. For the numerical solution, the model fits to all regions, including residuals, are presented in supplementary figures 13 and 14.

Table 3.

Summary of results from the pig experiments.

Pig experiment	Displacement experiment				Baseline-block experiment
	Numerical solution		Single step solution		Baseline-block experiment
	V_ND	∂^max	V_ND	∂^max	V_ND	Δ
1	7.47	40.7%	7.47	40.7%	7.46	16.2%
2	1.85	79.4%	2.15	80.4%	2.07	65.5%
3	2.41	86.4%	2.41	81.6%	–	–
4	3.96	81.6%	3.96	81.6%	–	–
5	3.44	75.1%	3.44	75.1%	–	–
6	5.61	86.2%	5.61	86.2%	–	–

Pig numbers correspond to those in Table 1. Results for the baseline-block experiment are calculated using the Lassen plot. Note that because the block scans were conducted after the displacement scans, the Lassen occupancies (Δ) correspond to a lower concentration of brivaracetam in plasma compared to the displacement scan occupancies (∂^max).

Comparison with Lassen plots

Experiment 1 and 2 each had a baseline scan before and a blocking scan after the displacement scan. These scans were analysed using the traditional 1TCM, and occupancies were estimated using the Lassen plot.^5,27 For both experiments, the Lassen occupancies in the block scans were lower than the estimates from the displacement scans (see Table 3). However, the plasma drug concentrations were also lower during the block scans, and the outcome of the Lassen-plot fit well with the dose-occupancy relationship estimated from the displacement scans (see Figure 4). The Lassen-V_NDs also showed good agreement with the ones calculated with the displacement model (see Table 3). For the high-dose pig scan, the Lassen plot returned a V_ND of 2.08, while the displacement model returned a V_ND of 1.85 with the numerical solution and 2.15 with the single step solution. For the low-dose pig scan, the Lassen-V_ND was 7.46, while the displacement model returned V_ND estimates of 7.47 with both solutions.

Figure 4.

Results from pig experiment. The two tiles show the estimated occupancies plotted against peak plasma brivaracetam concentrations for the numerical solution and single step solutions, respectively. In both, the Lassen-occupancies plotted against the peak plasma value during the block scans are also included.

Dose-occupancy model

Figure 4 shows occupancies ( $\partial^{\max})$ estimated from the displacement models for all pigs, plotted against the maximal plasma level of brivaracetam following its injection (C_briva). These occupancies could be well described by the E_max model, $\partial^{\max} = Δ_{\max} C_{briva} / (C_{briva} + {I C}_{50})$ , where IC₅₀ is the drug’s half maximal inhibitory concentration, and $Δ_{\max}$ is the maximal attainable occupancy for the population. The estimated values for $Δ_{\max}$ were 87.3% for the numerical solution, and 87.6% for the single step approach. The corresponding values for IC₅₀ were 1.26 µg/mL and 1.27 µg/mL for the numerical solution and single step, respectively.

Discussion

In this study, we present a pharmacokinetic model capable of describing PET time-activity curves after a pharmacological intervention. We have developed a generic and flexible model that allows for increasing occupancy, and that is incorporated into the standard PET compartmental models, to describe a displacement of the radiotracer during the scan. Because the differential equations are time-variant, we present two new approaches for quantification of PET data with arterial input functions. In the single step solution, the effect of the drug intervention is approximated to be instantaneous, and the system can thus be assumed to be time-invariant both before and after the effect of the intervention occurs. The “extended simplified reference tissue model” (ESRTM) is based on the same idea, although it relies on a reference tissue rather than an arterial input function.³⁶ The other solution is to solve the model differential equations using numerical methods. Here, we used the Euler Forward method together with a monotone, continuous and differentiable occupancy model.

The present results suggest that performing displacement scans is a viable alternative to the traditional two-scan setup to determine target engagement. We demonstrate the usefulness of the methods by displacing [¹¹C]UCB-J with brivaracetam in pigs. Experiments 1 and 2 showed that the displacement model resulted in consistent occupancy estimates to those obtained using the traditional Lassen plot. The estimated occupancies could be well described by the E_max model, which resulted in IC₅₀ estimates of 1.26 μg/mL for the single step method and 1.27 µg/mL for the numerical solution. The E_max model assumes that the drug concentration in plasma is constant but plasma brivaracetam changes rapidly following intravenous injection, and it is not entirely clear how to best map occupancies to such dynamic plasma levels. Because we saw a very rapid displacement of [¹¹C]UCB-J, we used peak plasma values, i.e. plasma values immediately following injection, as they align temporally to radiotracer displacement. In baseline-block experiments, where the drug is normally administered some time before the block scan, the plasma drug concentration will remain relatively constant during the scan. In these cases, the plasma concentration either at the start of the scan, the end of the scan, or the mean of the two, is often used in the E_max model.³⁷ For comparison, we re-ran the dose-occupancy analysis using mean plasma concentrations during each scan in place of the peak plasma values. This resulted in an IC₅₀ of 0.47 µg/mL, for both the numerical solution and single-step approximation. This is nearly identical to the brivaracetam IC₅₀ reported by Finnema and colleagues (0.46 µg/mL) from a [¹¹C]UCB-J baseline-block experiment in humans.³⁸

When performing drug development studies, it is common practice to use a range of doses to better characterise the dose-occupancy relationship. The displacement models presented here do not necessarily provide good estimates in the low occupancy ranges (∼25%, or lower). At 25% occupancy, there is also a large uncertainty in the V_ND estimate (supplementary figure 8), and in several cases (especially for the numerical solution) it hits the lower bound at V_ND = 0. Due to a strong positive correlation between ∂^max and V_ND (see supplementary tables 1 and 2, and supplementary figures 11 and 12), this leads to the apparent negative bias in occupancy that we see in Figure 2. Difficulties in determining low occupancies has also been reported with a baseline-block setup.^39–41 A possible solution is to fit multiple subjects simultaneously in a multilevel pharmacokinetic modelling framework, allowing the model to differentiate between displacement and normal scans.⁴² This could improve the occupancy estimates, even if normal and displacement scans were conducted in different research subjects. Such an approach could be particularly valuable if the displacement is small, e.g., when using a behavioural task to elicit neurotransmitter release rather than a pharmacological challenge.

The assumption of instantaneous occupancy (single step) that we have employed to allow the model to be solved analytically has already been shown to be a useful one for reference region quantification of displacement scans.³⁶ We emphasize that although the single step solution involves splitting TACs into two segments, each segment is not fitted independently. All rate constants are constrained to be constant throughout the scan, and they are estimated by fitting the entire TAC.

The objective of introducing a numerical solution, accounting for the time-course of occupancy, was to allow better quantification of occupancy for slow-acting drugs. Unexpectedly, the two approaches performed comparably across all experiments, even for the simulated slower drug (t_e = 30 min). In addition to the presented data, we simulated scans with a much slower drug (t_e = 120 min, see supplementary figure 5), where the final occupancy was reached 30 min after the last acquired data. Even in this case, we saw no advantage of the numerical solution over the single step simplification. In fact, both approaches performed poorly in this scenario. The poor performance for t_e = 120 min presents a limitation on the possible applications of our model, and likely even the displacement scan experimental setup itself. It is clear that the model will only be applicable to those drugs and modes of administration where ∂^max can be expected to be reached relatively fast.

The performance of the single step solution presented in this paper relies on using a relatively high time resolution (0.5 s frequency) in the convolution step. In our experience, this is a much shorter step size than what is usually used when solving compartment models. Consequently, the single step approach is computationally relatively heavy, and requires approximately 30 times longer run-time than the Euler Forward-based numerical solution.

Another advantage of the numerical approach is that it allows flexibility. In this study, we have used a monotonically increasing function to explain the time course of occupancy. These assumptions appear to be reasonable for the [¹¹C]UCB-J pig scans with brivaracetam intervention. However, depending on factors like the drug, radiotracer, and experimental design, in some cases it might be preferable to use a different type of function to describe the occupancy. For instance, the numerical approach allows for an occupancy model where both drug uptake and washout happen during the scan. We chose the occupancy model in (equation (1)) because it is a continuous and differentiable function that allows some key parameters to be estimated. In reality, we expect that the increase in occupancy in the time following a drug intervention is at first rapid, and then slows down as fewer binding sites remain available. In a recent study, Naganawa and colleagues present a similar displacement model for [¹¹C]UCB-J. In that study, the rate constants defining the drug’s uptake and clearance in tissue were estimated, together with the radioligand’s rate constants.⁴³ While this model more accurately reflects the underlying competition at the SV2A binding sites, the authors show that parameter identifiability becomes challenging with a model of that complexity. The approach presented in the current study is thus a pragmatic solution to derive occupancy estimates from a single displacement scan.

We are confident that the simulated TACs have a level of noise that is realistic to [¹¹C]UCB-J. The level of noise could however vary for different radiotracers. In supplementary figures 6 and 7 we show histograms of ∂^max estimates at different levels of noise. With increasing noise, the precision of the parameter estimation is reduced, but the estimates remain unbiased. With no noise added, both solutions to the model consistently return the true ∂^max value.

In both solutions to the model, the time of the intervention (t_b) is treated as a known parameter and is not fitted. In the presented results, the models were solved with the true t_b values. In reality, it might be difficult to identify the exact moment when the drug reaches its target. We therefore applied the models to some of the simulated data with wrong values for t_b (1 and 5 minutes before and after the true t_b). For both solutions to the model, the ∂^max and V_ND were generally unaffected by the different values for t_b (supplementary figure 16).

A limitation of the models is that they only consider a change in available binding sites. Pharmacological interventions may also affect perfusion, which could influence some model parameters. Further work is needed to develop models that can account for other changes than a reduction in B_avail induced by the pharmacological challenge, like some of the existing reference region based methods do.^14,19 Also, similar to available methods for baseline-block scans, the model does not account for specific binding of the radiotracer to the target, assuming that it is only present in tracer doses. Depending on the specific activity of the radiotracer, this could lead to some bias.

Although we derived displacement versions of both the 1TCM and 2TCM (see supplementary material, section A.), we only considered a tracer that can be described by 1TC kinetics. Additional work is needed to evaluate the performance of the 2TC displacement models, as well as reference tissue implementations.

A limitation of the simulation experiments is that, for the numerical solution, the same model is used both to simulate the data and solve it. This could offer an unfair advantage to the numerical solution over the single step approximation, but our pig experiments confirm that the two approaches perform well, and they are in agreement with the Lassen plot outcome. We also limited our case to a drug that after intravenous injection shows a very immediate interaction with the target. Future studies must show if the two methods perform equally well for more slow-acting drugs. For solving the proposed displacement model, both with the numerical approach and the analytical approximation, it is necessary to pool data from several brain regions. This is standard for methods of estimating occupancies in the absence of a reference region.^5,27,39–41 Yet, observed data does not always align with these assumptions, and spatially varying values for both V_ND and occupancy have been reported.^44–46 When using the displacement model it is therefore important to carefully consider which regions are appropriate to group together.

In conclusion, drug displacement PET scans constitute a promising alternative to determine occupancy, compared to baseline and follow-up studies. The kinetic models presented here enable estimation of occupancy from a single displacement scan, thereby obviating the need for two consecutive scans. This allows the number of scans required for target engagement studies to be substantially reduced, leading to lower radiation exposure and experimental costs, while also limiting the variation of biological and experimental factors. To facilitate the implementation of these models in other research centres, the MATLAB code is freely available for download at https://github.com/Gjertrud/ISI.

Supplemental Material

sj-pdf-1-jcb-10.1177_0271678X231168591 - Supplemental material for Kinetic models for estimating occupancy from single-scan PET displacement studies

Supplemental material, sj-pdf-1-jcb-10.1177_0271678X231168591 for Kinetic models for estimating occupancy from single-scan PET displacement studies by Gjertrud Louise Laurell, Pontus Plavén-Sigray, Annette Johansen, Nakul Ravi Raval, Arafat Nasser, Clara Aabye Madsen, Jacob Madsen, Hanne Demant Hansen, Lene Lundgaard Donovan, Gitte Moos Knudsen, Adriaan A Lammertsma, R Todd Ogden, Claus Svarer and Martin Schain in Journal of Cerebral Blood Flow & Metabolism

Footnotes

Funding

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This project has received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie grant agreement No 796759. GLL was supported by Rigshospitalet’s Research Council (R197-A8956). PPS was supported by the Swedish Society for Medicine and the Lundbeck Foundation (R303-2018-3263). AJ was supported by The Independent Research Fund Denmark. NRR was funded by the European Union’s Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie grant agreement No 813528. HDH was supported by the Lundbeck Foundation (R293-2018-738). RTO and AAL were supported by visiting professorships from the Lundbeck Foundation (R317-2019-787 and R290-2018-756).

Acknowledgements

The authors would like to thank the veterinarians and staff at the Department of Experimental Medicine, University of Copenhagen for their continued assistance with animal experiments.

Declaration of conflicting interests

The author(s) declared the following potential conflicts of interest with respect to the research, authorship, and/or publication of this article: GMK has received honoraria as a speaker for Sage Biogen and H. Lundbeck, and as a consultant for Sanos. MS has received compensation from Roche as a key opinion leader, and is an employee and owns stock options in Antaros Medical AB.

Authors’ contributions

MS, PPS, AAL, RTO, GMK and CS developed the theoretical framework. GLL and MS implemented the models, designed the simulations, and performed all calculations. AJ, NRR, CAM, LLD and HDH planned and conducted all PET experiments. AN analyzed all blood and plasma samples. JM led all radioligand syntheses. GLL and MS drafted the manuscript, and all authors provided critical feedback and helped shape the final version.

ORCID iDs

Annette Johansen

Nakul Ravi Raval

Clara Aabye Madsen

Hanne Demant Hansen

Claus Svarer

Supplemental material

Supplemental material for this article is available online.

References

Wouters

McKee

Luyten

Estimated research and development investment needed to bring a new medicine to market, 2009–2018. JAMA 2020; 323: 844–853.

Gunn

Rabiner

EA.

Imaging in central nervous system drug discovery. Semin Nucl Med 2017; 47: 89–98.

Innis

Cunningham

Delforge

, et al. Consensus nomenclature for in vivo imaging of reversibly binding radioligands. J Cereb Blood Flow Metab 2007; 27: 1533–1539.

Slifstein

Laruelle

Models and methods for derivation of in vivo neuroreceptor parameters with PET and SPECT reversible radiotracers. Nucl Med Biol 2001; 28: 595–608.

Cunningham

Rabiner

Slifstein

, et al. Measuring drug occupancy in the absence of a reference region: the lassen plot re-visited. J Cereb Blood Flow Metab 2010; 30: 46–50.

Steinberg

Rubin-Falcone

Galfalvy

, et al. Cortisol stress response and in vivo PET imaging of human brain serotonin 1A receptor binding. Int J Neuropsychopharmacol 2019; 22: 329–338.

da Cunha-Bang

Ettrup

Mc Mahon

, et al. Measuring endogenous changes in serotonergic neurotransmission with [(11)C]cimbi-36 positron emission tomography in humans. Transl Psychiatry 2019; 9: 134.

Mc Mahon

Andersen

Madsen

, et al. Seasonal difference in brain serotonin transporter binding predicts symptom severity in patients with seasonal affective disorder. Brain 2016; 139: 1605–1614.

Matheson

Schain

Almeida

, et al. Diurnal and seasonal variation of the brain serotonin system in healthy male subjects. NeuroImage 2015; 112: 225–231.

10.

Morris

Endres

Schmidt

, et al. Kinetic modeling in positron emission tomography. In: Wernick MN, Aarsvold JN (eds) Emission tomography: The Fundamentals of PET and SPECT. San Diego, California, USA: Elsevier Academic Press, 2004, pp.499–540.

11.

Endres

Kolachana

Saunders

, et al. Kinetic modeling of [11C]raclopride: combined PET-microdialysis studies. J Cereb Blood Flow Metab 1997; 17: 932–942.

12.

Endres

Carson

RE.

Assessment of dynamic neurotransmitter changes with bolus or infusion delivery of neuroreceptor ligands. J Cereb Blood Flow Metab 1998; 18: 1196–1210.

13.

Morris

Fisher

Alpert

, et al. In vivo imaging of neuromodulation using positron emission tomography: optimal ligand characteristics and task length for detection of activation. Hum Brain Mapp 1995; 3: 35–55.

14.

Alpert

Badgaiyan

Livni

, et al. A novel method for noninvasive detection of neuromodulatory changes in specific neurotransmitter systems. Neuroimage 2003; 19: 1049–1060.

15.

Yoder

Wang

Morris

ED.

Change in binding potential as a quantitative index of neurotransmitter release is highly sensitive to relative timing and kinetics of the tracer and the endogenous ligand. J Nucl Med 2004; 45: 903–911.

16.

Morris

Yoder

Wang

, et al. ntPET: a new application of PET imaging for characterizing the kinetics of endogenous neurotransmitter release. Mol Imaging 2005; 4: 473–489.

17.

Morris

Normandin

Schiffer

WK.

Initial comparison of ntPET with microdialysis measurements of methamphetamine-induced dopamine release in rats: support for estimation of dopamine curves from PET data. Mol Imaging Biol 2008; 10: 67–73.

18.

Normandin

Morris

ED.

Temporal resolution of ntPET using either arterial or reference region-derived plasma input functions. Conf Proc IEEE Eng Med Biol Soc 2006; 2006: 2005–2008.

19.

Normandin

Schiffer

Morris

ED.

A linear model for estimation of neurotransmitter response profiles from dynamic PET data. Neuroimage 2012; 59: 2689–2699.

20.

Kim

Sullivan

Wang

, et al. Voxelwise lp-ntPET for detecting localized, transient dopamine release of unknown timing: sensitivity analysis and application to cigarette smoking in the PET scanner. Hum Brain Mapp 2014; 35: 4876–4891.

21.

Wang

Kim

Cosgrove

, et al. A framework for designing dynamic lp-ntPET studies to maximize the sensitivity to transient neurotransmitter responses to drugs: application to dopamine and smoking. Neuroimage 2017; 146: 701–714.

22.

Liu

Zakiniaeiz

Cosgrove

, et al. Toward whole-brain dopamine movies: a critical review of PET imaging of dopamine transmission in the striatum and cortex. Brain Imaging Behav 2019; 13: 314–322.

23.

Irace

Mérida

Redouté

, et al. Bayesian estimation of the ntPET model in single-scan competition PET studies. Front Physiol 2020; 11: 498.

24.

Nabulsi

Mercier

Holden

, et al. Synthesis and preclinical evaluation of 11C-UCB-J as a PET tracer for imaging the synaptic vesicle glycoprotein 2A in the brain. J Nucl Med 2016; 57: 777–784.

25.

Yin

Goudriaan

Lantinga

, et al. A flexible sigmoid function of determinate growth. Ann Bot 2003; 91: 361–371.

26.

Lammertsma

Hume

SP.

Simplified reference tissue model for PET receptor studies. Neuroimage 1996; 4: 153–158.

27.

Lassen

Bartenstein

Lammertsma

, et al. Benzodiazepine receptor quantification in vivo in humans using [11C]flumazenil and PET: application of the steady-state principle. J Cereb Blood Flow Metab 1995; 15: 152–165.

28.

Finnema

Nabulsi

Mercier

, et al. Kinetic evaluation and test-retest reproducibility of [(11)C]UCB-J, a novel radioligand for positron emission tomography imaging of synaptic vesicle glycoprotein 2A in humans. J Cereb Blood Flow Metab 2018; 38: 2041–2052.

29.

Logan

Fowler

Volkow

, et al. A strategy for removing the bias in the graphical analysis method. J Cereb Blood Flow Metab 2001; 21: 307–320.

30.

Varga

Szabo

Modified regression model for the Logan plot. J Cereb Blood Flow Metab 2002; 22: 240–244.

31.

Raval

Gudmundsen

Juhl

, et al. Synaptic density and neuronal metabolic function measured by positron emission tomography in the unilateral 6-OHDA rat model of Parkinson’s disease. Front Synaptic Neurosci 2021; 13: 715811.

32.

Sureau

Reader

Comtat

, et al. Impact of image-space resolution modeling for studies with the high-resolution research tomograph. J Nucl Med 2008; 49: 1000–1008.

33.

Keller

Svarer

Sibomana

Attenuation correction for the HRRT PET-scanner using transmission scatter correction and total variation regularization. IEEE Trans Med Imaging 2013; 32: 1611–1621.

34.

Villadsen

Hansen

Jørgensen

, et al. Automatic delineation of brain regions on MRI and PET images from the pig. J Neurosci Methods 2018; 294: 51–58.

35.

Gillings

A restricted access material for rapid analysis of [(11)C]-labeled radiopharmaceuticals and their metabolites in plasma. Nucl Med Biol 2009; 36: 961–965.

36.

Zhou

Chen

M-K

Endres

, et al. An extended simplified reference tissue model for the quantification of dynamic PET with amphetamine challenge. Neuroimage 2006; 33: 550–563.

37.

Takano

Varrone

Gulyás

, et al. Guidelines to PET measurements of the target occupancy in the brain for drug development. Eur J Nucl Med Mol Imaging 2016; 43: 2255–2262.

38.

Finnema

Rossano

Naganawa

, et al. A single-center, open-label positron emission tomography study to evaluate brivaracetam and levetiracetam synaptic vesicle glycoprotein 2A binding in healthy volunteers. Epilepsia 2019; 60: 958–967.

39.

Schain

Zanderigo

Ogden

RT.

Likelihood estimation of drug occupancy for brain PET studies. NeuroImage 2018; 178: 255–265.

40.

Naganawa

Gallezot

J-D

Rossano

, et al. Quantitative PET imaging in drug development: estimation of target occupancy. Bull Math Biol 2019; 81: 3508–3541.

41.

Laurell

Plavén-Sigray

Svarer

, et al. Designing drug occupancy studies with PET neuroimaging: Sample size, occupancy ranges and analytical methods. NeuroImage 2022; 263: 119620.

42.

Chen

Goldsmith

Ogden

RT.

Nonlinear mixed-effects models for PET data. IEEE Trans Biomed Eng 2019; 66: 881–891.

43.

Naganawa

Gallezot

J-D

Finnema

, et al. Drug characteristics derived from kinetic modeling: combined 11C-UCB-J human PET imaging with levetiracetam and brivaracetam occupancy of SV2A. EJNMMI Res 2022; 12: 71.

44.

Svensson

Schain

Plavén-Sigray

, et al. Validity and reliability of extrastriatal [11C]raclopride binding quantification in the living human brain. NeuroImage 2019; 202: 116143.

45.

de Laat

Morris

ED.

A local-neighborhood Lassen plot filter for creating occupancy and non-displaceable binding images. J Cereb Blood Flow Metab 2021; 41: 1379–1389.

46.

Rossano

Toyonaga

Finnema

, et al. Assessment of a white matter reference region for (11)C-UCB-J PET quantification. J Cereb Blood Flow Metab 2020; 40: 1890–1901.

Supplementary Material

Please find the following supplemental material available below.

For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.

For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.

0.00 MB

1.43 MB