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]     

Role of Partitioning-defective 1/Microtubule Affinity-regulating Kinases in the Morphogenetic Activity of Helicobacter pylori CagA

Helicobacter pylori CagA plays a key role in gastric carcinogenesis. Upon delivery into gastric epithelial cells, CagA binds and deregulates SHP-2 phosphatase, a bona fide oncoprotein, thereby causing sustained ERK activation and impaired focal adhesions. CagA also binds and inhibits PAR1b/MARK2, one of the four members of the PAR1 family of kinases, to elicit epithelial polarity defect. In nonpolarized gastric epithelial cells, CagA induces the hummingbird phenotype, an extremely elongated cell shape characterized by a rear retraction defect. This morphological change is dependent on CagA-deregulated SHP-2 and is thus thought to reflect the oncogenic potential of CagA. In this study, we investigated the role of the PAR1 family of kinases in the hummingbird phenotype. We found that CagA binds not only PAR1b but also other PAR1 isoforms, with order of strength as follows: PAR1b > PAR1d ≥ PAR1a > PAR1c. Binding of CagA with PAR1 isoforms inhibits the kinase activity. This abolishes the ability of PAR1 to destabilize microtubules and thereby promotes disassembly of focal adhesions, which contributes to the hummingbird phenotype. Consistently, PAR1 knockdown potentiates induction of the hummingbird phenotype by CagA. The morphogenetic activity of CagA was also found to be augmented through inhibition of non-muscle myosin II. Because myosin II is functionally associated with PAR1, perturbation of PAR1-regulated myosin II by CagA may underlie the defect of rear retraction in the hummingbird phenotype. Our findings reveal that CagA systemically inhibits PAR1 family kinases and indicate that malfunctioning of microtubules and myosin II by CagA-mediated PAR1 inhibition cooperates with deregulated SHP-2 in the morphogenetic activity of CagA.Infection with Helicobacter pylori strains bearing cagA (cytotoxin-associated gene A)-positive strains is the strongest risk factor for the development of gastric carcinoma, the second leading cause of cancer-related death worldwide (1–3). The cagA gene is located within a 40-kb DNA fragment, termed the cag pathogenicity island, which is specifically present in the genome of cagA-positive H. pylori strains (4–6). In addition to cagA, there are ∼30 genes in the cag pathogenicity island, many of which encode a bacterial type IV secretion system that delivers the cagA-encoded CagA protein into gastric epithelial cells (7–10). Upon delivery into gastric epithelial cells, CagA is localized to the plasma membrane, where it undergoes tyrosine phosphorylation at the C-terminal Glu-Pro-Ile-Tyr-Ala motifs by Src family kinases or c-Abl kinase (11–14). The C-terminal Glu-Pro-Ile-Tyr-Ala-containing region of CagA is noted for the structural diversity among distinct H. pylori isolates. Oncogenic potential of CagA has recently been confirmed by a study showing that systemic expression of CagA in mice induces gastrointestinal and hematological malignancies (15).When expressed in gastric epithelial cells, CagA induces morphological transformation termed the hummingbird phenotype, which is characterized by the development of one or two long and thin protrusions resembling the beak of the hummingbird. It has been thought that the hummingbird phenotype is related to the oncogenic action of CagA (7, 16–19). Pathophysiological relevance for the hummingbird phenotype in gastric carcinogenesis has recently been provided by the observation that infection with H. pylori carrying CagA with greater ability to induce the hummingbird phenotype is more closely associated with gastric carcinoma (20–23). Elevated motility of hummingbird cells (cells showing the hummingbird phenotype) may also contribute to invasion and metastasis of gastric carcinoma.In host cells, CagA interacts with the SHP-2 phosphatase, C-terminal Src kinase, and Crk adaptor in a tyrosine phosphorylation-dependent manner (16, 24, 25) and also associates with Grb2 adaptor and c-Met in a phosphorylation-independent manner (26, 27). Among these CagA targets, much attention has been focused on SHP-2 because the phosphatase has been recognized as a bona fide oncoprotein, gain-of-function mutations of which are found in various human malignancies (17, 18, 28). Stable interaction of CagA with SHP-2 requires CagA dimerization, which is mediated by a 16-amino acid CagA-multimerization (CM)² sequence present in the C-terminal region of CagA (29). Upon complex formation, CagA aberrantly activates SHP-2 and thereby elicits sustained ERK MAP kinase activation that promotes mitogenesis (30). Also, CagA-activated SHP-2 dephosphorylates and inhibits focal adhesion kinase (FAK), causing impaired focal adhesions. It has been shown previously that both aberrant ERK activation and FAK inhibition by CagA-deregulated SHP-2 are involved in induction of the hummingbird phenotype (31).Partitioning-defective 1 (PAR1)/microtubule affinity-regulating kinase (MARK) is an evolutionally conserved serine/threonine kinase originally isolated in C. elegans (32–34). Mammalian cells possess four structurally related PAR1 isoforms, PAR1a/MARK3, PAR1b/MARK2, PAR1c/MARK1, and PAR1d/MARK4 (35–37). Among these, PAR1a, PAR1b, and PAR1c are expressed in a variety of cells, whereas PAR1d is predominantly expressed in neural cells (35, 37). These PAR1 isoforms phosphorylate microtubule-associated proteins (MAPs) and thereby destabilize microtubules (35, 38), allowing asymmetric distribution of molecules that are involved in the establishment and maintenance of cell polarity.In polarized epithelial cells, CagA disrupts the tight junctions and causes loss of apical-basolateral polarity (39, 40). This CagA activity involves the interaction of CagA with PAR1b/MARK2 (19, 41). CagA directly binds to the kinase domain of PAR1b in a tyrosine phosphorylation-independent manner and inhibits the kinase activity. Notably, CagA binds to PAR1b via the CM sequence (19). Because PAR1b is present as a dimer in cells (42), CagA may passively homodimerize upon complex formation with the PAR1 dimer via the CM sequence, and this PAR1-directed CagA dimer would form a stable complex with SHP-2 through its two SH2 domains.Because of the critical role of CagA in gastric carcinogenesis (7, 16–19), it is important to elucidate the molecular basis underlying the morphogenetic activity of CagA. In this study, we investigated the role of PAR1 isoforms in induction of the hummingbird phenotype by CagA, and we obtained evidence that CagA-mediated inhibition of PAR1 kinases contributes to the development of the morphological change by perturbing microtubules and non-muscle myosin II.     

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.     

