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.

Sequential Monte Carlo for Bayesian sequentially designed experiments for discrete data

In this paper we present a sequential Monte Carlo algorithm for Bayesian sequential experimental design applied to generalised non-linear models for discrete data. The approach is computationally convenient in that the information of newly observed data can be incorporated through a simple re-weighting step. We also consider a flexible parametric model for the stimulus-response relationship together with a newly developed hybrid design utility that can produce more robust estimates of the target stimulus in the presence of substantial model and parameter uncertainty. The algorithm is applied to hypothetical clinical trial or bioassay scenarios. In the discussion, potential generalisations of the algorithm are suggested to possibly extend its applicability to a wide variety of scenarios

Practical stability of approximating discrete-time filters with respect to model mismatch

This paper establishes practical stability results for an important range of approximate discrete-time filtering problems involving mismatch between the true system and the approximating filter model. Practical stability is established in the sense of an asymptotic bound on the amount of bias introduced by the model approximation. Our analysis applies to a wide range of estimation problems and justifies the common practice of approximating intractable infinite dimensional nonlinear filters by simpler computationally tractable filters.

Frame-rate stereopsis using non-parametric transforms and programmable logic

A frame-rate stereo vision system, based on non-parametric matching metrics, is described. Traditional metrics, such as normalized cross-correlation, are expensive in terms of logic. Non-parametric measures require only simple, parallelizable, functions such as comparators, counters and exclusive-or, and are thus very well suited to implementation in reprogrammable logic.

Synthesis of discrete allyl-functional polyethylene glycol for efficient hydrophilization of dihydr (polydimethylsiloxane)

Amphiphilic poly(ethylene glycol)-block-pol (dimethylsiloxane)-block-poly(ethylene glycol)(PEG-block-PDMS block-PEG) triblock copolymers have been successfully prepared via hydrosilylation using discrete and polydisperse PEG of various chain lengths. Facile synthesis of discrete PEG (dPEG) is achieved via systematic tosylation and etherification of lower glycols. Amphiphilicity of the dPEG block-PDMS-block-dPEG triblock copolymer is illustrated by dynamic light scattering (DLS) and measurement of the critical micelle concentration (CMC).

Discrete loading of the foot during hurdle drills: A pilot investigation

Introduction: Lower limb function in hurdling is patently asymmetrical. The lead limb undertakes the preparatory and landing steps while the trail limb contends with the hurdle and recovery steps. Discrete loading profiles of these steps will reflect the asymmetrical function and may provide useful insight into injury mechanisms. A pilot study was undertaken to determine the loading profiles of the hurdle, landing and recovery steps of elite male hurdlers. Equivalent data for steps between hurdles, where the running action is more symmetrical, were used for the purpose of comparison, simultaneously minimising the confounding effect of speed. Methodology: In-shoe pressures were recorded (FScan, 200 Hz) for four elite male hurdlers while they completed a routine hurdle drill at a self-selected fast but sub-race speed. The drill comprised of three consecutive hurdles. Data for the hurdle, landing and recovery steps of the first and second hurdles, along with data for the running steps between hurdles 1 and 2, and 2 and 3, were used for the purpose of analysis. Peak pressures within 1cm2 masks were determined for the hallux, first, central and fifth metatarsals (T1, M1, M2–4 and M5 respectively). Peak pressure (kPa) and loading duration (ms) for the hurdle, landing and recovery steps are reported as a percentage of the respective limb-matched values for between-hurdle steps. Results/discussion: For between-hurdle steps, T1, M1 and M2–4 peak pressures were 312/357, 356/306 and 362/368 kPa, lead/trail limbs respectively. For the hurdle, landing and recovery steps, pressures at T1 and M1 increased. For T1 the increases were in the order of 17%, 36% and 8% (hurdle, landing and recovery steps, respectively) while the corresponding increases at M1 were 7%, 54% and 20%. Pressures at M2–4 were similar for all steps, while M5 loaded erratically. For the between-hurdle steps, the loading durations at T1, M1 and M2–4, were 160/162, 170/142 and 190/191 ms, respectively. For the landing step, loading duration decreased for T1, M1and M2–4 (−8%, −19% and −18%, respectively). In the hurdle step, loading duration decreased for the metatarsals but not for T1. Conclusions: The hurdling action leads to regional pressure increases that act for shorter durations in comparison to the between-hurdle running steps. These changes are most notable at the first metatarsal, a common site of foot injury.

Lattice-free models of cell invasion : discrete simulations and travelling waves

Invasion waves of cells play an important role in development, disease and repair. Standard discrete models of such processes typically involve simulating cell motility, cell proliferation and cell-to-cell crowding effects in a lattice-based framework. The continuum-limit description is often given by a reaction–diffusion equation that is related to the Fisher–Kolmogorov equation. One of the limitations of a standard lattice-based approach is that real cells move and proliferate in continuous space and are not restricted to a predefined lattice structure. We present a lattice-free model of cell motility and proliferation, with cell-to-cell crowding effects, and we use the model to replicate invasion wave-type behaviour. The continuum-limit description of the discrete model is a reaction–diffusion equation with a proliferation term that is different from lattice-based models. Comparing lattice based and lattice-free simulations indicates that both models lead to invasion fronts that are similar at the leading edge, where the cell density is low. Conversely, the two models make different predictions in the high density region of the domain, well behind the leading edge. We analyse the continuum-limit description of the lattice based and lattice-free models to show that both give rise to invasion wave type solutions that move with the same speed but have very different shapes. We explore the significance of these differences by calibrating the parameters in the standard Fisher–Kolmogorov equation using data from the lattice-free model. We conclude that estimating parameters using this kind of standard procedure can produce misleading results.

Order tracking for discrete-random separation in variable speed conditions

The transmission path from the excitation to the measured vibration on the surface of a mechanical system introduces a distortion both in amplitude and in phase. Moreover, in variable speed conditions, the amplification/attenuation and the phase shift, due to the transfer function of the mechanical system, varies in time. This phenomenon reduces the effectiveness of the traditionally tachometer based order tracking, compromising the results of a discrete-random separation performed by a synchronous averaging. In this paper, for the first time, the extent of the distortion is identified both in the time domain and in the order spectrum of the signal, highlighting the consequences for the diagnostics of rotating machinery. A particular focus is given to gears, providing some indications on how to take advantage of the quantification of the disturbance to better tune the techniques developed for the compensation of the distortion. The full theoretical analysis is presented and the results are applied to an experimental case.

A Human identification technique using images of the iris and wavelet transform

A new approach for recognizing the iris of the human eye is presented. Zero-crossings of the wavelet transform at various resolution levels are calculated over concentric circles on the iris, and the resulting one-dimensional (1-D) signals are compared with model features using different dissimilarity functions.