Distribution kinetics of solutes in the isolated in-situ perfused rat head using the multiple indicator dilution technique and a physiological two-barrier model

The purpose of this study was to determine the pharmacokinetics of [C-14]diclofenac, [C-14]salicylate and [H-3]clonidine using a single pass rat head perfusion preparation. The head was perfused with 3-[N-morpholino] propane-sulfonic acid-buffered Ringer's solution. Tc-99m-red blood cells and a drug were injected in a bolus into the internal carotid artery and collected from the posterior facial vein over 28 min. A two-barrier stochastic organ model was used to estimate the statistical moments of the solutes. Plasma, interstitial and cellular distribution volumes for the solutes ranged from 1.0 mL (diclofenac) to 1.6 mL (salicylate), 2.0 mL (diclofenac) to 4.2 mL (water) and 3.9 mL (salicylate) to 20.9 mL (diclofenac), respectively. A comparison of these volumes to water indicated some exclusion of the drugs from the interstitial space and salicylate from the cellular space. Permeability-surface area (PS) products calculated from plasma to interstitial fluid permeation clearances (CLPI) (range 0.02-0.40 mL s(-1)) and fractions of solute unbound in the perfusate were in the order: diclofenac>salicylate >clonidine>sucrose (from 41.8 to 0.10 mL s(-1)). The slow efflux of diclofenac, compared with clonidine and salicylate, may be related to its low average unbound fraction in the cells. This work accounts for the tail of disposition curves in describing pharmacokinetics in the head.

Cyclosporin A pretreatment in a rat model of warm ischaemia/reperfusion injury

Background/Aims: These studies investigated the role of apoptosis following ischaemia/reperfusion (I/R) injury to the liver and the effect of pretreatment with Cyclosporin A. Methods: Male Sprague-Dawley rats received 30 min of warm ischaemia followed by a period of reperfusion of 6 h. Rats were given olive oil or Cyclosporin A (30 mg/kg p.o.) the day before surgery. Neutrophil numbers were assessed in haematoxylin-eosin-stained sections of liver. In situ staining of sections using TdT-mediated dUTP-fluoreseein nick-end labelling was carried out to determine the extent of apoptosis, followed by electron microscopy. Semi-quantitative polymerase chain reaction (PCR) analysis of the transcript for Fas antigen was performed. Results and Conclusions: High levels of apoptosis were observed in I/R injury, which were greatly ameliorated in Cyclosporin A-pretreated groups. PCR analysis indicated a reduction in the level of expression of Fas transcript in Cyclosporin A-treated rats. Histological analysis showed a significant increase in the number of neutrophils infiltrating I/R-injured tissue (62 +/- 10.69, it = 16), which was markedly reduced by Cyclosporin A pretreatment (16 +/- 7, n = 6, P < 0.05). These results indicate a role of parenchymal apoptosis in the pathogenesis of I/R injury, which occurs in association with neutrophil infiltration, both of which can be significantly reduced by Cyclosporin A pretreatment. (C) 2002 European Association for the Study of the Liver. Published by Elsevier Science B.V. All rights reserved.

A new level of conotoxin diversity, a non-native disulfide bond connectivity in alpha-conotoxin AuIB reduces structural definition but increases biological activity

alpha-Conotoxin AuIB and a disulfide bond variant of AuIB have been synthesized to determine the role of disulfide bond connectivity on structure and activity. Both of these peptides contain the 15 amino acid sequence GCCSYPPCFATNPDC, with the globular (native) isomer having the disulfide connectivity Cys(2-8 and 3-15) and the ribbon isomer having the disulfide connectivity Cys(2-15 and 3-8). The solution structures of the peptides were determined by NAIR spectroscopy, and their ability to block the nicotinic acetylcholine receptors on dissociated neurons of the rat parasympathetic ganglia was examined. The ribbon disulfide isomer, although having a less well defined structure, is surprisingly found to have approximately 10 times greater potency than the native peptide. To our knowledge this is the first demonstration of a non-native disulfide bond isomer of a conotoxin exhibiting greater biological activity than the native isomer.

Simulation of the load-unload response ratio and critical sensitivity in the lattice solid model

