Multiscale analysis of protein functions and stochastic modelling of gene transcriptional regulatory networks

Genomic and proteomic analyses have attracted a great deal of interests in biological research in recent years. Many methods have been applied to discover useful information contained in the enormous databases of genomic sequences and amino acid sequences. The results of these investigations inspire further research in biological fields in return. These biological sequences, which may be considered as multiscale sequences, have some specific features which need further efforts to characterise using more refined methods. This project aims to study some of these biological challenges with multiscale analysis methods and stochastic modelling approach. The first part of the thesis aims to cluster some unknown proteins, and classify their families as well as their structural classes. A development in proteomic analysis is concerned with the determination of protein functions. The first step in this development is to classify proteins and predict their families. This motives us to study some unknown proteins from specific families, and to cluster them into families and structural classes. We select a large number of proteins from the same families or superfamilies, and link them to simulate some unknown large proteins from these families. We use multifractal analysis and the wavelet method to capture the characteristics of these linked proteins. The simulation results show that the method is valid for the classification of large proteins. The second part of the thesis aims to explore the relationship of proteins based on a layered comparison with their components. Many methods are based on homology of proteins because the resemblance at the protein sequence level normally indicates the similarity of functions and structures. However, some proteins may have similar functions with low sequential identity. We consider protein sequences at detail level to investigate the problem of comparison of proteins. The comparison is based on the empirical mode decomposition (EMD), and protein sequences are detected with the intrinsic mode functions. A measure of similarity is introduced with a new cross-correlation formula. The similarity results show that the EMD is useful for detection of functional relationships of proteins. The third part of the thesis aims to investigate the transcriptional regulatory network of yeast cell cycle via stochastic differential equations. As the investigation of genome-wide gene expressions has become a focus in genomic analysis, researchers have tried to understand the mechanisms of the yeast genome for many years. How cells control gene expressions still needs further investigation. We use a stochastic differential equation to model the expression profile of a target gene. We modify the model with a Gaussian membership function. For each target gene, a transcriptional rate is obtained, and the estimated transcriptional rate is also calculated with the information from five possible transcriptional regulators. Some regulators of these target genes are verified with the related references. With these results, we construct a transcriptional regulatory network for the genes from the yeast Saccharomyces cerevisiae. The construction of transcriptional regulatory network is useful for detecting more mechanisms of the yeast cell cycle.

Influence of preexercise muscle glycogen content on transcriptional activity of metabolic and myogenic genes in well-trained humans

To determine whether pre-exercise muscle glycogen content influences the transcription of several early-response genes involved in the regulation of muscle growth, seven male strength-trained subjects performed one-legged cycling exercise to exhaustion to lower muscle glycogen levels (Low) in one leg compared with the leg with normal muscle glycogen (Norm) and then the following day completed a unilateral bout of resistance training (RT). Muscle biopsies from both legs were taken at rest, immediately after RT, and after 3 h of recovery. Resting glycogen content was higher in the control leg (Norm leg) than in the Low leg (435 ± 87 vs. 193 ± 29 mmol/kg dry wt; P < 0.01). RT decreased glycogen content in both legs (P < 0.05), but postexercise values remained significantly higher in the Norm than the Low leg (312 ± 129 vs. 102 ± 34 mmol/kg dry wt; P < 0.01). GLUT4 (3-fold; P < 0.01) and glycogenin mRNA abundance (2.5-fold; not significant) were elevated at rest in the Norm leg, but such differences were abolished after exercise. Preexercise mRNA abundance of atrogenes was also higher in the Norm compared with the Low leg [atrogin: ?14-fold, P < 0.01; RING (really interesting novel gene) finger: ?3-fold, P < 0.05] but decreased for atrogin in Norm following RT (P < 0.05). There were no differences in the mRNA abundance of myogenic regulatory factors and IGF-I in the Norm compared with the Low leg. Our results demonstrate that 1) low muscle glycogen content has variable effects on the basal transcription of select metabolic and myogenic genes at rest, and 2) any differences in basal transcription are completely abolished after a single bout of heavy resistance training. We conclude that commencing resistance exercise with low muscle glycogen does not enhance the activity of genes implicated in promoting hypertrophy.

Temporal tracking of mineralization and transcriptional developments of shell formation during the early life history of pearl oyster Pinctada maxima

