Quantification of human brain PDE4 occupancy by GSK356278: A [ 11 C](R)-rolipram PET study

Abstract

We characterized the relationship between the plasma concentration of the phospodiesterase (PDE)-4 inhibitor GSK356278 and occupancy of the PDE4 enzyme in the brain of healthy volunteers, using the positron emission tomography (PET) tracer [¹¹C](R)-rolipram. To this end, PET scans were acquired in eight male volunteers before and at 3 and 8 h after a single 14 mg oral dose of GSK356278. A metabolite-corrected arterial input function was used in conjunction with the dynamic PET emission data to estimate volumes of distribution (V_T) from a two-tissue compartment model. The administration of GSK356278 reduced [¹¹C](R)-rolipram whole brain V_T by 17% at 3 h post-dose (p = 0.01) and by 4% at 8 h post-dose. The mean plasma C_max was 42.3 ng/ml, leading to a PDE4 occupancy of 48% at T_max. The in vivo affinity of GSK356278 was estimated as EC₅₀ = 46 ± 3.6 ng/ml. We present the first report of a direct estimation of PDE4 blockade in the living human brain. In vivo affinity of GSK356278 for the PDE4, estimated in this early phase study, was combined with GSK356278 safety and tolerability data to decide on a therapeutic dose for future clinical development.

Introduction

Rolipram is a selective inhibitor of the enzyme phospodiesterase (PDE)-4, a member of the PDE enzyme family, which hydrolyses the second messenger cyclic adenosine monophosphate (cAMP). The carbon-11-labelled (R) enantiomer of rolipram has been demonstrated to be suitable for the in vivo evaluation of PDE4 availability and activity with positron emission tomography (PET).^1,2 [¹¹C](R)-rolipram binding is sensitive to pharmacological challenges in rat ³ and porcine brain. ⁴ However, to this point, no human blocking studies with PDE4 modulators have been published. GSK356278 is a potent, CNS penetrant inhibitor of cAMP hydrolytic activity ⁵ that binds to the high-affinity rolipram binding site (HARBS) with a pIC₅₀ of 8.6 and is nonselective for the PDE4A, B and D isoforms. ⁶ In a model of anxiety in common marmosets, the therapeutic index for GSK356278 was >10 versus <1 for rolipram, ⁶ which may be explained by the faster HARBS dissociation rate compared to rolipram. The present study was designed to evaluate the relationship between the plasma GSK356278 concentration and the occupancy of the brain PDE4 enzyme in healthy male subjects. These data were to be utilized to optimize the dose range for future Phase 1 and 2 studies.

Methods

This was an open-label study in eight healthy male volunteers with a mean (±standard deviation, SD) age of 34.4 ± 10.7. This study was approved by the Edinburgh Independent Ethics Committee for Medical Research and the Administration of Radioactive Substances Advisory Committee (ARSAC) and conducted in accordance with ICH GCP and the 2008 revision of the Declaration of Helsinki. PET scanning was performed at Imanova, Centre for Imaging Sciences (Hammersmith Hospital, London, UK). All subjects gave written informed consent, tested negative for drugs in urine, and were free from clinically significant illness or disease as determined by their medical history, a full physical examination, blood laboratory tests and electrocardiography.

Radiopharmaceutical preparation

[¹¹C](R)-rolipram was prepared by O-[¹¹C] methylation of its desmethyl analog (1 mg) in dimethylformamide (350 μl) using [¹¹C]methyl iodide in the presence of tetrabutylammonium hydroxide (0.5 M, 4 μl). The reaction mixture was heated to 85℃ for 5 min, diluted with high-performance liquid chromatography (HPLC) mobile phase (1 ml) and loaded onto semi-preparative HPLC for purification (Agilent Eclipse XDB C18, 5 µm, 250 × 9.4 mm). The retention time of [¹¹C](R)-rolipram was 5.4 min through isocratic elution with a mixture of 10 mM ammonium formate buffer/acetonitrile (60/40, v/v) at a flow rate of 6 ml/min. The desired fraction was collected, diluted with water (20 ml) and loaded onto a SepPak® Classic C18 cartridge for reformulation. Following an initial wash with sterile water (10 ml), [¹¹C](R)-rolipram was eluted off the cartridge with ethanol (1 ml) and subsequently diluted with saline (10 ml). In a final step, the resulting formulation solution was filtered through a 0.22 µm sterile filter into a sterile, pyrogen-free vial. Typical synthesis from a 55 μA, 30 min beam provided 3.35 ± 1.4 GBq (uncorrected) of [¹¹C](R)-rolipram in a total synthesis time of 35 min with a specific activity of 142 ± 133 GBq/µmol. Quality control showed that doses were obtained with radiochemical purities >99% and that the final product for injection was sterile and pyrogen free.

