Modelling acidosis and the cell cycle in multicellular tumour spheroids

A partial differential equation model is developed to understand the effect that nutrient and acidosis have on the distribution of proliferating and quiescent cells and dead cell material (necrotic and apopotic) within a multicellular tumour spheroid. The rates of cell quiescence and necrosis depend upon the local nutrient and acid concentrations and quiescent cells are assumed to consume less nutrient and produce less acid than proliferating cells. Analysis of the differences in nutrient consumption and acid production by quiescent and proliferating cells shows low nutrient levels do not necessarily lead to increased acid concentration via anaerobic metabolism. Rather, it is the balance between proliferating and quiescent cells within the tumour which is important; decreased nutrient levels lead to more quiescent cells, which produce less acid than proliferating cells. We examine this effect via a sensitivity analysis which also includes a quantification of the effect that nutrient and acid concentrations have on the rates of cell quiescence and necrosis.

Classifying general nonlinear force laws in cell-based models via the continuum limit

though discrete cell-based frameworks are now commonly used to simulate a whole range of biological phenomena, it is typically not obvious how the numerous different types of model are related to one another, nor which one is most appropriate in a given context. Here we demonstrate how individual cell movement on the discrete scale modeled using nonlinear force laws can be described by nonlinear diffusion coefficients on the continuum scale. A general relationship between nonlinear force laws and their respective diffusion coefficients is derived in one spatial dimension and, subsequently, a range of particular examples is considered. For each case excellent agreement is observed between numerical solutions of the discrete and corresponding continuum models. Three case studies are considered in which we demonstrate how the derived nonlinear diffusion coefficients can be used to (a) relate different discrete models of cell behavior; (b) derive discrete, intercell force laws from previously posed diffusion coefficients, and (c) describe aggregative behavior in discrete simulations.

Methodology To model the energy and greenhouse gas emissions of electronic software distributions

A new electronic software distribution (ESD) life cycle analysis (LCA)methodology and model structure were constructed to calculate energy consumption and greenhouse gas (GHG) emissions. In order to counteract the use of high level, top-down modeling efforts, and to increase result accuracy, a focus upon device details and data routes was taken. In order to compare ESD to a relevant physical distribution alternative,physical model boundaries and variables were described. The methodology was compiled from the analysis and operational data of a major online store which provides ESD and physical distribution options. The ESD method included the calculation of power consumption of data center server and networking devices. An in-depth method to calculate server efficiency and utilization was also included to account for virtualization and server efficiency features. Internet transfer power consumption was analyzed taking into account the number of data hops and networking devices used. The power consumed by online browsing and downloading was also factored into the model. The embedded CO2e of server and networking devices was proportioned to each ESD process. Three U.K.-based ESD scenarios were analyzed using the model which revealed potential CO2e savings of 83% when ESD was used over physical distribution. Results also highlighted the importance of server efficiency and utilization methods.

Upper-ocean heat budget and ocean eddy transport in the south-east Pacific in a high-resolution coupled model

We present an analysis of the oceanic heat advection and its variability in the upper 500 m in the southeastern tropical Pacific (100W–75W, 25S–10S) as simulated by the global coupled model HiGEM, which has one of the highest resolutions currently used in long-term integrations. The simulated climatology represents a temperature advection field arising from transient small-scale (<450 km) features, with structures and transport that appear consistent with estimates based on available observational data for the mooring at 20S, 85W. The transient structures are very persistent (>4 months), and in specific locations they generate an important contribution to the local upper-ocean heat budget, characterised by scales of a few hundred kilometres, and periods of over a year. The contribution from such structures to the local, long-term oceanic heat budget however can be of either sign, or vanishing, depending on the location; and, although there appears some organisation in preferential areas of activity, the average over the entire region is small. While several different mechanisms may be responsible for the temperature advection by transients, we find that a significant, and possibly dominant, component is associated with vortices embedded in the large-scale, climatological salinity gradient associated with the fresh intrusion of mid-latitude intermediate water which penetrates north-westward beneath the tropical thermocline

Riluzole neuroprotection in a parkinson’s disease model involves suppression of reactive astrocytosis but not GLT-1 regulation

