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.

Non-Equilibrium Reaction Rates in the Macroscopic Chemistry Method for DSMC Calculations

The Direct Simulation Monte Carlo (DSMC) method is used to simulate the flow of rarefied gases. In the Macroscopic Chemistry Method (MCM) for DSMC, chemical reaction rates calculated from local macroscopic flow properties are enforced in each cell. Unlike the standard total collision energy (TCE) chemistry model for DSMC, the new method is not restricted to an Arrhenius form of the reaction rate coefficient, nor is it restricted to a collision cross-section which yields a simple power-law viscosity. For reaction rates of interest in aerospace applications, chemically reacting collisions are generally infrequent events and, as such, local equilibrium conditions are established before a significant number of chemical reactions occur. Hence, the reaction rates which have been used in MCM have been calculated from the reaction rate data which are expected to be correct only for conditions of thermal equilibrium. Here we consider artificially high reaction rates so that the fraction of reacting collisions is not small and propose a simple method of estimating the rates of chemical reactions which can be used in the Macroscopic Chemistry Method in both equilibrium and non-equilibrium conditions. Two tests are presented: (1) The dissociation rates under conditions of thermal non-equilibrium are determined from a zero-dimensional Monte-Carlo sampling procedure which simulates ‘intra-modal’ non-equilibrium; that is, equilibrium distributions in each of the translational, rotational and vibrational modes but with different temperatures for each mode; (2) The 2-D hypersonic flow of molecular oxygen over a vertical plate at Mach 30 is calculated. In both cases the new method produces results in close agreement with those given by the standard TCE model in the same highly nonequilibrium conditions. We conclude that the general method of estimating the non-equilibrium reaction rate is a simple means by which information contained within non-equilibrium distribution functions predicted by the DSMC method can be included in the Macroscopic Chemistry Method.

Satellite protection and drag reduction using a purging gas flow

We have used the DSMC method to determine contamination (impingement of atmospheric molecules) and the aerodynamic forces on a cold satellite when a protective “purge gas” is ejected from a sting protruding ahead of the satellite. Forward ejection of the purge gas provides the greatest protection for a given mass of purge gas and the aerodynamic drag can be significantly reduced, thus compensating for the backward reaction from the forward ejection. If the purge gas is ejected backward from the sting (towards the satellite) the ejection provides thrust and the net retarding force can be reduced to zero. Contamination can be reduced and the mass of purging gas is less than the mass of conventional rocket propellant required to maintain the orbit of an unprotected satellite.

Should the health community promote smokeless tobacco (snus) as a harm reduction measure?

Background to the debate: The tobacco control community is divided on whether or not to inform the public that using oral, smokeless tobacco (Swedish snus) is less hazardous to health than smoking tobacco. Proponents of 'harm reduction' point to the Swedish experience. Snus seems to be widely used as an alternative to cigarettes in Sweden, say these proponents, contributing to the low overall prevalence of smoking and smoking-related disease. Harm reduction proponents thus argue that the health community should actively inform inveterate cigarette smokers of the benefits of switching to snus. However, critics of harm reduction say that snus has its own risks, that no form of tobacco should ever be promoted, and that Sweden’s experience is likely to be specific to that culture and not transferable to other settings. Critics also remain deeply suspicious that the tobacco industry will use snus marketing as a 'gateway' to promote cigarettes. In the interests of promoting debate, the authors (who are collaborators on a research project on the future of tobacco control) have agreed to outline the strongest arguments for and against promoting Swedish snus as a form of harm reduction.

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.

Extended GCD and Hermite normal form algorithms via lattice basis reduction

Extended gcd calculation has a long history and plays an important role in computational number theory and linear algebra. Recent results have shown that finding optimal multipliers in extended gcd calculations is difficult. We present an algorithm which uses lattice basis reduction to produce small integer multipliers x(1), ..., x(m) for the equation s = gcd (s(1), ..., s(m)) = x(1)s(1) + ... + x(m)s(m), where s1, ... , s(m) are given integers. The method generalises to produce small unimodular transformation matrices for computing the Hermite normal form of an integer matrix.

