Structure-hepatic disposition relationships for cationic drugs in isolated perfused rat livers: Transmembrane exchange and cytoplasmic binding process

This work studied the structure-hepatic disposition relationships for cationic drugs of varying lipophilicity using a single-pass, in situ rat liver preparation. The lipophilicity among the cationic drugs studied in this work is in the following order: diltiazem. propranolol. labetalol. prazosin. antipyrine. atenolol. Parameters characterizing the hepatic distribution and elimination kinetics of the drugs were estimated using the multiple indicator dilution method. The kinetic model used to describe drug transport (the two-phase stochastic model) integrated cytoplasmic binding kinetics and belongs to the class of barrier-limited and space-distributed liver models. Hepatic extraction ratio (E) (0.30-0.92) increased with lipophilicity. The intracellular binding rate constant (k(on)) and the equilibrium amount ratios characterizing the slowly and rapidly equilibrating binding sites (K-S and K-R) increase with the lipophilicity of drug (k(on) : 0.05-0.35 s(-1); K-S : 0.61-16.67; K-R : 0.36-0.95), whereas the intracellular unbinding rate constant (k(off)) decreases with the lipophilicity of drug (0.081-0.021 s(-1)). The partition ratio of influx (k(in)) and efflux rate constant (k(out)), k(in)/k(out), increases with increasing pK(a) value of the drug [from 1.72 for antipyrine (pK(a) = 1.45) to 9.76 for propranolol (pK(a) = 9.45)], the differences in k(in/kout) for the different drugs mainly arising from ion trapping in the mitochondria and lysosomes. The value of intrinsic elimination clearance (CLint), permeation clearance (CLpT), and permeability-surface area product (PS) all increase with the lipophilicity of drug [CLint (ml . min(-1) . g(-1) of liver): 10.08-67.41; CLpT (ml . min(-1) . g(-1) of liver): 10.80-5.35; PS (ml . min(-1) . g(-1) of liver): 14.59-90.54]. It is concluded that cationic drug kinetics in the liver can be modeled using models that integrate the presence of cytoplasmic binding, a hepatocyte barrier, and a vascular transit density function.

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.

Studies of solute self-association by sedimentation equilibrium: allowance for effects of thermodynamic non-ideality beyond the consequences of nearest-neighbor interactions

A sedimentation equilibrium study of a-chymotrypsin self-association in acetate-chloride buffer, pH 4.1 I 0.05, has been used to illustrate determination of a dimerization constant under conditions where thermodynamic non-ideality is manifested beyond the consequences of nearest-neighbor interactions. Because the expressions for the experimentally determinable interaction parameters comprise a mixture of equilibrium constant and excluded volume terms, the assignment of reasonable magnitudes to the relevant virial coefficients describing non-associative cluster formation is essential for the evaluation of a reliable estimate of the dimerization constant. Determination of these excluded volume parameters by numerical integration over the potential-of-mean-force is shown to be preferable to their calculation by approximate analytical solutions of the integral for this relatively small enzyme monomer with high net charge (+ 10) under conditions of low ionic strength (0.05 M). (C) 2001 Elsevier Science B.V. All rights reserved.

Skin-friction measurements in a supersonic combustor with crossflow fuel injection

Shock-tunnel experiments have been performed to measure the effect on skin-friction drag in a supersonic combustor of flow disturbances induced by hydrogen fuel injection transverse to the airstream. Constant-area, circular cross section combustors of lengths varying up to 0.52 m were employed. The experiments were done at a stagnation enthalpy of 7.2 MJ . kg(-1) and a Mach number of 4.3, with a boundary layer that was turbulent downstream of the 0.14-m station in the combustors. Combustor skin-friction drag was measured by a method based on the stress wave force balance, the method being validated by agreement between fuel-off skin-friction drag measurements and predictions using existing skin-friction theories. When fuel was injected, it was found that the drag remained at fuel-off values. Thus, the streamwise vortices and other flow disturbances induced by the fuel injection, mixing, and combustion, which are expected to be present in a scramjet combustor, did not influence the skin-friction drag of the combustors.

Three-dimensional orthotropic viscoelastic finite element model of a human ligament

