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.     

Splitting the BLOSUM Score into Numbers of Biological Significance

Mathematical tools developed in the context of Shannon information theory were used to analyze the meaning of the BLOSUM score, which was split into three components termed as the BLOSUM spectrum (or BLOSpectrum). These relate respectively to the sequence convergence (the stochastic similarity of the two protein sequences), to the background frequency divergence (typicality of the amino acid probability distribution in each sequence), and to the target frequency divergence (compliance of the amino acid variations between the two sequences to the protein model implicit in the BLOCKS database). This treatment sharpens the protein sequence comparison, providing a rationale for the biological significance of the obtained score, and helps to identify weakly related sequences. Moreover, the BLOSpectrum can guide the choice of the most appropriate scoring matrix, tailoring it to the evolutionary divergence associated with the two sequences, or indicate if a compositionally adjusted matrix could perform better.[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]     

Metal Deficiency Increases Aberrant Hydrophobicity of Mutant Superoxide Dismutases That Cause Amyotrophic Lateral Sclerosis

The mechanisms by which mutant variants of Cu/Zn-superoxide dismutase (SOD1) cause familial amyotrophic lateral sclerosis are not clearly understood. Evidence to date suggests that altered conformations of amyotrophic lateral sclerosis mutant SOD1s trigger perturbations of cellular homeostasis that ultimately cause motor neuron degeneration. In this study we correlated the metal contents and disulfide bond status of purified wild-type (WT) and mutant SOD1 proteins to changes in electrophoretic mobility and surface hydrophobicity as detected by 1-anilinonaphthalene-8-sulfonic acid (ANS) fluorescence. As-isolated WT and mutant SOD1s were copper-deficient and exhibited mobilities that correlated with their expected negative charge. However, upon disulfide reduction and demetallation at physiological pH, both WT and mutant SOD1s underwent a conformational change that produced a slower mobility indicative of partial unfolding. Furthermore, although ANS did not bind appreciably to the WT holoenzyme, incubation of metal-deficient WT or mutant SOD1s with ANS increased the ANS fluorescence and shifted its peak toward shorter wavelengths. This increased interaction with ANS was greater for the mutant SOD1s and could be reversed by the addition of metal ions, especially Cu²⁺, even for SOD1 variants incapable of forming the disulfide bond. Overall, our findings support the notion that misfolding associated with metal deficiency may facilitate aberrant interactions of SOD1 with itself or with other cellular constituents and may thereby contribute to neuronal toxicity.The sequence of events by which more than 100 mutations in the gene encoding Cu/Zn-superoxide dismutase (SOD1)³ cause familial forms of amyotrophic lateral sclerosis (ALS) is unknown. Studies of purified SOD1 proteins and cellular or rodent models of SOD1-linked ALS suggest that impaired metal ion binding or misfolding of mutant SOD1 proteins in the cellular environment may be related to their toxicity (1–10). Available evidence suggests that partially unfolded mutant SOD1 species could contribute to motor neuron death by promoting abnormal interactions that produce cellular dysfunction (11–16).In previous studies we characterized physicochemical properties of 14 different biologically metallated ALS SOD1 mutants (17) and demonstrated altered thermal stabilities of these mutants compared with wild-type (WT) SOD1 (18). These “as-isolated” SOD1 proteins, which contain variable amounts of copper and zinc, were broadly grouped into two classes based on their ability to incorporate and retain metal ions with high affinity. WT-like SOD1 mutants retain the ability to bind copper and zinc ions and exhibit dismutase activity similar to the normal enzyme, whereas metal binding region (MBR) mutants are significantly deficient in copper and/or zinc (17, 19). We also observed that ALS-associated SOD1 mutants were more susceptible than the WT enzyme to reduction of the intrasubunit disulfide bond between Cys-57 and Cys-146 (20). The significance of these results is that even WT-like mutants, which exhibit a nearly normal backbone structure (21–23), may be vulnerable to destabilizing influences in vivo. Our group and others subsequently showed that the mutant SOD1 proteins share a susceptibility to increased hydrophobicity under conditions that reduce disulfide bonds and/or chelate metal ions (5) and that similar hydrophobic species exist in tissue lysates from mutant SOD1 transgenic mice (4–6). One consequence of such hydrophobic exposure could be the facilitation of abnormal interactions between the mutant enzymes and other cellular constituents (e.g. chaperones, mitochondrial components, or other targets), which might influence pathways leading to motor neuron death (15, 16, 24–27).Accumulating evidence suggests that metal deficiency of SOD1 is an important factor that can influence SOD1 aggregation or neurotoxicity (4, 28–33), but the metal-deficient states of SOD1 that are most relevant to ALS remain unclear. Zinc-deficient, copper-replete SOD1 species, which can be produced in vitro by adding copper to SOD1 that has been stripped of its metal ions at acidic pH, were shown to be toxic to motor neurons in culture (28). However, it has not been shown that zinc-deficient, copper-replete SOD1 is produced in vivo as a consequence of ALS mutations, and loading of copper into SOD1 by the copper chaperone for SOD1 (CCS) is not required for toxicity (34, 35). Furthermore, the MBR mutants have a disrupted copper site and have been found to be severely deficient in both zinc and copper (17, 30), yet expression of these SOD1s still produces motor neuron disease (1, 2, 30, 34, 36, 37).When recombinant human SOD1 was overexpressed in insect cells, we instead observed zinc-replete but copper-deficient species for most WT-like mutants, probably because the capacity of the copper-loading mechanism was exceeded (17). These preparations indicate that zinc can be efficiently incorporated into many WT-like mutants in vivo, and much of it is retained after purification. Furthermore, these copper-deficient biologically metallated proteins may be useful reagents to assess the influence of copper binding upon other properties of SOD1 mutants that may be relevant to their neurotoxicity.We previously observed that reduction of the Cys-57—Cys-146 disulfide bond facilitates the ability of metal chelators to alter the electrophoretic mobility and to increase the hydrophobicity of SOD1 mutants (5). This is consistent with the known properties of this linkage to stabilize the dimeric interface, to orient Arg-143 via a hydrogen bond from the carbonyl oxygen of Cys-57 to Arg-143-NH₂, and to prevent metal ion loss (38–40). However, it remains unclear whether the Cys-57—Cys-146 bond is required to prevent abnormal SOD1 hydrophobic exposure or whether the aberrant conformational change primarily results from metal ion loss. Ablation of the disulfide bond by the experimental (non-ALS) mutants C57S and C146S provides useful reagents to test the relative influence of the disulfide bond and copper binding upon SOD1 properties.In this study we sought to correlate the consequences of copper deficiency, copper and zinc deficiency, and disulfide reduction upon the hydrodynamic behavior and surface hydrophobicity of WT and representative mutant SOD1 enzymes (Fig. 1A). We quantitated the metal contents of as-isolated SOD1 proteins, detected changes in conformation or metal occupancy using native PAGE to assess their electrophoretic mobility, a measure of global conformational change, and correlated these changes to hydrophobic exposure using 1-anilinonaphthalene-8-sulfonic acid (ANS), which is very sensitive to local conformational changes. ANS is a small amphipathic dye (Fig. 1B) that has been used as a sensitive probe to detect hydrophobic pockets on protein surfaces (41–44). Free ANS exhibits only weak fluorescence that is maximal near 520 nm, but when ANS binds to a hydrophobic site in a partially or fully folded protein, the fluorescence peak increases in amplitude and shifts to a shorter wavelength (42). ANS also has an anionic sulfonate group that can interact with cationic groups (e.g. Arg or Lys residues) through ion-pair formation which may be further strengthened by hydrophobic interactions (43–46).Open in a separate windowFIGURE 1.A, WT SOD1 structure showing the position of the C57-C146 intrasubunit disulfide bond (S–S, yellow), bound copper and zinc ions, and ALS mutant residues. The residues altered in A4V, G85R, G93A, D124V, and S134N SOD1s are indicated as green spheres. The backbone of the β-barrel core and the loops is shown in a rainbow color, from blue at the amino terminus to red at the carboxyl terminus. The figure was generated using PyMOL (84) and PDB entry 1HL5 (22). B, chemical structure of ANS fluorophore.To evaluate further the importance of metal ion binding, we measured spectral changes related to the binding of cobalt and copper to the same SOD1 proteins. We observed that as-isolated WT-like mutants containing zinc could interact with copper ions to produce an electrophoretic mobility and decreased hydrophobicity resembling that of the fully metalated holo-WT SOD1. In contrast, we saw no evidence for copper binding to MBR mutants in a manner that alters their hydrodynamic properties or their hydrophobicity. Our data suggest that binding of both copper and zinc are important determinants of SOD1 conformation and that perturbation of such binding may be relevant to the ALS disease process.     

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]     

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]     

