927 resultados para jump kinetics
Resumo:
Delays are an important feature in temporal models of genetic regulation due to slow biochemical processes, such as transcription and translation. In this paper, we show how to model intrinsic noise effects in a delayed setting by either using a delay stochastic simulation algorithm (DSSA) or, for larger and more complex systems, a generalized Binomial τ-leap method (Bτ-DSSA). As a particular application, we apply these ideas to modeling somite segmentation in zebra fish across a number of cells in which two linked oscillatory genes (her1 and her7) are synchronized via Notch signaling between the cells.
Resumo:
Recently the application of the quasi-steady-state approximation (QSSA) to the stochastic simulation algorithm (SSA) was suggested for the purpose of speeding up stochastic simulations of chemical systems that involve both relatively fast and slow chemical reactions [Rao and Arkin, J. Chem. Phys. 118, 4999 (2003)] and further work has led to the nested and slow-scale SSA. Improved numerical efficiency is obtained by respecting the vastly different time scales characterizing the system and then by advancing only the slow reactions exactly, based on a suitable approximation to the fast reactions. We considerably extend these works by applying the QSSA to numerical methods for the direct solution of the chemical master equation (CME) and, in particular, to the finite state projection algorithm [Munsky and Khammash, J. Chem. Phys. 124, 044104 (2006)], in conjunction with Krylov methods. In addition, we point out some important connections to the literature on the (deterministic) total QSSA (tQSSA) and place the stochastic analogue of the QSSA within the more general framework of aggregation of Markov processes. We demonstrate the new methods on four examples: Michaelis–Menten enzyme kinetics, double phosphorylation, the Goldbeter–Koshland switch, and the mitogen activated protein kinase cascade. Overall, we report dramatic improvements by applying the tQSSA to the CME solver.
Resumo:
Biologists are increasingly conscious of the critical role that noise plays in cellular functions such as genetic regulation, often in connection with fluctuations in small numbers of key regulatory molecules. This has inspired the development of models that capture this fundamentally discrete and stochastic nature of cellular biology - most notably the Gillespie stochastic simulation algorithm (SSA). The SSA simulates a temporally homogeneous, discrete-state, continuous-time Markov process, and of course the corresponding probabilities and numbers of each molecular species must all remain positive. While accurately serving this purpose, the SSA can be computationally inefficient due to very small time stepping so faster approximations such as the Poisson and Binomial τ-leap methods have been suggested. This work places these leap methods in the context of numerical methods for the solution of stochastic differential equations (SDEs) driven by Poisson noise. This allows analogues of Euler-Maruyuma, Milstein and even higher order methods to be developed through the Itô-Taylor expansions as well as similar derivative-free Runge-Kutta approaches. Numerical results demonstrate that these novel methods compare favourably with existing techniques for simulating biochemical reactions by more accurately capturing crucial properties such as the mean and variance than existing methods.
Resumo:
Self-segregation and compartimentalisation are observed experimentally to occur spontaneously on live membranes as well as reconstructed model membranes. It is believed that many of these processes are caused or supported by anomalous diffusive behaviours of biomolecules on membranes due to the complex and heterogeneous nature of these environments. These phenomena are on the one hand of great interest in biology, since they may be an important way for biological systems to selectively localize receptors, regulate signaling or modulate kinetics; and on the other, they provide an inspiration for engineering designs that mimick natural systems. We present an interactive software package we are developing for the purpose of simulating such processes numerically using a fundamental Monte Carlo approach. This program includes the ability to simulate kinetics and mass transport in the presence of either mobile or immobile obstacles and other relevant structures such as liquid-ordered lipid microdomains. We also present preliminary simulation results regarding the selective spatial localization and chemical kinetics modulating power of immobile obstacles on the membrane, obtained using the program.
Resumo:
The performance and electron recombination kinetics of dye-sensitized solar cells based on TiO2 films consisting of one-dimensional nanorod arrays (NR-DSSCs) which are sensitized with dye N719, C218 and D205 respectively have been studied. It has been found that the best efficiency is obtained with the dye C218 based NR-DSSCs, benefiting from a 40% higher short-circuit photocurrent density. However, the open circuit photovoltage of the N719 based cell is 40 mV higher than that of the organic dye C218 and D205 based devices. Investigation of the electron recombination kinetics of the NR-DSSCs has revealed that the effective electron lifetime, τn, of the N719 based NR-DSSC is the lowest whereas the τn of the C218 based NR-DSSC is the highest among the three dyes. The higher Voc with the N719 based NR-DSSC is originated from the more negative energy level of the conduction band of the TiO2 film. In addition, in comparison to the DSSCs with conventional nanocrystalline particles based TiO2 films, the NR-DSSCs have shown over two orders of magnitude higher τn when employing N719 as the sensitizer. Nevertheless, the τn of the DSSCs with the C218 based nanorod arrays is only ten-fold higher than the that of the nanoparticles based devices. The remarkable characteristic of the dye C218 in suppressing the electron recombination of DSSCs is discussed.
Resumo:
This study of photocatalytic oxidation of phenol over titanium dioxide films presents a method for the evaluation of true reaction kinetics. A flat plate reactor was designed for the specific purpose of investigating the influence of various reaction parameters, specifically photocatalytic film thickness, solution flow rate (1–8 l min−1), phenol concentration (20, 40 and 80 ppm), and irradiation intensity (70.6, 57.9, 37.1and 20.4 W m−2), in order to further understand their impact on the reaction kinetics. Special attention was given to the mass transfer phenomena and the influence of film thickness. The kinetics of phenol degradation were investigated with different irradiation levels and initial pollutant concentration. Photocatalytic degradation experiments were performed to evaluate the influence of mass transfer on the reaction and, in addition, the benzoic acid method was applied for the evaluation of mass transfer coefficient. For this study the reactor was modelled as a batch-recycle reactor. A system of equations that accounts for irradiation, mass transfer and reaction rate was developed to describe the photocatalytic process, to fit the experimental data and to obtain kinetic parameters. The rate of phenol photocatalytic oxidation was described by a Langmuir–Hinshelwood type law that included competitive adsorption and degradation of phenol and its by-products. The by-products were modelled through their additive effect on the solution total organic carbon.
Resumo:
Velocity jump processes are discrete random walk models that have many applications including the study of biological and ecological collective motion. In particular, velocity jump models are often used to represent a type of persistent motion, known as a “run and tumble”, which is exhibited by some isolated bacteria cells. All previous velocity jump processes are non-interacting, which means that crowding effects and agent-to-agent interactions are neglected. By neglecting these agent-to-agent interactions, traditional velocity jump models are only applicable to very dilute systems. Our work is motivated by the fact that many applications in cell biology, such as wound healing, cancer invasion and development, often involve tissues that are densely packed with cells where cell-to-cell contact and crowding effects can be important. To describe these kinds of high cell density problems using a velocity jump process we introduce three different classes of crowding interactions into a one-dimensional model. Simulation data and averaging arguments lead to a suite of continuum descriptions of the interacting velocity jump processes. We show that the resulting systems of hyperbolic partial differential equations predict the mean behavior of the stochastic simulations very well.
Resumo:
Fibroblasts and their activated phenotype, myofibroblasts, are the primary cell types involved in the contraction associated with dermal wound healing. Recent experimental evidence indicates that the transformation from fibroblasts to myofibroblasts involves two distinct processes: the cells are stimulated to change phenotype by the combined actions of transforming growth factor β (TGFβ) and mechanical tension. This observation indicates a need for a detailed exploration of the effect of the strong interactions between the mechanical changes and growth factors in dermal wound healing. We review the experimental findings in detail and develop a model of dermal wound healing that incorporates these phenomena. Our model includes the interactions between TGFβ and collagenase, providing a more biologically realistic form for the growth factor kinetics than those included in previous mechanochemical descriptions. A comparison is made between the model predictions and experimental data on human dermal wound healing and all the essential features are well matched.
Resumo:
Even though titanium dioxide photocatalysis has been promoted as a leading green technology for water purification, many issues have hindered its application on a large commercial scale. For the materials scientist the main issues have centred the synthesis of more efficient materials and the investigation of degradation mechanisms; whereas for the engineers the main issues have been the development of appropriate models and the evaluation of intrinsic kinetics parameters that allow the scale up or re-design of efficient large-scale photocatalytic reactors. In order to obtain intrinsic kinetics parameters the reaction must be analysed and modelled considering the influence of the radiation field, pollutant concentrations and fluid dynamics. In this way, the obtained kinetic parameters are independent of the reactor size and configuration and can be subsequently used for scale-up purposes or for the development of entirely new reactor designs. This work investigates the intrinsic kinetics of phenol degradation over titania film due to the practicality of a fixed film configuration over a slurry. A flat plate reactor was designed in order to be able to control reaction parameters that include the UV irradiance, flow rates, pollutant concentration and temperature. Particular attention was paid to the investigation of the radiation field over the reactive surface and to the issue of mass transfer limited reactions. The ability of different emission models to describe the radiation field was investigated and compared to actinometric measurements. The RAD-LSI model was found to give the best predictions over the conditions tested. Mass transfer issues often limit fixed film reactors. The influence of this phenomenon was investigated with specifically planned sets of benzoic acid experiments and with the adoption of the stagnant film model. The phenol mass transfer coefficient in the system was calculated to be km,phenol=8.5815x10-7Re0.65(ms-1). The data obtained from a wide range of experimental conditions, together with an appropriate model of the system, has enabled determination of intrinsic kinetic parameters. The experiments were performed in four different irradiation levels (70.7, 57.9, 37.1 and 20.4 W m-2) and combined with three different initial phenol concentrations (20, 40 and 80 ppm) to give a wide range of final pollutant conversions (from 22% to 85%). The simple model adopted was able to fit the wide range of conditions with only four kinetic parameters; two reaction rate constants (one for phenol and one for the family of intermediates) and their corresponding adsorption constants. The intrinsic kinetic parameters values were defined as kph = 0.5226 mmol m-1 s-1 W-1, kI = 0.120 mmol m-1 s-1 W-1, Kph = 8.5 x 10-4 m3 mmol-1 and KI = 2.2 x 10-3 m3 mmol-1. The flat plate reactor allowed the investigation of the reaction under two different light configurations; liquid and substrate side illumination. The latter of particular interest for real world applications where light absorption due to turbidity and pollutants contained in the water stream to be treated could represent a significant issue. The two light configurations allowed the investigation of the effects of film thickness and the determination of the catalyst optimal thickness. The experimental investigation confirmed the predictions of a porous medium model developed to investigate the influence of diffusion, advection and photocatalytic phenomena inside the porous titania film, with the optimal thickness value individuated at 5 ìm. The model used the intrinsic kinetic parameters obtained from the flat plate reactor to predict the influence of thickness and transport phenomena on the final observed phenol conversion without using any correction factor; the excellent match between predictions and experimental results provided further proof of the quality of the parameters obtained with the proposed method.
Resumo:
Cell invasion involves a population of cells that migrate along a substrate and proliferate to a carrying capacity density. These two processes, combined, lead to invasion fronts that move into unoccupied tissues. Traditional modelling approaches based on reaction–diffusion equations cannot incorporate individual–level observations of cell velocity, as information propagates with infinite velocity according to these parabolic models. In contrast, velocity jump processes allow us to explicitly incorporate individual–level observations of cell velocity, thus providing an alternative framework for modelling cell invasion. Here, we introduce proliferation into a standard velocity–jump process and show that the standard model does not support invasion fronts. Instead, we find that crowding effects must be explicitly incorporated into a proliferative velocity–jump process before invasion fronts can be observed. Our observations are supported by numerical and analytical solutions of a novel coupled system of partial differential equations, including travelling wave solutions, and associated random walk simulations.
Resumo:
Motor unit number estimation (MUNE) is a method which aims to provide a quantitative indicator of progression of diseases that lead to loss of motor units, such as motor neurone disease. However the development of a reliable, repeatable and fast real-time MUNE method has proved elusive hitherto. Ridall et al. (2007) implement a reversible jump Markov chain Monte Carlo (RJMCMC) algorithm to produce a posterior distribution for the number of motor units using a Bayesian hierarchical model that takes into account biological information about motor unit activation. However we find that the approach can be unreliable for some datasets since it can suffer from poor cross-dimensional mixing. Here we focus on improved inference by marginalising over latent variables to create the likelihood. In particular we explore how this can improve the RJMCMC mixing and investigate alternative approaches that utilise the likelihood (e.g. DIC (Spiegelhalter et al., 2002)). For this model the marginalisation is over latent variables which, for a larger number of motor units, is an intractable summation over all combinations of a set of latent binary variables whose joint sample space increases exponentially with the number of motor units. We provide a tractable and accurate approximation for this quantity and also investigate simulation approaches incorporated into RJMCMC using results of Andrieu and Roberts (2009).
Resumo:
Purpose: To investigate the effects of an acute multinutrient supplement on game-based running performance, peak power output, anaerobic by-products, hormonal profiles, markers of muscle damage, and perceived muscular soreness before, immediately after, and 24 h following competitive rugby union games. Methods: Twelve male rugby union players ingested either a comprehensive multinutrient supplement (SUPP), [RE-ACTIVATE:01], or a placebo (PL) for 5 d. Participants then performed a competitive rugby union game (with global positioning system tracking), with associated blood draws and vertical jump assessments pre, immediately post and 24 h following competition. Results: SUPP ingestion resulted in moderate to large effects for augmented 1st half very high intensity running (VHIR) mean speed (5.9 ± 0.4 vs 4.8 ± 2.3 m·min–1; d= 0.93). Further, moderate increases in 2nd half VHIR distance (137 ± 119 vs 83 ± 89 m; d= 0.73) and VHIR mean speed (5.9 ± 0.6 v 5.3 ± 1.7 m·min–1; d= 0.56) in SUPP condition were also apparent. Postgame aspartate aminotransferase (AST; 44.1 ± 11.8 vs 37.0 ± 3.2 UL; d= 1.16) and creatine kinase (CK; 882 ± 472 vs. 645 ± 123 UL; d= 0.97) measures demonstrated increased values in the SUPP condition, while AST and CK values correlated with 2nd half VHIR distance (r= –0.71 and r= –0.76 respectively). Elevated C-reactive protein (CRP) was observed postgame in both conditions; however, it was significantly blunted with SUPP (P= .05). Conclusions: These findings suggest SUPP may assist in the maintenance of VHIR during rugby union games, possibly via the buffering qualities of SUPP ingredients. However, correlations between increased work completed at very high intensities and muscular degradation in SUPP conditions, may mask any anticatabolic properties of the supplement.
Resumo:
In this paper, we analyse the impact of a (small) heterogeneity of jump type on the most simple localized solutions of a 3-component FitzHugh–Nagumo-type system. We show that the heterogeneity can pin a 1-front solution, which travels with constant (non-zero) speed in the homogeneous setting, to a fixed, explicitly determined, distance from the heterogeneity. Moreover, we establish the stability of this heterogeneous pinned 1-front solution. In addition, we analyse the pinning of 1-pulse, or 2-front, solutions. The paper is concluded with simulations in which we consider the dynamics and interactions of N-front patterns in domains with M heterogeneities of jump type (N = 3, 4, M ≥ 1).