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.

Induction of anti-chlamydial mucosal immunity by transcutaneous immunization is enhanced by topical application of GM-CSF

Transcutaneous immunization (TCI) involves the direct application of antigen plus adjuvant to skin, taking advantage of the large numbers of Langerhans cells and other resident skin dendritic cells, that process antigen then migrate to draining lymph nodes where immune responses are initiated. We have used this form of immunization to protect mice against genital tract and respiratory tract chlamydial infection. Protection was associated with local antibody responses in the vagina, uterus and lung as well as strong Th1 responses in the lymph nodes draining the reproductive tract and lungs respectively. In this study we show that topical application of GM-CSF to skin enhances the numbers and activation status of epidermal dendritic cells. Topical application of GM-CSF also increased the immune responses elicited by TCI. GM-CSF supplementation greatly increased cytokine (IFNgamma and IL-4) gene expression in lymph node and splenic cells compared to cells from animals immunized without GM-CSF. IgG responses in serum, uterine lavage and bronchoalveolar lavage and IgA responses in vaginal lavage were also increased by topical application of GM-CSF. The studies show that TCI induces protection against genital and respiratory tract chlamydial infections and that topical application of cytokines such as GM-CSF can enhance TCI-induced antibody and cell-mediated immunity.

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.

Oscillatory regulation of Hes1: Discrete stochastic delay modelling and simulation

Discrete stochastic simulations are a powerful tool for understanding the dynamics of chemical kinetics when there are small-to-moderate numbers of certain molecular species. In this paper we introduce delays into the stochastic simulation algorithm, thus mimicking delays associated with transcription and translation. We then show that this process may well explain more faithfully than continuous deterministic models the observed sustained oscillations in expression levels of hes1 mRNA and Hes1 protein.

Comparison of intranasal and transcutaneous immunization for induction of protective immunity against chlamydia muridarum respiratory tract infection

Chlamydia pneumoniae causes a range of respiratory infections including bronchitis, pharyngitis and pneumonia. Infection has also been implicated in exacerbation/initiation of asthma and chronic obstructive pulmonary disease (COPD) and may play a role in atherosclerosis and Alzheimer's disease. We have used a mouse model of Chlamydia respiratory infection to determine the effectiveness of intranasal (IN) and transcutaneous immunization (TCI) to prevent Chlamydia lung infection. Female BALB/c mice were immunized with chlamydial major outer membrane protein (MOMP) mixed with cholera toxin and CpG oligodeoxynucleotide adjuvants by either the IN or TCI routes. Serum and bronchoalveolar lavage (BAL) were collected for antibody analysis. Mononuclear cells from lung-draining lymph nodes were stimulated in vitro with MOMP and cytokine mRNA production determined by real time PCR. Animals were challenged with live Chlamydia and weighed daily following challenge. At day 10 (the peak of infection) animals were sacrificed and the numbers of recoverable Chlamydia in lungs determined by real time PCR. MOMP-specific antibody-secreting cells in lung tissues were also determined at day 10 post-infection. Both IN and TCI protected animals against weight loss compared to non-immunized controls with both immunized groups gaining weight by day 10-post challenge while controls had lost 6% of body weight. Both immunization protocols induced MOMP-specific IgG in serum and BAL while only IN immunization induced MOMP-specific IgA in BAL. Both immunization routes resulted in high numbers of MOMP-specific antibody-secreting cells in lung tissues (IN > TCI). Following in vitro re-stimulation of lung-draining lymph node cells with MOMP; IFNγ mRNA increased 20-fold in cells from IN immunized animals (compared to non-immunized controls) while IFNγ levels increased 6- to 7-fold in TCI animals. Ten days post challenge non-immunized animals had >7000 IFU in their lungs, IN immunized animals <50 IFU and TCI immunized animals <1500 IFU. Thus, both intranasal and transcutaneous immunization protected mice against respiratory challenge with Chlamydia. The best protection was obtained following IN immunization and correlated with IFNγ production by mononuclear cells in lung-draining LN and MOMP-specific IgA in BAL.

Occurrence of LINE, gypsy-like, and copia-like retrotransposons in the clonally propagated sweet potato (Ipomoea batatas L.)

