Pin1 Catalyzes Conformational Changes of Thr-187 in p27Kip1 and Mediates Its Stability through a Polyubiquitination Process

The cis-trans peptidylprolyl isomerase Pin1 plays a critical role in regulating a subset of phosphoproteins by catalyzing conformational changes on the phosphorylated Ser/Thr-Pro motifs. The phosphorylation-directed ubiquitination is one of the major mechanisms to regulate the abundance of p27^Kip1. In this study, we demonstrate that Pin1 catalyzes the cis-trans conformational changes of p27^Kip1 and further mediates its stability through the polyubiquitination mechanism. Our results show that the phosphorylated Thr-187-Pro motif in p27^Kip1 is a key Pin1-binding site. In addition, NMR analyses show that this phosphorylated Thr-187-Pro site undergoes conformational change catalyzed by Pin1. Moreover, in Pin1 knock-out mouse embryonic fibroblasts, p27^Kip1 has a shorter lifetime and displays a higher degree of polyubiquitination than in Pin1 wild-type mouse embryonic fibroblasts, suggesting that Pin1 plays a critical role in regulating p27^Kip1 degradation. Additionally, Pin1 dramatically reduces the interaction between p27^Kip1 and Cks1, possibly via isomerizing the cis-trans conformation of p27^Kip1. Our study thus reveals a novel regulatory mechanism for p27^Kip1 stability and sheds new light on the biological function of Pin1 as a general regulator of protein stability.Cellular differentiation and cell cycle inhibition are tightly controlled via sensitive molecular mechanisms. p27^Kip1, a member of the Cip/Kip family, is an essential cell cycle inhibitor that functions largely during the G₀/G₁ phase where it promotes the assembly of the cyclin D1-CDK4 complex and inhibits the kinase activity of the cyclin E-CDK2 complex in the G₁-S phase (1–4). Several review articles have elegantly summarized and discussed the detailed cellular functions of p27^Kip1 (1–6). p27^Kip1 is also a phosphoprotein with multiple Ser/Thr phosphorylation sites, including Ser-10, Ser-178, and Thr-187, followed by a proline residue. Hence, these motifs are potential substrate sites for proline-directed kinases (5, 6). Compared with Ser-178, which has not yet been well studied, the phosphorylation of Ser-10 and Thr-187 has been well characterized to be important for the regulation of p27^Kip1 function. For instance, Ser-10 has been found to be the major phosphorylation site of p27^Kip1 (7) and to play an important role in regulating cell migration (8–10), although the regulation of Ser-10 phosphorylation is still not completely defined (11, 12).In contrast to Ser-10 and Thr-178, Thr-187 is the best characterized phosphorylation site on p27^Kip1 and is known to regulate the complex formation of p27^Kip1-cyclin E-CDK2 (1–2). In addition, it is also widely accepted that Thr-187 plays a crucial role in determining the abundance of mature p27^Kip1 proteins. The phosphorylation of Thr-187 directs p27^Kip1 to an SCF^Skp2 ubiquitin ligase complex (consisting of Skp2-Skp1-Cks1-Cul1-Roc1), which in turn promotes the polyubiquitination and degradation of p27^Kip1 (13, 14). The crystal structure of the Skp1-Skp2-Cks1-p27^Kip1 phosphopeptide complex shows that p27^Kip1 binds both Cks1 and Skp2 and that the C terminus of Skp2 and Cks1 forms the substrate recognition core of the SCF complex (15). Furthermore, the structure of this complex has revealed that the phosphorylation of Thr-187 in p27^Kip1 is recognized by the phosphate-binding site of Cks1, indicating that Cks1 is not only a facilitator but also an indispensable component in p27^Kip1 degradation machinery (15).Pin1 is a unique peptidyl-prolyl isomerase (PPIase)² that recognizes only the phosphorylated Ser/Thr motif preceding a proline residue (16). In addition, Pin1 is very prominent in isomerizing the cis-trans conformation of prolyl-peptidyl bonds in its substrates, resulting in either the modification of their function (e.g. c-Jun (17), β-catenin (18), Bax (19), and Notch1 (20)) or modulation of their stability (e.g. cyclin D1 (21), p53 (22, 23), and NF-κB (24)). Loss of Pin1 in mice results in several phenotypes similar to those of cyclin D1-null mice (21) and neuronal degenerative phenotypes (25–28), suggesting the conformational changes mediated by Pin1 may be crucial for the normal functioning of cells. Additionally, Pin1 also plays important roles in cancer and other cellular events, which have been extensively discussed in several recent review articles (29–33).In this study, we show that Pin1 binds to p27^Kip1, mainly through the phosphorylated Thr-187-Pro motif, and causes subsequent prolyl isomerization of this cell cycle protein. Moreover, we also find that Pin1 can protect p27^Kip1 from degradation. Importantly, we demonstrate that by catalyzing conformational changes in p27^Kip1, Pin1 hinders its association with Cks1, resulting in a reduction of polyubiquitination of p27^Kip1 and protecting its degradation by SCF^Skp2 complexes. Our results suggest that the cis-trans isomerization catalyzed by Pin1 represents a novel regulatory mechanism during post-phosphorylation of proteins and polyubiquitination-directed degradation pathways.     

