Efficient on the fly calculation of time correlation functions in computer simulations

Time correlation functions yield profound information about the dynamics of a physical system and hence are frequently calculated in computer simulations. For systems whose dynamics span a wide range of time, currently used methods require significant computer time and memory. In this paper, we discuss the multiple-tau correlator method for the efficient calculation of accurate time correlation functions on the fly during computer simulations. The multiple-tau correlator is efficacious in terms of computational requirements and can be tuned to the desired level of accuracy. Further, we derive estimates for the error arising from the use of the multiple-tau correlator and extend it for use in the calculation of mean-square particle displacements and dynamic structure factors. The method described here, in hardware implementation, is routinely used in light scattering experiments but has not yet found widespread use in computer simulations.

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.

Generation of candidate human influenza vaccine strains in cell culture: rehearsing the European response to an H7N1 pandemic threat

Background: Although H5N1 avian influenza viruses pose the most obvious imminent pandemic threat, there have been several recent zoonotic incidents involving transmission of H7 viruses to humans. Vaccines are the primary public health defense against pandemics, but reliance on embryonated chickens eggs to propagate vaccine and logistic problems posed by the use of new technology may slow our ability to respond rapidly in a pandemic situation. Objectives: We sought to generate an H7 candidate vaccine virus suitable for administration to humans whose generation and amplification avoided the use of eggs. Methods: We generated a suitable H7 vaccine virus by reverse genetics. This virus, known as RD3, comprises the internal genes of A/Puerto Rico/8/34 with surface antigens of the highly pathogenic avian strain A/Chicken/Italy/13474/99 (H7N1). The multi-basic amino acid site in the HA gene, associated with high pathogenicity in chickens, was removed. Results: The HA modification did not alter the antigenicity of the virus and the resultant single basic motif was stably retained following several passages in Vero and PER. C6 cells. RD3 was attenuated for growth in embryonated eggs, chickens, and ferrets. RD3 induced an antibody response in infected animals reactive against both the homologous virus and other H7 influenza viruses associated with recent infection by H7 viruses in humans. Conclusions: This is the first report of a candidate H7 vaccine virus for use in humans generated by reverse genetics and propagated entirely in mammalian tissue culture. The vaccine has potential use against a wide range of H7 strains.

Class II phosphoinositide 3-kinase defines a novel signaling pathway in cell migration

The lipid products of phosphoinositide 3-kinase (PI3K) are involved in many cellular responses such as proliferation, migration, and survival. Disregulation of PI3K-activated pathways is implicated in different diseases including cancer and diabetes. Among the three classes of PI3Ks, class I is the best characterized, whereas class II has received increasing attention only recently and the precise role of these isoforms is unclear. Similarly, the role of phosphatidylinositol-3-phosphate (PtdIns-3-P) as an intracellular second messenger is only just beginning to be appreciated. Here, we show that lysophosphatidic acid (LPA) stimulates the production of PtdIns-3-P through activation of a class II PI3K (PI3K-C2β). Both PtdIns-3-P and PI3K-C2β are involved in LPA-mediated cell migration. This study is the first identification of PtdIns-3-P and PI3K-C2β as downstream effectors in LPA signaling and demonstration of an intracellular role for a class II PI3K. Defining this novel PI3K-C2β- PtdIns-3-P signaling pathway may help clarify the process of cell migration and may shed new light on PI3K-mediated intracellular events.

Functionalised amyloid fibrils for roles in cell adhesion

We describe experiments designed to explore the possibility of using amyloid fibrils as new nanoscale biomaterials for promoting and exploiting cell adhesion, migration and differentiation in vitro. We created peptides that add the biological cell adhesion sequence (RGD) or a control sequence (RAD) to the C-terminus of an 11-residue peptide corresponding to residues 105-115 of the amyloidogenic protein transthyretin. These peptides readily self-assemble in aqueous solution to form amyloid fibrils, and X-ray fibre diffraction shows that they possess the same strand and sheet spacing in the characteristic cross-beta structure as do fibrils formed by the parent peptide. We report that the fibrils containing the RGD sequence are bioactive and that these fibrils interact specifically with cells via the RGD group displayed on the fibril surface. As the design of such functionalized fibrils can be systematically altered, these findings suggest that it will be possible to generate nanomaterials based on amyloid fibrils that are tailored to promote interactions with a wide variety of cell types. (c) 2007 Elsevier Ltd. All rights reserved.

Differences in cell morphology, lipid and apo B secretory capacity in caco-2 cells following long term treatment with saturated and monounsaturated fatty acids

