Acute induction of anxiety in humans by delta-9-tetrahydrocannabinol related to amygdalar cannabinoid-1 (CB1) receptors.

Use of Cannabis, the most widely used illicit drug worldwide, is associated with acute anxiety, and anxiety disorders following regular use. The precise neural and receptor basis of these effects have not been tested in man. Employing a combination of functional MRI (fMRI) and positron emission tomography (PET), we investigated whether the effects of delta-9-tetrahydrocannabinol (delta-9-THC), the main psychoactive ingredient of cannabis, on anxiety and on amygdala response while processing fearful stimuli were related to local availability of its main central molecular target, cannabinoid-1 (CB1) receptors in man. Fourteen healthy males were studied with fMRI twice, one month apart, following an oral dose of either delta-9-THC (10 mg) or placebo, while they performed a fear-processing task. Baseline availability of the CB1 receptor was studied using PET with [11C]MePPEP, a CB1 inverse agonist radioligand. Relative to the placebo condition, delta-9-THC induced anxiety and modulated right amygdala activation while processing fear. Both these effects were positively correlated with CB1 receptor availability in the right amygdala. These results suggest that the acute effects of cannabis on anxiety in males are mediated by the modulation of amygdalar function by delta-9-THC and the extent of these effects are related to local availability of CB1 receptors.

Test-retest reproducibility of quantitative binding measures of [11C]Ro15-4513, a PET ligand for GABAA receptors containing alpha5 subunits.

INTRODUCTION: Alteration of γ-aminobutyric acid "A" (GABAA) receptor-mediated neurotransmission has been associated with various neurological and psychiatric disorders. [11C]Ro15-4513 is a PET ligand with high affinity for α5-subunit-containing GABAA receptors, which are highly expressed in limbic regions of the human brain (Sur et al., 1998). We quantified the test-retest reproducibility of measures of [11C]Ro15-4513 binding derived from six different quantification methods (12 variants). METHODS: Five healthy males (median age 40 years, range 38-49 years) had a 90-min PET scan on two occasions (median interval 12 days, range 11-30 days), after injection of a median dose of 441 MegaBequerels of [11C]Ro15-4513. Metabolite-corrected arterial plasma input functions (parent plasma input functions, ppIFs) were generated for all scans. We quantified regional binding using six methods (12 variants), some of which were region-based (applied to the average time-activity curve within a region) and others were voxel-based: 1) Models requiring arterial ppIFs - regional reversible compartmental models with one and two tissue compartments (2kbv and 4kbv); 2) Regional and voxelwise Logan's graphical analyses (Logan et al., 1990), which required arterial ppIFs; 3) Model-free regional and voxelwise (exponential) spectral analyses (SA; (Cunningham and Jones, 1993)), which also required arterial ppIFs; 4) methods not requiring arterial ppIFs - voxelwise standardised uptake values (Kenney et al., 1941), and regional and voxelwise simplified reference tissue models (SRTM/SRTM2) using brainstem or alternatively cerebellum as pseudo-reference regions (Lammertsma and Hume, 1996; Gunn et al., 1997). To compare the variants, we sampled the mean values of the outcome parameters within six bilateral, non-reference grey matter regions-of-interest. Reliability was quantified in terms of median absolute percentage test-retest differences (MA-TDs; preferentially low) and between-subject coefficient of variation (BS-CV, preferentially high), both compounded by the intraclass correlation coefficient (ICC). These measures were compared between variants, with particular interest in the hippocampus. RESULTS: Two of the six methods (5/12 variants) yielded reproducible data (i.e. MA-TD <10%): regional SRTMs and voxelwise SRTM2s, both using either the brainstem or the cerebellum; and voxelwise SA. However, the SRTMs using the brainstem yielded a lower median BS-CV (7% for regional, 7% voxelwise) than the other variants (8-11%), resulting in lower ICCs. The median ICCs across six regions were 0.89 (interquartile range 0.75-0.90) for voxelwise SA, 0.71 (0.64-0.84) for regional SRTM-cerebellum and 0.83 (0.70-0.86) for voxelwise SRTM-cerebellum. The ICCs for the hippocampus were 0.89 for voxelwise SA, 0.95 for regional SRTM-cerebellum and 0.93 for voxelwise SRTM-cerebellum. CONCLUSION: Quantification of [11C]Ro15-4513 binding shows very good to excellent reproducibility with SRTM and with voxelwise SA which, however, requires an arterial ppIF. Quantification in the α5 subunit-rich hippocampus is particularly reliable. The very low expression of the α5 in the cerebellum (Fritschy and Mohler, 1995; Veronese et al., 2016) and the substantial α1 subunit density in this region may hamper the application of reference tissue methods.

