Search PubMed⌕ Search

Biomedical subjects

L J Aarons

Publications and source records attributed to L J Aarons.

At least 19 recordsLinked to original sources

Optimal sampling strategies for early pharmacodynamic measures in tuberculosis.

OBJECTIVES: To evaluate whether methodological optimization of serial sputum colony counting (SSCC) studies, a potentially important component in the drug development process for tuberculosis, could significantly improve their power. METHODS: Simulations were carried out using a model derived from a large SSCC dataset. Variance inflation factors (VIFs) were calculated for model parameters, focusing on the elimination rate constant likely to reflect 'sterilizing' activity and sampling schemes were optimized relative to a scheme of daily sampling during the initial phase of therapy. Corresponding sample sizes required for SSCC studies using different schemes were also computed. RESULTS: Published sampling schemes lacked efficiency with respect to the 'sterilizing' phase. Pragmatic optimized schemes yielding greatest precision were achieved using eleven sampling points around a skeleton of 0, 2, 7, 14 and 56 days. The standard error of the 'sterilizing' rate constant was reduced more than 4-fold, and sample size for realistic treatment effects was effectively halved. Even schemes with a restricted duration of sampling to avoid high proportions of missing data and those with fewer sampling points still achieved significant gains in precision. Sensitivity analysis suggested that such schemes should continue to perform well over the immediately foreseeable range of improvements in therapy. CONCLUSIONS: Methodological improvements in the design of SSCC studies could make them a powerful tool in Phase II development of anti-tuberculosis agents.

Antitubercular Agents↗

Lumping of whole-body physiologically based pharmacokinetic models.

Lumping is a common pragmatic approach aimed at the reduction of whole-body physiologically based pharmacokinetic (PBPK) model dimensionality and complexity. Incorrect lumping is equivalent to model misspecification with all the negative consequences to the subsequent model implementation. Proper lumping should guarantee that no useful information about the kinetics of the underlying processes is lost. To enforce this guarantee, formal standard lumping procedures and techniques need to be defined and implemented. This study examines the lumping process from a system theory point of view, which provides a formal basis for the derivation of principles and standard procedures of lumping. The lumping principle in PBPK modeling is defined as follows: Only tissues with identical model specification, and occupying identical positions in the system structure should be lumped together at each lumping iteration. In order to lump together parallel tissues, they should have similar or close time constants. In order to lump together serial tissues, they should equilibrate very rapidly with one another. The lumping procedure should include the following stages: (i) tissue specification conversion (when tissues with different model specifications are to be lumped together); (ii) classification of the tissues into classes with significantly different kinetics, according to the basic principle of lumping above; (iii) calculation of the parameters of the lumped compartments; (iv) simulation of the lumped system; (v) lumping of the experimental data; and (vi) verification of the lumped model. The use of the lumping principles and procedures to be adopted is illustrated with an example of a commonly implemented whole-body physiologically based pharmacokinetic model structure to characterize the pharmacokinetics of a homologous series of barbiturates in the rat.

Algorithms↗

Physiologically based pharmacokinetic modeling of a homologous series of barbiturates in the rat: a sensitivity analysis.