Retrotransposons are a class of transposable elements that represent a major fraction of the repetitive DNA of most eukaryotes. Their abundance stems from their expansive replication strategies. We screened and isolated sequence fragments of long terminal repeat (LTR), gypsy-like reverse transcriptase (rt) and gypsy-like envelope (env) domains, and two partial sequences of non-LTR retrotransposons, long interspersed element (LINE), in the clonally propagated allohexaploid sweet potato (Ipomoea batatas (L.) Lam.) genome. Using dot-blot hybridization, these elements were found to be present in the ~1597 Mb haploid sweet potato genome with copy numbers ranging from ~50 to ~4100 as observed in the partial LTR (IbLtr-1) and LINE (IbLi-1) sequences, respectively. The continuous clonal propagation of sweet potato may have contributed to such a multitude of copies of some of these genomic elements. Interestingly, the isolated gypsy-like env and gypsy-like rt sequence fragments, IbGy-1 (~2100 copies) and IbGy-2 (~540 copies), respectively, were found to be homologous to the Bagy-2 cDNA sequences of barley (Hordeum vulgare L.). Although the isolated partial sequences were found to be homologous to other transcriptionally active elements, future studies are required to determine whether they represent elements that are transcriptionally active under normal and (or) stressful conditions.

Equitable remedies for elder financial abuse in inter vivos transactions

Elder law is a growing area of legal practice due to the increasing numbers and proportions of older people in Australian society. The older generation has greater access to financial resources by way of retirement capital than ever before. Coupled with the current generation’s high level of debt and an increasing dependence on inheritances to meet these debts, this has created an environment in which the potential for elder financial abuse is increasing. This article examines how equitable remedies can be used as an avenue of redress for elder financial abuse. The effectiveness of these remedies, and in particular the prospect of a costs order being awarded against the perpetrator of the abuse in successful claims, may act as a deterrent and assist in preventing elder financial abuse from occurring.

Systemic inflammatory responses to maximal versus submaximal lengthening contractions of the elbow flexors

We compared changes in markers of muscle damage and systemic inflammation after submaximal and maximal lengthening muscle contractions of the elbow flexors. Using a cross-over design, 10 healthy young men not involved in resistance training completed a submaximal trial (10 sets of 60 lengthening contractions at 10% maximum isometric strength, 1 min rest between sets), followed by a maximal trial (10 sets of three lengthening contractions at 100% maximum isometric strength, 3 min rest between sets). Lengthening contractions were performed on an isokinetic dynamometer. Opposite arms were used for the submaximal and maximal trials, and the trials were separated by a minimum of two weeks. Blood was sampled before, immediately after, 1 h, 3 h, and 1-4 d after each trial. Total leukocyte and neutrophil numbers, and the serum concentration of soluble tumor necrosis factor-alpha receptor 1 were elevated after both trials (P < 0.01), but there were no differences between the trials. Serum IL-6 concentration was elevated 3 h after the submaximal contractions (P < 0.01). The concentrations of serum tumor necrosis factor-alpha, IL-1 receptor antagonist, IL-10, granulocyte-colony stimulating factor and plasma C-reactive protein remained unchanged following both trials. Maximum isometric strength and range of motion decreased significantly (P < 0.001) after both trials, and were lower from 1-4 days after the maximal contractions compared to the submaximal contractions. Plasma myoglobin concentration and creatine kinase activity, muscle soreness and upper arm circumference all increased after both trials (P < 0.01), but were not significantly different between the trials. Therefore, there were no differences in markers of systemic inflammation, despite evidence of greater muscle damage following maximal versus submaximal lengthening contractions of the elbow flexors.

Look before you leap : a confidence-based method for selecting species criticality while avoiding negative populations in τ-leaping

The stochastic simulation algorithm was introduced by Gillespie and in a different form by Kurtz. There have been many attempts at accelerating the algorithm without deviating from the behavior of the simulated system. The crux of the explicit τ-leaping procedure is the use of Poisson random variables to approximate the number of occurrences of each type of reaction event during a carefully selected time period, τ. This method is acceptable providing the leap condition, that no propensity function changes “significantly” during any time-step, is met. Using this method there is a possibility that species numbers can, artificially, become negative. Several recent papers have demonstrated methods that avoid this situation. One such method classifies, as critical, those reactions in danger of sending species populations negative. At most, one of these critical reactions is allowed to occur in the next time-step. We argue that the criticality of a reactant species and its dependent reaction channels should be related to the probability of the species number becoming negative. This way only reactions that, if fired, produce a high probability of driving a reactant population negative are labeled critical. The number of firings of more reaction channels can be approximated using Poisson random variables thus speeding up the simulation while maintaining the accuracy. In implementing this revised method of criticality selection we make use of the probability distribution from which the random variable describing the change in species number is drawn. We give several numerical examples to demonstrate the effectiveness of our new method.

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.

