Metabolic Flux Increases Glycoprotein Sialylation: Implications for Cell Adhesion and Cancer Metastasis*

This study reports a global glycoproteomic analysis of pancreatic cancer cells that describes how flux through the sialic acid biosynthetic pathway selectively modulates a subset of N-glycosylation sites found within cellular proteins. These results provide evidence that sialoglycoprotein patterns are not determined exclusively by the transcription of biosynthetic enzymes or the availability of N-glycan sequons; instead, bulk metabolic flux through the sialic acid pathway has a remarkable ability to increase the abundance of certain sialoglycoproteins while having a minimal impact on others. Specifically, of 82 glycoproteins identified through a mass spectrometry and bioinformatics approach, ∼31% showed no change in sialylation, ∼29% exhibited a modest increase, whereas ∼40% experienced an increase of greater than twofold. Increased sialylation of specific glycoproteins resulted in changes to the adhesive properties of SW1990 pancreatic cancer cells (e.g. increased CD44-mediated adhesion to selectins under physiological flow and enhanced integrin-mediated cell mobility on collagen and fibronectin). These results indicate that cancer cells can become more aggressively malignant by controlling the sialylation of proteins implicated in metastatic transformation via metabolic flux.

The surfaces of mammalian cells are covered with a dense layer of carbohydrates, collectively known as the glycocalyx, that influence many aspects of the interaction between a cell and its microenvironment. To date, the biosynthesis of cell surface displayed glycans has been thought to be controlled largely by individual glycosyltransferases based on the assumption that flux through the metabolic pathways that supply activated nucleotide sugar donors (the substrates for these enzymes) is not a limiting factor. For example, this premise has been used in mathematical models of sialylation (1), where concentrations of CMP-Neu5Ac in the lumen of the Golgi were assumed to be much higher than the K m of sialyltransferases (2). In the past several years, the idea that nucleotide sugars, exemplified by CMP-Neu5Ac (shown in Fig.  1A), are not a limiting or controlling factor in glycosylation has garnered one major and unambiguous exception, notably that changes in flux through the hexosamine biosynthetic pathway (HBP) 1 can alter UDP-GlcNAc levels with profound consequences on the branching of N-linked glycans (3). By changing the valence of these glycoconjugates, flux through the HBP can alter the galectin lattice and affect a host of downstream biological events including cancer progression (4); cell differentiation and proliferation (3); and autoimmunity, metabolic syndromes, and aging (5).
In this report, we demonstrate that the HBP is not unique in its ability to control surface glycoproteins via bulk metabolic flux. In particular, counter to earlier assumptions that flux through the sialic acid pathway does not significantly alter the sialylation of individual glycans (2), analysis of two "high-demand" sialoglycans (i.e. polysialylated NCAM (6) and podocalyxin (7)) suggested that fluctuations in the intracellular concentrations of sialic acid and the corresponding supply of CMP-Neu5Ac critically affected their production. In the current report, we used a global cell level approach to investigate whether these two examples were outliers or whether metabolic flux controls the surface display of sialic acid with fine resolution across a wide range of N-linked glycoproteins. We found that the sialylation of certain N-linked glycoproteins increased dramatically (e.g. by five-to eightfold) when flux through the sialic acid path-way was enhanced by exogenously supplied substrate whereas there was negligible effect on other glycoproteins. Importantly, these changes altered the adhesive behavior of cancer cells in a manner consistent with the glycoproteomic analysis. Together, these findings expand the role of metabolic flux in controlling glycosylation beyond the HBP and provide a foundation to explore the hypothesis that cancer cells modulate their metastatic potential and malignant progression via changes to bulk metabolic flux through the sialic acid biosynthetic pathway. Cell Culture-The human pancreatic carcinoma cell line SW1990 was obtained from the American Type Culture Collection (Manassas, VA) and cultured in Dulbecco's modified Eagle's medium supplemented with 10% fetal calf serum and 1.0% (v/v) of a 100 ϫ stock solution of penicillin and streptomycin (Invitrogen, Carlsbad, CA). Prior to cell lysis, pancreatic cells were detached from culture flasks using Enzyme Free Cell Dissociation Media (20 min at 37°C; Chemicon, Phillipsburg, NJ). CHO cells stably transfected with full-length E-selectin (CHO-E) or with phosphatidylinositol glycan-linked extracellular domain of P-selectin (CHO-P) were kindly donated by Affymax (Palo Alto, CA) and cultured as described previously (8). Synthesis of 1,3,4-O-Bu 3 ManNAc followed a previously reported procedure (9).
Periodate Resorcinal Assay-Total and bound sialic acid were quantified by the periodate resorcinol method. Briefly, treated and nontreated SW1990 cells were harvested with Enzyme Free Cell Dissociation Media and counted. The cells (2.25 ϫ 10 6 ) were washed twice with 1.0 ml D-PBS (Dulbecco's phosphate-buffered saline), pelleted, resuspended in 200 l D-PBS and divided into two 100 l aliquots. Cell lysates were obtained by four freeze/thaw cycles; the lysates were then analyzed by the method of Jourdian and coworkers (10) adapted for a 96-well plate format (11). Oxidation of aliquots of each sample was performed in parallel at 0°C and 37°C to measure total sialic acid or glycoconjugate-bound sialic acid (note that CMPsialic acid, although soluble, is detected in the what has conventionally been termed the "glycoconjugate-bound" fraction in this assay), respectively. Sialic acid concentrations were calculated by comparison with a standard curve (0 -250 M sialic acid) and expressed as the fold change from the control.
Lectin and Antibody Analysis of SW1990 Cells-Surface sialic acid was analyzed by labeling with fluorescein-conjugated RCA, FITCconjugated MAA, and fluorescein-conjugated Sambucus nigra agglutinin (SNA) lectins. In brief, 500,000 SW1990 cells were harvested as described above, washed with D-PBS, and resuspended in 100 l D-PBS. Lectin (4.0 l at 0.1 mg/ml) was added and the cell-lectin suspension was incubated for 1.0 h. Excess lectin was removed by washes with 1.0 ml D-PBS and the cells were analyzed by flow cytometry as described previously (12). Monoclonal antibody binding to determine Le X , sLe X , and sLe A expression was performed by washing 500,000 cells with D-PBS and incubating them with antihuman CD15 FITC (eBioscience Inc., San Diego, CA), 10 g/ml of mAb PE-labeled HECA 452 (PharMingen, San Diego, CA) in 100 l D-PBS with 0.1% BSA, or anti-sialyl Lewis A (Millipore, Billerica, MA), respectively. Background fluorescence for each antibody was determined using isotype-matched Ig. The samples were measured in a flow cytometer (BD Biosciences, San Diego, CA) and the data were analyzed using the FlowJo software (Treestar, Inc., San Carlos, CA).
Western Analysis of CD44 -Proteins (12 g) were resolved by SDS-PAGE and transferred electrophoretically onto a nitrocellulose membrane. The membrane was blocked with 5% nonfat milk/0.1% TBS-Tween 20 at room temperature for 2.0 h and then probed with anti-human CD44 antibody (515) at 1:1000 at 4°C overnight, followed by three washes with 0.1% TBS-Tween 20. HRP-conjugated antimouse IgG antibody was added at 1:2000 and incubated at room temperature for 1.0 h, followed by three washes with 0.1% TBS-Tween 20. The signal was visualized using SuperSignal Substrate (Pierce, Rockford, IL).
Mass Spectroscopy and Glyocopeptide Identification-Replicates of SW1990 cells (9 ϫ 10 6 ) with or without treatment with 1,3,4-O-Bu 3 ManNAc were harvested and lysed using radioimmunoprecipitation assay buffer (Sigma-Aldrich) according to the manufacturer's instructions. The protein concentration was determined using the BCA assay (Pierce); 1.4 mg of protein was used in following experiment for each cell treatment condition. Urea (8.0 M) and tris(2-carboxyethyl) phosphine hydrochloride (5.0 mM) were added to each sample to the indicated final concentrations followed by incubation at 60°C for 2.0 h. Iodoacetamide (10 mM final concentration) was added to each sample and they were incubated for 30 min at RT in the dark. Each solution was diluted eightfold with 100 mM KH 2 PO 4 (pH 8.0). Trypsin (Promega, 30 g) was added to each sample followed by incubation at 37°C overnight with shaking. Silver staining was used to monitor the completeness of trypsin digestion as indicated by the disappearance of high MW protein bands and the appearance of lower MW peptide bands (Ͻ10 kDa). After digestion, the samples were centrifuged at 16,110 ϫ g for 5.0 min to remove particulate matter.
The peptides were cleaned by elution through a C18 column and divided into two equal parts, N-linked total glycopeptides were then isolated from one sample using the solid phase extraction of glycopeptides method (13)(14)(15). This method was modified to extract Nlinked sialoglycopeptides from the second aliquots of each sample. Briefly, the modified solid phase extraction of glycopeptidesmethod used 1.0 mM precooled sodium periodate at 0°C for 15 min to oxidize sialic acid selectively. The resulting sialoglycopeptides were covalently conjugated to a solid support via hydrazide chemistry whereas nonsialylated glycopeptides were removed by washing prior to release of N-linked glycopeptides from solid support by PNGase F. The resulting sialoglycopeptides were concentrated by elution through C18 columns, dried, and resuspended in 20 l of 0.4% acetic acid.
The peptides (5.0 l/sample) were labeled with isobaric tag for relative and absolute quantitation (iTRAQ) 8plex (AB SCIEX) according to the manufacturer's instructions. The formerly N-linked total glycopeptides from control cells were labeled in duplicate by iTRAQ with 113 and 114 whereas those from 1,3,4-O-Bu 3 ManNAc-treated cells were labeled with 115 and 116. The formerly N-linked sialoglycopeptides from control cells were labeled in duplicate by iTRAQ 117 and 118, whereas those from1,3,4-O-Bu 3 ManNAc-treated cells were labeled in duplicate with 119 and 121 separately. The labeled peptides were then mixed, cleaned using strong cation exchange columns, and resuspended in 15 l of 0.4% acetic acid.
iTRAQ-labeled peptides (6 l) were analyzed by liquid chromatography-tandem MS (LC-MS/MS) using an LTQ-Orbitrap velos (Thermo-Fisher, Waltham, MA) coupled with a 15 cm ϫ 75 m C18 column (5 m particles with 100 angstrom pore size). A nanoAquity UPLC at 300 nl/min with a 90-min linear acetonitrile gradient (from 5-32% B over 90 min; A ϭ 0.1% formic acid in water, B ϭ 0.1% formic acid in acetonitrile) was used. Top 10 data dependent with exclusion for 20 s was set. The samples were run with HCD fragmentation at a normalized collision energy of 45 and an isolation width of 1.2 Da. Monoisotopic precursor selection was enabled and the dynamic exclusion was set to 30 s with a repeat count of 1 and Ϯ 10 ppm mass window. The source voltage was 2.0 kDa. A lock mass of the polysiloxane peak at 371.10123 was used to correct the mass in MS and MS/MS. Target values in MS were 1e6 ions at a resolution setting of 30,000 and in MS2 1e5 ions at a resolution setting of 7500.
MS/MS spectra were searched with MASCOT (version 2.2.0) using Proteome Discoverer (version 1.0) (Thermo Fisher) against human subdatabase of NCBI Reference Sequence (RefSeq) (version 40, released at April 16, 2010) containing 29,704 sequences. Integration window tolerance was set at 20 ppm for peak integration. The peptide cutoff score was set at 30. For this database search, the precursor mass tolerance and fragment mass tolerance was set at 15 ppm and 0.05 Da respectively, fixed modifications were set as iTRAQ 8plex labeling (N-term and K) and carbamidomethylation (C), and other database-searching parameters were set as flexible modification as follows: deamidation (NQ) and oxidation (M). Semi-tryptic end and 1 missed cleavage site was allowed. The False Discovery Rate was set at 0.01 to eliminate low-probability protein identifications. The data generated by mass spectrometry for the global N-glycopeptides and N-sialoglycopeptides analysis may be downloaded from Proteome-Commons.org Tranche using the following hash: dcLmbIAQlBZvBYIgiJjGG5om/IK/fq7eϩ5TpB3TAg7PcWhk6CJzg-TuePjUH7HJktURuizj6JD7Ack9TJDkqr8y9WGGQAAAAAAAACtgϭϭ. The hash may be used to prove exactly what files were published as part of this manuscript's data set, and the hash may also be used to check that the data has not changed since publication. For single peptide identification, the matched spectra can be found in supplemental Table S2.
Kyoto Encyclopedia of Genes and Genomes Pathway Analysis of Identified Glycoproteins-The DAVID Bioinformatics Resources 6.7 tool (National Institute of Allergy and Infectious Diseases, NIAID, NIH) was used to analyze the glycoproteins and sialoglycoproteins identified in the present study and assign them into categories compiled by the KEGG database. All analysis parameters were used in the default setting.
Evaluation of the Cut-off of the Protein Abundance Ratio for Protein Changes-To correct for any systematic errors of protein ratio introduced by sample handling and to determine the appropriate cutoff for protein changes, the distribution of abundance ratios of glycosylation changes was generated. Because the majority of proteins were not expressed differently in two cell states, we normalized the ratio based on the distribution of the protein abundance ratios from two cell states. Proteins that fell outside the normal distribution from the abundance ratio of two cell states were considered as altered proteins. The threshold to select protein changes was based on the ratio distribution of two cell states. The mean and standard deviation of the ratio from the two cell states were calculated, and the abundance of proteins with an abundance ratio outside of one standard deviation from the mean were flagged as altered.
Flow-based Adhesion Assays-To simulate the physiological shear environment of the vasculature, pancreatic carcinoma cells in D-PBS containing Ca ϩ2 /Mg ϩ2 , 0.1% bovine serum albumin were perfused over immobilized E-, or L-selectin-coated dishes at prescribed wall shear stresses using a parallel plate flow chamber (250-m channel depth, 5.0-mm channel width). Adhesion was quantified by perfusing cells at 1.0 ϫ 10 6 /ml and recording for 2 or 5 min. Average rolling velocities were computed as the distance traveled by the centroid of the translating cell divided by the time interval at the given wall shear stress (16).
Blot Rolling Assays-SW1990 whole cell lysate was prepared by membrane disruption using 2.0% Nonidet P-40 followed by differential centrifugation. Blots of immunopurified CD44 from treated and untreated SW1990 whole cell lysate were stained with anti-CD44 (2C5) or anti-CD44 (515) or HECA-452 mAbs and rendered translucent by immersion in 90% D-PBS 10% glycerol. The blots were placed under a parallel plate flow chamber CHO transfectants expressing P-or E-selectin, resuspended at 1.0 ϫ 10 6 cells/ml in 90% D-PBS/10% glycerol, and perfused at the shear stress of 0.5 dynes/ cm 2 (16). Molecular weight markers were used as guides to aid placement of the flow chamber over stained bands of interest. The number of interacting cells per lane was averaged over five 10ϫ fields of view (0.55 mm 2 each) within each stained region. Nonspecific adhesion was assessed by perfusing 5.0 mM EDTA in the flow medium.
Surface Plasmon Resonance Binding Assay-Surface plasmon resonance binding studies were performed using a BIAcore 3000 system (BIAcore Inc., Piscataway, NJ). Recombinant human E-, L-, and P-selectin-Fc chimera (500 g/ml in 50 mM sodium acetate, pH 4.75; R&D Systems) were covalently immobilized via amine coupling to a sensor chip (CM5) to channels 2, 3, and 4 respectively as directed by the manufacturer (BIAcore). IgG was immobilized to Channel 1 and used as control. To protect the binding sites, CaCl 2 , MgCl 2 , fucose, and galactose were added to the immobilization buffer at a final concentration of 1.0 mM each. The immobilization of the selectins resulted in an average of 2500 resonance units. The binding assays were performed in PBS, with 1.0 mM CaCl 2 and 1.0 mM MgCl 2 , pH 7.4 at 25°C with a flow rate of 10 l/min. The interactions of the same concentration of immunoprecipitated CD44 from 1,3,4-O-Bu 3 ManNAc treated and untreated control SW1990 cells were monitored and the signal response was recorded and corrected for nonspecific binding to the control channel. Data was analyzed using BIAevaluation software (V4.1, BIAcore).
Wound Healing Assays-Culture plates were coated with 20 g/ml collagen-1 or 20 g/ml fibronectin (BD Biosciences, Bedford, MA) in D-PBS. Cells were plated and 4.0 l from a 50 mM stock solution of 1,3,4-O-Bu 3 ManNAc was added to treated samples (to give a final analog concentration of 100 M) and 4.0 l of 100% ethanol was added to the controls. Cells were incubated for 24 h at 37°C to allow the formation of a confluent monolayer that was scratched by carefully dragging a p200 pipette tip across the cells. Fresh medium without FBS, and with or without 1,3,4-O-Bu 3 ManNAc was added, and the cells were incubated at 37°C on a live-cell microscope where pictures of wound closure were taken every 20 min for 24 h. The wound area was measured using the Nikon Imaging Software (NIS) for each time point.
Modeled Representations of Identified Glycopeptides in the Overall Structures of CD44 and Integrin ␣6-For CD44, the identified glycan was modeled onto the NMR structure of CD44 (PDB ID:1UU), using the program GlyProt (http://www.glycosciences.de/glyprot/) (17). Because a corresponding structure is not available for integrin ␣6, the structure of the integrin ␣V (PBD ID: 3ije) was used to model the subunit repeat of the propeller where the identified glycopeptide was located using a comparative protein structural modeling approach (18). The GlycoProt program was used to N-glycosylate the model subunit with a simplified N-glycan structure. Visualization of structures was performed with the Chimera software (19).
Statistical Analysis-Data are expressed as the mean S.E. for at least three independent experiments. Statistical significance of differences between means was determined by analysis of variance. RESULTS

1,3,4-O-Bu 3 ManNAc Increased Intracellular and Surface
Sialic Acid-Pancreatic cancer SW1990 cells incubated with 1,3,4-O-Bu 3 ManNAc experienced an increase in total and glycoconjugate-bound sialic acid ( Fig. 2A) consistent with previous results that cells incubated with ManNAc analogs have greatly increased levels of intracellular sialic acid and more modestly elevated surface display of this monosaccharide (12,20,21). Lectin binding analysis (Fig. 2B) confirmed that the surface display of sialic acid increased insofar as the number of "empty" penultimate galactose/GalNAc residues decreased (as measured by RCA binding) whereas a corresponding increase in the surface display of ␣2,6-linked (SNA) and ␣2,3-linked (MAA) sialic acids occurred. Similarly, antibody binding assays (Fig. 2C) showed that Lewis X (Le X ) levels decreased in analog-treated cells whereas sialyl Lewis X (sLe X ) and sialyl Lewis A (sLe A ) levels increased. These results demonstrate that the overall level of glycans on the cell surface was not changing; instead the fraction of sialylated Nglycans was altered by analog-driven flux through the sialic acid pathway. It is noteworthy that the analysis of nonsialylated epitopes (e.g. CD15 (Fig. 2G)) indicated that analogenhanced flux through the sialic pathway did not simply saturate all possible sites of sialic acid attachment because if it did, CD15 levels would be reduced to the levels observed in the isotype control. Once we established that bulk flux through the sialic acid pathway in SW1990 cells altered surface sialic acid display, we used mass spectrometry-based methods to identify the specific sites of glycosylation associated with these changes and then investigated whether these changes resulted in altered cell behavior, as described below.
Glycoproteomic Analysis of Metabolic Flux-driven Changes to N-Glycan Sialylation-Glycoproteomic methods were next used to perform a global analysis of sialylation changes that occurred in 1,3,4-O-Bu 3 ManNAc treated SW1990 cells. More specifically, to identify individual proteins whose sialylation status was impacted by analog-driven flux, the formerly Nlinked total glycopeptides and sialoglycopeptides were isolated from cells that had been pretreated or not with 1,3,4-O-Bu 3 ManNAc, labeled with iTRAQ reagents, and analyzed by mass spectrometry (13)(14)(15)22). Two proteins, CD44 and integrin ␣6, were used as examples to illustrate this process, as shown in Fig. 3. When conducted on a global scale to identify cell wide changes with a cutoff of a 1% false discovery rate, 105 unique glycopeptides (containing the NXS/T motif) were identified representing 82 unique glycoproteins; the ratio of the peptide level and protein level in each sample obtained from analog-treated cells to control cells was calculated and is given in supplemental Table S1. An important finding was that the relative levels of peptides isolated with or without pretreatment with 1,3,4-O-Bu 3 ManNAc was approximately the same when they were captured via their N-linked glycans. By contrast, the levels of certain peptides isolated after pretreatment with 1,3,4-O-Bu 3 ManNAc varied significantly from control samples when captured via the sialic acid specific strategy (Table I). Together, these results verified that this analog gave rise to sialic acid specific changes while leaving N-glycan levels unchanged.
Several noteworthy trends observed from the glycoproteomic analyses are summarized in Fig. 4. First, histograms were used to compare the number of peptides with different abundance ratios. Because of errors introduced by analytical procedures, such as sample handling and quantification process, the ratios may be shifted slightly, therefore glycopeptides that were distributed within one S.D. (0.34) of the mean (1.07) were considered to be unchanged whereas those that fell beyond one S.D. of the normal distribution curve (Ͻ0.73 or Ͼ1.41) were considered to be changed. Based on these criteria, the ratio distribution of total N-linked glycoproteins identified from 1,3,4-O-Bu 3 ManNAc treated cells compared with untreated controls fit the curve of a normal distribution that almost exclusively fell within the "unchanged" group (i.e. between 0.73 and 1.41, Fig. 4A). By contrast, the ratio distribu-tion of sialoglycopeptides was shifted to the right after treatment of 1,3,4-O-Bu 3 ManNAc indicating increased sialylation (Fig. 4A). Because not all sialoglycopeptides exhibited an increase and the ones that did had a distinct biomodal distribution, we divided these samples into three groups for further analysis: Group 1, statistically unchanged sialylation; Group 2, modestly increased sialylation; and Group 3, substantially increased sialylation.
Pathway Analysis Revealed that Group 2 "Modestly Increased" Sialoglycoproteins Include Adhesion Molecules (e.g. CD44 and Integrin ␣6) Implicated in Metastasis-David pathway analysis was used to assign the glycopeptides identified from the mass spectrometry data to pathways delineated by the KEGG database. This analysis revealed that the total sets of proteins identified from both of the N-glycan and sialic acidspecific capture methods fit into 13 pathway classes that were distributed very similar to each other (Fig. 4B). The comparison of these two data sets revealed that each profile consisted of 12 categories (listed in Fig. 4C) with the main difference being a complete absence of proteins in the hypertrophic cardiomyopathy class (xiii) in the N-linked set and an absence of hypertrophic myopathy (ix) proteins in the sialylated set.
A more detailed analysis of the three groups of sialylated glycoproteins (i.e. Group 1, 2, and 3 from Fig. 4A) revealed that Groups 1 and 3 had a fairly broad distribution of pathway categories (Fig. 4D). By contrast, Group 2 proteins were restricted to 4 of the 13 categories with two of these categories (v, cell adhesion molecules; and xii, renin angiogensis system) completely absent from the unchanged (Group 1) and highly increased (Group 3) sialoglycoproteins.

Flux-driven Sialylation of CD44 and Integrin ␣6 Enhances
Adhesion Associated with Metastasis-We reasoned that the "Group 2" proteins might be restricted to a modest, but real, increase in sialylation because they were involved in critical cellular processes very sensitive to changes in this PTM. The identification of CD44 and integrin ␣6 in this category (Table I) provided an opportunity to test whether any of the sialic acid-specific changes identified in the glycoproteomic data had a meaningful impact on cell behavior because of these molecules involvement in cancer progression and metastasis (23)(24)(25). Metastasis is a multistep process in which cancerous cells separate from the primary tissue and enter the circulatory system where they interact with various host cells before they migrate and lodge in the target organ and form secondary metastatic colonies. Mounting evidence suggests a role for selectins with ligands such as CD44 in the initial interactions between host cells and tumor cells in vasculature and integrins such as integrin ␣6 in tumor cell migration through the extracellular matrix. Accordingly, we used several in vitro   Fig.  3 for CD44 and integrin ␣6). B, Pathway classes of proteins represented by both of the classes of glycopeptides (e.g. analog-treated and control) shown in (A); the classes i to xiii are listed in (C). Because virtually all of the proteins identified from the total N-linked glycans fell within one S.D. of the mean (i.e. between 0.73 and 1.41) in a bell shaped distribution, pathway analysis was performed only for the entire group. D, Additional pathway analysis, however, was conducted for the sialoglycoproteins based on the groups of sialopeptides shown in A: specifically, Group 1 included proteins that fell within one S.D. of the mean; Group 2 included proteins that experienced a modest but statistically significant increase in sialylation upon analog treatment; and Group 3 included proteins with sialylation increases of greater than ϳtwofold. a Amino acids indicated in lower case at both ends of the sequence represent the iTRAQ labeled N-termini and Lys, lower case n in the nXT/S motif represents the formerly glycosylated Asp and delaminated after PNGase F release, lower case m represents the oxidation on M.
b Pep Avg TG S/C specifies the ratio of total glycopeptides of sample/control quantified at peptide level. c Pep Avg SG S/C specifies the ratio of sialoglycopeptides of sample/control quantified at peptide level. d Pro Avg TG S/C specifies the ratio of total glycopeptides of sample/control quantified at protein level. e Pro Avg SG S/C specifies the ratio of sialoglycopeptides of sample/control quantified at protein level. and control SW1990 cells to immobilized selectins under physiological flow conditions. The analog-treated cells adhered strongly to the E-selectin coated surface, and consequently rolled significantly slower than control cells at 0.5 dynes/cm 2 (Fig. 5A). A smaller, but still statistically significant, decrease was noted for L-selectin-dependent rolling under the same shear conditions (Fig. 5B) consistent with prior work showing that CD44 possesses ligand activity for both of these selectins (27)(28)(29). Importantly, the site density of CD44 expressed in the analog-treated cells did not increase, as assessed by flow cytometry (Fig. 5C) and Western blot analysis (Fig. 5D) implicating the increased sialylation of these molecules as the sole factor contributing to the increased adhesion to selectins observed in the 1,3,4-O-Bu 3 ManNAc treated cells.
Although sialylated epitopes on CD44 are linked firmly to cancer progression, additional factors nonetheless could have contributed to the increased adhesion of 1,3,4-O-Bu 3 ManNAc treated cells to E-and L-selectin. Therefore, to gain additional evidence that the observed increase in sialylation of CD44 was causatively linked to changes in the adhesive behavior of analog-treated SW1990 cells, a blot rolling assay was performed in which selectin-expressing CHO cells were perfused over SDS-PAGE-resolved immunoprecipitated CD44 from 1,3,4-O-Bu 3 ManNAc treated or untreated SW1990 cells. As shown in Fig. 5E, immunopurified CD44 from analog-treated cells exhibited a markedly increased capacity to capture E-and P-selectin-transfected CHO cells under dynamic flow conditions compared with control cells; this binding was completely ablated by sialidase treatment. In line with the flow-based rolling assays, a smaller difference was observed with L-selectin-expressing lymphocytes (data not shown). To further verify that CD44 from 1,3,4-O-Bu 3 ManNAc-treated SW1990 cells interacted more efficiently with E-selectin relative to controls, we used surface plasmon resonance to show that E-selectin rapidly and avidly bound to CD44 isolated from analog-treated SW1990 cells (Fig. 5F). In contrast, only modestly increased affinity of CD44 isolated from 1,3,4-O-Bu 3 ManNAc treated cells was observed for Land P-selectins, indicating that sialylation of this glycoprotein preferentially enhanced adhesion to E-selectin. Taken together, this data indicates that changes in CD44 sialylation and not transcriptional control or other impacts to glycosylation were primarily responsible for the enhanced adhesion of analog-treated pancreatic cancer cells to E-selectin.
Moving to integrins, in general there have been several examples reported where hypersialylation of these cell adhesion molecules contributes to cancer progression by increasing cell motility (24,25) through the extracellular matrix ECM after selectin-mediated extravasation from the vascu- lature. More specifically, however, although integrin ␣6 has long been known to play a role in the motility of pancreatic cancer cells (30), the precise involvement of sialic acids in this process remains unclear. To address this issue and test whether the increased sialylation of integrin ␣6 observed in 1,3,4-O-Bu 3 ManNAc treated cells influenced cell migration through the ECM, we used a wound-healing assay where the mobility of SW1990 cells on ECM substrates known to be integrin ligands was monitored. As shown in Fig. 6, the migration of analog-treated cells increased slightly on collagen type I (Fig. 6A), strongly on fibronectin (Fig. 6B), whereas no difference was seen on a control, BSA-coated surface (Fig. 6C). The function blocking GoH3 antibody to integrin ␣6 significantly impeded cell migration of 1,3,4-O-Bu 3 ManNAc-treated cells verifying that this integrin played a key role in the increasing the motility of the analog-treated cells. Finally, flow cytometry analysis with GoH3 showed that there was no difference in the protein levels of integrin ␣6 in the analog-treated cells (Fig. 6D), thereby verifying that the observed changes in adhesion were specific to changes in sialic acid display. DISCUSSION In the past, two methods have been used to alter metabolic flux through the sialic acid pathway (Fig. 1A) and both have had limitations. One has been the manipulation of GNE, an enzyme that regulates flux into the sialic acid pathway (31,32). GNE, however, is a multifunctional protein that modulates the transcription of enzymes involved in glycosylation (33) and it also interacts directly with cell adhesion molecules (34) in multiple ways that could affect the endpoints evaluated in this study. Alternately, the sialic acid pathway can be supplemented with exogenous ManNAc but the high levels needed (often in the tens of millimolar (11)) are problematic whereas more efficiently utilized peracetylated analogs (35) have significant and often deleterious side effects (36). The current study avoided these problems by employing the "high flux" tributanoylated analog 1,3,4-O-Bu 3 ManNAc that allows intracellular sialic acid levels to be elevated to high levels (9) with negligible cytotoxicity or perturbation to gene regulation (37,38). In this study, treatment of SW1990 cells with 1,3,4-O-Bu 3 ManNAc increased total sialic acid levels several fold and impacted glycoconjugate-bound and surface sialylation more modestly (e.g. by Յ twofold) when measured using lectins or antibodies at cell-level resolution (Fig. 2).
Having confirmed that globally increased surface sialylation occurred in 1,3,4-O-Bu 3 ManNAc-treated SW1990 cells, we conducted a glycoproteomic analysis (Fig. 3) that revealed important insights into the global cellular impact of bulk flux through the sialic acid pathway. First, flux driven changes to intracellular sialic acid levels quite remarkably did not uniformly increase the sialylation of all surface glycans but instead selectively tuned the sialylation status of individual glycoproteins ( Fig. 4A) with ϳ31% of these molecules remaining FIG. 6. Wound healing, cell migration assays. SW1990 cells were treated (or not) with 100 M 1,3,4-O-Bu 3 ManNAc (Cpd 1) and plated on culture plates pretreated with 20 g/ml (A) collagen I, (B) fibronectin, or (C) BSA. After 24 h, wounds were created and new media without FBS but with (or without) fresh 1 was added. For each condition, 10 g/ml of the integrin ␣6 functional blocking antibody (GoH3) was added to adjacent wells. The mobility of analog-treated cells was evaluated in comparison with control cells that had not been treated with analog by measuring the accumulated distance traveled per hour with the data given as a fold increase of analog-treated cells compared with control samples. D, Verification that protein levels of integrin ␣6 did not change was obtained by flow cytometry. statistically unchanged in 1,3,4,-O-Bu 3 ManNAc-treated cells while others (ϳ40%) experienced an increase of 200 -800% or more. An even more interesting group was the subset of ϳ30% of the glycoproteins that experienced a small (e.g. ϳ1.5to twofold) but statistically significant increase in sialylation after treatment with 1,3,4-O-Bu 3 ManNAc. David pathway analysis revealed that these proteins were limited to four categories (v, cell adhesion molecules; vi, lysosome; xi, hematopoietic stem cell lineages, and xii, renin angiogensis system, Fig. 4D).
An intriguing explanation for the fairly restricted subset of "Group 2" glycoproteins ( Fig. 4D) that experienced tightly regulated increases in sialylation was that the activities of these proteins are unusually sensitive to sialylation. Accordingly, although cells are able to increase their sialylation thereby moving them out of "Group 1", tight regulation is required to prevent their activity from being modified too severely, thus preventing them from moving into "Group 3." Alternately, there may be no need to increase sialylation too greatly because even modest increases can modulate biolog- FIG. 7. Representation of the identified glycopeptides aFnSTLPTmAQmEk in CD44 and eINSLnLTESHnSR and anHSGAVVLk in integrin ␣6. A, A cartoon representation of CD44 is shown with blue lollipops illustrating the positions of putative N-linked glycans. Successive "zoomed in" depictions show the location of the aFnSTLPTmZQEk glycopeptide within a computationally generated surface illustration of the HA binding domain of CD44. B, A cartoon representation of the integrin ␣6␤1 complex (bottom) employs blue lollipops to illustrate the positions of putative N-linked glycans except for the red and green lollipops, which represent the actual N-glycans identified by mass spectrometry in this study. A modeled depiction of the integrin propeller subunit repeat containing the anHSGAVVLk glycopeptide is shown along with a zoomed in view of this glycopeptide with a representative N-glycan attached and shown using a sticks format. The green lollipop in the Calf-2 region represents the eINSLnLTESHnSR glycopeptide that also experienced selectively enhanced sialylation; sufficient structural information, however, is not available to further model this site. ical activities. We obtained experimental evidence for the latter premise by demonstrating that the sialylation status of CD44 and integrin ␣6, which are members of the cell adhesion category of tightly regulated proteins (i.e. "v" in the David pathway analysis, Fig. 4), contributed to changes in cancer cell mobility consistent with the known role of sialic acid in these processes. For example, increased sialylation of CD44 resulted in enhanced binding to selectins (Fig. 5); cancer cells exploit selectin-mediated adhesion to exit the vascular during metastasis. Another way sialylation assists metastasis is through integrin-mediated adhesion that facilitates cell migration through the ECM. We measured this end point via a wound healing assay wherein analog-treated SW1990 cells showed an increased ability to migrate across ECM components (Fig. 6). It is noteworthy that the increased mobility of analog-treated cells on ECM components complements enhanced selectin-mediated adhesion to facilitate two distinct facets of metastasis, namely the (1) tethering and rolling and (2) firm adhesion steps of the extravasation process.
Of the multiple potential sites of N-glycosylation for CD44 and integrin ␣6 that may have affected the activity of these molecules, we identified one site for CD44 and two for integrin ␣6 that were selectively hypersialated in cells treated with 1,3,4-O-Bu 3 ManNAc. Two of these three changes occurred in domains of the host proteins previously implicated in adhesion and provide the new information that sialic acid attached to a specific site of N-glycosylation is an important determinant of adhesion. These two sites are the hyaluronan binding domain of CD44 (Fig. 7A) (39) and the propeller domain of an integrin (Fig. 7B) (40). Although the impact of the increased sialylation observed for the third glycopeptide, located in the Calf-2 domain of integrin ␣6 (Fig. 7B) on adhesion is less clear, precedent that N-glycans present in the hinge region of an integrin can stabilize the open, activated form (41) suggests that allosteric effects of this third glycopeptide theoretically also could modulate integrin-mediated adhesion.
In conclusion, this paper utilized 1,3,4-O-Bu 3 ManNAc, a recently developed molecular tool for manipulating flux through the sialic acid pathway (9, 42) that avoids pitfalls of previous genetic and small molecule-based approaches, to demonstrate that intracellular sialic acid levels selectively tune the sialylation status of individual surface glycans. This finding conclusively demonstrated that metabolic flux can determine the display and biological activity of sialic acid, an important cell surface carbohydrate. This work also supports the hypothesis that metabolic flux can alter the metastatic potential of cancer cells in a glycan-dependent manner reminiscent of how the activity of N-acetylgalactosyltransferases MGAT4/5 and the subsequent branching of N-glycans depends on flux through the HBP (3,4,43), which is the current exemplar of how metabolic flux-driven changes can modulate cell surface glycosylation.