The suitability of the caco-2 cell line as a model for studying the long term impact of dietary fatty acids on intestinal lipid handling and chylomicron production was examined. Chronic supplementation of caco-2 cells with palmitic acid (PA) resulted in a lower triacylglycerol secretion than oleic acid (OA). This was coupled with a detrimental effect of PA, but not OA, on transepithelial electrical resistance (TER) measurements, suggesting a loss of structural integrity across the cell monolayer. Addition of OA reversed the adverse effects of PA and stearic acid on TER and increased the ability of cells to synthesise and accumulate lipid, but did not normalise the secretion of lipids by caco-2 cells. Increasing amounts of OA and decreasing amounts of PA in the incubation media markedly improved the ability of cells to synthesise apolipoprotein B and secrete lipids. Real time RT-PCR revealed a down regulation of genes involved in lipoprotein synthesis following PA than OA. Electron microscopy showed adverse effects of PA on cellular morphology consistent with immature enterocytes such as stunted microvilli and poor tight junction formation. In conclusion, previously reported differences in lipoprotein secretion by caco-2 cells supplemented with saturated fatty acids (SFA) and OA may partly reflect early cytotoxic effects of SFA on cellular integrity and function. (C) 2007 Elsevier B.V. All rights reserved.

The role of 80K/MARCKS, a specific substrate of protein kinase C, in cell growth and tumour progression

Since its discovery more than a decade ago [Wu et al., 1982; Rozengurt et al., 1983], the 80-87 kDa myristoylated a lanine-rich C-kinase substrate (80K/MARCKS) protein has attracted a great deal of attention from researchers interested in cell growth and tumour progression. However, despite its ubiquitous distribution, a definitive functional role for 80K/MARCKS has not been found. The purpose of this review is to describe the properties, distribution and regulation of 80K/MARCKS and to discuss some of the most recent findings, both from our laboratory and from others, that have suggested a functional role for this protein in modulating cell growth and tumour progression. Furthermore, I will present data from our laboratory that implicates 80K/MARCKS as a novel tumour suppressor in cells of melanocyte origin.

Pathogenic LRRK2 mutations do not alter gene expression in cell model systems or human brain tissue

Point mutations in LRRK2 cause autosomal dominant Parkinson's disease. Despite extensive efforts to determine the mechanism of cell death in patients with LRRK2 mutations, the aetiology of LRRK2 PD is not well understood. To examine possible alterations in gene expression linked to the presence of LRRK2 mutations, we carried out a case versus control analysis of global gene expression in three systems: fibroblasts isolated from LRRK2 mutation carriers and healthy, non-mutation carrying controls; brain tissue from G2019S mutation carriers and controls; and HEK293 inducible LRRK2 wild type and mutant cell lines. No significant alteration in gene expression was found in these systems following correction for multiple testing. These data suggest that any alterations in basal gene expression in fibroblasts or cell lines containing mutations in LRRK2 are likely to be quantitatively small. This work suggests that LRRK2 is unlikely to play a direct role in modulation of gene expression, although it remains possible that this protein can influence mRNA expression under pathogenic cicumstances.

Interaction between attaching and effacing Escherichia coli serotypes O157:H7 and O26:K60 in cell culture

Ruminants harbour both O157:H7 and non-O157 Attaching Effacing Escherichia coli (AEEC) strains but to date only nonO157 AEEC have been shown to induce attaching effacing lesions in naturally infected animals. However, O157 may induce lesions in deliberate oral inoculation studies and persistence is considered dependent upon the bacterially encoded locus for enterocyte effacement. In concurrent infections in ruminants it is unclear whether non-O157 AEEC contribute either positively or negatively to the persistence of E. coli O157:H7. To investigate this, and prior to animal studies, E. coli O157:H7 NCTC 12900, a non-toxigenic strain that persists in conventionally reared sheep, and non-toxigenic AEEC O26:K60 isolates of sheep origin were tested for adherence to Hep-2 tissue culture alone and in competition one with another. Applied together, both strains adhered in similar numbers but lower than when either was applied separately. Pre-incubation of tissue culture with either one strain reduced significantly (P < 0.05) the extent of adherence of the strain that was applied second. It was particularly noticeable that AEEC O26 when applied first reduced adherence and inhibited microcolony formation, as demonstrated by confocal microscopy, of E. coli 01 57:H7. The possibility that prior colonisation of a ruminant by non-O157 AEEC such as O26 may antagonise O157 colonisation and persistence in ruminants is discussed. (C) 2004 Elsevier B.V. All rights reserved.

A novel alternative to cryopreservation for the short-term storage of stem cells for use in cell therapy using alginate encapsulation

