Development of a semimechanistic model to describe the pharmacokinetics of gentamicin in patients receiving Hemodialysis

The aim of this study was to ascertain the most suitable dosing schedule for gentamicin in patients receiving hemodialysis. We developed a model to describe the concentrationtime course of gentamicin in patients receiving hemodialysis. Using the model, an optimal dosing schedule was evaluated. Various dosing regimens were compared in their ability to achieve maximum concentration (C-max, >= 8 mg/L) and area under the concentration time-curve (AUC >= 70 mg(.)h/L and <= 120 mg(.)h/L per 24 hours). The model was evaluated by comparing model predictions against real data collected retrospectively. Simulations from the model confirmed the benefits of predialysis dosing. The mean optimal dose was 230 mg administered immediately before dialysis. The model was found to have good predictive performance when simulated data were compared to data observed in real patients. In summary, a model was developed that describes gentamicin pharmacokinetics in patients receiving hemodialysis. Predialysis dosing provided a superior pharmacokinetic profile than did postdialysis dosing.

The cost of uniqueness in groundwater model calibration

Calibration of a groundwater model requires that hydraulic properties be estimated throughout a model domain. This generally constitutes an underdetermined inverse problem, for which a Solution can only be found when some kind of regularization device is included in the inversion process. Inclusion of regularization in the calibration process can be implicit, for example through the use of zones of constant parameter value, or explicit, for example through solution of a constrained minimization problem in which parameters are made to respect preferred values, or preferred relationships, to the degree necessary for a unique solution to be obtained. The cost of uniqueness is this: no matter which regularization methodology is employed, the inevitable consequence of its use is a loss of detail in the calibrated field. This, ill turn, can lead to erroneous predictions made by a model that is ostensibly well calibrated. Information made available as a by-product of the regularized inversion process allows the reasons for this loss of detail to be better understood. In particular, it is easily demonstrated that the estimated value for an hydraulic property at any point within a model domain is, in fact, a weighted average of the true hydraulic property over a much larger area. This averaging process causes loss of resolution in the estimated field. Where hydraulic conductivity is the hydraulic property being estimated, high averaging weights exist in areas that are strategically disposed with respect to measurement wells, while other areas may contribute very little to the estimated hydraulic conductivity at any point within the model domain, this possibly making the detection of hydraulic conductivity anomalies in these latter areas almost impossible. A study of the post-calibration parameter field covariance matrix allows further insights into the loss of system detail incurred through the calibration process to be gained. A comparison of pre- and post-calibration parameter covariance matrices shows that the latter often possess a much smaller spectral bandwidth than the former. It is also demonstrated that, as all inevitable consequence of the fact that a calibrated model cannot replicate every detail of the true system, model-to-measurement residuals can show a high degree of spatial correlation, a fact which must be taken into account when assessing these residuals either qualitatively, or quantitatively in the exploration of model predictive uncertainty. These principles are demonstrated using a synthetic case in which spatial parameter definition is based oil pilot points, and calibration is Implemented using both zones of piecewise constancy and constrained minimization regularization. (C) 2005 Elsevier Ltd. All rights reserved.

Type III secretion contact-dependent model for the intracellular development of Chlamydia

The medically significant genus Chlamydia is a class of obligate intracellular bacterial pathogens that replicate within vacuoles in host eukaryotic cells termed inclusions. Chlamydia's developmental cycle involves two forms; an infectious extracellular form, known as an elementary body (EB), and a non-infectious form, known as the reticulate body (RB), that replicates inside the vacuoles of the host cells. The RB surface is covered in projections that are in intimate contact with the inclusion membrane. Late in the developmental cycle, these reticulate bodies differentiate into the elementary body form. In this paper, we present a hypothesis for the modulation of these developmental events involving the contact-dependent type III secretion (TTS) system. TTS surface projections mediate intimate contact between the RB and the inclusion membrane. Below a certain number of projections, detachment of the RB provides a signal for late differentiation of RB into EB. We use data and develop a mathematical model investigating this hypothesis. If the hypothesis proves to be accurate, then we have shown that increasing the number of inclusions per host cell will increase the number of infectious progeny EB until some optimal number of inclusions. For more inclusions than this optimum, the infectious yield is reduced because of spatial restrictions. We also predict that a reduction in the number of projections on the surface of the RB (and as early as possible during development) will significantly reduce the burst size of infectious EB particles. Many of the results predicted by the model can be tested experimentally and may lead to the identification of potential targets for drug design. © Society for Mathematical Biology 2006.