Using [(11)C]Ro15 4513 PET to characterise GABA-benzodiazepine receptors in opiate addiction: Similarities and differences with alcoholism.

The importance of the GABA-benzodiazepine receptor complex and its subtypes are increasingly recognised in addiction. Using the α1/α5 benzodiazepine receptor PET radioligand [(11)C]Ro15 4513, we previously showed reduced binding in the nucleus accumbens and hippocampus in abstinent alcohol dependence. We proposed that reduced [(11)C]Ro15 4513 binding in the nucleus accumbens was a marker of addiction whilst the reduction in hippocampus and positive relationship with memory was a consequence of chronic alcohol abuse. To examine this further we assessed [(11)C]Ro15 4513 binding in another addiction, opiate dependence, and used spectral analysis to estimate contributions of α1 and α5 subtypes to [(11)C]Ro15 4513 binding in opiate and previously acquired alcohol-dependent groups. Opiate substitute maintained opiate-dependent men (n=12) underwent an [(11)C]Ro15 4513 PET scan and compared with matched healthy controls (n=13). We found a significant reduction in [(11)C]Ro15 4513 binding in the nucleus accumbens in the opiate-dependent compared with the healthy control group. There was no relationship between [(11)C]Ro15 4513 binding in the hippocampus with memory. We found that reduced [(11)C]Ro15 4513 binding was associated with reduced α5 but not α1 subtypes in the opiate-dependent group. This was also seen in an alcohol-dependent group where an association between memory performance and [(11)C]Ro15 4513 binding was primarily driven by α5 and not α1 subtype. We suggest that reduced α5 levels in the nucleus accumbens are associated with addiction since we have now shown this in dependence to two pharmacologically different substances, alcohol and opiates.

Test-retest reproducibility of cannabinoid-receptor type 1 availability quantified with the PET ligand [¹¹C]MePPEP.