Efficient transport of stem/progenitor cells without affecting their survival and function is a key factor in any practical cell-based therapy. However, the current approach using liquid nitrogen for the transfer of stem cells requires a short delivery time window is technically challenging and financially expensive. The present study aims to use semipermeable alginate hydrogels (crosslinked by strontium) to encapsulate, store, and release stem cells, to replace the conventional cryopreservation method for the transport of therapeutic cells within world-wide distribution time frame. Human mesenchymal stem cell (hMSC) and mouse embryonic stem cells (mESCs) were successfully stored inside alginate hydrogels for 5 days under ambient conditions in an air-tight environment (sealed cryovial). Cell viability, of the cells extracted from alginate gel, gave 74% (mESC) and 80% (hMSC) survival rates, which compared favorably to cryopreservation. More importantly, the subsequent proliferation rate and detection of common stem cell markers (both in mRNA and protein level) from hMSCs and mESCs retrieved from alginate hydrogels were also comparable to (if not better than) results gained following cryopreservation. In conclusion, this new and simple application of alginate hydrogel encapsulation may offer a cheap and robust alternative to cryopreservation for the transport and storage of stem cells for both clinical and research purposes.

Parameterised diabatic processes in numerical simulations of an extratropical cyclone

The parameterisation of diabatic processes in numerical models is critical for the accuracy of weather forecasts and for climate projections. A novel approach to the evaluation of these processes in models is introduced in this contribution. The approach combines a suite of on-line tracer diagnostics with off-line trajectory calculations. Each tracer tracks accumulative changes in potential temperature associated with a particular parameterised diabatic process in the model. A comparison of tracers therefore allows the identification of the most active diabatic processes and their downstream impacts. The tracers are combined with trajectories computed using model-resolved winds, allowing the various diabatic contributions to be tracked back to their time and location of occurrence. We have used this approach to investigate diabatic processes within a simulated extratropical cyclone. We focus on the warm conveyor belt, in which the dominant diabatic contributions come from large-scale latent heating and parameterised convection. By contrasting two simulations, one with standard convection parameterisation settings and another with reduced parameterised convection, the effects of parameterised convection on the structure of the cyclone have been determined. Under reduced parameterised convection conditions, the large-scale latent heating is forced to release convective instability that would otherwise have been released by the convection parameterisation. Although the spatial distribution of precipitation depends on the details of the split between parameterised convection and large-scale latent heating, the total precipitation amount associated with the cyclone remains largely unchanged. For reduced parameterised convection, a more rapid and stronger latent heating episode takes place as air ascends within the warm conveyor belt.

Vertical structure and physical processes of the Madden-Julian oscillation: exploring key model physics in climate simulations

Aimed at reducing deficiencies in representing the Madden-Julian oscillation (MJO) in general circulation models (GCMs), a global model evaluation project on vertical structure and physical processes of the MJO was coordinated. In this paper, results from the climate simulation component of this project are reported. It is shown that the MJO remains a great challenge in these latest generation GCMs. The systematic eastward propagation of the MJO is only well simulated in about one-fourth of the total participating models. The observed vertical westward tilt with altitude of the MJO is well simulated in good MJO models, but not in the poor ones. Damped Kelvin wave responses to the east of convection in the lower troposphere could be responsible for the missing MJO preconditioning process in these poor MJO models. Several process-oriented diagnostics were conducted to discriminate key processes for realistic MJO simulations. While large-scale rainfall partition and low-level mean zonal winds over the Indo-Pacific in a model are not found to be closely associated with its MJO skill, two metrics, including the low-level relative humidity difference between high and low rain events and seasonal mean gross moist stability, exhibit statistically significant correlations with the MJO performance. It is further indicated that increased cloud-radiative feedback tends to be associated with reduced amplitude of intraseasonal variability, which is incompatible with the radiative instability theory previously proposed for the MJO. Results in this study confirm that inclusion of air-sea interaction can lead to significant improvement in simulating the MJO.

Analysis of convectively-generated gravity waves in mesoscale model simulations and wind profiler observations

The characteristics of convectively-generated gravity waves during an episode of deep convection near the coast of Wales are examined in both high resolution mesoscale simulations [with the (UK) Met Oce Unified Model] and in observations from a Mesosphere-Stratosphere-Troposphere (MST) wind profiling Doppler radar. Deep convection reached the tropopause and generated vertically propagating, high frequency waves in the lower stratosphere that produced vertical velocity perturbations O(1 m/s). Wavelet analysis is applied in order to determine the characteristic periods and wavelengths of the waves. In both the simulations and observations, the wavelet spectra contain several distinct preferred scales indicated by multiple spectral peaks. The peaks are most pronounced in the horizontal spectra at several wavelengths less than 50 km. Although these peaks are most clear and of largest amplitude in the highest resolution simulations (with 1 km horizontal grid length), they are also evident in coarser simulations (with 4 km horizontal grid length). Peaks also exist in the vertical and temporal spectra (between approximately 2.5 and 4.5 km, and 10 to 30 minutes, respectively) with good agreement between simulation and observation. Two-dimensional (wavenumber-frequency) spectra demonstrate that each of the selected horizontal scales contains peaks at each of preferred temporal scales revealed by the one- dimensional spectra alone.