Global parameter identification of stochastic reaction networks from single trajectories

We consider the problem of inferring the unknown parameters of a stochastic biochemical network model from a single measured time-course of the concentration of some of the involved species. Such measurements are available, e.g., from live-cell fluorescence microscopy in image-based systems biology. In addition, fluctuation time-courses from, e.g., fluorescence correlation spectroscopy provide additional information about the system dynamics that can be used to more robustly infer parameters than when considering only mean concentrations. Estimating model parameters from a single experimental trajectory enables single-cell measurements and quantification of cell--cell variability. We propose a novel combination of an adaptive Monte Carlo sampler, called Gaussian Adaptation, and efficient exact stochastic simulation algorithms that allows parameter identification from single stochastic trajectories. We benchmark the proposed method on a linear and a non-linear reaction network at steady state and during transient phases. In addition, we demonstrate that the present method also provides an ellipsoidal volume estimate of the viable part of parameter space and is able to estimate the physical volume of the compartment in which the observed reactions take place.Comment: Article in print as a book chapter in Springer's "Advances in Systems Biology

Global parameter estimation methods for stochastic biochemical systems

<p>Abstract</p> <p>Background</p> <p>The importance of stochasticity in cellular processes having low number of molecules has resulted in the development of stochastic models such as chemical master equation. As in other modelling frameworks, the accompanying rate constants are important for the end-applications like analyzing system properties (e.g. robustness) or predicting the effects of genetic perturbations. Prior knowledge of kinetic constants is usually limited and the model identification routine typically includes parameter estimation from experimental data. Although the subject of parameter estimation is well-established for deterministic models, it is not yet routine for the chemical master equation. In addition, recent advances in measurement technology have made the quantification of genetic substrates possible to single molecular levels. Thus, the purpose of this work is to develop practical and effective methods for estimating kinetic model parameters in the chemical master equation and other stochastic models from single cell and cell population experimental data.</p> <p>Results</p> <p>Three parameter estimation methods are proposed based on the maximum likelihood and density function distance, including probability and cumulative density functions. Since stochastic models such as chemical master equations are typically solved using a Monte Carlo approach in which only a finite number of Monte Carlo realizations are computationally practical, specific considerations are given to account for the effect of finite sampling in the histogram binning of the state density functions. Applications to three practical case studies showed that while maximum likelihood method can effectively handle low replicate measurements, the density function distance methods, particularly the cumulative density function distance estimation, are more robust in estimating the parameters with consistently higher accuracy, even for systems showing multimodality.</p> <p>Conclusions</p> <p>The parameter estimation methodologies described in this work have provided an effective and practical approach in the estimation of kinetic parameters of stochastic systems from either sparse or dense cell population data. Nevertheless, similar to kinetic parameter estimation in other modelling frameworks, not all parameters can be estimated accurately, which is a common problem arising from the lack of complete parameter identifiability from the available data.</p

Melatonin the "light of night" in human biology and adolescent idiopathic scoliosis

Melatonin "the light of night" is secreted from the pineal gland principally at night. The hormone is involved in sleep regulation, as well as in a number of other cyclical bodily activities and circadian rhythm in humans. Melatonin is exclusively involved in signalling the 'time of day' and 'time of year' (hence considered to help both clock and calendar functions) to all tissues and is thus considered to be the body's chronological pacemaker or 'Zeitgeber'. The last decades melatonin has been used as a therapeutic chemical in a large spectrum of diseases, mainly in sleep disturbances and tumours and may play a role in the biologic regulation of mood, affective disorders, cardiovascular system, reproduction and aging. There are few papers regarding melatonin and its role in adolescent idiopathic scoliosis (AIS). Melatonin may play a role in the pathogenesis of scoliosis (neuroendocrine hypothesis) but at present, the data available cannot clearly support this hypothesis. Uncertainties and doubts still surround the role of melatonin in human physiology and pathophysiology and future research is needed

Optimization in computational systems biology

Optimization aims to make a system or design as effective or functional as possible. Mathematical optimization methods are widely used in engineering, economics and science. This commentary is focused on applications of mathematical optimization in computational systems biology. Examples are given where optimization methods are used for topics ranging from model building and optimal experimental design to metabolic engineering and synthetic biology. Finally, several perspectives for future research are outlined

Predators and parasitoids of egg sacs of the widow spiders, Latrodectus geometricus and Latrodectus hesperus (Araneae: Theridiidae) in southern California

Volume: 40Start Page: 209End Page: 21