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.

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.

Magnetization transfer imaging for polymer gel dosimetry

Off-resonance RF pre-saturation was used to obtain contrast in MRI images of polymer gel dosimeters irradiated to doses up to 50 Gy. Two different polymer gel dosimeters composed of 2-hydroxyethyl-acryl ate or methacrylic acid monomers mixed with N, N'-methylene-bisacrylamide (BIS), dispersed in an aqueous gelatin matrix were evaluated. Radiation-induced polymerization of the co-monomers generates a fast-relaxing insoluble polymer. Saturation of the polymer using off-resonance Gaussian RF pulses prior to a spin-echo read-out with a short echo time leads to contrast that is dependent on the absorbed dose. This contrast is attributed to magnetization transfer (MT) between free water and the polymer, and direct saturation of water was found to be negligible under the prevailing experimental conditions. The usefulness of MT imaging was assessed by computing the dose resolution obtained with this technique. We found a low value of dose resolution over a wide range of doses could be obtained with a single experiment. This is an advantage over multiple spin echo (MSE) experiments using a single echo spacing where an optimal dose resolution is achieved over only very limited ranges of doses. The results suggest MT imaging protocols may be developed into a useful tool for polymer gel dosimetry.

Formal methods within a totally functional approach to programming

Taking functional programming to its extremities in search of simplicity still requires integration with other development (e.g. formal) methods. Induction is the key to deriving and verifying functional programs, but can be simplified through packaging proofs with functions, particularly folds, on data (structures). Totally Functional Programming avoids the complexities of interpretation by directly representing data (structures) as platonic combinators - the functions characteristic to the data. The link between the two simplifications is that platonic combinators are a kind of partially-applied fold, which means that platonic combinators inherit fold-theoretic properties, but with some apparent simplifications due to the platonic combinator representation. However, despite observable behaviour within functional programming that suggests that TFP is widely-applicable, significant work remains before TFP as such could be widely adopted.

Determination of the force constant of a single-beam gradient trap by measurement of backscattered light

A single-beam gradient trap could potentially be used to hold a stylus for scanning force microscopy. With a view to development of this technique, we modeled the optical trap as a harmonic oscillator and therefore characterized it by its force constant. We measured force constants and resonant frequencies for 1-4-mu m-diameter polystyrene spheres in a single-beam gradient trap using measurements of back-scattered light. Force constants were determined with both Gaussian and doughnut laser modes, with powers of 3 and 1 mW, respectively. Typical values for spring constants were measured to be between 10(-6) and 4 x 10(-6) N/m. The resonant frequencies of trapped particles were measured to be between 1 and 10 kHz, and the rms amplitudes of oscillations were estimated to be around 40 nm. Our results confirm that the use of the doughnut mode for single-beam trapping is more efficient in the axial direction. (C) 1996 Optical Society of America.

On the validity of the dispersion model of hepatic drug elimination when intravascular transit time densities are long-tailed

The dispersion model with mixed boundary conditions uses a single parameter, the dispersion number, to describe the hepatic elimination of xenobiotics and endogenous substances. An implicit a priori assumption of the model is that the transit time density of intravascular indicators is approximated by an inverse Gaussian distribution. This approximation is limited in that the model poorly describes the tail part of the hepatic outflow curves of vascular indicators. A sum of two inverse Gaussian functions is proposed as ail alternative, more flexible empirical model for transit time densities of vascular references. This model suggests that a more accurate description of the tail portion of vascular reference curves yields an elimination rate constant (or intrinsic clearance) which is 40% less than predicted by the dispersion model with mixed boundary conditions. The results emphasize the need to accurately describe outflow curves in using them as a basis for determining pharmacokinetic parameters using hepatic elimination models. (C) 1997 Society for Mathematical Biology.

Relative dispersion of intravascular transit times during isolated human limb perfusions for recurrent melanoma

Aims We have characterized the relative dispersion of vascular and extravascular markers in the limbs of three patients undergoing isolated limb perfusions with the cytotoxic melphalan for recurrent malignant melanoma both before and after melphalan dosing. Methods A bolus of injectate containing [Cr-51] labelled red blood cells, [C-14]-sucrose and [H-3]-water was injected into an iliac or femoral artery and outflow samples collected at 1 s intervals by a fraction collector. The radioactivity due to each isotype was analysed by either gamma [Cr-51] or beta [C-14 and H-3] counting. The moments of the outflow fraction-time profiles were estimated by a nonparametric (numerical integration) method and a parametric model (sum of two inverse Gaussian functions). Results The availability, mean transit time and normalised variance (CV2) obtained for labelled red blood cells, sucrose and water were similar before and after melphalan dosing and with the two methods of calculation but varied between the patients. Conclusions The vascular space is not well-stirred but characterized by a CV2 similar that reported previously for in situ rat hind limb and rat liver perfusions. A flow-limited blood-tissue exchange was observed for the permeating indicators. Administration of melphalan did not influence the distribution characteristics of the indicators.

Topological charge and angular momentum of light beams carrying optical vortices

We analyze the properties of light beams carrying phase singularities, or optical vortices. The transformations of topological charge during free-space propagation of a light wave, which is a combination of a Gaussian beam and a multiple charged optical vortex within a Gaussian envelope, are studied both in theory and experiment. We revise the existing knowledge about topological charge conservation, and demonstrate possible scenarios where additional vortices appear or annihilate during free propagation of such a combined beam. Coaxial interference of optical vortices is also analyzed, and the general rule for angular-momentum density distribution in a combined beam is established. We show that, in spite of any variation in the number of vortices in a combined beam, the total angular momentum is constant during the propagation. [S1050-2947(97)09910-1].

