Study in non-human primates
The study was approved by the Animal Research Ethical Committee of the Northern Stockholm Region. Three adult female cynomolgus monkeys (mean weight 6.1 kg) were included. The non-human primates (NHPs) were owned by the Centre for Psychiatry Research, Department of Clinical Neuroscience, Karolinska Institutet, and housed in the Astrid Fagraeus Laboratory, Karolinska Institutet, Solna, Sweden.
Study design
The study in NHP comprised of 8 experimental days. Two PET-measurements with [11C]VC-002 were performed on each day. A baseline PET-measurement was followed by a second measurement 2.5h later either at baseline conditions (test-retest) or after i.v. administration of tiotropium (pretreatment). Two of the monkeys (NHP1 and NHP2) participated in the test-retest measurements. All three monkeys (NHP1-3) participated in two pretreatment sessions each. In the pretreatment measurements (n=6), tiotropium was infused over 15 min starting 20 min before the [11C]VC-002 injection (for details see supplement Figure, Supplemental Digital Content 1). The six doses of tiotropium varied from 0.03 to 1.0 μg/kg. The highest doses were expected to induce saturating levels of tiotropium binding in lungs. Such conditions were feasible as higher doses were well tolerated in NHP.
Measurement of tiotropium plasma concentration
At each PET measurement with active drug, venous blood samples were collected for determination of tiotropium plasma concentration. The samples were drawn before tiotropium administration at approximately 35 min before [11C]VC-002 injection, and after tiotropium administration at 0.5, 1, 5, 15, 30 and 60 min after [11C]VC-002 injection. The plasma concentration of tiotropium was measured by solid phase extraction, followed by high-performance liquid chromatography and atmospheric pressure chemical ionization tandem mass spectrometry17. The lower limit of quantification (LLOQ) was 5 pM. The area under the plasma concentration curve for the time of the PET measurement (90 min) was calculated and divided by the duration to obtain the average plasma concentration during PET.
PET measurements
General anaesthesia was induced by intramuscular injection of ketamine hydrochloride (approximately 10 mg/kg) and after endotracheal intubation maintained by administration of a mixture of sevoflurane (2–8%), oxygen, and medical air. Head immobilization and safety monitoring procedures (body temperature, heart rate, blood pressure, fluid balance) were performed as described previously15.
The radioligand [11C]VC-002 was prepared at the PET center at Karolinska University Hospital as previously described13,15. In each PET measurement, a sterile physiological, phosphate-buffered (pH 7.4) saline (PBS) solution containing [11C]VC-002, in a volume not exceeding 5mL, was injected as a bolus into a sural vein during 5 s. PET data acquisition started at time of the bolus injection. In all measurements the radiochemical purity of [11C]VC-002 exceeded 99% at time of injection. The mean radioactivity injected was 154 MBq (SD ±9 MBq, range 139 – 172 MBq, N=16 measurements). The molar activity at time of injection was 689 GBq/µmol (SD ±814 GBq/µmol, range 129 – 3195 GBq/µmol) corresponding to an injected mass of 0.25 μg (SD ±18 μg, range 0.025 – 0.57 μg).
PET measurements were conducted using the high-resolution research tomograph (HRRT) (Siemens Molecular Imaging). Radioactivity in lungs was measured continuously over 63 min using an imaging protocol described in detail previously15. In short, a transmission scan of 6 min using a single 137Cs source was performed immediately before [11C]VC-002 injection. Data were acquired continuously in list mode for 63 min after i.v. injection of [11C]VC-002. Images were reconstructed for a series of time frames (9 × 10 s, 2 × 15 s, 3 × 20 s, 4 × 30 s, 4 × 1 min, 4 × 3 min, and 7 × 6 min).
Blood sampling for the measurement of radioactivity
Venous blood samples (1–3 mL) were obtained manually at 1, 2, 3, 5, 15, 30, 45 and 60 min for the measurement of radioactivity in whole blood and plasma using a well-counter. The procedure has previously been described in detail15.
Region of interest definition
Regions of interest (ROIs) were manually delineated for the lungs (pooled) and the arch of the thoracic aorta on the summation (0-60 min integral) and the early (0-2 min) PET images, respectively. The ROIs were applied to the series of PET images to generate time-activity curves (TACs).
Study in humans
Subjects and study design
The study was approved by the Regional Ethics Committee in Stockholm and the Radiation Protection Committee at the Karolinska University Hospital, Stockholm. Written informed consent was obtained from each subject.
Seven male subjects, 20-50 years of age, were recruited and studied at the PET Centre, Department of Clinical Neuroscience, Karolinska Institutet, Stockholm, Sweden. Subjects were healthy according to medical history, clinical examination, and routine laboratory blood and urine tests. No medications were allowed at time of the study.
The PET-measurements were carried out on a whole-body GE Discovery 710 PET/CT-system at the Department of Nuclear Medicine, Karolinska University Hospital, Solna. Each subject participated in two or three PET examinations that were performed on separate days.
In part 1, subjects H1-H3 first participated in a baseline (BL) measurement. A second PET measurement was performed 7-12 days later and 2 hours after inhalation of a therapeutic dose of tiotropium (18 µg).
In part 2, each of four subjects H4-H7 participated in three PET measurements. The first PET (BL) was followed by two PET measurements, each one after pretreatment with tiotropium (18 µg). The pretreatment PET measurements started 30 min after inhalation of tiotropium in subjects H4-H5 and 2 hours after inhalation in subjects H6-H7 (for details see supplement Figure, Supplemental Digital Content 1). In both study parts, adverse events were monitored on experimental days and up to one week after the last measurement via a follow-up telephone call.
Measurement of tiotropium plasma concentration
In part 1 (subjects H1-H3) the concentration of tiotropium in plasma was not measured. In part 2 (subjects H4-H7), 3 venous blood samples were obtained at start, middle time and end of the PET measurement, i.e. 30–90 or 120–180 min post inhalation of tiotropium. A standard LC-MS/MS analysis of tiotropium concentration in venous plasma was performed 17. The lower limit of quantification (LLOQ) was 1.3 pM. The average concentration during PET acquisition was obtained by calculating the area under the plasma concentration curve and dividing it by the duration of PET.
PET/CT measurements
PET examinations, radiochemistry, blood sampling, image processing, and quantification were performed essentially as described in detail in our previous publication on the test-retest reliability of [11C]VC-002 binding in the human lung16. The radiochemical purity of [11C]VC-002 exceeded 99% at time of injection in all measurements. The mean injected radioactivity was 220 MBq (SD ±37 MBq, range 165-296 MBq, N=18) and the molar activity was 398 GBq/µmol (SD ±250 GBq/µmol, range 51-955 GBq/µmol). The corresponding mass of the radioligand injected was 0.34 μg (SD ±0.34, range 0.08-1.24 μg).
At time of imaging examination each subject was positioned in the PET/CT system head-first, supine with the chest located within the 16 cm field of view. Initially, a low-dose CT was obtained for the chest. After i.v. injection of [11C]VC-002, radioactivity was acquired in list mode for 63 min providing a reconstructed 4D PET image with 33 timeframes with equivalent frame timings as for the NHP study.
The PET image was evaluated and, if possible, corrected for inter-frame subject movements as described previously16. In short, the reconstructed 4D PET image was frame-by-frame realigned to have a consistent subject position throughout. In case of mismatch of the PET frame and the CT scan used for attenuation correction, the PET image was re-reconstructed to minimize this effect of a movement artifact, by ensuring the correct alignment of the attenuation map for each time frame. The final corrected image was quality checked for the presence of potential residual artifacts related to severe intraframe subject movement.
Blood sampling for measuring [11C]VC-002 radioactivity concentration
Before PET an arterial cannula was inserted into one of the radial arteries. After injection of [11C]VC-002, 3 arterial blood samples (2 ml each) were drawn manually at 10, 25, and 45 min for the purpose of determining the average plasma-to-whole-blood ratio.
Regions of interest delineation
ROIs for the lungs and the arch of the aorta were delineated automatically as described previously16. In short, the method relied on the co-registered CT images to identify voxels in the body with a density below water based on the CT Hounsfield unit value. The binary image of such voxels was then refined using the summation PET image and image processing morphological operations to arrive at identifying a ROI covering both lungs. Voxels falling within the volume of the arch of the aorta were identified based on the voxels’ early time course of high radioactivity after iv. injection. In addition, the search was restricted to the mediastinal space, and distinguished between venous (e.g. vena cava) and arterial voxels based on the relative timing of peak radioactivity. Finally, the ROIs were applied to the series of PET images to obtain TACs for the whole lung (bilateral) and aorta, respectively.
Quantification of radioligand binding and receptor occupancy
Quantification of [11C]VC-002 binding in lungs
The total plasma radioactivity concentration was used as the input function for quantification of binding as previously described15,16. To derive this curve, the TAC for the aorta was multiplied with the average plasma-to-whole-blood radioactivity ratio calculated from the drawn blood samples as described above.
The parameter used to express [11C]VC-002 binding was the total volume of distribution (VT), which is an index of total binding in tissue, i.e. the sum of non-specific and receptor-specific binding. VT is numerically equivalent to the tissue partition coefficient, i.e. the ratio of tissue to plasma concentration in case of proper steady state conditions.
The quantification was carried out for each volume element (voxel) using data-driven estimation of parametric images based on compartmental theory (DEPICT), as described previously16,18. The range of exponents, employed in the calculation of the table of kinetic basis functions used by DEPICT, was between 0.0136 1/min and 0.6 1/min (as per 18). The main output of DEPICT was voxel-wise parametric images of VT in the lungs.
Quantification of receptor occupancy
Receptor occupancy calculation is straightforward when the availability of a reference region with negligible binding provides an estimate of non-specific binding and thus allows for differentiation between specific and non-specific binding in the target region. However, in pulmonary imaging there is no reference region and the specific binding component of a radioligand such as [11C]VC-002 cannot be directly quantified16. The estimation of receptor occupancy following drug administration was instead estimated using the Lassen plot, an established indirect approach19,20.
The indirect Lassen approach is based on the total binding (VT) values and relies on the assumptions that various parts of the lungs have different levels of mAChR expression, and thus specific [11C]VC-002 binding; and that the background of non-specific binding as well as drug exposure at the target is similar across the organ. Based on these assumptions, lung VT data from baseline and pretreatment conditions were entered into a graphical evaluation (Lassen plot), which provided estimates of receptor occupancy as well as the level of non-specific (non-displaceable) binding, VND19–21.
The scatter points for the Lassen plot analysis were obtained from the parametric images of VT in lung tissue, following a downsampling procedure to reduce voxel-wise variation of VT estimates, and, consequently, the spread (variance) of the scatter points in the plot. In detail, the voxel count in the parametric images was reduced by a factor of 4 along each dimension, i.e. a 64-fold reduction in the total number of voxels, by calculating the mean VT estimate for each adjacent 4x4x4 cube of voxels. Then the voxel-wise VT estimates within the lungs at baseline and at pretreatment were extracted and projected into the Lassen plot. For this purpose, the lung ROI mask, obtained at the original PET image resolution, was first eroded to discard voxels at the edge of the lungs (within a 2-voxel margin from the surface), and then downsampled in two steps. Each downsampling step reduced the volume size by a factor of 2 along each dimension, i.e. a 8-fold reduction in the total count of voxels, and then retained voxels as “within-lung” which had a downsampled value >0.5 (i.e. >50% probability for lung). The final downsampled lung ROI mask had the same resolution as the downsampled VT parametric image and could thus be applied for the extraction of voxel-wise lung VT values.
Lassen’s original graphical approach entailed plotting the difference in VT between baseline and pretreatment conditions vs. the baseline VT values across lungs in a scatter plot. Receptor occupancy was estimated by performing a constrained, weighted linear fit that gave more weight to scatter points in the Lassen plot with a higher local kernel density. A modified variant of the Lassen approach was used in one subject (H6), which, in lack of usable baseline data, relied on VT data from only one pair of experiments following drug inhalation as well as the measured average drug concentrations in plasma during PET21.
The presently obtained receptor occupancy values were compared to the previously reported test-retest study in humans16. For that purpose, the test-retest dataset was re-analyzed using the same settings as in the present study, i.e. using only the first 63 min of the acquired data and the same DEPICT settings as specified above.
Assessment of the plasma exposure–occupancy relationship
The occupancy values obtained from the Lassen plot were entered into an analysis of the curvilinear relationship between drug plasma concentration and pulmonary receptor occupancy. In detail, a weighted non-linear least square curve fit was performed using the following equation:
where Occ is occupancy, Occmax is the maximum level of occupancy achievable by tiotropium at the target binding sites detected by [11C]VC-002, Cp is the average plasma tiotropium concentration during PET and Ki is the apparent inhibition constant, i.e. the plasma concentration when Occ equals half of Occmax.
The model was fitted to the occupancy values estimated obtained from the Lassen plot analysis, yielding estimates of Ki. The maximal occupancy, Occmax, was either fixed at 100% or, alternatively, estimated by the fit. The R-squared values of the linear fit in the Lassen plot occupancy estimation were used as weighting factors when fitting the occupancy model, i.e. in effect giving more weight to occupancy values that were more certain. The model was fitted in separate analyses to NHP (i.v. administration-based) and human (inhalation based) data, respectively.
Statistics
Processing and computations were performed using the Matlab, version R2014b (www.mathworks.com). The estimation of occupancy using weighted linear least squares provided statistical assessment of the goodness-of-fit, such as the coefficient of determination, R2, the standard error and a 95% confidence interval for the coefficients. A simple t-test was separately performed on the inter-individual mean occupancy values for each group (pretreatment or test-retest) and species (NHP or human). Statistical assessment was performed on the results of the exposure–occupancy model fits: 1) to ascertain the goodness of fit (among others calculating the standard error and 95% confidence interval of fitted parameters), 2) to compare alternative models of the NHP data (by using an F-test, Akaike’s information criterion values, AIC, and adjusted R-squared values), and 3) to compare i.v. administration based (NHP) and inhalation based (human) Ki estimates (by using two-tailed t-tests). In all analyses, the statistical significance (alpha level) was set at p < 0.05.