The Load-Unload Response Ratio (LURR) method is an intermediate-term earthquake prediction approach that has shown considerable promise. It involves calculating the ratio of a specified energy release measure during loading and unloading where loading and unloading periods are determined from the earth tide induced perturbations in the Coulomb Failure Stress on optimally oriented faults. In the lead-up to large earthquakes, high LURR values are frequently observed a few months or years prior to the event. These signals may have a similar origin to the observed accelerating seismic moment release (AMR) prior to many large earthquakes or may be due to critical sensitivity of the crust when a large earthquake is imminent. As a first step towards studying the underlying physical mechanism for the LURR observations, numerical studies are conducted using the particle based lattice solid model (LSM) to determine whether LURR observations can be reproduced. The model is initialized as a heterogeneous 2-D block made up of random-sized particles bonded by elastic-brittle links. The system is subjected to uniaxial compression from rigid driving plates on the upper and lower edges of the model. Experiments are conducted using both strain and stress control to load the plates. A sinusoidal stress perturbation is added to the gradual compressional loading to simulate loading and unloading cycles and LURR is calculated. The results reproduce signals similar to those observed in earthquake prediction practice with a high LURR value followed by a sudden drop prior to macroscopic failure of the sample. The results suggest that LURR provides a good predictor for catastrophic failure in elastic-brittle systems and motivate further research to study the underlying physical mechanisms and statistical properties of high LURR values. The results provide encouragement for earthquake prediction research and the use of advanced simulation models to probe the physics of earthquakes.

Simulation of the influence of rate- and state-dependent friction on the macroscopic behavior of complex fault zones with the lattice solid model

In order to understand the earthquake nucleation process, we need to understand the effective frictional behavior of faults with complex geometry and fault gouge zones. One important aspect of this is the interaction between the friction law governing the behavior of the fault on the microscopic level and the resulting macroscopic behavior of the fault zone. Numerical simulations offer a possibility to investigate the behavior of faults on many different scales and thus provide a means to gain insight into fault zone dynamics on scales which are not accessible to laboratory experiments. Numerical experiments have been performed to investigate the influence of the geometric configuration of faults with a rate- and state-dependent friction at the particle contacts on the effective frictional behavior of these faults. The numerical experiments are designed to be similar to laboratory experiments by DIETERICH and KILGORE (1994) in which a slide-hold-slide cycle was performed between two blocks of material and the resulting peak friction was plotted vs. holding time. Simulations with a flat fault without a fault gouge have been performed to verify the implementation. These have shown close agreement with comparable laboratory experiments. The simulations performed with a fault containing fault gouge have demonstrated a strong dependence of the critical slip distance D-c on the roughness of the fault surfaces and are in qualitative agreement with laboratory experiments.

Genes induced by growth hormone in a model of adipogenic differentiation

A substantial number of GH regulated genes have been reported in mature hepatocytes. but genes involved in GH-initiated cell differentiation have not yet been identified. Here we have studied a, ell-characterised model of GH-dependent differentiation, adipogenesis of 3T3-F442A preadipocytes, to identify genes rapidly induced by GH. Using the suppression subtractive hybridisation technique, we have identified eight genes induced within 60 min of GH treatment, and verified these by northern analysis. Six were identifiable as Stat 2. Stat 3, thrombospondin-1. oncostatin M receptor beta chain. a DEAD box RNA helicase. and muscleblind. a developmental transcription factor. Bioinformatic approaches assigned one of the two remaining unknown genes as a novel 436 residue serine,threonine kinase. As each of the identified genes hake important developmental roles. they may be important in initiating GH-induced adipogenesis. (C) 2002 Elsevier Science Ireland Ltd. All rights reserved.

Role of PGF(2 alpha) and oxytocin in parturition in the brushtail possum (Trichosurus vulpecula)

Maturation of the fetal pituitary and adrenal glands allows the secretion of cortisol, which in turn leads to an increase in prostaglandin and mesotocin production. The production of prostaglandin and mesotocin results in an increase in uterine contractions and initiates birth in marsupials. The major metabolite of PGF(2alpha), 13,14-dihydro-15-keto-prostaglandin F-2alpha (PGFM), has been found in the plasma of the possum at the time of birth and administration of PGF(2alpha) to female possums induced the adoption of the birth position. Evidence that mesotocin is an integral hormone of birth in the tammar wallaby indicates that both PGF(2alpha) and mesotocin or oxytocin are required for marsupial birth. The presence of PGF(2alpha) receptors in the uterus and corpus luteum of the possum, and the in vitro uterine responsiveness to PGF(2alpha) or oxytocin, were examined. PGF(2alpha) receptors were not observed in possum uteri and the inability of PGF(2alpha) to cause contractions indicates that PGF(2alpha) is not involved directly in contraction of the uterus at parturition. The presence of oxytocin and mesotocin receptors in the uterus of possoms and the ability of oxytocin to induce uterine contraction in vitro supports the view that mesotocin is required for expulsion of the young from the uterus. Low numbers of PGF(2alpha) receptors were found in the possum corpus luteum at birth, indicating an involvement of PGF(2alpha) in regression of the corpus luteum.