Positron emission tomography data acquisition

In total, 22 [¹¹C](R)-rolipram PET scans were acquired on a Siemens Biograph HiRez XVI PET-CT scanner (Siemens Healthcare, Erlangen, Germany) of which eight scans were acquired at baseline (PET 1) and seven scans at approximately 3 and 8 h post-dose each (PET 2 and PET 3 respectively). A CT scan of the head was acquired for attenuation and model-based scatter correction. Subjects were injected with an intravenous bolus of the radioligand, and dynamic emission data were acquired continuously for 90 min. PET data were reconstructed using filtered backprojection with corrections for attenuation, dead time, random coincidences and scatter. Dynamic data were binned into 26 frames (8 × 15 s, 3 × 1, 5 × 2, 5 × 5 and 5 × 10 min). A continuous sampling system (ABSS Allogg, Mariefred, Sweden) was used to measure whole blood activity for the first 15 min (sampled at 5 ml/min). Discrete blood samples were acquired to measure blood and plasma radioactivity concentration at 5, 10, 15, 20, 30, 40, 50, 70 and 90 min and to determine with HPLC the fraction of radioactivity corresponding to intact parent [¹¹C](R)-rolipram in arterial plasma. Plasmafree fraction (f_P) was determined by ultrafiltration. The three discrete blood samples at 5, 10 and 15 min post-injection were used to calibrate the continuous blood data. The continuous and discrete datasets were used to form a whole-blood activity curve, covering the duration of the scan. Radioactivity concentrations in discrete plasma samples were divided by the corresponding whole-blood samples to form plasma-over-blood (POB) data, and a constant POB model was fitted (i.e. average of POB values). This POB value was then multiplied by the whole blood curve to generate a total plasma curve. Parent fraction data were fitted to a sigmoid model f = ((1−(t³/(t³+ 10^a)))^b+ c)/(1 + c), where t is time and a, b and c are fitted parameters. The resulting fitted parent fraction profile was multiplied by the total plasma curve and then smoothed post-peak using a tri-exponential fit to derive the required parent plasma input function. For each scan, a time delay was fitted and applied to the input function to account for any temporal delay between blood sample measurement and the tomographic measurements of the tissue data.

Image analysis

T1-weighted magnetic resonance images (MRI) were acquired to aid in the definition of the anatomic regions of interest (ROIs) using a Siemens Magnetom Trio 3T MRI scanner at Imanova. Dynamic PET images were registered to each individual subject’s MRI, and corrected for motion using a frame-to-frame registration process with a mutual information cost function. Regional analysis was facilitated by using a human brain atlas containing ROIs that had been manually delineated on a T1-MR image template according to strict anatomical criteria. ⁷ The target ROIs included the caudate, putamen, thalamus, hippocampus, frontal cortex, parietal lobe, temporal lobe, occipital pole and cerebellum. For each subject, the template and the corresponding atlas were individually warped to the subject’s anatomical MRI which was previously linearly registered to the PET images and used to generate time activity curves (TACs).

Compartmental modelling and specific signal estimation

The individual metabolite-corrected arterial input function and the whole-blood activity curve were used in conjunction with the dynamic PET data to estimate regional volumes of distribution (V_T) using a two-tissue compartment model. The blood volume component (V_B) was fixed at 5%.

