Context-dependent Bcl-2/Bak Interactions Regulate Lymphoid Cell Apoptosis

The release of cytochrome c from mitochondria, which leads to activation of the intrinsic apoptotic pathway, is regulated by interactions of Bax and Bak with antiapoptotic Bcl-2 family members. The factors that regulate these interactions are, at the present time, incompletely understood. Recent studies showing preferences in binding between synthetic Bcl-2 homology domain 3 and antiapoptotic Bcl-2 family members in vitro have suggested that the antiapoptotic proteins Mcl-1 and Bcl-x_L, but not Bcl-2, restrain proapoptotic Bak from inducing mitochondrial membrane permeabilization and apoptosis. Here we show that Bak protein has a much higher affinity than the 26-amino acid Bak Bcl-2 homology domain 3 for Bcl-2, that some naturally occurring Bcl-2 allelic variants have an affinity for full-length Bak that is only 3-fold lower than that of Mcl-1, and that endogenous levels of these Bcl-2 variants (which are as much as 40-fold more abundant than Mcl-1) restrain part of the Bak in intact lymphoid cells. In addition, we demonstrate that Bcl-2 variants can, depending on their affinity for Bak, substitute for Mcl-1 in protecting cells. Thus, the ability of Bcl-2 to protect cells from activated Bak depends on two important contextual variables, the identity of the Bcl-2 present and the amount expressed.The release of cytochrome c from mitochondria, which leads to activation of the intrinsic apoptotic pathway, is regulated by Bcl-2 family members (1–5). This group of proteins consists of three subgroups: Bax and Bak, which oligomerize upon death stimulation to form a putative pore in the outer mitochondrial membrane, thereby allowing efflux of cytochrome c and other mitochondrial intermembrane space components; Bcl-2, Bcl-x_L, Mcl-1, and other antiapoptotic homologs, which antagonize the effects of Bax and Bak; and BH3-only proteins² such as Bim, Bid, and Puma, which are proapoptotic Bcl-2 family members that share only limited homology with the other two groups in a single 15-amino acid domain (the BH3 domain, see Ref. 6). Although it is clear that BH3-only proteins serve as molecular sensors of various stresses and, when activated, trigger apoptosis (3, 6–11), the mechanism by which they do so remains incompletely understood. One current model suggests that BH3-only proteins trigger apoptosis solely by binding and neutralizing antiapoptotic Bcl-2 family members, thereby causing them to release the activated Bax and Bak that are bound (reviewed in Refs. 9 and 10; see also Refs. 12 and 13), whereas another current model suggests that certain BH3-only proteins also directly bind to and activate Bax (reviewed in Ref. 3; see also Refs. 14–17). Whichever model turns out to be correct, both models agree that certain antiapoptotic Bcl-2 family members can inhibit apoptosis, at least in part, by binding and neutralizing activated Bax and Bak before they permeabilize the outer mitochondrial membrane (13, 18, 19).Much of the information about the interactions between pro- and antiapoptotic Bcl-2 family members has been derived from the study of synthetic peptides corresponding to BH3 domains. In particular, these synthetic peptides have been utilized as surrogates for the full-length proapoptotic proteins during structure determinations (20–22) as well as in functional studies exploring the effect of purified BH3 domains on isolated mitochondria (14, 23) and on Bax-mediated permeabilization of lipid vesicles (15).Recent studies using these same peptides have suggested that interactions of the BH3 domains of Bax, Bak, and the BH3-only proteins with the “BH3 receptors” of the antiapoptotic Bcl-2 family members are not all equivalent. Surface plasmon resonance, a technique that is widely used to examine the interactions of biomolecules under cell-free conditions (24–26), has demonstrated that synthetic BH3 peptides of some BH3-only family members show striking preferences, with the Bad BH3 peptide binding to Bcl-2 and Bcl-x_L but not Mcl-1, and the Noxa BH3 peptide binding to Mcl-1 but not Bcl-2 or Bcl-x_L (27). Likewise, the Bak BH3 peptide exhibits selectivity, with high affinity for Bcl-x_L and Mcl-1 but not Bcl-2 (12). The latter results have led to a model in which Bcl-x_L and Mcl-1 restrain Bak and inhibit Bak-dependent apoptosis, whereas Bcl-2 does not (10).Because the Bak protein contains multiple recognizable domains in addition to its BH3 motif (28, 29), we compared the binding of Bak BH3 peptide and Bak protein to Bcl-2. Surface plasmon resonance demonstrated that Bcl-2 binds Bak protein with much higher affinity than the Bak 26-mer BH3 peptide. Further experiments demonstrated that the K_D for Bak differs among naturally occurring Bcl-2 sequence variants but is only 3-fold higher than that of Mcl-1 in some cases. In light of previous reports that Bcl-2 overexpression contributes to neoplastic transformation (30–33) and drug resistance (34–36) in lymphoid cells, we also examined Bcl-2 expression and Bak binding in a panel of neoplastic lymphoid cell lines. Results of these experiments demonstrated that Bcl-2 expression varies among different lymphoid cell lines but is up to 40-fold more abundant than Mcl-1. In lymphoid cell lines with abundant Bcl-2, Bak is detected in Bcl-2 as well as Mcl-1 immunoprecipitates; and Bak-dependent apoptosis induced by Mcl-1 down-regulation can be prevented by Bcl-2 overexpression. Collectively, these observations shed new light on the role of Bcl-2 in binding and neutralizing Bak.     

