I have over 15 years of experience in computational biology. My research goal is to predict the dynamic behaviour of the living cell by computer simulation of the genome scale network models representing experimental data on interaction between molecules.
I am convinced that we can fully exploit information about full genomic sequence of human and other organisms only if we use legacy of molecular biology data to build predictive mechanistic models of genotype-phenotype relationship. Due to the number of molecular components in the cell and non-linearity of their interactions this goal can only be achieved by computer simulation. The successful computer simulation of the molecular cell biology will enable prediction of the individual genetic differences on the trajectories of major diseases providing foundation for predictive and personalized medicine of the future. Likewise, industrial biotechnology is being revolutionized by increasing ability to computer simulate the effects of genetic engineering in commercial cell lines and therefore rationally design industrial fermentation processes.I developed novel Quasi Steady State Petri Net (QSSPN) algorithm (Bioinformatics, doi: 10.1093/bioinformatics/btt552) integrating Petri nets and constraint-based analysis to predict the feasibility of qualitative dynamic behaviours in models of gene regulation, signalling and whole-cell metabolism. We presented the first dynamic simulations including regulatory mechanisms and a genome scale metabolic network in human cell, using bile acid homeostasis in human hepatocytes as a case study. QSSPN simulations reproduced experimentally determined qualitative dynamic behaviours and permitted mechanistic analysis of genotype-phenotype relationships. (QSSPN Project website: http://sysbio3.fhms.surrey.ac.uk/qsspn/index.html)
I have performed computational part of the project leading to the first reconstruction of the Genome Scale Metabolic Reaction Network of Mycobacterium tuberculosis, causative agent of Tuberculosis disease (Genome Biology, 2007). The tools developed for this project motivated have been matured into SurreyFBA software recently published by my group (Bioinformatics, 2011). I have also been working on analysis of gene expression data in the context of genome scale metabolic networks (PLoS Computational Biology, 2011) and development of software for web based computation with FBA models (BMC Bioinformatics, 2011). Industrial biotechnology is an important application area for genome scale metabolic modeling; I worked on FBA simulations in the context of bioprocess feed development for antibiotic production in Streptomyces coelicolor (Metabolic Engineering, 2008).
I have been modeling stochastic effects in molecular interaction network dynamics for 10 years. I have constructed detailed model of prokaryotic gene expression and investigated dependence between accuracy of gene expression and transcription and translation initiation rates (J. Biol. Chem, 2001). This work has also lead to the publication of STOCKS software for stochastic simulation of molecular interaction network (Bioinformatics, 2002). Subsequently, we have developed Maximal Timestep Method, a hybrid algorithm enabling stochastic simulation of systems with reaction rates varying by many orders of magnitude. The method has been applied to investigate propagation of gene expression noise to the level of metabolic processes leading to epigenetically inherited changes in single cell physiology (Biophysical Journal 2004). More recently, I was working the influence of RNA regulators on gene expression noise (Biophysical Journal 2009) and constructed stochastic kinetic model of Two Component System Signalling (Molecular Biosystems 2010).
In have past bioinformatics experience in the field of homology modeling of protein structure (Nucl. Acids. Research 2003, Nature Immunology 2003), regulatory sequence analysis (J. Biol. Chem 2005) and annotation of genome sequences (Nature 2004). I did my PhD in the area of Biophysics and worked on the agent-based simulations of protein crystal growth (Biophysical Journal 1997). I have also performed molecular dynamics simulations and analysed light scaterring spectra (J. Phys. Chem. 1999).
1. BMS3072: Systems Biology: Genomes in Action
2. BMS1023: Numeracy skills and Statistics
4. MSc Coourses: Statistics
Professor of Systems Biology
Module organiser for BMS3072
Find me on campus Room: 19 AY 04
Dynamic simulation of genome-scale molecular interaction networks will enable the mechanistic prediction of genotype-phenotype relationships. Despite advances in quantitative biology, full parameterization of whole-cell models is not yet possible. Simulation methods capable of using available qualitative data are required to develop dynamic whole-cell models through an iterative process of modelling and experimental validation.
One of the most challenging problems in microbiology is to understand how a small fraction of microbes that resists killing by antibiotics can emerge in a population of genetically identical cells, the phenomenon known as persistence or drug tolerance. Its characteristic signature is the biphasic kill curve, whereby microbes exposed to a bactericidal agent are initially killed very rapidly but then much more slowly. Here we relate this problem to the more general problem of understanding the emergence of distinct growth phenotypes in clonal populations. We address the problem mathematically by adopting the framework of the phenomenon of so-called weak ergodicity breaking, well known in dynamical physical systems, which we extend to the biological context. We show analytically and by direct stochastic simulations that distinct growth phenotypes can emerge as a consequence of slow-down of stochastic fluctuations in the expression of a gene controlling growth rate. In the regime of fast gene transcription, the system is ergodic, the growth rate distribution is unimodal, and accounts for one phenotype only. In contrast, at slow transcription and fast translation, weakly non-ergodic components emerge, the population distribution of growth rates becomes bimodal, and two distinct growth phenotypes are identified. When coupled to the well-established growth rate dependence of antibiotic killing, this model describes the observed fast and slow killing phases, and reproduces much of the phenomenology of bacterial persistence. The model has major implications for efforts to develop control strategies for persistent infections.
Streptomycetes sense and respond to the stress of phosphate starvation via the two-component PhoR-PhoP signal transduction system. To identify the in vivo targets of PhoP we have undertaken a chromatin-immunoprecipitation-on-microarray analysis of wild-type and phoP mutant cultures and, in parallel, have quantified their transcriptomes. Most (ca. 80%) of the previously in vitro characterized PhoP targets were identified in this study among several hundred other putative novel PhoP targets. In addition to activating genes for phosphate scavenging systems PhoP was shown to target two gene clusters for cell wall/extracellular polymer biosynthesis. Furthermore PhoP was found to repress an unprecedented range of pathways upon entering phosphate limitation including nitrogen assimilation, oxidative phosphorylation, nucleotide biosynthesis and glycogen catabolism. Moreover, PhoP was shown to target many key genes involved in antibiotic production and morphological differentiation, including afsS, atrA, bldA, bldC, bldD, bldK, bldM, cdaR, cdgA, cdgB and scbR-scbA. Intriguingly, in the PhoP-dependent cpk polyketide gene cluster, PhoP accumulates substantially at three specific sites within the giant polyketide synthase-encoding genes. This study suggests that, following phosphate limitation, Streptomyces coelicolor PhoP functions as a 'master' regulator, suppressing central metabolism, secondary metabolism and developmental pathways until sufficient phosphate is salvaged to support further growth and, ultimately, morphological development.
Phenotypic differences of genetically identical cells under the same environmental conditions have been attributed to the inherent stochasticity of biochemical processes. Various mechanisms have been suggested, including the existence of alternative steady states in regulatory networks that are reached by means of stochastic fluctuations, long transient excursions from a stable state to an unstable excited state, and the switching on and off of a reaction network according to the availability of a constituent chemical species. Here we analyse a detailed stochastic kinetic model of two-component system signalling in bacteria, and show that alternative phenotypes emerge in the absence of these features. We perform a bifurcation analysis of deterministic reaction rate equations derived from the model, and find that they cannot reproduce the whole range of qualitative responses to external signals demonstrated by direct stochastic simulations. In particular, the mixed mode, where stochastic switching and a graded response are seen simultaneously, is absent. However, probabilistic and equation-free analyses of the stochastic model that calculate stationary states for the mean of an ensemble of stochastic trajectories reveal that slow transcription of either response regulator or histidine kinase leads to the coexistence of an approximate basal solution and a graded response that combine to produce the mixed mode, thus establishing its essential stochastic nature. The same techniques also show that stochasticity results in the observation of an all-or-none bistable response over a much wider range of external signals than would be expected on deterministic grounds. Thus we demonstrate the application of numerical equation-free methods to a detailed biochemical reaction network model, and show that it can provide new insight into the role of stochasticity in the emergence of phenotypic diversity.
BACKGROUND: Neisseria meningitidis is an important human commensal and pathogen that causes several thousand deaths each year, mostly in young children. How the pathogen replicates and causes disease in the host is largely unknown, particularly the role of metabolism in colonization and disease. Completed genome sequences are available for several strains but our understanding of how these data relate to phenotype remains limited. RESULTS: To investigate the metabolism of N. meningitidis we generated and selected a representative Tn5 library on rich medium, a minimal defined medium and in human serum to identify genes essential for growth under these conditions. To relate these data to a systems-wide understanding of the pathogen's biology we constructed a genome-scale metabolic network: Nmb_iTM560. This model was able to distinguish essential and non-essential genes as predicted by the global mutagenesis. These essentiality data, the library and the Nmb_iTM560 model are powerful and widely applicable resources for the study of meningococcal metabolism and physiology. We demonstrate the utility of these resources by predicting and demonstrating metabolic requirements on minimal medium such as a requirement for PEP carboxylase, and by describing the nutritional and biochemical status of N. meningitidis when grown in serum, including a requirement for both the synthesis and transport of amino acids. CONCLUSIONS: This study describes the application of a genome scale transposon library combined with an experimentally validated genome-scale metabolic network of N. meningitidis to identify essential genes and provide novel insight to the pathogen's metabolism both in vitro and during infection.
Constraint-based approaches facilitate the prediction of cellular metabolic capabilities, based, in turn on predictions of the repertoire of enzymes encoded in the genome. Recently, genome annotations have been used to reconstruct genome scale metabolic reaction networks for numerous species, including Homo sapiens, which allow simulations that provide valuable insights into topics, including predictions of gene essentiality of pathogens, interpretation of genetic polymorphism in metabolic disease syndromes and suggestions for novel approaches to microbial metabolic engineering. These constraint-based simulations are being integrated with the functional genomics portals, an activity that requires efficient implementation of the constraint-based simulations in the web-based environment.
A large number of cDNA inserts were sequenced from a high-quality library of chicken bursal lymphocyte cDNAs. Comparisons to public gene databases indicate that the cDNA collection represents more than 2,000 new, full-length transcripts. This resource defines the structure and the coding potential of a large fraction of B-cell specific and housekeeping genes whose function can be analyzed by disruption in the chicken DT40 B-cell line.
Page Owner: bss1ak
Page Created: Friday 31 July 2009 11:38:48 by t00345
Last Modified: Monday 28 April 2014 10:04:08 by kj0008
Assembly date: Wed Aug 27 21:49:25 BST 2014
Content ID: 11033