BACKGROUND: Endocannabinoids are involved in normal cognition, and dysfunction in cannabinoid-receptor-mediated neurotransmission has been suggested in a variety of neurological and psychiatric pathologies. The type 1 cannabinoid receptor (CB1) is widely expressed in the human central nervous system. The objective of this study was to quantify the test-retest reproducibility of measures of the PET ligand [(11)C]MePPEP in order to assess the stability of CB1-receptor quantification in humans in vivo. METHODS: Fifteen healthy subjects (eight females; median age 32 years, range 25 to 65 years) had a 90-minute PET scan on two occasions after injection of a median dose of [(11)C]MePPEP of 364 MBq. Metabolite-corrected arterial plasma input functions were obtained for all scans. Eight ROIs, reflecting different levels of receptor densities/concentrations, were defined automatically: hippocampus, anterior cingulate gyrus, inferior frontal gyrus, caudate nucleus, globus pallidus, nucleus accumbens, thalamus, and pons. We used seven quantification methods: reversible compartmental models with one and two tissue classes, two and four rate constants, and a variable blood volume term (2kbv; 4kbv); model-free (spectral) analyses with and without regularisation, including one with voxel-wise quantification; the simplified reference tissue model (SRTM) with pons as a pseudo-reference region; and modified standard uptake values (mSUVs) calculated for the period of ~30-60 min after injection. Percentage test-retest change and between-subject variability were both assessed, and test-retest reliability was quantified by the intraclass correlation coefficient (ICC). The ratio of binding estimates pallidum:pons served as an indicator of a method's ability to reflect binding heterogeneity. RESULTS: Neither the SRTM nor the 4kbv model produced reliable measures, with ICCs around zero. Very good (>0.75) or excellent (>0.80) ICCs were obtained with the other methods. The most reliable were spectral analysis parametric maps (average across regions±standard deviation 0.83±0.03), rank shaping regularised spectral analysis (0.82±0.05), and the 2kbv model (0.82±0.09), but mSUVs were also reliable for most regions (0.79±0.13). Mean test-retest changes among the five well-performing methods ranged from 12±10% for mSUVs to 16% for 2kbv. Intersubject variability was high, with mean between-subject coefficients of variation ranging from 32±13% for mSUVs to 45% for 2kbv. The highest pallidum:pons ratios of binding estimates were achieved by mSUV (4.2), spectral analysis-derived parametric maps (3.6), and 2kbv (3.6). CONCLUSION: Quantification of CB1 receptor availability using [(11)C]MePPEP shows good to excellent reproducibility with several kinetic models and model-free analyses, whether applied on a region-of-interest or voxelwise basis. Simple mSUV measures were also reliable for most regions, but do not allow fully quantitative interpretation. [(11)C]MePPEP PET is well placed as a tool to investigate CB1-receptor mediated neurotransmission in health and disease.

PET imaging: implications for the future of therapy monitoring with PET/CT in oncology.

Among the methods based on molecular imaging, the measure of the tracer uptake variation between a baseline and follow-up scan with the SUV and [(18)F]FDG-PET/CT is a very powerful tool for assessing response to treatment in oncology. However, the development of new targeted therapeutics and tissue pharmacokinetic evaluation of existing ones are increasingly requiring therapy monitoring with alternative tracers and indicators. In parallel, the potential predictive and prognostic value of other image-derived parameters, such as tumour volume and textural features, relating to tumoral heterogeneity, has recently emerged from several works.

Imaging of cellular proliferation in liver metastasis by [18F]fluorothymidine positron emission tomography: effect of therapy.

Although [(18)F]fluorothymidine positron emission tomography (FLT-PET) permits estimation of tumor thymidine kinase-1 expression, and thus, cell proliferation, high physiological uptake of tracer in liver tissue can limit its utility. We evaluated FLT-PET combined with a temporal-intensity information-based voxel-clustering approach termed kinetic spatial filtering (FLT-PET(KSF)) for detecting drug response in liver metastases. FLT-PET and computed tomography data were collected from patients with confirmed breast or colorectal liver metastases before, and two weeks after the first cycle of chemotherapy. Changes in tumor FLT-PET and FLT-PET(KSF) variables were determined. Visual distinction between tumor and normal liver was seen in FLT-PET(KSF) images. Of the 33 metastases from 20 patients studied, 26 were visible after kinetic filtering. The net irreversible retention of the tracer (Ki; from unfiltered data) in the tumor, correlated strongly with tracer uptake when the imaging variable was an unfiltered average or maximal standardized uptake value, 60 min post-injection (SUV(60,av): r = 0.9, SUV(60,max): r = 0.7; p < 0.0001 for both) and occurrence of high intensity voxels derived from FLT-PET(KSF) (r = 0.7, p < 0.0001). Overall, a significant reduction in the imaging variables was seen in responders compared to non-responders; however, the two week time point selected for imaging was too early to allow prediction of long term clinical benefit from chemotherapy. FLT-PET and FLT-PET(KSF) detected changes in proliferation in liver metastases.

Evaluation of limited blood sampling population input approaches for kinetic quantification of [18F]fluorothymidine PET data.

