Search

Abstract Background Rust fungi are biotrophic basidiomycete plant pathogens that cause major diseases on plants and trees world-wide, affecting agriculture and forestry. Their biotrophic nature precludes many established molecular genetic manipulations and lines of research. The generation of genomic resources for these microbes is leading to novel insights into biology such as interactions with the hosts and guiding directions for breakthrough research in plant pathology. Results To support gene discovery and gene model verification in the genome of the wheat leaf rust fungus, Puccinia triticina (Pt), we have generated Expressed Sequence Tags (ESTs) by sampling several life cycle stages. We focused on several spore stages and isolated haustorial structures from infected wheat, generating 17,684 ESTs. We produced sequences from both the sexual (pycniospores, aeciospores and teliospores) and asexual (germinated urediniospores) stages of the life cycle. From pycniospores and aeciospores, produced by infecting the alternate host, meadow rue (Thalictrum speciosissimum), 4,869 and 1,292 reads were generated, respectively. We generated 3,703 ESTs from teliospores produced on the senescent primary wheat host. Finally, we generated 6,817 reads from haustoria isolated from infected wheat as well as 1,003 sequences from germinated urediniospores. Along with 25,558 previously generated ESTs, we compiled a database of 13,328 non-redundant sequences (4,506 singlets and 8,822 contigs). Fungal genes were predicted using the EST version of the self-training GeneMarkS algorithm. To refine the EST database, we compared EST sequences by BLASTN to a set of 454 pyrosequencing-generated contigs and Sanger BAC-end sequences derived both from the Pt genome, and to ESTs and genome reads from wheat. A collection of 6,308 fungal genes was identified and compared to sequences of the cereal rusts, Puccinia graminis f. sp. tritici (Pgt) and stripe rust, P. striiformis f. sp. tritici (Pst), and poplar leaf rust Melampsora species, and the corn smut fungus, Ustilago maydis (Um). While extensive homologies were found, many genes appeared novel and species-specific; over 40% of genes did not match any known sequence in existing databases. Focusing on spore stages, direct comparison to Um identified potential functional homologs, possibly allowing heterologous functional analysis in that model fungus. Many potentially secreted protein genes were identified by similarity searches against genes and proteins of Pgt and Melampsora spp., revealing apparent orthologs. Conclusions The current set of Pt unigenes contributes to gene discovery in this major cereal pathogen and will be invaluable for gene model verification in the genome sequence.

Abstract Background Ustilago maydis is the basidiomycete fungus responsible for common smut of corn and is a model organism for the study of fungal phytopathogenesis. To aid in the annotation of the genome sequence of this organism, several expressed sequence tag (EST) libraries were generated from a variety of U. maydis cell types. In addition to utility in the context of gene identification and structure annotation, the ESTs were analyzed to identify differentially abundant transcripts and to detect evidence of alternative splicing and anti-sense transcription. Results Four cDNA libraries were constructed using RNA isolated from U. maydis diploid teliospores (U. maydis strains 518 × 521) and haploid cells of strain 521 grown under nutrient rich, carbon starved, and nitrogen starved conditions. Using the genome sequence as a scaffold, the 15,901 ESTs were assembled into 6,101 contiguous expressed sequences (contigs); among these, 5,482 corresponded to predicted genes in the MUMDB (MIPS Ustilago maydis database), while 619 aligned to regions of the genome not yet designated as genes in MUMDB. A comparison of EST abundance identified numerous genes that may be regulated in a cell type or starvation-specific manner. The transcriptional response to nitrogen starvation was assessed using RT-qPCR. The results of this suggest that there may be cross-talk between the nitrogen and carbon signalling pathways in U. maydis. Bioinformatic analysis identified numerous examples of alternative splicing and anti-sense transcription. While intron retention was the predominant form of alternative splicing in U. maydis, other varieties were also evident (e.g. exon skipping). Selected instances of both alternative splicing and anti-sense transcription were independently confirmed using RT-PCR. Conclusion Through this work: 1) substantial sequence information has been provided for U. maydis genome annotation; 2) new genes were identified through the discovery of 619 contigs that had previously escaped annotation; 3) evidence is provided that suggests the regulation of nitrogen metabolism in U. maydis differs from that of other model fungi, and 4) Alternative splicing and anti-sense transcription were identified in U. maydis and, amid similar observations in other basidiomycetes, this suggests these phenomena may be widespread in this group of fungi. These advances emphasize the importance of EST analysis in genome annotation.

