Generalized binomial Tau-leap method for biochemical kinetics incorporating both delay and intrinsic noise

The delay stochastic simulation algorithm (DSSA) by Barrio et al. [Plos Comput. Biol.2, 117–E (2006)] was developed to simulate delayed processes in cell biology in the presence of intrinsic noise, that is, when there are small-to-moderate numbers of certain key molecules present in a chemical reaction system. These delayed processes can faithfully represent complex interactions and mechanisms that imply a number of spatiotemporal processes often not explicitly modeled such as transcription and translation, basic in the modeling of cell signaling pathways. However, for systems with widely varying reaction rate constants or large numbers of molecules, the simulation time steps of both the stochastic simulation algorithm (SSA) and the DSSA can become very small causing considerable computational overheads. In order to overcome the limit of small step sizes, various τ-leap strategies have been suggested for improving computational performance of the SSA. In this paper, we present a binomial τ- DSSA method that extends the τ-leap idea to the delay setting and avoids drawing insufficient numbers of reactions, a common shortcoming of existing binomial τ-leap methods that becomes evident when dealing with complex chemical interactions. The resulting inaccuracies are most evident in the delayed case, even when considering reaction products as potential reactants within the same time step in which they are produced. Moreover, we extend the framework to account for multicellular systems with different degrees of intercellular communication. We apply these ideas to two important genetic regulatory models, namely, the hes1 gene, implicated as a molecular clock, and a Her1/Her 7 model for coupled oscillating cells.

Graphene-like nano-sheets for surface acoustic wave gas sensor applications

The gas sensing properties of graphene-like nano-sheets deposited on 36° YX lithium tantalate (LiTaO3) surface acoustic wave (SAW) transducers are reported. The thin graphene-like nano-sheets were produced via the reduction of graphite oxide which was deposited on SAW interdigitated transducers (IDTs). Their sensing performance was assessed towards hydrogen (H2) and carbon monoxide (CO) in a synthetic air carrier gas at room temperature (25 °C) and 40 °C. Raman and X-ray photoelectron spectroscopy (XPS) revealed that the deposited graphite oxide (GO) was not completely reduced creating small, graphitic nanocrystals ∼2.7 nm in size. © 2008 Elsevier B.V.

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.

Investigation of phenol degradation : True reaction kinetics on fixed film titanium dioxide photocatalyst

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.

Ethane-1,2-diaminium 4,4'-sulfonyldibenzoate

In the title compound 2(C2H10N2+) C14H8O6S2-, both the ethylenediaminium cations and the 4,4'-sulfonyldibenzoate dianions have crystallographic two-fold rotational symmetry and are interlinked by aminium N-H...O(carboxyl) hydrogen-bonding associations giving sheets which extend along (101) and are further linked down b forming a three-dimensional structure.

Poly[tetra-mu-aqua-diaquatetrakis[mu-(E)-2-nitrocinnamato]tetrarubidium]

In the structure of the title compound [Rb4(C9H6NO4)4(H~2~O)6]n, the asymmetric unit comprises four rubidium complex cations, two of which have an RbO7 coordination polyhedron with a monocapped distorted octahedral stereochemistry and two of which have a distorted RbO6 octahedral coordination. The bonding about both the seven-coordinate centres is similar, comprising one monodentate water molecule together with three bridging water molecules and three carboxylate O-atom donors, two of which are bridging. The environments about the six-coordinate cations are also similar, comprising a monodentate nitro O-atom donor, a bridging water molecule and four bridging carboxylate O-atom donors [overall Rb-O range, 2.849(2)-3.190(2)A]. The coordination leads to a two-dimensional polymeric structure extending parallel to (001), which is stabilized by interlayer water O-H...O hydrogen-bonding associations to water, carboxyl and nitro O-atom acceptors, together with weak inter-ring pi--pi interactions [minimum ring centroid separation = 3.5319(19)A].

catena-Poly[[diaqualithium]-mu-[rac-cis-(2-carboxycyclohexane-1-carboxylato-kO1:O2]]

In the structure of the title compound, [C8H11LiO4(H2O)2]n the distorted tetrahadral LiO4 coordination sphere comprises two water molecules and two carboxyl O-donors from separate bridging cis-2-carboxycyclohexane-1-carboxylate monoanions [Li-O range, 1.887(4)-1.946(3)A], giving chain substructures which extend along (010). Water-water and water-carboxyl O-H...O hydrogen bonds stabilize these chain structures and provide inter-chain links, resulting in a two-dimensional layered structure extending across (011).

Synthesis and thermal stability of hydrotalcites based upon gallium

Hydrotalcites based upon gallium as a replacement for aluminium in hydrotalcite over a Mg/Al ratio of 2:1 to 4:1 were synthesised. The d(003) spacing varied from 7.83 A ° for the 2:1 hydrotalcite to 8.15 A ° for the 3:1 gallium containing hydrotalcite. A comparison is made with the Mg Al hydrotalcite in which the d(003) spacing for the Mg/Al hydrotalcite varied from 7.62 A ° for the 2:1Mg hydrotalcite to 7.98 A ° for the 4:1 hydrotalcite. The thermal stability of the gallium containing hydrotalcite was determined using thermogravimetric analysis. Four mass loss steps at 77, 263–280,485 and 828 degrees C with mass losses of 10.23, 21.55, 5.20 and 7.58% are attributed to dehydration, dehydroxylation and decarbonation. The thermal stability of the galliumcontaining hydrotalcite is slightly less than the aluminium hydrotalcite.