BACKGROUND: Quantification of kinetic parameters of positron emission tomography (PET) imaging agents normally requires collecting arterial blood samples which is inconvenient for patients and difficult to implement in routine clinical practice. The aim of this study was to investigate whether a population-based input function (POP-IF) reliant on only a few individual discrete samples allows accurate estimates of tumour proliferation using [18F]fluorothymidine (FLT). METHODS: Thirty-six historical FLT-PET data with concurrent arterial sampling were available for this study. A population average of baseline scans blood data was constructed using leave-one-out cross-validation for each scan and used in conjunction with individual blood samples. Three limited sampling protocols were investigated including, respectively, only seven (POP-IF7), five (POP-IF5) and three (POP-IF3) discrete samples of the historical dataset. Additionally, using the three-point protocol, we derived a POP-IF3M, the only input function which was not corrected for the fraction of radiolabelled metabolites present in blood. The kinetic parameter for net FLT retention at steady state, Ki, was derived using the modified Patlak plot and compared with the original full arterial set for validation. RESULTS: Small percentage differences in the area under the curve between all the POP-IFs and full arterial sampling IF was found over 60 min (4.2%-5.7%), while there were, as expected, larger differences in the peak position and peak height.A high correlation between Ki values calculated using the original arterial input function and all the population-derived IFs was observed (R2 = 0.85-0.98). The population-based input showed good intra-subject reproducibility of Ki values (R2 = 0.81-0.94) and good correlation (R2 = 0.60-0.85) with Ki-67. CONCLUSIONS: Input functions generated using these simplified protocols over scan duration of 60 min estimate net PET-FLT retention with reasonable accuracy.

Characterisation of the contribution of the GABA-benzodiazepine α1 receptor subtype to [(11)C]Ro15-4513 PET images.

This positron emission tomography (PET) study aimed to further define selectivity of [(11)C]Ro15-4513 binding to the GABARα5 relative to the GABARα1 benzodiazepine receptor subtype. The impact of zolpidem, a GABARα1-selective agonist, on [(11)C]Ro15-4513, which shows selectivity for GABARα5, and the nonselective benzodiazepine ligand [(11)C]flumazenil binding was assessed in humans. Compartmental modelling of the kinetics of [(11)C]Ro15-4513 time-activity curves was used to describe distribution volume (V(T)) differences in regions populated by different GABA receptor subtypes. Those with low α5 were best fitted by one-tissue compartment models; and those with high α5 required a more complex model. The heterogeneity between brain regions suggested spectral analysis as a more appropriate method to quantify binding as it does not a priori specify compartments. Spectral analysis revealed that zolpidem caused a significant V(T) decrease (~10%) in [(11)C]flumazenil, but no decrease in [(11)C]Ro15-4513 binding. Further analysis of [(11)C]Ro15-4513 kinetics revealed additional frequency components present in regions containing both α1 and α5 subtypes compared with those containing only α1. Zolpidem reduced one component (mean±s.d.: 71%±41%), presumed to reflect α1-subtype binding, but not another (13%±22%), presumed to reflect α5. The proposed method for [(11)C]Ro15-4513 analysis may allow more accurate selective binding assays and estimation of drug occupancy for other nonselective ligands.

A [11C]Ro15 4513 PET study suggests that alcohol dependence in man is associated with reduced α5 benzodiazepine receptors in limbic regions.

Preclinical evidence suggests the α5 subtype of the GABA-benzodiazepine receptor is involved in some of the actions of alcohol and in memory. The positron emission tomography (PET) tracer, [(11)C]Ro15 4513 shows relative selectivity in labelling the α5 subtype over the other GABA-benzodiazepine receptor subtypes in limbic regions of the brain. We used this tracer to investigate the distribution of α5 subtype availability in human alcohol dependence and its relationship to clinical variables. Abstinent (>6 weeks) alcohol-dependent men and healthy male controls underwent an [(11)C]Ro15 4513 PET scan. We report [(11)C]Ro15 4513 brain uptake for 8 alcohol-dependent men and 11 healthy controls. We found a significant reduction in [(11)C]Ro15 4513 binding in the nucleus accumbens, parahippocampal gyri, right hippocampus and amygdala in the alcohol-dependent compared with the healthy control group. Levels of [(11)C]Ro15 4513 binding in both hippocampi were significantly and positively associated with performance on a delayed verbal memory task in the alcohol-dependent but not the control group. We speculate that the reduced limbic [(11)C]Ro15 4513 binding seen here results from the effects of alcohol, though we cannot currently distinguish whether they are compensatory in nature or evidence of brain toxicity.

