The mechanical properties of Saccharomyces cerevisiae

Cell-wall mechanical properties play an integral part in the growth and form of Saccharomyces cerevisiae, In contrast to the tremendous knowledge on the genetics of S. cerevisiae, almost nothing is known about its mechanical properties. We have developed a micromanipulation technique to measure the force required to burst single cells and have recently established a mathematical model to extract the mechanical properties of the cell wall from such data, Here we determine the average surface modulus of the S, cerevisiae cell wall to be 11.1 +/- 0.6 N/m and 12.9 +/- 0.7 N/m in exponential and stationary phases, respectively, giving corresponding Young's moduli of 112 +/- 6 MPa and 107 +/- 6 MPa, This result demonstrates that yeast cell populations strengthen as they enter stationary phase by increasing wall thickness and hence the surface modulus, without altering the average elastic properties of the cell-wall material. We also determined the average breaking strain of the cell wall to be 82% +/- 3% in exponential phase and 80% +/- 3% in stationary phase, This finding provides a failure criterion that can be used to predict when applied stresses (e,g,, because of fluid flow) will lead to wall rupture, This work analyzes yeast compression experiments in different growth phases by using engineering methodology.

Structured modeling of recombinant protein production in batch and fed-batch culture of baculovirus-infected insect cells

The infection of insect cells with baculovirus was described in a mathematical model as a part of the structured dynamic model describing whole animal cell metabolism. The model presented here is capable of simulating cell population dynamics, the concentrations of extracellular and intracellular viral components, and the heterologous product titers. The model describes the whole processes of viral infection and the effect of the infection on the host cell metabolism. Dynamic simulation of the model in batch and fed-batch mode gave good agreement between model predictions and experimental data. Optimum conditions for insect cell culture and viral infection in batch and fed-batch culture were studied using the model.

Finite element modelling of heat transfer through permeable cracks in hydrothermal systems with upward throughflow

We use the finite element method to model the heat transfer phenomenon through permeable cracks in hydrothermal systems with upward throughflow. Since the finite element method is an approximate numerical method, the method must be validated before it is used to soh,e any new, kind of problem. However, the analytical solution, which can be used to validate the finite element method and other numerical methods, is rather limited in the literature, especially, for the problem considered here. Keeping this in mind, we have derived analytical solutions for the temperature distribution along the vertical axis of a crack in a fluid-saturated porous layer. After the finite element method is validated by comparing the numerical solution with the analytical solution for the same benchmark problem, it is used to investigate the pore-fluid flow and heat transfer in layered hydrothermal systems with vertical permeable cracks. The related analytical and numerical results have demonstrated that vertical cracks are effective and efficient members to transfer heat energy from the bottom section to the top section in hydrothermal systems with upward throughflow.

Three-dimensional electrical impedance, tomography: A topology optimization approach

Electrical impedance tomography is a technique to estimate the impedance distribution within a domain, based on measurements on its boundary. In other words, given the mathematical model of the domain, its geometry and boundary conditions, a nonlinear inverse problem of estimating the electric impedance distribution can be solved. Several impedance estimation algorithms have been proposed to solve this problem. In this paper, we present a three-dimensional algorithm, based on the topology optimization method, as an alternative. A sequence of linear programming problems, allowing for constraints, is solved utilizing this method. In each iteration, the finite element method provides the electric potential field within the model of the domain. An electrode model is also proposed (thus, increasing the accuracy of the finite element results). The algorithm is tested using numerically simulated data and also experimental data, and absolute resistivity values are obtained. These results, corresponding to phantoms with two different conductive materials, exhibit relatively well-defined boundaries between them, and show that this is a practical and potentially useful technique to be applied to monitor lung aeration, including the possibility of imaging a pneumothorax.

Diel vertical migration: Modelling light-mediated mechanisms

