Defining the Structural Basis of Human Plasminogen Binding by Streptococcal Surface Enolase

The flesh-eating bacterium group A Streptococcus (GAS) binds and activates human plasminogen, promoting invasive disease. Streptococcal surface enolase (SEN), a glycolytic pathway enzyme, is an identified plasminogen receptor of GAS. Here we used mass spectrometry (MS) to confirm that GAS SEN is octameric, thereby validating in silico modeling based on the crystal structure of Streptococcus pneumoniae α-enolase. Site-directed mutagenesis of surface-located lysine residues (SEN^{K252 + 255A}, SEN^K304A, SEN^K334A, SEN^K344E, SEN^K435L, and SEN^Δ434–435) was used to examine their roles in maintaining structural integrity, enzymatic function, and plasminogen binding. Structural integrity of the GAS SEN octamer was retained for all mutants except SEN^K344E, as determined by circular dichroism spectroscopy and MS. However, ion mobility MS revealed distinct differences in the stability of several mutant octamers in comparison with wild type. Enzymatic analysis indicated that SEN^K344E had lost α-enolase activity, which was also reduced in SEN^K334A and SEN^Δ434–435. Surface plasmon resonance demonstrated that the capacity to bind human plasminogen was abolished in SEN^{K252 + 255A}, SEN^K435L, and SEN^Δ434–435. The lysine residues at positions 252, 255, 434, and 435 therefore play a concerted role in plasminogen acquisition. This study demonstrates the ability of combining in silico structural modeling with ion mobility-MS validation for undertaking functional studies on complex protein structures.Streptococcus pyogenes (group A Streptococcus, GAS)⁸ is a common bacterial pathogen, causing over 700 million human disease episodes each year (1). These range from serious life-threatening invasive diseases including necrotizing fasciitis and streptococcal toxic shock-like syndrome to non-invasive infections like pharyngitis and pyoderma. Invasive disease, in combination with postinfection immune sequelae including rheumatic heart disease and acute poststreptococcal glomerulonephritis, account for over half a million deaths each year (1). Although a resurgence of GAS invasive infections has occurred in western countries since the mid-1980s, disease burden is much greater in developing countries and indigenous populations of developed nations, where GAS infections are endemic (2–4).GAS is able to bind human plasminogen and activate the captured zymogen to the serine protease plasmin (5–17). The capacity of GAS to do this plays a critical role in virulence and invasive disease initiation (3, 17–19). The plasminogen activation system in humans is an important and highly regulated process that is responsible for breakdown of extracellular matrix components, dissolution of blood clots, and cell migration (20, 21). Plasminogen is a 92-kDa zymogen that circulates in human plasma at a concentration of 2 μm (22). It consists of a binding region of five homologous triple loop kringle domains and an N-terminal serine protease domain that flank the Arg⁵⁶¹–Val⁵⁶² site (23), where it is cleaved by tissue plasminogen activator and urokinase plasminogen activator to yield the active protease plasmin (20, 23). GAS also has the ability to activate human plasminogen by secreting the virulence determinant streptokinase. Streptokinase forms stable complexes with plasminogen or plasmin, both of which exhibit plasmin activity (20, 24). Activation of plasminogen by the plasmin(ogen)-streptokinase complex circumvents regulation by the host plasminogen activation inhibitors, α₂-antiplasmin and α₂-macroglobulin (11, 20). GAS can bind the plasmin(ogen)-streptokinase complex and/or plasmin(ogen) directly via plasmin(ogen) receptors at the bacterial cell surface (6). These receptors include the plasminogen-binding group A streptococcal M-like protein (PAM) (25), the PAM-related protein (19), glyceraldehyde-3-phosphate dehydrogenase (GAPDH; also known as streptococcal plasmin receptor, Plr, or streptococcal surface dehydrogenase) (9, 26), and streptococcal surface enolase (SEN or α-enolase) (27). Interactions with these GAS receptors occurs via lysine-binding sites within the kringle domains of plasminogen (6).In addition to its ability to bind human plasminogen, SEN is primarily the glycolytic enzyme that converts 2-phosphoglycerate to phosphoenolpyruvate (27–29). SEN is abundantly expressed in the cytosol of most bacterial species but has also been identified as a surface-located protein in GAS and other bacteria including pneumococci, despite lacking classical cell surface protein motifs such as a signal sequence, membrane-spanning domain, or cell-wall anchor motif (27, 28, 30, 31). The interaction between SEN and plasminogen is reported to be facilitated by the two C-terminal lysine residues at positions 434 and 435 (27, 32). In contrast, an internal binding motif containing lysines at positions 252 and 255 in the closely related α-enolase of Streptococcus pneumoniae has been shown to play a pivotal role in the acquisition of plasminogen in this bacterial species (33). The octameric pneumococcal α-enolase structure consists of a tetramer of dimers. Hence, potential binding sites could be buried in the interface between subunits. In fact, the crystal structure of S. pneumoniae α-enolase revealed that the two C-terminal lysine residues are significantly less exposed than the internal plasminogen-binding motif (34).In this study, we constructed an in silico model of GAS SEN, based on the pneumococcal octameric α-enolase crystal structure, and validated this model using ion mobility (IM) mass spectrometry (MS). Site-directed mutagenesis followed by structural and functional analyses revealed that Lys³⁴⁴ plays a crucial role in structural integrity and enzymatic function. Furthermore, we demonstrate that the plasminogen-binding motif residues Lys²⁵² and Lys²⁵⁵ and the C-terminal Lys⁴³⁴ and Lys⁴³⁵ residues are located adjacently in the GAS SEN structure and play a concerted role in the binding of human plasminogen.     

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.     