An overview of APSIM, a model designed for farming systems simulation

The Agricultural Production Systems Simulator (APSIM) is a modular modelling framework that has been developed by the Agricultural Production Systems Research Unit in Australia. APSIM was developed to simulate biophysical process in farming systems, in particular where there is interest in the economic and ecological outcomes of management practice in the face of climatic risk. The paper outlines APSIM's structure and provides details of the concepts behind the different plant, soil and management modules. These modules include a diverse range of crops, pastures and trees, soil processes including water balance, N and P transformations, soil pH, erosion and a full range of management controls. Reports of APSIM testing in a diverse range of systems and environments are summarised. An example of model performance in a long-term cropping systems trial is provided. APSIM has been used in a broad range of applications, including support for on-farm decision making, farming systems design for production or resource management objectives, assessment of the value of seasonal climate forecasting, analysis of supply chain issues in agribusiness activities, development of waste management guidelines, risk assessment for government policy making and as a guide to research and education activity. An extensive citation list for these model testing and application studies is provided. Crown Copyright (C) 2002 Published by Elsevier Science B.V. All rights reserved.

A mathematical model of concentrate solids content for the wet drum magnetic separator

Low concentrate density from wet drum magnetic separators in dense medium circuits can cause operating difficulties due to inability to obtain the required circulating medium density and, indirectly, high medium solids losses. The literature is almost silent on the processes controlling concentrate density. However, the common name for the region through which concentrate is discharged-the squeeze pan gap-implies that some extrusion process is thought to be at work. There is no model of magnetics recovery in a wet drum magnetic separator, which includes as inputs all significant machine and operating variables. A series of trials, in both factorial experiments and in single variable experiments, was done using a purpose built rig which featured a small industrial scale (700 mm lip length, 900 turn diameter) wet drum magnetic separator. A substantial data set of 191 trials was generated in this work. The results of the factorial experiments were used to identify the variables having a significant effect on magnetics recovery. It is proposed, based both on the experimental observations of the present work and on observations reported in the literature, that the process controlling magnetic separator concentrate density is one of drainage. Such a process should be able to be defined by an initial moisture, a drainage rate and a drainage time, the latter being defined by the volumetric flowrate and the volume within the drainage zone. The magnetics can be characterised by an experimentally derived ultimate drainage moisture. A model based on these concepts and containing adjustable parameters was developed. This model was then fitted to a randomly chosen 80% of the data, and validated by application to the remaining 20%. The model is shown to be a good fit to data over concentrate solids content values from 40% solids to 80% solids and for both magnetite and ferrosilicon feeds. (C) 2003 Elsevier Science B.V. All rights reserved.

A mathematical model of recovery of dense medium magnetics in the wet drum magnetic separator

Loss of magnetic medium solids from dense medium circuits is a substantial contributor to operating cost. Much of this loss is by way of wet drum magnetic separator effluent. A model of the separator would be useful for process design, optimisation and control. A review of the literature established that although various rules of thumb exist, largely based on empirical or anecdotal evidence, there is no model of magnetics recovery in a wet drum magnetic separator which includes as inputs all significant machine and operating variables. A series of trials, in both factorial experiments and in single variable experiments, was therefore carried out using a purpose built rig which featured a small industrial scale (700 mm lip length, 900 mm diameter) wet drum magnetic separator. A substantial data set of 191 trials was generated in the work. The results of the factorial experiments were used to identify the variables having a significant effect on magnetics recovery. Observations carried out as an adjunct to this work, as well as magnetic theory, suggests that the capture of magnetic particles in the wet drum magnetic separator is by a flocculation process. Such a process should be defined by a flocculation rate and a flocculation time; the latter being defined by the volumetric flowrate and the volume within the separation zone. A model based on this concept and containing adjustable parameters was developed. This model was then fitted to a randomly chosen 80% of the data, and validated by application to the remaining 20%. The model is shown to provide a satisfactory fit to the data over three orders of magnitude of magnetics loss. (C) 2003 Elsevier Science BY. All rights reserved.