Cold weather operability is sometimes a limiting factor in the use of biodiesel blends for transportation. Regional temperature variability can therefore influence biodiesel adoption, with potential economic and environmental implications. This study assesses present and future biodiesel cold weather operability limits in North America according to temperature data from weather stations, atmospheric reanalysis, and global climate models with highest resolution over Ontario, Canada. Future temperature projections using the RCP8.5 climate change scenario show increases in the potential duration for certain seasonal fuel blends. For example, biodiesel blends whose cloud point temperature is -9 °C may expand their duration by 3-7% in North America for nonwinter seasons according to projections for 2040. Cloud point specifications among supply orbits in Ontario increase up to +6 °C during nonwinter seasons, with most increases observed in Fall and Spring. In winter, however, the modeling suggests no change in Ontario cloud point specifications because the coldest temperatures by mid-century are not significantly warmer than the past climate normal according to our climate simulations. This study provides a quantitative analysis on biodiesel usage scenarios under a changing climate, including Ontario region geographic temperature clusters that could prove useful for biodiesel blend-related decision-making.

A key motivation behind the development and adoption of industrial biotechnology is the reduction of negative environmental impacts. However, accurately assessing these impacts remains a formidable task. Environmental impacts of industrial biotechnology may be significant across a number of categories that include, but may not be limited to, nonrenewable resource depletion, water withdrawals and consumption, climate change, and natural land transformation/occupation. In this chapter, we highlight some key environmental issues across two broad areas: (a) processes that use biobased feedstocks and (b) industrial activity that is supported by biological processes. We also address further issues in accounting for related environmental impacts such as geographic and temporal scope, co-product management, and uncertainty and variability in impacts. Case studies relating to (a) lignocellulosic ethanol, (b) biobased plastics, and (c) enzyme use in the detergent industry are then presented, which illustrate more specific applications. Finally, emerging trends in the area of environmental impacts of biotechnology are discussed.

Kt/V is a nondimensional number and a scaling parameter that has, with arbitrary definitions, been recast as a measure of dialysis by Gotch and Lysaght. This editorial discusses the concept of nondimensional numbers within the context of dialysis measurement, modeling, and medical evidence. It concludes that Gotch's Kt/V, Lysaght's Kt/V, and standardized Kt/V are not well suited to measure dialysis. An ideal dialysis measure would be proportional to toxins cleared by the kidney, portable with regard to dialysis modality, practical (largely devoid of calculations) and reflective of the pathology/physiology., Competing Interests: Declaration of Conflicting Interests: The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Lignocellulosic ethanol has potential for lower life cycle greenhouse gas emissions compared to gasoline and conventional grain-based ethanol. Ethanol production 'pathways' need to meet economic and environmental goals. Numerous life cycle assessments of lignocellulosic ethanol have been published over the last 15 years, but gaps remain in understanding life cycle performance due to insufficient data, and model and methodological issues. We highlight key aspects of these issues, drawing on literature and a case study of corn stover ethanol. Challenges include the complexity of feedstock/ecosystems and market-mediated aspects and the short history of commercial lignocellulosic ethanol facilities, which collectively have led to uncertainty in GHG emissions estimates, and to debates on LCA methods and the role of uncertainty in decision making., (Copyright © 2016 Elsevier Ltd. All rights reserved.)

