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.

Analysis of multicomponent adsorption kinetics on activated carbon

An integrated mathematical model for the kinetics of multicomponent adsorption on microporous carbon was developed. Transport in this bidisperse solid is represented by balance equations in the macropore and micropore phases, in which gas-phase diffusion dominates the mass transfer in the macropores, with the phenomenological diffusivities represented by the generalized Maxwell-Stefan (GMS) formulation. Viscous flow also contributes to the macropore fluxes and is included in the MS expressions. Diffusion of the adsorbed phase controls the mass transfer in the micro ore phase, p which is also described in a similar way by the MS method. The adsorption isotherms are represented by a new heterogeneous modified vacancy solution theory formulation of adsorption, which has proved to be a robust method for adsorption on activated carbons. The model is applied to the coadsorption and codesorption of C2H6 and C3H8 on Ajax and Norit carbon, as well as the displacement on Ajax carbon. The effect of the viscous flow in the macropore phase is not significant for the cases studied. The model accurately predicts the overshoot behavior and rollup of C2H6 during coadsorption. The prediction for the heavier compound C3H8 is always satisfactory, though at higher C3H8 mole fraction, the overshoot extent of C2H6 is overpredicted, possibly due to neglect of heat effects.

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.

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.

Hepatic disposition and metabolite kinetics of a homologous series of diflunisal esters