The genetics of alcohol intake and of alcohol dependence

Background: Because alcohol has multiple dose-dependent consequences, it is important to understand the causes of individual variation in the amount of alcohol used. The aims of this study were to assess the long-term repeatability and genetic or environmental causes of variation in alcohol intake and to estimate the degree of overlap with causes of susceptibility to alcohol dependence. Methods: Data were used from three studies conducted between 1980 and 1995 on volunteer adult male and female Australian twin subjects. In each study, alcohol intake was reported both as quantity X frequency and as past-week data. Repeatability was calculated as correlations between occasions and between measures, and the effects of genes and environment were estimated by multivariate model fitting to the twin pair repeated measures of alcohol use. Relationships between mean alcohol use and the lifetime history of DSM-III-R alcohol dependence were tested by bivariate model fitting. Results: Repeatability of the alcohol intake measures was between 0.54 and 0.85, with the highest repeatability between measures within study and the lowest repeatability between the first and last studies. Reported alcohol consumption was mainly affected by genetic factors affecting all times of study and by nonshared environmental factors (including measurement error) unique to each time of study. Genes that affect alcohol intake do affect alcohol dependence, but genetic effects unique to dependence are also significant; environmental effects are largely unique to either intake and dependence. Conclusions: Nearly all the repeatable component of variation in alcohol intake is due to genetic effects. Genes affecting intake also affect dependence risk, but there are other genes that affect dependence alone. Studies aiming to identify genes that affect alcohol use disorders need to test loci and candidate genes against both phenotypes.

Monarchs across the Pacific: the Columbus hypothesis revisited

The 'Columbus hypothesis' suggests that the annual north-south return migration of Danaus plexippus in North America is a very recently evolved behaviour, less than 200 years old. This hypothesis rests, in part, on an analysis of the 19th century spread of the monarch across the Pacific that assumes a continuous east to west movement and is based predominantly on one publication. We review all the contemporary literature and present new analysis of the data. The movement of the monarch across the Pacific in the second half of the 19th century is best explained by a model which involves no more than three spot introductions, directly or indirectly aided by human movement, followed by natural spread of the monarch across island groups. Contemporary records refer to 'boom' and 'bust' population cycles on newly settled islands, which may have led to high rates of monarch movement. We see no evidence in the records to suggest an east to west sweep by monarch populations as suggested by the Columbus hypothesis. (C) 2004 The Linnean Society of London.

Which version of the Eysenck Personality Profiler is best? 6-, 12- or 20-items per scale

Data provided by 400 first year undergraduate students were analysed to develop two short forms of the Eysenck Personality Profiler (EPP) in which each of the 22 primary scales is assessed by a 6-item and a 12-item version instead of the usual 20-item per scale measure. In comparison with the 6-item per scale measure, the 12-item version retains more of the characteristics of the long version and seems a good compromise between quality of data and administration time. (C) 2004 Elsevier Ltd. All rights reserved.

Tractable molecular theory of transport of Lennard-Jones fluids in nanopores

We present here a tractable theory of transport of simple fluids in cylindrical nanopores, which is applicable over a wide range of densities and pore sizes. In the Henry law low-density region the theory considers the trajectories of molecules oscillating between diffuse wall collisions, while at higher densities beyond this region the contribution from viscous flow becomes significant and is included through our recent approach utilizing a local average density model. The model is validated by means of equilibrium as well nonequilibrium molecular dynamics simulations of supercritical methane transport in cylindrical silica pores over a wide range of temperature, density, and pore size. The model for the Henry law region is exact and found to yield an excellent match with simulations at all conditions, including the single-file region of very small pore size where it is shown to provide the density-independent collective transport coefficient. It is also shown that in the absence of dispersive interactions the model reduces to the classical Knudsen result, but in the presence of such interactions the latter model drastically overpredicts the transport coefficient. For larger micropores beyond the single-file region the transport coefficient is reduced at high density because of intermolecular interactions and hindrance to particle crossings leading to a large decrease in surface slip that is not well represented by the model. However, for mesopores the transport coefficient increases monotonically with density, over the range studied, and is very well predicted by the theory, though at very high density the contribution from surface slip is slightly overpredicted. It is also seen that the concept of activated diffusion, commonly associated with diffusion in small pores, is fundamentally invalid for smooth pores, and the apparent activation energy is not simply related to the minimum pore potential or the adsorption energy as generally assumed. (C) 2004 American Institute of Physics.