Mcl-1 Degradation during Hepatocyte Lipoapoptosis

The mechanisms of free fatty acid-induced lipoapoptosis are incompletely understood. Here we demonstrate that Mcl-1, an anti-apoptotic member of the Bcl-2 family, was rapidly degraded in hepatocytes in response to palmitate and stearate by a proteasome-dependent pathway. Overexpression of a ubiquitin-resistant Mcl-1 mutant in Huh-7 cells attenuated palmitate-mediated Mcl-1 loss and lipoapoptosis; conversely, short hairpin RNA-targeted knockdown of Mcl-1 sensitized these cells to lipoapoptosis. Palmitate-induced Mcl-1 degradation was attenuated by the novel protein kinase C (PKC) inhibitor rottlerin. Of the two human novel PKC isozymes, PKCδ and PKCθ, only activation of PKCθ was observed by phospho-immunoblot analysis. As compared with Jurkat cells, a smaller PKCθ polypeptide and mRNA were expressed in hepatocytes consistent with an alternative splice variant. Short hairpin RNA-mediated knockdown of PKCθ reduced Mcl-1 degradation and lipoapoptosis. Likewise, genetic deletion of Pkcθ also attenuated Mcl-1 degradation and cytotoxicity by palmitate in primary hepatocytes. During treatment with palmitate, rottlerin inhibited phosphorylation of Mcl-1 at Ser¹⁵⁹, a phosphorylation site previously implicated in Mcl-1 turnover. Consistent with these results, an Mcl-1 S159A mutant was resistant to degradation and improved cell survival during palmitate treatment. Collectively, these results implicate PKCθ-dependent destabilization of Mcl-1 as a mechanism contributing to hepatocyte lipoapoptosis.Current evidence suggests that hepatic steatosis is present in up to 30% of the American population (1). A subset of these individuals develop severe hepatic lipotoxicity, a syndrome referred to as NASH² (2), which can progress to cirrhosis and its chronic sequela (3, 4). A major risk factor for hepatic lipotoxicity is insulin resistance (5–7), resulting in excessive lipolysis within peripheral adipose tissue with release of high levels of free fatty acids (FFA) to the circulation. Circulating FFA are taken up by the liver via fatty acid transporter 5 and CD36 (8–10), and the bulk of hepatic neutral fat is derived from re-esterification of circulating FFA (8). Current concepts indicate that FFA, and not their esterified product (triglyceride), mediate hepatic lipotoxicity (11, 12). Elevated serum FFA correlate with liver disease severity (13–15), and therapies that enhance insulin sensitivity ameliorate hepatic lipotoxicity, in part, by decreasing plasma FFA (16). Hepatic FFA also accumulate in experimental steatohepatitis, further supporting a role for these nutrients in hepatic lipotoxicity (17). Saturated FFA are more strongly implicated in hepatic lipotoxicity than unsaturated FFA (18, 19). Saturated FFA induce hepatocyte apoptosis (20, 21), a cardinal feature of nonalcoholic fatty liver disease (22), and serum biomarkers of apoptosis are useful for identifying hepatic lipotoxicity (23). Thus, FFA-mediated lipotoxicity occurs, in part, by apoptosis.Apoptosis is regulated by members of the Bcl-2 protein family (24). These proteins can be categorized into three subsets as follows: the guardians or anti-apoptotic members of this family, which include Bcl-2, A1, Mcl-1, Bcl-x_L, and Bcl-w; the multidomain executioners or proapoptotic members of this family, which include Bax and Bak; and the messengers or biosensors of cell death, which share only the third Bcl-2 homology domain and are referred to as BH3-only proteins. This last group of proteins includes Bid, Bim, Bmf, Puma, Noxa, Hrk, Bad, and Bik. We have previously reported that cytotoxic FFA induce Bim expression by a FoxO3a-dependent mechanism that contributes, in part, to lipoapoptosis by activating Bax (20, 21). However, Bax activation can be held in check by anti-apoptotic members of the Bcl-2 family suggesting their function may also be dysregulated during FFA-mediated cytotoxicity.Bcl-2 is not expressed in hepatocytes at the protein level (25), whereas Bcl-w and Bfl-1/A1 knock-out mice have no liver phenotype (26–28). However, both potent anti-apoptotic proteins Bcl-x_L and Mcl-1 are expressed by hepatocytes and exhibit a liver phenotype in knock-out mice (29, 30), whereas up-regulation of Mcl-1 renders hepatocytes resistant to apoptosis (31–33). It has also been posited that cellular elimination of Mcl-1 is a critical step in certain proapoptotic cascades (34, 35). Mcl-1 is unique among Bcl-2 proteins in that it has a short half-life, 30–120 min in most cell types, due to the presence of two sequences rich in proline, glutamic acid, serine, and threonine, which target the protein for rapid degradation by the proteasome (36). Proteasomal degradation of Mcl-1 is promoted by ubiquitination, which in turn is regulated by various kinase cascades (36). Despite its potential importance, a role for Mcl-1 in regulating hepatocyte FFA-mediated lipoapoptosis remains unexplored.Given that FFA induce insulin resistance (37), the kinases potentially regulating lipoapoptosis are likely those also identified in insulin resistance syndromes, especially the novel PKC isoforms PKCδ and PKCθ (38). The novel PKC isoforms are activated by diacylglycerol, which rises in the presence of FFA (39–41), and diacylglycerol levels are significantly increased in NASH (42). A role for PKCδ in apoptosis has not been described. PKCθ has recently been shown to be activated by endoplasmic reticulum stress in liver cells (43) and lipids in vivo (44, 45). Furthermore, PKCθ has also been implicated in apoptosis of Jurkat cells, neuroblastoma cells, and myeloid leukemia cells (46, 47). However, neither its role in mediating lipoapoptosis nor modulating levels/activity of Bcl-2 proteins has been examined.This study addresses the role of Mcl-1 and PKCθ in FFA-induced lipoapoptosis. We identify a pathway that involves PKCθ-dependent proteasomal degradation of Mcl-1. Using inhibitors of various steps along this pathway, along with Mcl-1 mutants that are resistant to proteasomal degradation or Ser¹⁵⁹ phosphorylation, our studies implicate Mcl-1 degradation via a PKCθ-dependent process as a critical step in lipoapoptosis.     

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]     