Development of a high-level gene expression system for the production of bioplastics in sugarcane

Despite various approaches, the production of biodegradable plastics such as polyhydroxybutyrate (PHB) in transgenic plants has met with limited success due largely to low expression levels. Even in the few instances where high levels of protein expression have been reported, the transgenic plants have been stunted indicating PHB is phytotoxic (Poirier 2002). This PhD describes the application of a novel virus-based gene expression technology, termed InPAct („In Plant Activation.), for the production of PHB in tobacco and sugarcane. InPAct is based on the rolling circle replication mechanism by which circular ssDNA viruses replicate and provides a system for controlled, high-level gene expression. Based on these features, InPAct was thought to represent an ideal system to enable the controlled, high-level expression of the three phb genes (phbA, phbB and phbC) required for PHB production in sugarcane at a preferred stage of plant growth. A Tobacco yellow dwarf virus (TbYDV)-based InPAct-phbA vector, as well as linear vectors constitutively expressing phbB and phbC were constructed and different combinations were used to transform tobacco leaf discs. A total of four, eight, three and three phenotypically normal tobacco lines were generated from discs transformed with InPAct-phbA, InPAct-phbA + p1300-TaBV P-phbB/phbC- 35S T, p1300-35S P-phbA-NOS T + p1300-TaBV P-phbB/phbC-35S T and InPAct-GUS, respectively. To determine whether the InPAct cassette could be activated in the presence of the TbYDV Rep, leaf samples from the eight InPActphbA + p1300-TaBV P-phbB/phbC-35S T plants were agroinfiltrated with p1300- TbYDV-Rep/RepA. Three days later, successful activation was indicated by the detection of episomes using both PCR and Southern analysis. Leaf discs from the eight InPAct-phbA + p1300-TaBV P-phbB/phbC-35S T transgenic plant lines were agroinfiltrated with p1300-TbYDV-Rep/RepA and leaf tissue was collected ten days post-infiltration and examined for the presence of PHB granules. Confocal microscopy and TEM revealed the presence of typical PHB granules in five of the eight lines, thus demonstrating the functionality of InPActbased PHB production in tobacco. However, analysis of leaf extracts by HPLC failed to detect the presence of PHB suggesting only very low level expression levels. Subsequent molecular analysis of three lines revealed low levels of correctly processed mRNA from the catalase intron contained within the InPAct cassette and also the presence of cryptic splice sites within the intron. In an attempt to increase expression levels, new InPAct-phb cassettes were generated in which the castorbean catalase intron was replaced with a synthetic intron (syntron). Further, in an attempt to both increase and better control Rep/RepA-mediated activation of InPAct cassettes, Rep/RepA expression was placed under the control of a stably integrated alc switch. Leaf discs from a transgenic tobacco line (Alc ML) containing 35S P-AlcR-AlcA P-Rep/RepA were supertransformed with InPAct-phbAsyn or InPAct-GUSsyn using Agrobacterium and three plants (lines) were regenerated for each construct. Analysis of the RNA processing of the InPAct-phbAsyn cassette revealed highly efficient and correct splicing of the syntron, thus supporting its inclusion within the InPAct system. To determine the efficiency of the alc switch to activate InPAct, leaf material from the three Alc ML + InPAct-phbAsyn lines was either agroinfiltrated with 35S P-Rep/RepA or treated with ethanol. Unexpectedly, episomes were detected not only in the infiltrated and ethanol treated samples, but also in non-treated samples. Subsequent analysis of transgenic Alc ML + InPAct-GUS lines, confirmed that the alc switch was leaky in tissue culture. Although this was shown to be reversible once plants were removed from the tissue culture environment, it made the regeneration of Alc ML + InPAct-phbsyn plant lines extremely difficult, due to unintentional Rep expression and therefore high levels of phb expression and phytotoxic PHB production. Two Alc ML + InPAct-phbAsyn + p1300-TaBV P-phbB/phbC-35S T transgenic lines were able to be regenerated, and these were acclimatised, alcohol-treated and analysed. Although episome formation was detected as late as 21 days post activation, no PHB was detected in the leaves of any plants using either microscopy or HPLC, suggesting the presence of a corrupt InPAct-phbA cassette in both lines. The final component of this thesis involved the application of both the alc switch and the InPAct systems to sugarcane in an attempt to produce PHB. Initial experiments using transgenic Alc ML + InPAct-GUS lines indicated that the alc system was not functional in sugarcane under the conditions tested. The functionality of the InPAct system, independent of the alc gene switch, was subsequently examined by bombarding the 35S Rep/RepA cassette into leaf and immature leaf whorl cells derived from InPAct-GUS transgenic sugarcane plants. No GUS expression was observed in leaf tissue, whereas weak and irregular GUS expression was observed in immature leaf whorl tissue derived from two InPAct- GUS lines and two InPAct-GUS + 35S P-AlcR-AlcA P-GUS lines. The most plausible reason to explain the inconsistent and low levels of GUS expression in leaf whorls is a combination of low numbers of sugarcane cells in the DNA replication-conducive S-phase and the irregular and random nature of sugarcane cells bombarded with Rep/RepA. This study details the first report to develop a TbYDV-based InPAct system under control of the alc switch to produce PHB in tobacco and sugarcane. Despite the inability to detect quantifiable levels of PHB levels in either tobacco or sugarcane, the findings of this study should nevertheless assist in the further development of both the InPAct system and the alc system, particularly for sugarcane and ultimately lead to an ethanol-inducible InPAct gene expression system for the production of bioplastics and other proteins of commercial value in plants.

