Modeling nicotine arterial-venous differences to predict arterial concentrations and input based on venous measurements: application to smokeless tobacco and nicotine gum

Sampling Site Has a Critical Impact on Physiologically Based Pharmacokinetic Modeling

It has been shown that arterial (central) and venous (peripheral) plasma drug concentrations can be very different. While pharmacokinetic studies typically measure drug concentrations from the peripheral vein such as the arm vein, physiologically based pharmacokinetic (PBPK) models generally output simulated concentrations from the central venous compartment that physiologically represents the right atrium, a merge of the superior and inferior vena cava. In this study, a physiologically based peripheral forearm sampling site model was developed and verified using nicotine, ketamine, lidocaine, and fentanyl as model drugs. This verified model allows output of simulated peripheral venous concentrations that can be meaningfully compared with observed pharmacokinetic data from the arm vein. The generalized effect of PBPK model sampling site on simulation output was investigated. Drugs and metabolites with large volumes of distribution showed considerable concentration discrepancy between the simulated central venous compartment and the peripheral arm vein after intravenous or oral administration, resulting in significant differences in values for C _max and time taken to reach C _max (t _max ) In addition, the simulated central venous metabolite profile showed an unexpected profile that was not observed in the peripheral arm vein. Using fentanyl as a model compound, we show that using the wrong sampling site in PBPK models can lead to biased model evaluation and subsequent erroneous model parameter optimization. Such an error in model parameters along with the discrepant sampling site could dramatically mislead the pharmacokinetic prediction in unstudied clinical scenarios, affecting the assessment of drug safety and efficacy. Overall, this study shows that PBPK model publications should specify the model sampling sites and match them with those employed in clinical studies. SIGNIFICANCE STATEMENT: Our study shows that sampling from the central venous compartment (right atrium) during physiologically based pharmacokinetic model development gives rise to biased model evaluation and erroneous model parameterization when observed data are collected from the peripheral arm vein. This can lead to a clinically significant error in predictions of plasma concentration-time profiles in unstudied scenarios. To address this error, we developed and verified a novel peripheral sampling site model to simulate arm vein drug concentrations that can be applied to different drug dosing scenarios.

Are Physiologically Based Pharmacokinetic Models Reporting the Right C(max)? Central Venous Versus Peripheral Sampling Site

Physiologically based pharmacokinetic (PBPK) models can over-predict maximum plasma concentrations (C_max) following intravenous administration. A proposed explanation is that invariably PBPK models report the concentration in the central venous compartment, rather than the site where the samples are drawn. The purpose of this study was to identify and validate potential corrective models based on anatomy and physiology governing the blood supply at the site of sampling and incorporate them into a PBPK platform. Four models were developed and scrutinised for their corrective potential. All assumed the peripheral sampling site concentration could be described by contributions from surrounding tissues and utilised tissue-specific concentration-time profiles reported from the full-PBPK model within the Simcyp Simulator. Predicted concentrations for the peripheral site were compared to the observed C_max. The models results were compared to clinical data for 15 studies over seven compounds (alprazolam, imipramine, metoprolol, midazolam, omeprazole, rosiglitazone and theophylline). The final model utilised tissue concentrations from adipose, skin, muscle and a contribution from artery. Predicted C_max values considering the central venous compartment can over-predict the observed values up to 10-fold whereas the new sampling site predictions were within 2-fold of observed values. The model was particularly relevant for studies where traditional PBPK models over-predict early time point concentrations. A successful corrective model for C_max prediction has been developed, subject to further validation. These models can be enrolled as built-up modules into PBPK platforms and potentially account for factors that may affect the initial mixing of the blood at the site of sampling.

Arterial input function derived from pairwise correlations between PET-image voxels

A metabolite corrected arterial input function is a prerequisite for quantification of positron emission tomography (PET) data by compartmental analysis. This quantitative approach is also necessary for radioligands without suitable reference regions in brain. The measurement is laborious and requires cannulation of a peripheral artery, a procedure that can be associated with patient discomfort and potential adverse events. A non invasive procedure for obtaining the arterial input function is thus preferable. In this study, we present a novel method to obtain image-derived input functions (IDIFs). The method is based on calculation of the Pearson correlation coefficient between the time-activity curves of voxel pairs in the PET image to localize voxels displaying blood-like behavior. The method was evaluated using data obtained in human studies with the radioligands [(11)C]flumazenil and [(11)C]AZ10419369, and its performance was compared with three previously published methods. The distribution volumes (VT) obtained using IDIFs were compared with those obtained using traditional arterial measurements. Overall, the agreement in VT was good (∼3% difference) for input functions obtained using the pairwise correlation approach. This approach performed similarly or even better than the other methods, and could be considered in applied clinical studies. Applications to other radioligands are needed for further verification.