Delineation of a Carcinogenic Helicobacter pylori Proteome

Helicobacter pylori is the strongest known risk factor for gastric adenocarcinoma, yet only a fraction of infected persons ever develop cancer. The extensive genetic diversity inherent to this pathogen has precluded comprehensive analyses of constituents that mediate carcinogenesis. We previously reported that in vivo adaptation of a non-carcinogenic H. pylori strain endowed the output derivative with the ability to induce adenocarcinoma, providing a unique opportunity to identify proteins selectively expressed by an oncogenic H. pylori strain. Using a global proteomics DIGE/MS approach, a novel missense mutation of the flagellar protein FlaA was identified that affects structure and function of this virulence-related organelle. Among 25 additional differentially abundant proteins, this approach also identified new proteins previously unassociated with gastric cancer, generating a profile of H. pylori proteins to use in vaccine development and for screening persons infected with strains most likely to induce severe disease.Helicobacter pylori is a Gram-negative bacterial species that selectively colonizes gastric epithelium and induces an inflammatory response within the stomach that persists for decades (1, 2). Biological costs incurred by the long term relationship between H. pylori and humans include an increased risk for distal gastric adenocarcinoma (3–8), and eradication of this pathogen significantly decreases cancer risk among infected individuals without premalignant lesions (9). However, only a fraction of colonized persons ever develop neoplasia, and enhanced cancer risk is related to H. pylori strain differences, inflammatory responses governed by host genetic diversity, and/or specific interactions between host and microbial determinants (10).H. pylori strains are remarkably diverse (11–15), and the genetic composition of strains can change over time within an individual colonized stomach (16, 17). Despite this diversity, several genetic loci have been identified that augment disease risk. The cag pathogenicity island encodes a type IV bacterial secretion system, and the product of the terminal gene in this island, CagA, is translocated into host epithelial cells by the cag secretion system following adherence (18–20). Within the host cell, CagA undergoes Src- and Abl-dependent tyrosine phosphorylation (21) and activates the eukaryotic phosphatase SHP-2, leading to dephosphorylation of host cell proteins and cellular morphological changes (19–21). CagA also dysregulates β-catenin signaling (22, 23) and apical-junctional complexes (24), events linked to increased cell motility and oncogenic transformation in several models (25, 26). Another H. pylori constituent linked to gastric cancer is the cytotoxin VacA, encoded by the gene vacA, which is present in virtually all H. pylori strains (27). In vitro, VacA induces the formation of intracellular vacuoles (27) and can induce apoptosis (28), and vacuolating activity is significantly associated with the presence of the cag pathogenicity island (3).Approximately 20% of H. pylori bind to gastric epithelial cells in vivo (29), and sequence analysis has revealed that the H. pylori genome contains an unusually high number of ORFs relative to its genome size that are predicted to encode outer membrane proteins (15). BabA, a member of a family of highly conserved outer membrane proteins and encoded by the strain-specific gene babA2, binds the Lewis^b histo-blood group antigen on gastric epithelial cells (30, 31), and H. pylori babA2⁺ strains are associated with an increased risk for gastric cancer (30). However, not all persons infected with cag⁺ babA2⁺ toxigenic strains develop gastric cancer, indicating that additional H. pylori constituents are important in carcinogenesis.We recently identified a strain of H. pylori, 7.13, that reproducibly induces gastric cancer in two rodent models of gastritis, Mongolian gerbils and hypergastrinemic INS-GAS mice (22). This strain was derived via in vivo adaptation of a clinical H. pylori strain, B128, which induces inflammation, but not cancer, in rodent gastric mucosa. The oncogenic 7.13 phenotype is not due to an enhanced ability of strain 7.13 to colonize as there were no significant differences in gastric colonization density or efficiency between strains B128 and 7.13 as assessed by either quantitative culture or histology. However, carcinogenic strain 7.13 binds more avidly to gastric epithelial cells in vitro than does strain B128, suggesting that the two strains may variably express different outer membrane proteins.To define proteins that may mediate the development of H. pylori-induced gastric cancer, we performed two-dimensional (2D)¹ DIGE coupled with MS to identify differentially abundant membrane-associated and cytosolic proteins from non-carcinogenic H. pylori strain B128 and its carcinogenic derivative, strain 7.13 (22). DIGE/MS is a well established proteomics technology based on conventional 2D gel protein separations whereby prelabeling samples with spectrally resolvable fluorescent dyes and multiplexing samples onto a series of gels that contain a mixture of all experimental samples (internal standard) provide quantitative data on abundance changes for thousands of intact proteins from multiple experimental conditions, each measured in replicate for statistical confidence (32–36). Techniques including DIGE/MS have recently been utilized to robustly define differences in protein abundance profiles between bacterial strains and to compare expression patterns of proteins harvested from bacteria maintained under different growth conditions (37, 38).Utilizing DIGE/MS, we detected and identified 26 proteins with statistically significant differences between strains B128 and 7.13, including a novel cysteine-to-arginine mutation in the H. pylori flagellar protein FlaA. We demonstrate that this FlaA mutation results in structural and functional aberrations. Application of this technique to two genetically related bacterial strains that induce distinct phenotypes also identified several novel candidate H. pylori virulence factors, providing a framework for studies targeting the pathogenesis of microbially induced cancer.     

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]     

