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.

Extended kalman filter for estimation of parameters in nonlinear state-space models of biochemical networks.

It is system dynamics that determines the function of cells, tissues and organisms. To develop mathematical models and estimate their parameters are an essential issue for studying dynamic behaviors of biological systems which include metabolic networks, genetic regulatory networks and signal transduction pathways, under perturbation of external stimuli. In general, biological dynamic systems are partially observed. Therefore, a natural way to model dynamic biological systems is to employ nonlinear state-space equations. Although statistical methods for parameter estimation of linear models in biological dynamic systems have been developed intensively in the recent years, the estimation of both states and parameters of nonlinear dynamic systems remains a challenging task. In this report, we apply extended Kalman Filter (EKF) to the estimation of both states and parameters of nonlinear state-space models. To evaluate the performance of the EKF for parameter estimation, we apply the EKF to a simulation dataset and two real datasets: JAK-STAT signal transduction pathway and Ras/Raf/MEK/ERK signaling transduction pathways datasets. The preliminary results show that EKF can accurately estimate the parameters and predict states in nonlinear state-space equations for modeling dynamic biochemical networks.

Altered transmission of HOX and apoptotic SNPs identify a potential common pathway for clubfoot.

Clubfoot is a common birth defect that affects 135,000 newborns each year worldwide. It is characterized by equinus deformity of one or both feet and hypoplastic calf muscles. Despite numerous study approaches, the cause(s) remains poorly understood although a multifactorial etiology is generally accepted. We considered the HOXA and HOXD gene clusters and insulin-like growth factor binding protein 3 (IGFBP3) as candidate genes because of their important roles in limb and muscle morphogenesis. Twenty SNPs from the HOXA and HOXD gene clusters and 12 SNPs in IGFBP3 were genotyped in a sample composed of non-Hispanic white and Hispanic multiplex and simplex families (discovery samples) and a second sample of non-Hispanic white simplex trios (validation sample). Four SNPs (rs6668, rs2428431, rs3801776, and rs3779456) in the HOXA cluster demonstrated altered transmission in the discovery sample, but only rs3801776, located in the HOXA basal promoter region, showed altered transmission in both the discovery and validation samples (P = 0.004 and 0.028). Interestingly, HOXA9 is expressed in muscle during development. An SNP in IGFBP3, rs13223993, also showed altered transmission (P = 0.003) in the discovery sample. Gene-gene interactions were identified between variants in HOXA, HOXD, and IGFBP3 and with previously associated SNPs in mitochondrial-mediated apoptotic genes. The most significant interactions were found between CASP3 SNPS and variants in HOXA, HOXD, and IGFBP3. These results suggest a biologic model for clubfoot in which perturbation of HOX and apoptotic genes together affect muscle and limb development, which may cause the downstream failure of limb rotation into a plantar grade position.

Xenobiotic Metabolism Genes and Clubfoot

Idiopathic or isolated clubfoot is a common orthopedic birth defect that affects approximately 135,000 children worldwide. It is characterized by equinus, varus and adductus deformities of the ankle and foot. Correction of clubfoot involves months of serial manipulations, castings and bracing, with surgical correction needed in forty percent of cases. Multifactorial etiology has been suggested in numerous studies with both environmental and genetic factors playing an etiologic role. Maternal smoking during pregnancy is the only common environmental factor that has consistently been shown to increase the risk for clubfoot. Moreover, a positive family history of clubfoot and maternal smoking increases the risk of clubfoot twenty fold. These findings suggest that genetic variation in smoking metabolism genes may increase susceptibility to clubfoot. Based on this reasoning, we interrogated eight candidate genes, chosen based on their involvement in phase 1 and 2 cigarette smoke metabolism. Twenty-two SNPs and two null alleles in eight genes (CYP1A1, CYP1A2, CYP1B1, CYP2A6, EPHX1, NAT2, GSTM1 and GSTT1) were genotyped in a dataset composed of nonHispanic white and Hispanic multiplex and simplex families. Only one SNP in CYP1A1, rs1048943, had significantly altered transmission in the aggregate and multiplex NHW datasets (p=0.003 and p=0.009). Perturbation of CYP1A1 by rs1048943 polymorphism causes an increase in the amount of harmful, adduct forming metabolic intermediates. A significant gene interaction between EPHX1 and NAT2 was also found (p=0.007). This interaction may affect the metabolism of harmful metabolic intermediates. Additionally, marginal interactions were found for other xenobiotic genes and these interactions may play a contributory role in clubfoot. Importantly, for CYP1A2, significant maternal (p=0.03; RR=1.24; 95% CI: 1.04-1.44) and fetal (p=0.01; RR=1.33; 95% CI: 1.13-1.54) genotypic effects were identified suggesting that both maternal and fetal genotypes impact normal limb development. No association was found for maternal smoking status and tobacco metabolism genes. Together, these results suggest that xenobiotic metabolism genes may play a contributory role in the etiology of clubfoot regardless of maternal smoking status and may impact foot development through perturbation of tobacco metabolic pathways.

