Nonlinear breathing modes at a defect

Recent molecular dynamics (MD) simulations of Cubero et al (1999) of a DNA duplex containing the 'rogue' base difluorotoluene (F) in place of a thymine (T) base show that breathing events can occur on the nanosecond timescale, whereas breathing events in a normal DNA duplex take place on the microsecond timescale. The main aim of this paper is to analyse a nonlinear Klein-Gordon lattice model of the DNA duplex including both nonlinear interactions between opposing bases and a defect in the interaction at one lattice site; each of which can cause localisation of energy. Solutions for a breather mode either side of the defect are derived using multiple-scales asymptotics and are pieced together across the defect to form a solution which includes the effects of the nonlinearity and the defect. We consider defects in the inter-chain interactions and in the along chain interactions. In most cases we find in-phase breather modes and/or out-of-phase breather modes, with one case displaying a shifted mode.

Towards a predictive model of Ca²⁺ puffs

We investigate key characteristics of Ca²⁺ puffs in deterministic and stochastic frameworks that all incorporate the cellular morphology of IP[subscript]3 receptor channel clusters. In a first step, we numerically study Ca²⁺ liberation in a three dimensional representation of a cluster environment with reaction-diffusion dynamics in both the cytosol and the lumen. These simulations reveal that Ca²⁺ concentrations at a releasing cluster range from 80 µM to 170 µM and equilibrate almost instantaneously on the time scale of the release duration. These highly elevated Ca²⁺ concentrations eliminate Ca²⁺ oscillations in a deterministic model of an IP[subscript]3R channel cluster at physiological parameter values as revealed by a linear stability analysis. The reason lies in the saturation of all feedback processes in the IP[subscript]3R gating dynamics, so that only fluctuations can restore experimentally observed Ca²⁺ oscillations. In this spirit, we derive master equations that allow us to analytically quantify the onset of Ca²⁺ puffs and hence the stochastic time scale of intracellular Ca²⁺ dynamics. Moving up the spatial scale, we suggest to formulate cellular dynamics in terms of waiting time distribution functions. This approach prevents the state space explosion that is typical for the description of cellular dynamics based on channel states and still contains information on molecular fluctuations. We illustrate this method by studying global Ca²⁺ oscillations.