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.

The structure of the mineral leogangite Cu10(OH)6(SO4)(AsO4)4•8H2O–implications for arsenic accumulation and removal

The objective of this research is to determine the molecular structure of the mineral leogangite. The formation of the types of arsenosulphate minerals offers a mechanism for arsenate removal from soils and mine dumps. Raman and infrared spectroscopy have been used to characterise the mineral. Observed bands are assigned to the stretching and bending vibrations of (SO4)2- and (AsO4)3- units, stretching and bending vibrations of hydrogen bonded (OH)- ions and Cu2+-(O,OH) units. The approximate range of O-H...O hydrogen bond lengths is inferred from the Raman spectra. Raman spectra of leogangite from different origins differ in that some spectra are more complex, where bands are sharp and the degenerate bands of (SO4)2- and (AsO4)3- are split and more intense. Lower wavenumbers of  H2O bending vibration in the spectrum may indicate the presence of weaker hydrogen bonds compared with those in a different leogangite samples. The formation of leogangite offers a mechanism for the removal of arsenic from the environment.

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.

The molecular structure of the mineral sarmientite Fe2(AsO4,SO4)2(OH)6•5H2O : implications for arsenic accumulation and removal

Sarmientite is an environmental mineral; its formation in soils enables the entrapment and immobilisation of arsenic. The mineral sarmientite is often amorphous making the application of X-ray diffraction difficult. Vibrational spectroscopy has been applied to the study of sarmientite. Bands are attributed to the vibrational units of arsenate, sulphate, hydroxyl and water. Raman bands at 794, 814 and 831 cm−1 are assigned to the ν3 (AsO4)3− antisymmetric stretching modes and the ν1 symmetric stretching mode is observed at 891 cm−1. Raman bands at 1003 and 1106 cm−1 are attributed to vibrations. The Raman band at 484 cm−1 is assigned to the triply degenerate (AsO4)3− bending vibration. The high intensity Raman band observed at 355 cm−1 (both lower and upper) is considered to be due to the (AsO4)3−ν2 bending vibration. Bands attributed to water and OH stretching vibrations are observed.

Vibrational spectroscopy of the multi anion mineral zykaite Fe4(AsO4)(SO4)(OH)•15H2O – implications for arsenate removal

Some minerals are colloidal and are poorly diffracting . Vibrational spectroscopy offers one of the few methods for the assessment of the structure of these types of minerals. Among this group of minerals is zykaite with formula Fe4(AsO4)(SO4)(OH)•15H2O. The objective of this research is to determine the molecular structure of the mineral zykaite using vibrational spectroscopy. Raman and infrared bands are attributed to the AsO43-, SO42- and water stretching vibrations. The sharp band at 3515 cm-1 is assigned to the stretching vibration of the OH units. This mineral offers a mechanism for the formation of more crystalline minerals such as scorodite and bukovskyite. Arsenate ions can be removed from aqueous systems through the addition of ferric compounds such as ferric chloride. This results in the formation of minerals such as zykaite and pitticite (Fe3+,AsO4,SO4,H2O).

The molecular structure of the multianion mineral hidalgoite PbAl3(AsO4)(SO4)(OH)6 : implications for arsenic removal from soils

The objective of this research is to determine the molecular structure of the mineral hidalgoite PbAl3(AsO4)(SO4)(OH)6 using vibrational spectroscopy. The mineral is found in old mine sites. Observed bands are assigned to the stretching and bending vibrations of (SO4)2- and (AsO4)3- units, stretching and bending vibrations of hydrogen bonded (OH)- ions and Al3+-(O,OH) units. The approximate range of O-H...O hydrogen bond lengths is inferred from the Raman and infrared spectra. Values of 2.6989 Å, 2.7682 Å, 2.8659 Å were obtained. The formation of hidalgoite may offer a mechanism for the removal of arsenic from the environment.