Quantification of ligand PET studies using a reference region with a displaceable fraction: application to occupancy studies with [(11)C]-DASB as an example.

This paper aims to build novel methodology for the use of a reference region with specific binding for the quantification of brain studies with radioligands and positron emission tomography (PET). In particular: (1) we introduce a definition of binding potential BP(D)=DVR-1 where DVR is the volume of distribution relative to a reference tissue that contains ligand in specifically bound form, (2) we validate a numerical methodology, rank-shaping regularization of exponential spectral analysis (RS-ESA), for the calculation of BP(D) that can cope with a reference region with specific bound ligand, (3) we demonstrate the use of RS-ESA for the accurate estimation of drug occupancies with the use of correction factors to account for the specific binding in the reference. [(11)C]-DASB with cerebellum as a reference was chosen as an example to validate the methodology. Two data sets were used; four normal subjects scanned after infusion of citalopram or placebo and further six test-retest data sets. In the drug occupancy study, the use of RS-ESA with cerebellar input plus corrections produced estimates of occupancy very close the ones obtained with plasma input. Test-retest results demonstrated a tight linear relationship between BP(D) calculated either with plasma or with a reference input and high reproducibility.

Further evaluation of quantum chemical methods for the prediction of non-specific binding of positron emission tomography tracers.

The non-specific binding of candidate positron emission tomography (PET) radiotracers causes resulting PET images to have poor contrast and is a key determinant for the success or failure of imaging drugs. Non-specific binding is thought to arise when radiotracers bind to cell membranes and moieties other than their intended target. Our previous preliminary work has proposed the use of the drug-lipid interaction energy descriptor to predict the level of non-specific binding in vivo using a limited set of ten well known PET radiotracers with kinetic modelling data taken from the literature. This work validates and extends the use of the drug-lipid interaction energy descriptor using a new set of twenty-two candidate PET radiotracers with non-specific binding data recently collected at the same imaging centre with consistent methodology. As with the previous set of radiotracers, a significant correlation is found between the quantum chemical drug-lipid interaction energy and in vivo non-specific binding experimental values. In an effort to speed up the calculation process, several semi-empirical quantum chemical methods were assessed for their ability to reproduce the ab initio results. However no single semi-empirical method was found to consistently reproduce the level of correlation achieved with ab initio quantum chemical methods.

[18F]-3'Deoxy-3'-fluorothymidine positron emission tomography and breast cancer response to docetaxel.

PURPOSE: To establish biomarkers indicating clinical response to taxanes, we determined whether early changes in [(18)F]-3'deoxy-3'-fluorothymidine positron emission tomography (FLT-PET) can predict benefit from docetaxel therapy in breast cancer. EXPERIMENTAL DESIGN: This was a prospective unblinded study in 20 patients with American Joint Committee on Cancer (AJCC) stage II-IV breast cancer unresponsive to first-line chemotherapy or progressing on previous therapy. Individuals underwent a baseline dynamic FLT-PET scan followed by a scan 2 weeks after initiating the first or second cycle of docetaxel. PET variables were compared with anatomic response midtherapy (after 3 cycles). RESULTS: Average and maximum tumor standardized uptake values at 60 minutes (SUV(60,av) and SUV(60,max)) normalized to body surface area ranged between 1.7 and 17.0 and 5.6 and 26.9 × 10(-5) m(2)/mL, respectively. Docetaxel treatment resulted in a significant decrease in FLT uptake (P = 0.0003 for SUV(60,av) and P = 0.0002 for SUV(60,max)). Reduction in tumor SUV(60,av) was associated with target lesion size changes midtherapy (Pearson R for SUV(60,av) = 0.64; P = 0.004) and predicted midtherapy target lesion response (0.85 sensitivity and 0.80 specificity). Decreases in SUV(60,av) in responders were due, at least in part, to reduced net intracellular trapping of FLT (rate constant, K(i)). Docetaxel significantly reduced K(i) by 51.1% (±28.4%, P = 0.0009). CONCLUSION: Changes in tumor proliferation assessed by FLT-PET early after initiating docetaxel chemotherapy can predict lesion response midtherapy with good sensitivity warranting prospective trials to assess the ability to stop therapy in the event of non-FLT-PET response.