This paper aims to comprehensively distinguish among the merits of different vehicles using a common primary energy source. In this study, we consider compressed natural gas (CNG) use directly in conventional vehicles (CV) and hybrid electric vehicles (HEV), and natural gas-derived electricity (NG-e) use in plug-in battery electric vehicles (BEV). This study evaluates the incremental life cycle air emissions (climate change and human health) impacts and life cycle ownership costs of non-plug-in (CV and HEV) and plug-in light-duty vehicles. Replacing a gasoline CV with a CNG CV, or a CNG CV with a CNG HEV, can provide life cycle air emissions impact benefits without increasing ownership costs; however, the NG-e BEV will likely increase costs (90% confidence interval: $1000 to $31 000 incremental cost per vehicle lifetime). Furthermore, eliminating HEV tailpipe emissions via plug-in vehicles has an insignificant incremental benefit, due to high uncertainties, with emissions cost benefits between -$1000 and $2000. Vehicle criteria air contaminants are a relatively minor contributor to life cycle air emissions impacts because of strict vehicle emissions standards. Therefore, policies should focus on adoption of plug-in vehicles in nonattainment regions, because CNG vehicles are likely more cost-effective at providing overall life cycle air emissions impact benefits.

Our study evaluates life cycle energy use and GHG emissions of lignocellulosic ethanol and bioelectricity use in U.S. light-duty vehicles. The well-to-pump, pump-to-wheel, and vehicle cycle stages are modeled. All ethanol (E85) and bioelectricity pathways have similar life cycle fossil energy use (~ 100 MJ/100 vehicle kilometers traveled (VKT)) and net GHG emissions (~5 kg CO2eq./100 VKT), considerably lower (65-85%) than those of reference gasoline and U.S. grid-electricity pathways. E85 use in a hybrid vehicle and bioelectricity use in a fully electric vehicle also have similar life cycle biomass and total energy use (~ 350 and ~450 MJ/100 VKT, respectively); differences in well-to-pump and pump-to-wheel efficiencies can largely offset each other. Our energy use and net GHG emissions results contrast with findings in literature, which report better performance on these metrics for bioelectricity compared to ethanol. The primary source of differences in the studies is related to our development of pathways with comparable vehicle characteristics. Ethanol or vehicle electrification can reduce petroleum use, while bioelectricity may displace nonpetroleum energy sources. Regional characteristics may create conditions under which either ethanol or bioelectricity may be the superior option; however, neither has a clear advantage in terms of GHG emissions or energy use.

Background: Non-dimensional analysis is a powerful approach that can be applied to multivariate problems to better understand their behaviour and interpret complex interactions of variables. It is has not been rigorously applied to the parameters that define renal dialysis treatments and may provide insight into the planning of hemodialysis treatments., Methods: Buckingham's non-dimensional approach was applied to the parameters that define hemodialysis treatments. Non-dimensional groups were derived with knowledge of a mass transfer model and independent of it. Using a mass transfer model, the derived non-dimensional groups were plotted to develop an understanding of key relationships governing hemodialysis and toxin profiles in patients with end-stage renal disease., Results: Three non-dimensional groups are sufficient to describe hemodialysis, if there is no residual renal function (RRF). The non-dimensional groups found represent (1) the number of half-lives that characterize the mass transfer, (2) the toxin concentration divided by the rise in toxin concentration without dialysis for the cycle time (the inverse of the dialysis frequency), and (3) the ratio of dialysis time to the cycle time. If there is RRF, one additional non-dimensional group is needed (the ratio between cycle time and intradialytic elimination rate constant). Alternate non-dimensional groups can be derived from the four unique groups., Conclusions: Physical interpretation of the non-dimensional groups allows for greater insight into the parameters that determine dialysis effectiveness. This technique can be applied to any toxin and facilitates a greater understanding of dialysis treatment options. Quantitative measures of dialysis adequacy should be based on dimensional variables.

In this article, the mechanisms of thermoinactivation of glucose isomerase (GI) from Streptomyces rubiginosus (in soluble and immobilized forms) were investigated, particularly the contributions of thiol oxidation of the enzyme's cysteine residue and a "Maillard-like" reaction between the enzyme and sugars in high fructose corn syrup (HFCS). Soluble GI (SGI) was successfully immobilized on silica gel (13.5 microm particle size), with an activity yield between 20 and 40%. The immobilized GI (IGI) has high enzyme retention on the support during the glucose isomerization process. In batch reactors, SGI (half-life =145 h) was more stable than IGI (half-life =27 h) at 60 degrees C in HFCS, whereas at 80 degrees C, IGI (half-life =12 h) was more stable than SGI (half-life =5.2 h). IGI was subject to thiol oxidation at 60 degrees C, which contributed to the enzyme's deactivation. IGI was subject to thiol oxidation at 80 degrees C, but this did not contribute to the deactivation of the enzyme. SGI did not undergo thiol oxidation at 60 degrees C, but at 80 degrees C SGI underwent severe precipitation and thiol oxidation, which caused the enzyme to deactivate. Experimental results show that immobilization suppresses the destabilizing effect of thiol oxidation on GI. A "Maillard-like" reaction between SGI and the sugars also caused SGI thermoinactivation at 60, 70, and 80 degrees C, but had minimal effect on IGI. At 60 and 80 degrees C, IGI had higher thermostability in continuous reactors than in batch reactors, possibly because of reduced contact with deleterious compounds in HFCS.