Printor, a Novel TorsinA-interacting Protein Implicated in Dystonia Pathogenesis

Early onset generalized dystonia (DYT1) is an autosomal dominant neurological disorder caused by deletion of a single glutamate residue (torsinA ΔE) in the C-terminal region of the AAA⁺ (ATPases associated with a variety of cellular activities) protein torsinA. The pathogenic mechanism by which torsinA ΔE mutation leads to dystonia remains unknown. Here we report the identification and characterization of a 628-amino acid novel protein, printor, that interacts with torsinA. Printor co-distributes with torsinA in multiple brain regions and co-localizes with torsinA in the endoplasmic reticulum. Interestingly, printor selectively binds to the ATP-free form but not to the ATP-bound form of torsinA, supporting a role for printor as a cofactor rather than a substrate of torsinA. The interaction of printor with torsinA is completely abolished by the dystonia-associated torsinA ΔE mutation. Our findings suggest that printor is a new component of the DYT1 pathogenic pathway and provide a potential molecular target for therapeutic intervention in dystonia.Early onset generalized torsion dystonia (DYT1) is the most common and severe form of hereditary dystonia, a movement disorder characterized by involuntary movements and sustained muscle spasms (1). This autosomal dominant disease has childhood onset and its dystonic symptoms are thought to result from neuronal dysfunction rather than neurodegeneration (2, 3). Most DYT1 cases are caused by deletion of a single glutamate residue at positions 302 or 303 (torsinA ΔE) of the 332-amino acid protein torsinA (4). In addition, a different torsinA mutation that deletes amino acids Phe³²³–Tyr³²⁸ (torsinA Δ323–328) was identified in a single family with dystonia (5), although the pathogenic significance of this torsinA mutation is unclear because these patients contain a concomitant mutation in another dystonia-related protein, ϵ-sarcoglycan (6). Recently, genetic association studies have implicated polymorphisms in the torsinA gene as a genetic risk factor in the development of adult-onset idiopathic dystonia (7, 8).TorsinA contains an N-terminal endoplasmic reticulum (ER)³ signal sequence and a 20-amino acid hydrophobic region followed by a conserved AAA⁺ (ATPases associated with a variety of cellular activities) domain (9, 10). Because members of the AAA⁺ family are known to facilitate conformational changes in target proteins (11, 12), it has been proposed that torsinA may function as a molecular chaperone (13, 14). TorsinA is widely expressed in brain and multiple other tissues (15) and is primarily associated with the ER and nuclear envelope (NE) compartments in cells (16–20). TorsinA is believed to mainly reside in the lumen of the ER and NE (17–19) and has been shown to bind lamina-associated polypeptide 1 (LAP1) (21), lumenal domain-like LAP1 (LULL1) (21), and nesprins (22). In addition, recent evidence indicates that a significant pool of torsinA exhibits a topology in which the AAA⁺ domain faces the cytoplasm (20). In support of this topology, torsinA is found in the cytoplasm, neuronal processes, and synaptic terminals (2, 3, 15, 23–26) and has been shown to bind cytosolic proteins snapin (27) and kinesin light chain 1 (20). TorsinA has been proposed to play a role in several cellular processes, including dopaminergic neurotransmission (28–31), NE organization and dynamics (17, 22, 32), and protein trafficking (27, 33). However, the precise biological function of torsinA and its regulation remain unknown.To gain insights into torsinA function, we performed yeast two-hybrid screens to search for torsinA-interacting proteins in the brain. We report here the isolation and characterization of a novel protein named printor (protein interactor of torsinA) that interacts selectively with wild-type (WT) torsinA but not the dystonia-associated torsinA ΔE mutant. Our data suggest that printor may serve as a cofactor of torsinA and provide a new molecular target for understanding and treating dystonia.     