Ligaments undergo finite strain displaying hyperelastic behaviour as the initially tangled fibrils present straighten out, combined with viscoelastic behaviour (strain rate sensitivity). In the present study the anterior cruciate ligament of the human knee joint is modelled in three dimensions to gain an understanding of the stress distribution over the ligament due to motion imposed on the ends, determined from experimental studies. A three dimensional, finite strain material model of ligaments has recently been proposed by Pioletti in Ref. [2]. It is attractive as it separates out elastic stress from that due to the present strain rate and that due to the past history of deformation. However, it treats the ligament as isotropic and incompressible. While the second assumption is reasonable, the first is clearly untrue. In the present study an alternative model of the elastic behaviour due to Bonet and Burton (Ref. [4]) is generalized. Bonet and Burton consider finite strain with constant modulii for the fibres and for the matrix of a transversely isotropic composite. In the present work, the fibre modulus is first made to increase exponentially from zero with an invariant that provides a measure of the stretch in the fibre direction. At 12% strain in the fibre direction, a new reference state is then adopted, after which the material modulus is made constant, as in Bonet and Burton's model. The strain rate dependence can be added, either using Pioletti's isotropic approximation, or by making the effect depend on the strain rate in the fibre direction only. A solid model of a ligament is constructed, based on experimentally measured sections, and the deformation predicted using explicit integration in time. This approach simplifies the coding of the material model, but has a limitation due to the detrimental effect on stability of integration of the substantial damping implied by the nonlinear dependence of stress on strain rate. At present, an artificially high density is being used to provide stability, while the dynamics are being removed from the solution using artificial viscosity. The result is a quasi-static solution incorporating the effect of strain rate. Alternate approaches to material modelling and integration are discussed, that may result in a better model.

A study of the oxidation of two aromatic polyimide fluoropolymers containing phenylphosphine oxide by a water plasma and gamma-radiolysis in air

The oxidation of two fluorinated polyimides containing phenylphosphine oxide units, TOR-RC and TOR-RC ODPA, have been studied at 300 K for treatment by a water plasma and gamma -radiolysis in air. The changes in the O 1s/C 1s ratios obtained from x-ray photoelectron spectroscopy (XPS) analysis showed that for exposure to the water plasma the ratio increases at short exposure times and then levels to a constant value. Evidence for the formation of phosphate species was also obtained from the XPS analyses. Similar observations were made for gamma -radiolysis of the polymers in air. The polymers containing phenylphosphine oxide were found to be more resistant to oxidation in the water plasma than Kapton(R). Radiolysis of the polymers in air to high doses were also accompanied by a red shift in the visible absorption spectra.

Free-radical bulk polymerization of styrene: ESR and near-infrared spectroscopic study of the entire conversion range

The free radical polymerization of styrene in bulk was monitored by ESR and FT near-infrared spectroscopy at 70°C for a series of concentrations of the initiator, dimethyl 2,2′-azobis(isobutyrate). In order to obtain detailed kinetic information over the intire conversion range, and the gel effect range in particular, conversion and free radical concentration data points were accumulated with exceptionally short time intervals. The polystyrene radical concentration ([St•]) went through a sharp maximum at the gel effect, a feature that has hitherto escaped observation due to the rapid concentration changes in the gel effect range relative to the data point time intervals of previous studies. Temperature measurements throughout the polymerization were employed to calculate that a temperature increase was not the cause of the [St•] maximum, which thus appeares to be a genuine feature of the gel effect of this system under isothermal conditions. The propagation rate constant (kp) as a function of monomer conversion exhibited a marked dependence on initiator concentration at high monomer conversion; the sharp decrease in kp with increasing conversion was shifted to higher conversions with increasing initiator concentration.

The behaviour of phenylphosphine oxide containing polyimides in a water plasma following gamma-radiolysis

The surface oxidation of two polyimides containing fluorinated phenylphosphine oxide units, TOR-RC and TOR-RC ODPA, have been studied by (XPS) spectroscopy following gamma -radiolysis under vacuum or in air and subsequent treatment in a water plasma. The changes in the O 1s/C 1s ratios obtained from (XPS) analysis showed that on exposure to the water plasma the ratio increases and then levels to a constant value which is similar to that found for exposure to the plasma without prior gamma -radiation treatment. Evidence for the formation of phosphate species was also obtained from the (XPS) analyses. (C) 2001 Elsevier Science Ltd. All rights reserved.