Fractional dynamics pharmacokinetics-pharmacodynamic models

While an increasing number of fractional order integrals and differential equations applications have been reported in the physics, signal processing, engineering and bioengineering literatures, little attention has been paid to this class of models in the pharmacokinetics-pharmacodynamic (PKPD) literature. One of the reasons is computational: while the analytical solution of fractional differential equations is available in special cases, it this turns out that even the simplest PKPD models that can be constructed using fractional calculus do not allow an analytical solution. In this paper, we first introduce new families of PKPD models incorporating fractional order integrals and differential equations, and, second, exemplify and investigate their qualitative behavior. The families represent extensions of frequently used PK link and PD direct and indirect action models, using the tools of fractional calculus. In addition the PD models can be a function of a variable, the active drug, which can smoothly transition from concentration to exposure, to hyper-exposure, according to a fractional integral transformation. To investigate the behavior of the models we propose, we implement numerical algorithms for fractional integration and for the numerical solution of a system of fractional differential equations. For simplicity, in our investigation we concentrate on the pharmacodynamic side of the models, assuming standard (integer order) pharmacokinetics.

Exposure to domestic violence associated with adult smoking in India: a population based study

OBJECTIVE

To investigate the relation between domestic violence and tobacco use among adults in India.

DESIGN

Multilevel cross sectional analyses of a nationally representative population based sample from the 1998-9 Indian national family health survey.

PARTICIPANTS

278,977 individuals aged 15 or older; and 89,092 ever married women aged 15-49.

MAIN OUTCOME

Dichotomous variables for smoking and chewing tobacco.

RESULTS

Women who reported being abused more than one year ago and those who reported being abused in the past year were more likely to smoke and chew tobacco than women who have never experienced domestic violence. Compared to individuals who lived in homes where no abuse was reported, those who lived in homes where a woman reported experiencing domestic violence were more likely to smoke and chew tobacco.

CONCLUSION

Domestic violence is associated with higher odds of smoking and chewing tobacco in India. Efforts to control tobacco use need to consider the larger psychosocial circumstances within which individuals who practise such harmful health behaviours reside.

Effect of smokeless tobacco product marketing and use on population harm from tobacco use policy perspective for tobacco-risk reduction

This article presents policy perspectives on the marketing of smokeless tobacco products to reduce population harm from tobacco use. Despite consensus that smokeless tobacco products as sold in the United States are less dangerous than cigarettes, there is no consensus on how to proceed. Diverse factions have different policy concerns. While the tobacco industry is exempted from U.S. Food and Drug Administration (FDA) oversight, the pharmaceutical industry whose nicotine replacement therapy (NRT) medicines compete with smokeless tobacco as noncombustible nicotine-delivery systems are regulated by the FDA. Some public health experts support smokeless tobacco use to reduce population harm from tobacco; other public health experts oppose promoting smokeless tobacco for harm reduction. Adult consumers can freely purchase currently-marketed smokeless tobacco products and even more-deadly cigarettes. Concerns with and advantages of smokeless tobacco products are discussed. In that noncombustible medicinal nicotine-delivery systems have been proven to be effective smoking-cessation aids, smokeless tobacco, as another source of psychoactive doses of nicotine, could be used similarly, in a dose-response fashion as a smoking-cessation aid (consistent with FDA principles for evaluating generic versions of drugs). Price measures should be used on tobacco products to make costs to consumers proportional to product health risks (which would make smokeless tobacco much cheaper than cigarettes), and smokeless tobacco should be encouraged as an option for smoking cessation in adult smokers, particularly for those who have failed to stop smoking using NRT or other methods.

Possible role of ammonia on the deposition, retention, and absorption of nicotine in humans while smoking

This perspective presents an overview of the properties of tobacco smoke aerosol and the possible effect of ammonia on the deposition location, retention and the amount and rate of nicotine absorption during cigarette smoking. Three main mechanisms describe the absorption of smoke constituents: (A) gas-phase constituents deposit directly; (B) particles deposit and the constituents then diffuse through the particle into the biological buffer and then into the tissue; and (C) particulate phase constituents evaporate from the particles and then deposit from the gas phase. Nicotine from smoking deposits and is absorbed predominately in the lungs. When particles deposit on the lung-blood interfaces, nicotine is absorbed rapidly, regardless of the acid-base nature of the particles. This is due to the buffering capacity of the lung-blood interfaces and the small mass of nicotine per puff distributed over a large number of particles depositing onto a huge lung surface. The composition of both tobacco smoke aerosol particles and the gas phase are time dependent. Ammonia in mainstream smoke evaporates faster from particles than nicotine. It is, therefore, unlikely that ammonia can significantly affect the volatility of MS smoke nicotine from particles in the smoke aerosol. It is certain that no single measurement of tobacco or of smoke, especially one made under equilibrium conditions, can adequately characterize the time-dependent properties of mainstream smoke aerosol. Thus, the fraction of nonprotonated freebase nicotine in trapped, aged smoke particulate matter has not been shown to be a useful predictor of the amount or total rate of nicotine uptake in human smokers. Similarly, "smoke pH" and "pHeff" are not useful practical parameters for providing understanding or predictability of tobacco smoke chemistry or nicotine bioavailability.