Studies were conducted on maltodextrin saccharification and on simultaneous saccharification and fermentation (SSF) with various commercial glucoamylases. In kinetics studies, none of the glucoamylases were able to completely convert maltodextrin into glucose. Typically, about 85% conversion was obtained, and glucose yields were about 75%. Typically, the kinetics were biphasic, with 1 h of rapid conversion, then a significant reduction in rate. Data were consistent with strong product inhibition and/or enzyme inactivation. Some glucoamylases followed first-order kinetics, initially slower at dextrin conversion, but eventually achieving comparable conversion and glucose concentrations. Most of the glucoamylases were more active at 55 degrees C than at 35 degrees C, but pH had little effect on activity. Screening studies in an SSF system demonstrated little difference between the glucoamylases, with a few exceptions. Subsequent targeted studies showed clear differences in performance, depending on the fermentation temperature and yeast used, suggesting that these are key parameters that would guide the selection of a glucoamylase.

The performance of cellulase and amylase immobilized on siliceous supports was investigated. Enzyme uptake onto the support depended on the enzyme source and immobilization conditions. For amylase, the uptake ranged between 20 and 60%, and for cellulase, 7-10%. Immobilized amylase performance was assessed by batch kinetics in 100-300 g/L of corn flour at 65 degrees C. Depending on the substrate and enzyme loading, between 40 and 60% starch conversion was obtained. Immobilized amylase was more stable than soluble amylase. Enzyme samples were preincubated in a water bath at various temperatures, then tested for activity. At 105 degrees C, soluble amylase lost approximately 55% of its activity, compared with approximately 30% loss for immobilized amylase. The performance of immobilized cellulase was evaluated from batch kinetics in 10 g/L of substrate (shredded wastepaper) at 55 degrees C. Significant hydrolysis of the wastepaper was also observed, indicating that immobilization does not preclude access to and hydrolysis of insoluble cellulose.

The degradation of phenol by tyrosinase immobilized on chemically modified sodium aluminosilicate (NaA), calcium aluminosilicate (CaA), and silica gel was studied. Phenol conversion by immobilized tyrosinase ranged between approximately 15% and 60%, depending upon the initial phenol concentration, pH, and enzyme loading. Tyrosinase immobilized on CaA and on NaA could be re-used repeatedly without any decrease in performance. However, in studies at pH 8.0, significant enzyme inhibition was observed, since phenol conversion was rapid for approximately 20 min, then reached a plateau. The inhibition was reversible; activity was restored upon placing the immobilized enzyme in fresh substrate. Reducing the pH to 6.8 from 8.0 led to higher conversion of phenol, and decreased the inhibition of the immobilized enzyme.

The inactivation of the catecholase activity of mushroom tyrosinase was investigated under nonaqueous conditions. The enzyme was immobilized on glass beads, and assays were conducted in chloroform, toluene, amyl acetate, isopropyl ether, and butanol. The reaction components were pre-equilibrated for 2 weeks with a saturated salt solution at a water activity of 0.90. The initial reaction velocity varied between 1.3 x 10(3) mol product/((mol enzyme)(min)) in toluene and 8.7 x 10(3) mol product/((mol enzyme)(min)) in amyl acetate. The turnover number varied between 8.1 x 10(3) mol product/mol enzyme in toluene and 7.2 x 10(4) mol product/mol enzyme in amyl acetate. In each solvent, the tyrosinase reaction inactivation parameters were represented by a probabilistic model. Changes in the probability of inactivation were followed throughout the course of the reaction using a second model which relates the reaction velocity to the amount of product formed. These models reveal that the inactivation rate of tyrosinase decreases as the reaction progresses, and that the inactivation kinetics are independent of the quinone concentration in toluene, chloroform, butanol, and amyl acetate. Significant effects of quinone concentration were, however, observed in isopropyl ether. The likelihood of inactivation of the enzyme was found to be greatest toward the beginning of the reaction. In the latter phase of the reaction, inactivation probability was less and tended to remain constant until the completion of the reaction., (Copyright 1999 John Wiley & Sons, Inc.)