A Sonic Hedgehog Missense Mutation Associated with Holoprosencephaly Causes Defective Binding to GAS1

Holoprosencephaly (HPE) is a common birth defect predominantly affecting the forebrain and face and has been linked to mutations in the sonic hedgehog (SHH) gene. HPE is genetically heterogeneous, and clinical presentation represents a spectrum of phenotypes. We have previously shown that Gas1 encodes a cell-autonomous Hedgehog signaling enhancer. Combining cell surface binding, in vitro activity, and explant culture assays, we provide evidence that SHH contains a previously unknown unique binding surface for its interaction with GAS1 and that this surface is also important for maximal signaling activity. Within this surface, the Asn-115 residue of human SHH has been documented to associate with HPE when mutated to lysine (N115K). We provide evidence that HPE associated with this mutation can be mechanistically explained by a severely reduced binding of SHH to GAS1, and we predict a similar result if a mutation were to occur at Tyr-80. Our data should encourage future searches for mutations in GAS1 as possible modifiers contributing to the wide spectrum of HPE.Holoprosencephaly (HPE)² is a developmental defect of the brain and face estimated to affect 1 in 250 conceptuses (1). Clinical presentation represents a spectrum of phenotypes, ranging from the most severe (alobar), where embryos have cyclopia and the prosencephalon fails to divide into hemispheres, to relatively mild defects (microform HPE) such as maxillary central incisor fusion, midfacial hypoplasia and clefting, and the presence of a single nostril (2). The use of mice as a model has proven invaluable for investigating the molecular and genetic causes of HPE. We have previously reported that microform HPE develops in growth arrest-specific gene 1 (Gas1) mutant mice (3, 4). Additionally, we determined that the 37-kDa, cell surface-presented, and glycosylphosphatidylinositol-anchored GAS1 protein binds to the secreted cell-cell signaling protein Sonic hedgehog (SHH) and that it functions as a cell-autonomous enhancer of SHH signaling activity (3, 5, 6). Consistently, the Gas1 mutant phenotype is more severe when an allele of Shh is removed, supporting a genetic interaction between the two genes (3, 4). Given the strong evidence that mutations in Shh can cause HPE in mice and humans (7–11), we investigated the hypothesis that some of these mutations cause defective SHH signaling due to a failed interaction with GAS1.Here we identify specific residues on SHH that are required for maximal binding to GAS1 and show, in both cell culture and explant culture assays, that these mutant SHH proteins have decreased signaling activity due to their defective interaction with GAS1. Significantly, one of these mutations has been associated with autosomal dominant HPE in a human family (9). These results lead us to propose that human embryos carrying this mutation may develop HPE due to a failed GAS1-SHH protein interaction.     

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]     