Reproducibility of 18F-FDG and 3'-deoxy-3'-18F-fluorothymidine PET tumor volume measurements.

UNLABELLED: The objective of this study was to establish the repeatability and reproducibility limits of several volume-related PET image-derived indices-namely tumor volume (TV), mean standardized uptake value, total glycolytic volume (TGV), and total proliferative volume (TPV)-relative to those of maximum standardized uptake value (SUV(max)), commonly used in clinical practice. METHODS: Fixed and adaptive thresholding, fuzzy C-means, and fuzzy locally adaptive Bayesian methodology were considered for TV delineation. Double-baseline (18)F-FDG (17 lesions, 14 esophageal cancer patients) and 3'-deoxy-3'-(18)F-fluorothymidine ((18)F-FLT) (12 lesions, 9 breast cancer patients) PET scans, acquired at a mean interval of 4 d and before any treatment, were used for reproducibility evaluation. The repeatability of each method was evaluated for the same datasets and compared with manual delineation. RESULTS: A negligible variability of less than 5% was measured for all segmentation approaches in comparison to manual delineation (5%-35%). SUV(max) reproducibility levels were similar to others previously reported, with a mean percentage difference of 1.8% +/- 16.7% and -0.9% +/- 14.9% for the (18)F-FDG and (18)F-FLT lesions, respectively. The best TV, TGV, and TPV reproducibility limits ranged from -21% to 31% and -30% to 37% for (18)F-FDG and (18)F-FLT images, respectively, whereas the worst reproducibility limits ranged from -90% to 73% and -68% to 52%, respectively. CONCLUSION: The reproducibility of estimating TV, mean standardized uptake value, and derived TGV and TPV was found to vary among segmentation algorithms. Some differences between (18)F-FDG and (18)F-FLT scans were observed, mainly because of differences in overall image quality. The smaller reproducibility limits for volume-derived image indices were similar to those for SUV(max), suggesting that the use of appropriate delineation tools should allow the determination of tumor functional volumes in PET images in a repeatable and reproducible fashion.

Kinetic filtering of [(18)F]Fluorothymidine in positron emission tomography studies.

[(18)F]Fluorothymidine (FLT) is a cell proliferation marker that undergoes predominantly hepatic metabolism and therefore shows a high level of accumulation in the liver, as well as in rapidly proliferating tumours. Furthermore, the tracer's uptake is substantial in other organs including the heart. We present a nonlinear kinetic filtering technique which enhances the visualization of tumours imaged with FLT positron emission tomography (FLT-PET). A classification algorithm to isolate cancerous tissue from healthy organs was developed and validated using 29 scan data from patients with locally advanced or metastatic breast cancer. A large reduction in signal from the liver and heart of 80% was observed following application of the kinetic filter, whilst the majority of signal from both primary tumours and metastases was retained. A scan acquisition time of 60 min has been shown to be sufficient to obtain the necessary kinetic data. The algorithm extends utility of FLT-PET imaging in oncology research.

Further evaluation of the carbon11-labeled D(2/3) agonist PET radiotracer PHNO: reproducibility in tracer characteristics and characterization of extrastriatal binding.