Analysis of Free Energy Signals Arising from Nucleotide Hybridization Between rRNA and mRNA Sequences during Translation in Eubacteria

A decoding algorithm is tested that mechanistically models the progressive alignments that arise as the mRNA moves past the rRNA tail during translation elongation. Each of these alignments provides an opportunity for hybridization between the single-stranded, -terminal nucleotides of the 16S rRNA and the spatially accessible window of mRNA sequence, from which a free energy value can be calculated. Using this algorithm we show that a periodic, energetic pattern of frequency 1/3 is revealed. This periodic signal exists in the majority of coding regions of eubacterial genes, but not in the non-coding regions encoding the 16S and 23S rRNAs. Signal analysis reveals that the population of coding regions of each bacterial species has a mean phase that is correlated in a statistically significant way with species () content. These results suggest that the periodic signal could function as a synchronization signal for the maintenance of reading frame and that codon usage provides a mechanism for manipulation of signal phase.[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32]     

Phosphorylation Dependence of Hsp27 Multimeric Size and Molecular Chaperone Function

The molecular chaperone Hsp27 exists as a distribution of large oligomers that are disassembled by phosphorylation at Ser-15, -78, and -82. It is controversial whether the unphosphorylated Hsp27 or the widely used triple Ser-to-Asp phospho-mimic mutant is the more active molecular chaperone in vitro. This question was investigated here by correlating chaperone activity, as measured by the aggregation of reduced insulin or α-lactalbumin, with Hsp27 self-association as monitored by analytical ultracentrifugation. Furthermore, because the phospho-mimic is generally assumed to reproduce the phosphorylated molecule, the size and chaperone activity of phosphorylated Hsp27 were compared with that of the phospho-mimic. Hsp27 was triply phosphorylated by MAPKAP-2 kinase, and phosphorylation was tracked by urea-PAGE. An increasing degree of suppression of insulin or α-lactalbumin aggregation correlated with a decreasing Hsp27 self-association, which was the least for phosphorylated Hsp27 followed by the mimic followed by the unphosphorylated protein. It was also found that Hsp27 added to pre-aggregated insulin did not reverse aggregation but did inhibit these aggregates from assembling into even larger aggregates. This chaperone activity appears to be independent of Hsp27 phosphorylation. In conclusion, the most active chaperone of insulin and α-lactalbumin was the Hsp27 (elongated) dimer, the smallest Hsp27 subunit observed under physiological conditions. Next, the Hsp27 phospho-mimic is only a partial mimic of phosphorylated Hsp27, both in self-association and in chaperone function. Finally, the efficient inhibition of insulin aggregation by Hsp27 dimer led to the proposal of two models for this chaperone activity.Oligomeric heat shock protein 27 (Hsp27)² is a ubiquitous mammalian protein with a variety of functions in health and disease (1–8). These functions include ATP-independent chaperone activity in response to environmental stress, e.g. heat shock and oxidative stress, control of apoptosis, and regulation of actin cytoskeleton dynamics. Hsp27 is a member of the α-crystallin small heat shock protein family of which αB-crystallin is the archetype. These proteins are characterized by an α-crystallin domain of 80–90 residues consisting of roughly eight β-strands that form an intermolecular β-sheet interaction interface within a dimer, the basic building subunit of the oligomer (2, 4, 9–11).Hsp27 is in equilibrium between high molecular weight oligomers and much lower molecular weight multimers. It has been reported that unphosphorylated Hsp27 includes predominantly a distribution of high molecular species ranging in size from 12-mer to 35-mer (12–19). Phosphorylation of Hsp27 at serines 15, 78, and 82 by the p38-activated MAPKAP-2 kinase (20–22) or the use of the triple Ser-to-Asp phospho-mimic results in a major shift in the equilibrium toward much smaller multimers (23) and in an alteration of its function (1, 3, 6, 7, 24, 25). The size distribution of the smaller species has been reported to be between monomer and tetramer (12–16, 18, 19).Small heat shock proteins, including Hsp27, behave as ATP-independent molecular chaperones during cellular heat shock. They bind partially unfolded proteins and prevent their aggregation until the proteins can be refolded by larger ATP-dependent chaperones or are digested (7, 8, 26). This function includes the up-regulation and/or phosphorylation of Hsp27.It is not entirely clear what the role of Hsp27 size and phosphorylation state plays in its heat shock function because there are conflicting results in the literature. Some in vitro studies concluded that the unphosphorylated oligomeric Hsp27 (or the murine isoform Hsp25) protects proteins against aggregation better than does the phosphorylation mimic (13, 19, 27), whereas others found no difference (16, 28, 29), and still other studies found that the mimic protects better than does the unphosphorylated wild type (27, 30, 31). In-cell studies found that phosphorylation of Hsp27 was essential for thermo-protection of actin filaments (32), and the Hsp27 phosphorylation mimic decreased inclusion body formation better than did unphosphorylated Hsp27 (33). This study was undertaken to investigate the molecular chaperone function of Hsp27 by correlating chaperone activity with Hsp27 size and by comparing fully phosphorylated Hsp27 with its phospho-mimic.     

