Differential dynamic properties of scleroderma fibroblasts in response to perturbation of environmental stimuli.

Diseases are believed to arise from dysregulation of biological systems (pathways) perturbed by environmental triggers. Biological systems as a whole are not just the sum of their components, rather ever-changing, complex and dynamic systems over time in response to internal and external perturbation. In the past, biologists have mainly focused on studying either functions of isolated genes or steady-states of small biological pathways. However, it is systems dynamics that play an essential role in giving rise to cellular function/dysfunction which cause diseases, such as growth, differentiation, division and apoptosis. Biological phenomena of the entire organism are not only determined by steady-state characteristics of the biological systems, but also by intrinsic dynamic properties of biological systems, including stability, transient-response, and controllability, which determine how the systems maintain their functions and performance under a broad range of random internal and external perturbations. As a proof of principle, we examine signal transduction pathways and genetic regulatory pathways as biological systems. We employ widely used state-space equations in systems science to model biological systems, and use expectation-maximization (EM) algorithms and Kalman filter to estimate the parameters in the models. We apply the developed state-space models to human fibroblasts obtained from the autoimmune fibrosing disease, scleroderma, and then perform dynamic analysis of partial TGF-beta pathway in both normal and scleroderma fibroblasts stimulated by silica. We find that TGF-beta pathway under perturbation of silica shows significant differences in dynamic properties between normal and scleroderma fibroblasts. Our findings may open a new avenue in exploring the functions of cells and mechanism operative in disease development.

Differential dynamic properties of scleroderma fibroblasts in response to perturbation of environmental stimuli.

Diseases are believed to arise from dysregulation of biological systems (pathways) perturbed by environmental triggers. Biological systems as a whole are not just the sum of their components, rather ever-changing, complex and dynamic systems over time in response to internal and external perturbation. In the past, biologists have mainly focused on studying either functions of isolated genes or steady-states of small biological pathways. However, it is systems dynamics that play an essential role in giving rise to cellular function/dysfunction which cause diseases, such as growth, differentiation, division and apoptosis. Biological phenomena of the entire organism are not only determined by steady-state characteristics of the biological systems, but also by intrinsic dynamic properties of biological systems, including stability, transient-response, and controllability, which determine how the systems maintain their functions and performance under a broad range of random internal and external perturbations. As a proof of principle, we examine signal transduction pathways and genetic regulatory pathways as biological systems. We employ widely used state-space equations in systems science to model biological systems, and use expectation-maximization (EM) algorithms and Kalman filter to estimate the parameters in the models. We apply the developed state-space models to human fibroblasts obtained from the autoimmune fibrosing disease, scleroderma, and then perform dynamic analysis of partial TGF-beta pathway in both normal and scleroderma fibroblasts stimulated by silica. We find that TGF-beta pathway under perturbation of silica shows significant differences in dynamic properties between normal and scleroderma fibroblasts. Our findings may open a new avenue in exploring the functions of cells and mechanism operative in disease development.

Identification of diabetic retinopathy genes through a genome-wide association study among Mexican-Americans from Starr County, Texas

To identify genetic susceptibility loci for severe diabetic retinopathy, 286 Mexican-Americans with type 2 diabetes from Starr County, Texas completed detailed physical and ophthalmologic examinations including fundus photography for diabetic retinopathy grading. 103 individuals with moderate-to-severe non-proliferative diabetic retinopathy or proliferative diabetic retinopathy were defined as cases for this study. DNA samples extracted from study subjects were genotyped using the Affymetrix GeneChip® Human Mapping 100K Set, which includes 116,204 single nucleotide polymorphisms (SNPs) across the whole genome. Single-marker allelic tests and 2- to 8-SNP sliding-window Haplotype Trend Regression implemented in HelixTreeTM were first performed with these direct genotypes to identify genes/regions contributing to the risk of severe diabetic retinopathy. An additional 1,885,781 HapMap Phase II SNPs were imputed from the direct genotypes to expand the genomic coverage for a more detailed exploration of genetic susceptibility to diabetic retinopathy. The average estimated allelic dosage and imputed genotypes with the highest posterior probabilities were subsequently analyzed for associations using logistic regression and Fisher's Exact allelic tests, respectively. To move beyond these SNP-based approaches, 104,572 directly genotyped and 333,375 well-imputed SNPs were used to construct genetic distance matrices based on 262 retinopathy candidate genes and their 112 related biological pathways. Multivariate distance matrix regression was then used to test hypotheses with genes and pathways as the units of inference in the context of susceptibility to diabetic retinopathy. This study provides a framework for genome-wide association analyses, and implicated several genes involved in the regulation of oxidative stress, inflammatory processes, histidine metabolism, and pancreatic cancer pathways associated with severe diabetic retinopathy. Many of these loci have not previously been implicated in either diabetic retinopathy or diabetes. In summary, CDC73, IL12RB2, and SULF1 had the best evidence as candidates to influence diabetic retinopathy, possibly through novel biological mechanisms related to VEGF-mediated signaling pathway or inflammatory processes. While this study uncovered some genes for diabetic retinopathy, a comprehensive picture of the genetic architecture of diabetic retinopathy has not yet been achieved. Once fully understood, the genetics and biology of diabetic retinopathy will contribute to better strategies for diagnosis, treatment and prevention of this disease.^