Background Riluzole is a neuroprotective drug used in the treatment of motor neurone disease. Recent evidence suggests that riluzole can up-regulate the expression and activity of the astrocyte glutamate transporter, GLT-1. Given that regulation of glutamate transport is predicted to be neuroprotective in Parkinson's disease, we tested the effect of riluzole in parkinsonian rats which had received a unilateral 6-hydroxydopamine injection into the median forebrain bundle. Results Rats were treated with intraperitoneal riluzole (4 mg/kg or 8 mg/kg), 1 hour before the lesion then once daily for seven days. Riluzole produced a modest but significant attenuation of dopamine neurone degeneration, assessed by suppression of amphetamine-induced rotations, preservation of tyrosine hydroxylase positive neuronal cell bodies in the substantia nigra pars compacta and attenuation of striatal tyrosine hydroxylase protein loss. Seven days after 6-hydroxydopamine lesion, reactive astrocytosis was observed in the striatum, as determined by increases in expression of glial fibrillary acidic protein, however the glutamate transporter, GLT-1, which is also expressed in astrocytes was not regulated by the lesion. Conclusions The results confirm that riluzole is a neuroprotective agent in a rodent model of parkinson’s disease. Riluzole administration did not regulate GLT-1 levels but significantly reduced GFAP levels, in the lesioned striatum. Riluzole suppression of reactive astrocytosis is an intriguing finding which might contribute to the neuroprotective effects of this drug.

Age related changes in speed and mechanism of adult skeletal muscle stem cell migration

Skeletal muscle undergoes a progressive age-related loss in mass and function. Preservation of muscle mass depends in part on satellite cells, the resident stem cells of skeletal muscle. Reduced satellite cell function may contribute to the age-associated decrease in muscle mass. Here we focused on characterising the effect of age on satellite cell migration. We report that aged satellite cells migrate at less than half the speed of young cells. In addition, aged cells show abnormal membrane extension and retraction characteristics required for amoeboid based cell migration. Aged satellite cells displayed low levels of integrin expression. By deploying a mathematical model approach to investigate mechanism of migration, we have found that young satellite cells move in a random ‘memoryless’ manner whereas old cells demonstrate superdiffusive tendencies. Most importantly, we show that nitric oxide, a key regulator of cell migration, reversed the loss in migration speed and reinstated the unbiased mechanism of movement in aged satellite cells. Finally we found that although Hepatocyte Growth Factor increased the rate of aged satellite cell movement it did not restore the memoryless migration characteristics displayed in young cells. Our study shows that satellite cell migration, a key component of skeletal muscle regeneration, is compromised during aging. However, we propose clinically approved drugs could be used to overcome these detrimental changes.

Influence of substrate on corneal epithelial cell viability within ocular surface models

Corneal tissue engineering has improved dramatically over recent years. It is now possible to apply these technological advancements to the development of superior in vitro ocular surface models to reduce animal testing. We aim to show the effect different substrates can have on the viability of expanded corneal epithelial cells and that those which more accurately mimic the stromal surface provide the most protection against toxic assault. Compressed collagen gel as a substrate for the expansion of a human epithelial cell line was compared against two well-known substrates for modeling the ocular surface (polycarbonate membrane and conventional collagen gel). Cells were expanded over 10 days at which point cell stratification, cell number and expression of junctional proteins were assessed by electron microscopy, immunohistochemistry and RT-PCR. The effect of increasing concentrations of sodium lauryl sulphate on epithelial cell viability was quantified by MTT assay. Results showed improvement in terms of stratification, cell number and tight junction expression in human epithelial cells expanded upon either the polycarbonate membrane or compressed collagen gel when compared to a the use of a conventional collagen gel. However, cell viability was significantly higher in cells expanded upon the compressed collagen gel. We conclude that the more naturalistic composition and mechanical properties of compressed collagen gels produces a more robust corneal model.

Galleria mellonella as an infection model for Campylobacter jejuni virulence

Larvae of Galleria mellonella (Greater Wax Moth) have been shown to be susceptible to Campylobacter jejuni infection and our study characterizes this infection model. Following infection with C. jejuni human isolates, bacteria were visible in the haemocoel and gut of challenged larvae, and there was extensive damage to the gut. Bacteria were found in the extracellular and cell-associated fraction in the haemocoel, and it was shown that C. jejuni can survive in insect cells. Finally, we have used the model to screen a further 67 C. jejuni isolates belonging to different MLST types. Isolates belonging to ST257 were the most virulent in the Galleria model, whereas those belonging to ST21 were the least virulent.