The production of L-DOPA immobilized on chemically modified nylon 6,6 membranes was studied in a batch reactor. Tyrosinase was immobilized on nylon using glutaraldehyde as a crosslinking agent. The effects of membrane pore size and glutaraldehyde concentration upon enzyme uptake and L-DOPA production were investigated. Enzyme uptake was unaffected by glutaraldehyde concentration; approximately 70% uptake was observed when 25% w/v (group 1), 5% (group 2), and 3% (group 3) glutaraldehyde were used, indicating that glutaraldehyde was in excess. Similarly, uptake was the same for membranes with 0.20 and 10 microm pore sizes.Membranes produced using different levels of glutaraldehyde exhibited dramatically different capacities for L-DOPA production, despite the fact that enzyme uptake was equivalent. Membranes from groups 2 and 3 (5% and 3% glutaraldehyde) produced L-DOPA at a rate of 1.70 mg L(-1) h(-1) over 170 h in a 500-mL batch reactor. However, no free L-DOPA was detected when group 1 membranes were used. Experimental evidence suggests that L-DOPA was produced, but remained bound to these membranes via excess glutaraldehyde left over from the immobilization process. Membrane pore size also effected L-DOPA production; less production was observed when 10-microm membranes were used, despite equivalent enzyme uptake. The observed difference in production may be due to differences in the pore density on the two types of membranes which could affect the access of the substrate to the immobilized enzyme.The results of these studies indicate that tyrosinase can be effectively immobilized on nylon 6,6. L-DOPA production was optimal when 0.20-microm-pore-size membranes were activated with 3-5% glutaraldehyde. Stability studies indicated a 20% reduction in activity over 14 days when the immobilized enzyme was used under turnover conditions.

Controlled release ocular inserts have been found to increase the amount of drug which is absorbed into the aqueous humour when compared to eyedrops. Systemic absorption following delivery using a controlled release insert has been found to be dependent on the release rate of the insert. The objective of this study was to determine if ocular inserts affect drug absorption into other ocular tissues such as the conjunctiva and iris-ciliary body. Ocular absorption studies were performed using albino rabbits and ethylene-vinyl acetate controlled release devices containing timolol maleate. A compartmental model previously developed to simulate ocular absorption following eyedrop administration was modified and used to simulate these experiments. The conjunctival absorption coefficient calculated by the model and the AUC of the conjunctiva per mumol of delivered drug were found to be 2.7 and 42 times higher, respectively, for the ocular insert as compared to eyedrop administration. The increased conjunctiva absorption was likely the result of reduced tear mixing, which caused a high local concentration of timolol between the insert and the conjunctiva. The AUC of the iris-ciliary body per mumol of delivered drug was found to be 24 times higher for the ocular inserts as compared to eyedrop administration. The AUC of the iris-ciliary body was found to be 1.4 times higher than the AUC of the aqueous humour for eyedrop administration, but 9 times greater for delivery via the ocular inserts. Thus, the increased absorption into the iris-ciliary body and aqueous humour observed for ocular inserts is partially the result of an increase in the amount of drug which enters these tissues via penetration across the conjunctiva and sclera.

Studies were conducted in rabbits to determine the ocular distribution and elimination of cyclosporine, with the objective of developing a comprehensive pharmacokinetic model. Following a bolus dose into the anterior chamber, drug levels were measured in the aqueous humor, cornea, iris/ciliary body, lens, sclera, and conjunctiva. Cyclosporine was rapidly eliminated from the aqueous, but drug levels in ocular tissues persisted for in excess of 48 hours, particularly in the cornea and iris/ciliary body. The terminal elimination half life from these tissues was 45 hr and 30 hr, respectively, providing evidence that these tissues could act as a reservoir for the drug. It was found that a compartmental model accurately described the experimental data. A single compartment was used for each of the tissues and fluids sampled, except for the cornea, which was subdivided into two compartments, representing its tissue and aqueous regions.