Predicting brain concentrations of drug using positron emission tomography and venous input: modeling of arterial-venous concentration differences

OBJECTIVE

In a positron emission tomography (PET) study, the concentrations of the labeled drug (radiotracer) are often different in arterial and venous plasma, especially immediately following administration. In a PET study, the transfer of the drug from plasma to brain is usually described using arterial plasma concentrations, whereas venous sampling is standard in clinical pharmacokinetic studies of new drug candidates. The purpose of the study was to demonstrate the modeling of brain drug kinetics based on PET data in combination with venous blood sampling and an arterio-venous transform (T(av)).

METHODS

Brain kinetics (C(br)) was described as the convolution of arterial plasma kinetics (C(ar)) with an arterial-to-brain impulse response function (T(br)). The arterial plasma kinetics was obtained as venous plasma kinetics (C(ve)) convolved with the inverse of the arterio-venous transform (T(av) (-1)). The brain kinetics was then given by C(br)=C(ve)*T(av) (-1)*T(br). This concept was applied on data from a clinical PET study in which both arterial and venous plasma sampling was done in parallel to PET measurement of brain drug kinetics. The predictions of the brain kinetics based on an arterial input were compared with predictions using a venous input with and without an arterio-venous transform.

RESULTS

The venous based models for brain distribution, including a biexponential arterio-venous transform, performed comparably to models based on arterial data and better than venous based models without the transform. It was also shown that three different brain regions with different shaped concentration curves could be modeled with a common arterio-venous transform together with an individual brain distribution model.

CONCLUSION

We demonstrated the feasibility of modeling brain drug kinetics based on PET data in combination with venous blood sampling and an arterio-venous transform. Such a model can in turn be used for the calculation of brain kinetics resulting from an arbitrary administration mode by applying this model on venous plasma pharmacokinetics. This would be an important advantage in the development of drugs acting in the brain, and in other circumstances when the effect is likely to be closer related to the brain than the plasma concentration.

Blood-brain barrier penetration of zolmitriptan--modelling of positron emission tomography data

Positron emission tomography (PET) with the drug radiolabelled allows a direct measurement of brain or other organ kinetics, information which can be essential in drug development. Usually, however, a PET-tracer is administered intravenously (i.v.), whereas the therapeutic drug is mostly given orally or by a different route to the PET-tracer. In such cases, a recalculation is needed to make the PET data representative for the alternative administration route. To investigate the blood-brain barrier penetration of a drug (zolmitriptan) using dynamic PET and by PK modelling quantify the brain concentration of the drug after the nasal administration of a therapeutic dose. [11C]Zolmitriptan at tracer dose was administered as a short i.v. infusion and the brain tissue and venous blood kinetics of [11C]zolmitriptan was measured by PET in 7 healthy volunteers. One PET study was performed before and one 30 min after the administration of 5 mg zolmitriptan as nasal spray. At each of the instances, the brain radioactivity concentration after subtraction of the vascular component was determined up to 90 min after administration and compared to venous plasma radioactivity concentration after correction for radiolabelled metabolites. Convolution methods were used to describe the relationship between arterial and venous tracer concentrations, respectively between brain and arterial tracer concentration. Finally, the impulse response functions derived from the PET studies were applied on plasma PK data to estimate the brain zolmitriptan concentration after a nasal administration of a therapeutic dose. The studies shows that the PET data on brain kinetics could well be described as the convolution of venous tracer kinetics with an impulse response including terms for arterial-to-venous plasma and arterial-to-brain impulse responses. Application of the PET derived impulse responses on the plasma PK from nasal administration demonstrated that brain PK of zolmitriptan increased with time, achieving about 0.5 mg/ml at 30 min and close to a maximum of 1.5 mg/ml after 2 hr. A significant brain concentration was observed already after 5 min. The data support the notation of a rapid brain availability of zolmitriptan after nasal administration.

Intraportal nicotine infusion in rats decreases hepatic blood flow through endothelin-1 and both endothelin A and endothelin B receptors

