Folding and modulation of the helical conformation of Glycophorin A by point mutations.

Transmembrane helix folding and self-association play important roles in biological signaling and transportation pathways across biomembranes. With molecular simulations, studies to explore the structural biochemistry of this process have been limited to focusing on individual fragments of this process - either helix formation or dimerization. While at an atomistic resolution, it can be prohibitive to access long spatio-temporal scales, at the coarse grained (CG) level, current methods either employ additional constraints to prevent spontaneous unfolding or have a low resolution on sidechain beads that restricts the study of dimer disruption caused by mutations. To address these research gaps, in this work, we apply our recent, in-house developed CG model (ProMPT) to study the folding and dimerization of Glycophorin A (GpA) and its mutants in the presence of Dodecyl-phosphocholine (DPC) micelles. Our results first validate the two-stage model that folding and dimerization are independent events for transmembrane helices and found a positive correlation between helix folding and DPC-peptide contacts. The wild type (WT) GpA is observed to be a right-handed dimer with specific GxxxG contacts, which agrees with experimental findings. Specific point mutations reveal several features responsible for the structural stability of GpA. While the T87L mutant forms anti-parallel dimers due to an absence of T87 interhelical hydrogen bonds, a slight loss in helicity and a hinge-like feature at the GxxxG region develops for the G79L mutant. We note that the local changes in the hydrophobic environment, affected by the point mutation, contribute to the development of this helical bend. This work presents a holistic overview of the structural stability of GpA in a micellar environment, while taking secondary structural fluctuations into account. Moreover, it presents opportunities for applications of computationally efficient CG models to study conformational alterations of transmembrane proteins that have physiological relevance.

A Balance between Transmembrane-Mediated ER/Golgi Retention and Forward Trafficking Signals in Glycophorin-Anion Exchanger-1 Interaction.

Anion exchanger-1 (AE1) is the main erythroid Cl-/HCO3- transporter that supports CO2 transport. Glycophorin A (GPA), a component of the AE1 complexes, facilitates AE1 expression and anion transport, but Glycophorin B (GPB) does not. Here, we dissected the structural components of GPA/GPB involved in glycophorin-AE1 trafficking by comparing them with three GPB variants-GPBhead (lacking the transmembrane domain [TMD]), GPBtail (mainly the TMD), and GP.Mur (glycophorin B-A-B hybrid). GPB-derived GP.Mur bears an O-glycopeptide that encompasses the R18 epitope, which is present in GPA but not GPB. By flow cytometry, AE1 expression in the control erythrocytes increased with the GPA-R18 expression; GYP.Mur+/+ erythrocytes bearing both GP.Mur and GPA expressed more R18 epitopes and more AE1 proteins. In contrast, heterologously expressed GPBtail and GPB were predominantly localized in the Golgi apparatus of HEK-293 cells, whereas GBhead was diffuse throughout the cytosol, suggesting that glycophorin transmembrane encoded an ER/Golgi retention signal. AE1 coexpression could reduce the ER/Golgi retention of GPB, but not of GPBtail or GPBhead. Thus, there are forward-trafficking and transmembrane-driven ER/Golgi retention signals encoded in the glycophorin sequences. How the balance between these opposite trafficking signals could affect glycophorin sorting into AE1 complexes and influence erythroid anion transport remains to be explored.

Effects of cholesterol and PIP2 on interactions between glycophorin A and Band 3 in lipid bilayers.

In the erythrocyte membrane, the interactions between glycophorin A (GPA) and Band 3 are associated strongly with the biological function of the membrane and several blood disorders. In this work, using coarse-grained molecular-dynamics simulations, we systematically investigate the effects of cholesterol and phosphatidylinositol-4,5-bisphosphate (PIP2) on the interactions of GPA with Band 3 in the model erythrocyte membranes. We examine the dynamics of the interactions of GPA with Band 3 in different lipid bilayers on the microsecond time scale and calculate the binding free energy between GPA and Band 3. The results indicate that cholesterols thermodynamically favor the binding of GPA to Band 3 by increasing the thickness of the lipid bilayer and by producing an effective attraction between the proteins due to the depletion effect. Cholesterols also slow the kinetics of the binding of GPA to Band 3 by reducing the lateral mobility of the lipids and proteins and may influence the binding sites between the proteins. The anionic PIP2 lipids prefer binding to the surface of the proteins through electrostatic attraction between the PIP2 headgroup and the positively charged residues on the protein surface. Ions in the solvent facilitate PIP2 aggregation, which promotes the binding of GPA to Band 3.