Ensemble forecasts of a flood-producing storm: comparison of the influence of model-state perturbations and parameter modifications

High-resolution ensemble simulations (Δx = 1 km) are performed with the Met Office Unified Model for the Boscastle (Cornwall, UK) flash-flooding event of 16 August 2004. Forecast uncertainties arising from imperfections in the forecast model are analysed by comparing the simulation results produced by two types of perturbation strategy. Motivated by the meteorology of the event, one type of perturbation alters relevant physics choices or parameter settings in the model's parametrization schemes. The other type of perturbation is designed to account for representativity error in the boundary-layer parametrization. It makes direct changes to the model state and provides a lower bound against which to judge the spread produced by other uncertainties. The Boscastle has genuine skill at scales of approximately 60 km and an ensemble spread which can be estimated to within ∼ 10% with only eight members. Differences between the model-state perturbation and physics modification strategies are discussed, the former being more important for triggering and the latter for subsequent cell development, including the average internal structure of convective cells. Despite such differences, the spread in rainfall evaluated at skilful scales is shown to be only weakly sensitive to the perturbation strategy. This suggests that relatively simple strategies for treating model uncertainty may be sufficient for practical, convective-scale ensemble forecasting.

Effects of fermentation products of pro- and prebiotics on trans-epithelial electrical resistance in an in vitro model of the colon

Evidence from in vivo and in vitro studies suggests that the consumption of pro- and prebiotics may inhibit colon carcinogenesis; however, the mechanisms involved have, thus far, proved elusive. There are some indications from animal studies that the effects are being exerted during the promotion stage of carcinogenesis. One feature of the promotion stage of colorectal cancer is the disruption of tight junctions, leading to a loss of integrity across the intestinal barrier. We have used the Caco-2 human adenocarcinoma cell line as a model for the intestinal epithelia. Trans-epithelial electrical resistance measurements indicate Caco-2 monolayer integrity, and we recorded changes to this integrity following exposure to the fermentation products of selected probiotics and prebiotics, in the form of nondigestible oligosaccharides (NDOs). Our results indicate that NDOs themselves exert varying, but generally minor, effects upon the strength of the tight junctions, whereas the fermentation products of probiotics and NDOs tend to raise tight junction integrity above that of the controls. This effect was bacterial species and oligosaccharide specific. Bifidobacterium Bb 12 was particularly effective, as were the fermentation products of Raftiline and Raftilose. We further investigated the ability of Raftilose fermentations to protect against the negative effects of deoxycholic acid (DCA) upon tight junction integrity. We found protection to be species dependent and dependent upon the presence of the fermentation products in the media at the same time as or after exposure to the DCA. Results suggest that the Raftilose fermentation products may prevent disruption of the intestinal epithelial barrier function during damage by tumor promoters.

In vitro assessment of the bioaccessibility of brominated flame retardants in indoor dust using a colon extended model of the human gastrointestinal tract

An in vitro colon extended physiologically based extraction test (CEPBET) which incorporates human gastrointestinal tract (GIT) parameters (including pH and chemistry, solid-to-fluid ratio, mixing and emptying rates) was applied for the first time to study the bioaccessibility of brominated flame retardants (BFRs) from the 3 main GIT compartments (stomach, small intestine and colon) following ingestion of indoor dust. Results revealed the bioaccessibility of γ-HBCD (72%) was less than that for α- and β-isomers (92% and 80% respectively) which may be attributed to the lower aqueous solubility of the γ-isomer (2 μg L−1) compared to the α- and β-isomers (45 and 15 μg L−1 respectively). No significant change in the enantiomeric fractions of HBCDs was observed in any of the studied samples. However, this does not completely exclude the possibility of in vivo enantioselective absorption of HBCDs, as the GIT cell lining and bacterial flora – which may act enantioselectively – are not included in the current CE-PBET model. While TBBP-A was almost completely (94%) bioaccessible, BDE-209 was the least (14%) bioaccessible of the studied BFRs. Bioaccessibility of tri-hepta BDEs ranged from 32–58%. No decrease in the bioaccessibility with increasing level of bromination was observed in the studied PBDEs.

Adaptation-induced collective dynamics of a single-cell protozoan