Modeling molecular transport in slit pores

We examine the transport of methane in microporous carbon by performing equilibrium and nonequilibrium molecular dynamics simulations over a range of pore sizes, densities, and temperatures. We interpret these simulation results using two models of the transport process. At low densities, we consider a molecular flow model, in which intermolecular interactions are neglected, and find excellent agreement between transport diffusion coefficients determined from simulation, and those predicted by the model. Simulation results indicate that the model can be applied up to fluid densities of the order to 0.1-1 nm(-3). Above these densities, we consider a slip flow model, combining hydrodynamic theory with a slip condition at the solid-fluid interface. As the diffusion coefficient at low densities can be accurately determined by the molecular flow model, we also consider a model where the slip condition is supplied by the molecular flow model. We find that both density-dependent models provide a useful means of estimating the transport coefficient that compares well with simulation. (C) 2004 American Institute of Physics.

Thermosetting Blends of Polybenzoxazine and Poly(ε-caprolactone): Phase Behavior and Intermolecular Specific Interactions

Polybenzoxazine (PBA-a)/poly(epsilon-caprolactone) (PCL) blends were prepared by an in situ curing reaction of benzoxazine (BA-a) in the presence of PCL. Before curing, the benzoxazine (BA-a)/PCL blends are miscible, which was evidenced by the behaviors of single and composition-dependant glass transition temperature and equilibrium melting point depression. However, the phase separation induced by polymerization was observed after curing at elevated temperature. It was expected that after curing, the PBA-a/PCL blends would be miscible since the phenolic hydroxyls in the PBA-a molecular backbone have the potential to form inter- molecular hydrogen-bonding interactions with the carbonyls of PCL and thus would fulfil the miscibility of the blends. The resulting morphology of the blends prompted an investigation of the status of association between PBA-a and PCL under the curing conditions. Although Fourier-transform infrared spectroscopy (FT-IR) showed that there were intermolecular hydrogen-bonding interactions between PBA-a and PCL at room temperature, especially for the PCL-rich blends, the results of variable temperature FT-IR spectroscopy by the model compound indicate that the phenolic hydroxyl groups could not form efficient intermolecular hydrogen-bonding interactions at elevated temperatures, i.e., the phenolic hydroxyl groups existed mainly in the non-associated form in the system during curing. The results are valuable to understand the effect of curing temperature on the resulting morphology of the thermosetting blends. SEM micrograph of the dichloromethane-etched fracture surface of a 90:10 PBA-a PCL blend showing a heterogeneous morphology.

Increased site 1 affinity improves biopotency of porcine growth hormone - Evidence against diffusion dependent receptor dimerization

Based on phage display optimization studies with human growth hormone (GH), it is thought that the biopotency of GH cannot be increased. This is proposed to be a result of the affinity of the first receptor for hormone far exceeding that which is required to trap the hormone long enough to allow diffusion of the second receptor to form the ternary complex, which initiates signaling. We report here that despite similar site 1 kinetics to the hGH/hGH receptor interaction, the potency of porcine GH for its receptor can be increased up to 5-fold by substituting hGH residues involved in site 1 binding into pGH. Based on extensive mutations and BIAcore studies, we show that the higher potency and site 1 affinity of hGH for the pGHR is primarily a result of a decreased off-rate associated with residues in the extended loop between helices 1 and 2 that interact with the two key tryptophans Trp(104) and Trp(169) in the receptor binding hot spot. Our mutagenic analysis has also identified a second determinant (Lys(165)), which in addition to His(169), restricts the ability of non-primate hormones to activate hGH receptor. The increased biopotency of GH that we observe can be explained by a model for GH receptor activation where subunit alignment is critical for effective signaling.

Adsorbate transport in nanopores

We present a tractable theory of transport of simple fluids in cylindrical nanopores, considering trajectories of molecules between diffuse wall collisions at low-density, and including viscous flow contributions at higher densities. The model is validated through molecular dynamics simulations of supercritical methane transport, over a wide range of conditions. We find excellent agreement between model and simulation at low to medium densities. However, at high densities the model tends to over-predict the transport behaviour, due to a large decrease in surface slip that is not well represented by the model. It is also seen that the concept of activated diffusion, commonly associated with diffusion in small pores, is fundamentally invalid for smooth pores.

Prediction of high-pressure adsorption equilibrium of supercritical gases using density functional theory