Engineered red blood cells carrying PCSK9 inhibitors persistently lower LDL and prevent obesity.

Low plasma levels of Proprotein Convertase Subtilisin/Kexin 9 (PCSK9) are associated with decreased low-density lipoprotein (LDL) cholesterol and a reduced risk of cardiovascular disease. PCSK9 binds to the epidermal growth factor-like repeat A (EGFA) domain of LDL receptors (LDLR), very low-density lipoprotein receptors (VLDLR), apolipoprotein E receptor 2 (ApoER2), and lipoprotein receptor-related protein 1 (LRP1) and accelerates their degradation, thus acting as a key regulator of lipid metabolism. Antibody and RNAi-based PCSK9 inhibitor treatments lower cholesterol and prevent cardiovascular incidents in patients, but their high-cost hampers market penetration. We sought to develop a safe, long-term and one-time solution to treat hyperlipidemia. We created a cDNA encoding a chimeric protein in which the extracellular N- terminus of red blood cells (RBCs) specific glycophorin A was fused to the LDLR EGFA domain and introduced this gene into mouse bone marrow hematopoietic stem and progenitor cells (HSPCs). Following transplantation into irradiated mice, the animals produced RBCs with the EGFA domain (EGFA-GPA RBCs) displayed on their surface. These animals showed significantly reduced plasma PCSK9 (66.5% decrease) and reduced LDL levels (40% decrease) for as long as 12 months post-transplantation. Furthermore, the EGFA- GPA mice remained lean for life and maintained normal body weight under a high-fat diet. Hematopoietic stem cell gene therapy can generate red blood cells expressing an EGFA-glycophorin A chimeric protein as a practical and long-term strategy for treating chronic hyperlipidemia and obesity.

Addressing the Excessive Aggregation of Membrane Proteins in the MARTINI Model.

The MARTINI model is a widely used coarse-grained force field popular for its capacity to represent a diverse array of complex biomolecules. However, efforts to simulate increasingly realistic models of membranes, involving complex lipid mixtures and multiple proteins, suggest that membrane protein aggregates are overstabilized by the MARTINI v2.2 force field. In this study, we address this shortcoming of the MARTINI model. We determined the free energy of dimerization of four transmembrane protein systems using the nonpolarizable MARTINI model. Comparison with experimental FRET-based estimates of the dimerization free energy was used to quantify the significant overstabilization of each protein homodimer studied. To improve the agreement between simulation and experiment, a single uniform scaling factor, α, was used to enhance the protein-lipid Lennard-Jones interaction. A value of α = 1.04-1.045 was found to provide the best fit to the dimerization free energies for the proteins studied while maintaining the specificity of contacts at the dimer interface. To further validate the modified force field, we performed a multiprotein simulation using both MARTINI v2.2 and the reparameterized MARTINI model. While the original MARTINI model predicts oligomerization of protein into a single aggregate, the reparameterized MARTINI model maintains a dynamic equilibrium between monomers and dimers as predicted by experimental studies. The proposed reparameterization is an alternative to the standard MARTINI model for use in simulations of realistic models of a biological membrane containing diverse lipids and proteins.

Atomistic mechanism of transmembrane helix association.

Transmembrane helix association is a fundamental step in the folding of helical membrane proteins. The prototypical example of this association is formation of the glycophorin dimer. While its structure and stability have been well-characterized experimentally, the detailed assembly mechanism is harder to obtain. Here, we use all-atom simulations within phospholipid membrane to study glycophorin association. We find that initial association results in the formation of a non-native intermediate, separated by a significant free energy barrier from the dimer with a native binding interface. We have used transition-path sampling to determine the association mechanism. We find that the mechanism of the initial bimolecular association to form the intermediate state can be mediated by many possible contacts, but seems to be particularly favoured by formation of non-native contacts between the C-termini of the two helices. On the other hand, the contacts which are key to determining progression from the intermediate to the native state are those which define the native binding interface, reminiscent of the role played by native contacts in determining folding of globular proteins. As a check on the simulations, we have computed association and dissociation rates from the transition-path sampling. We obtain results in reasonable accord with available experimental data, after correcting for differences in native state stability. Our results yield an atomistic description of the mechanism for a simple prototype of helical membrane protein folding.

Non-synonymous amino acid alterations in PfEBA-175 modulate the merozoite ligand's ability to interact with host's Glycophorin A receptor.