CELL POLARITY REGULATES ORGAN GROWTH THROUGH THE HIPPO PATHWAY

Defects in apical-basal cell polarity and abnormal expression of cell polarity determinants are linked to human cancer. Loss of polarity is highly correlated with malignancy. In Drosophila, perturbation of apical-basal polarity, including overexpressing the apical determinant Crumbs, can lead to uncontrolled tissue growth. Cells mutant for the basolateral determinant scribble overproliferate and can form neoplastic tumors. Interestingly, scribble mutant clones that arise in wild-type tissues are eliminated and therefore do not manifest their tumorigenic potential. However, the mechanisms by which cell polarity coordinates with growth control pathways in developing organs to achieve appropriate organ size remain obscure. To investigate the function of apical determinants in growth regulation, I investigated the mechanism by which the apical determinant Crumbs affects growth in Drosophila imaginal discs. I found that crumbs gain and loss of function cause overgrowth and induction of Hippo target genes. In addition, Crumbs is required for the proper localization of Expanded, an upstream component of the Hippo pathway. Furthermore, we uncoupled the cell polarity and growth control function of Crb through structure-functional analysis. Taken together, our data identify a role of Crb in growth regulation specifically through modulation of the Hippo pathway. To further explore the role of polarity in growth control, I investigated how cells mutant for basolateral determinants are eliminated by using patches of cells mutant for scribble (scribble mutant clones) as a model system. We found that competitive cell-cell interactions eliminate tumorigenic scribble cells by modulation of the Hippo pathway. The regulation of Hippo signaling is required and sufficient to restrain the tumorous growth of scribble mutant cells. Artificially increasing the relative fitness of scribble mutant cells unleashes their tumorigenic potential. Therefore, we have identified a novel tumor-suppression mechanism that depends on signaling between normal and tumorigenic cells. These data identify evasion of cell competition as a critical step toward malignancy and illustrate a role for wild-type tissue in eliminating abnormal cells and preventing the formation of tumors.

Gene Discovery in Nonsyndromic Cleft Lip with or without Cleft Palate

Nonsyndromic cleft lip with or without cleft palate (NSCLP), a common, complex orofacial birth defect that affects approximately 4,000 newborns each year in the United States, is caused by both genetic and environmental factors. Orofacial clefts affect the mouth and nose, causing severe deformity of the face, which require medical, dental and speech therapies. Despite having substantial genetic liability, less than 25% of the genetic contribute to NSCLP has been identified. The studies described in this thesis were performed to identify genes that contribute to NSCLP and to demonstrate the role of these genes in normal craniofacial development. Using genome scan and candidate gene approaches, novel associations with NSCLP were identified. These include MYH9 (7 SNPs, 0.009≤p<0.05), Wnt3A (4 SNPs, 0.001≤p≤0.005), Wnt11 (2 SNPs, 0.001≤p≤0.01) and CRISPLD2 (4 SNPs, 0.001≤p<0.05). The most interesting findings were for CRISPLD2. This gene is expressed in the fused mouse palate at E17.5. In zebrafish, crispld2 localized to the craniofacial region by one day post fertilization. Morpholino knockdown of crispld2 resulted in a lower survival rates and altered neural crest cell (NCC) clustering. Because NCCs form the tissues that populate the craniofacies, this NCC abnormality resulted in cartilage abnormalities of the jaw including fewer ceratobranchial cartilages forming the lower jaw (three pairs compared to five) and broader craniofacies compared to wild-type zebrafish. These findings suggest that the CRISPLD2 gene plays an important role in normal craniofacial development and perturbation of this gene in humans contributes to orofacial clefting. Overall, these results are important because they contribute to our understanding of normal craniofacial development and orofacial clefting etiology, information that can be used to develop better methods to diagnose, counsel and potentially treat NSCLP patients.

The functional and structural interactions between v-{\it mos\/} proteins and cell cycle control elements

The v-mos gene of Moloney murine sarcoma virus (Mo-MuSv) encodes a serine/threonine protein kinase capable of inducing cellular transformation. The c-mos protein is an important cell cycle regulator that functions during meiotic cell division cycles in germ cells. The overall function of c-mos in controlling meiosis is becoming better understood but the role of v-mos in malignant transformation of cells is largely unknown.^ In this study, v-mos protein was shown to be phosphorylated by M phase kinase in vitro and in vivo. The kinase activity and neoplastic transforming ability of v-mos is positively regulated by the phosphorylation. Together with the earlier finding of activation of M phase kinase by c-mos, these results raise the possibility of mutual regulation between M phase kinase and mos kinases.^ In addition to its functional interaction with the M phase kinase, the v-mos protein was shown to be present in the same protein complex with a cyclin-dependent kinase (cdk). In addition, an antibody that recognizes the cdk proteins was shown to co-precipitate the v-mos proteins in the interphase and mitotic cells transformed by p85$\sp{\rm gag-mos}$. Cdk proteins have been shown to be associated with nonmitotic cyclins which are potential oncogenes. The perturbation of cdk kinase or the activation of non-mitotic cyclins as oncogenes by v-mos could contribute directly to v-mos induced cellular transformation. v-mos proteins were also shown to interact with tubulin and vimentin, the essential components of microtubules and type IV intermediate filaments, respectively. The organizations of both microtubules and intermediate filaments are cell cycle-regulated. These results suggest that the v-mos kinase could be directly involved in inducing morphological changes typically seen in transformed cells.^ The interactions between the v-mos protein and these cell cycle control elements in regards to v-mos induced neoplastic transformation are discussed in detail in the text. ^

