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.

A software tool for simulating spatial separation and kinetics in biological membranes

Self-segregation and compartimentalisation are observed experimentally to occur spontaneously on live membranes as well as reconstructed model membranes. It is believed that many of these processes are caused or supported by anomalous diffusive behaviours of biomolecules on membranes due to the complex and heterogeneous nature of these environments. These phenomena are on the one hand of great interest in biology, since they may be an important way for biological systems to selectively localize receptors, regulate signaling or modulate kinetics; and on the other, they provide an inspiration for engineering designs that mimick natural systems. We present an interactive software package we are developing for the purpose of simulating such processes numerically using a fundamental Monte Carlo approach. This program includes the ability to simulate kinetics and mass transport in the presence of either mobile or immobile obstacles and other relevant structures such as liquid-ordered lipid microdomains. We also present preliminary simulation results regarding the selective spatial localization and chemical kinetics modulating power of immobile obstacles on the membrane, obtained using the program.

A fibrocontractive mechanochemical model of dermal wound closure incorporating realistic growth factor kinetics

Fibroblasts and their activated phenotype, myofibroblasts, are the primary cell types involved in the contraction associated with dermal wound healing. Recent experimental evidence indicates that the transformation from fibroblasts to myofibroblasts involves two distinct processes: the cells are stimulated to change phenotype by the combined actions of transforming growth factor β (TGFβ) and mechanical tension. This observation indicates a need for a detailed exploration of the effect of the strong interactions between the mechanical changes and growth factors in dermal wound healing. We review the experimental findings in detail and develop a model of dermal wound healing that incorporates these phenomena. Our model includes the interactions between TGFβ and collagenase, providing a more biologically realistic form for the growth factor kinetics than those included in previous mechanochemical descriptions. A comparison is made between the model predictions and experimental data on human dermal wound healing and all the essential features are well matched.

Clinical strategies for the alleviation of contractures from a predictive mathematical model of dermal repair

Hypertrophic scars arise when there is an overproduction of collagen during wound healing. These are often associated with poor regulation of the rate of programmed cell death(apoptosis) of the cells synthesizing the collagen or by an exuberant inflammatory response that prolongs collagen production and increases wound contraction. Severe contractures that occur, for example, after a deep burn can cause loss of function especially if the wound is over a joint such as the elbow or knee. Recently, we have developed a morphoelastic mathematical model for dermal repair that incorporates the chemical, cellular and mechanical aspects of dermal wound healing. Using this model, we examine pathological scarring in dermal repair by first assuming a smaller than usual apoptotic rate for myofibroblasts, and then considering a prolonged inflammatory response, in an attempt to determine a possible optimal intervention strategy to promote normal repair, or terminate the fibrotic scarring response. Our model predicts that in both cases it is best to apply the intervention strategy early in the wound healing response. Further, the earlier an intervention is made, the less aggressive the intervention required. Finally, if intervention is conducted at a late time during healing, a significant intervention is required; however, there is a threshold concentration of the drug or therapy applied, above which minimal further improvement to wound repair is obtained.

A continuum model of the growth of engineered epidermal skin substitutes

We present a porous medium model of the growth and deterioration of the viable sublayers of an epidermal skin substitute. It consists of five species: cells, intracellular and extracellular calcium, tight junctions, and a hypothesised signal chemical emanating from the stratum corneum. The model is solved numerically in Matlab using a finite difference scheme. Steady state calcium distributions are predicted that agree well with the experimental data. Our model also demonstrates epidermal skin substitute deterioration if the calcium diffusion coefficient is reduced compared to reported values in the literature.

A new geological model of Cenozoic New Zealand suggests a complex phylogeographic history for New Zealand terrestrial taxa

The new model of North Island Cenozoic palaeogeography developed by Kamp et al. has a range of important implications for the evolution of New Zealand terrestrial taxa over the past 30 Ma. Key aspects include the prolonged isolation of the biota on the North Island landmass from the larger and more diverse greater South Island, and the founding of North Island taxa from the potentially unusual ecosystem of a small island around Northland. The prolonged period of isolation is expected to have generated deep phylogenetic splits within taxa present on both islands, and an important current aim should be to identify such signals in surviving endemics to start building a picture of the historical phylogeography, and inferred ecology of both islands through the Cenozoic. Given the potential differences in founding terrestrial species and climatic conditions, it seems likely that the ecology may have been very diferent between the North and South Islands. New genetic data from the 10 or so species of extinct moa suggest that the radiation of moa was much more recent than previously suggested, and reveals a complex pattern that is inferred to result from the interplay of the Cenozoic biogeography, marine barriers, and glacial cycles.