Hepatocyte [H-3]-palmitate uptake: effect of albumin surface charge modification

The role of plasma proteins on the cellular uptake of lipophilic substrates has perplexed investigators for many years. We tested the hypothesis that an ionic interaction between the protein-ligand complex and hepatocyte surface may be responsible for supplying more ligand to the cell for uptake. The surface-charged groups on albumin were modified to yield proteins having a range of isoelectric points (ALB, ALBs, ALBm, ALBe had values of 4.8-5.0, 4.5-4.7, 3.0-3.5, 8.4-8.6, respectively). [H-3]-Palmitate uptake studies were performed with adult rat hepatocyte suspensions using similar unbound ligand fractions in the presence of the different binding proteins. Mass spectrometry, isoelectric focusing (pI), and heptane : water partitioning were used to determine protein molecular weight, pI, and protein-palmitate equilibrium binding constant, respectively. Hepatocyte [H-3]-palmitate clearance in the presence of ALBs and ALBm were significantly lower (p < 0.05) than ALB, whereas [H-3]-palmitate clearance in the presence of ALBe was significantly higher (p < 0.05) than ALB. The data were consistent with the notion that ionic interactions between extracellular protein-ligand complexes and the hepatocyte surface facilitate the uptake of long-chain fatty acids.

Time-dependent master equation simulation of complex elementary reactions in combustion: Application to the reaction of (CH2)-C-1 with C2H2 from 300-2000 K

Computational simulations of the title reaction are presented, covering a temperature range from 300 to 2000 K. At lower temperatures we find that initial formation of the cyclopropene complex by addition of methylene to acetylene is irreversible, as is the stabilisation process via collisional energy transfer. Product branching between propargyl and the stable isomers is predicted at 300 K as a function of pressure for the first time. At intermediate temperatures (1200 K), complex temporal evolution involving multiple steady states begins to emerge. At high temperatures (2000 K) the timescale for subsequent unimolecular decay of thermalized intermediates begins to impinge on the timescale for reaction of methylene, such that the rate of formation of propargyl product does not admit a simple analysis in terms of a single time-independent rate constant until the methylene supply becomes depleted. Likewise, at the elevated temperatures the thermalized intermediates cannot be regarded as irreversible product channels. Our solution algorithm involves spectral propagation of a symmetrised version of the discretized master equation matrix, and is implemented in a high precision environment which makes hitherto unachievable low-temperature modelling a reality.

Application of Petrov-Galerkin methods to transient boundary value problems in chemical engineering: adsorption with steep gradients in bidisperse solids

Petrov-Galerkin methods are known to be versatile techniques for the solution of a wide variety of convection-dispersion transport problems, including those involving steep gradients. but have hitherto received little attention by chemical engineers. We illustrate the technique by means of the well-known problem of simultaneous diffusion and adsorption in a spherical sorbent pellet comprised of spherical, non-overlapping microparticles of uniform size and investigate the uptake dynamics. Solutions to adsorption problems exhibit steep gradients when macropore diffusion controls or micropore diffusion controls, and the application of classical numerical methods to such problems can present difficulties. In this paper, a semi-discrete Petrov-Galerkin finite element method for numerically solving adsorption problems with steep gradients in bidisperse solids is presented. The numerical solution was found to match the analytical solution when the adsorption isotherm is linear and the diffusivities are constant. Computed results for the Langmuir isotherm and non-constant diffusivity in microparticle are numerically evaluated for comparison with results of a fitted-mesh collocation method, which was proposed by Liu and Bhatia (Comput. Chem. Engng. 23 (1999) 933-943). The new method is simple, highly efficient, and well-suited to a variety of adsorption and desorption problems involving steep gradients. (C) 2001 Elsevier Science Ltd. All rights reserved.

Solution techniques for transport problems involving steep concentration gradients: application to noncatalytic fluid solid reactions