Light is generally regarded as the most likely cue used by zooplankton to regulate their vertical movements through the water column. However, the way in which light is used by zooplankton as a cue is not well understood. In this paper we present a mathematical model of diel vertical migration which produces vertical distributions of zooplankton that vary in space and time. The model is used to predict the patterns of vertical distribution which result when animals are assumed to adopt one of three commonly proposed mechanisms for vertical swimming. First, we assume zooplankton tend to swim towards a preferred intensity of light. We then assume zooplankton swim in response to either the rate of change in light intensity or the relative rate of change in light intensity. The model predicts that for all three mechanisms movement is fastest at sunset and sunrise and populations are primarily influenced by eddy diffusion at night in the absence of a light stimulus. Daytime patterns of vertical distribution differ between the three mechanisms and the reasons for the predicted differences are discussed. Swimming responses to properties of the light field are shown to be adequate for describing diel vertical migration where animals congregate in near surface waters during the evening and reside at deeper depths during the day. However, the model is unable to explain how some populations halt their ascent before reaching surface waters or how populations re-congregate in surface waters a few hours before sunrise, a phenomenon which is sometimes observed in the held. The model results indicate that other exogenous or endogenous factors besides light may play important roles in regulating vertical movement.

The kinetics of baculovirus adsorption to insect cells in suspension culture

The influence of various culture parameters on the attachment of a recombinant baculovirus to suspended insect cells was examined under normal culture conditions. These parameters included cell density, multiplicity of infection, and composition of the cell growth medium. It was found that the fractional rate of virus attachment was independent of the multiplicity of infection but dependent on the cell density. A first order mathematical model was used to simulate the adsorption kinetics and predict the efficiency of virus attachment under the various culture conditions. This calculated efficiency of virus attachment was observed to decrease at high cell densities, which was attributed to cell clumping. It was also observed that virus attachment was more efficient in Sf900II serum free medium than it was in IPL-41 serum-supplemented medium. This effect was attributed to the protein in serum which may coat the cells and so inhibit adsorption. A general discussion relating the observations made in-these experiments to the kinetics of recombinant baculovirus adsorption to suspended insect cells is presented.

Odorant-induced currents in intact patches from rat olfactory receptor neurons: Theory and experiment

Odorant-induced currents in mammalian olfactory receptor neurons have proved difficult to obtain reliably using conventional whole-cell recording. By using a mathematical model of the electrical circuit of the patch and rest-of-cell, we demonstrate how cell-attached patch measurements can be used to quantitatively analyze responses to odorants or a high (100 mM) K+ solution. High K+ induced an immediate current flux from cell to pipette, which was modeled as a depolarization of similar to 52 mV, close to that expected from the Nernst equation (56 mV), and no change in the patch conductance. By contrast, a cocktail of cAMP-stimulating odorants induced a current flux from pipette into cell following a significant (4-10 s) delay. This was modeled as an average patch conductance increase of 36 pS and a depolarization of 13 mV, Odorant-induced single channels had a conductance of 16 pS. In cells bathed with no Mg2+ and 0.25 mM Ca2+, odorants induced a current flow from cell to pipette, which was modeled as a patch conductance increase of similar to 115 pS and depolarization of similar to 32 mV, All these results are consistent with cAMP-gated cation channels dominating the odorant response, This approach, which provides useful estimates of odorant-induced voltage and conductance changes, is applicable to similar measurements in any small cells.

The Potential Impact of Using Donations After Cardiac Death on the Liver Transplantation Program and Waiting List in the State of Sao Paulo, Brazil

Liver transplantation was first performed at the University of Sao Paulo School of Medicine in 1968. Since then, the patient waiting list for liver transplantation has increased at a rate of 150 new cases per month. Liver transplantation itself rose 1.84-fold (from 160 to 295) from 1988 to 2004. However, the number of patients on the liver waiting list jumped 2.71-fold (from 553 to 1500). Consequently, the number of deaths on the liver waiting list moved to a higher level, from 321 to 671, increasing 2.09-fold. We have applied a mathematical model to analyze the potential impact of using a donation after cardiac death (DCD) policy on our liver transplantation program and on the waiting list. Five thousand one hundred people died because of accidents and other violent causes in our state in 2004; of these, only 295 were donors of liver grafts that were transplanted. The model assumed that 5% of these grafts would have been DCD. We found a relative reduction of 27% in the size of the liver transplantation waiting list if DCD had been used by assuming that 248 additional liver transplants would have been performed annually. In conclusion, the use of DCD in our transplantation program would reduce the pressure on our liver transplantation waiting list, reducing it by at least 27%. On the basis of this model, the projected number of averted deaths is about 41,487 in the next 20 years. Liver Transpl 14:1732-1736, 2008. (C) 2008 AASLD.

