Pharmacokinetic-pharmacodynamic modelling of QT interval prolongation following citalopram overdoses

Aims To develop a pharmacokinetic-pharmacodynamic model describing the time-course of QT interval prolongation after citalopram overdose and to evaluate the effect of charcoal on the relative risk of developing abnormal QT and heart-rate combinations. Methods Plasma concentrations and electrocardiograph (ECG) data from 52 patients after 62 citalopram overdose events were analysed in WinBUGS using a Bayesian approach. The reported doses ranged from 20 to 1700 mg and on 17 of the events a single dose of activated charcoal was administered. The developed pharmacokinetic-pharmacodynamic model was used for predicting the probability of having abnormal combinations of QT-RR, which was assumed to be related to an increased risk for torsade de pointes (TdP). Results The absolute QT interval was related to the observed heart rate with an estimated individual heart-rate correction factor [alpha = 0.36, between-subject coefficient of variation (CV) = 29%]. The heart-rate corrected QT interval was linearly dependent on the predicted citalopram concentration (slope = 40 ms l mg(-1), between-subject CV = 70%) in a hypothetical effect-compartment (half-life of effect-delay = 1.4 h). The heart-rate corrected QT was predicted to be higher in women than in men and to increase with age. Administration of activated charcoal resulted in a pronounced reduction of the QT prolongation and was shown to reduce the risk of having abnormal combinations of QT-RR by approximately 60% for citalopram doses above 600 mg. Conclusion Citalopram caused a delayed lengthening of the QT interval. Administration of activated charcoal was shown to reduce the risk that the QT interval exceeds a previously defined threshold and therefore is expected to reduce the risk of TdP.

Modelling of two-stage anaerobic digestion using the IWA Anaerobic Digestion Model No. 1 (ADM1)

The aim of the study presented was to implement a process model to simulate the dynamic behaviour of a pilot-scale process for anaerobic two-stage digestion of sewage sludge. The model implemented was initiated to support experimental investigations of the anaerobic two-stage digestion process. The model concept implemented in the simulation software package MATLAB(TM)/Simulink(R) is a derivative of the IWA Anaerobic Digestion Model No.1 (ADM1) that has been developed by the IWA task group for mathematical modelling of anaerobic processes. In the present study the original model concept has been adapted and applied to replicate a two-stage digestion process. Testing procedures, including balance checks and 'benchmarking' tests were carried out to verify the accuracy of the implementation. These combined measures ensured a faultless model implementation without numerical inconsistencies. Parameters for both, the thermophilic and the mesophilic process stage, have been estimated successfully using data from lab-scale experiments described in literature. Due to the high number of parameters in the structured model, it was necessary to develop a customised procedure that limited the range of parameters to be estimated. The accuracy of the optimised parameter sets has been assessed against experimental data from pilot-scale experiments. Under these conditions, the model predicted reasonably well the dynamic behaviour of a two-stage digestion process in pilot scale. (C) 2004 Elsevier Ltd. All rights reserved.

Matching conditional moments in PDF modelling of nonpremixed combustion

The rate of generation of fluctuations with respect to the scalar values conditioned on the mixture fraction, which significantly affects turbulent nonpremixed combustion processes, is examined. Simulation of the rate in a major mixing model is investigated and the derived equations can assist in selecting the model parameters so that the level of conditional fluctuations is better reproduced by the models. A more general formulation of the multiple mapping conditioning (MMC) model that distinguishes the reference and conditioning variables is suggested. This formulation can be viewed as a methodology of enforcing certain desired conditional properties onto conventional mixing models. Examples of constructing consistent MMC models with dissipation and velocity conditioning as well as of combining MMC with large eddy simulations (LES) are also provided. (c) 2005 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

Mechanistic modelling of pilling. Part II: Individual-fibre computational model

A one-dimensional computational model of pilling of a fibre assembly has been created. The model follows a set of individual fibres, as free ends and loops appear as fuzz and arc progressively withdrawn from the body of the assembly, and entangle to form pills, which eventually break off or are pulled out. The time dependence of the computation is given by ticks, which correspond to cycles of a wear and laundering process. The movement of the fibres is treated as a reptation process. A set of standard values is used as inputs to the computation. Predictions arc given of the change with a number Of cycles of mass of fuzz, mass of pills, and mass removed from the assembly. Changes in the standard values allow sensitivity studies to be carried out.

Mechanistic modelling of pilling. Part I: Detailing of mechanisms

Mechanistic models of pilling are discussed in general terms, and a framework for pilling simulations is thereby created. A fundamental flaw in earlier models of pilling is revealed. A more comprehensive model of fibre diffusion and withdrawal from the fabric is proposed, and this is solved in general terms to find the rate of fuzz growth. Fuzz wear-off and entanglement into pills are discussed. Fibre fatigue is introduced, and it is demonstrated that this potentially increases the rate of withdrawal of anchor fibres.

Exploratory modelling of grinding pressure within a compressed particle bed

