Stiffly accurate Runge-Kutta methods for stiff stochastic differential equations

In this paper we discuss implicit methods based on stiffly accurate Runge-Kutta methods and splitting techniques for solving Stratonovich stochastic differential equations (SDEs). Two splitting techniques: the balanced splitting technique and the deterministic splitting technique, are used in this paper. We construct a two-stage implicit Runge-Kutta method with strong order 1.0 which is corrected twice and no update is needed. The stability properties and numerical results show that this approach is suitable for solving stiff SDEs. (C) 2001 Elsevier Science B.V. All rights reserved.

A retrospective analysis of mycophenolic acid and cyclosporin concentrations with acute rejection in renal transplant recipients

Objectives: Although monitoring of cyclosporin (CsA) is standard clinical practice postrenal transplantation. mycophenolic acid (MPA) concentrations are not routinely measured. There is evidence that a relationship exists between MPA area under the concentration-time curve (AUC) and rejection. In this study, a retrospective analysis was undertaken of 27 adult renal transplant recipients. Methods: Patients received CsA and MPA therapy and had a four-point MPA AUC investigation. The relationship between MPA AUC performed in the first week after transplantation, as well as median trough cyclosporin concentrations, and clinical outcomes in the first month posttransplant were evaluated. Results: A total of 12 patients experienced biopsy proven rejection (44.4%) and 4 patients had gastrointestinal adverse events (14.8%). A statistically significant relationship was observed between the incidence of biopsy proven rejection and both MPA AUC (p = 0.02) and median trough CsA concentration (p = 0.008). No relationship between trough MPA concentration and rejection was observed (p = 0.21). Only 3 of 11 (27%) patients with an MPA AUC > 30 mg.h/L and a median trough CsA > 175 mug/L experienced acute rejection, compared with a 56% incidence of rejection for the remaining 16 patients. Patients who experienced adverse gastrointestinal events had significantly lower MPA AUC (p = 0.04), but median trough CsA concentrations were not significantly different (p = 0.24). Further, 3 of these 4 patients had rejection episodes. Conclusions: in addition to standard CsA monitoring, we propose further investigation of the use of a 4-point sampling strategy to predict MPA AUC in the first week posttransplant, which may facilitate optimization of mycophenolate mofetil dose at a rime when patients are most vulnerable to acute rejection. (C) 2001 The Canadian Society of Clinical Chemists. All rights reserved.

Crystallization and preliminary X-ray diffraction studies of FHA domains of Dun1 and Rad53 protein kinases

Forkhead-associated (FHA) domains are modular protein–protein interaction domains of ~130 amino acids present in numerous signalling proteins. FHA-domain-dependent protein interactions are regulated by phosphorylation of target proteins and FHA domains may be multifunctional phosphopeptide-recognition modules. FHA domains of the budding yeast cell-cycle checkpoint protein kinases Dun1p and Rad53p have been crystallized. Crystals of the Dun1-FHA domain exhibit the symmetry of the space group P6122 or P6522, with unit-cell parameters a = b = 127.3, c = 386.3 Å; diffraction data have been collected to 3.1 Å resolution on a synchrotron source. Crystals of the N-terminal FHA domain (FHA1) of Rad53p diffract to 4.0 Å resolution on a laboratory X-ray source and have Laue-group symmetry 4/mmm, with unit-cell parameters a = b = 61.7, c = 104.3 Å.

Model based procedure for scale-up of wet overflow ball mills Part 1: Outline of the methodology

Bond's method for ball mill scale-up only gives the mill power draw for a given duty. This method is incompatible with computer modelling and simulation techniques. It might not be applicable for the design of fine grinding ball mills and ball mills preceded by autogenous and semi-autogenous grinding mills. Model-based ball mill scale-up methods have not been validated using a wide range of full-scale circuit data. Their accuracy is therefore questionable. Some of these methods also need expensive pilot testing. A new ball mill scale-up procedure is developed which does not have these limitations. This procedure uses data from two laboratory tests to determine the parameters of a ball mill model. A set of scale-up criteria then scales-up these parameters. The procedure uses the scaled-up parameters to simulate the steady state performance of full-scale mill circuits. At the end of the simulation, the scale-up procedure gives the size distribution, the volumetric flowrate and the mass flowrate of all the streams in the circuit, and the mill power draw.

Time evolution in the unimolecular master equation at low temperatures: full spectral solution with scalable iterative methods and high precision

A new method is presented to determine an accurate eigendecomposition of difficult low temperature unimolecular master equation problems. Based on a generalisation of the Nesbet method, the new method is capable of achieving complete spectral resolution of the master equation matrix with relative accuracy in the eigenvectors. The method is applied to a test case of the decomposition of ethane at 300 K from a microcanonical initial population with energy transfer modelled by both Ergodic Collision Theory and the exponential-down model. The fact that quadruple precision (16-byte) arithmetic is required irrespective of the eigensolution method used is demonstrated. (C) 2001 Elsevier Science B.V. All rights reserved.

Analytical validation for a series of marker compounds used to assess renal drug elimination processes