Algorithms for Finding Small Attractors in Boolean Networks

A Boolean network is a model used to study the interactions between different genes in genetic regulatory networks. In this paper, we present several algorithms using gene ordering and feedback vertex sets to identify singleton attractors and small attractors in Boolean networks. We analyze the average case time complexities of some of the proposed algorithms. For instance, it is shown that the outdegree-based ordering algorithm for finding singleton attractors works in time for , which is much faster than the naive time algorithm, where is the number of genes and is the maximum indegree. We performed extensive computational experiments on these algorithms, which resulted in good agreement with theoretical results. In contrast, we give a simple and complete proof for showing that finding an attractor with the shortest period is NP-hard.[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32]     

A Ca2+ Release-activated Ca2+ (CRAC) Modulatory Domain (CMD) within STIM1 Mediates Fast Ca2+-dependent Inactivation of ORAI1 Channels

STIM1 and ORAI1, the two limiting components in the Ca²⁺ release-activated Ca²⁺ (CRAC) signaling cascade, have been reported to interact upon store depletion, culminating in CRAC current activation. We have recently identified a modulatory domain between amino acids 474 and 485 in the cytosolic part of STIM1 that comprises 7 negatively charged residues. A STIM1 C-terminal fragment lacking this domain exhibits enhanced interaction with ORAI1 and 2–3-fold higher ORAI1/CRAC current densities. Here we focused on the role of this CRAC modulatory domain (CMD) in the fast inactivation of ORAI1/CRAC channels, utilizing the whole-cell patch clamp technique. STIM1 mutants either with C-terminal deletions including CMD or with 7 alanines replacing the negative amino acids within CMD gave rise to ORAI1 currents that displayed significantly reduced or even abolished inactivation when compared with STIM1 mutants with preserved CMD. Consistent results were obtained with cytosolic C-terminal fragments of STIM1, both in ORAI1-expressing HEK 293 cells and in RBL-2H3 mast cells containing endogenous CRAC channels. Inactivation of the latter, however, was much more pronounced than that of ORAI1. The extent of inactivation of ORAI3 channels, which is also considerably more prominent than that of ORAI1, was also substantially reduced by co-expression of STIM1 constructs missing CMD. Regarding the dependence of inactivation on Ca²⁺, a decrease in intracellular Ca²⁺ chelator concentrations promoted ORAI1 current fast inactivation, whereas Ba²⁺ substitution for extracellular Ca²⁺ completely abrogated it. In summary, CMD within the STIM1 cytosolic part provides a negative feedback signal to Ca²⁺ entry by triggering fast Ca²⁺-dependent inactivation of ORAI/CRAC channels.The Ca²⁺ release-activated Ca²⁺ (CRAC)⁵ channel is one of the best characterized store-operated entry pathways (1–7). Substantial efforts have led to identification of two key components of the CRAC channel machinery: the stromal interaction molecule 1 (STIM1), which is located in the endoplasmic reticulum and acts as a Ca²⁺ sensor (8–10), and ORAI1/CRACM1, the pore-forming subunit of the CRAC channel (11–13). Besides ORAI1, two further homologues named ORAI2 and ORAI3 belong to the ORAI channel family (12, 14).STIM1 senses endoplasmic reticulum store depletion primarily by its luminal EF-hand in its N terminus (8, 15), redistributes close to the plasma membrane, where it forms puncta-like structures, and co-clusters with ORAI1, leading to inward Ca²⁺ currents (12, 16–19). The STIM1 C terminus, located in the cytosol, contains two coiled-coil regions overlapping with an ezrin-radixin-moesin (ERM)-like domain followed by a serine/proline- and a lysine-rich region (2, 8, 20–22). Three recent studies have described the essential ORAI-activating region within the ERM domain, termed SOAR (Stim ORAI-activating region) (23), OASF (ORAI-activating small fragment) (24), and CAD (CRAC-activating domain) (25), including the second coiled coil domain and the following ∼55 amino acids. We and others have provided evidence that store depletion leads to a dynamic coupling of STIM1 to ORAI1 (26–28) that is mediated by a direct interaction of the STIM1 C terminus with ORAI1 C terminus probably involving the putative coiled-coil domain in the latter (27).Furthermore, different groups have proven that the C terminus of STIM1 is sufficient to activate CRAC as well as ORAI1 channels independent of store depletion (22–25, 27, 29). We have identified that OASF-(233–474) or shorter fragments exhibit further enhanced coupling to ORAI1 resulting in 3-fold increased constitutive Ca²⁺ currents. A STIM1 fragment containing an additional cluster of anionic amino acids C-terminal to position 474 displays weaker interaction with ORAI1 as well as reduced Ca²⁺ current comparable with that mediated by wild-type STIM1 C terminus. Hence, we have suggested that these 11 amino acids (474–485) act in a modulatory manner onto ORAI1; however, their detailed mechanistic impact within the STIM1/ORAI1 signaling machinery has remained so far unclear.In this study, we focused on the impact of this negative cluster on fast inactivation of STIM1-mediated ORAI Ca²⁺ currents. Lis et al. (30) have shown that all three ORAI homologues display distinct inactivation profiles, where ORAI2 and ORAI3 show a much more pronounced fast inactivation than ORAI1. Moreover, it has been reported (31) that different expression levels of STIM1 to ORAI1 affect the properties of CRAC current inactivation. Yamashita et al. (32) have demonstrated a linkage between the selectivity filter of ORAI1 and its Ca²⁺-dependent fast inactivation. Here we provide evidence that a cluster of acidic residues within the C terminus of STIM1 is involved in the fast inactivation of ORAI1 and further promotes that of ORAI3 and native CRAC currents.     