The V_T is the sum of the tracer’s specific (V_s) and non-displaceable (V_ND) binding. As there is no suitable (pseudo)reference region in the human brain, devoid of PDE4, it is not straightforward to estimate V_ND and calculate the non-displaceable binding potential (BP_ND). Target occupancy and V_ND may be estimated from the V_T data and the related plasma concentration of GSK356278 if a sufficiently large reduction in V_T is observed post-block. However, we were not able to acquire data following high levels of PDE4 blockade due to the safety limitations on the maximal dose of GSK356278 we could administer. Therefore, an estimate of the average brain binding potential ( ${\mathit{BP}}_{\text{ND}}^{\text{baseline}}$ ) was used from the literature, ² acquired in an experiment where the binding of [¹¹C](R)-rolipram and the inactive enantiomer [¹¹C](S)-rolipram in the human brain was examined. The specific binding of [¹¹C](R)-rolipram was estimated by making explicit assumptions that the (S) enantiomer demonstrates only non-displaceable binding, and the magnitude of this non-displaceable component is similar for both stereoisomers. For each subject in the current study, an individual estimate of V_ND was made using ${\text{BP}}_{\text{ND}}^{*}$  = 0.5 and the formula, ${\text{V}}_{\text{ND}}=\frac{\overline{{\text{V}}_{\text{T}}^{\text{baseline}}}}{1+{\text{BP}}_{\text{ND}}^{*}}$ , where $\overline{{\text{V}}_{\text{T}}^{\text{baseline}}}$ is the global average brain V_T before drug administration. In order to evaluate the sensitivity of our data to the assumption of the literature ${\text{BP}}_{\text{ND}}^{*}$ estimate, we repeated the process with ${\text{BP}}_{\text{ND}}^{*}$ of 0.25 and 1. For each ROI, average V_T ( $\overline{{\text{V}}_{\text{T}}})$ values per subject and scan were calculated as $\overline{{\text{V}}_{\text{T}}}=\frac{\sum _i{\alpha }_i\times V_{T_i}}{\sum _i{\alpha }_i}$ , where ${\alpha }_i$ is the volume (in mm³) of region i and $V_{T_i}$ is the total volume of distribution of region i. Finally, GSK356278 occupancy of PDE4 was calculated as $\text{occupancy}=1-\frac{{\text{BP}}_{\text{P}}^{\text{drug}}}{{\text{BP}}_{\text{P}}^{\text{baseline}}}$ , where ${\text{BP}}_{\text{P}}=\frac{K_1{k}_3}{k_2{k}_4}$ . Plasma GSK356278 concentrations ( $C_P^{\text{GSK}356278}$ ) were obtained from blood samples using mass-spectrometry with a lower limit of quantification of 0.3 ng/ml. Changes in [¹¹C](R)-rolipram specific binding following the administration of GSK356278 were related to plasma GSK356278 concentrations at the start of PET scanning using an E_max occupancy model. The model $\mathit{Occ}=\frac{C_P^{\text{GSK}356278}}{C_P^{\text{GSK}356278}+{\text{EC}}_{50}}$ was fitted to the data to obtain estimates of the half maximal effective concentration (EC₅₀). A paired sample one-tailed t-test with a significance level of 0.05 was used, under the assumption that population is normally distributed, to assess V_T changes at 3 h and 8 h compared to baseline. Considering the small sample size, the Wilcoxon signed rank test was also used to compare PET 1 with PET 2 and PET 3 with the critical value Wilcoxon W ≤ 3 for p ≤ 0.05. Six subjects completed both post-dose scans and two subjects completed only a single post-dose scan (one at 3 h and one at 8 h), therefore yielding N = 7 for both post-dose t-tests.

Results

The mean (±SD) injected dose for the 22 PET scans in this study was 270 ± 61 MBq with specific activity of 76 ± 56 GBq/µmol. After injection, the tracer readily entered the brain and showed widespread distribution. The images and the regional TACs for a representative subject are presented in Figure 1(a) and (b), respectively.

Figure 1.

(a) Anatomical MRI and PET-integrated activity from 30 to 90 min post [¹¹C](R)-rolipram injection at baseline, post-dose 1 (3 h) and post-dose 2 (8 h) for subject 2. Data for each scan have been normalized by the injected activity per liter (%ID/l). (b) Regional measured data (circles) and model fits (lines) for subject 2. For each region, the baseline (red), post-dose 1 (3 h, green) and post-dose 2 (8 h, blue) time activity curves are shown.