PTIP Regulates 53BP1 and SMC1 at the DNA Damage Sites

Although PTIP is implicated in the DNA damage response, through interactions with 53BP1, the function of PTIP in the DNA damage response remain elusive. Here, we show that RNF8 controls DNA damage-induced nuclear foci formation of PTIP, which in turn regulates 53BP1 localization to the DNA damage sites. In addition, SMC1, a substrate of ATM, could not be phosphorylated at the DNA damage sites in the absence of PTIP. The PTIP-dependent pathway is important for DNA double strand breaks repair and DNA damage-induced intra-S phase checkpoint activation. Taken together, these results suggest that the role of PTIP in the DNA damage response is downstream of RNF8 and upstream of 53BP1. Thus, PTIP regulates 53BP1-dependent signaling pathway following DNA damage.The DNA damage response pathways are signal transduction pathways with DNA damage sensors, mediators, and effectors, which are essential for maintaining genomic stability (1–3). Following DNA double strand breaks, histone H2AX at the DNA damage sites is rapidly phosphorylated by ATM/ATR/DNAPK (4–10), a family homologous to phosphoinositide 3-kinases (11, 12). Subsequently, phospho-H2AX (γH2AX) provides the platform for accumulation of a larger group of DNA damage response factors, such as MDC1, BRCA1, 53BP1, and the MRE11·RAD50·NBS1 complex (13, 14), at the DNA damage sites. Translocalization of these proteins to the DNA double strand breaks (DSBs)³ facilitates DNA damage checkpoint activation and enhances the efficiency of DNA damage repair (14, 15).Recently, PTIP (Pax2 transactivation domain-interacting protein, or Paxip) has been identified as a DNA damage response protein and is required for cell survival when exposed to ionizing radiation (IR) (1, 16–18). PTIP is a 1069-amino acid nuclear protein and has been originally identified in a yeast two-hybrid screening as a partner of Pax2 (19). Genetic deletion of the PTIP gene in mice leads to early embryonic lethality at embryonic day 8.5, suggesting that PTIP is essential for early embryonic development (20). Structurally, PTIP contains six tandem BRCT (BRCA1 carboxyl-terminal) domains (16–18, 21). The BRCT domain is a phospho-group binding domain that mediates protein-protein interactions (17, 22, 23). Interestingly, the BRCT domain has been found in a large number of proteins involved in the cellular response to DNA damages, such as BRCA1, MDC1, and 53BP1 (7, 24–29). Like other BRCT domain-containing proteins, upon exposure to IR, PTIP forms nuclear foci at the DSBs, which is dependent on its BRCT domains (16–18). By protein affinity purification, PTIP has been found in two large complexes. One includes the histone H3K4 methyltransferase ALR and its associated cofactors, the other contains DNA damage response proteins, including 53BP1 and SMC1 (30, 31). Further experiments have revealed that DNA damage enhances the interaction between PTIP and 53BP1 (18, 31).To elucidate the DNA damage response pathways, we have examined the upstream and downstream partners of PTIP. Here, we report that PTIP is downstream of RNF8 and upstream of 53BP1 in response to DNA damage. Moreover, PTIP and 53BP1 are required for the phospho-ATM association with the chromatin, which phosphorylates SMC1 at the DSBs. This PTIP-dependent pathway is involved in DSBs repair.     

