Browsing by Author "Rajasimha, Harsha Karur"

Now showing 1 - 2 of 2
  • Thumbnail Image
    Insights Into Mitochondrial Genetic and Morphologic Dynamics Gained by Stochastic Simulation
    Rajasimha, Harsha Karur (Virginia Tech, 2007-12-05)
    MtDNA mutations in mammalian cells are implicated in cellular ageing and encephalomyopathies, although mechanisms involved are not completely understood. The mitochondrial genetic bottleneck has puzzled biologists for a long time. Approximate models of genetic bottleneck proposed in the literature do not accurately model underlying biology. Recent studies indicate mitochondrial morphology changes during cellular aging in culture. In particular, the rates of mitochondrial fission and fusion are shown to be in tight balance, though this rate decreases with age. Some proteins involved in mitochondrial morphology maintenance are implicated in apoptosis. Hence, mitochondrial genetic and morphologic dynamics are critical to the life and death of cells. By working closely with experimental collaborators and by utilizing data derived from literature, we have developed stochastic simulation models of mitochondrial genetic and morphologic dynamics. Hypotheses from the mitochondrial genetic dynamics model include: (1) the decay of mtDNA heteroplasmy in blood is exponential and not linear as reported in literature. (2) Blood heteroplasmy measurements are a good proxy for the blood stem cell heteroplasmy. (3) By analyzing our simulation results in tandem with published longitudinal clinical data, we propose for the first time, a way to correct for the patient's age in the analysis of heteroplasmy data. (4) We develop a direct model of the genetic bottleneck process during mouse embryogenesis. (5) Partitioning of mtDNA into daughter cells during blastocyst formation and relaxed replication of mtDNA during the exponential growth phase of primordial germ cells leads to the variation in heteroplasmy inherited by offspring from the same mother. (6) We develop a “simulation control” for experimental studies on mtDNA heteroplasmy variation in cell cultures. Hypothesis from the mitochondrial morphologic dynamics model: (7) A cell adjusts the mitochondrial fusion rate to compensate for the fluctuations in the fission rate, but not vice versa. A deterministic model for this control is proposed. Contributions: extensible simulation models of mitochondrial genetic and morphologic dynamics to aide in the powerful analysis of published and new experimental data. Our results have direct relevance to cell biology and clinical diagnosis. The work also illustrates scientific success by tight integration of theory with practice.
  • Thumbnail Image
    PathMeld: A Methodology for The Unification of Metabolic Pathway Databases
    Rajasimha, Harsha Karur (Virginia Tech, 2004-12-15)
    A biological pathway database is a database that describes biochemical pathways, reactions, enzymes that catalyze the reactions, and the substrates that participate in these reactions. A pathway genome database (PGDB) integrates pathway information with information about the complete genome of various sequenced organisms. Two of the popular PGDBs available today are the Kyoto Encyclopedia of Genes and Genomes (KEGG) and MetaCyc. The proliferation of biological databases in general raises several questions for the life scientist. Which of these databases is most accurate, most current, or most comprehensive? Do they have a standard format? Do they complement each other? Overall, which database should be used for what purpose? If more than one database is deemed relevant, it is desirable to have a unified database containing information from all the shortlisted databases. There is no standard methodology yet for integrating biological pathway databases and, to the best of our knowledge, no commercial software that can perform such integration tasks. While XML based pathway data exchange standards such as BioPAX and SBML are emerging, these do not address the basic problems such as inconsistent nomenclature and substrate matching between databases in the unification of pathway databases. Here, we present the PathMeld methodology to unify KEGG and MetaCyc databases starting from their flat files. Individual PGDBs are transformed into a unified schema that we design. With individual PGDBs in the common unified schema, the key to the PathMeld methodology is to find the entity correspondences between the KEGG and MetaCyc substrates. We present a heuristic driven approach for one-to-one mapping of the substrates between KEGG and MetaCyc. Using the exact name and chemical formula match criteria, 82.6% of the substrates in MetaCyc were matched accurately to corresponding substrates in KEGG. The substrate names in the MetaCyc database contain html tags and non-characters such as , , , , &, and $. The MetaCyc chemical formula are stored in lisp format in the database while KEGG stores them as continuous strings. Hence, we subject MetaCyc chemical formulae to transformation into KEGG format to make them directly comparable. Applying pre-processing to transform MetaCyc substrate names and formulae improved substrate matching by 2%. To investigate how many of the remaining 17:4% substrates are indeed absent from KEGG, we employ a standard UNIX based approximate string matching tool called agrep. The resulting matches are curated into four mutually exlusive groups: 3:83% are correct matches, 3:17% are close matches, and 7:45% are incorrect matches. 3:68% of MetaCyc substrate names are not matched at all. This shows that 11:13% of MetaCyc substrate names are absent in KEGG. We note some of the implementation issues we solved. First, parsing only one flat file to populate one database table is not sufficient. Second, intermediate database tables are needed. Third, transformation of substrate names and chemical formula from one of the component databases is required for comparison. Fourth, a biochemist's intervention is needed in evaluating the approximate substrate matches from agrep. In conclusion, the PathMeld methodology successfully uni¯es KEGG and MetaCyc °at ¯le databases into a uni¯ed PostgreSQL database. Matching substrates between databases is the key issue in the uni¯cation process. About 83% of the substrate correspondences can be computationally achieved, while the remaining 17% substrates require approximate matching and manual curation by a biochemist. We presented several di®erent techniques for substrate matching and showed that about 10% of the MetaCyc substrates do not match and hence are absent from KEGG.