Molluscan larval ontogeny is a highly conserved process comprising three principal developmental stages. A characteristic unique to each of these stages is shell design, termed prodissoconch I, prodissoconch II and dissoconch. These shells vary in morphology, mineralogy and microstructure. The discrete temporal transitions in shell biomineralization between these larval stages are utilized in this study to investigate transcriptional involvement in several distinct biomineralization events. Scanning electron microscopy and X-ray diffraction analysis of P. maxima larvae and juveniles collected throughout post-embryonic ontogenesis, document the mineralogy and microstructure of each shelled stage as well as establishing a timeline for transitions in biomineralization. P. maxima larval samples most representative of these biomineralization distinctions and transitions were analyzed for differential gene expression on the microarray platform PmaxArray 1.0. A number of transcripts are reported as differentially expressed in correlation to the mineralization events of P. maxima larval ontogeny. Some of those isolated are known shell matrix genes while others are novel; these are discussed in relation to potential shell formation roles. This interdisciplinary investigation has linked the shell developments of P. maxima larval ontogeny with corresponding gene expression profiles, furthering the elucidation of shell biomineralization.

Transcriptional regulation of IRS5/DOK4 expression in non-small-cell lung cancer cells

The insulin-receptor substrate family plays important roles in cellular growth, signaling, and survival. Two new members of this family have recently been isolated: IRS5/Dok4 and IRS6/Dok5. This study examines the expression of IRS5/DOK4 in a panel of lung cancer cell lines and tumor specimens. The results demonstrate that expression of IRS5/DOK4 is frequently altered with both elevated and decreased expression in non-small-cell lung cancer (NSCLC) tumor specimens. The altered expression of IRS5/DOK4 observed in tumor samples is not due to aberrant methylation. In vitro cell culture studies demonstrate that treatment of NSCLC cell lines with the histone deacetylase inhibitor trichostatin A (TSA) upregulates IRS5/DOK4. This finding indicates that expression is regulated epigenetically at the level of chromatin remodeling. Chromatin immunoprecipitation experiments confirm that the IRS5/DOK4 promoter has enhanced histone hyperacetylation following treatments with TSA. Finally, hypoxia was demonstrated to downregulate IRS5/DOK4 expression. This expression was restored by TSA. The clinical relevance of altered IRS5/DOK4 expression in NSCLC requires fur ther evaluation.

C. elegans RNA-dependent RNA polymerases rrf-1 and ego-1 silence Drosophila transgenes by differing mechanisms

Drosophila possesses the core gene silencing machinery but, like all insects, lacks the canonical RNA-dependent RNA polymerases (RdRps) that in C. elegans either trigger or enhance two major small RNA-dependent gene silencing pathways. Introduction of two different nematode RdRps into Drosophila showed them to be functional, resulting in differing silencing activities. While RRF-1 enhanced transitive dsRNA-dependent silencing, EGO-1 triggered dsRNA-independent silencing, specifically of transgenes. The strain w; da-Gal4; UAST-ego-1, constitutively expressing ego-1, is capable of silencing transgene including dsRNA hairpin upon a single cross, which created a powerful tool for research in Drosophila. In C. elegans, EGO-1 is involved in transcriptional gene silencing (TGS) of chromosome regions that are unpaired during meiosis. There was no opportunity for meiotic interactions involving EGO-1 in Drosophila that would explain the observed transgene silencing. Transgene DNA is, however, unpaired during the pairing of chromosomes in embryonic mitosis that is an unusual characteristic of Diptera, suggesting that in Drosophila, EGO-1 triggers transcriptional silencing of unpaired DNA during embryonic mitosis. © 2012 Springer Basel.

Transcriptional regulation of flavonoid biosynthesis in nectarine (Prunus persica) by a set of R2R3 MYB transcription factors