A two-tissue compartment model using the metabolite-corrected plasma input function described the kinetics of [¹¹C](R)-rolipram well in all ROIs. The plasmafree fractions (f_P) were similar in PET 1, 2 and 3 with mean (±SD) values of 7.4 ± 3.4%, 7.5 ± 2.6% and 8.5 ± 2.9%, respectively. [¹¹C](R)-rolipram V_T estimates are presented in Table 1 and Figure 2. The baseline [¹¹C](R)-rolipram V_T values seen in our study were consistent with those seen in previous human studies. ⁸ Intersubject variability (coefficient of variation) in the baseline V_T values (N = 8) was 20%. Oral administration of 14 mg of GSK356278 led to a mean (±SD) reduction in V_T of 17.2 ± 14.1% at approximately 3 h post-dose and 4.1 ± 19.1% at approximately 8 h post-dose. The reduction in V_T was statistically significant at PET 2 (p = 0.012, Wilcoxon W = 1) but not at PET 3 (p = 0.271, W = 6). The magnitude of the decrease was similar across the nine ROIs in the brain (range 16–20% for PET 2 and 0–8% for PET 3).

Table 1.

[¹¹C](R)-rolipram volumes of distribution (V_T), plasma GSK356278 concentration and PDE4 occupancy at post-dose scans.

	Whole-brain							ROI VT
Subject	PET time post-dose (h)	Global VT	VT reduction post-dose	Plasma GSK356278 (ng/ml)	VND estimates and PDE4 occupancy			Caudate	Putamen	Thalamus	Hippo- campus	Frontal cortex	Parietal lobe	Occipital pole	Temporal lobe	Cerebellum
					${\text{BP}}_{\text{ND}}^{*}$  = 0.25	${\text{BP}}_{\text{ND}}^{*}$ = 0.5	${\text{BP}}_{\text{ND}}^{*}$ = 1
1	Pre	0.54			0.43	0.37	0.27	0.53	0.6	0.58	0.51	0.55	0.52	0.56	0.53	0.48
	2.9	0.44	18.0%	48.5	88%	54%	36%	0.43	0.50	0.46	0.43	0.44	0.42	0.47	0.45	0.39
	7.7	0.51	5.4%	12.2	25%	16%	11%	0.49	0.59	0.52	0.5	0.52	0.52	0.53	0.5	0.43
2	Pre	0.80			0.64	0.56	0.40	0.76	0.91	0.85	0.74	0.89	0.79	0.82	0.8	0.68
	2.9	0.65	18.7%	47.2	95%	59%	39%	0.62	0.73	0.68	0.64	0.7	0.62	0.69	0.66	0.55
	8	0.76	5.0%	12.0	26%	16%	10%	0.7	0.85	0.8	0.74	0.82	0.77	0.8	0.78	0.62
3	Pre	0.74			0.59	0.50	0.37	0.68	0.85	0.79	0.73	0.75	0.7	0.76	0.72	0.64
	2.6	0.77	−4.8%	43.8	0	0	0	0.76	0.9	0.79	0.72	0.8	0.76	0.81	0.76	0.64
	7.6	0.98	−33.5%	11.5	0	0	0	0.93	1.08	1.04	0.98	1.01	0.97	1.02	0.96	0.85
4	Pre	0.68			0.54	0.46	0.34	0.65	0.81	0.7	0.68	0.72	0.66	0.69	0.63	0.57
	2.6	0.60	12.6%	30.1	57%	37%	25%	0.56	0.69	0.63	0.59	0.62	0.56	0.63	0.57	0.49
	7.5	0.62	8.5%	18.6	44%	27%	18%	0.66	0.69	0.65	0.59	0.66	0.6	0.62	0.61	0.51
5	Pre	0.83			0.67	0.57	0.42	0.82	0.96	0.91	0.86	0.85	0.8	0.86	0.82	0.63
	6.6	0.77	7.2%	16.0	30%	20%	14%	0.8	0.88	0.82	0.81	0.8	0.73	0.78	0.78	0.57
6	Pre	0.74			0.60	0.50	0.37	0.75	0.91	0.8	0.76	0.73	0.68	0.78	0.71	0.58
	2.5	0.61	17.8%	44.4	73%	50%	35%	0.58	0.76	0.65	0.63	0.63	0.56	0.65	0.59	0.46
