Stochastic simulation of chemical reactions in spatially complex media

Discrete stochastic simulations, via techniques such as the Stochastic Simulation Algorithm (SSA) are a powerful tool for understanding the dynamics of chemical kinetics when there are low numbers of certain molecular species. However, an important constraint is the assumption of well-mixedness and homogeneity. In this paper, we show how to use Monte Carlo simulations to estimate an anomalous diffusion parameter that encapsulates the crowdedness of the spatial environment. We then use this parameter to replace the rate constants of bimolecular reactions by a time-dependent power law to produce an SSA valid in cases where anomalous diffusion occurs or the system is not well-mixed (ASSA). Simulations then show that ASSA can successfully predict the temporal dynamics of chemical kinetics in a spatially constrained environment.

Synthesis and characterisation of metal oxyhydroxide and oxide nanomaterials

In this work, a range of nanomaterials have been synthesised based on metal oxyhydroxides MO(OH), where M=Al, Co, Cr, etc. Through a self-assembly hydrothermal route, metal oxyhydroxide nanomaterials with various morphologies were successfully synthesised: one dimensional boehmite (AlO(OH)) nanofibres, zero dimensional indium hydroxide (In(OH)3) nanocubes and chromium oxyhydroxide (CrO(OH)) nanoparticles, as well as two dimensional cobalt hydroxide and oxyhydroxide (Co(OH)2 & CoO(OH)) nanodiscs. In order to control the synthetic nanomaterial morphology and growth, several factors were investigated including cation concentration, temperature, hydrothermal treatment time, and pH. Metal ion doping is a promising technique to modify and control the properties of materials by intentionally introducing impurities or defects into the material. Chromium was successfully applied as a dopant for fabricating doped boehmite nanofibres. The thermal stability of the boehmite nanofibres was enhanced by chromium doping, and the photoluminescence property was introduced to the chromium doped alumina nanofibres. Doping proved to be an efficient method to modify and functionalize nanomaterials. The synthesised nanomaterials were fully characterised by X-ray diffraction (XRD), transmission electron microscopy (TEM) combined with selected area electron diffraction (SAED), scanning electron microscopy (SEM), BET specific surface area analysis, X-ray photoelectron spectroscopy (XPS) and thermo gravimetric analysis (TGA). Hot-stage Raman and infrared emission spectroscopy were applied to study the chemical reactions during dehydration and dehydroxylation. The advantage of these techniques is that the changes in molecular structure can be followed in situ and at the elevated temperatures.

Diels-Alder reactions as an efficient route to high purity cyclic polymers

A simple and efficient route for the synthesis of cyclic polymer systems is presented. Linear furan protected α-maleimide-ω-cyclopentadienyl functionalized precursors (poly(methyl methacrylate) and poly(tert-butyl acrylate)) were synthesized via atom transfer radical polymerization (ATRP) and subsequent substitution of the bromine end-group with cyclopentadiene. Upon heating at high dilution, deprotection of the dieneophile occurs followed by an intramolecular Diels–Alder reaction yielding a high purity cyclic product.

Modeling intrinsic noise and delays in chemical kinetics of coupled autoregulated oscillating cells

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.

Stochastic chemical kinetics and the quasi-steady-state assumption : application to the stochastic simulation algorithm and chemical master equation

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.

A Krylov-based finite state projection algorithm for solving the chemical master equation arising in the discrete modelling of biological systems

Biochemical reactions underlying genetic regulation are often modelled as a continuous-time, discrete-state, Markov process, and the evolution of the associated probability density is described by the so-called chemical master equation (CME). However the CME is typically difficult to solve, since the state-space involved can be very large or even countably infinite. Recently a finite state projection method (FSP) that truncates the state-space was suggested and shown to be effective in an example of a model of the Pap-pili epigenetic switch. However in this example, both the model and the final time at which the solution was computed, were relatively small. Presented here is a Krylov FSP algorithm based on a combination of state-space truncation and inexact matrix-vector product routines. This allows larger-scale models to be studied and solutions for larger final times to be computed in a realistic execution time. Additionally the new method computes the solution at intermediate times at virtually no extra cost, since it is derived from Krylov-type methods for computing matrix exponentials. For the purpose of comparison the new algorithm is applied to the model of the Pap-pili epigenetic switch, where the original FSP was first demonstrated. Also the method is applied to a more sophisticated model of regulated transcription. Numerical results indicate that the new approach is significantly faster and extendable to larger biological models.

Accelerated leap methods for simulating discrete stochastic chemical kinetics

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.