Cost-effectiveness analysis of a hypothetical hepatitis C vaccine compared to antiviral therapy

We propose a mathematical model to simulate the dynamics of hepatitis C virus (HCV) infection in the state of Sao Paulo, Brazil. We assumed that a hypothetical vaccine, which cost was taken to be the initial cost of the vaccine against hepatitis B exists and it is introduced in the model. We computed its cost-effectiveness compared with the anti-HCV therapy. The calculated basic reproduction number was 1.20. The model predicts that without intervention a steady state exists with an HCV prevalence of 3%, in agreement with the Current epidemiological data. Starting from this steady state three interventions were simulated: indiscriminate vaccination, selective vaccination and anti-HCV therapy. Selective vaccination proved to be the strategy with the best cost-effectiveness ratio, followed by indiscriminate vaccination and anti-HCV therapy.

Autopsy-Proven Determinants of Immediate Fire Death in Lungs

In immediate fire deaths, pulmonary injury may be the main source of mortality, being important to document the histologic findings for the purpose of excluding other modes of death, such as from asphyxia with no gross findings. In this context, a group of morphologic determinants have been targeted with useful makers of pulmonary injury. To facilitate the determination of whether an individual was deceased before the start of a fire and validate the importance of parenchymal alterations in pulmonary injury in fire deaths, we studied lungs in victims of fire (N = 28) and suffocation (N = 40), creating a mathematical model using cluster analysis. For this purpose, a semiquantitative analysis of the distal parenchyma was performed to evaluate the amount of bronchiolar dilatation, overinsufflation (ductal and alveolar), collapse (ductal and alveolar), passive congestion, alveolar edema, and hemorrhage (interstitial and alveolar). These 7 histologic determinants were useful to discriminate fire (bronchiolar dilatation, ductal overinsuflation, alveolar overinsuflation, alveolar hemorrhage) from suffocation lung injuries (alveolar collapse, congestion, and edema). We conclude that these determinants should be included in the routine of forensic pathology.

The risk of chikungunya fever in a dengue-endemic area

Background. Chikungunya, an alphavirus of the Togaviridae family, causes a febrile disease transmitted to humans by the bite of infected Aedes mosquitoes. This infection is reaching endemic levels in many Southeast Asian countries. Symptoms include sudden onset of fever, chills, headache, nausea, vomiting, joint pain with or without swelling, low back pain, and rash. According to the World Health Organization, there are 2 billion people living in Aedes-infested areas. In addition, traveling to these areas is popular, making the potential risk of infections transmitted by the bite of infected Aedes mosquitoes very high. Methods. We proposed a mathematical model to estimate the risk of acquiring chikungunya fever in an Aedes-infested area by taking the prevalence of dengue fever into account. The basic reproduction number for chikungunya fever R-0chik can be written as a function of the basic reproduction number of dengue R-0dengue by calculating the ratio R-0chik/R-0dengue. From R-0chik, we estimated the force of infection and the risk of acquiring the disease both for local residents of a dengue-endemic area and for travelers to this area. Results. We calculated that R-0chik is 64.4% that of R-0dengue. The model was applied to a hypothetical situation, namely, estimating the individual risk of acquiring chikungunya fever in a dengue-endemic area, both for local inhabitants (22% in steady state) and for visiting travelers (from 0.31% to 1.23% depending on the time spent in the area). Conclusions. The method proposed based on the output of a dynamical model is innovative and provided an estimation of the risk of infection, both for local inhabitants and for visiting travelers.

Modelling the control strategies against dengue in Singapore

Notified cases of dengue infections in Singapore reached historical highs in 2004 (9459 cases) and 2005 (13 817 cases) and the reason for such all increase is still to be established. We apply a mathematical model for dengue infection that takes into account the seasonal variation in incidence, characteristic of dengue fever, and which mimics the 2004-2005 epidemics in Singapore. We simulated a set of possible control strategies and confirmed the intuitive belief that killing adult mosquitoes is the most effective strategy to control an ongoing epidemic. On the other hand, the control of immature forms was very efficient ill preventing the resurgence of dengue epidemics. Since the control of immature forms allows the reduction of adulticide, it seems that the best strategy is to combine both adulticide and larvicide control measures during an outbreak, followed by the maintenance of larvicide methods after the epidemic has subsided. In addition, the model showed that the mixed strategy of adulticide and larvicide methods introduced by the government seems to be very effective in reducing the number of cases in the first weeks after the start of control.

