Some novel techniques of parameter estimation for the dynamical models in biological systems

Inverse problems based on using experimental data to estimate unknown parameters of a system often arise in biological and chaotic systems. In this paper, we consider parameter estimation in systems biology involving linear and non-linear complex dynamical models, including the Michaelis–Menten enzyme kinetic system, a dynamical model of competence induction in Bacillus subtilis bacteria and a model of feedback bypass in B. subtilis bacteria. We propose some novel techniques for inverse problems. Firstly, we establish an approximation of a non-linear differential algebraic equation that corresponds to the given biological systems. Secondly, we use the Picard contraction mapping, collage methods and numerical integration techniques to convert the parameter estimation into a minimization problem of the parameters. We propose two optimization techniques: a grid approximation method and a modified hybrid Nelder–Mead simplex search and particle swarm optimization (MH-NMSS-PSO) for non-linear parameter estimation. The two techniques are used for parameter estimation in a model of competence induction in B. subtilis bacteria with noisy data. The MH-NMSS-PSO scheme is applied to a dynamical model of competence induction in B. subtilis bacteria based on experimental data and the model for feedback bypass. Numerical results demonstrate the effectiveness of our approach.

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.

Biological and environmental determinants

- describe the complex web of determinants as part of broad causal pathways that affect health - identify and discuss the range of physical, biological and environmental determinants that impact on health - suggest why it is important to the practice of public health that you understand how determinants contribute to health - understand the complexity of health and illness and the multifaceted role of health determinants - relate determinants of health to public health activity and realise the need for multisectoral action and multiple approaches when working to improve health

Biological responses of human bone marrow mesenchymal stem cells to Sr-M-Si (M = Zn, Mg) silicate bioceramics

Strontium (Sr), Zinc (Zn), magnesium (Mg), and silicon (Si) are reported to be essential trace elements for the growth and mineralization of bone. We speculated that the combination of these bioactive elements in bioceramics may be effective to regulate the osteogenic property of boneforming cells. In this study, two Sr-containing silicate bioceramics, Sr2ZnSi2O7 (SZS) and Sr2MgSi2O7 (SMS), were prepared. The biological response of human bone marrow mesenchymal stem cells (BMSCs) to the two bioceramics (in the forms of powders and dense ceramic bulks) was systematically studied. In powder form, the effect of powder extracts on the viability and alkaline phosphatase (ALP) activity of BMSCs was investigated. In ceramic disc form, both direct and indirect coculture of BMSCs with ceramic discs were used to investigate their biological response, including attachment, proliferation, ALP activity, and bone-related genes expression. Beta-tricalcium phosphate (b-TCP) and akermanite (Ca2MgSi2O7, CMS) were used as control materials. The results showed that the Sr, Zn, and Si (or Sr, Mg, and Si)-containing ionic products from SZS and SMS powders enhanced ALP activity of BMSCs, compared to those from b-TCP. Both SZS and SMS ceramic discs supported the growth of BMSCs, and most importantly, significantly enhanced the ALP activity and bone-related genes expression of BMSCs as compared to b-TCP. The results suggest that the specific combination of bioactive ions (Sr, Zn, Si, e.g.) in bioceramics is a viable way to improve the biological performance of biomaterials, and the form of materials and surface properties were nonnegligible factors to influence cell response.

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.