7	Pre	0.85			0.68	0.58	0.42	0.73	1.02	0.87	0.86	0.86	0.83	0.89	0.85	0.7
	3.2	0.48	43.2%	25.2	195%	130%	88%	0.42	0.55	0.5	0.48	0.51	0.47	0.51	0.5	0.38
	7.7	0.58	31.4%	7.8	141%	94%	63%	0.5	0.69	0.61	0.59	0.59	0.56	0.61	0.59	0.48
8	Pre	1.07			0.86	0.74	0.54	0.93	1.3	1.1	1.08	1.13	1.1	1.12	1.04	0.85
	3.8	0.91	15.2%	31.4	69%	46%	31%	0.8	1.07	0.95	0.88	0.95	0.93	0.96	0.93	0.71
	7.8	1.02	4.7%	14.7	28%	17%	11%	1.02	1.16	1.07	1.02	1.06	1.03	1.06	1.01	0.77
N = 8	Pre	0.78						% VT reduction per ROI
N = 7	2.9	0.64	17.2%	38.7	54%			17	19	18	18	17	18	16	16	20
N = 7	7.6	0.75	4.1%	13.3	27%			0	8	5	4	5	4	5	3	7

V_T: volume of distribution; V_ND: non-displaceable volume of distribution; ${\text{BP}}_{\text{ND}}^{*}$ : estimated non-displaceable binding potential.

[¹¹C](R)-rolipram binding in all eight subjects expressed as V_T for the whole brain (global) and nine regions of interest (ROI). Plasma GSK356278 concentration was measured at the start of PET. The bottom three rows show the means for all subjects, i.e. N = 8 for baseline scans and N = 7 for post-dose scans.

Figure 2.

Global brain volumes of distribution (V_T) for baseline, 3 and 8 h post-dose scans in all eight subjects. Subjects 5 and 6 completed one post-dose scan.

The individual plasma GSK356278 concentration–time curves are presented in Figure 3. The mean C_max was 42.3 ng/ml. Figure 4 shows plasma GSK356278–PDE4 occupancy plots using three different estimates of ${\text{BP}}_{\text{ND}}^{*}$ . Visual inspection suggested that subjects 3 and 7 may be outliers, with ΔV_T (and occupancies) for these subjects being notably different from the mean ΔV_T of the sample as a whole. Plasmafree fractions were similar between the three PET scans in subject 3 (range 6.3–9.3%) and subject 7 (range 11.8–12.0%)_. The modified Thompson Tau test was used to determine if the ΔV_T for these subjects were outliers. The method takes into account the sample mean, SD and N, and provides a statistically determined rejection value. ΔV_T of subject 7 for both post-dose scans exceeded the rejection region of 1.7 SD. Subsequently, the test was applied to ΔV_T with N = 6 and identified both post-dose scans of subject 3 as outliers. The test on N = 5 did not identify additional outliers. We therefore examined the group ΔV_T again without subjects 3 and 7. The reduction (±SD) in average V_T for N = 5 was 16.4 ± 2.5% and 6.1 ± 1.7% for 3 h and 8 h data respectively (p < 0.001, W = 0). The estimated mean PDE4 occupancy was 49 ± 8% at 3 h, and 19 ± 5% at 8 h. The in vivo affinity of GSK356278, excluding subjects 3 and 7, was estimated as EC₅₀ = 46 ± 3.6 ng/ml, leading to an estimated PDE4 occupancy of 48% at plasma T_max.

Figure 3.

Individual plasma GSK356278 concentration-time curves. The dotted lines show the average time of the post-dose scans.

Figure 4.

Model fits of the PET occupancy data and plasma GSK356278 concentration with estimated EC₅₀ in all subjects examined (top row) and for the dataset excluding subjects 3 and 7 (bottom row). Column 2 represents the most likely value of the non-displaceable binding potential (BP_ND) while the other two columns are intended to show the effect of a different estimate of BP_ND on the EC₅₀.

Discussion