[(11)C]-(+)-PHNO is a new dopamine D(2/3) receptor agonist radiotracer which has been successfully used to measure D(2/3) receptor availability in experimental animals and man. Here we report in vivo evaluation in the rat of the biodistribution, metabolism, specificity, selectivity, and dopamine sensitivity of carbon11-labeled PHNO ([(11)C]-3-PHNO) produced by an alternative radiochemical synthesis method. [(11)C]-3-PHNO showed rapid metabolism and clearance from most peripheral organs and tissues. [(11)C]-3-PHNO, but not its polar metabolite, readily crossed the blood-brain barrier and showed high levels of uptake in the D(2/3)-rich striatum. Pretreatment with unlabeled PHNO and the D(2/3) receptor antagonist raclopride indicated that binding in the striatum was specific and selective to D(2/3) receptors. PET studies in anesthetized rats revealed significant reductions in [(11)C]-3-PHNO binding in the striatum following amphetamine administration, indicating sensitivity to increases in endogenous dopamine concentrations. D(2/3) antagonist pretreatment additionally indicated moderate levels of [(11)C]-3-PHNO specific binding in several extrastriatal brain areas-most notably the olfactory bulbs and tubercles, thalamus, and hypothalamus. Of particular interest, approximately 30% of [(11)C]-3-PHNO signal in the cerebellum-a region often used as a "low-binding" reference region for PET quantification-was attributable to specific signal. These data demonstrate that [(11)C]-3-PHNO shows similar tracer characteristics to [(11)C]-(+)-PHNO, but additionally indicate that radiolabeled PHNO may be used to estimate D(2/3) receptor availability in select extrastriatal brain regions with PET.

A new model for prediction of drug distribution in tumor and normal tissues: pharmacokinetics of temozolomide in glioma patients.

Difficulties in direct measurement of drug concentrations in human tissues have hampered the understanding of drug accumulation in tumors and normal tissues. We propose a new system analysis modeling approach to characterize drug distribution in tissues based on human positron emission tomography (PET) data. The PET system analysis method was applied to temozolomide, an important alkylating agent used in the treatment of brain tumors, as part of standard temozolomide treatment regimens in patients. The system analysis technique, embodied in the convolution integral, generated an impulse response function that, when convolved with temozolomide plasma concentration input functions, yielded predicted normal brain and brain tumor temozolomide concentration profiles for different temozolomide dosing regimens (75-200 mg/m(2)/d). Predicted peak concentrations of temozolomide ranged from 2.9 to 6.7 microg/mL in human glioma tumors and from 1.8 to 3.7 microg/mL in normal brain, with the total drug exposure, as indicated by the tissue/plasma area under the curve ratio, being about 1.3 in tumor compared with 0.9 in normal brain. The higher temozolomide exposures in brain tumor relative to normal brain were attributed to breakdown of the blood-brain barrier and possibly secondary to increased intratumoral angiogenesis. Overall, the method is considered a robust tool to analyze and predict tissue drug concentrations to help select the most rational dosing schedules.

Ab initio computational study of positron emission tomography ligands interacting with lipid molecule for the prediction of nonspecific binding.

Nonspecific binding is a poorly understood biological phenomenon of relevance in the study of small molecules interactions in vivo and in drug development. Nonspecific binding is thought to be correlated in part to a molecule's lipophilicity, typically estimated by measuring (or calculating) octanol-water partition coefficient. This is, however, a gross simplification of a complex phenomenon. In this article, we present a computational method whose aim is to help identify positron emission tomography (PET) ligands with low nonspecific binding characteristics by investigating the molecular basis of ligand-membrane interaction. We considered a set consisting of 10 well-studied central nervous system PET radiotracers acting on a variety of molecular targets. Quantum mechanical calculations were used to estimate the strength of the interaction between each drug molecule and one phospholipid molecule commonly present in mammalian membranes. The results indicate a correlation between the computed drug-lipid interaction energy and the in vivo nonspecific distribution volume relative to the free tracer plasma concentration, calculated using standard compartmental modeling for the analysis of PET data. Significantly, the drugs whose interaction with the lipid molecule more favorably possessed, in general, a higher nonspecific binding value, whereas for the drugs taken in consideration in this study, the water-octanol partition coefficient, log P, did not show good predictive power of the nonspecific binding. This study also illustrates how ab initio chemical methods may offer meaningful and unbiased insights for the understanding of the underlying chemical mechanisms in biological systems.