Deletion of the Chloride Transporter Slc26a7 Causes Distal Renal Tubular Acidosis and Impairs Gastric Acid Secretion

SLC26A7 (human)/Slc26a7 (mouse) is a recently identified chloride-base exchanger and/or chloride transporter that is expressed on the basolateral membrane of acid-secreting cells in the renal outer medullary collecting duct (OMCD) and in gastric parietal cells. Here, we show that mice with genetic deletion of Slc26a7 expression develop distal renal tubular acidosis, as manifested by metabolic acidosis and alkaline urine pH. In the kidney, basolateral Cl⁻/HCO₃⁻ exchange activity in acid-secreting intercalated cells in the OMCD was significantly decreased in hypertonic medium (a normal milieu for the medulla) but was reduced only mildly in isotonic medium. Changing from a hypertonic to isotonic medium (relative hypotonicity) decreased the membrane abundance of Slc26a7 in kidney cells in vivo and in vitro. In the stomach, stimulated acid secretion was significantly impaired in isolated gastric mucosa and in the intact organ. We propose that SLC26A7 dysfunction should be investigated as a potential cause of unexplained distal renal tubular acidosis or decreased gastric acid secretion in humans.The collecting duct segment of the distal kidney nephron plays a major role in systemic acid base homeostasis by acid secretion and bicarbonate absorption. The acid secretion occurs via H⁺-ATPase and H-K-ATPase into the lumen and bicarbonate is absorbed via basolateral Cl⁻/HCO₃⁻ exchangers (1–4). The tubules, which are located within the outer medullary region of the kidney collecting duct (OMCD),² have the highest rate of acid secretion among the distal tubule segments and are therefore essential to the maintenance of acid base balance (2).The gastric parietal cell is the site of generation of acid and bicarbonate through the action of cytosolic carbonic anhydrase II (5, 6). The intracellular acid is secreted into the lumen via gastric H-K-ATPase, which works in conjunction with a chloride channel and a K⁺ recycling pathway (7–10). The intracellular bicarbonate is transported to the blood via basolateral Cl⁻/HCO₃⁻ exchangers (11–14).SLC26 (human)/Slc26 (mouse) isoforms are members of a conserved family of anion transporters that display tissue-specific patterns of expression in epithelial cells (15–24). Several SLC26 members can function as chloride/bicarbonate exchangers. These include SLC26A3 (DRA), SLC26A4 (pendrin), SLC26A6 (PAT1 or CFEX), SLC26A7, and SLC26A9 (25–31). SLC26A7 and SLC26A9 can also function as chloride channels (32–34).SLC26A7/Slc26a7 is predominantly expressed in the kidney and stomach (28, 29). In the kidney, Slc26a7 co-localizes with AE1, a well-known Cl⁻/HCO₃⁻ exchanger, on the basolateral membrane of (acid-secreting) A-intercalated cells in OMCD cells (29, 35, 36) (supplemental Fig. 1). In the stomach, Slc26a7 co-localizes with AE2, a major Cl⁻/HCO₃⁻ exchanger, on the basolateral membrane of acid secreting parietal cells (28). To address the physiological function of Slc26a7 in the intact mouse, we have generated Slc26a7 ko mice. We report here that Slc26a7 ko mice exhibit distal renal tubular acidosis and impaired gastric acidification in the absence of morphological abnormalities in kidney or stomach.     