Vinyl Sulfones as Antiparasitic Agents and a Structural Basis for Drug Design

Cysteine proteases of the papain superfamily are implicated in a number of cellular processes and are important virulence factors in the pathogenesis of parasitic disease. These enzymes have therefore emerged as promising targets for antiparasitic drugs. We report the crystal structures of three major parasite cysteine proteases, cruzain, falcipain-3, and the first reported structure of rhodesain, in complex with a class of potent, small molecule, cysteine protease inhibitors, the vinyl sulfones. These data, in conjunction with comparative inhibition kinetics, provide insight into the molecular mechanisms that drive cysteine protease inhibition by vinyl sulfones, the binding specificity of these important proteases and the potential of vinyl sulfones as antiparasitic drugs.Sleeping sickness (African trypanosomiasis), caused by Trypanosoma brucei, and malaria, caused by Plasmodium falciparum, are significant, parasitic diseases of sub-Saharan Africa (1). Chagas'' disease (South American trypanosomiasis), caused by Trypanosoma cruzi, affects approximately, 16–18 million people in South and Central America. For all three of these protozoan diseases, resistance and toxicity to current therapies makes treatment increasingly problematic, and thus the development of new drugs is an important priority (2–4).T. cruzi, T. brucei, and P. falciparum produce an array of potential target enzymes implicated in pathogenesis and host cell invasion, including a number of essential and closely related papain-family cysteine proteases (5, 6). Inhibitors of cruzain and rhodesain, major cathepsin L-like papain-family cysteine proteases of T. cruzi and T. brucei rhodesiense (7–10) display considerable antitrypanosomal activity (11, 12), and some classes have been shown to cure T. cruzi infection in mouse models (11, 13, 14).In P. falciparum, the papain-family cysteine proteases falcipain-2 (FP-2)⁶ and falcipain-3 (FP-3) are known to catalyze the proteolysis of host hemoglobin, a process that is essential for the development of erythrocytic parasites (15–17). Specific inhibitors, targeted to both enzymes, display antiplasmodial activity (18). However, although the abnormal phenotype of FP-2 knock-outs is “rescued” during later stages of trophozoite development (17), FP-3 has proved recalcitrant to gene knock-out (16) suggesting a critical function for this enzyme and underscoring its potential as a drug target.Sequence analyses and substrate profiling identify cruzain, rhodesain, and FP-3 as cathepsin L-like, and several studies describe classes of small molecule inhibitors that target multiple cathepsin L-like cysteine proteases, some with overlapping antiparasitic activity (19–22). Among these small molecules, vinyl sulfones have been shown to be effective inhibitors of a number of papain family-like cysteine proteases (19, 23–27). Vinyl sulfones have many desirable attributes, including selectivity for cysteine proteases over serine proteases, stable inactivation of the target enzyme, and relative inertness in the absence of the protease target active site (25). This class has also been shown to have desirable pharmacokinetic and safety profiles in rodents, dogs, and primates (28, 29). We have determined the crystal structures of cruzain, rhodesain, and FP-3 bound to vinyl sulfone inhibitors and performed inhibition kinetics for each enzyme. Our results highlight key areas of interaction between proteases and inhibitors. These results help validate the vinyl sulfones as a class of antiparasitic drugs and provide structural insights to facilitate the design or modification of other small molecule inhibitor scaffolds.     

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]     