Background Flavonoids such as anthocyanins, flavonols and proanthocyanidins, play a central role in fruit colour, flavour and health attributes. In peach and nectarine (Prunus persica) these compounds vary during fruit growth and ripening. Flavonoids are produced by a well studied pathway which is transcriptionally regulated by members of the MYB and bHLH transcription factor families. We have isolated nectarine flavonoid regulating genes and examined their expression patterns, which suggests a critical role in the regulation of flavonoid biosynthesis. Results In nectarine, expression of the genes encoding enzymes of the flavonoid pathway correlated with the concentration of proanthocyanidins, which strongly increases at mid-development. In contrast, the only gene which showed a similar pattern to anthocyanin concentration was UDP-glucose-flavonoid-3-O-glucosyltransferase (UFGT), which was high at the beginning and end of fruit growth, remaining low during the other developmental stages. Expression of flavonol synthase (FLS1) correlated with flavonol levels, both temporally and in a tissue specific manner. The pattern of UFGT gene expression may be explained by the involvement of different transcription factors, which up-regulate flavonoid biosynthesis (MYB10, MYB123, and bHLH3), or repress (MYB111 and MYB16) the transcription of the biosynthetic genes. The expression of a potential proanthocyanidin-regulating transcription factor, MYBPA1, corresponded with proanthocyanidin levels. Functional assays of these transcription factors were used to test the specificity for flavonoid regulation. Conclusions MYB10 positively regulates the promoters of UFGT and dihydroflavonol 4-reductase (DFR) but not leucoanthocyanidin reductase (LAR). In contrast, MYBPA1 trans-activates the promoters of DFR and LAR, but not UFGT. This suggests exclusive roles of anthocyanin regulation by MYB10 and proanthocyanidin regulation by MYBPA1. Further, these transcription factors appeared to be responsive to both developmental and environmental stimuli.

The binding of nuclear factors to the as-1 element in the CaMV 35S promoter is affected by cytosine methylation in vitro

Transcriptional gene silencing (TCS) is often associated with an increased level of cytosine methylation in the affected promoters. The effect of methylation of the cauliflower mosaic virus (CaMV) 35S promoter sequence on its binding to factors present in the nuclei was analyzed by electrophoretic mobility shift assays using extracts of petunia flowers. Specific DNA-protein interactions were detected in the region of the CaMV 35S promoter that contains the as-1 element and the region between -345 and -208. The binding of protein factor(s) to the as-1 element was influenced by cytosine methylation, whereas the binding to the region between -345 and -208 was unaffected. The results suggest that cytosine methylation of the as-1 element potentially affects the activity of the CaMV 35S promoter. © Georg Thieme Verlag KG Stuttgart.

A rapid transcriptional activation is induced by the dormancy-breaking chemical hydrogen cyanamide in kiwifruit (Actinidia deliciosa) buds

Budbreak in kiwifruit (Actinidia deliciosa) can be poor in locations that have warm winters with insufficient winter chilling. Kiwifruit vines are often treated with the dormancy-breaking chemical hydrogen cyanamide (HC) to increase and synchronize budbreak. This treatment also offers a tool to understand the processes involved in budbreak. A genomics approach is presented here to increase our understanding of budbreak in kiwifruit. Most genes identified following HC application appear to be associated with responses to stress, but a number of genes appear to be associated with the reactivation of growth. Three patterns of gene expression were identified: Profile 1, an HC-induced transient activation; Profile 2, an HC-induced transient activation followed by a growth-related activation; and Profile 3, HC- and growth-repressed. One group of genes that was rapidly up-regulated in response to HC was the glutathione S-transferase (GST) class of genes, which have been associated with stress and signalling. Previous budbreak studies, in three other species, also report up-regulated GST expression. Phylogenetic analysis of these GSTs showed that they clustered into two sub-clades, suggesting a strong correlation between their expression and budbreak across species.

Staurosporine augments EGF-mediated EMT in PMC42-LA cells through actin depolymerisation, focal contact size reduction and Snail1 induction : a model for cross-modulation

Background A feature of epithelial to mesenchymal transition (EMT) relevant to tumour dissemination is the reorganization of actin cytoskeleton/focal contacts, influencing cellular ECM adherence and motility. This is coupled with the transcriptional repression of E-cadherin, often mediated by Snail1, Snail2 and Zeb1/δEF1. These genes, overexpressed in breast carcinomas, are known targets of growth factor-initiated pathways, however it is less clear how alterations in ECM attachment cross-modulate to regulate these pathways. EGF induces EMT in the breast cancer cell line PMC42-LA and the kinase inhibitor staurosporine (ST) induces EMT in embryonic neural epithelial cells, with F-actin de-bundling and disruption of cell-cell adhesion, via inhibition of aPKC. Methods PMC42-LA cells were treated for 72 h with 10 ng/ml EGF, 40 nM ST, or both, and assessed for expression of E-cadherin repressor genes (Snail1, Snail2, Zeb1/δEF1) and EMT-related genes by QRT-PCR, multiplex tandem PCR (MT-PCR) and immunofluorescence +/- cycloheximide. Actin and focal contacts (paxillin) were visualized by confocal microscopy. A public database of human breast cancers was assessed for expression of Snail1 and Snail2 in relation to outcome. Results When PMC42-LA were treated with EGF, Snail2 was the principal E-cadherin repressor induced. With ST or ST+EGF this shifted to Snail1, with more extreme EMT and Zeb1/δEF1 induction seen with ST+EGF. ST reduced stress fibres and focal contact size rapidly and independently of gene transcription. Gene expression analysis by MT-PCR indicated that ST repressed many genes which were induced by EGF (EGFR, CAV1, CTGF, CYR61, CD44, S100A4) and induced genes which alter the actin cytoskeleton (NLF1, NLF2, EPHB4). Examination of the public database of breast cancers revealed tumours exhibiting higher Snail1 expression have an increased risk of disease-recurrence. This was not seen for Snail2, and Zeb1/δEF1 showed a reverse correlation with lower expression values being predictive of increased risk. Conclusion ST in combination with EGF directed a greater EMT via actin depolymerisation and focal contact size reduction, resulting in a loosening of cell-ECM attachment along with Snail1-Zeb1/δEF1 induction. This appeared fundamentally different to the EGF-induced EMT, highlighting the multiple pathways which can regulate EMT. Our findings add support for a functional role for Snail1 in invasive breast cancer.

