Modeling of hepatic elimination and organ distribution kinetics with the extended convection-dispersion model

The conventional convection-dispersion (also called axial dispersion) model is widely used to interrelate hepatic availability (F) and clearance (Cl) with the morphology and physiology of the liver and to predict effects such as changes in liver blood flow on F and Cl. An extended form of the convection-dispersion model has been developed to adequately describe the outflow concentration-time profiles for vascular markers at both short and long times after bolus injections into perfused livers. The model, based on flux concentration and a convolution of catheters and large vessels, assumes that solute elimination in hepatocytes follows either fast distribution into or radial diffusion in hepatocytes. The model includes a secondary vascular compartment, postulated to be interconnecting sinusoids. Analysis of the mean hepatic transit time (MTT) and normalized variance (CV2) of solutes with extraction showed that the discrepancy between the predictions of MTT and CV2 for the extended and conventional models are essentially identical irrespective of the magnitude of rate constants representing permeability, volume, and clearance parameters, providing that there is significant hepatic extraction. In conclusion, the application of a newly developed extended convection-dispersion model has shown that the unweighted conventional convection-dispersion model can be used to describe the disposition of extracted solutes and, in particular, to estimate hepatic availability and clearance in booth experimental and clinical situations.

Modeling of subcutaneous absorption kinetics of infusion solutions in the elderly using technetium

Absorption kinetics of solutes given with the subcutaneous administration of fluids is ill-defined. The gamma emitter, technitium pertechnetate, enabled estimates of absorption rate to be estimated independently using two approaches. In the first approach, the counts remaining at the site were estimated by imaging above the subcutaneous administration site, whereas in the second approach, the plasma technetium concentration-time profiles were monitored up to 8 hr after technetium administration. Boluses of technetium pertechnetate were given both intravenously and subcutaneously on separate occasions with a multiple dosing regimen using three doses on each occasion. The disposition of technetium after iv administration was best described by biexponential kinetics with a V-ss of 0.30 +/- 0.11 L/kg and a clearance of 30.0 +/- 13.1 ml/min. The subcutaneous absorption kinetics was best described as a single exponential process with a half-life of 18.16 +/- 3.97 min by image analysis and a half-life of 11.58 +/- 2.48 min using plasma technetium time data. The bioavailability of technetium by the subcutaneous route was estimated to be 0.96 +/- 0.12. The absorption half-life showed no consistent change with the duration of the subcutaneous infusion. The amount remaining at the absorption site with time was similar when analyzed using image analysis, and plasma concentrations assuming multiexponential disposition kinetics and a first-order absorption process. Profiles of fraction remaining at the absorption sire generated by deconvolution analysis, image analysis, and assumption of a constant first-order absorption process were similar. Slowing of absorption from the subcutaneous administration site is apparent after the last bolus dose in three of the subjects and can De associated with the stopping of the infusion. In a fourth subject, the retention of technetium at the subcutaneous site is more consistent with accumulation of technetium near the absorption site as a result of systemic recirculation.

Diffusion modeling of percutaneous absorption kinetics: 2. Finite vehicle volume and solvent deposited solids

The diffusion model for percutaneous absorption is developed for the specific case of delivery to the skin being limited by the application of a finite amount of solute. Two cases are considered; in the first, there is an application of a finite donor (vehicle) volume, and in the second, there are solvent-deposited solids and a thin vehicle with a high partition coefficient. In both cases, the potential effect of an interfacial resistance at the stratum corneum surface is also considered. As in the previous paper, which was concerned with the application of a constant donor concentration, clearance limitations due to the viable eqidermis, the in vitro sampling rate, or perfusion rate in vivo are included. Numerical inversion of the Laplace domain solutions was used for simulations of solute flux and cumulative amount absorbed and to model specific examples of percutaneous absorption of solvent-deposited solids. It was concluded that numerical inversions of the Laplace domain solutions for a diffusion model of the percutaneous absorption, using standard scientific software (such as SCIENTIST, MicroMath Scientific software) on modern personal computers, is a practical alternative to computation of infinite series solutions. Limits of the Laplace domain solutions were used to define the moments of the flux-time profiles for finite donor volumes and the slope of the terminal log flux-time profile. The mean transit time could be related to the diffusion time through stratum corneum, viable epidermal, and donor diffusion layer resistances and clearance from the receptor phase. Approximate expressions for the time to reach maximum flux (peak time) and maximum flux were also derived. The model was then validated using reported amount-time and flux-time profiles for finite doses applied to the skin. It was concluded that for very small donor phase volume or for very large stratum corneum-vehicle partitioning coefficients (e.g., for solvent deposited solids), the flux and amount of solute absorbed are affected by receptor conditions to a lesser extent than is obvious for a constant donor constant donor concentrations. (C) 2001 Wiley-Liss, Inc. and the American Pharmaceutical Association J Pharm Sci 90:504-520, 2001.