Sensitivity analysis studies the effects of the inherent variability and uncertainty in model parameters on the model outputs and may be a useful tool at all stages of the pharmacokinetic modeling process. The present study examined the sensitivity of a whole-body physiologically based pharmacokinetic (PBPK) model for the distribution kinetics of nine 5-n-alkyl-5-ethyl barbituric acids in arterial blood and 14 tissues (lung, liver, kidney, stomach, pancreas, spleen, gut, muscle, adipose, skin, bone, heart, brain, testes) after i.v. bolus administration to rats. The aims were to obtain new insights into the model used, to rank the model parameters involved according to their impact on the model outputs and to study the changes in the sensitivity induced by the increase in the lipophilicity of the homologues on ascending the series. Two approaches for sensitivity analysis have been implemented. The first, based on the Matrix Perturbation Theory, uses a sensitivity index defined as the normalized sensitivity of the 2-norm of the model compartmental matrix to perturbations in its entries. The second approach uses the traditional definition of the normalized sensitivity function as the relative change in a model state (a tissue concentration) corresponding to a relative change in a model parameter. Autosensitivity has been defined as sensitivity of a state to any of its parameters; cross-sensitivity as the sensitivity of a state to any other states' parameters. Using the two approaches, the sensitivity of representative tissue concentrations (lung, liver, kidney, stomach, gut, adipose, heart, and brain) to the following model parameters: tissue-to-unbound plasma partition coefficients, tissue blood flows, unbound renal and intrinsic hepatic clearance, permeability surface area product of the brain, have been analyzed. Both the tissues and the parameters were ranked according to their sensitivity and impact. The following general conclusions were drawn: (i) the overall sensitivity of the system to all parameters involved is small due to the weak connectivity of the system structure; (ii) the time course of both the auto- and cross-sensitivity functions for all tissues depends on the dynamics of the tissues themselves, e.g., the higher the perfusion of a tissue, the higher are both its cross-sensitivity to other tissues' parameters and the cross-sensitivities of other tissues to its parameters; and (iii) with a few exceptions, there is not a marked influence of the lipophilicity of the homologues on either the pattern or the values of the sensitivity functions. The estimates of the sensitivity and the subsequent tissue and parameter rankings may be extended to other drugs, sharing the same common structure of the whole body PBPK model, and having similar model parameters. Results show also that the computationally simple Matrix Perturbation Analysis should be used only when an initial idea about the sensitivity of a system is required. If comprehensive information regarding the sensitivity is needed, the numerically expensive Direct Sensitivity Analysis should be used.

Animals↗

Quantitative structure-pharmacokinetics relationships: I. Development of a whole-body physiologically based model to characterize changes in pharmacokinetics across a homologous series of barbiturates in the rat.

As part of an overall program to develop a framework for evaluating the contribution of structural and physicochemical properties to pharmacokinetics, the distribution kinetics of nine 5-n-alkyl-5-ethyl barbituric acids in arterial blood and 14 tissues (lung, liver, kidney, stomach, pancreas, spleen, gut, muscle, adipose, skin, bone, heart, brain, testes) was examined after i.v. bolus administration in rats. The barbituric acids studied form a true homologous series; therefore any differences in pharmacokinetics, noted between congeners, can be directly linked to the increase in lipophilicity, resulting from the addition of a methylene group. A whole-body physiologically based pharmacokinetic model has been developed, assuming most of the tissues to be well-stirred compartments. Brain and testes, in which distribution for the lower homologues was permeability rate-limited, were represented by two compartments. For each homologue, the model parameters have been optimized, using the tissue concentration-time data. The initial distribution processes in the system were very rapid, making it quite stiff, and essentially over before the first samples were taken. A progressively increasing redistribution from lean tissues into adipose on ascending the homologous series was observed, characterized by a tendency for a progressive decrease in the magnitude of the concentration-time profiles for some of the lean and well-perfused tissues, an increase in the adipose concentration-time profile, and an increase in the time to reach the maximum adipose concentration. A shift from permeability rate limitation to perfusion rate limitation of the distribution processes for brain and testes, as well as an increase in the intrinsic hepatic clearance and decrease in the renal clearance with the increase of lipophilicity of the homologues, were quantified. An increase in the total unbound volume of distribution on ascending the homologous series was also observed. Muscle was found to be the major drug depot at steady state, accounting for approximately 50% of the total unbound volume of distribution, regardless of the lipophilicity of the homologue; the unbound volume of distribution of adipose increases more than 10-fold with the increase of lipophilicity.

Adipose Tissue↗

Application of a linear recirculation model to drug targeting.

Current interest in drug targeting has inspired theoretical considerations of its potential and problems. Previously, drug targeting has been considered in terms of more or less elaborate compartmental models. The present paper shows how an equivalent analysis of the potential advantage of drug targeting may be derived with the minimum reliance on a specific model. A linear recirculation model is used to describe the drug concentration profile at some target site and in the rest of the body. Equations for the AUCs of drug and of a drug-carrier conjugate can then be derived. These AUCs are used to define a drug targeting index (DTI), a measure of drug targeting selectivity previously derived from a specific model. It is shown that the DTI can be defined solely in terms of extraction ratios for elimination of free drug, when release of drug is confined to the target site. The expression for DTI is shown to be equivalent to that previously derived from several more model-dependent approaches.

Drug Administration Routes↗

Mean residence time for drugs subject to enterohepatic cycling.