Pathway analysis of genome-wide association study data in primary biliary cirrhosisa

Genome-wide association studies (GWAS) have successfully identified several genetic loci associated with inherited predisposition to primary biliary cirrhosis (PBC), the most common autoimmune disease of the liver. Pathway-based tests constitute a novel paradigm for GWAS analysis. By evaluating genetic variation across a biological pathway (gene set), these tests have the potential to determine the collective impact of variants with subtle effects that are individually too weak to be detected in traditional single variant GWAS analysis. To identify biological pathways associated with the risk of development of PBC, GWAS of PBC from Italy (449 cases and 940 controls) and Canada (530 cases and 398 controls) were independently analyzed. The linear combination test (LCT), a recently developed pathway-level statistical method was used for this analysis. For additional validation, pathways that were replicated at the P <0.05 level of significance in both GWAS on LCT analysis were also tested for association with PBC in each dataset using two complementary GWAS pathway approaches. The complementary approaches included a modification of the gene set enrichment analysis algorithm (i-GSEA4GWAS) and Fisher's exact test for pathway enrichment ratios. Twenty-five pathways were associated with PBC risk on LCT analysis in the Italian dataset at P<0.05, of which eight had an FDR<0.25. The top pathway in the Italian dataset was the TNF/stress related signaling pathway (p=7.38×10 -4, FDR=0.18). Twenty-six pathways were associated with PBC at the P<0.05 level using the LCT in the Canadian dataset with the regulation and function of ChREBP in liver pathway (p=5.68×10-4, FDR=0.285) emerging as the most significant pathway. Two pathways, phosphatidylinositol signaling system (Italian: p=0.016, FDR=0.436; Canadian: p=0.034, FDR=0.693) and hedgehog signaling (Italian: p=0.044, FDR=0.636; Canadian: p=0.041, FDR=0.693), were replicated at LCT P<0.05 in both datasets. Statistically significant association of both pathways with PBC genetic susceptibility was confirmed in the Italian dataset on i-GSEA4GWAS. Results for the phosphatidylinositol signaling system were also significant in both datasets on applying Fisher's exact test for pathway enrichment ratios. This study identified a combination of known and novel pathway-level associations with PBC risk. If functionally validated, the findings may yield fresh insights into the etiology of this complex autoimmune disease with possible preventive and therapeutic application.^

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.

SYNTHESIS AND BIOLOGICAL EVALUATION OF ANTITUMOR ALDOPHOSPHAMIDE ACETALS

Although the major metabolic pathways of cyclophosphamide are well established, the mechanism of antitumor drug selectivity is highly controversial. However, it is widely accepted that aldophosphamide, one of the primary metabolites, plays a crucial role in drug selectivity. In an attempt to gain a better understanding of the mechanism of selectivity of cyclophosphamide, a series of aldophosphamide analogs have been synthesized.^ The new analogs, unlike aldophosphamide, are relatively stable in neutral solution; however, they are converted rapidly to aldehydo intermediates in the presence of carboxylate esterase. Due to structural differences, these analogs may be classified into three different groups, arbitrarily designated as A, B, C, depending upon the facility with which the intermediate aldehydes form 4-hydroxy cyclic tautomers. The half-life of the aldehydo/4-hydroxy cyclic tautomeric mixture is longer for bis(acetoxy)aldophosphamide acetal I (a representative of group A), shorter for the n-ethyl analog III (B), and shortest for the N,N-dimethyl analog IV (C). The ratio of aldophosphamide: 4-hydroxycyclophosphamide at pseudoequilibrium is 1: 4 for compound I, 1: 2 for compound III and 0: 1 for compound IV. The therapeutic efficacy of these compounds are group A $>$ group B $>$ group C. It is apparent that the equilibrium position between the aldehydo and 4-hydroxy cyclic tautomers, which determines their stability, is a crucial determinant of both the cytotoxicity and antitumor selectivity. These findings, taken in conjunction with the aldehyde dehydrogenase selectivity hypothesis, may provide an explanation for the unique antitumor activity of cyclophosphamide. ^