In this paper, we present the results of the prediction of the high-pressure adsorption equilibrium of supercritical. gases (Ar, N-2, CH4, and CO2) on various activated carbons (BPL, PCB, and Norit R1 extra) at various temperatures using a density-functional-theory-based finite wall thickness (FWT) model. Pore size distribution results of the carbons are taken from our recent previous work 1,2 using this approach for characterization. To validate the model, isotherms calculated from the density functional theory (DFT) approach are comprehensively verified against those determined by grand canonical Monte Carlo (GCMC) simulation, before the theoretical adsorption isotherms of these investigated carbons calculated by the model are compared with the experimental adsorption measurements of the carbons. We illustrate the accuracy and consistency of the FWT model for the prediction of adsorption isotherms of the all investigated gases. The pore network connectivity problem occurring in the examined carbons is also discussed, and on the basis of the success of the predictions assuming a similar pore size distribution for accessible and inaccessible regions, it is suggested that this is largely related to the disordered nature of the carbon.

Subjective and predicted sleepiness while driving in young adults

Sleepiness is a significant contributor to car crashes and sleepiness related crashes have higher mortality and morbidity than other crashes. Young adult drivers are at particular risk for sleepiness related car crashes. It has been suggested that this is because young adults are typically sleepier than older adults because of chronic sleep loss, and more often drive at times of increased risk of acute sleepiness. This prospective study aimed to determine the relationship between predicted and perceived sleepiness while driving in 47 young-adult drivers over a 4-week period. Sleepiness levels were predicted by a model incorporating known circadian and sleep factors influencing alertness, and compared to subjective ratings of sleepiness during 25 18 driving episodes. Results suggested that young drivers frequently drive while at risk of crashing, at times of predicted sleepiness (>7% of episodes) and at times they felt themselves to be sleepy (>23% of episodes). A significant relationship was found between perceived and predicted estimates of sleepiness. However, the participants nonetheless drove at these times. The results of this study may help preventative programs to specifically target factors leading to increased sleepiness when driving (particularly time of day), and to focus interventions to stop young adults from driving when they feel sleepy. (c) 2005 Elsevier Ltd. All rights reserved.

Divergent selection and the evolution of signal traits and mating preferences

Mating preferences are common in natural populations, and their divergence among populations is considered an important source of reproductive isolation during speciation. Although mechanisms for the divergence of mating preferences have received substantial theoretical treatment, complementary experimental tests are lacking. We conducted a laboratory evolution experiment, using the fruit fly Drosophila serrata, to explore the role of divergent selection between environments in the evolution of female mating preferences. Replicate populations of D. serrata were derived from a common ancestor and propagated in one of three resource environments: two novel environments and the ancestral laboratory environment. Adaptation to both novel environments involved changes in cuticular hydrocarbons, traits that predict mating success in these populations. Furthermore, female mating preferences for these cuticular hydrocarbons also diverged among populations. A component of this divergence occurred among treatment environments, accounting for at least 17.4% of the among- population divergence in linear mating preferences and 17.2% of the among-population divergence in nonlinear mating preferences. The divergence of mating preferences in correlation with environment is consistent with the classic by- product model of speciation in which premating isolation evolves as a side effect of divergent selection adapting populations to their different environments.

Thermal strain in the mushy zone related to hot tearing

A volume-averaged two-phase model addressing the main transport phenomena associated with hot tearing in an isotropic mushy zone during solidification of metallic alloys has recently been presented elsewhere along with a new hot tearing criterion addressing both inadequate melt feeding and excessive deformation at relatively high solid fractions. The viscoplastic deformation in the mushy zone is addressed by a model in which the coherent mush is considered as a porous medium saturated with liquid. The thermal straining of the mush is accounted for by a recently developed model taking into account that there is no thermal strain in the mushy zone at low solid fractions because the dendrites then are free to move in the liquid, and that the thermal strain in the mushy zone tends toward the thermal strain in the fully solidified material when the solid fraction tends toward one. In the present work, the authors determined how variations in the parameters of the constitutive equation for thermal strain influence the hot tearing susceptibility calculated by the criterion. It turns out that varying the parameters in this equation has a significant effect on both liquid pressure drop and viscoplastic strain, which are key parameters in the hot tearing criterion. However, changing the parameters in this constitutive equation will result in changes in the viscoplastic strain and the liquid pressure drop that have opposite effects on the hot tearing susceptibility. The net effect on the hot tearing susceptibility is thus small.