Stochastic simulation for spatial modelling of dynamic process in a living cell

One of the fundamental motivations underlying computational cell biology is to gain insight into the complicated dynamical processes taking place, for example, on the plasma membrane or in the cytosol of a cell. These processes are often so complicated that purely temporal mathematical models cannot adequately capture the complex chemical kinetics and transport processes of, for example, proteins or vesicles. On the other hand, spatial models such as Monte Carlo approaches can have very large computational overheads. This chapter gives an overview of the state of the art in the development of stochastic simulation techniques for the spatial modelling of dynamic processes in a living cell.

Linear models for endocytic transformations from live cell imaging

Endocytosis is the process by which cells internalise molecules including nutrient proteins from the extracellular media. In one form, macropinocytosis, the membrane at the cell surface ruffles and folds over to give rise to an internalised vesicle. Negatively charged phospholipids within the membrane called phosphoinositides then undergo a series of transformations that are critical for the correct trafficking of the vesicle within the cell, and which are often pirated by pathogens such as Salmonella. Advanced fluorescent video microscopy imaging now allows the detailed observation and quantification of these events in live cells over time. Here we use these observations as a basis for building differential equation models of the transformations. An initial investigation of these interactions was modelled with reaction rates proportional to the sum of the concentrations of the individual constituents. A first order linear system for the concentrations results. The structure of the system enables analytical expressions to be obtained and the problem becomes one of determining the reaction rates which generate the observed data plots. We present results with reaction rates which capture the general behaviour of the reactions so that we now have a complete mathematical model of phosphoinositide transformations that fits the experimental observations. Some excellent fits are obtained with modulated exponential functions; however, these are not solutions of the linear system. The question arises as to how the model may be modified to obtain a system whose solution provides a more accurate fit.

Phenomenological modeling of cell-to-cell and beat-to-beat variability in isolated Guinea Pig ventricular myocytes

Experimental action potential (AP) recordings in isolated ventricular myoctes display significant temporal beat-to-beat variability in morphology and duration. Furthermore, significant cell-to-cell differences in AP also exist even for isolated cells originating from the same region of the same heart. However, current mathematical models of ventricular AP fail to replicate the temporal and cell-to-cell variability in AP observed experimentally. In this study, we propose a novel mathematical framework for the development of phenomenological AP models capable of capturing cell-to-cell and temporal variabilty in cardiac APs. A novel stochastic phenomenological model of the AP is developed, based on the deterministic Bueno-Orovio/Fentonmodel. Experimental recordings of AP are fit to the model to produce AP models of individual cells from the apex and the base of the guinea-pig ventricles. Our results show that the phenomenological model is able to capture the considerable differences in AP recorded from isolated cells originating from the location. We demonstrate the closeness of fit to the available experimental data which may be achieved using a phenomenological model, and also demonstrate the ability of the stochastic form of the model to capture the observed beat-to-beat variablity in action potential duration.

An ultraviolet selective photodetector based on a nanocrystalline TiO2 photoelectrochemical cell

Solar ultraviolet (UV) radiation causes a range of skin disorders as well as affecting vision and the immune system. It also inhibits development of plants and animals. UV radiation monitoring is used routinely in some locations in order to alert the population to harmful solar radiation levels. There is ongoing research to develop UV-selective-sensors [1–3]. A personal, inexpensive and simple UV-selective-sensor would be desirable to measure UV intensity exposure. A prototype of such a detector has been developed and evaluated in our laboratory. It comprises a sealed two-electrode photoelectrochemical cell (PEC) based on nanocrystalline TiO2. This abundant semiconducting oxide, which is innocuous and very sta-ble, is the subject of intense study at present due to its application in dye sensitized solar cells (DSSC) [4]. Since TiO2 has a wide band gap (EG = 3.0 eV for rutile and EG = 3.2 eV for anatase), it is inher-ently UV-selective, so that UV filters are not required. This further reduces the cost of the proposed photodetector in comparison with conventional silicon detectors. The PEC is a semiconductor–electrolyte device that generates a photovoltage when it is illuminated and a corresponding photocur-rent if the external circuit is closed. The device does not require external bias, and the short circuit current is generally a linear function of illumination intensity. This greatly simplifies the elec-trical circuit needed when using the PEC as a photodetector. DSSC technology, which is based on a PEC containing nanocrystalline TiO2 sensitized with a ruthenium dye, holds out the promise of solar cells that are significantly cheaper than traditional silicon solar cells. The UV-sensor proposed in this paper relies on the cre-ation of electron–hole pairs in the TiO2 by UV radiation, so that it would be even cheaper than a DSSC since no sensitizer dye is needed. Although TiO2 has been reported as a suitable material for UV sensing [3], to the best of our knowledge, the PEC configuration described in the present paper is a new approach. In the present study, a novel double-layer TiO2 structure has been investigated. Fabrication is based on a simple and inexpensive technique for nanostructured TiO2 deposition using microwave-activated chemical bath deposition (MW-CBD) that has been reported recently [5]. The highly transparent TiO2 (anatase) films obtained are densely packed, and they adhere very well to the transparent oxide (TCO) substrate [6]. These compact layers have been studied as contacting layers in double-layer TiO2 structures for DSSC since improvement of electron extraction at the TiO2–TCO interface is expected [7]. Here we compare devices incorporating a single mesoporous nanocrystalline TiO2 structure with devices based on a double structure in which a MW-CBD film is situated between the TCO and the mesoporous nanocrystalline TiO2 layer. Besides improving electron extraction, this film could also help to block recombination of electrons transferred to the TCO with oxidized species in the electrolyte, as has been reported in the case of DSSC for compact TiO2 films obtained by other deposition tech-niques [8,9]. The two types of UV-selective sensors were characterized in detail. The current voltage characteristics, spectral response, inten-sity dependence of short circuit current and response times were measured and analyzed in order to evaluate the potential of sealed mesoporous TiO2-based photoelectrochemical cells (PEC) as low cost personal UV-photodetectors.