Regulation of Microtubule Dynamic Instability in Vitro by Differentially Phosphorylated Stathmin

Stathmin is an important regulator of microtubule polymerization and dynamics. When unphosphorylated it destabilizes microtubules in two ways, by reducing the microtubule polymer mass through sequestration of soluble tubulin into an assembly-incompetent T2S complex (two α:β tubulin dimers per molecule of stathmin), and by increasing the switching frequency (catastrophe frequency) from growth to shortening at plus and minus ends by binding directly to the microtubules. Phosphorylation of stathmin on one or more of its four serine residues (Ser¹⁶, Ser²⁵, Ser³⁸, and Ser⁶³) reduces its microtubule-destabilizing activity. However, the effects of phosphorylation of the individual serine residues of stathmin on microtubule dynamic instability have not been investigated systematically. Here we analyzed the effects of stathmin singly phosphorylated at Ser¹⁶ or Ser⁶³, and doubly phosphorylated at Ser²⁵ and Ser³⁸, on its ability to modulate microtubule dynamic instability at steady-state in vitro. Phosphorylation at either Ser¹⁶ or Ser⁶³ strongly reduced or abolished the ability of stathmin to bind to and sequester soluble tubulin and its ability to act as a catastrophe factor by directly binding to the microtubules. In contrast, double phosphorylation of Ser²⁵ and Ser³⁸ did not affect the binding of stathmin to tubulin or microtubules or its catastrophe-promoting activity. Our results indicate that the effects of stathmin on dynamic instability are strongly but differently attenuated by phosphorylation at Ser¹⁶ and Ser⁶³ and support the hypothesis that selective targeting by Ser¹⁶-specific or Ser⁶³-specific kinases provides complimentary mechanisms for regulating microtubule function.Stathmin is an 18-kDa ubiquitously expressed microtubule-destabilizing phosphoprotein whose activity is modulated by phosphorylation of its four serine residues, Ser¹⁶, Ser²⁵, Ser³⁸, and Ser⁶³ (1–7). Several classes of kinases have been identified that phosphorylate stathmin, including kinases associated with cell growth and differentiation such as members of the mitogen-activated protein kinase (MAPK)2 family, cAMP-dependent protein kinase (1–5, 8–11), and kinases associated with cell cycle regulation such as cyclin-dependent kinase 1 (3, 12–14). Phosphorylation of stathmin is required for cell cycle progression through mitosis and for proper assembly/function of the mitotic spindle (3, 13–16). Inhibition of stathmin phosphorylation produces strong mitotic phenotypes characterized by disassembly and disorganization of mitotic spindles and abnormal chromosome distributions (3, 13–14).Stathmin is known to destabilize microtubules in two ways. One is by binding to soluble tubulin and forming a stable complex that cannot polymerize into microtubules, consisting of one molecule of stathmin and two molecules of tubulin (T2S complex) (17–24). Addition of stathmin to microtubules in equilibrium with soluble tubulin results in sequestration of the tubulin and a reduction in the level of microtubule polymer (17–18, 22, 25–28). In addition to reducing the amount of assembled polymer, tubulin sequestration by stathmin has been shown to increase the switching frequency at microtubule plus ends from growth to shortening (called the catastrophe frequency) as the microtubules relax to a new steady state (17, 29). The second way is by binding directly to microtubules (27–30). The direct binding of stathmin to microtubules increases the catastrophe frequency at both ends of the microtubules and considerably more strongly at minus ends than at plus ends (27). Consistent with its strong catastrophe-promoting activity at minus ends, stathmin increases the treadmilling rate of steady-state microtubules in vitro (27). These results have led to the suggestion that stathmin might be an important cellular regulator of minus-end microtubule dynamics (27).Phosphorylation of stathmin diminishes its ability to regulate microtubule polymerization (3, 14, 25–26). Phosphorylation of Ser¹⁶ or Ser⁶³ appears to be more critical than phosphorylation of Ser²⁵ and Ser³⁸ for the ability of stathmin to bind to soluble tubulin and to inhibit microtubule assembly in vitro (3, 25). Inhibition of stathmin phosphorylation induces defects in spindle assembly and organization (3, 14) suggesting that not only soluble tubulin-microtubule levels are regulated by phosphorylation of stathmin, but the dynamics of microtubules could also be regulated in a phosphorylation-dependent manner.It is not known how phosphorylation at any of the four serine residues of stathmin affects its ability to regulate microtubule dynamics and, specifically, its ability to increase the catastrophe frequency at plus and minus ends due to its direct interaction with microtubules. Thus, we determined the effects of stathmin individually phosphorylated at either Ser¹⁶ or Ser⁶³ and doubly phosphorylated at both Ser²⁵ and Ser³⁸ on dynamic instability at plus and minus ends in vitro at microtubule polymer steady state and physiological pH (pH 7.2). We find that phosphorylation of Ser¹⁶ strongly reduces the direct catastrophe-promoting activity of stathmin at plus ends and abolishes it at minus ends, whereas phosphorylation of Ser⁶³ abolishes the activity at both ends. The effects of phosphorylation of individual serines correlated well with stathmin''s reduced abilities to form stable T2S complexes, to inhibit microtubule polymerization, and to bind to microtubules. In contrast, double phosphorylation of Ser²⁵ and Ser³⁸ did not alter the ability of stathmin to modulate dynamic instability at the microtubule ends, its ability to form a stable T2S complex, or its ability to bind to microtubules. The data further support the hypotheses that phosphorylation of stathmin on either Ser¹⁶ or Ser⁶³ plays a critical role in regulating microtubule polymerization and dynamics in cells.     

The Wavelet-Based Cluster Analysis for Temporal Gene Expression Data

A variety of high-throughput methods have made it possible to generate detailed temporal expression data for a single gene or large numbers of genes. Common methods for analysis of these large data sets can be problematic. One challenge is the comparison of temporal expression data obtained from different growth conditions where the patterns of expression may be shifted in time. We propose the use of wavelet analysis to transform the data obtained under different growth conditions to permit comparison of expression patterns from experiments that have time shifts or delays. We demonstrate this approach using detailed temporal data for a single bacterial gene obtained under 72 different growth conditions. This general strategy can be applied in the analysis of data sets of thousands of genes under different conditions.[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29]