The pathological outcome of malaria due to Plasmodium falciparum infection depends largely on erythrocyte invasion by blood-stage merozoites which employ a cascade of interactions occurring between parasite ligands and RBC receptors. In a previous study exploring the genetic diversity of region-II of PfEBA-175, a ligand that plays a crucial part in parasite's RBC entry through Glycophorin A (GPA) receptor, we demonstrated that F2 domain of region-II underwent positive selection in Indian P. falciparum population through the accumulation of non-synonymous polymorphisms. Here, we examine the functional impact of two highly prevalent non-synonymous alterations in F2, namely Q584E & E592A, using a battery of molecular, biophysical and in-silico techniques. Application of circular dichroism, FTIR, fluorescence spectroscopy reveals that secondary and three-dimensional folding of recombinant-F2 protein carrying 584E and 592A residues (F2-Mut) differs significantly from that carrying 584Q and 592E (F2-3D7). A comparison of spectroscopic and thermodynamic parameters shows that F2-Mut is capable of forming a complex with GPA with higher efficiency compared to F2-3D7. In silico docking predicts both artemisinin and artesunate possess the capacity of slipping into the GPA binding crevices of PfEBA-175 and disrupt PfEBA-GPA association. However, the estimated affinity of artesunate towards PfEBA-175 with 584E and 592A residues is higher than that of artemisinin. Thermodynamic parameters computed using isotherms are concordant with this in-silico prediction. Together, our data suggest that the presence of amino acid alterations in F2 provide structural and functional stability favoring PfEBA-GPA interaction and artesunate can efficiently disrupt the interaction between GPA and PfEBA-175 even carrying altered amino acid residues. The present study alerts the malaria research community by presenting evidence that the parasite is gaining evolutionary fitness by cultivating genetic alterations in many of its proteins.

Non-uniform distribution of myosin-mediated forces governs red blood cell membrane curvature through tension modulation.

The biconcave disk shape of the mammalian red blood cell (RBC) is unique to the RBC and is vital for its circulatory function. Due to the absence of a transcellular cytoskeleton, RBC shape is determined by the membrane skeleton, a network of actin filaments cross-linked by spectrin and attached to membrane proteins. While the physical properties of a uniformly distributed actin network interacting with the lipid bilayer membrane have been assumed to control RBC shape, recent experiments reveal that RBC biconcave shape also depends on the contractile activity of nonmuscle myosin IIA (NMIIA) motor proteins. Here, we use the classical Helfrich-Canham model for the RBC membrane to test the role of heterogeneous force distributions along the membrane and mimic the contractile activity of sparsely distributed NMIIA filaments. By incorporating this additional contribution to the Helfrich-Canham energy, we find that the RBC biconcave shape depends on the ratio of forces per unit volume in the dimple and rim regions of the RBC. Experimental measurements of NMIIA densities at the dimple and rim validate our prediction that (a) membrane forces must be non-uniform along the RBC membrane and (b) the force density must be larger in the dimple than the rim to produce the observed membrane curvatures. Furthermore, we predict that RBC membrane tension and the orientation of the applied forces play important roles in regulating this force-shape landscape. Our findings of heterogeneous force distributions on the plasma membrane for RBC shape maintenance may also have implications for shape maintenance in different cell types.

An unusual variant glycophorin expressing protease-resistant M antigen encoded by the GYPB-E(2-4)-B hybrid gene.

BACKGROUND AND OBJECTIVES: MNS is a highly polymorphic blood group comprising 49 antigens recognized by International Society of Blood Transfusion, some of which may have been generated by genomic recombination among the closely linked genes GYPA, GYPB and GYPE. The GYPE gene has an almost identical sequence to GYPA*01 allele in exon 2 (99% homology), which accounts for M antigen. We investigated an unusual glycophorin molecule with protease-resistant M antigen. METHODS: Blood samples were screened by an automated blood typing system (PK7300) using bromelain-treated red blood cells (RBCs) and murine monoclonal anti-M. The M-positive RBC samples were analysed by immunoblotting using anti-M as the primary antibody. GYPA, GYPB and GYPE genes were analysed by polymerase chain reaction (PCR), cloning and sequencing using reticulocyte mRNA and genomic DNA. RESULTS: Serological tests and immunoblotting revealed that 103 of the 193 009 individuals (0·0534%) expressed protease-resistant M-active glycophorin having a molecular weight of 20 kDa. All the 103 individuals were S+ s- or S- s+. When reticulocyte mRNA from the individuals with M-active glycophorin (20 kDa) was examined by PCR and cloning followed by sequencing, a novel GYPE-B hybrid transcript was identified. Long-range PCR and sequencing using genomic DNA revealed that the individuals had a GYPB-E(2-4)-B hybrid gene. This hybrid gene was predicted to encode a 59-amino-acid mature glycoprotein that expresses no S or s antigens CONCLUSIONS: The prevalence of the M-active glycophorin (20 kDa) in the Japanese population is 0·0534%. This glycophorin is predicted to be a 59 amino acids polypeptide encoded by the novel GYPB-E(2-4)-B hybrid gene.