Effect of cold water immersion on repeated cycling performance and limb blood flow

The purpose of the present study was to compare the effects of cold water immersion (CWI) and active recovery (ACT) on resting limb blood flow, rectal temperature and repeated cycling performance in the heat. Ten subjects completed two testing sessions separated by 1 week; each trial consisted of an initial all-out 35-min exercise bout, one of two 15-min recovery interventions (randomised: CWI or ACT), followed by a 40-min passive recovery period before repeating the 35-min exercise bout. Performance was measured as the change in total work completed during the exercise bouts. Resting limb blood flow, heart rate, rectal temperature and blood lactate were recorded throughout the testing sessions. There was a significant decline in performance after ACT (mean (SD) −1.81% (1.05%)) compared with CWI where performance remained unchanged (0.10% (0.71%)). Rectal temperature was reduced after CWI (36.8°C (1.0°C)) compared with ACT (38.3°C (0.4°C)), as was blood flow to the arms (CWI 3.64 (1.47) ml/100 ml/min; ACT 16.85 (3.57) ml/100 ml/min) and legs (CW 4.83 (2.49) ml/100 ml/min; ACT 4.83 (2.49) ml/100 ml/min). Leg blood flow at the end of the second exercise bout was not different between the active (15.25 (4.33) ml/100 ml/min) and cold trials (14.99 (4.96) ml/100 ml/min), whereas rectal temperature (CWI 38.1°C (0.3°C); ACT 38.8°C (0.2°C)) and arm blood flow (CWI 20.55 (3.78) ml/100 ml/min; ACT 23.83 (5.32) ml/100 ml/min) remained depressed until the end of the cold trial. These findings indicate that CWI is an effective intervention for maintaining repeat cycling performance in the heat and this performance benefit is associated with alterations in core temperature and limb blood flow.

Oxygen delivery through high-flow nasal cannulae increase end-expiratory lung volume and reduce respiratory rate in post-cardiac surgical patients

Background: High-flow nasal cannulae (HFNC) create positive oropharyngeal airway pressure but it is unclear how their use affects lung volume. Electrical impedance tomography (EIT) allows assessment of changes in lung volume by measuring changes in lung impedance. Primary objectives were to investigate the effects of HFNC on airway pressure (Paw) and end-expiratory lung volume (EELV), and to identify any correlation between the two. Secondary objectives were to investigate the effects of HFNC on respiratory rate (RR), dyspnoea, tidal volume and oxygenation; and the interaction between body mass index (BMI) and EELV. Methods: Twenty patients prescribed HFNC post-cardiac surgery were investigated. Impedance measures, Paw, PaO2/FiO2 ratio, RR and modified Borg scores were recorded first on low flow oxygen (nasal cannula or Hudson face mask) and then on HFNC. Results: A strong and significant correlation existed between Paw and end-expiratory lung impedance (EELI) (r=0.7, p<0.001). Compared with low flow oxygen, HFNC significantly increased EELI by 25.6% (95% CI 24.3, 26.9) and Paw by 3.0 cmH2O (95% CI 2.4, 3.7). RR reduced by 3.4 breaths per minute (95% CI 1.7, 5.2) with HFNC use, tidal impedance variation increased by 10.5% (95% CI 6.1, 18.3) and PaO2/FiO2 ratio improved by 30.6 mmHg (95% CI 17.9, 43.3). HFNC improved subjective dyspnoea scoring (p=0.023). Increases in EELI were significantly influenced by BMI, with larger increases associated with higher BMIs (p<0.001). Conclusions: This study suggests that HFNC improve dyspnoea and oxygenation by increasing both EELV and tidal volume, and are most beneficial in patients with higher BMIs.