Renal drug elimination is determined by glomerular filtration, tubular secretion, and tubular reabsorption. Changes in the integrity of these processes influence renal drug clearance, and these changes may not be detected by conventional measures of renal function such as creatinine clearance. The aim of the current study was to examine the analytic issues needed to develop a cocktail of marker drugs (fluconazole, rac-pindolol, para-aminohippuric acid, sinistrin) to measure simultaneously the mechanisms contributing to renal clearance. High-performance liquid chromatographic methods of analysis for fluconazole, pindolol, para-aminohippuric acid, and creatinine and an enzymatic assay for sinistrin were developed or modified and then validated to allow determination of each of the compounds in both plasma and urine in the presence of all other marker drugs. A pilot clinical study in one volunteer was conducted to ensure that the assays were suitable for quantitating all the marker drugs to the sensitivity and specificity needed to allow accurate determination of individual renal clearances. The performance of all assays (plasma and urine) complied with published validation criteria. All standard curves displayed linearity over the concentration ranges required, with coefficients of correlation greater than 0.99. The precision of the interday and intraday variabilities of quality controls for each marker in plasma and urine were all less than 11.9% for each marker. Recoveries of markers (and internal standards) in plasma and urine were all at least 90%. All markers investigated were shown to be stable when plasma or urine was frozen and thawed. For all the assays developed, there were no interferences from other markers or endogenous substances. In a pilot clinical study, concentrations of all markers could be accurately and reproducibly determined for a sufficient duration of time after administration to calculate accurate renal clearance for each marker. This article presents details of the analytic techniques developed for measuring concentrations of marker drugs for different renal elimination processes administered as a single dose to define the processes contributing to renal drug elimination.

Outcomes of cephalosporin treatment for serious infection due to apparently susceptible organisms producing extended-spectrum b-Lactamases: Implications for the clinical microbiology laboratory

Although extended-spectrum beta-lactamases (ESBLs) hydrolyze cephalosporin antibiotics, some ESBL-producing organisms are not resistant to all cephalosporins when tested in vitro. Some authors have suggested that screening klebsiellae or Escherichia coli for ESBL production is not clinically necessary, and when most recently surveyed the majority of American clinical microbiology laboratories did not make efforts to detect ESBLs, We performed a prospective, multinational study of Klebsiella pneumoniae bacteremia and identified 10 patients who were treated for ESBL-producing K. pneumoniae bacteremia with cephalosporins and whose infecting organisms were not resistant in vitro to the utilized cephalosporin. In addition, we reviewed 26 similar cases of severe infections which had previously been reported. Of these 36 patients, 4 had to be excluded from analysis. Of the remaining 32 patients, 100% (4 of 4) patients experienced clinical failure when MICs of the cephalosporin used for treatment were in the intermediate range and 54% (15 of 28) experienced failure when MICs of the cephalosporin used for treatment were in the susceptible range, Thus, it is clinically important to detect ESBL production by klebsiellae or E, coli even when cephalosporin MICs are in the susceptible range (less than or equal to 8 mug/ml) and to report ESBL-producing organisms as resistant to aztreonam and all cephalosporins (with the exception of cephamycins).

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.

On surrogate methods for detecting lateral gene transfer

Surrogate methods for detecting lateral gene transfer are those that do not require inference of phylogenetic trees. Herein I apply four such methods to identify open reading frames (ORFs) in the genome of Escherichia coli K12 that may have arisen by lateral gene transfer. Only two of these methods detect the same ORFs more frequently than expected by chance, whereas several intersections contain many fewer ORFs than expected. Each of the four methods detects a different non-random set of ORFs. The methods may detect lateral ORFs of different relative ages; testing this hypothesis will require rigorous inference of trees. (C) 2001 Federation of European Microbiological Societies. Published by Elsevier Science BN. All rights reserved.

Probabilistic modelling of CO2 corrosion laboratory data using neural networks

The present paper addresses two major concerns that were identified when developing neural network based prediction models and which can limit their wider applicability in the industry. The first problem is that it appears neural network models are not readily available to a corrosion engineer. Therefore the first part of this paper describes a neural network model of CO2 corrosion which was created using a standard commercial software package and simple modelling strategies. It was found that such a model was able to capture practically all of the trends noticed in the experimental data with acceptable accuracy. This exercise has proven that a corrosion engineer could readily develop a neural network model such as the one described below for any problem at hand, given that sufficient experimental data exist. This applies even in the cases when the understanding of the underlying processes is poor. The second problem arises from cases when all the required inputs for a model are not known or can be estimated with a limited degree of accuracy. It seems advantageous to have models that can take as input a range rather than a single value. One such model, based on the so-called Monte Carlo approach, is presented. A number of comparisons are shown which have illustrated how a corrosion engineer might use this approach to rapidly test the sensitivity of a model to the uncertainities associated with the input parameters. (C) 2001 Elsevier Science Ltd. All rights reserved.

Effect of hemochromatosis genotype and lifestyle factors on iron and red cell indices in a community population

Background: Heterozygotes for the C282Y mutation of the HFE gene may have altered hematology indices and higher iron stores than wild-type subjects. Methods: We performed a cross-sectional analysis of 1488 females and 1522 males 20-79 years of age drawn from the Busselton (Australia) population study to assess the effects of HFE genotype, age, gender, and lifestyle on serum iron and hematology indices. Results: Male C282Y heterozygotes had increased transferrin saturation compared with the wild-type genotype. Neither male nor female heterozygotes had significantly increased ferritin values compared with the wild-type genotype. Younger (20-29 years) wild-type males, but not heterozygous males, had significantly lower ferritin values than wild-type males in the older age groups. Compound heterozygous subjects had increased means for serum iron, transferrin saturation, corpuscular volume, and corpuscular hemoglobin compared with the wild-type genotype, and the males also had increased ferritin values (medians 323 vs 177 mug/L; P = 0.003). In both male and female wild-type subjects, an increased body mass index was associated with decreased serum iron and transferrin saturation and increased ferritin values. There was a significant increase in ferritin concentrations in both genders with increasing frequency of red meat consumption above a baseline of 1-2 times per week and alcohol intakes >10 g/day. Conclusions: Male C282Y heterozygotes had significantly increased transferrin saturation values. Compound heterozygous (C282Y/H63D) subjects formed a separate category of C282Y heterozygotes in whom both iron and red cell indices were significantly increased compared with the wild-type genotype. (C) 2001 American Association for Clinical Chemistry.