Parallel implementation of stochastic simulation for large scale cellular processes

Experimental and theoretical studies have shown the importance of stochastic processes in genetic regulatory networks and cellular processes. Cellular networks and genetic circuits often involve small numbers of key proteins such as transcriptional factors and signaling proteins. In recent years stochastic models have been used successfully for studying noise in biological pathways, and stochastic modelling of biological systems has become a very important research field in computational biology. One of the challenge problems in this field is the reduction of the huge computing time in stochastic simulations. Based on the system of the mitogen-activated protein kinase cascade that is activated by epidermal growth factor, this work give a parallel implementation by using OpenMP and parallelism across the simulation. Special attention is paid to the independence of the generated random numbers in parallel computing, that is a key criterion for the success of stochastic simulations. Numerical results indicate that parallel computers can be used as an efficient tool for simulating the dynamics of large-scale genetic regulatory networks and cellular processes

Higher order spectral analysis of Chua's circuit

Higher order spectral analysis is used to investigate nonlinearities in time series of voltages measured from a realization of Chua's circuit. For period-doubled limit cycles, quadratic and cubic nonlinear interactions result in phase coupling and energy exchange between increasing numbers of triads and quartets of Fourier components as the nonlinearity of the system is increased. For circuit parameters that result in a chaotic Rossler-type attractor, bicoherence and tricoherence spectra indicate that both quadratic and cubic nonlinear interactions are important to the dynamics. When the circuit exhibits a double-scroll chaotic attractor the bispectrum is zero, but the tricoherences are high, consistent with the importance of higher-than-second order nonlinear interactions during chaos associated with the double scroll.

Regenerateengageexperiment in 2011! - News from Spatial Science Discipline at QUT

Abstract: Goals and potential impacts of QUT corporate Blueprint3 framework, university has made significant investments in physical infrastructure, and investments to improve staff profiles, particularly in relation to science, technology, engineering, and mathematics (STEM) disciplines. The most significant physical change to the Faculty’s infrastructure has seen new workshop and teaching and research spaces located in Science and Technology precinct under construction. Also includes Alumni news, input and output numbers Spatial Science discussion, Work Integrated Learning (WIL) in 2011, some key teaching administrative dates in 2011.

Spatial Science News from QUT April 2011

Multiple awards for Spatial/Surveying lecturer, raising entry quality and commencing numbers at QUT, Gardens Point rapt in promise of things to come, STEM building progress.