An analytical solution for diffusion and nonlinear uptake of oxygen in a spherical cell

We develop a new analytical solution for a reactive transport model that describes the steady-state distribution of oxygen subject to diffusive transport and nonlinear uptake in a sphere. This model was originally reported by Lin (Journal of Theoretical Biology, 1976 v60, pp449–457) to represent the distribution of oxygen inside a cell and has since been studied extensively by both the numerical analysis and formal analysis communities. Here we extend these previous studies by deriving an analytical solution to a generalized reaction-diffusion equation that encompasses Lin’s model as a particular case. We evaluate the solution for the parameter combinations presented by Lin and show that the new solutions are identical to a grid-independent numerical approximation.

Natural convection flow in a porous enclosure with localized heating from below

Natural convection flow in a two-dimensional fluid saturated porous enclosure with localized heating from below, symmetrical cooling from the sides and the top and rest of the bottom walls are insulated, has been investigated numerically. Darcy’s law for porous media along with the energy equation based on the 1st law of thermodynamics has been considered. Implicit finite volume method with TDMA solver is used to solve the governing equations. Localized heating is simulated by a centrally located isothermal heat source on the bottom wall, and four different values of the dimensionless heat source length, 1/5, 2/5, 3/5 and 4/5 are considered. The effect of heat source length and the Rayleigh number on streamlines and isotherms are presented, as well as the variation of the local rate of heat transfer in terms of the local Nusselt number from the heated wall. Finally, the average Nusselt number at the heated part of the bottom wall has been shown against Rayleigh number for the non-dimensional heat source length.

Data flow analysis of embedded program expressions

Data flow analysis techniques can be used to help assess threats to data confidentiality and integrity in security critical program code. However, a fundamental weakness of static analysis techniques is that they overestimate the ways in which data may propagate at run time. Discounting large numbers of these false-positive data flow paths wastes an information security evaluator's time and effort. Here we show how to automatically eliminate some false-positive data flow paths by precisely modelling how classified data is blocked by certain expressions in embedded C code. We present a library of detailed data flow models of individual expression elements and an algorithm for introducing these components into conventional data flow graphs. The resulting models can be used to accurately trace byte-level or even bit-level data flow through expressions that are normally treated as atomic. This allows us to identify expressions that safely downgrade their classified inputs and thereby eliminate false-positive data flow paths from the security evaluation process. To validate the approach we have implemented and tested it in an existing data flow analysis toolkit.

Melanoma cell invasiveness is regulated by miR-211 suppression of the BRN2 transcription factor

To identify microRNAs potentially involved in melanomagenesis, we compared microRNA expression profiles between melanoma cell lines and cultured melanocytes. The most differentially expressed microRNA between the normal and tumor cell lines was miR-211. We focused on this pigment-cell-enriched miRNA as it is derived from the microphthalmia-associated transcription factor (MITF)-regulated gene, TRPM1 (melastatin). We find that miR-211 expression is greatly decreased in melanoma cells and melanoblasts compared to melanocytes. Bioinformatic analysis identified a large number of potential targets of miR-211, including POU3F2 (BRN2). Inhibition of miR-211 in normal melanocytes resulted in increased BRN2 protein, indicating that endogenous miR-211 represses BRN2 in differentiated cells. Over-expression of miR-211 in melanoma cell lines changed the invasive potential of the cells in vitro through directly targeting BRN2 translation. We propose a model for the apparent non-overlapping expression levels of BRN2 and MITF in melanoma, mediated by miR-211 expression.

Human neural cell interactions with orientated electrospun nanofibers in vitro

