Abstract
The histamine H3 receptor is implicated in the pathophysiology of several central nervous system disorders. N-methyl-6-(3-cyclobutyl-2,3,4,5-tetrahydro-1H-benzo[d]azepin-7-yloxy)-nicotamide (GSK189254) is a highly potent, selective, and brain-penetrant H3 receptor antagonist. Previous studies in the pig using PET have shown that 11C-GSK189254 uptake in H3-rich regions of the brain can be blocked by the selective H3 antagonist ciproxifan. The purpose of the present study was to evaluate 11C-GSK189254 as a PET radioligand for human studies and to determine the dose–receptor occupancy relationship of GSK189254 in the human brain. Methods: Dynamic PET scans were obtained in healthy subjects over 90 min after intravenous administration of approximately 370 MBq of 11C-GSK189254. Blood samples were taken throughout the scans to derive the arterial plasma parent input function. Each subject was scanned twice, either with tracer alone (test–retest) or before and after a single oral dose of GSK189254 (10–100 μg). Data were analyzed by compartmental analysis, and regional receptor-occupancy estimates were obtained by graphical analysis of changes in the total volumes of distribution (VT) of the radioligand. Results: 11C-GSK189254 readily entered the brain; its regional brain distribution reflected the known distribution of H3 receptors, with high binding in the caudate and putamen, intermediate binding in cortical regions, and low binding in the cerebellum. GSK189254 displayed a high receptor affinity, and a marked reduction in VT was apparent at all the doses tested. The oral dose equaling 50% occupancy of the available receptor sites (ED50) was estimated as 4.33 μg. Additional data on plasma pharmacokinetics after oral dosing and the plasma free fraction gave a corresponding estimate of the free concentration of GSK189254 required to occupy 50% of the available receptor sites (EC50) (0.011 nM). The test–retest data showed reductions in regional VT on the second scan in all subjects. A nonlinear compartmental analysis of this effect demonstrated that this reduction was consistent with carryover of a tracer mass dose effect with an estimated in vivo apparent dissociation constant of 0.010 nM, close to the independent estimate of the plasma EC50. Conclusion: 11C-GSK189254 can be used to quantify H3 receptor availability in humans in vivo using PET but requires high specific activity; the possibility of tracer mass dose effects should be carefully analyzed.
Central histaminergic neurones originate almost exclusively in the tuberomammillary nucleus of the posterior hypothalamus and project to almost all regions of the mammalian brain (1). Histamine has a global coordinating function in the brain (2), particularly in relation to the regulation of the sleep–wake cycle (3), and mediates its biologic effects through activation of 4 transmembrane G-protein–coupled receptors (H1, H2, H3, and H4) (4,5). H3 receptors are heterogeneously expressed in the mammalian brain, with particularly high densities in the basal ganglia (6). They are located presynaptically on histaminergic neurones, where they act as autoreceptors regulating the synthesis and release of histamine (7). In addition, H3 receptors act as heteroreceptors on nonhistaminergic neurones, controlling the release of other neurotransmitters, including acetylcholine and dopamine (8,9).
The H3 receptor has been implicated in a range of disorders including cognitive deficits, obesity, diabetes mellitus, and sleep disorders (3,10,11). Considerable interest in the pharmacology of this receptor exists, and several pharmaceutical companies are actively researching and developing H3 agonists and antagonists as drug candidates for the treatment of cognitive disorders (12).
Since the discovery of the H3 receptor more than 25 y ago, several ligands, radiolabeled with either 3H or 125I, have been developed as probes for investigation of its physiology, pharmacology, and distribution within the brain, including 3H-(R)α-methylhistamine and 125I-iodoproxyfan (13–15). Until recently, there has been a lack of suitable radioligands for imaging of the H3 receptor in vivo in the central nervous system of humans using PET. A few compounds previously investigated as potential radioligands for PET include the iodoproxyfan analog 18F-fluoroproxyfan (16) and the thioperamide analog 18F-VUF5000 (17). Although 18F-fluoroproxyfan may still prove to be potentially useful to the imaging of H3 receptors in the human brain, other compounds, including 18F-VUF5000, do not show the appropriate characteristics necessary for successful imaging agents. 18F-VUF5000, for instance, shows a lack of brain penetration.
Hamill et al. (18) have described 2 new H3 inverse agonist PET radioligands that look promising for imaging and quantifying H3 receptors. In addition, Plisson et al. (19) have described the radiosynthesis and preclinical evaluation of 11C-N-methyl-6-(3-cyclobutyl-2,3,4,5-tetrahydro-1H-benzo[d]azepin-7-yloxy)-nicotamide (11C-GSK189254) as a novel PET radioligand for the H3 receptor. GSK189254 is a novel, highly potent, selective, and brain-penetrant H3 receptor antagonist, with high affinity for both human and rat H3 receptors (−log(10) Ki, where Ki is the inhibition constant, 9.59–9.90 and 8.51–9.17, respectively) (20). Previous studies in the pig have shown that 11C-GSK189254 is brain-penetrant, has suitable kinetics for quantification, and demonstrates uptake (which can be blocked by the H3 antagonist ciproxifan) in H3-rich regions of the brain.
Here we describe the evaluation of 11C-GSK189254 in the human brain using PET. Healthy volunteers underwent 2 scans, either under test–retest conditions or before and after an oral dose of GSK189254 (10–100 μg), to assess specific binding of 11C-GSK189254 and determine the dose–H3 receptor occupancy relationship of GSK189254 in the human brain.
MATERIALS AND METHODS
The study was performed at the Centre for Addiction and Mental Health, Toronto, on behalf of GlaxoSmithKline.
Radiochemistry
6-[(3-cyclobutyl-2,3,4,5-tetrahydro-1H-3-benzo[d]azepin-7-yl)oxy]-3-nicotamide (GSK185071B) and GSK189254 were synthesized at GlaxoSmithKline. 11C-GSK189254 was prepared by 11C-methylation of the carboxamide precursor GSK185071B with 11C-iodomethane (21) inside a high-performance liquid chromatography (HPLC) sample loop using the previously described loop method (22). Briefly, the normethyl precursor (0.5 mg) was dissolved in a mixture of dimethylformamide (80 μL) and aqueous potassium hydroxide (2N, 2 μL) and reacted on the loop with trapped 11C-iodomethane for 4 min at room temperature. After purification by HPLC (250 × 10 mm; 20/80 acetonitrile/0.1N ammonium formate, pH 4; 7 mL/min) (Luna C18; Phenomenex) and formulation in buffered saline, 11C-GSK189254 was obtained with 15%−20% radiochemical yield (uncorrected for decay, from 11C-CO2) at 27 min after the end of bombardment.
Study Subjects
Data from 21 healthy volunteers (16 men and 5 women) with a mean age of 43 y (range, 35–67 y) are reported for this study. Each subject underwent a physical examination, including a 12-lead electrocardiogram, vital-sign assessment, and blood analyses for substances of abuse. Each subject provided written informed consent, and the study was approved by the local research ethics board and Health Canada. On each day of PET, the subjects were cannulated in the forearm or antecubital vein for tracer administration and in the contralateral radial artery for blood sampling.
PET with 11C-GSK189254
Each subject was scanned twice on a Biograph HiRez XVI PET camera system (Siemens Molecular Imaging), the performance characteristics of which are described by Brambilla et al. (23). On each scanning occasion, a scout view was acquired for the accurate positioning of the subject, followed by a low-dose (effective dose, 0.2 mSv) CT scan. A 511-keV attenuation image was generated from this CT image, which was then forward-projected to generate the attenuation correction for the emission PET sinograms. After the transmission scan, a 32-bit, list-mode 90-min emission scan was acquired. 11C-GSK189254 (247–433 MBq) was administered by intravenous bolus injection via a cannula. The specific activity at the time of injection was 18.0–79.5 GBq/μmol. Because GSK189254 has a molecular weight of 351.5, the mass dose associated with these injections was 1.63–4.82 μg. The emission list-mode data were rebinned into a series of 3-dimensional sinograms (span, 11; ring difference, 27). Frames were defined as 1 × background frame (start of scan to time of injection, ∼1 min), 8 × 15 s, 3 × 60 s, 5 × 120 s, 5 × 300 s, and 5 × 600 s. For each 3-dimensional sinogram, the data were normalized, and attenuation and scatter were corrected before Fourier rebinning was applied to convert the 3-dimensional sinograms into 2-dimensional sinograms. The 2-dimensional sinograms were then reconstructed into image space using a 2-dimensional filtered backprojection algorithm, with a ramp filter at a Nyquist cutoff frequency. After reconstruction, a gaussian filter (5 mm full width at half maximum) and calibration factor were applied.
Six of the subjects underwent 2 baseline scans (test–retest), with approximately 2 h between scan starts. Fifteen subjects received 1 baseline scan and an oral dose of GSK189254 (10–100 μg) at either 4 or 24 h before a second PET scan. Subjects also underwent 3-dimensional volumetric MRI using a fast spoiled gradient-recalled echo sequence acquired in the axial plane, to provide an anatomic image for coregistration with their PET data. MRI scans were obtained using a Signa Excite HD 1.5-T scanner system (GE Healthcare).
Blood Sampling and Metabolite Analysis
Arterial blood samples were collected continuously at a rate of 5 mL/min during the first 10 min after 11C-GSK189254 administration. Further discrete samples were collected throughout each PET scan for measurement of the whole-blood radioactivity, plasma radioactivity, and parent fraction. For radiolabeled metabolite analyses, HPLC was performed by minor modifications of the method described by Hilton et al. (24). Briefly, plasma samples (4, 7, 10, 15, 30, 45, 60, and 75 min after 11C-GSK189254 administration) were loaded onto a 5-mL HPLC injector loop (Valco) and injected onto a small capture column (4.6 × 20 mm) packed in-house with OASIS HLB (30 μm; Waters Corp.). The capture column was eluted with 1% aqueous acetonitrile (2 mL/min) for 3–4 min and then backflushed (25% acetonitrile/75% H2O + 0.1N ammonium formate, pH 4, 2 mL/min) onto a 10-μm Luna C18 column (250 × 4.6 mm; Phenomenex). Both column effluents were monitored through a flow detector (Flow-Count; Bioscan) operated in coincidence mode. All radioactivity data were corrected for physical decay and integrated using a personal computer. The recovery of radioactivity in the HPLC system was greater than 97% in all cases.
Kinetic Analysis of PET Data
A metabolite-corrected plasma input function was generated from the blood and plasma time–activity profiles and metabolite data. For each subject, regions of interest (ROIs) (cerebellum, insula cortex, anterior cingulate cortex, frontal cortex, parietal cortex, occipital cortex, caudate, putamen, hippocampus, amygdala, and thalamus) were defined on the MR image for that subject using a region template (25). The MR image was coregistered to each PET image for that subject using SPM2 (26). The ROIs were subsequently applied to the dynamic PET data to generate tissue time–activity curves. For presentation purposes, tissue radioactivity concentrations were normalized for the radioactive dose injected and for the subject's body mass, giving standardized uptake values ([kBq/mL of tissue]/[kBq/g of body mass]). Data are presented as mean ± SD.
A 2-tissue-compartment model was applied to the plasma and tissue data, with rate constants K1 and k2 describing the exchange of tracer between plasma and tissue and rate constants k3 and k4 describing the exchange of tracer between the 2 tissue compartments. High initial estimates of K1 indicated that there were significant dispersion effects in the measured blood activity, attributable to the length of the sampling cannula (∼1.25 m). Delay and dispersion terms were therefore incorporated in the model (Appendix). Because of the slow kinetics and correspondingly high estimates of the total volumes of distribution (VT) for this radioligand, it was also necessary to apply constraints to the model using a shared-parameter approach. Thus, individual estimates of K1 and k3 were obtained for each ROI, whereas estimates of the K1/k2 ratio and k4 were obtained globally across all ROIs for each scan (Appendix). Time stability was analyzed by sequentially examining datasets with decreasing scan durations, allowing for an assessment of what scanning duration was required to obtain a stable estimate of the VT in each region. Results are expressed in terms of the VT (the equilibrium partition coefficient between tissue and plasma, VT = K1/k2(1 + k3/k4)) for each ROI (27).
Because there was no suitable reference region within the brain to quantify free plus nonspecific binding (VND) (27), estimates of fractional occupancy (Occ) and VND for each subject were obtained by graphical analysis according to the occupancy plot described by Equation 1.Eq. 1where and are the VTs across all regions obtained at baseline and after drug administration, respectively (28).
Estimates of the oral dose of GSK189254 required for 50% occupancy of the available receptor sites (ED50) were obtained by fitting an Emax model (Eq. 2) to a plot of fractional occupancy versus oral dose (D) across all subjects.Eq. 2where Emax is the apparent maximal attainable occupancy.
Data were also available from a previous pharmacokinetic study of the concentration of the cold drug in plasma after oral administration, together with an estimate of the free fraction of the drug in plasma. These data were used to convert the ED50 to the equivalent free concentration of GSK189254 required to occupy 50% of the available receptor sites (EC50), assuming passive transport of GSK189254 across the blood–brain barrier.
After the estimation of the EC50, the test–retest data were analyzed further to investigate the possibility of the mass dose associated with the tracer injection affecting occupancy measurements. In brief, regional data from corresponding test and retest scans were fitted in sequence to a parallel 2-tissue-compartment model, with equations describing the time courses of radioactivity in nondisplaceable and specifically bound compartments together with the concentration time courses of cold drug in the 2 compartments (Appendix).
RESULTS
Representative images from the baseline scans are shown in Figure 1. After intravenous administration of 11C-GSK189254, radioactivity readily entered the brain, accumulating in each ROI including the cerebellum. The distribution was heterogeneous, with highest concentrations in the caudate and putamen. Representative time–activity curves from a baseline scan are shown in Figure 2, together with a corresponding plasma parent input function. The parent fraction of 11C-GSK189254 in plasma was 86% ± 5% after 75 min. The kinetics of the uptake of 11C-GSK189254 in the brain were slow in that there was no obvious washout phase in the regional time–activity curves within the time frame of the baseline scans. Attempts to fit compartmental models to high-binding individual regions alone frequently gave unrealistic estimates of VT because of large errors in the off-rate constant (k4). Fits were therefore obtained by simultaneous analysis of all regions, applying the constraints described in the “Materials and Methods” section, with global estimates of the K1/k2 ratio, k4, and dispersion rate constant (kdisp) for each scan.
Test–retest results are shown in Table 1. In each region studied, there was a significant reduction in VT in the second scan relative to the first. The significance of this observation will be commented on later.
Figures 1 and 2 show the reduction in the overall brain uptake of 11C-GSK189254 after oral administration of GSK189254. A marked reduction in VT was apparent at all the doses tested, even doses as low as 10 μg (Table 2). This effect was seen in all regions, including the cerebellum, which was therefore unsuitable for use as a reference region devoid of specific binding.
The overall estimate of kdisp was 0.025 ± 0.016 s−1. Simulations (results not shown) demonstrated a negligible propagation of error from kdisp to VT.
The effects of scan length on the stability of estimates of VT are illustrated in Figure 3. There was no apparent bias in low-, intermediate-, or high-binding regions, but as expected for this slow ligand the increase in variance was more marked in high-binding regions, indicating a need for maximal scanning times.
The estimate of the K1/k2 ratio for the scans reported in Table 2 was 4.4 ± 3.4. However, because this estimate is not a stable macroparameter (29) in the absence of a suitable reference region, estimates of receptor occupancy and VND were therefore obtained from occupancy plots as described in the “Materials and Methods” section.
Figure 4 shows occupancy plots for each subject who received GSK189254 at 4 h before the second scan. Occupancies were high across all doses including the lowest oral dose of 10 μg. High occupancies were also seen 24 h after an oral dose of 25 μg (67.6 ± 12.2, n = 5), consistent with the persistence of the drug in the systemic circulation. The overall estimate of VND was 10.54 ± 1.78.
The relationship between dosage and receptor occupancy at 4 h for GSK189254 is shown in Figure 5. A least-squares fit of Equation 2 to these data gave an estimate of ED50 = 0.055 ± 0.014 μg/kg of body weight. The mean body weight of these subjects was 79.4 kg, corresponding to an oral dose ED50 of 4.33 μg per subject. In-house data were available on the plasma concentration time course of GSK189254 in human volunteers after a range of oral doses from 25 to 500 μg per subject. These data were normalized to the ED50, assuming linearity, giving a mean plasma concentration over the scanning period (4–5.5 h) of 0.031 nM. A 0.36 value for the plasma free fraction (fP) (in-house data) gives an EC50 of 0.011 nM.
This observation of high occupancies at low oral doses suggested that a decrease in VT in retest data relative to test data could be due to a tracer mass dose effect with carryover from the test to the retest scan, even at the low mass doses (<5 μg) associated with the injection of 11C-GSK189254. The test–retest data were therefore reanalyzed using a model taking account of this carryover, formulated in terms of a regional binding capacity (Bmax) and apparent dissociation constant (KD) with occupancy of receptor sites varying over time (Appendix). Although the test–retest part of the study was clearly not originally designed as a saturation experiment, it was nonetheless possible in this way to obtain an independent estimate of the KD by simultaneously fitting multiple regions across the test and retest scans. Satisfactory fits were obtained in 5 of 6 subjects, with 1 outlier. A typical fit is shown in Figure 6, and the results are summarized in Table 3.
The estimation of KD in Table 3 (0.0095 nM) used the mean value of VND (10.54) and fP (0.36) (Appendix). This estimate of KD agrees well with that of the plasma EC50 obtained from the analysis of the occupancy studies (0.011 nM). This agreement is consistent with the possibility that the reduction in VT between test and retest (Table 2) was due to a tracer mass dose effect. Estimates of Bmax for low (cerebellum), high (putamen), and intermediate (cortex) regions are also shown in Table 3.
DISCUSSION
The present study was designed to investigate 11C-GSK189254 in the human brain and to determine the dose–H3 receptor occupancy relationship of GSK189254. All scans were performed during normal waking hours when the levels of histamine in the brain are known to be stable (30).
After intravenous administration of 11C-GSK189254 in humans, the regional distribution of radioactivity in the brain was consistent with the known distribution of H3 receptors (6). Two key issues had to be addressed in the analysis. The first was the slow off-rate of the radioligand from the receptor. The constraints that were applied to the 2-tissue-compartment model relating tissue and plasma data assume that regional differences in 11C-GSK189254 kinetics are due to differences in delivery and available receptor density but that nondisplaceable binding and receptor off-rates are the same across regions. The resultant baseline estimates of VT were high in all regions, particularly in the putamen (98 ± 50) and caudate (74 ± 27). The stability of these estimates was investigated by successive removal of late–time frame data from the analyses. No apparent bias in low-, intermediate-, or high-binding regions was observed, indicating the suitability of the modeling approach.
The second issue was the lack of a suitable reference region devoid of specific binding. Estimates of VND, required for the calculation of occupancy after oral dosing, were therefore derived from occupancy plots. Linear plots, indicating occupancy to be the same across all ROIs after oral dosing of GSK189254, were obtained for all subjects, with consistent estimates of VND being obtained from the intercepts.
Given the estimate of VND, together with a value for the fP (from which the tissue free fraction could be obtained), it was then possible to obtain 2 essentially independent estimates of the in vivo affinity of GSK189254 for its receptor sites, assuming passive diffusion. The first, using a conventional estimate of the oral ED50, together with plasma pharmacokinetic data, gave an in vivo plasma EC50 of 0.011 nM. The second, from an analysis of a carryover effect in the test–retest data, together with the known tracer specific activities, gave a tissue KD of 0.0095 nM. If there was also carryover in the postdose scans, simulations (not shown) indicated that this carryover would be negligible at high occupancies and that the occupancies quoted in this study would be slightly underestimated.
The affinities measured in vivo using PET, equivalent to a negative log of the dissociation constant (pK) of about 11, is an order of magnitude higher than that measured in vitro for the human H3 receptor (pK = 9.59–9.90) (20). In contrast, in a previous study in the pig, the estimate of the in vivo pK using PET (9.9) was similar to the in vitro estimate, both being an order of magnitude lower than that seen here in humans, with the tracer kinetics being more obviously reversible in vivo (19). The estimated receptor densities are, however, similar in the 2 species. The reason for this interspecies difference does not appear to be that GSK189254 is a P-glycoprotein substrate in the pig, because brain uptake of this compound in vivo at tracer doses was not affected by cyclosporine A (C. Salinas, unpublished data, 2009). Differences in the concentration of endogenous histamine are unlikely to contribute to the observed interspecies difference, and we therefore do not have an explanation for this difference in the affinity of GSK189254.
CONCLUSION
In general, care must be taken when using high-affinity PET radioligands to maintain a negligible level of receptor occupancy by the radioligand itself during the scan acquisition. Typically for an ideal PET radioligand, the mass dose associated with its injection should not occupy more than about 5% of the available receptors. This mass dose limitation needs to be considered in future study designs if 11C-GSK189254 is to be used to measure H3 receptor occupancy of other compounds in humans. Suitable time intervals between PET acquisitions will be required to allow for the washout of the mass dose associated with baseline scans before postdose scans. The ED50 estimated here would correspond to a dose of no more than 0.003 μg/kg of body weight to achieve less than 5% occupancy. Alternatively, a higher tracer mass dose effect could be accounted for by modeling as shown here, particularly if relative biases between baseline and posttreatment scans are minimized by the administration of similar tracer mass doses in the 2 scans.
APPENDIX
A 2-Tissue-Compartment Model Incorporating Dispersion of Input Function
Preliminary investigation of the data showed that a reversible 2-tissue-compartment (4k) model was necessary to fit low-binding regions but gave unreasonable estimates of K1, indicating a dispersion effect. The operational equation of a catenary 2-tissue-compartment model with a true arterial plasma input function (CA(t)) is given by (29):Eq. 1Awhere t is time; CT(t) the total tissue activity; K1 the plasma-to-tissue clearance; and Φ1, Φ2, θ1, and θ2 are functions of the rate constants k2, k3, and k4, describing the exchange of tracer between plasma and tissue compartments (27). ⊗denotes convolution.
Assuming a simple exponential dispersion function (31), then the measured input function (CM (t)) is given by:Eq. 2ASubstitution of CA(t) by CM(t) in Equation 1A then gives:Eq. 3Awhere Φ3 and Φ4 are functions of k2, k3, k4, and kdisp. The particular formulation of these functions is equivalent to those of a full reference-tissue model treating CM(t) as the reference input (29), with VND equaling 1 and turnover rate constant equaling kdisp.
Sharing of Common Parameters Across Regions Within Scans
The above model involves 5 identifiable parameters, K1, k2, k3, k4, and kdisp. Estimating these parameters independently for each ROI gave acceptable fits and residuals (results not shown) but frequently with unreasonable numeric estimates of the individual parameters (e.g., negative values) or outlying large estimates of VT for high-binding regions. Constraints were therefore applied by minimizing the residual sum of squares simultaneously across all regions, with individual estimates of K1 and k3 for each region but with common (global) estimates of K1/k2, k4, and kdisp.
Estimation of Tracer Mass Dose Effect and of Carryover Between Test–Retest Scans
The aim of the modeling was to obtain estimates of KD and Bmax from test–retest data, using a simple competition model.
Parallel 2-tissue-compartment models were set up (32) describing the kinetics of 11C-GSK189254 and of unlabeled GSK189254, coupling the test–retest scans and taking into account the different specific activities of the injected doses and the carryover of mass dose from the first to the second scan. Carryover of radioactivity between scans was assumed to be negligible, because of decay, but the initial-condition concentrations of unlabeled GSK189254 in tissue compartments at the start of the retest scan were calculated by extrapolation of the fits to the test scan. This was achieved by applying a single exponential extrapolation to the tail of the input functions over the approximately 30-min period between the end of the first scan and the beginning of the second and calculating the corresponding compartmental concentrations as initial conditions for the second scan. The input functions were precorrected for dispersion with the mean estimate of kdisp obtained from fits to individual scans using Equation 3A.
k3 varies in time throughout the 2 scans according to the following equation:Eq. 4Awhere DS(t) is the concentration of the unlabeled GSK189254 in the specifically bound compartment. Bmax, fND, kon, and k4 (=koff) are time-invariant.
The model was parameterized in terms of K1, VND, koff, , and Bmax, whereEq. 5AEquations were integrated numerically. VND was fixed (see main text), and the remaining parameters were estimated using a Levenberg–Marquardt algorithm. Under the assumption of passive diffusion, an estimate of fND was derived from the measured fP:Eq. 6AThis then allows for the true KD to be derived asEq. 7AIndividual estimates of K1 were obtained for each region and scan, and estimates of Bmax were obtained for each region across test and retest scans. KD and koff were treated as common parameters (see above) within subjects across regions and across the 2 scans.
Acknowledgments
We thank staff of CAMH who carried out the acquisition of the data and provided care for the subjects while they underwent the procedures.
Footnotes
-
COPYRIGHT © 2010 by the Society of Nuclear Medicine, Inc.
References
- Received for publication October 20, 2009.
- Accepted for publication January 12, 2010.