Studies were conducted to evaluate the efficacy of direct injections of cyclosporine (CsA) into the anterior chamber for the prevention of corneal allograft rejection in Dutch Belted rabbits. The mean survival time (MST) of grafts progressively increased from 50 to 89 days as the CsA concentration in the dose was increased from 1 to 10 mg/mL. Injection of 30 microL of 20 mg/mL CsA in olive oil prolonged graft survival to beyond 125 days without any signs of rejection. By comparison, the MST of allografts in control animals which received no therapy was 32 +/- 5 days, and the MST in animals administered a placebo of olive oil only was 31 +/- 4 days. The observed concentration dependence of the MST on CsA concentration is likely related to the time over which the drug delivery rate provides sufficient drug to achieve a therapeutic concentration in the aqueous humor; these studies suggest that the minimum delivery rate to the anterior chamber is between 200 and 325 ng/day. The efficacy of CsA was due to local delivery, and was likely not a systemic effect, because CsA was not detected in the systemic circulation at any time. This indicates that direct delivery of CsA to the eye can be useful in prolonging corneal graft survival, while minimizing systemic side effects. Separate experiments revealed that episodes of advanced rejection could not be reversed by a 30 microL dose of 20 mg/mL CsA to the anterior chamber, indicating the importance of avoiding long periods of subtherapeutic dosing.

The primary route into the eye for many drugs is transcorneal permeation. A better understanding of the mechanisms involved in transcorneal permeation could lead to improvements in drug dosage forms or the development of drug delivery devices which enhance the ocular bioavailability of drugs. A corneal permeation model has been developed which can be used to study the mechanisms involved in corneal permeation. The model uses five compartments in series to simulate the tear film, epithelium, stroma, endothelium and aqueous humour. These tissues were assumed to be adequately represented by plane sheet barriers of physiological thickness. The tear film was assumed to be perfectly mixed and the stroma completely stagnant. Due to inadequate knowledge of the hydrodynamics of the aqueous humour, both stagnant and perfectly mixed extremes were studied. The four routes of drug loss which were considered the most significant and therefore included in the model were lacrimal drainage, conjunctival absorption, aqueous drainage and iris-ciliary body absorption. The equilibrium that can exist between the ionic and non-ionic forms of a drug was found to be an important step in the mechanism of transcorneal permeation. Including the equilibrium condition in the model resulted in aqueous humour drug levels that were over 50 times higher than the levels predicted by a model which did not use the equilibrium mechanism. A relationship between the lipophilicity of each of the two drug forms and its permeability in each layer of the cornea was used in the model. The model was used to predict aqueous humour drug concentrations resulting from a constant release of timolol into the tear film or from the application of timolol, levobunolol and pilocarpine eyedrops. The model produced transient aqueous humour drug levels that closely followed experimental in vivo data from literature. Using the model, it was also possible to predict the amount of instilled drug that is lost through each of the four elimination routes of the eye.

The uptake and metabolism of lipophilic compounds by the liver were studied by administering a model compound, lidocaine, to the isolated rat liver. Lidocaine was continuously infused into the liver until steady state was reached. Subsequent step changes in the inlet concentration were used to obtain information on rates of cellular uptake and release and to assess the extent of mixing within the organ. A simple heterogeneous model combining mass transfer and enzyme reactions was required to simulate the effluent levels of lidocaine and two primary metabolites, monoethylglycinexylidide and 3-hydroxylidocaine. The rate constants for uptake and release of lidocaine were 1200 and 46 min-1, respectively. The rate-limiting step was intracellular reaction, with a rate constant of 0.49 min-1. Although the rate of lidocaine uptake was fast, it was 50 times slower than the rate of facilitated uptake of galactose, a fact suggesting passive transport of lidocaine between the tissue and the vasculature. The rates of mass transfer of lidocaine and its metabolites differed, but the ratios of the rate of uptake to the rate of release were the same. The results suggested that all three species had an affinity for the cellular region of the liver; concentrations in tissue were approximately five times greater than concentrations in effluent. Because of the large capacity of the organ for uptake of lidocaine and its metabolites, concentrations from washout experiments were controlled by linear mass transfer from the tissue.(ABSTRACT TRUNCATED AT 250 WORDS)