Anthocyanin synthesis, growth and nutrient uptake in suspension cultures of strawberry cells

Strawberry (Fragaria ananassa cv. Shikinari) cell suspension cultures carried out in shake flasks for 18 d were closely examined for cell growth, anthocyanin synthesis and the development of pigmented cells in relation to the uptake of carbohydrate, extracellular PO4, NO3, NH4, and calcium. Cell viability, extracellular anthocyanin content, pH and electrical conductivity of the broth were also monitored. The specific growth rate of strawberry cells at exponential phase was 0.27 and 0.28 d(-1) based on fresh and dry weight, respectively. Anthocyanin synthesis was observed to increase continuously to a maximum value of 0.86 mg/g fresh cell weight (FCW) at day 6, and was partially growth-associated. Anthocyanin synthesis was linearly related to the increase in pigmented cell ratio, which increased with time and reached a maximum value of ca. 70% at day 6 due to reduction in cell viability and depletion of substrate. Total carbohydrate uptake was closely associated with increase in cell growth, and glucose was utilized in preference to fructose. Nitrate and ammonia were consumed until 9 d of culture, but phosphate was completely absorbed within 4 d. Calcium was assimilated throughout the growth cycle. After 9 d, cell lysis was observed which resulted in the leakage of intracellular substances and a concomitant pH rise. Anthocyanin was never detected in the broth although the broth became darkly pigmented during the lysis period. This suggests that anthocyanin was synthesized only by viable pigmented cells, and degraded rapidly upon cell death and lysis. Based on the results of kinetic analysis, a model was developed by incorporating governing equations for the ratio of pigmented cells into a Bailey and Nicholson's model. This was verified by comparison with the experimental data. The results suggest Bat the model satisfactorily describes the strawberry cell culture process, and may thus be used for process optimization.

Geochemical mass-balance and oxygen-isotope constrains on silcrete formation and its paleoclimatic implications in Southern Australia

Two geographically distinct silcrete associations are present in southern Australia, inland and eastern; these were sampled in central South Australia and central Victoria, respectively, At each site, both silicified and immediately adjacent unsilicified parent material were collected. Analytical data from these pairs were used to construct isocons, assuming Zr immobility, and to calculate the volume change and amount of silica introduced during silicification, These results, together with whole-rock oxygen isotope compositions, were used to determine the delta(18)O of th, introduced silica, The results show that the eastern silcretes in central Victoria are probably linked genetically to the associated basalts, weathering of which supplied the introduced silica, This conclusion is based on the close spatial connection between the two, as well as the substantial amount of introduced silica in the silcretes (greater than in the inland silcretes), resulting in volume increases in some eastern silcretes, Oxygen isotopic calculations for the silcretes indicate that the silica precipitated from groundwaters at temperatures slightly higher than present conditions. Silcrete formation apparently occurred during the Miocene and Pliocene (basalts in Victoria younger than Pliocene lack associated silcrete) and may reflect the much wetter climate in southeastern Australia at that time. The inland silcretes of central South Australia can be divided into pedogenic (the most common) and groundwater varieties. The pedogenic silcretes, which show typical soil features like columnar and nodular textures, contain moderate amounts of introduced silica that precipitated by evaporation from saline groundwaters, For the groundwater silcretes, which have massive textures and formed at or close to the water table, insufficient data are available to determine the mode of formation. The inland pedogenic silcretes have probably been farming from the Eocene-Miocene to the present, implying that conditions of seasonally high evaporation have occurred in central Australia during this time period. Thus silcrete formation depends on a complex interplay between climate and silica supply, and it is impossible to generalize that the presence of silcrete is indicative of a particular climate. Likewise, the elemental composition of silcretes, particularly Ti content, is not necessarily of climatic significance, Nevertheless, detailed geochemical and oxygen isotopic studies of a silcrete and its parent material can elucidate the mechanisms of silcrete formation, and if evaporation is indicated as a major factor in silcrete formation, then the climate at the time was likely to have been at least seasonally arid.