Detection and phenotyping of extracellular vesicles by size exclusion chromatography coupled with on-line fluorescence detection.

New methods for quantifying extracellular vesicles (EVs) in complex biofluids are critically needed. We report the development of a new technology combining size exclusion chromatography (SEC), a commonly used EV purification technique, with fluorescence detection of specifically labelled EVs. The resulting platform, Flu-SEC, demonstrates a linear response to concentration of specific EVs and could form the basis of a system with phenotyping capability. Flu-SEC was validated using red blood cell derived EVs (REVs), which provide an ideal EV model with monodisperse size distribution and high EV concentration. Microfluidic Resistive Pulse Sensing (MRPS) was used to accurately determine the size distribution and concentration of REVs. Anti-CD235a antibody, specific to glycophorin A, and the more general wheat germ agglutinin (WGA), were selected to label REVs. The results show the quantitative power of Flu-SEC: a highly linear fluorescence response over a wide range of concentrations. Moreover, the Flu-SEC technique reports the ratio of EV-bound and free-antibody molecules, an important metric for determining optimal labelling conditions for other applications. Flu-SEC represents an orthogonal tool to single-particle fluorescent methods such as flow cytometry and fluorescent NTA, for the quantification and phenotyping of EVs.

Accelerating Membrane Simulations with Hydrogen Mass Repartitioning.

The time step of atomistic molecular dynamics (MD) simulations is determined by the fastest motions in the system and is typically limited to 2 fs. An increasingly popular approach is to increase the mass of the hydrogen atoms to ∼3 amu and decrease the mass of the parent atom by an equivalent amount. This approach, known as hydrogen-mass repartitioning (HMR), permits time steps up to 4 fs with reasonable simulation stability. While HMR has been applied in many published studies to date, it has not been extensively tested for membrane-containing systems. Here, we compare the results of simulations of a variety of membranes and membrane-protein systems run using a 2 fs time step and a 4 fs time step with HMR. For pure membrane systems, we find almost no difference in structural properties, such as area-per-lipid, electron density profiles, and order parameters, although there are differences in kinetic properties such as the diffusion constant. Conductance through a porin in an applied field, partitioning of a small peptide, hydrogen-bond dynamics, and membrane mixing show very little dependence on HMR and the time step. We also tested a 9 Å cutoff as compared to the standard CHARMM cutoff of 12 Å, finding significant deviations in many properties tested. We conclude that HMR is a valid approach for membrane systems, but a 9 Å cutoff is not.

Interaction of the human erythrocyte Band 3 anion exchanger 1 (AE1, SLC4A1) with lipids and glycophorin A: Molecular organization of the Wright (Wr) blood group antigen.

The Band 3 (AE1, SLC4A1) membrane protein is found in red blood cells and in kidney where it functions as an electro-neutral chloride/bicarbonate exchanger. In this study, we have used molecular dynamics simulations to provide the first realistic model of the dimeric membrane domain of human Band 3 in an asymmetric lipid bilayer containing a full complement of phospholipids, including phosphatidylinositol 4,5-bisphosphate (PIP2) and cholesterol, and its partner membrane protein Glycophorin A (GPA). The simulations show that the annular layer in the inner leaflet surrounding Band 3 was enriched in phosphatidylserine and PIP2 molecules. Cholesterol was also enriched around Band 3 but also at the dimer interface. The interaction of these lipids with specific sites on Band 3 may play a role in the folding and function of this anion transport membrane protein. GPA associates with Band 3 to form the Wright (Wr) blood group antigen, an interaction that involves an ionic bond between Glu658 in Band 3 and Arg61 in GPA. We were able to recreate this complex by performing simulations to allow the dimeric transmembrane portion of GPA to interact with Band 3 in a model membrane. Large-scale simulations showed that the GPA dimer can bridge Band 3 dimers resulting in the dynamic formation of long strands of alternating Band 3 and GPA dimers.