The liver is, by nature, heterogeneous. It contains a complex vascular network for blood flow and a stationary phase consisting of enzymes within parenchymal cells. Several physiological processes, therefore, may combine to give observed ranges in drug elimination. Net changes in concentration are a consequence of a series of steps: uptake of substrate into liver cells, enzymatic reactions within the cells, release of metabolites and unconverted substrate from the cells into the sinusoids, and the net flow of the perfusing medium in the vasculature. In addition, substrate binding to proteins in the blood and in the liver can influence hepatic elimination. An understanding of each of these processes is necessary to fully comprehend the overall process of drug elimination, and these processes must be accounted for, either individually or by grouping and approximation, if a model for drug elimination is to be developed. Existing models of hepatic elimination may be classified according to their treatment of mixing within the vasculature and whether or not the model explicitly accounts for mass transfer between the heterogeneous phases of the liver. Four major classes may be defined: 1. Nonparametric homogeneous models, which assume that either complete mixing or no mixing occurs within the vasculature of the organ. 2. Homogeneous mixing models, which allow for a range of mixing phenomena. 3. Heterogeneous micromixing models, which allow for mass transport between the cells and vasculature and describe mixing within the vasculature on a microscopic level. 4. Heterogeneous compartmental models, which also describe interphase mass transfer but assume complete mixing on a microscopic level, and therefore use a time and spatially averaged approach to model mixing. The utility of these models of hepatic elimination will be critically assessed based upon (1) their ability to account for the influence of the aforementioned physiological processes upon elimination; (2) the data requirements of the model, in addition to its mathematical complexity and ease of use; and (3) the range of compounds and metabolites which may be described using the model.

The steady state kinetics of lidocaine and its metabolites were modeled using nonlinear elimination pathways for multiple enzymes. The main metabolites, monoethylglycinexylidide and 3-hydroxy-lidocaine, were infused in the absence of lidocaine to measure the kinetic parameters for secondary elimination. Data from continuous perfusion of lidocaine in the isolated perfused rat liver at concentrations ranging from 9.6 to 278 microM (N = 16) were used to calculate the kinetic parameters for formation of the main metabolites. The elimination of lidocaine in the liver was approximated by the well stirred model. The whole liver study gave higher elimination rates than were predicted from microsomal studies. The major pathways for elimination of lidocaine in the rat were deethylation and hydroxylation, and subsequent elimination along these pathways accounted for the poor material balance at low dosage levels. The observed competitive inhibition of hydroxylation was in agreement with the predictions of the model.

The effects of lidocaine on hepatic enzyme activity were studied using the isolated perfused rat liver. The in vivo liver activity was examined by infusing lidocaine via the jugular vein, followed by organ isolation and drug perfusion 24 h later. The liver was studied in vitro by perfusing the organ with lidocaine until steady state was reached, then allowing the drug and metabolites to wash out of the organ, followed by a second infusion of lidocaine to probe enzyme activity. In both types of experiments, pretreatment with lidocaine caused a reduction in deethylation, and led to a more rapid attainment of steady state. The experimental concentration-time profiles and literature data were successfully described by a mathematical model.

Three models for the elimination of lidocaine in the isolated perfused rat liver were used to simulate the time course of both lidocaine and its metabolites in a single-pass perfusion system: well stirred, parallel tube, and a two-compartment model to test the effects of enzyme heterogeneity. All models included multiple enzymes and multiple metabolic pathways, as well as varying degrees of tissue binding. Although the well stirred and parallel tube models gave qualitatively different results, neither model predicted that the concentration of monoethylglycinexylidide would pass through a maximum and then decline to a lower steady state value, as observed in continuous perfusion experiments. Although each of the three models tested would give reasonable agreement with steady state observations, the test of the time-dependent behavior of both lidocaine and monoethylglycinexylidide was more discriminating. Each model gave characteristic predictions for the time course of the metabolites.