Studies of MEKK2 activation mechanisms and its biological role in LPS/TLR4 signaling

MEKK2 is an evolutionarily conserved mitogen-activated protein kinase (MAPK) kinase kinase (MAP3K) that controls the MAPK and IKK-NF-κB pathways. The MAPK and IKK pathways are intracellular signaling networks that are crucial for the Toll-like receptor (TLR) mediated innate immunity, cellular stress and many other physiological responses. Members of the MAP3K family are central to the activation of these processes. However, the molecular mechanisms underlying stimuli-mediated MAP3K activation remain largely unknown. In this study, we identified a key phosphoserine residue, Ser-519 in MEKK2, and its equivalent site Ser-526 in MEKK3 within their activation loop whose phosphorylation are essential for their optimal activation. Mutation of this regulatory serine to an alanine severely impaired MEKK2 activation and MEKK2 signaling to its downstream targets. To demonstrate that physiological stimuli induce this serine phosphorylation, we generated an antibody that specifically recognizes the phosphorylated serine residue. We found that many, but not all, of the MAPK agonists, including the TLR ligands, growth factors, cytokines and cellular stresses, induced this regulatory serine phosphorylation in MEKK2, suggesting an involvement of MEKK2 in the activation of the MAPK cascade leading to different cellular responses. We further investigated the specific role of MEKK2 in LPS/TLR4 signaling by using MEKK2−/− mice. We found that MEKK2 was selectively required for LPS-induced ERK1/2 activation, but not JNK, p38 or NF-κB activation. We also found that MEKK2 was involved in TLR4 dependent induction of proinflammatory cytokines and LPS-induced septic shock. In conclusion, we identified a key regulatory serine residue in the activation loop of MEKK2 whose phosphorylation is a key sensor of receptor- and cellular stress-mediated signals. We also demonstrated that MEKK2 is crucial for TLR4-mediated innate immunity. ^

p21 activated kinase 5 links multiple GTPase signaling pathways at mitochondria

The p21-activated kinase 5 (PAK5) is a serine/threonine protein kinase associated with the group 2 subfamily of PAKs. Although our understanding about PAK5 is very limited, it is receiving increasing interest due to its tissue specific expression pattern and important signaling properties. PAK5 is highly expressed in brain. Its overexpression induces neurite outgrowth in neuroblastoma cells and promotes survival in fibroblasts. ^ The serine/threonine protein kinase Raf-1 is an essential mediator of Ras-dependent signaling that controls the ERK/MAPK pathway. In contrast to PAK5, Raf-1 has been the subject of intensive investigation. However due to the complexity of its activation mechanism, the biological inputs controlling Raf-1 activation are not fully understood. ^ PAKs 1-3 are the known kinases responsible for phosphorylation of Raf-1 on serine 338, which is a crucial phosphorylation site for Raf-1 activation. However, dominant negative versions of these kinases do not block EGF-induced Raf-1 activation, indicating that other kinases may regulate the phosphorylation of Raf-1 on serine 338. ^ This thesis work was initiated to test whether the group 2 PAKs 4, 5 and 6 are responsible for EGF-induced Raf-1 activation. We found that PAK5, and to a lesser extent PAK4, can activate Raf-1 in cells. Our studies thereafter focused on PAK5. With the progress of our study we found that PAK5 does not significantly stimulate serine 338 phosphorylation of Triton X-100 soluble Raf-1. PAK5, however, constitutively and specifically associates with Raf-1 and targets it to a Triton X-100 insoluble, mitochondrial compartment, where PAK5 phosphorylates serine 338 of Raf-1. We further demonstrated that endogenous PAK5 and Raf-1 colocalize in Hela cells at the mitochondrial outer membrane. In addition, we found that the mitochondria-targeting of PAK5 is determined by its C-terminal kinase domain plus the upstream proximal region, and facilitated by the N-terminal p21 binding domain. We also demonstrated that Rho GTPases Cdc42 and RhoD associate with and regulate the subcellular localization of PAK5. Taken together, this work suggests that the mitochondria-targeting of PAK5 may link Ras and Rho GTPase-mediated signaling pathways, and sheds light on aspects of PAK5 signaling that may be important for regulating neuronal homeostasis. ^