Balancing Force Field Protein-Lipid Interactions To Capture Transmembrane Helix-Helix Association.

Atomistic simulations have recently been shown to be sufficiently accurate to reversibly fold globular proteins and have provided insights into folding mechanisms. Gaining similar understanding from simulations of membrane protein folding and association would be of great medical interest. All-atom simulations of the folding and assembly of transmembrane protein domains are much more challenging, not least due to very slow diffusion within the lipid bilayer membrane. Here, we focus on a simple and well-characterized prototype of membrane protein folding and assembly, namely the dimerization of glycophorin A, a homodimer of single transmembrane helices. We have determined the free energy landscape for association of the dimer using the CHARMM36 force field. We find that the native structure is a metastable state, but not stable as expected from experimental estimates of the dissociation constant and numerous experimental structures obtained under a variety of conditions. We explore two straightforward approaches to address this problem and demonstrate that they result in stable dimers with dissociation constants consistent with experimental data.

Genetic structure of two erythrocyte binding antigens of Plasmodium falciparum reveals a contrasting pattern of selection.

Erythrocyte binding antigens 175 (EBA-175) and 140 (EBA-140) play key roles in erythrocyte invasion by binding to glycophorin A (GPA) and C (GPC) respectively in human malaria. Since antigenic variation in malaria endemic region is a major barrier to development of effective vaccine, we explore the nature and pattern of sequence diversity of these two vaccine candidates in Kolkata, India. Population genetic parameters based on parasite sequences representing region II of Pfeba-175 and Pfeba-140 genes were estimated using DnaSP V.5.10 and MEGA version 6.0. A novel molecular docking approach was implemented to assess the binding affinities of Kolkata Pfeba-175 variants with GPA. P. falciparum Kolkata isolates experienced a recent population expansion as documented by negative Tajima's D, Fu & Li's statistics, unimodal mismatch distribution and star-like median-joining network for both loci. Positive selection seemed to play a major role in shaping the diversity of Pfeba-175 (dN/dS=2.45, and McDonald-Kreitman P-value=0.04) with successive accumulation of Q584K/E, E592A and R664S deriving high frequency haplotypes designated here as F2KH3 and F2KH1. In silico molecular docking demonstrated that polypeptides encoded by F2KH1 and F2KH3 were capable of engaging the parasite ligand into energetically favorable interaction with GPA. Our data demonstrated emergence of Pfeba-175 sequences harboring selectively advantageous nonsynonymous substitutions on Pf3D7 sequence background in the Kolkata parasite isolates. A contrasting pattern of Pf3D7-centric expansion of parasite sequences was noted for Pfeba-140. Together, this study provides a firm genetic and biological support favoring a dominant role of EBA-175 in erythrocyte invasion.

Combination of Cα-H Hydrogen Bonds and van der Waals Packing Modulates the Stability of GxxxG-Mediated Dimers in Membranes.

The GxxxG motif is frequently found at the dimerization interface of a transmembrane structural motif called GASright, which is characterized by a short interhelical distance and a right-handed crossing angle between the helices. In GASright dimers, such as glycophorin A (GpA), BNIP3, and members of the ErbB family, the backbones of the helices are in contact, and they invariably display networks of 4 to 8 weak hydrogen bonds between Cα-H carbon donors and carbonyl acceptors on opposing helices (Cα-H···O═C hydrogen bonds). These networks of weak hydrogen bonds at the helix-helix interface are presumably stabilizing, but their energetic contribution to dimerization has yet to be determined experimentally. Here, we present a computational and experimental structure-based analysis of GASright dimers of different predicted stabilities, which show that a combination of van der Waals packing and Cα-H hydrogen bonding predicts the experimental trend of dimerization propensities. This finding provides experimental support for the hypothesis that the networks of Cα-H hydrogen bonds are major contributors to the free energy of association of GxxxG-mediated dimers. The structural comparison between groups of GASright dimers of different stabilities reveals distinct sequence as well as conformational preferences. Stability correlates with shorter interhelical distances, narrower crossing angles, better packing, and the formation of larger networks of Cα-H hydrogen bonds. The identification of these structural rules provides insight on how nature could modulate stability in GASright and finely tune dimerization to support biological function.

Identification of lectin counter-receptors on cell membranes by proximity labeling.