Tyrosine phosphorylation mediates ConA-induced membrane type 1-matrix metalloproteinase expression and matrix metalloproteinase-2 activation in MDA-MB-231 human breast carcinoma cells

ConA-induced cell surface activation of pro-matrix metalloproteinase-2 (pro-MMP-2) by MDA-MB-231 human breast cancer cells is apparently mediated by up-regulation of membrane type 1 MMP (MT1-MMP) through transcriptional and posttranscriptional mechanisms. Here, we have explored the respective roles of cell surface clustering and protein tyrosine phosphorylation in the ConA- induction effects. Treatment with succinyl-ConA, a variant lacking significant clusterability, partially stimulated MT1-MMP mRNA and protein levels but did not induce MMP-2 activation, suggesting that clustering contributes to the transcriptional regulation by ConA but appears to be critical for the nontranscriptional component. We further found that genistein, an inhibitor of tyrosine phosphorylation, blocked ConA-induced pro-MMP-2 activation and ConA-induced MT1-MMP mRNA level in a dose-dependent manner, implicating tyrosine phosphorylation in the transcriptional aspect. This was confirmed by the dose-dependent promotion of pro-MMP-2 activation by sodium orthovanadate in the presence of suboptimal concentrations of ConA (7.5 μg/ml), with optimal effects seen at 25 μg/g orthovanadate. Genistein did not inhibit the ConA potentiation of MMP-2 activation in MCF-7 cells, in which transfected MT1-MMP is driven by a heterologous promoter, supporting the major implication of phosphotyrosine in the transcriptional component of ConA regulation. These data describe a major signaling event upstream of MT1- MMP induction by ConA and set the stage for further analysis of the nontranscriptional component.

Type I collagen abrogates the clathrin-mediated internalization of membrane type 1 matrix metalloproteinase (MT1-MMP) via the MT1-MMP hemopexin domain

Type I collagen (Col I)-stimulated matrix metalloproteinase-2 (MMP-2) activation via membrane type 1 MMP (MT1-MMP) involves both a transcriptional increase in MT1-MMP expression and a nontranscriptional response mediated by preexisting MT1-MMP. In order to identify which MT1-MMP domains were required for the nontranscriptional response, MCF-7 cells that lack endogenous MT1-MMP were transfected with either wild type or domain mutant MT1-MMP constructs. We observed that mutant constructs lacking the MT1-MMP cytoplasmic tail were able to activate MMP-2 in response to Col I but not a construct lacking the MT1-MMP hemopexin domain. Col I did not alter total MT1-MMP protein levels; nor did it appear to directly induce MT1-MMP oligomerization. Col I did, however, redistribute preexisting MT1-MMP to the cell periphery compared with unstimulated cells that displayed amore diffuse staining pattern. In addition, Col I blocked the internalization of MT1-MMP in a dynamin-dependent manner via clathrin-coated pit-mediated endocytosis. This mechanism of impaired internalization is different from that reported for concanavalin A, since it is not mediated by the cytoplasmic tail of MT1-MMP but rather by the hemopexin domain. In summary, upon Col I binding to its cell surface receptor, MT1-MMP internalization via clathrin-coated pit-mediated endocytosis is impaired through interactions with the hemopexin domain, thereby regulating its function and ability to activate MMP-2.

Functional role for senataxin, defective in ataxia oculomotor apraxia type 2, in transcriptional regulation

