Juan
Wei
*ab,
Dimitris
Papanastasiou
c,
Mariangela
Kosmopoulou
c,
Athanasios
Smyrnakis
c,
Pengyu
Hong
d,
Nafisa
Tursumamat
a,
Joshua A.
Klein
b,
Chaoshuang
Xia
b,
Yang
Tang
be,
Joseph
Zaia
b,
Catherine E.
Costello
be and
Cheng
Lin
*b
aShanghai Jiao Tong University, 800 Dongchuan Road, Shanghai, 200240, China
bCenter for Biomedical Mass Spectrometry, Boston University Chobanian & Avedisian School of Medicine, Boston, MA 02118, USA. E-mail: chenglin@bu.edu
cFasmatech Science and Technology, 15310 Athens, Greece
dDepartment of Computer Science, Brandeis University, Waltham, MA 02454, USA
eDepartment of Chemistry, Boston University, Boston, MA 02215, USA
First published on 2nd June 2023
Comprehensive de novo glycan sequencing remains an elusive goal due to the structural diversity and complexity of glycans. Present strategies employing collision-induced dissociation (CID) and higher energy collisional dissociation (HCD)-based multi-stage tandem mass spectrometry (MSn) or MS/MS combined with sequential exoglycosidase digestions are inherently low-throughput and difficult to automate. Compared to CID and HCD, electron transfer dissociation (ETD) and electron capture dissociation (ECD) each generate more cross-ring cleavages informative about linkage positions, but electronic excitation dissociation (EED) exceeds the information content of all other methods and is also applicable to analysis of singly charged precursors. Although EED can provide extensive glycan structural information in a single stage of MS/MS, its performance has largely been limited to FTICR MS, and thus it has not been widely adopted by the glycoscience research community. Here, the effective performance of EED MS/MS was demonstrated on a hybrid Orbitrap-Omnitrap QE-HF instrument, with high sensitivity, fragmentation efficiency, and analysis speed. In addition, a novel EED MS2-guided MS3 approach was developed for detailed glycan structural analysis. Automated topology reconstruction from MS2 and MS3 spectra could be achieved with a modified GlycoDeNovo software. We showed that the topology and linkage configurations of the Man9GlcNAc2 glycan can be accurately determined from first principles based on one EED MS2 and two CID-EED MS3 analyses, without reliance on biological knowledge, a structure database or a spectral library. The presented approach holds great promise for autonomous, comprehensive and de novo glycan sequencing.
Major limitations of collision-based MS/MS arise from the slow-heating nature of collisional activation, and may be overcome by employing alternative, radical-driven ion activation methods, such as free radical activated glycan sequencing,21,22 ultraviolet photodissociation,23–25 charge transfer dissociation,26 and various electron-activated dissociation (ExD) methods.15,27–34 Among them, electronic excitation dissociation (EED) MS/MS has recently emerged as a powerful tool for structural glycomics. EED can produce detailed structural information in a single stage of MS/MS analysis, and has been successfully implemented with on-line LC and ion mobility separation for effective characterization of glycan mixtures.34–40 The power of an integrated approach that combines EED MS/MS with on-line chromatographic separation and software-assisted spectral interpretation was demonstrated in a recent study that employed porous graphitic carbon (PGC)-LC-EED-MS/MS.40 Candidate topologies of 18 oligomannose glycans released from RNase B were reconstructed from their EED MS2 spectra by GlycoDeNovo,41 a de novo glycan sequencing software, without reliance on any database. Putative topologies were consistently ranked as the top structures by the software.40 Many linkages could also be determined de novo, although in some cases, biological insights, such as the glycan biosynthetic rules, were needed to fully define the structure.
To date, glycan structural elucidation by EED MS/MS has only been demonstrated on Fourier transform-ion cyclotron resonance (FTICR) MS instruments. Broad adoption of EED MS/MS by the glycoscience research community has been hindered by the limited accessibility and high operating cost of FTICR MS. Recently, electromagnetostatic cells that may be installed in non-ICR instruments have become commercially available, allowing protein sequencing by electron-based dissociation on Orbitrap, time-of-flight (TOF), and ion trap mass spectrometers.42–48 However, effective ExD characterization of glycans on these instruments has yet to be fully realized, in part because ion-electron interaction occurs during analyte transfer through the ExD cell without trapping, and this limits the ExD efficiency, particularly for higher-energy ExD processes. Improved electron activated dissociation (EAD) efficiency may be achieved by trapping the precursor ions for more extensive interactions with electrons.49
In this study, we took advantage of a new ExD cell design on a prototype Omnitrap instrument.50Fig. 1a shows the schematic of an Omnitrap platform, consisting of a system of linear ion trap segments, connected to the rear of the HCD cell of a QE-HF Orbitrap instrument via a transfer hexapole. The segmented ion trap offers flexibilities for ion storage and manipulation, enabling high performance MSn analysis by utilizing a wide range of fragmentation methods, including CID and ExD. Interfacing the Omnitrap with an Orbitrap allows subsequent detection of fragment ions with high mass accuracy and resolving power, and this is beneficial for MS/MS-based glycan structural analysis as isobaric fragments are frequently produced by glycans. The principle of operation for ExD experiments on a QE-Omnitrap instrument is illustrated in Fig. 1b. ExD analysis is carried out in Q5 with an orthogonally mounted electron source consisting of a tantalum disc cathode and electrostatic focusing lenses. Electrons are injected into Q5 when a positive potential is applied to the X-electrodes and blocked when a negative potential is applied. The E2 split electrode provides a deflection pulse to ensure electron injection only during the correct RF phase. Unlike conventional ion traps that employ sinusoidal RF waveforms, rectangular waveforms are applied to the Omnitrap electrodes. The rectangular waveform produces a “static’ trapping field, allowing precise control of the electron energy. Time-controlled ion-electron interaction is possible due to analyte trapping in Q5 during the entire ExD event, and this should lead to improved ExD efficiencies.
Here, we report, for the first time, efficient detailed glycan structural characterization by EED MS/MS on a hybrid QE-Omnitrap instrument. To address the challenges in de novo glycan sequencing by MS2, we developed a novel EED MS2-guided MS3 strategy to differentiate iso-topologies and to eliminate ambiguities in linkage assignment among different antennae. Additionally, the GlycoDeNovo software was modified to facilitate automated interpretation of both MS2 and MS3 spectra. The potential of this approach for de novo glycan sequencing is demonstrated for characterization of the canonical Man9 chitobiose structure (Man9GlcNAc2) and the biantennary G2 glycan.
A comparative study was performed on a 12T solariX FTICR MS equipped with a hollow cathode dispenser (Bruker Daltonics, Bremen, Germany). Samples were loaded into a pulled fused silica capillary and directly infused into the mass spectrometer by nano-electrospray. Targeted ions were isolated by the front-end quadrupole and accumulated in an external hexapole collision cell for up to 400 ms. For EED MS/MS analysis, precursor ions were irradiated by 18–20 eV electrons for 100–300 ms. For MS3 analysis, quadrupole-selected precursors were accumulated and fragmented by CID in the collision cell, and fragments were sent to the ICR cell and selected by an in-cell sweep isolation for EED MS3 analysis. Mass spectra were acquired with a low-mass cutoff of m/z 200, and a 512k or 1 M data size, resulting in a transient length of 288 ms or 577 ms, and an estimated resolving power of 66 K or 130 K at m/z 400, respectively. Five and forty spectra were averaged for EED MS/MS and CID-EED MS3 analyses, respectively.
Fig. 2 shows the EED MS2 spectra of the doubly-sodiated Man9GlcNAc2 precursor (m/z 1218.1042) acquired on the QE-Omnitrap and on the FTICR MS, with all assigned fragments listed in ESI Tables S1 and S2,† respectively. EED produced very similar fragmentation patterns on these two instruments, while the EED efficiency achieved on QE-Omnitrap was significantly higher than on FTICR MS, producing on average 5–10 times higher relative fragment abundance and signal-to-noise ratio. With a higher S/N ratio, 31 more peaks were assigned in the QE-Omnitrap EED MS2 spectrum than in the corresponding FTICR spectrum, representing a 25% increase in the number of assigned peaks. Notably, the Omnitrap spectrum of the five-fold diluted sample solution was acquired with a shorter ion injection time (5 ms vs. 400 ms) and electron irradiation time (50 ms vs. 200 ms) than the ICR spectrum. The shorter spectral acquisition time and higher EED efficiency achieved on QE-Omnitrap make it far more suitable for high-throughput LC-MS/MS analysis, allowing characterization of a higher number of glycoforms, including those present in lower abundance. The higher EED efficiency achieved on the QE-Omnitrap instrument may have benefited from collisional cooling and focusing in the higher-pressure Q5 that improves the electron–ion interaction, whereas in the ultra-high vacuum of the ICR cell, axial excitation and magnetron expansion can lead to poorer ion-electron overlap.54
Fig. 2 EED MS2 spectra and cleavage maps of reduced- and permethylated Man9GlcNAc2 ([M + 2Na]2+ at m/z 1218.1042) acquired on (a) QE-Omnitrap MS and (b) FTICR MS, respectively. Each EED MS2 spectrum was the summed result of 5 scans. A complete series of C-ions and Y/Z/1,5X triplets are labeled in red and blue, respectively. Linkage-diagnostic cross-ring and internal fragments are labeled in green. Asterisk indicates singly charged fragments with two sodium atoms. A complete list of assigned fragments can be found in ESI Tables S1 and S2.† |
On either instrument platform, EED of the Man9GlcNAc2 glycan produced a complete series of C-type ions (C1/C1‡, C2/C2‡, C3β/C3β‡, C3α/C3α‡, C4‡, C5‡, where the symbol ‡ indicates loss of two hydrogens), and Y/Z/1,5X triplets. The presence of these sequence ions allows de novo reconstruction of the Man9GlcNAc2 candidate topologies by GlycoDeNovo. The canonical structure, listed in bold in Fig. 3a with its topology encircled in Fig. 3b, is one of the five candidates with the highest IonClassifier score. The IonClassifier was established by a machine learning approach that identifies common spectral context of a certain type of fragment in a training data set, defined as a collection of its related neighboring peaks, represented by their respective mass shifts and relative abundances. The IonClassifier can be used to evaluate the accuracy of an assignment by examining the spectral context of the fragment.41 The IonClassifier score of a candidate topology is the sum of the IonClassifier scores of all its supporting peaks, namely, glycosidic fragments consistent with the topology. Here, the five top-ranked topologies are indistinguishable at the MS2 level, as they produce the same set of supporting peaks, including non-reducing-end glycosidic fragments with compositions of Hex, Hex2, Hex3, Hex5, Hex9, and Hex9HexNAc. Previously, the correct Man9GlcNAc2 topology was assigned based on prior glycan biosynthetic knowledge.40 Among the five top-ranked candidates, only topology 2, with a linear tri-hexose branch and a branched penta-hexose branch, can be derived from the structure of the tetradecasaccharide N-linked glycan precursor (Fig. 3c).
Fig. 3 (a) Top-ranked candidate topologies reconstructed by GlycoDeNovo from the EED MS2 spectrum of deutero-reduced and permethylated Man9GlcNAc2 in Fig. 2; (b) graphic representations of the five top-ranked isotopolomers of Man9GlcNAc2; (c) SNFG representation of the N-glycan precursor. |
Once a topology is defined, most linkages can be independently determined based on linkage-diagnostic fragments without further reference to the precursor structure. For example, the presence of 3,5A4 (m/z 1145.5566), 0,4A4 (m/z 1117.5256), and C4/Z3β (m/z 1247.5889) ions not only places the pentasaccharide branch to the 6-antenna, but also localizes the tri-hexose branch to the 3-antenna, since a high-abundance C/Z ion is associated with the loss of the C3-substituent from a Z ion. The observation of 1,3A4 (m/z 723.3408) and 0,2X2 (m/z 618.3319) ions also corroborates with the assignment of the tri-hexose to the 3-antenna. Within the 6-antenna, the presence of 0,4A3α (m/z 505.2253), 3,5A3α (m/z 533.2568), C3α/Z4α′′ (m/z 635.2886), and 0,2X3α (m/z 728.8606) ions is consistent with a penta-hexose structure consisting of two di-hexose branches, at the 6- and 3-positions of the branching mannose, respectively, although the 1→6 linkage assignment is not conclusive, as the observed 0,4A3α and 3,5A3α ions can theoretically be generated from the linear tri-hexose 3-antenna as well. Similarly, the 1→2 linkage between the second and third mannose residues on the 3-antenna cannot be assigned definitively, since its diagnostic 1,3A3β (m/z 519.2411) ion is isomeric to the 1,3A3α′′ at the 6-antenna. Finally, while the presence of a 1,3A2 ion (m/z 315.1410) and the absence of 0,2X4, 3,5A2, and 0,4A2 ions strongly suggest the existence of at least one 1→2-linked non-reducing-end mannose residue, its location (or their locations) cannot be unambiguously determined. ESI Scheme S1b† shows a hypothetical structure that can produce the same group of cross-ring and internal fragments as the true structure (ESI Scheme S1a†). Detailed discussions on the general EED mechanism and formation of specific diagnostic fragments have been presented elsewhere.36,37,40,55
The analysis above clearly illustrates the challenges of MS2-based de novo glycan sequencing. Even with complete series of glycosidic fragments and cross-ring fragments, it is sometimes not possible to differentiate iso-topologies or accurately determine linkages among different branches. Additional structural details may be revealed by performing sequential tandem MS analysis on MS2 fragment(s) of interest. Presently, there are several ways to choose MS2 fragments as precursors for later stages of MSn analysis, each with its merits and limitations. The most straightforward method is to select product ion(s) of the highest abundance. Although such an approach allows autonomous precursor selection, the highest-abundance fragments do not always produce the most structurally informative MS3 spectra. Alternatively, precursors may be selected by expert users, aided by existing knowledge of the glycan fragmentation behavior and/or biological insights.56–58 This is typically performed during direct infusion analyses on ion trap instruments, whose low mass resolving power can negatively impact the accuracy of the fragment assignment. Further, without chromatographic separation, each MS2 precursor may consist of multiple isomeric structures, and each MS2 fragment chosen for further MSn analysis can be a mixture of isobaric structures. Consequently, confident structural assignment requires investigation of many gas-phase disassembly pathways, often involving deeper level of MSn. This process is difficult to automate, and the need to perform MSn>3 limits both the sensitivity and throughput. Recently, Sun and coworkers developed a glycan intelligent precursor selection (GIPS) strategy to guide MSn experiments.59,60 GIPS utilizes a statistical model to calculate the distinguishing power of fragments, which dictates the selection of precursor(s) for next stage of tandem MS analyses. The GIPS approach has only been demonstrated for the identification of the glycan branching patterns, but not for linkage determination. Moreover, when calculating the distinguishing power, it only considers structures present in the existing glycan database, and thus it cannot identify new structures.
The MSn method may also be combined with spectral library search for glycan structural elucidation, as utilized in logically derived sequence (LODES)/MSn.20,61,62 The LODES/MSn approach employs a built-in logical procedure to successively break down an oligosaccharide sequence into its monosaccharide and disaccharide components. The CID spectra of these smaller fragments can then be searched against a spectral database for structural assignment. The LODES/MSn approach can provide detailed structural information, provided that the corresponding spectral library of disaccharide standards is available. However, it still requires investigation of many MSn pathways, and is not readily compatible with on-line LC-MS analysis. Multiple LC injections were needed to fully define the structure of glycans as small as trisaccharides.62 For larger oligosaccharides, such as released N-linked glycans, knowledge on their biosynthesis was utilized to reduce the number of MSn pathways that need to be examined. For complex mixtures, it was necessary to perform multi-dimensional LC fractionation off-line before LODES/MSn analysis.20
Here, precursor(s) for MS3 analyses were selected based on their projected differentiating power on top-ranked topology candidates identified by GlycoDeNovo from the EED MS2 spectrum. For Man9GlcNAc2, the top five topologies (Fig. 3b) can be uniquely differentiated by performing MS3 analysis on the CID-generated Hex3 and Hex5 fragments (Fig. 4a). EED of a Hex5 fragment, 0,4A4 at m/z 1117.5259, produced glycosidic fragments with either one or two Hex units (C1 ion at m/z 259.1162, B2 at m/z 445.2045, and C2 at m/z 463.2153), but not with three or four Hex units (Fig. 3b), hence eliminating topologies 1, 3, and 5 from consideration. The remaining two topologies, 2 and 4, can be distinguished by EED analysis of a Hex3 fragment, B3β at m/z 649.3046. The presence of two Hex2 fragments (B2 at m/z 445.2048 and C2 at m/z 463.2162) in the EED spectrum of B3β (Fig. 3c) clearly established topology 2 as the only structure consistent with the MS3 results. The same conclusion could be reached by applying a modified GlycoDeNovo algorithm for analysis of the MS3 spectra. In this case, the software treated the mass difference between the intact glycan and the MS2 fragment as a reducing-end modification. GlycoDeNovo successfully identified a penta-saccharide with two di-hexose branches and a linear tri-hexose structure as the top-ranked topologies for the Hex5 and Hex3 substructures, respectively. Thus, we showed that the correct topology of the Man9GlcNAc2 glycan can be defined by just two MS3 analyses, autonomously from the first principles, with no need for reference spectral libraries, or biological knowledge. Here, the applicability of EED MS/MS to analysis of singly charged precursor is essential for the CID-EED MS3 workflow, as both the 0,4A4 and B3β fragments selected for MS3 analysis were singly charged, and therefore could not be characterized by charge-reducing MS/MS methods, such as ETD and ECD.
Fig. 4 (a) CID MS2 spectra of deutero-reduced- and permethylated Man9GlcNAc2 ([M + 2Na]2+ at m/z 1218.1042) acquired on the QE-Omnitrap system. Fragments selected for EED MS3 are labeled. (b) and (c) CID-EED MS3 spectra and cleavage maps of the 0,4A4 (m/z 1117.5259) and B3β (m/z 649.3046) ions, respectively, acquired on QE-Omnitrap. Linkage-diagnostic cross-ring and internal fragments are labeled in green. Peaks marked by asterisks are background peaks with mass defects that are not expected from glycan fragments. Complete lists of assigned fragments can be found in ESI Tables S3 and S4.† |
It is worth noting that a given substructure may be present in several fragments. For MS2 fragments with the same differentiating power, the choice is generally made based on the fragment ion abundance and the ability to achieve a clean isolation. Sometimes, more than one MS2 fragment can generate similarly informative MS3 spectra. For example, the Y3β2+ ion at m/z 904.9469 also retained the Man5 substructure without interference from the Man3 branch. EED of Y3β2+ (ESI Fig. S1 and Table S5†) similarly produced non-reducing-end glycosidic fragments with compositions of Hex1 or Hex2, but not Hex3 or Hex4, consistent with the branched Hex5 substructure assignment based on the EED spectrum of the 0,4A4 ion.
As discussed earlier, once topology 2 is established as the correct structure, the Hex5 and Hex3 branches can be accurately located to the 6- and 3-antennae, respectively, on the basis of the characteristic cross-ring and abundant C/Z internal fragments observed in the EED MS2 spectrum. Within the 6-antenna, one of the di-hexose branches can be assigned to the 3-position based on the presence of a high-abundance C3/Z4α′′ ion at m/z 635.2886, but the location of the other di-hexose branch cannot be confidently assigned due to potential interference by fragments from the linear tri-hexose 3-antenna. With MS3, the Hex5 substructure can be isolated for accurate linkage analysis. EED MS3 of the 0,4A4 ion generated both the 0,4A3 ion at m/z 505.2262 and the 3,5A3 ion at m/z 533.2575, thus unambiguously assigning the second di-hexose branch to the 6-position. The presence of the 1,3A2 ion at m/z 315.1416 and a C‡/Y (Hex) ion at m/z 243.0844, together with the lack of a 0,2X4β/0,4A4 ion at m/z 751.3359, placed the terminal mannose residues to the 2-position. Similarly, the linear tri-hexose branch at the 3-antenna was isolated in the B3β fragment, whose EED MS3 analysis generated 1,3A2 at m/z 315.1415, 1,3A3 at m/z 519.2409, C‡/Y (Hex) ion at m/z 243.0842, but 0,2X4β at m/z 283.1149, supporting the 1→2 linkages between the mannose residues. Collectively, the linkages within both antennae could be determined de novo by the CID-EED MS3 analysis.
With the Omnitrap-Orbitrap system, the MS3 experiment can be performed in several different configurations. The CID-EED sequence was chosen here as CID typically generates a smaller number of fragments in higher abundance than EED. Thus, the CID spectrum is usually not as congested as the EED spectrum, making it easier to achieve a clean isolation of the high-abundance MS2 fragment of interest. For Man9GlcNAc2, isolation of the 0,4A4 fragment at m/z 1117.5256 from EED MS2 can be complicated by co-isolation of the doubly charged 3,5A6 ion at m/z 1115.0348 (Fig. 2a and Table S1†), whereas the interfering 3,5A6 ion was absent in the CID spectrum of Man9GlcNAc2 (Fig. 4a). On the other hand, EED is the preferred final-stage dissociation method, producing MS3 spectra with much higher structural information content than CID. The EED and CID MS3 cleavage maps of the 0,4A4 ion are shown in Fig. 4b and ESI Fig. S2,† respectively. It is evident that many linkage-diagnostic fragments, such as 1,3A and C‡/Y (Hex) ions, were absent in the CID spectrum (Fig. S2†). Finally, although it is possible to perform CID-EED MS3 analysis on a hybrid Qh-FTICR MS instrument, selection of the MS3 precursor can only be achieved with in-cell isolation, resulting in a much lower ion count for the precursor of interest. Consequently, even with more than 10 times the ion injection time, the CID-EED MS3 spectrum of the 0,4A4 ion on FTICR MS (ESI Fig. S3†) was of much lower quality than that obtained on the Omnitrap-Orbitrap system (Fig. 4b).
While the discussion thus far has been focused on the Man9GlcNAc2 structure, the workflow presented here can be effectively applied to characterize other glycan structures, including complex-type glycans. Note that the Man9GlcNAc2 structure was chosen as the primary example not because of its apparent structural simplicity, but rather, because it represents one of the hardest cases for isotopolomer differentiation. With only one type of monosaccharide residue (hexose) present in all branches, many structures can produce the same set of glycosidic fragments and are therefore impossible to differentiate at the MS2 level, even with the IonClassifier. Two additional MS3 analyses were required to unambiguously define the canonical structure. On the other hand, complex-type glycans consist of a larger variety of monosaccharide building blocks that have different masses and are less likely to produce isotopologies. True topologies of many complex-type glycans can be identified at the MS2 level, though there are cases where MS3 is necessary. Fig. 5 shows one such case for characterization of the biantennary G2 glycan. At the MS2 level, topologies 1 and 2 (T1 and T2) would produce the same set of glycosidic fragments (Hex, HexHexNAc, Hex2HexNAc, etc.), and were co-ranked as the top candidate by the IonClassifier based on the EED MS/MS spectrum of the deutero-reduced and permethylated G2 glycan (ESI Fig. S4†). These two topologies could be differentiated based on the HCD-EED MS3 spectrum of the Y4 ion (ESI Fig. S5†), where the presence of HexHexNAc fragments and their complementary ions supported T1 as the correct topology. Lists of assigned peaks for the EED spectra of the G2 glycan and its Y4 fragment can be found in ESI Tables S6 and S7.†
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d3sc00870c |
This journal is © The Royal Society of Chemistry 2023 |