Hemodynamic and ventilatory effects of manual respiratory physiotherapy techniques of chest clapping, vibration, and shaking in an animal model

Chest clapping, vibration, and shaking were studied in 10 physiotherapists who applied these techniques on an anesthetized animal model. Hemodynamic variables (such as heart rate, blood pressure, pulmonary artery pressure, and right atrial pressure) were measured during the application of these techniques to verify claims of adverse events. In addition, expired tidal volume and peak expiratory flow rate were measured to ascertain effects of these techniques. Physiotherapists in this study applied chest clapping at a rate of 6.2 +/- 0.9 Hz, vibration at 10.5 +/- 2.3 Hz, and shaking at 6.2 +/- 2.3 Hz. With the use of these rates, esophageal pressure swings of 8.8 +/- 5.0, 0.7 +/- 0.3, and 1.4 +/- 0.7 mmHg resulted from clapping, vibration, and shaking respectively. Variability in rates and forces generated by these techniques was 80% of variance in shaking force (P = 0.003). Application of these techniques by physiotherapists was found to have no significant effects on hemodynamic and most ventilatory variables in this study. From this study, we conclude that chest clapping, vibration, and shaking 1) can be consistently performed by physiotherapists; 2) are significantly related to physiotherapists' characteristics, particularly clinical experience; and 3) caused no significant hemodynamic effects.

Impulse-response function of splanchnic circulation with model-independent constraints: theory and experimental validation

Modeling physiological processes using tracer kinetic methods requires knowledge of the time course of the tracer concentration in blood supplying the organ. For liver studies, however, inaccessibility of the portal vein makes direct measurement of the hepatic dual-input function impossible in humans. We want to develop a method to predict the portal venous time-activity curve from measurements of an arterial time-activity curve. An impulse-response function based on a continuous distribution of washout constants is developed and validated for the gut. Experiments with simultaneous blood sampling in aorta and portal vein were made in 13 anesthetized pigs following inhalation of intravascular [O-15] CO or injections of diffusible 3-O[ C-11] methylglucose (MG). The parameters of the impulse-response function have a physiological interpretation in terms of the distribution of washout constants and are mathematically equivalent to the mean transit time ( T) and standard deviation of transit times. The results include estimates of mean transit times from the aorta to the portal vein in pigs: (T) over bar = 0.35 +/- 0.05 min for CO and 1.7 +/- 0.1 min for MG. The prediction of the portal venous time-activity curve benefits from constraining the regression fits by parameters estimated independently. This is strong evidence for the physiological relevance of the impulse-response function, which includes asymptotically, and thereby justifies kinetically, a useful and simple power law. Similarity between our parameter estimates in pigs and parameter estimates in normal humans suggests that the proposed model can be adapted for use in humans.

An EM-based Semi-Parametric Mixture Model Approach to the Regression Analysis of Competing-Risks Data

We consider a mixture model approach to the regression analysis of competing-risks data. Attention is focused on inference concerning the effects of factors on both the probability of occurrence and the hazard rate conditional on each of the failure types. These two quantities are specified in the mixture model using the logistic model and the proportional hazards model, respectively. We propose a semi-parametric mixture method to estimate the logistic and regression coefficients jointly, whereby the component-baseline hazard functions are completely unspecified. Estimation is based on maximum likelihood on the basis of the full likelihood, implemented via an expectation-conditional maximization (ECM) algorithm. Simulation studies are performed to compare the performance of the proposed semi-parametric method with a fully parametric mixture approach. The results show that when the component-baseline hazard is monotonic increasing, the semi-parametric and fully parametric mixture approaches are comparable for mildly and moderately censored samples. When the component-baseline hazard is not monotonic increasing, the semi-parametric method consistently provides less biased estimates than a fully parametric approach and is comparable in efficiency in the estimation of the parameters for all levels of censoring. The methods are illustrated using a real data set of prostate cancer patients treated with different dosages of the drug diethylstilbestrol. Copyright (C) 2003 John Wiley Sons, Ltd.