The hepatic disposition and metabolite kinetics of a homologous series of diflunisal O-acyl esters (acetyl, butanoyl, pentanoyl, anti hexanoyl) were determined using a single-pass perfused in situ rat liver preparation. The experiments were conducted using 2% BSA Krebs-Henseleit buffer (pH 7.4), and perfusions were performed at 30 mL/min in each liver. O-Acyl esters of diflunisal and pregenerated diflunisal were injected separately into the portal vein. The venous outflow samples containing the esters and metabolite diflunisal were analyzed by high performance liquid chromatography (HPLC). The normalized outflow concentration-time profiles for each parent ester and the formed metabolite, diflunisal, were analyzed using statistical moments analysis and the two-compartment dispersion model. Data (presented as mean +/- standard error for triplicate experiments) was compared using ANOVA repeated measures, significance level P < 0.05. The hepatic availability (AUC'), the fraction of the injected dose recovered in the outflowing perfusate, for O-acetyldiflunisal (C2D = 0.21 +/- 0.03) was significantly lower than the other esters (0.34-0.38). However, R-N/f(u), the removal efficiency number R-N divided by the unbound fraction in perfusate f(u), which represents the removal efficiency of unbound ester by the liver, was significantly higher for the most lipophilic ester (O-hexanoyldiflunisal, C6D = 16.50 +/- 0.22) compared to the other members of the series (9.57 to 11.17). The most lipophilic ester, C6D, had the largest permeability surface area (PS) product (94.52 +/- 38.20 mt min-l g-l liver) and tissue distribution value VT (35.62 +/- 11.33 mL g(-1) liver) in this series. The MTT of these O-acyl esters of diflunisal were not significantly different from one another. However, the metabolite diflunisal MTTs tended to increase with the increase in the parent ester lipophilicity (11.41 +/- 2.19 s for C2D to 38.63 +/- 9.81 s for C6D). The two-compartment dispersion model equations adequately described the outflow profiles for the parent esters and the metabolite diflunisal formed from the O-acyl esters of diflunisal in the liver.

Kinetic analysis of vascular marker distribution in perfused rat livers after regeneration following partial hepatectomy

Background/Aims: Liver clearance models are based on information (or assumptions) on solute distribution kinetics within the microvasculatory system, The aim was to study albumin distribution kinetics in regenerated livers and in livers of normal adult rats, Methods: A novel mathematical model was used to evaluate the distribution space and the transit time dispersion of albumin in livers following regeneration after a two-thirds hepatectomy compared to livers of normal adult rats. Outflow curves of albumin measured after bolus injection in single-pass perfused rat livers were analyzed by correcting for the influence of catheters and fitting a long-tailed function to the data. Results: The curves were well described by the proposed model. The distribution volume and the transit time dispersion of albumin observed in the partial hepatectomy group were not significantly different from livers of normal adult rats. Conclusions: These findings suggest that the distribution space and the transit time dispersion of albumin (CV2) is relatively constant irrespective of the presence of rapid and extensive repair. This invariance of CV2 implies, as a first approximation, a similar degree of intrasinusoidal mixing, The finding that a sum of two (instead of one) inverse Gaussian densities is an appropriate empirical function to describe the outflow curve of vascular indicators has consequences for an improved prediction of hepatic solute extraction.

Thermogravimetric analysis of carbon deposition over Ni/gamma-Al2O3 catalysts in carbon dioxide reforming of methane

Carbon formation on Ni/gamma-Al2O3 catalysts and its kinetics during methane reforming with carbon dioxide was studied in the temperature range of 500-700 degrees C using a thermogravimetric analysis technique. The activation energies of methane cracking, carbon gasification in CO2, as well as carbon deposition in CO2-CH4 reforming were obtained. The results show that the activation energy for carbon gasification is larger than that of carbon formation in methane cracking and that the activation energy of coking in CO2-CH4 reforming is also larger than that of methane decomposition to carbon. The dependencies of coking rate on partial pressures of CH4 and CO2 indicate that methane decomposition is the main route for carbon deposition. A mechanism and kinetic model for carbon deposition is proposed.

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.

Functional analysis, using in vitro mutagenesis, of amino acids located in the phenylalanine hydroxylase active site

The 3-dimensionaI structure determination of rat phenylalanine hydroxylase (PAH) has identified potentially important amino acids lining the active site cleft with the majority of these having hydrophobic side-chains including several with aromatic side chains. Here we have analyzed the effect on rat PAH enzyme kinetics of in vitro mutagenesis of a number of these amino acids lining the PAH active site. Mutation of F299, Y324, F331, and Y343 caused a significant decrease in enzyme activity but no change in the K-m for substrate or cofactor. me conclude that these aromatic residues are essential for activity but are not significantly involved in binding of the substrate or cofactor. in contrast the PAH mutant, S349T, showed an 18-fold increase in K-m for phenylalanine, showing the first functional evidence that this residue was binding at or near the phenylalanine binding site. This confirms the recently published model for the binding of phenylalanine to the PAH active site that postulated S349 interacts with the amino group on the main chain of the phenylalanine molecule. This result differs with that found for the equivalent mutation (S395T), in the closely related tyrosine hydroxylase, which had no effect on substrate K-m, showing that while the architecture of the two active sites are very similar the amino acids that bind to the respective substrates are different. (C) 2000 Academic Press.

Metabolic and kinetic analysis of poly(3-hydroxybutyrate) production by recombinant Escherichia coli

A quantitatively repeatable protocol was developed for poly(3-hydroxybutyrate) (PHB) production by Escherichia coli XL1-Blue (pSYL107). Two constant-glucose fed-batch fermentations of duration 25 h were carried out in a 5-L bioreactor, with the measured oxygen volumetric mass-transfer coefficient (k(L)a) held constant at 1.1 min(-1). All major consumption and production rates were quantified. The intracellular concentration profiles of acetyl-CoA (300 to 600 mug.g RCM-1) and 3-hydroxy-butyryl-CoA (20 to 40 mug.g RCM-1) were measured, which is the first time this has been performed for E. coli during PHB production. The kinetics of PHB production were examined and likely ranges were established for polyhydroxyalkanoate (PHA) enzyme activity and the concentration of pathway metabolites. These measured and estimated values are quite similar to the available literature estimates for the native PHB producer Ralstonia eutropha. Metabolic control analysis performed on the PHB metabolic pathway showed that the PHB flux was highly sensitive to acetyl-CoA/CoA ratio (response coefficient 0.8), total acetyl-CoA + CoA concentration (response coefficient 0.7), and pH (response coefficient -1.25). It was less sensitive (response coefficient 0.25) to NADPH/NADP ratio. NADP(H) concentration (NADPH + NADP) had a negligible effect. No single enzyme had a dominant flux control coefficient under the experimental conditions examined (0.6, 0.25, and 0.15 for 3-ketoacyl-CoA reductase, PHA synthase, and 3-ketothiolase, respectively). In conjunction with metabolic flux analysis, kinetic analysis was used to provide a metabolic explanation for the observed fermentation profile. In particular, the rapid onset of PHB production was shown to be caused by oxygen limitation, which initiated a cascade of secondary metabolic events, including cessation of TCA cycle flux and an increase in acetyl-CoA/CoA ratio. (C) 2001 John Wiley & Sons. Inc. Biotechnol Bioeng 74: 70-80, 2001.

A computational analysis of substrate binding strength by phosphorylase kinase and protein kinase A

Protein kinases exhibit various degrees of substrate specificity. The large number of different protein kinases in the eukaryotic proteomes makes it impractical to determine the specificity of each enzyme experimentally. To test if it were possible to discriminate potential substrates from non-substrates by simple computational techniques, we analysed the binding enthalpies of modelled enzyme-substrate complexes and attempted to correlate it with experimental enzyme kinetics measurements. The crystal structures of phosphorylase kinase and cAMP-dependent protein kinase were used to generate models of the enzyme with a series of known peptide substrates and non-substrates, and the approximate enthalpy of binding assessed following energy minimization. We show that the computed enthalpies do not correlate closely with kinetic measurements, but the method can distinguish good substrates from weak substrates and non-substrates. Copyright (C) 2002 John Wiley Sons, Ltd.

Altered kinetics and benzodiazepine sensitivity of a GABA(A) receptor subunit mutation [gamma(2)(R43Q)] found in human epilepsy

The gamma-aminobutyric acid type A (GABA(A)) receptor mediates fast inhibitory synaptic transmission in the CNS. Dysfunction of the GABA(A) receptor would be expected to cause neuronal hyperexcitability, a phenomenon linked with epileptogenesis. We have investigated the functional consequences of an arginine-to-glutamine mutation at position 43 within the GABA(A) gamma(2)-subunit found in a family with childhood absence epilepsy and febrile seizures. Rapid-application experiments performed on receptors expressed in HEK-293 cells demonstrated that the mutation slows GABA(A) receptor deactivation and increases the rate of desensitization, resulting in an accumulation of desensitized receptors during repeated, short applications. In Xenopus laevis oocytes, two-electrode voltage-clamp analysis of steady-state currents obtained from alpha(1)beta(2)gamma(2) or alpha(1)beta(2)gamma(2)(R43Q) receptors did not reveal any differences in GABA sensitivity. However, differences in the benzodiazepine pharmacology of mutant receptors were apparent. Mutant receptors expressed in oocytes displayed reduced sensitivity to diazepam and flunitrazepam but not the imiclazopyricline zolpidem. These results provide evidence of impaired GABA(A) receptor function that could decrease the efficacy of transmission at inhibitory synapses, possibly generating a hyperexcitable neuronal state in thalamocortical networks of epileptic patients possessing the mutant subunit.

Kinetics of photoacclimation in corals

Traditional models describing the relationship between photosynthesis (P) and irradiance (I) do not account for photoacclimation to short-term variation in irradiance. Here we develop and test a model that predicts the rate of photosynthesis under fluctuating irradiances at the scale of days to weeks. Using oxygen respirometry, we measured the rates of change in the P-I model parameters P-max (maximum rate of gross photosynthesis) and I-k (sub-saturation irradiance) of the photo-symbiotic coral Turbinaria mesenterina (Lamarck) following large and small increases and decreases in growth irradiance. We analyse the behaviour of the dynamic P-I model in turbid-water conditions using a dataset of 3-month continuous irradiance as the input variable. In response to upward or downward changes in experimental growth irradiance, I-k values decreased or increased exponentially, reaching new and stable levels within 5-10 days. I-k responded 4 times stronger than P-max to changes in growth irradiance. The kinetics of I-k did not show hysteresis, and changed in similar ways when irradiance was increased or decreased in small or large amounts. This suggests that mechanisms associated with photo-protection during increases in irradiance, and the maximisation of photosynthetic efficiency during decreases in irradiance, are equally potent. On the scale of months, the dynamic P-I model did not predict higher rates of photosynthesis than the static P-I model, but buffered the variation in photosynthesis during periods of reduced irradiance. Fourier analysis indicated that the kinetics of I-k closely matches the main periodicities in daily irradiance (1-2 weeks). The recorded kinetics of photoacclimation in the Turbinaria-zooxanthella symbiosis is comparable to that of free-living phytoplankton and faster than that of higher plants.