Structure and dynamics of phospholipid bilayers using recently developed general all-atom force fields.

Two fully hydrated pure-species phospholipids bilayers, 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) and 1,2-dioleoyl-sn-glycero-3-phosphorylcholine (DOPC), in the fluid phase and explicit solvent have been studied using molecular dynamics simulation. Atom interactions were modeled using recently developed force fields based on AMBER with full atomistic details. Several representative liquid phase properties for the structure and dynamics of lipids with different length of hydrocarbon chains and different level of saturation have been reproduced without artificially biasing the system in order to match experimental data. In particular, as the new GAFF (General Amber Force Field) has not been explicitly developed to reproduce lipid characteristics and is naturally compatible with standard AMBER nucleic acids and proteins parameters, it is here proven a promising tool to study mixed lipid-protein processes as protein activity dependence on membrane composition, permeation of solute across membranes, and other cellular processes.

Efficient and precise solvation free energies via alchemical adiabatic molecular dynamics.

A new molecular dynamics method for calculating free energies associated with transformations of the thermodynamic state or chemical composition of a system (also known as alchemical transformations) is presented. The new method extends the adiabatic dynamics approach recently introduced by Rosso et al. [J. Chem. Phys. 116, 4389 (2002)] and is based on the use of an additional degree of freedom, lambda, that is used as a switching parameter between the potential energy functions that characterize the two states. In the new method, the coupling parameter lambda is introduced as a fictitious dynamical variable in the Hamiltonian, and a system of switching functions is employed that leads to a barrier in the lambda free energy profile between the relevant thermodynamic end points. The presence of such a barrier, therefore, enhances sampling in the end point (lambda = 0 and lambda = 1) regions which are most important for computing relevant free energy differences. In order to ensure efficient barrier crossing, a high temperature T(lambda) is assigned to lambda and a fictitious mass m(lambda) is introduced as a means of creating an adiabatic separation between lambda and the rest of the system. Under these conditions, it is shown that the lambda free energy profile can be directly computed from the adiabatic probability distribution function of lambda without any postprocessing or unbiasing of the output data. The new method is illustrated on two model problems and in the calculation of the solvation free energy of amino acid side-chain analogs in TIP3P water. Comparisons to previous work using thermodynamic integration and free energy perturbation show that the new lambda adiabatic free energy dynamics method results in very precise free energy calculations using significantly shorter trajectories.

Mapping the backbone dihedral free-energy surfaces in small peptides in solution using adiabatic free-energy dynamics.

Ramachandran surfaces for the alanine di- and tripeptides in gas phase and solution are mapped out using the recently introduced adiabatic free-energy dynamics (AFED) approach introduced by Rosso et al. (J. Chem. Phys. 2002, 116, 4389) as applied to the CHARMM22 force field. It is shown that complete surfaces can be mapped out with an order of magnitude of greater efficiency with the AFED approach than they can using the popular umbrella sampling method. In the alanine dipeptide, it is found, in agreement with numerous other studies using the CHARMM22 force field, that the lowest free-energy structure is the extended beta conformation, (phi, psi) = (-81, 81), while in solution, the extended beta, (phi, psi) = (-81, 153) and right-handed alpha-helical, (phi, psi) = (-81, 63) conformations are nearly isoenergetic. In solution, a secondary minimum at (phi, psi) = (63, -81), corresponding to a C(7)ax conformation, occurs approximately 2.3 kcal/mol above the global free-energy minimum. The alanine tripeptide, a system that has received considerably less attention in the literature, is found to exhibit a similar structure to the alanine dipeptide with the extended beta conformation being the free-energy minimum in the gas phase and the beta and right-handed alpha-helical conformations being isoenergetic in solution. These studies indicate that the AFED method can be a powerful tool for studying multidimensional free-energy surfaces in complex systems.