A compartmental model of hepatic disposition kinetics: 1. Model development and application to linear kinetics

The conventional convection-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. The extension of this model to include nonlinear kinetics and zonal heterogeneity of the liver is not straightforward and requires numerical solution of partial differential equation, which is not available in standard nonlinear regression analysis software. In this paper, we describe an alternative compartmental model representation of hepatic disposition (including elimination). The model allows the use of standard software for data analysis and accurately describes the outflow concentration-time profile for a vascular marker after bolus injection into the liver. In an evaluation of a number of different compartmental models, the most accurate model required eight vascular compartments, two of them with back mixing. In addition, the model includes two adjacent secondary vascular compartments to describe the tail section of the concentration-time profile for a reference marker. The model has the added flexibility of being easy to modify to model various enzyme distributions and nonlinear elimination. Model predictions of F, MTT, CV2, and concentration-time profile as well as parameter estimates for experimental data of an eliminated solute (palmitate) are comparable to those for the extended convection-dispersion model.

Commentary: Using the convection-dispersion model and transit time density functions in the analysis of organ distribution kinetics

The convection-dispersion model and its extended form have been used to describe solute disposition in organs and to predict hepatic availabilities. A range of empirical transit-time density functions has also been used for a similar purpose. The use of the dispersion model with mixed boundary conditions and transit-time density functions has been queried recently by Hisaka and Sugiyanaa in this journal. We suggest that, consistent with soil science and chemical engineering literature, the mixed boundary conditions are appropriate providing concentrations are defined in terms of flux to ensure continuity at the boundaries and mass balance. It is suggested that the use of the inverse Gaussian or other functions as empirical transit-time densities is independent of any boundary condition consideration. The mixed boundary condition solutions of the convection-dispersion model are the easiest to use when linear kinetics applies. In contrast, the closed conditions are easier to apply in a numerical analysis of nonlinear disposition of solutes in organs. We therefore argue that the use of hepatic elimination models should be based on pragmatic considerations, giving emphasis to using the simplest or easiest solution that will give a sufficiently accurate prediction of hepatic pharmacokinetics for a particular application. (C) 2000 Wiley-Liss Inc. and the American Pharmaceutical Association J Pharm Sci 89:1579-1586, 2000.

Modulational Instability in Periodic Quadratic Nonlinear Materials

We investigate the modulational instability of plane waves in quadratic nonlinear materials with linear and nonlinear quasi-phase-matching gratings. Exact Floquet calculations, confirmed by numerical simulations, show that the periodicity can drastically alter the gain spectrum but never completely removes the instability. The low-frequency part of the gain spectrum is accurately predicted by an averaged theory and disappears for certain gratings. The high-frequency part is related to the inherent gain of the homogeneous non-phase-matched material and is a consistent spectral feature.

On the nonlinear Neumann problem with critical and supercritical nonlinearities

We investigate the solvability of the Neumann problem (1.1) involving a critical Sobolev exponent. In the first part of this work it is assumed that the coeffcients Q and h are at least continuous. Moreover Q is positive on overline Omega and lambda > 0 is a parameter. We examine the common effect of the mean curvature and the shape of the graphs of the coeffcients Q and h on the existence of low energy solutions. In the second part of this work we consider the same problem with Q replaced by - Q. In this case the problem can be supercritical and the existence results depend on integrability conditions on Q and h.

Nodal solutions for nonlinear eigenvalue problems

We are concerned with determining values of, for which there exist nodal solutions of the boundary value problems u" + ra(t) f(u) = 0, 0 < t < 1, u(O) = u(1) = 0. The proof of our main result is based upon bifurcation techniques.

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.

An improved novel method for solving nonlinear boundary value problems

A modified formula for the integral transform of a nonlinear function is proposed for a class of nonlinear boundary value problems. The technique presented in this paper results in analytical solutions. Iterations and initial guess, which are needed in other techniques, are not required in this novel technique. The analytical solutions are found to agree surprisingly well with the numerically exact solutions for two examples of power law reaction and Langmuir-Hinshelwood reaction in a catalyst pellet.

The kinetics of NO and N2O reduction over coal chars in fluidised-bed combustion