With increasing industry interest in high pressure roll grinding (HPGR) technology, there is a strong incentive for improved understanding of the nature of grinding pressure that exists in the interior of a compressed particle bed. This corresponds to the crushing region of the HPGR. The relationship between applied pressure (stress) to the particle bed and induced pressure (stress) within particles and at contact points between particles is of particular interest. A detailed parametric investigation is beyond the scope of this exploratory paper. However, this exploratory investigation does suggest some interesting behaviour. The compressed particle bed within an 80 turn diameter piston has been modelled using Particle Flow Code for three dimensions. PFC3D is a discrete element code. The total number of simulated particles was 1225 and 2450 for two beds of different thickness. Particle diameters were uniformly distributed between 4 and 4.5 mm. The results of the simulations show that stress intensity within the simulated particle beds and within the observed particles increased with increase of the applied stress. The intensity of the average vertical stress in the selected particles tended to be comparable with the intensity of the pressure applied to the surface of particle bed and was only occasionally higher. However, the stress at contact points between particles could be several times higher. In a real crusher, such high stress amplification at contacts will quickly decrease due to local crushing and a resultant increase the size of the contact area. Therefore, its significance is likely to be relatively small in an industrial context. The modelling results also suggest that failure within the particle bed will progress from the crushing surface towards the depth of the bed. (c) 2006 Elsevier Ltd. All rights reserved.

Modelling of entrainment in industrial flotation cells: Water recovery and degree of entrainment

Entrainment in flotation can be considered as a two-step process, including the transfer of the suspended solids in the top of the pulp region just below the pulp-froth interface to the froth phase and the transfer of the entrained particles in the froth phase to the concentrate. Both steps have a strong classification characteristic. The degree of entrainment describes the classification effect of the drainage process in the froth phase. This paper briefly reviews two existing models of degree of entrainment. Experimental data were collected from an Outokumpu 3 m(3) tank cell in the Xstrata Mt. Isa Mines copper concentrator. The data are fitted to the models and the effect of cell operating conditions including air rate and froth height on the degree of entrainment is examined on a size-by-size basis. It is found that there is a strong correlation between the entrainment and the water recovery, which is close to lineal. for the fines. The degree of entrainment decreases with increase in particle size. Within the normal range of cell operating conditions, few particles coarser than 50 mu m are recovered by entrainment. In general, the degree of entrainment increases with increase in the ail rate and decreases with increase in the froth height. Air rate and froth height strongly interact with each other and affect the entrainment process mainly via changes in the froth retention time, the froth structure and froth properties. As a result, other mechanisms such as entrapment may become important in recovering the coarse entrained particles. (c) 2005 Elsevier Ltd. All rights reserved.

Multi-level zero-inflated Poisson regression modelling of correlated count data with excess zeros

Count data with excess zeros relative to a Poisson distribution are common in many biomedical applications. A popular approach to the analysis of such data is to use a zero-inflated Poisson (ZIP) regression model. Often, because of the hierarchical Study design or the data collection procedure, zero-inflation and lack of independence may occur simultaneously, which tender the standard ZIP model inadequate. To account for the preponderance of zero counts and the inherent correlation of observations, a class of multi-level ZIP regression model with random effects is presented. Model fitting is facilitated using an expectation-maximization algorithm, whereas variance components are estimated via residual maximum likelihood estimating equations. A score test for zero-inflation is also presented. The multi-level ZIP model is then generalized to cope with a more complex correlation structure. Application to the analysis of correlated count data from a longitudinal infant feeding study illustrates the usefulness of the approach.

Network topology and the evolution of dynamics in an artificial genetic regulatory network model created by whole genome duplication and divergence

Topological measures of large-scale complex networks are applied to a specific artificial regulatory network model created through a whole genome duplication and divergence mechanism. This class of networks share topological features with natural transcriptional regulatory networks. Specifically, these networks display scale-free and small-world topology and possess subgraph distributions similar to those of natural networks. Thus, the topologies inherent in natural networks may be in part due to their method of creation rather than being exclusively shaped by subsequent evolution under selection. The evolvability of the dynamics of these networks is also examined by evolving networks in simulation to obtain three simple types of output dynamics. The networks obtained from this process show a wide variety of topologies and numbers of genes indicating that it is relatively easy to evolve these classes of dynamics in this model. (c) 2006 Elsevier Ireland Ltd. All rights reserved.

On the modelling of wave breaking and set-up on coral reefs

An extended refraction-diffraction equation [Massel, S.R., 1993. Extended refraction-diffraction equation for surface waves. Coastal Eng. 19, 97-126] has been applied to predict wave transformation and breaking as well as wave-induced set-up on two-dimensional reef profiles of various shapes. A free empirical coefficient alpha in a formula for the average rate of energy dissipation [epsilon(b)] = (alpha rho g omega/8 pi)(root gh/C)(H-3/h) in the modified periodic bore model was found to be a function of the dimensionless parameter F-c0 = (g(1.25)H(0)(0.5)T(2.5))/h(r)(1.75), proposed by Gourlay [Gourlayl M.R., 1994. Wave transformation on a coral reef. Coastal Eng. 23, 17-42]. The applicability of the developed model has been demonstrated for reefs of various shapes subjected to various incident wave conditions. Assuming proposed relationships of the coefficient alpha and F-c0, the model provides results on wave height attenuation and set-up elevation which compare well with experimental data. (C) 2000 Elsevier Science B.V. All rights reserved.

Stochastic modelling of the invasion process of nematodes in fly larvae

Stochastic models based on Markov birth processes are constructed to describe the process of invasion of a fly larva by entomopathogenic nematodes. Various forms for the birth (invasion) rates are proposed. These models are then fitted to data sets describing the observed numbers of nematodes that have invaded a fly larval after a fixed period of time. Non-linear birthrates are required to achieve good fits to these data, with their precise form leading to different patterns of invasion being identified for three populations of nematodes considered. One of these (Nemasys) showed the greatest propensity for invasion. This form of modelling may be useful more generally for analysing data that show variation which is different from that expected from a binomial distribution.