Lectin-glycan interactions play important roles in many biological systems, but the nature of glycoprotein counter-receptors expressed on cell membranes is often poorly understood. To help overcome this problem, we developed a method based on proximity labeling technology. Using a peroxidase-coupled lectin, addition of H2O2 and tyramide-biotin substrates leads to generation of short-range biotin radicals that biotinylate proteins in the immediate vicinity of the bound lectin, which can subsequently be identified. As a proof-of-principle, sialoadhesin-horseradish peroxidase-human IgG1 Fc recombinant protein constructs were precomplexed with anti-Fc antibodies, bound to human erythrocytes and reacted with H2O2 and tyramide-SS-biotin. The erythrocyte membrane protein with strongest biotinylation was identified as glycophorin A, in agreement with early studies using lectin overlay and reglycosylation approaches. As a further test of the method, the plant lectin MAL II was conjugated with horseradish peroxidase and used in proximity labeling of human erythrocytes. Glycophorin A was again selectively labeled, which is consistent with previous reports that MAL II has high affinity for glycophorin. This method could be applied to other lectins to identify their membrane counter-receptors.

Resistance to malaria through structural variation of red blood cell invasion receptors.

The malaria parasite Plasmodium falciparum invades human red blood cells by a series of interactions between host and parasite surface proteins. By analyzing genome sequence data from human populations, including 1269 individuals from sub-Saharan Africa, we identify a diverse array of large copy-number variants affecting the host invasion receptor genes GYPA and GYPB We find that a nearby association with severe malaria is explained by a complex structural rearrangement involving the loss of GYPB and gain of two GYPB-A hybrid genes, which encode a serologically distinct blood group antigen known as Dantu. This variant reduces the risk of severe malaria by 40% and has recently increased in frequency in parts of Kenya, yet it appears to be absent from west Africa. These findings link structural variation of red blood cell invasion receptors with natural resistance to severe malaria.

Identifying ionic interactions within a membrane using BLaTM, a genetic tool to measure homo- and heterotypic transmembrane helix-helix interactions.

The assembly of integral membrane protein complexes is frequently supported by transmembrane domain (TMD) interactions. Here, we present the BLaTM assay that measures homotypic as well as heterotypic TMD-TMD interactions in a bacterial membrane. The system is based on complementation of ß-lactamase fragments genetically fused to interacting TMDs, which confers ampicillin resistance to expressing cells. We validated BLaTM by showing that the assay faithfully reports known sequence-specific interactions of both types. In a practical application, we used BLaTM to screen a focussed combinatorial library for heterotypic interactions driven by electrostatic forces. The results reveal novel patterns of ionizable amino acids within the isolated TMD pairs. Those patterns indicate that formation of heterotypic TMD pairs is most efficiently supported by closely spaced ionizable residues of opposite charge. In addition, TMD heteromerization can apparently be driven by hydrogen bonding between basic or between acidic residues.

Quantification of Membrane Protein Self-Association with a High-Throughput Compatible Fluorescence Assay.

For many membrane proteins, self-association serves both structural and functional roles. Studies of such association can be simplified by switching to micelles as the membrane-mimicking environment, but native interaction is not preserved in all detergents. The selection of suitable conditions for biochemical experiments would be greatly facilitated by a quantitative high-throughput assay. Here we showed that the fluorescence polarization reduction, which resulted from homo-Förster resonance energy transfer and was measured in a high-throughput compatible format, can be used to determine both association states and constants for membrane proteins in micelles.

Classification of Tandem Mass Spectra for Identification of N- and O-linked Glycopeptides.

Analysis of intact glycopeptides by mass spectrometry is essential to determining the microheterogeneity of protein glycosylation. Higher-energy collisional dissociation (HCD) fragmentation of glycopeptides generates mono- or disaccharide ions called oxonium ions that carry information about the structure of the fragmented glycans. Here, we investigated the link between glycan structures and the intensity of oxonium ions in the spectra of glycopeptides and utilized this information to improve the identification of glycopeptides in biological samples. Tandem spectra of glycopeptides from fetuin, glycophorin A, ovalbumin and gp120 tryptic digests were used to build a spectral database of N- and O-linked glycopeptides. Logistic regression was applied to this database to develop model to distinguish between the spectra of N- and O-linked glycopeptides. Remarkably, the developed model was found to reliably distinguish between the N- and O-linked glycopeptides using the spectral features of the oxonium ions using verification spectral set. Finally, the performance of the developed predictive model was evaluated in HILIC enriched glycopeptides extracted from human serum. The results showed that pre-classification of tandem spectra based on their glycosylation type improved the identification of N-linked glycopeptides. The developed model facilitates interpretation of tandem mass spectrometry data for assignment of glycopeptides.