Ataxia oculomotor apraxia type 2 (AOA2) is an autosomal recessive neurodegenerative disorder characterized by cerebellar ataxia and oculomotor apraxia. The gene mutated in AOA2, SETX, encodes senataxin, a putative DNA/RNA helicase which shares high homology to the yeast Sen1p protein and has been shown to play a role in the response to oxidative stress. To investigate further the function of senataxin, we identified novel senataxin-interacting proteins, the majority of which are involved in transcription and RNA processing, including RNA polymerase II. Binding of RNA polymerase II to candidate genes was significantly reduced in senataxin deficient cells and this was accompanied by decreased transcription of these genes, suggesting a role for senataxin in the regulation/modulation of transcription. RNA polymerase II-dependent transcription termination was defective in cells depleted of senataxin in keeping with the observed interaction of senataxin with poly(A) binding proteins 1 and 2. Splicing efficiency of specific mRNAs and alternate splice-site selection of both endogenous genes and artificial minigenes were altered in senataxin depleted cells. These data suggest that senataxin, similar to its yeast homolog Sen1p, plays a role in coordinating transcriptional events, in addition to its role in DNA repair.

Transcriptional pathway signatures predict MEK addiction and response to selumetinib (AZD6244)

Selumetinib (AZD6244, ARRY-142886) is a selective, non-ATP-competitive inhibitor of mitogen-activated protein/extracellular signal-regulated kinase kinase (MEK)-1/2. The range of antitumor activity seen preclinically and in patients highlights the importance of identifying determinants of response to this drug. In large tumor cell panels of diverse lineage, we show that MEK inhibitor response does not have an absolute correlation with mutational or phospho-protein markers of BRAF/MEK, RAS, or phosphoinositide 3-kinase (PI3K) activity. We aimed to enhance predictivity by measuring pathway output through coregulated gene networks displaying differential mRNA expression exclusive to resistant cell subsets and correlated to mutational or dynamic pathway activity. We discovered an 18-gene signature enabling measurement of MEK functional output independent of tumor genotype. Where the MEK pathway is activated but the cells remain resistant to selumetinib, we identified a 13-gene signature that implicates the existence of compensatory signaling from RAS effectors other than PI3K. The ability of these signatures to stratify samples according to functional activation of MEK and/or selumetinib sensitivity was shown in multiple independent melanoma, colon, breast, and lung tumor cell lines and in xenograft models. Furthermore, we were able to measure these signatures in fixed archival melanoma tumor samples using a single RT-qPCR-based test and found intergene correlations and associations with genetic markers of pathway activity to be preserved. These signatures offer useful tools for the study of MEK biology and clinical application of MEK inhibitors, and the novel approaches taken may benefit other targeted therapies.

Innate transcriptional networks activated in bladder in response to uropathogenic Escherichia coli drive diverse biological pathways and rapid synthesis of IL-10 for defense against bacterial urinary tract infection

Early transcriptional activation events that occur in bladder immediately following bacterial urinary tract infection (UTI) are not well defined. In this study, we describe the whole bladder transcriptome of uropathogenic Escherichia coli (UPEC) cystitis in mice using genome-wide expression profiling to define the transcriptome of innate immune activation stemming from UPEC colonization of the bladder. Bladder RNA from female C57BL/6 mice, analyzed using 1.0 ST-Affymetrix microarrays, revealed extensive activation of diverse sets of innate immune response genes, including those that encode multiple IL-family members, receptors, metabolic regulators, MAPK activators, and lymphocyte signaling molecules. These were among 1564 genes differentially regulated at 2 h postinfection, highlighting a rapid and broad innate immune response to bladder colonization. Integrative systems-level analyses using InnateDB (http://www.innatedb.com) bioinformatics and ingenuity pathway analysis identified multiple distinct biological pathways in the bladder transcriptome with extensive involvement of lymphocyte signaling, cell cycle alterations, cytoskeletal, and metabolic changes. A key regulator of IL activity identified in the transcriptome was IL-10, which was analyzed functionally to reveal marked exacerbation of cystitis in IL-10–deficient mice. Studies of clinical UTI revealed significantly elevated urinary IL-10 in patients with UPEC cystitis, indicating a role for IL-10 in the innate response to human UTI. The whole bladder transcriptome presented in this work provides new insight into the diversity of innate factors that determine UTI on a genome-wide scale and will be valuable for further data mining. Identification of protective roles for other elements in the transcriptome will provide critical new insight into the complex cascade of events that underpin UTI.