Pathway Semantics: An Algebraic Data Driven Algorithm to Generate Hypotheses about Molecular Patterns Underlying Disease Progression

The overarching goal of the Pathway Semantics Algorithm (PSA) is to improve the in silico identification of clinically useful hypotheses about molecular patterns in disease progression. By framing biomedical questions within a variety of matrix representations, PSA has the flexibility to analyze combined quantitative and qualitative data over a wide range of stratifications. The resulting hypothetical answers can then move to in vitro and in vivo verification, research assay optimization, clinical validation, and commercialization. Herein PSA is shown to generate novel hypotheses about the significant biological pathways in two disease domains: shock / trauma and hemophilia A, and validated experimentally in the latter. The PSA matrix algebra approach identified differential molecular patterns in biological networks over time and outcome that would not be easily found through direct assays, literature or database searches. In this dissertation, Chapter 1 provides a broad overview of the background and motivation for the study, followed by Chapter 2 with a literature review of relevant computational methods. Chapters 3 and 4 describe PSA for node and edge analysis respectively, and apply the method to disease progression in shock / trauma. Chapter 5 demonstrates the application of PSA to hemophilia A and the validation with experimental results. The work is summarized in Chapter 6, followed by extensive references and an Appendix with additional material.

Cellular dynamic simulator: an event driven molecular simulation environment for cellular physiology.

In this paper, we present the Cellular Dynamic Simulator (CDS) for simulating diffusion and chemical reactions within crowded molecular environments. CDS is based on a novel event driven algorithm specifically designed for precise calculation of the timing of collisions, reactions and other events for each individual molecule in the environment. Generic mesh based compartments allow the creation / importation of very simple or detailed cellular structures that exist in a 3D environment. Multiple levels of compartments and static obstacles can be used to create a dense environment to mimic cellular boundaries and the intracellular space. The CDS algorithm takes into account volume exclusion and molecular crowding that may impact signaling cascades in small sub-cellular compartments such as dendritic spines. With the CDS, we can simulate simple enzyme reactions; aggregation, channel transport, as well as highly complicated chemical reaction networks of both freely diffusing and membrane bound multi-protein complexes. Components of the CDS are generally defined such that the simulator can be applied to a wide range of environments in terms of scale and level of detail. Through an initialization GUI, a simple simulation environment can be created and populated within minutes yet is powerful enough to design complex 3D cellular architecture. The initialization tool allows visual confirmation of the environment construction prior to execution by the simulator. This paper describes the CDS algorithm, design implementation, and provides an overview of the types of features available and the utility of those features are highlighted in demonstrations.

THEORETICAL STUDIES ON THE METHODS OF RECONSTRUCTING PHYLOGENETIC TREES FROM DNA SEQUENCE DATA (MOLECULAR EVOLUTION, HOMINOID EVOLUTION, COMPUTER SIMULATION)