A physiologically realistic model of enterohepatic cycling (EHC) which includes separate liver and gallbladder compartments, discontinuous gallbladder emptying and first-order absorption from both an oral formulation and secreted bile (kapo and kab, respectively) has been developed. The effect of EHC on area under the first-moment curve (AUMC) of drug concentration in plasma and on parameters derived from the AUMC was investigated. Unlike AUC, AUMC is dependent on the time and time-course of gallbladder emptying, increasing as the interval between gallbladder emptying increases. Consequently, mean residence time (MRT) is also a time-dependent parameter. Analytical solutions for MRTiv and MRTpo were derived. Mean absorption time (MAT = MRTpo - MRTiv) is also time-dependent, contrary to findings previously published for a model of EHC with a continuous time lag. MAT is also dependent on kapo, kba and the hepatic extraction ratio. The difference between MRTpos for two formulations with unequal kapo values may deviate from the difference in the inverse of their absorption rate constants. Implications for design and interpretation of pharmacokinetic studies include (i) MAT values may be dominated by the time-course of recycling rather than the time-course of the initial absorption, depending on the extent of EHC and (ii) the unpredictable nature of the time of gallbladder emptying will contribute to intrasubject variability in derived parameters during crossover studies. Knowledge of the extent of EHC is invaluable in deciding whether modification of the in vitro release characteristics of an oral formulation will have any effect on the overall time-course of absorption in vivo. Techniques to monitor or control gallbladder emptying may be helpful for reducing variability in pharmacokinetic studies for compounds which are extensively cycled in bile.

Absorption↗

Interaction of mixed micelles formed from glycocholic acid and lecithin with the protein binding of various drugs.

Mixed micelles (MM) formed from glycocholic acid and lecithin are suited to solubilize lipophilic drugs for intravenous use. To test for possible drug-drug interactions, the protein binding of a series of agents known to bind to different sites on albumin (diazepam, warfarin, ketoprofen, frusemide, probenecid) and additionally (prazosin, quinidine, propranolol) or exclusively (disopyramide) to alpha 1-acid glycoprotein or to transcortin (prednisolone) was determined in the presence and absence of MM. Concentrations of MM, corresponding to the maximum possible plasma concentration achieved by injecting the highest clinical doses of MM into the systemic circulation, had little or no effect on the unbound fractions of drugs known to bind exclusively to albumin. Only at five times higher MM concentrations were the free fractions substantially increased (by up to 45%). Unbound fractions of drugs bound with high affinity but low capacity to alpha 1-acid glycoprotein were increased between 50-85% even at 'therapeutic' doses of MM. The present study suggests that drugs solubilized by MM should be given by slow injection or infusion to patients already receiving drugs which are highly bound to alpha 1-acid glycoprotein.

Glycocholic Acid↗

Estimation of area under the curve for drugs subject to enterohepatic cycling.

A physiologically realistic model is used to provide insight into the design of sampling protocols for accurate determination of AUC(0-infinity) for drugs subject to enterohepatic cycling. Through simulation of plasma concentration-time curves for such drugs it is found that more than one peak is predicted after oral and intravenous administration of a single dose of drug, the relative magnitude of peaks is dependent on the hepatic extraction ratio for both oral and intravenous drug administration, the percent of the AUC(0-infinity) in later time intervals is also a function of the hepatic extraction ratio, and present methods for the design of sampling protocols may not provide accurate estimates of AUC(0-infinity) (especially for highly extracted drugs), because peaks are only evident at later times after intravenous administration when plasma sampling is less frequent, much of the area occurs at later times, and the amount of drug in the sampling compartment after oral administration is much lower than that after intravenous administration of drug and could be incorrectly interpreted as low bioavailability if sampling is not carried out for a long period of time. The types of oral and intravenous profiles predicted for highly extracted drugs are exemplified by data for naltrexone in the monkey.

Administration, Oral↗

Pharmacokinetics of lidocaine and its deethylated metabolite: dose and time dependency studies in man.