Diffusion modeling of percutaneous absorption kinetics: 3. Variable diffusion and partition coefficients, consequences for stratum corneum depth profiles and desorption kinetics

Stratum corneum (SC) desorption experiments have yielded higher calculated steady-state fluxes than those obtained by epidermal penetration studies. A possible explanation of this result is a variable diffusion or partition coefficient across the SC. We therefore developed the diffusion model for percutaneous penetration and desorption to study the effects of either a variable diffusion coefficient or variable partition coefficient in the SC over the diffusion path length. Steady-state flux, lag time, and mean desorption time were obtained from Laplace domain solutions. Numerical inversion of the Laplace domain solutions was used for simulations of solute concentration-distance and amount penetrated (desorbed)-time profiles. Diffusion and partition coefficients heterogeneity were examined using six different models. The effect of heterogeneity on predicted flux from desorption studies was compared with that obtained in permeation studies. Partition coefficient heterogeneity had a more profound effect on predicted fluxes than diffusion coefficient heterogeneity. Concentration-distance profiles show even larger dependence on heterogeneity, which is consistent with experimental tape-stripping data reported for clobetasol propionate and other solutes. The clobetasol propionate tape-stripping data were most consistent with the partition coefficient decreasing exponentially for half the SC and then becoming a constant for the remaining SC. (C) 2004 Wiley-Liss, Inc.

Reactivity of chars and carbons: New insights through molecular modeling

A new model proposed for the gasification of chars and carbons incorporates features of the turbostratic nanoscale structure that exists in such materials. The model also considers the effect of initial surface chemistry and different reactivities perpendicular to the edges and to the faces of the underlying crystallite planes comprising the turbostratic structure. It may be more realistic than earlier models based on pore or grain structure idealizations when the carbon contains large amounts of crystallite matter. Shrinkage of the carbon particles in the chemically controlled regime is also possible due to the random complete gasification of crystallitic planes. This mechanism can explain observations in the literature of particle size reduction. Based on the model predictions, both initial surface chemistry and the number of stacked planes in the crystallites strongly influence the reactivity and particle shrinkage. Its test results agree well with literature data on the air-oxidation of Spherocarb and show that it accurately predicts the variation of particle size with conversion. Model parameters are determined entirely from rate measurements.

Modeling of columnar-to-equiaxed transition in solidified Al-Si alloys

A Cellular-Automaton Finite-Volume-Method (CAFVM) algorithm has been developed, coupling with macroscopic model for heat transfer calculation and microscopic models for nucleation and growth. The solution equations have been solved to determine the time-dependent constitutional undercooling and interface retardation during solidification. The constitutional undercooling is then coupled into the CAFVM algorithm to investigate both the effects of thermal and constitutional undercooling on columnar growth and crystal selection in the columnar zone, and formation of equiaxed crystals in the bulk liquid. The model cannot only simulate microstructures of alloys but also investigates nucleation mechanisms and growth kinetics of alloys solidified with various solute concentrations and solidification morphologies.

Race modeling: An effective cognitive strategy for the 100 m sprinter?

The purpose of this study was to investigate the effects of a specific cognitive race plan on 100 m sprint performance, Twelve elite sprinters (11 male and 1 female) performed 100 m time trials under normal (control) conditions and then under experimental conditions (use of race cues). In the experimental condition, participants were asked to think about specific thought content in each of three segments of the 100 m. A multiple baseline design was employed. A mean improvement of 0.26 s was found. Eleven of the 12 participants showed improvement using the specific cognitive race plan (p < .005). Participants also produced more consistent sprint performances when using the cues (p < .01). Subjective evaluations made by the participants unanimously supported the use of the race plan for optimizing sprint performance. Environmental conditions, effort, and practice effects were considered as possible influences on the results.

Surface stickiness of drops of carbohydrate and organic acid solutions during convective drying: Experiments and modeling

Drying kinetics of low molecular weight sugars such as fructose, glucose, sucrose and organic acid such as citric acid and high molecular weight carbohydrate such as maltodextrin (DE 6) were determined experimentally using single drop drying experiments as well as predicted numerically by solving the mass and heat transfer equations. The predicted moisture and temperature histories agreed with the experimental ones within 6% average relative (absolute) error and average difference of +/- 1degreesC, respectively. The stickiness histories of these drops were determined experimentally and predicted numerically based on the glass transition temperature (T-g) of surface layer. The model predicted the experimental observations with good accuracy. A nonsticky regime for these materials during spray drying is proposed by simulating a drop, initially 120 mum in diameter, in a spray drying environment.