This paper presents a comprehensive and critical review of the mechanisms and kinetics of NO and N2O reduction reaction with coal chars under fluidised-bed combustion conditions (FBC). The heterogeneous reactions of NO and N2O with char/carbon surface have been well recognised as the most important processes in reducing both NOx and N2O in situ FBC. Compared to NO-carbon reactions in FBC, the reactions of N2O with chars have been relatively less understood and studied. Beginning with the overall reaction schemes for both NO and N2O reduction, the paper extensively discusses the reaction mechanisms including the effects of active surface sites. Generally, NO- and N2O-carbon reactions follow a series of step reactions. However, questions remain concerning the role of adsorbed phases of NO and N2O, and the behaviour of different surface sites. Important kinetics factors such as the rate expressions, kinetics parameters as well as the effects of surface area and pore structure are discussed in detail. The main factors influencing the reduction of NO and N2O in FBC conditions are the chemical and physical properties of chars, and the operating parameters of FBC such as temperature, presence of CO, O-2 and pressure. It is shown that under similar conditions, N2O is more readily reduced on the char surface than NO. Temperature was found to be a very important parameter in both NO and N2O reduction. It is generally agreed that both NO- and N2O-carbon reactions follow first-order reaction kinetics with respect to the NO and N2O concentrations. The kinetic parameters for NO and N2O reduction largely depend on the pore structure of chars. The correlation between the char surface area and the reactivities of NO/N2O-char reactions is considered to be of great importance to the determination of the reaction kinetics. The rate of NO reduction by chars is strongly enhanced by the presence of CO and O-2, but these species may not have significant effects on the rate of N2O reduction. However, the presence of these gases in FBC presents difficulties in the study of kinetics since CO cannot be easily eliminated from the carbon surface. In N2O reduction reactions, ash in chars is found to have significant catalytic effects, which must be accounted for in the kinetic models and data evaluation. (C) 1997 Elsevier Science Ltd.

Hepatic structure-pharmacokinetic relationships: The hepatic disposition and metabolite kinetics of a homologous series of O-acyl derivatives of salicylic acid

1 The hepatic disposition and metabolite kinetics of a homologous series of O-acyl (acetyl, propionyl, butanoyl, pentanoyl, hexanoyl and octanoyl) esters of salicylic acid (C2SA, C3SA, C4SA, C5SA, C6SA and C8SA, respectively) was determined using a single-pass, in-sills rat liver preparation. 2 The hepatic venous outflow profiles for the parent esters and the generated metabolite, salicylic acid (SA) were analysed by HPLC. Non-parametric moments analysis was used to determine the area under the curve (AUC'), mean transit time (MTT) and normalized variance (CV2) for the parent esters and generated SA. 3 Pregenerated SA ([C-14]-salicylic acid) was injected into each liver with the parent ester to determine its distribution characteristics. 4 The overall recovery of ester plus metabolite was 89% of the ester dose injected and independent of the ester carbon number, suggesting that ester extraction was due to hepatic metabolism to salicylic acid. 5 The metabolite AUC' value increased directly with the lipophilicity of the parent ester (from 0.12 for C2SA to 0.95 for C8SA). By contrast, the parent AUC' decreased with the lipophilicity (from 0.85 for C2SA to zero for C8SA). The metabolite MTT value also showed a trend to increase with the lipophilicity of the parent ester (from 15.72 s for C3SA to 61.97 s for C8SA). However, the parent MTT value shows no significant change across the series. 6 The two-compartment dispersion model was used to derive the kinetic parameters for parent ester, pregenerated SA and generated SA. Consequently, these parameters were used to estimate the values of AUG', MITT and CV2 for the parent ester and metabolite. The moments values obtained using the two-compartment dispersion model show similar trends to the corresponding moments values obtained from the outflow profiles using a non-parametric approach. 7 The more lipophilic aspirin analogues are more confined to the portal circulation after oral administration than aspirin due to their more extensive hepatic elimination avoiding systemic prostacyclin inhibition. Given that aspirin's selectivity as an anti-thrombotic agent has been postulated to be due to selective anti-platelet effects in the portal circulation, the more lipophilic and highly extracted analogues are potentially more selective anti-thrombotic agents than aspirin.