A new fork for clinical application: Targeting forkhead transcription factors in cancer therapy

RAS-ERK-MAPK (Mitogen-activated protein kinase) pathway plays an essential role in proliferation, differentiation, and tumor progression. In this study, we showed that ERK downregulated FOXO3a through directly interacting with and phosphorylating FOXO3a at Serine 294, Serine 344, and Serine 425. ERK-phosphorylated FOXO3a was degraded by MDM2-mediated ubiquitin-proteosome pathway. FOXO3a phosphorylation and degradation consequently promoted cell proliferation and tumorigenesis. However, the non-phosphorylated FOXO3a mutant, which was resistant to the interaction and degradation by MDM2, resulted in inhibition of tumor formation. Forkhead O transcription factors (FOXOs) are important in the regulation of cellular functions including cell cycle arrest and cell death. Perturbation of FOXOs function leads to deregulated cell proliferation and cancer. Inactivation of FOXO proteins by activation of cell survival pathways, such as PI3K/AKT/IKK, is associated with tumorigenesis. Our study will further highlight FOXOs as new therapeutic targets in a broad spectrum of cancers. ^ Chemotherapeutic drug resistance is the most concerned problem in cancer therapy as resistance ultimately leads to treatment failure of cancer patients. In another study, we showed that blocking ERK activity with AZD6244, an established MEK1/2 inhibitor currently under human cancer clinical trials, enhances FOXO3a expression in various human cancer cell lines in vitro, and also in human colon cancer cell xenografts in vivo. Knocking down FOXO3a and its downstream gene Bim impaired AZD6244-induced growth suppression, whereas restoring activation of FOXO3a sensitized human cancer cell to AZD6244-induced growth arrest and apoptosis. More importantly, AZD6244-resistant cancer cells showed impaired endogenous FOXO3a nuclear translocation, reduced FOXO3a-Bim promoter association and significantly decreased Bim expression in response to AZD6244. AZD6244-resistant cancer cells can be sensitized to API-2 (an AKT inhibitor) and LY294002 (a PI3K inhibitor) in suppressing cell growth and colony formation, these inhibitors were known to enhance FOXO3a activity/nuclear translocation through inhibiting PI3K-AKT pathway. This study reveals novel molecular mechanism contributing to AZD6244-resistance through regulation of FOXO3a activity, further provides significant clinical implication of combining AZD6244 with PI3K/AKT inhibitors for sensitizing AZD6244-resistant cancer cells by activating FOXO3a. FOXO3a activation can be an essential pharmacological target and indicator to mediate and predict AZD6244 efficacy in clinical use. ^

Elucidating the IGFBP2 signaling pathway in glioma development and progression

Diffuse gliomas are highly lethal central nervous system malignancies which, unfortunately, are the most common primary brain tumor and also the least responsive to the very few therapeutic modalities currently available to treat them. IGFBP2 is a newly recognized oncogene that is operative in multiple cancer types, including glioma, and shows promise for a targeted therapeutic approach. Elevated IGFBP2 expression is present in high-grade glioma and correlates with poor survival. We have previously demonstrated that IGFBP2 induces glioma development and progression in a spontaneous glioma mouse model, which highlighted its significance and potential for future therapy. However, we did not yet know the key physiological pathways associated with this newly characterized oncogene. We first evaluated human glioma genomics data harnessed from the publicly available Rembrandt source to identify major pathways associated with IGFBP2 expression. Integrin and ILK, among other cell migration and invasion-related pathways, were the most prominently associated. We confirmed that these pathways are regulated by IGFBP2 in glioma cells lines, and demonstrated that 1) IGFBP2 activates integrin α5β1, leading to the activation of key pathways important in glioma; 2) IGFBP2 mediates cell migration pathways through ILK; and 3) IGFBP2 activates NF-kB via an integrin α5 interaction. We then sought to determine whether this was a physiologically active signaling pathway in vivo by assessing its ability to induce glioma progression in the RCAS/tv-a spontaneous glioma mouse model. We found that ILK is a key downstream mediator of IGFBP2 that is required for the induction of glioma progression. Most significantly, a genetic therapeutic approach revealed that perturbation of any point in the pathway thwarted tumor progression, providing strong evidence that targeting the key players could potentially produce a significant benefit for human glioma patients. The elucidation of this signaling pathway is a critical step, since efforts to create a small molecule drug targeting IGFBP2 have so far not been successful, but a number of inhibitors of the other pathway constituents, including ILK, integrin and NF-kB, have been developed.