Anisotropy-based performance analysis of linear discrete time invariant control systems

The anisotropic norm of a linear discrete-time-invariant system measures system output sensitivity to stationary Gaussian input disturbances of bounded mean anisotropy. Mean anisotropy characterizes the degree of predictability (or colouredness) and spatial non-roundness of the noise. The anisotropic norm falls between the H-2 and H-infinity norms and accommodates their loss of performance when the probability structure of input disturbances is not exactly known. This paper develops a method for numerical computation of the anisotropic norm which involves linked Riccati and Lyapunov equations and an associated special type equation.

Elastic moduli of model random three-dimensional closed-cell cellular solids

Most cellular solids are random materials, while practically all theoretical structure-property results are for periodic models. To be able to generate theoretical results for random models, the finite element method (FEM) was used to study the elastic properties of solids with a closed-cell cellular structure. We have computed the density (rho) and microstructure dependence of the Young's modulus (E) and Poisson's ratio (PR) for several different isotropic random models based on Voronoi tessellations and level-cut Gaussian random fields. The effect of partially open cells is also considered. The results, which are best described by a power law E infinity rho (n) (1<n<2), show the influence of randomness and isotropy on the properties of closed-cell cellular materials, and are found to be in good agreement with experimental data. (C) 2001 Acta Materialia Inc. Published by Elsevier Science Ltd. All rights reserved.

Considerations in the design of a mixed-method cluster evaluation of a community programme for 'at-risk' young people

This article discusses the design of a comprehensive evaluation of a community development programme for young people 'at-risk' of self-harming behaviour. It outlines considerations in the design of the evaluation and focuses on the complexities and difficulties associated with the evaluation of a community development programme. The challenge was to fulfil the needs of the funding body for a broad, outcome-focused evaluation while remaining close enough to the programme to accurately represent its activities and potential effects at a community level. Specifically, the strengths and limitations of a mixed-method evaluation plan are discussed with recommendations for future evaluation practice.

Caveolin and Ras function

Experimental protocols that allow confident assignment of signaling proteins to specific subdomains of the plasma membrane are essential for a full understanding of the complexities of signal transduction. This is especially relevant for Ras proteins, where the different membrane anchors of the Ras isoforms target them to functionally distinct microdomains that in turn allow quantitatively different signal outputs from otherwise highly homologous proteins. The methods outlined in this chapter, in addition to being invaluable in addressing Ras function, should also have wide utility in the study of many mammalian signal transduction pathways.

Interlaminar differences in the pyramidal cell phenotype in cortical areas 7m and STP (the superior temporal polysensory area) of the macaque monkey

Pyramidal neurones were injected with Lucifer Yellow in slices cut tangential to the surface of area 7m and the superior temporal polysensory area (STP) of the macaque monkey. Comparison of the basal dendritic arbors of supra- and infragranular pyramidal neurones (n=139) that were injected in the same putative modules in the different cortical areas revealed variation in their structure. Moreover, there were relative differences in dendritic morphology of supra- and infragranular pyramidal neurones in the two cortical areas. Shell analyses revealed that layer III pyramidal neurones in area STP had considerably higher peak complexity (maximum number of dendritic intersections per Shell circle) than those in layer V, whereas peak complexities were similar for supra- and infragranular pyramidal neurones in area 7m. In both cortical areas, the basal dendritic trees of layer m pyramidal neurones were characterized by a higher spine density than those in layer V. Calculations of the total number of dendritic spines in the average basal dendritic arbor revealed that layer V pyramidal neurones in area 7m had twice as many spines as cells in layer III. (4535 and 2294, respectively). A similar calculation for neurones in area STP revealed that layer III pyramidal neurones had approximately the same number of spines as cells in layer V (3585 and 3850 spines, respectively). Relative differences in the branching patterns of, and the number of spines in, the basal dendritic arbors of supra- and infragranular pyramidal neurones in the different cortical areas may allow for integration of different numbers of inputs, and different degrees of dendritic processing. These results support the thesis that intra-areal circuitry differs in different cortical areas.

Derivation of the probability distribution function for the local density of states of a disordered quantum wire via the replica trick and supersymmetry

We consider the statistical properties of the local density of states of a one-dimensional Dirac equation in the presence of various types of disorder with Gaussian white-noise distribution. It is shown how either the replica trick or supersymmetry can be used to calculate exactly all the moments of the local density of states.' Careful attention is paid to how the results change if the local density of states is averaged over atomic length scales. For both the replica trick and supersymmetry the problem is reduced to finding the ground state of a zero-dimensional Hamiltonian which is written solely in terms of a pair of coupled spins which are elements of u(1, 1). This ground state is explicitly found for the particular case of the Dirac equation corresponding to an infinite metallic quantum wire with a single conduction channel. The calculated moments of the local density of states agree with those found previously by Al'tshuler and Prigodin [Sov. Phys. JETP 68 (1989) 198] using a technique based on recursion relations for Feynman diagrams. (C) 2001 Elsevier Science B.V. All rights reserved.

SAR imaging of buried objects from MoM modelled scattered field

Synthetic aperture radar (SAR) images of resonant buried objects are modelled in the presence of ground surface clutter. The method of moments (MoM) is used to model scattered fields from a resonant buried conductor and clutter is modelled as a bivariant Gaussian distribution. A diffraction stack SAR imaging technique is applied to the ultra-wideband waveforms to give a bipolar signal image. A number of examples have been computed to illustrate the combined effects of SAR processing with resonant targets and clutter. SAR images of different targets show differences which may facilitate target identification. To maximise the peak signal-to-clutter ratio, an image correlation technique is applied and the results are shown.