Investigating the potential value of individual parameters of histological grading systems in a sheep model of cartilage damage : the modified Mankin method

A total histological grade does not necessarily distinguish between different manifestations of cartilage damage or degeneration. An accurate and reliable histological assessment method is required to separate normal and pathological tissue within a joint during treatment of degenerative joint conditions and to sub-classify the latter in meaningful ways. The Modified Mankin method may be adaptable for this purpose. We investigated how much detail may be lost by assigning one composite score/grade to represent different degenerative components of the osteoarthritic condition. We used four ovine injury models (sham surgery, anterior cruciate ligament/medial collateral ligament instability, simulated anatomic anterior cruciate ligament reconstruction and meniscal removal) to induce different degrees and potentially 'types' (mechanisms) of osteoarthritis. Articular cartilage was systematically harvested, prepared for histological examination and graded in a blinded fashion using a Modified Mankin grading method. Results showed that the possible permutations of cartilage damage were significant and far more varied than the current intended use that histological grading systems allow. Of 1352 cartilage specimens graded, 234 different manifestations of potential histological damage were observed across 23 potential individual grades of the Modified Mankin grading method. The results presented here show that current composite histological grading may contain additional information that could potentially discern different stages or mechanisms of cartilage damage and degeneration in a sheep model. This approach may be applicable to other grading systems.

Offsetting greenhouse gas emissions through biological carbon sequestration in North Eastern Australia

The Kyoto Protocol recognises trees as a sink of carbon and a valid means to offset greenhouse gas emissions and meet internationally agreed emissions targets. This study details biological carbon sequestration rates for common plantation species Araucaria cunninghamii (hoop pine), Eucalyptus cloeziana, Eucalyptus argophloia, Pinus elliottii and Pinus caribaea var hondurensis and individual land areas required in north-eastern Australia to offset greenhouse gas emissions of 1000tCO 2e. The 3PG simulation model was used to predict above and below-ground estimates of biomass carbon for a range of soil productivity conditions for six representative locations in agricultural regions of north-eastern Australia. The total area required to offset 1000tCO 2e ranges from 1ha of E. cloeziana under high productivity conditions in coastal North Queensland to 45ha of hoop pine in low productivity conditions of inland Central Queensland. These areas must remain planted for a minimum of 30years to meet the offset of 1000tCO 2e.

An improved model of tumour-immune system interactions

The immune system plays an important role in defending the body against tumours and other threats. Currently, mechanisms involved in immune system interactions with tumour cells are not fully understood. Here we develop a mathematical tool that can be used in aiding to address this shortfall in understanding. This paper de- scribes a hybrid cellular automata model of the interaction between a growing tumour and cells of the innate and specific immune system including the effects of chemokines that builds on previous models of tumour-immune system interactions. In particular, the model is focused on the response of immune cells to tumour cells and how the dynamics of the tumour cells change due to the immune system of the host. We present results and predictions of in silico experiments including simulations of Kaplan-Meier survival-like curves.

Review of 3D model development with application to groundwater studies in New Zealand 1996 – 2011

Three dimensional geological modelling techniques have been applied since 1996 with an aim to characterise the lithological and chronological units of New Zealand’s many diverse aquifers. Models of property-scattered data have also been applied to assess physical properties of aquifers and the distribution of groundwater chemistry, including groundwater age, to inform an understanding of groundwater systems. These models, fundamental to understanding groundwater recharge, flow and discharge have found many uses as outlined in this paper.

Application of a human bone engineering platform to an in vitro and in vivo breast cancer metastasis model