To our knowledge, this is the first study to explore the blockade of the human PDE4 enzyme in the brain in vivo. Baseline V_T values obtained in this study were consistent with those seen in previous human studies.^2,8 The between-subject variability (coefficient of variation) of the baseline V_T was 20%, comparable to earlier reports of 25%. ⁸ Oral administration of the PDE4 inhibitor GSK356278 led to a mean V_T reduction of 17% around 3 h post-dose compared to the baseline V_T and a reduction of 4% around 8 h post-dose, consistent with the hypothesis that GSK356278 enters the brain readily and binds to PDE4. V_T change in all subjects followed the plasma concentration of GSK356278, with ΔV_T PET 2 > than ΔV_T PET 3. The estimated relationship between the plasma concentration of GSK356278 and PDE4 occupancy indicates that occupancies of close to 50% of PDE4 are achieved at T_max following the administration of single oral doses of 14 mg of GSK356278. The available data provide no evidence for indirect pharmacokinetics for GSK356278 in the human brain, suggesting that the plasma EC₅₀ (46 ng/ml) estimated in this study can be used to calculate PDE4 occupancy following repeat dose administration.

The assessment of the relationship between PDE4 occupancy and plasma GSK356278 concentration was complicated by our inability to estimate the [¹¹C](R)-rolipram BP_ND directly from the study data, due to the relatively low levels of occupancy achieved. Our population estimate of BP_ND from the literature ignores between-subject variability in PDE4 expression. Although an error in the estimate of BP_ND would lead to an error in the estimated EC₅₀, a relatively large range of BP_ND estimates (0.25–1) produced a modest difference in GSK356278 EC₅₀ (21–84 ng/ml), indicating that our estimated EC₅₀ is relatively robust to variability in assumed ${\text{BP}}_{\text{ND}}^{*}$ . In the absence of higher levels of occupancy or associated estimates with the inactive enantiomer [¹¹C](S)-rolipram, it is not possible to confirm individual variability in ${\text{BP}}_{\text{ND}}^{*}$ directly. GSK356278 EC₅₀ of 46 ng/ml with a plasma C_max of 42.3 ng/ml leads to an estimated PDE4 occupancy of 48% at T_max, with a range of 34–67% (depending on the BP_ND estimate).

Subjects 3 and 7 were identified as outliers based on abnormal post-dose V_T compared to the baseline. Specifically, subject 3 showed an increase in post-dose V_T, whereas subject 7 showed an exceptionally large reduction, leading to occupancy at 3 h > 100%. The metabolite-corrected arterial input function and plasmafree fraction of [¹¹C](R)-rolipram for these subjects did not differ between scans, making it unlikely that changes in blood flow or plasma protein binding could explain these findings. Variability in estimated parameters is thus the most likely explanation for the findings in these subjects. Test–retest variability of V_T in healthy humans was shown to be approximately 19% in an earlier study of 12 subjects. ⁸

PDE enzyme activity dysfunction has been implicated in disease states such as asthma, ischemic stroke and CNS disorders. In a study of patients with major depressive disorder, [¹¹C](R)-rolipram binding was reduced by 18% compared to healthy control subjects and could be partially normalized with selective serotonin reuptake inhibitor treatment. ⁹ PDE4 selectively metabolizes cAMP in the brain to the inactive monophosphate and is therefore an important component of the cAMP cascade. The enzymatic activity of PDE4 is regulated by protein kinase A (PKA) via a feedback mechanism, with high concentrations of cAMP stimulating PKA to phosphorylate PDE4, thereby increasing its enzymatic activity and returning the concentration of cAMP to steady state. Unilateral injection of a PKA activator and inhibitor into the rat striatum was shown to significantly increase and decrease, respectively, the binding of [¹¹C](R)-rolipram as measured with PET. ³ Upregulation of the cAMP cascade through long-term pharmacological inhibition of the PDE4 enzyme is a promising therapeutic intervention for a range of conditions. In preclinical studies, GSK356278 was shown to improve performance in an object retrieval test in cynomolgus macaques, ⁶ consistent with the reported effects of rolipram in various tests of cognition. ¹⁰ Despite the possible benefits of brain-penetrant PDE4 inhibitors, clinical use has been limited by mechanism-dependent adverse events such as nausea and emesis. ¹¹ The highest dose in this study was limited to 14 mg (equivalent to ∼50% occupancy) by these adverse events in Phase 1 studies.

In conclusion, we present the first human report of PDE4 occupancy measured directly in the human brain with PET. Our data will be used in conjunction with the known plasma GSK356278 pharmacokinetics to determine optimal doses to be used in future clinical development.