The concentrations of lidocaine and of its deethylated metabolite, MEGX, were measured in blood following the intravenous administration of 50 and 100 mg lidocaine hydrochloride, the oral administration of 100, 300, and 500 mg lidocaine hydrochloride monohydrate, and the oral administration of 300 mg lidocaine hydrochloride monohydrate every 8 h for seven doses, to three healthy volunteers. The range of values for the parameters defining the disposition kinetics of lidocaine were: terminal half-life, 50-231 min; total clearance, 13-17 ml/min/kg; initial dilution space, 0.13-2.5 liters/kg; and volume of distribution at steady state, 0.6-4.5 liters/kg. Lidocaine absorption from solution was rapid, but due to presystemic hepatic metabolism, the availability was low, the range of average values lying between 0.19 and 0.38. No dose or time dependency in lidocaine and monoethylglycinexylidide pharmacokinetics following the single dose studies of lidocaine were noted. Effective hepatic blood flow, based on total clearance and availability measurements, was estimated to be 18-27 ml/min/kg. The concentrations of MEGX were approximately one-third of those of lidocaine following intravenous lidocaine and were comparable following oral lidocaine, but as predicted, the dose normalized area under the MEGX concentration-time curve was constant and independent of the route of administration of lidocaine. In two subjects, the blood concentrations of lidocaine and MEGX following multiple doses of oral lidocaine were those predicted from the single dose studies. In the third subject, the degree of accumulation of lidocaine was greater than predicted. The reasons and mechanism for this difference between subjects on multiple dosing remains unclear.

Administration, Oral↗

Interpretation of CO2 exhalation rate data from demethylation of aminopyrine and its metabolite monomethylaminoantipyrine.

1 Aminopyrine breath tests make use of the commercially available (N-dimethyl-[14C])-aminopyrine. A pharmacokinetic model has been proposed to relate 14CO2 exhalation rates (CER) to the demethylation of ([14C]-methyl)-aminopyrine (AP) and -monomethylaminoantipyrine (MAP). 2 computer simulations based on the model show that the shape of the CER-time profile is largely dependent on the ratio of the MAP to AP elimination rate constants. If this ratio equals 0.5 then the CERs decline in the monoexponential fashion. Ratios less than 0.5 result in concave biexponential curves whereas ratios greater than 0.5 result in convex curves. When demethylation is not complete for both compounds the transfer from biexponential to monoexponential behaviour will only occur at ratios greater than 0.5. 3 The resolution of concave biexponential CER-time profiles to give accurate estimates of AP and MAP elimination rate constants can only be achieved when the length of the experiment is adequate. The commonly employed 2 microCi tracer dose of aminopyrine is insufficient to monitor CER over the necessary time period to detect the proposed biexponential decline.

Aminopyridines↗

Kinetics of drug displacement interactions.

A simple model simulating the kinetics of drugs displacement kinetics is investigated. It is demonstrated that for highly bound, lowly cleared drugs, displacement interactions are transitory. Consequently, the kinetics of the interaction have to be considered as well as the in vitro interaction. It is possible to have a significant in vitro displacement interaction with no in vivo counterpart. Methods of moderating drug displacement by adjusting the rate and the timing of administration of the displacing agent are discussed.

Binding, Competitive↗

An in vitro study of drug displacement interactions: warfarin-salicylate and warfarin-phenylbutazone.

The binding interactions between warfarin-salicylate and warfarin-phenylbutazone in the presence of 4 g percent bovine serum albumin at 37 degrees C were studied using equilibrium dialysis. Methods of representing and analysing drug binding interactions are discussed. Scatchard plots, double reciprocal plots and the like are shown to be of no use in representing drug displacement interactions since they display only one drug and they can be potentially misleading. It is argued that a preferable method of analysing drug displacement interaction data is in terms of a stepwise multiple equilibria model. The numerical problems associated with fitting this kind of model to the data are discussed. A three-dimensional representation of the binding surface is proposed as a superior means of visualizing drug displacement interactions.

Animals↗

Problems associated with analysis and interpretation of small molecule/macromolecule binding data.

In the analysis of binding data, arbitrary transformations such as the Scatchard plot, may give misleading estimates of the binding parameters. The statistically correct approach is to determine values of K and n by non-linear regression of the actual dependent variable against the actual independent variable. In the case of the spectrophotometric titration method the dependent variable is the absorbance and the independent variable is the composition of the drug/macromolecule mixture. The method relies on an accurate estimate of the extinction coefficient of the bound drug and this is best treated as a parameter to be estimated in the regression analysis. In testing models by data fits alone it is emphasized that whilst a model may be rejected if it does not fit the data, a good fit does not ensure uniquieness and confirmatory, independent evidence must be sought.

DNA↗