Proinsulin and the Genetics of Diabetes Mellitus

Insulin plays a central role in the regulation of vertebrate metabolism. The hormone, the post-translational product of a single-chain precursor, is a globular protein containing two chains, A (21 residues) and B (30 residues). Recent advances in human genetics have identified dominant mutations in the insulin gene causing permanent neonatal-onset DM² (1–4). The mutations are predicted to block folding of the precursor in the ER of pancreatic β-cells. Although expression of the wild-type allele would in other circumstances be sufficient to maintain homeostasis, studies of a corresponding mouse model (5–7) suggest that the misfolded variant perturbs wild-type biosynthesis (8, 9). Impaired β-cell secretion is associated with ER stress, distorted organelle architecture, and cell death (10). These findings have renewed interest in insulin biosynthesis (11–13) and the structural basis of disulfide pairing (14–19). Protein evolution is constrained not only by structure and function but also by susceptibility to toxic misfolding.Insulin plays a central role in the regulation of vertebrate metabolism. The hormone, the post-translational product of a single-chain precursor, is a globular protein containing two chains, A (21 residues) and B (30 residues). Recent advances in human genetics have identified dominant mutations in the insulin gene causing permanent neonatal-onset DM² (1–4). The mutations are predicted to block folding of the precursor in the ER of pancreatic β-cells. Although expression of the wild-type allele would in other circumstances be sufficient to maintain homeostasis, studies of a corresponding mouse model (5–7) suggest that the misfolded variant perturbs wild-type biosynthesis (8, 9). Impaired β-cell secretion is associated with ER stress, distorted organelle architecture, and cell death (10). These findings have renewed interest in insulin biosynthesis (11–13) and the structural basis of disulfide pairing (14–19). Protein evolution is constrained not only by structure and function but also by susceptibility to toxic misfolding.     