A Combined Proteomics and Metabolomics Profiling of Gastric Cardia Cancer Reveals Characteristic Dysregulations in Glucose Metabolism

Gastric cardia cancer (GCC), which occurs at the gastric-esophageal boundary, is one of the most malignant tumors. Despite its high mortality and morbidity, the molecular mechanism of initiation and progression of this disease is largely unknown. In this study, using proteomics and metabolomics approaches, we found that the level of several enzymes and their related metabolic intermediates involved in glucose metabolism were deregulated in GCC. Among these enzymes, two subunits controlling pyruvic acid efflux, lactate dehydrogenase A (LDHA) and pyruvate dehydrogenase B (PDHB), were further analyzed in vitro. Either down-regulation of LDH subunit LDHA or overexpression of PDH subunit PDHB could force pyruvic acid into the Krebs cycle rather than the glycolysis process in AGS gastric cancer cells, which inhibited cell growth and cell migration. Our results reflect an important glucose metabolic signature, especially the dysregulation of pyruvic acid efflux in the development of GCC. Forced transition from glycolysis to the Krebs cycle had an inhibitory effect on GCC progression, providing potential therapeutic targets for this disease.Gastric cardia cancer (GCC),1 which occurs at the gastric-esophageal boundary, is one of the most malignant tumors. Despite the steadily falling incidence of gastric non-cardia cancer in the past two decades (1), the rate of GCC has risen rapidly, establishing gastric cancer as the second major cause of cancer-related deaths throughout the world (2). GCC has become a significant cause of mortality and morbidity both in the west (3–5) and in Asia (6, 7), especially in China (8). Although this cancer has become an important health problem worldwide, the its pathogenesis has not been well characterized (1). Most patients are diagnosed at an advanced stage, contributing to the high mortality rate of the disease.Systematic proteomics analysis has proved to be a powerful approach in a variety of human cancer research, including lung (9), esophagus (10), gastric (11), liver (12), breast (13), and brain cancer (14). Metabolomics, another new bio-omics technology recently introduced into cancer research (15), is the global analysis of the small metabolites produced by normal or pathologic cellular processes. Some metabolic intermediates have been identified as new cancer biomarkers (16).Using proteomics and metabolomics methods in this study, we found that a series of proteins and metabolic intermediates, mainly involved in glucose metabolism, were altered during the development of GCC. The high activity of anaerobic glycolysis and the impairment of aerobic respiration occurring in these cells recapitulated the Warburg effect (17). Further studies using a gastric cancer cell line demonstrated that the predominant anaerobic glycolysis was essential for tumor cells to sustain rapid proliferation, whereas forced transition from anaerobic glycolysis to aerobic respiration inhibited the growth of tumor cells. In conclusion, our study revealed the major metabolic alterations essential for the development of GCC and discovered a biomarker signature of GCC. Such a finding has the potential to improve early diagnosis and prognosis and helps to identify new therapeutic targets.     

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.