Phosphoproteomics of Klebsiella pneumoniae NTUH-K2044 Reveals a Tight Link between Tyrosine Phosphorylation and Virulence

Encapsulated Klebsiella pneumoniae is the predominant causative agent of pyogenic liver abscess, an emerging infectious disease that often complicates metastatic meningitis or endophthalmitis. The capsular polysaccharide on K. pneumoniae surface was determined as the key to virulence. Although the regulation of capsular polysaccharide biosynthesis is largely unclear, it was found that protein-tyrosine kinases and phosphatases are involved. Therefore, the identification and characterization of such kinases, phosphatases, and their substrates would advance our knowledge of the underlying mechanism in capsule formation and could contribute to the development of new therapeutic strategies. Here, we analyzed the phosphoproteome of K. pneumoniae NTUH-K2044 with a shotgun approach and identified 117 unique phosphopeptides along with 93 in vivo phosphorylated sites corresponding to 81 proteins. Interestingly, three of the identified tyrosine phosphorylated proteins, namely protein-tyrosine kinase (Wzc), phosphomannomutase (ManB), and undecaprenyl-phosphate glycosyltransferase (WcaJ), were found to be distributed in the cps locus and thus were speculated to be involved in the converging signal transduction of capsule biosynthesis. Consequently, we decided to focus on the lesser studied ManB and WcaJ for mutation analysis. The capsular polysaccharides of WcaJ mutant (WcaJY5F) were dramatically reduced quantitatively, and the LD₅₀ increased by 200-fold in a mouse peritonitis model compared with the wild-type strain. However, the capsular polysaccharides of ManB mutant (ManBY26F) showed no difference in quantity, and the LD₅₀ increased by merely 6-fold in mice test. Our study provided a clear trend that WcaJ tyrosine phosphorylation can regulate the biosynthesis of capsular polysaccharides and result in the pathogenicity of K. pneumoniae NTUH-K2044.Protein phosphorylation is one of the most biologically relevant and ubiquitous post-translational modifications in both eukaryotic and prokaryotic organisms. It is best known that protein phosphorylation is a reversible enzyme-catalyzed process that is controlled by various kinases and phosphatases. The aberrant functions often result in irregular protein phosphorylation and ultimately lead to serious disease states such as malignant transformation, immune disorders, and pathogenic infections in mammals (1, 2). Recently, accumulating evidences suggest that Ser/Thr/Tyr phosphorylations also contribute to regulate a diverse range of cellular responses and physiological processes in prokaryotes (1). Among them, tyrosine phosphorylation in encapsulated bacteria has been discovered to play key roles in capsular polysaccharide (CPS1; K antigen) biosynthesis, which leads to virulence (3, 4). This thick layer of exopolysaccharide on many pathogenic bacteria can act as a physical boundary to evade phagocytosis and complement-mediated killing and further inhibit complement activation of the host (1, 5, 6).In 1996, Acinetobacter johnsonii protein-tyrosine kinase (Ptk) was first discovered and categorized under the bacterial protein-tyrosine kinase (BY-kinase) family (1, 7, 8). Shortly after, its function in bacterial exopolysaccharide production and transport was characterized (1, 7, 8). From then on, many more bacterial tyrosine kinases such as Wzc of Escherichia coli (1, 9) and EpsB of Pseudomonas solanacearum (10, 11) were found to possess this conserved property; deletion of such tyrosine kinases will result in the loss of exopolysaccharide production (12). Therefore, several experiments were conducted to investigate the role of the downstream substrates of the tyrosine kinases in different strains of bacteria, and some targeted proteins were found to participate in the exopolysaccharide anabolism (13, 14). These findings demonstrated a direct relationship between bacterial tyrosine phosphorylation and exopolysaccharide biosynthesis that was directly reflected in the strain virulence.In the past, the functional roles of the critical components involved in protein phosphorylation were defined by basic biochemical and genetic approaches (1). However, there exists a salient gap between the growing number of identified protein-tyrosine kinases/phosphatases and the relative paucity of protein substrates characterized to date. Genomic sequence analyses and advanced high resolution/high accuracy MS systems with vastly improved phosphopeptide enrichment strategies are among the two key enabling technologies that allow a high efficiency identification of the scarcely detectable site-specific phosphorylations in bacterial systems (15). Mann et al. (16) were the first to initiate a systematic study of the phosphoproteome of B. subtilis in 2007 followed by similar site-specific phosphoproteomics analyses of E. coli (17), Lactococcus lactis (18), and Halobacterium salinarum (19). These pioneering works have since set the foundation in bacterial phosphoproteomics but have not been specifically carried out to address a particular biological issue of causal relevance to virulence or pathogenesis.Klebsiella pneumoniae is a Gram-negative, non-motile, facultative anaerobic, and rod-shaped bacterium. It is commonly found in water and soil (20) as well as on plants (21) and mucosal surfaces of mammals, such as human, horse, and swine (22, 23). It was demonstrated that CPS on the surface of K. pneumoniae is the prime factor of virulence and toxicity in causing pyogenic liver abscess (PLA), a common intra-abdominal infection with a high 10–30% mortality rate worldwide (24–29). There are also variations in virulence in regard to different capsular serotypes; K1 and K2 were found to be especially pathogenic in causing PLA in a mouse model (30) compared with other serotypes, which show little or no effect (31–34). The K. pneumoniae NTUH-K2044 (K2044) strain, encapsulated with K1 antigen (35), was isolated from clinical K. pneumoniae liver abscess patients. It has become an important emerging pathogen (36) because it usually complicates metastatic septic endophthalmitis and irreversible central nervous system infections independent of host underlying diseases (30, 34). The transmission rate is high (37), and it often rapidly leads to outbreaks of community-acquired infections, such as bacteremia, nosocomial pneumonia, and sepsis, common in immunocompromised individuals (38).In this study, we wanted to prove that the biosynthesis of CPS is mediated through tyrosine phosphorylation of a subset of proteins. An MS-based systematic phosphoproteomics analysis was conducted on K2044 to identify tyrosine phosphorylated proteins that are also associated with CPS biosynthesis. We further validated the relationship between tyrosine phosphorylation on those proteins and virulence of K2044 by site-directed mutagenesis, CPS quantification, serum killing, and mouse lethality assay.