Some efficient solution techniques for solving models of noncatalytic gas-solid and fluid-solid reactions are presented. These models include those with non-constant diffusivities for which the formulation reduces to that of a convection-diffusion problem. A singular perturbation problem results for such models in the presence of a large Thiele modulus, for which the classical numerical methods can present difficulties. For the convection-diffusion like case, the time-dependent partial differential equations are transformed by a semi-discrete Petrov-Galerkin finite element method into a system of ordinary differential equations of the initial-value type that can be readily solved. In the presence of a constant diffusivity, in slab geometry the convection-like terms are absent, and the combination of a fitted mesh finite difference method with a predictor-corrector method is used to solve the problem. Both the methods are found to converge, and general reaction rate forms can be treated. These methods are simple and highly efficient for arbitrary particle geometry and parameters, including a large Thiele modulus. (C) 2001 Elsevier Science Ltd. All rights reserved.

Fluorescence spectrum of a two-level atom driven by a multiple modulated field

We investigate the fluorescence spectrum of a two-level atom driven by a multiple amplitude-modulated field. The driving held is modeled as a polychromatic field composed of a strong central (resonant) component and a large number of symmetrically detuned sideband fields displaced from the central component by integer multiples of a constant detuning. Spectra obtained here differ qualitatively from those observed for a single pair of modulating fields [B. Blind, P.R. Fontana, and P. Thomann, J. Phys. B 13, 2717 (1980)]. In the case of a small number of the modulating fields, a multipeaked spectrum is obtained with the spectral features located at fixed frequencies that are independent of the number of modulating fields and their Rabi frequencies. As the number of the modulating fields increases, the spectrum ultimately evolves to the well-known Mellow triplet with the sidebands shifted from the central component by an effective Rabi frequency whose magnitude depends on the initial relative phases of the components of the driving held. For equal relative phases, the effective Rabi frequency of the driving field can be reduced to zero resulting in the disappearance of fluorescence spectrum, i.e., the atom can stop interacting with the field. When the central component and the modulating fields are 180 degrees out of phase, the spectrum retains its triplet structure with the sidebands located at frequencies equal to the sum of the Rabi frequencies of the component of the driving field. Moreover, we shaw that the frequency of spontaneous emission can be controlled and switched from one frequency to another when the Rabi frequency or initial phase of the modulating fields are varied.

Dynamical tunnelling of ultracold atoms

The divergence of quantum and classical descriptions of particle motion is clearly apparent in quantum tunnelling(1,2) between two regions of classically stable motion. An archetype of such nonclassical motion is tunnelling through an energy barrier. In the 1980s, a new process, 'dynamical' tunnelling(1-3), was predicted, involving no potential energy barrier; however, a constant of the motion (other than energy) still forbids classically the quantum-allowed motion. This process should occur, for example, in periodically driven, nonlinear hamiltonian systems with one degree of freedom(4-6). Such systems may be chaotic, consisting of regions in phase space of stable, regular motion embedded in a sea of chaos. Previous studies predicted(4) dynamical tunnelling between these stable regions. Here we observe dynamical tunnelling of ultracold atoms from a Bose-Einstein condensate in an amplitude-modulated optical standing wave. Atoms coherently tunnel back and forth between their initial state of oscillatory motion (corresponding to an island of regular motion) and the state oscillating 180 degrees out of phase with the initial state.

Tidal effects on groundwater dynamics in unconfined aquifers

The variation of seawater level resulting from tidal fluctuations is usually neglected in regional groundwater flow studies. Although the tidal oscillation is damped near the shoreline, there is a quasi-steady-slate rise in the mean water-table position, which may have an influence on regional groundwater flow. In this paper the effects of tidal fluctuations on groundwater hydraulics are investigated using a variably saturated numerical model that includes the effects of a realistic mild beach slope, seepage face and the unsaturated zone. In particular the impact of these factors on the velocity field in the aquifer is assessed. Simulations show that the tidal fluctuation has substantial consequences for the local velocity field in the vicinity of the exit face, which affects the nearshore migration of contaminant in coastal aquifers. An overheight in the water table as a result of the tidal fluctuation is observed anti this has a significant effect on groundwater discharge to the sea when the landward boundary condition is a constant water level. The effect of beach slope is very significant and simplifying the problem by considering a vertical beach face causes serious errors in predicting the water-table position and the groundwater flux. For media with a high effective capillary fringe, the moisture retained above the water table is important in determining the effects of the tidal fluctuations. Copyright (C) 2001 John Wiley & Sons, Ltd.