Aim: Electrospun nanofibers represent potent guidance substrates for nervous tissue repair. Development of nanofiber-based scaffolds for CNS repair requires, as a first step, an understanding of appropriate neural cell type-substrate interactions. Materials & methods: Astrocyte–nanofiber interactions (e.g., adhesion, proliferation, process extension and migration) were studied by comparing human neural progenitor-derived astrocytes (hNP-ACs) and a human astrocytoma cell line (U373) with aligned polycaprolactone (PCL) nanofibers or blended (25% type I collagen/75% PCL) nanofibers. Neuron–nanofiber interactions were assessed using a differentiated human neuroblastoma cell line (SH-SY5Y). Results & discussion: U373 cells and hNP-AC showed similar process alignment and length when associated with PCL or Type I collagen/PCL nanofibers. Cell adhesion and migration by hNP-AC were clearly improved by functionalization of nanofiber surfaces with type I collagen. Functionalized nanofibers had no such effect on U373 cells. Another clear difference between the U373 cells and hNP-AC interactions with the nanofiber substrate was proliferation; the cell line demonstrating strong proliferation, whereas the hNP-AC line showed no proliferation on either type of nanofiber. Long axonal growth (up to 600 µm in length) of SH-SY5Y neurons followed the orientation of both types of nanofibers even though adhesion of the processes to the fibers was poor. Conclusion: The use of cell lines is of only limited predictive value when studying cell–substrate interactions but both morphology and alignment of human astrocytes were affected profoundly by nanofibers. Nanofiber surface functionalization with collagen significantly improved hNP-AC adhesion and migration. Alternative forms of functionalization may be required for optimal axon–nanofiber interactions.

Cabergoline Decreases Alcohol Drinking and Seeking Behaviors Via Glial Cell Line-Derived Neurotrophic Factor

Background: Cabergoline is an ergotamine derivative that increases the expression of glial cell line-derived neurotrophic factor (GDNF) in vitro. We recently showed that GDNF in the ventral tegmental area (VTA) reduces the motivation to consume alcohol. We therefore set out to determine whether cabergoline administration decreases alcohol-drinking and -seeking behaviors via GDNF. Methods: Reverse transcription polymerase chain reaction (RT-PCR) and Enzyme-Linked ImmunoSorbent Assay (ELISA) were used to measure GDNF levels. Western blot analysis was used for phosphorylation experiments. Operant self-administration in rats and a two-bottle choice procedure in mice were used to assess alcohol-drinking behaviors. Instrumental performance tested during extinction was used to measure alcohol-seeking behavior. The [35S]GTPγS binding assay was used to assess the expression and function of the dopamine D2 receptor (D2R). Results: We found that treatment of the dopaminergic-like cell line SH-SY5Y with cabergoline and systemic administration of cabergoline in rats resulted in an increase in GDNF level and in the activation of the GDNF pathway. Cabergoline treatment decreased alcohol-drinking and -seeking behaviors including relapse, and its action to reduce alcohol consumption was localized to the VTA. Finally, the increase in GDNF expression and the decrease in alcohol consumption by cabergoline were abolished in GDNF heterozygous knockout mice. Conclusions: Together, these findings suggest that cabergoline-mediated upregulation of the GDNF pathway attenuates alcohol-drinking behaviors and relapse. Alcohol abuse and addiction are devastating and costly problems worldwide. This study puts forward the possibility that cabergoline might be an effective treatment for these disorders. © 2009 Society of Biological Psychiatry.

Mixed convection laminar flow along a horizontal plate subject to streamwise sinusoidal surface temperature

Mixed convection of a two-dimensional laminar incompressible flow along a horizontal flat plate with streamwise sinusoidal surface temperature has been numerically investigated for different values of Rayleigh number and Reynolds number for constant values of Prandtl number, amplitude and frequency of periodic temperature. The numerical scheme is based on the finite element method adapted to rectangular non-uniform mesh elements by a non-linear parametric solution algorithm. The fluid considered in this study is air. The results are obtained for the Rayleigh number and Reynolds number ranging from 102 to 104 and 1 to 100, respectively, with constant physical properties for the fluid medium considered. Velocity and temperature profiles, streamlines, isotherms, and average Nusselt numbers are presented to observe the effect of the investigating parameters on fluid flow and heat transfer characteristics. The present results show that the convective phenomena are greatly influenced by the variation of Rayleigh numbers and Reynolds number.

Flow regime transition criteria for two-phase flow at reduced gravity conditions

Flow regime transition criteria are of practical importance for two-phase flow analyses at reduced gravity conditions. Here, flow regime transition criteria which take the friction pressure loss effect into account were studied in detail. Criteria at reduced gravity conditions were developed by extending an existing model with various experimental datasets taken at microgravity conditions showed satisfactory agreement. Sample computations of the model were performed at various gravity conditions, such as 0.196, 1.62, 3.71, and 9.81 m/s2 corresponding to micro-gravity and lunar, Martian and Earth surface gravity, respectively. It was found that the effect of gravity on bubbly-slug and slug-annular (churn) transitions in a two-phase flow system was more pronounced at low liquid flow conditions, whereas the gravity effect could be ignored at high mixture volumetric flux conditions. While for the annular flow transitions due to flow reversal and onset of dropset entrainment, higher superficial gas velocity was obtained at higher gravity level.