Smoking has been demonstrated to aggravate liver injury. Nicotine, a major pharmacological component of tobacco smoke, affects a multitude of functions. Smoking and nicotine induce synthesis of endothelin (ET)-1. The effect of intraportal infusion of nicotine on hepatic circulation and an involvement of ET-1 and ET receptor in the action of nicotine were investigated in rats. Nicotine (0-100 microg/kg/h) was infused into the portal vein of urethane-anesthetized rats, and changes of hepatic blood flow were evaluated. Intraportal infusion of nicotine dose-dependently decreased hepatic blood flow and increased portal pressure without any alteration of heart rate or arterial blood pressure. This action of intraportal nicotine was completely abolished by pretreatment of ET-1 antibody. Either BQ485 (ET(A) receptor antagonist) or BQ788 (ET(B) receptor antagonist) partially reversed the effect of nicotine, and combination of BQ788 and BQ485 completely abolished it. These findings suggest that nicotine inhibits hepatic circulation through ET-1, and ET(A) and ET(B) receptor.

Physiologically based pharmacokinetic modeling of arterial - antecubital vein concentration difference

BACKGROUND

Modeling of pharmacokinetic parameters and pharmacodynamic actions requires knowledge of the arterial blood concentration. In most cases, experimental measurements are only available for a peripheral vein (usually antecubital) whose concentration may differ significantly from both arterial and central vein concentration.

METHODS

A physiologically based pharmacokinetic (PBPK) model for the tissues drained by the antecubital vein (referred to as "arm") is developed. It is assumed that the "arm" is composed of tissues with identical properties (partition coefficient, blood flow/gm) as the whole body tissues plus a new "tissue" representing skin arteriovenous shunts. The antecubital vein concentration depends on the following parameters: the fraction of "arm" blood flow contributed by muscle, skin, adipose, connective tissue and arteriovenous shunts, and the flow per gram of the arteriovenous shunt. The value of these parameters was investigated using simultaneous experimental measurements of arterial and antecubital concentrations for eight solutes: ethanol, thiopental, 99Tcm-diethylene triamine pentaacetate (DTPA), ketamine, D2O, acetone, methylene chloride and toluene. A new procedure is described that can be used to determine the arterial concentration for an arbitrary solute by deconvolution of the antecubital concentration. These procedures are implemented in PKQuest, a general PBPK program that is freely distributed http://www.pkquest.com.

RESULTS

One set of "standard arm" parameters provides an adequate description of the arterial/antecubital vein concentration for ethanol, DTPA, thiopental and ketamine. A significantly different set of "arm" parameters was required to describe the data for D2O, acetone, methylene chloride and toluene - probably because the "arm" is in a different physiological state.

CONCLUSIONS

Using the set of "standard arm" parameters, the antecubital vein concentration can be used to determine the whole body PBPK model parameters for an arbitrary solute without any additional adjustable parameters. Also, the antecubital vein concentration can be used to estimate the arterial concentration for an arbitrary input for solutes for which no arterial concentration data is available.

Physiologically based pharmacokinetic modeling of arterial - antecubital vein concentration difference

BACKGROUND

Modeling of pharmacokinetic parameters and pharmacodynamic actions requires knowledge of the arterial blood concentration. In most cases, experimental measurements are only available for a peripheral vein (usually antecubital) whose concentration may differ significantly from both arterial and central vein concentration.

METHODS

A physiologically based pharmacokinetic (PBPK) model for the tissues drained by the antecubital vein (referred to as "arm") is developed. It is assumed that the "arm" is composed of tissues with identical properties (partition coefficient, blood flow/gm) as the whole body tissues plus a new "tissue" representing skin arteriovenous shunts. The antecubital vein concentration depends on the following parameters: the fraction of "arm" blood flow contributed by muscle, skin, adipose, connective tissue and arteriovenous shunts, and the flow per gram of the arteriovenous shunt. The value of these parameters was investigated using simultaneous experimental measurements of arterial and antecubital concentrations for eight solutes: ethanol, thiopental, 99Tcm-diethylene triamine pentaacetate (DTPA), ketamine, D2O, acetone, methylene chloride and toluene. A new procedure is described that can be used to determine the arterial concentration for an arbitrary solute by deconvolution of the antecubital concentration. These procedures are implemented in PKQuest, a general PBPK program that is freely distributed http://www.pkquest.com.

RESULTS

One set of "standard arm" parameters provides an adequate description of the arterial/antecubital vein concentration for ethanol, DTPA, thiopental and ketamine. A significantly different set of "arm" parameters was required to describe the data for D2O, acetone, methylene chloride and toluene - probably because the "arm" is in a different physiological state.

CONCLUSIONS

Using the set of "standard arm" parameters, the antecubital vein concentration can be used to determine the whole body PBPK model parameters for an arbitrary solute without any additional adjustable parameters. Also, the antecubital vein concentration can be used to estimate the arterial concentration for an arbitrary input for solutes for which no arterial concentration data is available.