Breast cancer in its advanced stage has a high predilection to the skeleton. Currently, treatment options of breast cancer-related bone metastasis are restricted to only palliative therapeutic modalities. This is due to the fact that mechanisms regarding the breast cancer celI-bone colonisation as well as the interactions of breast cancer cells with the bone microenvironment are not fully understood, yet. This might be explained through a lack of appropriate in vitro and in vivo models that are currently addressing the above mentioned issue. Hence the hypothesis that the translation of a bone tissue engineering platform could lead to improved and more physiological in vitro and in vivo model systems in order to investigate breast cancer related bone colonisation was embraced in this PhD thesis. Therefore the first objective was to develop an in vitro model system that mimics human mineralised bone matrix to the highest possible extent to examine the specific biological question, how the human bone matrix influences breast cancer cell behaviour. Thus, primary human osteoblasts were isolated from human bone and cultured under osteogenic conditions. Upon ammonium hydroxide treatment, a cell-free intact mineralised human bone matrix was left behind. Analyses revealed a similar protein and mineral composition of the decellularised osteoblast matrix to human bone. Seeding of a panel of breast cancer cells onto the bone mimicking matrix as well as reference substrates like standard tissue culture plastic and collagen coated tissue culture plastic revealed substrate specific differences of cellular behaviour. Analyses of attachment, alignment, migration, proliferation, invasion, as well as downstream signalling pathways showed that these cellular properties were influenced through the osteoblast matrix. The second objective of this PhD project was the development of a human ectopic bone model in NOD/SCID mice using medical grade polycaprolactone tricalcium phosphate (mPCL-TCP) scaffold. Human osteoblasts and mesenchymal stem cells were seeded onto an mPCL-TCP scaffold, fabricated using a fused deposition modelling technique. After subcutaneous implantation in conjunction with the bone morphogenetic protein 7, limited bone formation was observed due to the mechanical properties of the applied scaffold and restricted integration into the soft tissue of flank of NOD/SCID mice. Thus, a different scaffold fabrication technique was chosen using the same polymer. Electrospun tubular scaffolds were seeded with human osteoblasts, as they showed previously the highest amount of bone formation and implanted into the flanks of NOD/SCID mice. Ectopic bone formation with sufficient vascularisation could be observed. After implantation of breast cancer cells using a polyethylene glycol hydrogel in close proximity to the newly formed bone, macroscopic communication between the newly formed bone and the tumour could be observed. Taken together, this PhD project showed that bone tissue engineering platforms could be used to develop an in vitro and in vivo model system to study cancer cell colonisation in the bone microenvironment.

The expression of plasminogen activator system in a rat model of periodontal wound healing

BACKGROUND: The plasminogen activator system has been proposed to play a role in proteolytic degradation of extracellular matrices in tissue remodeling, including wound healing. The aim of this study was to elucidate the presence of components of the plasminogen activator system during different stages of periodontal wound healing. METHODS: Periodontal wounds were created around the molars of adult rats and healing was followed for 28 days. Immunohistochemical analyses of the healing tissues and an analysis of the periodontal wound healing fluid by ELISA were carried out for the detection of tissue-type plasminogen activator (t-PA), urokinase-type plasminogen activator (u-PA), and 2 plasminogen activator inhibitors (PAI-1 and PAI-2). RESULTS: During the early stages (days 1 to 3) of periodontal wound healing, PAI-1 and PAI-2 were found to be closely associated with the deposition of a fibrin clot in the gingival sulcus. These components were strongly associated with the infiltrating inflammatory cells around the fibrin clot. During days 3 to 7, u-PA, PAI-1, and PAI-2 were associated with cells (particularly monocytes/macrophages, fibroblasts, and endothelial cells) in the newly formed granulation tissue. During days 7 to 14, a new attachment apparatus was formed during which PAI-1, PAI-2, and u-PA were localized in both periodontal ligament fibroblasts (PDL) and epithelial cells at sites where these cells were attaching to the root surface. In the periodontal wound healing fluid, the concentration for t-PA increased and peaked during the first week. PAI-2 had a similar expression to t-PA, but at a lower level over the entire wound-healing period. CONCLUSIONS: These findings indicate that the plasminogen activator system is involved in the entire process of periodontal wound healing, in particular with the formation of fibrin matrix on the root surface and its replacement by granulation tissue, as well as the subsequent formation of the attachment of soft tissue to the root surface during the later stages of wound repair.

Physical and linkage maps for Drosophila serrata, a model species for studies of clinal adaptation and sexual selection

Drosophila serrata is a member of the montium group, which contains more than 98 species and until recently was considered a subgroup within the melanogaster group. This Drosophila species is an emerging model system for evolutionary quantitative genetics and has been used in studies of species borders, clinal variation and sexual selection. Despite the importance of D. serrata as a model for evolutionary research, our poor understanding of its genome remains a significant limitation. Here, we provide a first-generation gene-based linkage map and a physical map for this species. Consistent with previous studies of other drosophilids we observed strong conservation of genes within chromosome arms homologous with D. melanogaster but major differences in within-arm synteny. These resources will be a useful complement to ongoing genome sequencing efforts and QTL mapping studies in this species

An expressed sequence tag (EST) library for Drosophila serrata, a model system for sexual selection and climatic adaptation studies

The native Australian fly Drosophila serrata belongs to the highly speciose montium subgroup of the melanogaster species group. It has recently emerged as an excellent model system with which to address a number of important questions, including the evolution of traits under sexual selection and traits involved in climatic adaptation along latitudinal gradients. Understanding the molecular genetic basis of such traits has been limited by a lack of genomic resources for this species. Here, we present the first expressed sequence tag (EST) collection for D. serrata that will enable the identification of genes underlying sexually-selected phenotypes and physiological responses to environmental change and may help resolve controversial phylogenetic relationships within the montium subgroup.