Dbf4-Cdc7 Phosphorylation of Mcm2 Is Required for Cell Growth

The Dbf4-Cdc7 kinase (DDK) is required for the activation of the origins of replication, and DDK phosphorylates Mcm2 in vitro. We find that budding yeast Cdc7 alone exists in solution as a weakly active multimer. Dbf4 forms a likely heterodimer with Cdc7, and this species phosphorylates Mcm2 with substantially higher specific activity. Dbf4 alone binds tightly to Mcm2, whereas Cdc7 alone binds weakly to Mcm2, suggesting that Dbf4 recruits Cdc7 to phosphorylate Mcm2. DDK phosphorylates two serine residues of Mcm2 near the N terminus of the protein, Ser-164 and Ser-170. Expression of mcm2-S170A is lethal to yeast cells that lack endogenous MCM2 (mcm2Δ); however, this lethality is rescued in cells harboring the DDK bypass mutant mcm5-bob1. We conclude that DDK phosphorylation of Mcm2 is required for cell growth.The Cdc7 protein kinase is required throughout the yeast S phase to activate origins (1, 2). The S phase cyclin-dependent kinase also activates yeast origins of replication (3–5). It has been proposed that Dbf4 activates Cdc7 kinase in S phase, and that Dbf4 interaction with Cdc7 is essential for Cdc7 kinase activity (6). However, it is not known how Dbf4-Cdc7 (DDK)² acts during S phase to trigger the initiation of DNA replication. DDK has homologs in other eukaryotic species, and the role of Cdc7 in activation of replication origins during S phase may be conserved (7–10).The Mcm2-7 complex functions with Cdc45 and GINS to unwind DNA at a replication fork (11–15). A mutation of MCM5 (mcm5-bob1) bypasses the cellular requirements for DBF4 and CDC7 (16), suggesting a critical physiologic interaction between Dbf4-Cdc7 and Mcm proteins. DDK phosphorylates Mcm2 in vitro with proteins purified from budding yeast (17, 18) or human cells (19). Furthermore, there are mutants of MCM2 that show synthetic lethality with DBF4 mutants (6, 17), suggesting a biologically relevant interaction between DBF4 and MCM2. Nevertheless, the physiologic role of DDK phosphorylation of Mcm2 is a matter of dispute. In human cells, replacement of MCM2 DDK-phosphoacceptor residues with alanines inhibits DNA replication, suggesting that Dbf4-Cdc7 phosphorylation of Mcm2 in humans is important for DNA replication (20). In contrast, mutation of putative DDK phosphorylation sites at the N terminus of Schizosaccharomyces pombe Mcm2 results in viable cells, suggesting that phosphorylation of S. pombe Mcm2 by DDK is not critical for cell growth (10).In budding yeast, Cdc7 is present at high levels in G₁ and S phase, whereas Dbf4 levels peak in S phase (18, 21, 22). Furthermore, budding yeast DDK binds to chromatin during S phase (6), and it has been shown that Dbf4 is required for Cdc7 binding to chromatin in budding yeast (23, 24), fission yeast (25), and Xenopus (9). Human and fission yeast Cdc7 are inert on their own (7, 8), but Dbf4-Cdc7 is active in phosphorylating Mcm proteins in budding yeast (6, 26), fission yeast (7), and human (8, 10). Based on these data, it has been proposed that Dbf4 activates Cdc7 kinase in S phase and that Dbf4 interaction with Cdc7 is essential for Cdc7 kinase activity (6, 9, 18, 21–24). However, a mechanistic analysis of how Dbf4 activates Cdc7 has not yet been accomplished. For example, the multimeric state of the active Dbf4-Cdc7 complex is currently disputed. A heterodimer of fission yeast Cdc7 (Hsk1) in complex with fission yeast Dbf4 (Dfp1) can phosphorylate Mcm2 (7). However, in budding yeast, oligomers of Cdc7 exist in the cell (27), and Dbf4-Cdc7 exists as oligomers of 180 and 300 kDa (27).DDK phosphorylates the N termini of human Mcm2 (19, 20, 28), human Mcm4 (10), budding yeast Mcm4 (26), and fission yeast Mcm6 (10). Although the sequences of the Mcm N termini are poorly conserved, the DDK sites identified in each study have neighboring acidic residues. The residues of budding yeast Mcm2 that are phosphorylated by DDK have not yet been identified.In this study, we find that budding yeast Cdc7 is weakly active as a multimer in phosphorylating Mcm2. However, a low molecular weight form of Dbf4-Cdc7, likely a heterodimer, has a higher specific activity for phosphorylation of Mcm2. Dbf4 or DDK, but not Cdc7, binds tightly to Mcm2, suggesting that Dbf4 recruits Cdc7 to Mcm2. DDK phosphorylates two serine residues of Mcm2, Ser-164 and Ser-170, in an acidic region of the protein. Mutation of Ser-170 is lethal to yeast cells, but this phenotype is rescued by the DDK bypass mutant mcm5-bob1. We conclude that DDK phosphorylation of Ser-170 of Mcm2 is required for budding yeast growth.