A numerical approach to modeling the catalytic voltammetry of surface-confined redox enzymes

A finite difference method for simulating voltammograms of electrochemically driven enzyme catalysis is presented. The method enables any enzyme mechanism to be simulated. The finite difference equations can be represented as a matrix equation containing a nonlinear sparse matrix. This equation has been solved using the software package Mathematica. Our focus is on the use of cyclic voltammetry since this is the most commonly employed electrochemical method used to elucidate mechanisms. The use of cyclic voltammetry to obtain data from systems obeying Michaelis-Menten kinetics is discussed, and we then verify our observations on the Michaelis-Menten system using the finite difference simulation. Finally, we demonstrate how the method can be used to obtain mechanistic information on a real redox enzyme system, the complex bacterial molybdoenzyme xanthine dehydrogenase.

Disposition kinetics of propranolol isomers in the perfused rat liver

The aim of this study was to define the determinants of the linear hepatic disposition kinetics of propranolol optical isomers using a perfused rat liver. Monensin was used to abolish the lysosomal proton gradient to allow an estimation of propranolol ion trapping by hepatic acidic vesicles. In vitro studies were used for independent estimates of microsomal binding and intrinsic clearance. Hepatic extraction and mean transit time were determined from outflow-concentration profiles using a nonparametric method. Kinetic parameters were derived from a physiologically based pharmacokinetic model. Modeling showed an approximate 34-fold decrease in ion trapping following monensin treatment. The observed model-derived ion trapping was similar to estimated theoretical values. No differences in ion-trapping values was found between R(+)- and S(-)- propranolol. Hepatic propranolol extraction was sensitive to changes in liver perfusate flow, permeability-surface area product, and intrinsic clearance. Ion trapping, microsomal and nonspecific binding, and distribution of unbound propranolol accounted for 47.4, 47.1, and 5.5% of the sequestration of propranolol in the liver, respectively. It is concluded that the physiologically more active S(-)- propranolol differs from the R(+)- isomer in higher permeability-surface area product, intrinsic clearance, and intracellular binding site values.

Enzyme deactivation in fixed bed reactors with michaelis-menten kinetics

Rate expression for enzyme poisoning which are consistent with a Michaelis-Menten main reaction are used to analyze the performance of a fixed bed reactor containing immobilized enzyme. When enzyme deactivation results from the irreversible bonding of a product molecule to an existing substrate-enzyme complex, it is shown that minimum enzyme activity can occur in the interior of the bed, well away from the ends. This suggests that bed sectioning techniques may enable direct evaluation of fundamental poisoning mechanisms.

Fixed bed reactors with catalyst poisoning - first order kinetics.

The long performance of an isothermal fixed bed reactor undergoing catalyst poisoning is theoretically analyzed using the dispersion model. First order reaction with dth order deactivation is assumed and the model equations are solved by matched asymptotic expansions for large Peclet number. Simple closed-form solutions, uniformly valid in time, are obtained.

A finite integral transform technique for solving the diffusion-reaction equation with Michaelis-Menten kinetics

An approximate analytical technique employing a finite integral transform is developed to solve the reaction diffusion problem with Michaelis-Menten kinetics in a solid of general shape. A simple infinite series solution for the substrate concentration is obtained as a function of the Thiele modulus, modified Sherwood number, and Michaelis constant. An iteration scheme is developed to bring the approximate solution closer to the exact solution. Comparison with the known exact solutions for slab geometry (quadrature) and numerically exact solutions for spherical geometry (orthogonal collocation) shows excellent agreement for all values of the Thiele modulus and Michaelis constant.

Substrate-inhibited kinetics with catalyst deactivation in an isothermal CSTR

Analytical expressions are developed for the time-dependent reactant concentration and catalyst activity in an isothermal CSTR with Langmuir-Hinshelwood kinetics of deactivation and reaction. Several parallel and series posioning mechanisms are considered for a reactor which, without poisoning, would operate at a unique steady state. The use of matched asymptotic expansions and abandonment of the usual initial-steady-state assumption give results, valid from startup to final loss of activity, whose accuracy can be improved systematically.

Substrate-inhibited kinetics with catalyst deactivation in an isothermal CSTR. II. Multiple pseudosteady states and reactor failure

Analytical expressions are derived for the time and magnitude of failure of an isothermal CSTR with substrate-inhibited kinetics, caused by slow catalyst deactivation under three types of parallel and series mechanisms. Reactors operating at high space velocity are found to be most susceptible to early failure and poisoning by product is more dangerous than by reactant. The magnitude of the jump across steady states depends solely on the Langmuir-Hinshelwood kinetic parameters and a detailed analysis of reactor behavior during the jump itself is given.