(1) A mathematical theory for computing the probabilities of various nucleotide configurations is developed, and the probability of obtaining the correct phylogenetic tree (model tree) from sequence data is evaluated for six phylogenetic tree-making methods (UPGMA, distance Wagner method, transformed distance method, Fitch-Margoliash's method, maximum parsimony method, and compatibility method). The number of nucleotides (m*) necessary to obtain the correct tree with a probability of 95% is estimated with special reference to the human, chimpanzee, and gorilla divergence. m* is at least 4,200, but the availability of outgroup species greatly reduces m* for all methods except UPGMA. m* increases if transitions occur more frequently than transversions as in the case of mitochondrial DNA. (2) A new tree-making method called the neighbor-joining method is proposed. This method is applicable either for distance data or character state data. Computer simulation has shown that the neighbor-joining method is generally better than UPGMA, Farris' method, Li's method, and modified Farris method on recovering the true topology when distance data are used. A related method, the simultaneous partitioning method, is also discussed. (3) The maximum likelihood (ML) method for phylogeny reconstruction under the assumption of both constant and varying evolutionary rates is studied, and a new algorithm for obtaining the ML tree is presented. This method gives a tree similar to that obtained by UPGMA when constant evolutionary rate is assumed, whereas it gives a tree similar to that obtained by the maximum parsimony tree and the neighbor-joining method when varying evolutionary rate is assumed. ^

STUDIES INTO THE MOLECULAR MECHANISMS REGULATING MUSCULAR ATROPHY RESULTING FROM HINDLIMB IMMOBILIZATION IN THE RAT (PROTEIN, SYNTHESIS, RNA)

The purpose of the work performed in this dissertation was to examine some of the possible regulatory mechanisms involved in the initiation of muscular atrophy during periods of decreased muscle utilization resulting from hindlimb immobilization in the rat. A 37% decrease in the rate of total muscle protein synthesis which has been observed to occur in the first 6 h of immobilization contributes significantly to the observed loss of protein during immobilization.^ The rates of cytochrome c and actin synthesis were determined in adult rat red vastus lateralis and gastrocnemius muscles, respectively, by the constant infusion and incorporation of ('3)H-tyrosine into protein. The fractional synthesis rates of both actin and cytochrome c were significantly decreased (P < 0.05) in the 6th h of hindlimb immobilization.^ RHA was extracted from adult rat gastrocnemius muscle by modification of the phenol: chloroform: SDS extraction procedures commonly used for preparation of RNA for hybridization analysis from other mammalian tissues. RNA content of rat gastrocnemius muscle, as determined by this method of extraction and its subsequent quantification by UV absorbance and orcinol assay, was significantly greater than the RNA content previously determined for adult rat gastrocnemius by other commonly employed methods.^ RNA extracted by this method from gastrocnemius muscles of control and 6h immobilized rats was subjected to "dot blot" hybridization to ('32)P-labelled probe from plasmid p749, containing a cDNA sequence complementary to (alpha)-actin mRNA and from rat skeletal muscle. (alpha)-Actin specific mRNA content as estimated by this procedure is not significantly decreased in rat gastrocnemius following 6h or hindlimb immobilization. However, (alpha)-actin specific mRNA content is significantly decreased (P < 0.05) in adult rat gastrocnemius (alpha)-actin specific mRNA is not decreased in adult rat gastrocnemius muscle following 6h of immobilization, a time when actin synthesis is significantly decreased, it is concluded that a change in (alpha)-actin specific mRNA content is not the initiating event responsible for the early decrease in actin synthesis observed in the 6th h of immobilization. ^

Characterization of healing tissue in a tooth extraction socket in a rabbit model

The histology of healing in a tooth extraction socket has been described in many studies. The focus of research in bone biology and healing is now centered on molecular events that regulate repair of injured tissue. Rapid progress in cellular and molecular biology has resulted in identification of many signaling molecules (growth factors and cytokines) associated with formation and repair of skeletal tissues. Some of these include members of the transforming growth factor-β superfamily (including the bone morphogenetic proteins), fibroblast growth factors, platelet derived growth factors and insulin like growth factors. ^ Healing of a tooth extraction socket is a complex process involving tissue repair and regeneration. It involves chemotaxis of appropriate cells into the wound, transformation of undifferentiated mesenchymal cells to osteoprogenitor cells, proliferation and differentiation of committed bone forming cells, extracellular matrix synthesis, mineralization of osteoid, maturation and remodeling of bone. Current data suggests that these cellular events are precisely controlled and regulated by specific signaling molecules. A plethora of cytokines; have been identified and studied in the past two decades. Some of these like transforming growth factor beta (TGF-β), vascular endothelial growth factor (VEGF), platelet derived growth factor (PDGF) and fibroblast growth factors (FGFs) are well conserved proteins involved in the initial response to injury and repair in soft and hard tissue. ^ The purpose of this study was to characterize the spatial and temporal localization of TGF-βl, VEGF, PDGF-A, FGF-2 and BMP-2, and secretory IgA in a tooth extraction socket model, and evaluate correlation of spatial and temporal changes of these growth factors to histological events. The results of this study showed positive correlation of histological events to spatial and temporal localization of TGF-β1, BMP-2, FGF-2, PDGF-A, and VEGF in a rabbit tooth extraction model. ^