Arachidonic Acid Stimulates Cell Adhesion through a Novel p38 MAPK-RhoA Signaling Pathway That Involves Heat Shock Protein 27

Rho GTPases are critical components of cellular signal transduction pathways. Both hyperactivity and overexpression of these proteins have been observed in human cancers and have been implicated as important factors in metastasis. We previously showed that dietary n-6 fatty acids increase cancer cell adhesion to extracellular matrix proteins, such as type IV collagen. Here we report that in MDA-MB-435 human melanoma cells, arachidonic acid activates RhoA, and inhibition of RhoA signaling with either C3 exoenzyme or dominant negative Rho blocked arachidonic acid-induced cell adhesion. Inhibition of the Rho kinase (ROCK) with either small molecule inhibitors or ROCK II-specific small interfering RNA (siRNA) blocked the fatty acid-induced adhesion. However, unlike other systems, inhibition of ROCK did not block the activation of p38 mitogen-activated protein kinase (MAPK); instead, Rho activation depended on p38 MAPK activity and the presence of heat shock protein 27 (HSP27), which is phosphorylated downstream of p38 after arachidonic acid treatment. HSP27 associated with p115RhoGEF in fatty acid-treated cells, and this association was blocked when p38 was inhibited. Furthermore, siRNA knockdown of HSP27 blocked the fatty acid-stimulated Rho activity. Expression of dominant negative p115-RhoGEF or p115RhoGEF-specific siRNA inhibited both RhoA activation and adhesion on type IV collagen, whereas a constitutively active p115RhoGEF restored the arachidonic acid stimulation in cells in which the p38 MAPK had been inhibited. These data suggest that n-6 dietary fatty acids stimulate a set of interactions that regulates cell adhesion through RhoA and ROCK II via a p38 MAPK-dependent association of HSP27 and p115RhoGEF.The ability of tumor cells to metastasize to secondary sites is a hallmark of neoplastic disease. Unfortunately, this propensity to spread is the primary cause of morbidity and death in cancer patients (1). Metastasis is clearly a highly regulated, multistep process that occurs in a spatiotemporal manner (2–4). To escape the restrictive compartment boundaries characteristic of adult tissue, separate intravasation and extravasation steps requiring alterations in co-adhesion, adhesion, invasion, and migration must occur. Execution of these biological processes, involving multiple proteins and cellular organelles, require highly coordinated cell signaling mechanisms.The Rho family of small GTPases regulates many facets of cytoskeletal rearrangements that facilitate cell attachment and migration (5–7). Rho GTPases act as molecular switches by changing from an inactive GDP-bound conformation to an active GTP-bound conformation, thereby regulating a signaling pathway. These proteins are directly regulated by Rho guanine nucleotide exchange factors (GEFs),² Rho GTPase activating proteins, and Rho GDP-dissociation inhibitors (8–12). RhoGEFs bind to the GTPase to catalyze the dissociation of GDP, allowing the binding of GTP and thereby promoting Rho activation (8). The RGS (regulators of G protein signaling) domain-containing RhoGEFs are a recently described family of GEFs. Currently, there are three members of this family, PDZ-RhoGEF, LARG, and p115RhoGEF (13–15), in which the RGS domains function as a heterotrimeric GTPase-activating domain (13, 15, 16). The RGS family of RhoGEFs has been shown to regulate Rho during several processes including cytoskeletal rearrangements, cell adhesion, and cancer progression (17–21).There is significant interplay between the activity of small GTPases and signaling derived from fatty acid metabolism (22–28). Linoleic acid, which is metabolized to arachidonic acid, is an n-6 polyunsaturated fatty acid that is present at high levels in most western diets (29). In animal models, diets high in n-6 polyunsaturated fatty acids have been shown to enhance tumor progression and metastasis (30, 31). Additionally, arachidonic acid is stored in cell membranes and is made available by phospholipases under conditions of increased inflammatory response (32). Arachidonic acid is further metabolized by cyclooxygenases (COX), lipoxygenases (LOX), and cytochrome P450 monooxygenases to yield bioactive products that have myriad effects on cells, and altered metabolism of arachidonic acid by COX, LOX, and P450 has been implicated in cancer progression (31, 33–36).We have studied mechanisms of cell adhesion using the MDA-MB-435 cells as a model of a highly metastatic human cancer cell line (37). These cells have been extensively studied for their ability to recapitulate the metastatic cascade in vivo and in vitro, although recent work indicates that the cells currently in use are most likely a human melanoma line (38). We initially observed that arachidonic acid (AA) enhanced adhesion of MDA-MB-435 cells to type IV collagen through specific integrin-mediated pathways (37). Exogenous AA led to the activation of mitogen-activated protein kinase (MAPK)-activated protein kinase 2 and the phosphorylation of heat shock protein 27 (HSP27) via a p38 MAPK-dependent process (39). Inhibition of p38 MAPK activation blocked cell adhesion as did function-blocking antibodies specific for subunits of the collagen receptor (40). More recently, we identified the key metabolite of AA (15-(S)- hydroxyeicosatetraenoic acid) and the upstream kinases (TAK1 and MKK6) that are responsible for activation of p38 MAPK in this system (41).In this study we investigated the role of Rho activation in the MDA-MB-435 cells after exposure to arachidonic acid. Several aspects of the regulation of Rho signaling in these cells provide insights into the cross-talk between important signaling pathways.     

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]     

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]     

Multipattern Consensus Regions in Multiple Aligned Protein Sequences and Their Segmentation

Decomposing a biological sequence into its functional regions is an important prerequisite to understand the molecule. Using the multiple alignments of the sequences, we evaluate a segmentation based on the type of statistical variation pattern from each of the aligned sites. To describe such a more general pattern, we introduce multipattern consensus regions as segmented regions based on conserved as well as interdependent patterns. Thus the proposed consensus region considers patterns that are statistically significant and extends a local neighborhood. To show its relevance in protein sequence analysis, a cancer suppressor gene called p53 is examined. The results show significant associations between the detected regions and tendency of mutations, location on the 3D structure, and cancer hereditable factors that can be inferred from human twin studies.[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]