FOR DEBATE A hypothesis for the 2007 dengue outbreak in Singapore

A previous mathematical model explaining dengue in Singapore predicted a reasonable outbreak of about 6500 cases for 2006 and a very mild outbreak with about 2000 cases for 2007. However, only 3051 cases were reported in 2006 while more than 7800 were reported in the first 44 weeks of 2007. We hypothesized that the combination of haze with other local sources of particulate matter had a significant impact on mosquito life expectancy, significantly increasing their mortality rate. To test the hypothesis a mathematical model based on the reproduction number of dengue fever and aimed at comparing the impact of several possible alternative control strategies was proposed. This model also aimed at contributing to the understanding of the causes of dengue resurgence in Singapore in the last decade. The model`s simulation demonstrated that an increase in mosquito mortality in 2006 and either a reduction in mortality or an increase in the carrying capacity of mosquitoes in 2007 explained the patterned observed in Singapore. Based on the model`s simulation we concluded that the fewer than expected number of dengue cases in Singapore in 2006 was caused by an increase in mosquito mortality due to the disproportionate haze affecting the country that year and that particularly favourable environmental conditions in 2007 propitiated mosquitoes with a lower mortality rate, which explains the greater than expected number of dengue cases in 2007. Whether our hypothesis is plausible or not should be debated further.

Effect of PaCO(2) variation on standard base excess value in critically ill patients

Purpose: The aim of this study was to investigate the impact of acute PaCO(2) temporal variation on the standard base excess (SBE) value in critically ill patients. Methods: A total of 265 patients were prospectively observed; 158 were allocated to the modeling group, and 107 were allocated to the validation group. Two models were developed in the modeling group (one including and one excluding PaCO(2) as a variable determinant of SBE), and both were tested in the validation group. Results: In the modeling group, the mathematical model including SIDai, SIG, L-lactate, albumin, phosphate, and PaCO(2) had a predictive superiority in comparison with the model without PaCO(2) (R(2) = 0.978 and 0.916, respectively). In the validation group, the results were confirmed with significant F change statistics (R(2) change = 0.059, P < .001) between the model with and without PaCO(2). A high correlation (R = 0.99, P < .001) and agreement (bias = -0.25 mEq/L, limits of agreement 95% = -0.72 to 0.22 mEq/L) were found between the model-predicted SBE value and the SBE calculated using the Van Slyke equation. Conclusions: Acute PaCO(2), temporal variation is related to SBE changes in critically ill patients. (C) 2009 Elsevier Inc. All rights reserved.

Modeling upper eyelid kinematics during spontaneous and reflex blinks

To test a mathematical model for measuring blinking kinematics. Spontaneous and reflex blinks of 23 healthy subjects were recorded with two different temporal resolutions. A magnetic search coil was used to record 77 blinks sampled at 200 Hz and 2 kHz in 13 subjects. A video system with low temporal resolution (30 Hz) was employed to register 60 blinks of 10 other subjects. The experimental data points were fitted with a model that assumes that the upper eyelid movement can be divided into two parts: an impulsive accelerated motion followed by a damped harmonic oscillation. All spontaneous and reflex blinks, including those recorded with low resolution, were well fitted by the model with a median coefficient of determination of 0.990. No significant difference was observed when the parameters of the blinks were estimated with the under-damped or critically damped solutions of the harmonic oscillator. On the other hand, the over-damped solution was not applicable to fit any movement. There was good agreement between the model and numerical estimation of the amplitude but not of maximum velocity. Spontaneous and reflex blinks can be mathematically described as consisting of two different phases. The down-phase is mainly an accelerated movement followed by a short time that represents the initial part of the damped harmonic oscillation. The latter is entirely responsible for the up-phase of the movement. Depending on the instantaneous characteristics of each movement, the under-damped or critically damped oscillation is better suited to describe the second phase of the blink. (C) 2010 Elsevier B.V. All rights reserved.