We investigate the behavior of a single-cell protozoan in a narrow tubular ring. This environment forces them to swim under a one-dimensional periodic boundary condition. Above a critical density, single-cell protozoa aggregate spontaneously without external stimulation. The high-density zone of swimming cells exhibits a characteristic collective dynamics including translation and boundary fluctuation. We analyzed the velocity distribution and turn rate of swimming cells and found that the regulation of the turing rate leads to a stable aggregation and that acceleration of velocity triggers instability of aggregation. These two opposing effects may help to explain the spontaneous dynamics of collective behavior. We also propose a stochastic model for the mechanism underlying the collective behavior of swimming cells.

Characteristics of high-resolution versions of the Met Office unified model for forecasting convection over the United Kingdom

With many operational centers moving toward order 1-km-gridlength models for routine weather forecasting, this paper presents a systematic investigation of the properties of high-resolution versions of the Met Office Unified Model for short-range forecasting of convective rainfall events. The authors describe a suite of configurations of the Met Office Unified Model running with grid lengths of 12, 4, and 1 km and analyze results from these models for a number of convective cases from the summers of 2003, 2004, and 2005. The analysis includes subjective evaluation of the rainfall fields and comparisons of rainfall amounts, initiation, cell statistics, and a scale-selective verification technique. It is shown that the 4- and 1-km-gridlength models often give more realistic-looking precipitation fields because convection is represented explicitly rather than parameterized. However, the 4-km model representation suffers from large convective cells and delayed initiation because the grid length is too long to correctly reproduce the convection explicitly. These problems are not as evident in the 1-km model, although it does suffer from too numerous small cells in some situations. Both the 4- and 1-km models suffer from poor representation at the start of the forecast in the period when the high-resolution detail is spinning up from the lower-resolution (12 km) starting data used. A scale-selective precipitation verification technique implies that for later times in the forecasts (after the spinup period) the 1-km model performs better than the 12- and 4-km models for lower rainfall thresholds. For higher thresholds the 4-km model scores almost as well as the 1-km model, and both do better than the 12-km model.

Hurricane-type vortices in a general circulation model

A very high resolution atmospheric general circulation model, T106-L19, has been used for the simulation of hurricanes in a multi-year numerical experiment. Individual storms as well as their geographical and seasonal distribution agree remarkably well with observations. In spite of the fact that only the thermal and dynamical structure of the storms have been used as criteria of their identification, practically all of them occur in areas where the sea surface temperature is higher or equal to 26 °C. There are considerable variations from year to year in the number of storms in spite of the fact that there are no interannual variations in the SST pattern. It is found that the number of storms in particular areas appear to depend on the intensity of the Hadley-Walker cell. The result is clearly resolution-dependant. At lower horizonal resolution, T42, for example, the intensity of the storms is significantly reduced and their overall structure is less realistic, including their vertical form and extent.

Response of the middle atmosphere to CO2 doubling: results from the Canadian Middle Atmosphere Model

The Canadian Middle Atmosphere Model (CMAM) has been used to examine the middle atmosphere response to CO2 doubling. The radiative-photochemical response induced by doubling CO2 alone and the response produced by changes in prescribed SSTs are found to be approximately additive, with the former effect dominating throughout the middle atmosphere. The paper discusses the overall response, with emphasis on the effects of SST changes, which allow a tropospheric response to the CO2 forcing. The overall response is a cooling of the middle atmosphere accompanied by significant increases in the ozone and water vapor abundances. The ozone radiative feedback occurs through both an increase in solar heating and a decrease in infrared cooling, with the latter accounting for up to 15% of the total effect. Changes in global mean water vapor cooling are negligible above ~30 hPa. Near the polar summer mesopause, the temperature response is weak and not statistically significant. The main effects of SST changes are a warmer troposphere, a warmer and higher tropopause, cell-like structures of heating and cooling at low and middlelatitudes in the middle atmosphere, warming in the summer mesosphere, water vapor increase throughout the domain, and O3 decrease in the lower tropical stratosphere. No noticeable change in upwardpropagating planetary wave activity in the extratropical winter–spring stratosphere and no significant temperature response in the polar winter–spring stratosphere have been detected. Increased upwelling in the tropical stratosphere has been found to be linked to changed wave driving at low latitudes.