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

Gene profiling of the erythro- and megakaryoblastic leukaemias induced by the Graffi murine retrovirus

<p>Abstract</p> <p>Background</p> <p>Acute erythro- and megakaryoblastic leukaemias are associated with very poor prognoses and the mechanism of blastic transformation is insufficiently elucidated. The murine Graffi leukaemia retrovirus induces erythro- and megakaryoblastic leukaemias when inoculated into NFS mice and represents a good model to study these leukaemias.</p> <p>Methods</p> <p>To expand our understanding of genes specific to these leukaemias, we compared gene expression profiles, measured by microarray and RT-PCR, of all leukaemia types induced by this virus.</p> <p>Results</p> <p>The transcriptome level changes, present between the different leukaemias, led to the identification of specific cancerous signatures. We reported numerous genes that may be potential oncogenes, may have a function related to erythropoiesis or megakaryopoiesis or have a poorly elucidated physiological role. The expression pattern of these genes has been further tested by RT-PCR in different samples, in a Friend erythroleukaemic model and in human leukaemic cell lines.</p> <p>We also screened the megakaryoblastic leukaemias for viral integrations and identified genes targeted by these integrations and potentially implicated in the onset of the disease.</p> <p>Conclusions</p> <p>Taken as a whole, the data obtained from this global gene profiling experiment have provided a detailed characterization of Graffi virus induced erythro- and megakaryoblastic leukaemias with many genes reported specific to the transcriptome of these leukaemias for the first time.</p

A Modifier Screen for Bazooka/PAR-3 Interacting Genes in the Drosophila Embryo Epithelium

The development and homeostasis of multicellular organisms depends on sheets of epithelial cells. Bazooka (Baz; PAR-3) localizes to the apical circumference of epithelial cells and is a key hub in the protein interaction network regulating epithelial structure. We sought to identify additional proteins that function with Baz to regulate epithelial structure in the Drosophila embryo.The baz zygotic mutant cuticle phenotype could be dominantly enhanced by loss of known interaction partners. To identify additional enhancers, we screened molecularly defined chromosome 2 and 3 deficiencies. 37 deficiencies acted as strong dominant enhancers. Using deficiency mapping, bioinformatics, and available single gene mutations, we identified 17 interacting genes encoding known and predicted polarity, cytoskeletal, transmembrane, trafficking and signaling proteins. For each gene, their loss of function enhanced adherens junction defects in zygotic baz mutants during early embryogenesis. To further evaluate involvement in epithelial polarity, we generated GFP fusion proteins for 15 of the genes which had not been found to localize to the apical domain previously. We found that GFP fusion proteins for Drosophila ASAP, Arf79F, CG11210, Septin 5 and Sds22 could be recruited to the apical circumference of epithelial cells. Nine of the other proteins showed various intracellular distributions, and one was not detected.Our enhancer screen identified 17 genes that function with Baz to regulate epithelial structure in the Drosophila embryo. Our secondary localization screen indicated that some of the proteins may affect epithelial cell polarity by acting at the apical cell cortex while others may act through intracellular processes. For 13 of the 17 genes, this is the first report of a link to baz or the regulation of epithelial structure

Comparisons between Chemical Mapping and Binding to Isoenergetic Oligonucleotide Microarrays Reveal Unexpected Patterns of Binding to the Bacillus subtilis RNase P RNA Specificity Domain†

ABSTRACT: Microarrays with isoenergetic pentamer and hexamer 20-O-methyl oligonucleotide probes with LNA (locked nucleic acid) and 2,6-diaminopurine substitutions were used to probe the binding sites on theRNase P RNA specificity domain of Bacillus subtilis. Unexpected binding patterns were revealed. Because of their enhanced binding free energies, isoenergetic probes can break short duplexes, merge adjacent loops, and/or induce refolding. This suggests new approaches to the rational design of short oligonucleotide therapeutics but limits the utility of microarrays for providing constraints for RNA structure determination. The microarray results are compared to results from chemical mapping experiments, which do provide constraints. Results from both types of experiments indicate that the RNase P RNA folds similarly in 1MNaþ and 10 mMMg2þ. Binding of RNA to RNA is important for many natural func-tions, includingproteinsynthesis (1,2), translationregulation (3,4), gene silencing (5, 6), metabolic regulation (7), RNAmodification (8, 9), etc. (10-13). Binding of oligonucleotides toRNAs is impor-tant for therapeutic approaches, such as siRNA, ribozymes, and antisense therapy (14, 15).Much remains to bediscovered, however, of the rules for predicting binding sites andpotential therapeutics