Multiplexed mass spectrometry of individual ions improves measurement of proteoforms and their complexes.

An Orbitrap-based ion analysis procedure determines the direct charge for numerous individual protein ions to generate true mass spectra. This individual ion mass spectrometry (I2MS) method for charge detection enables the characterization of highly complicated mixtures of proteoforms and their complexes in both denatured and native modes of operation, revealing information not obtainable by typical measurements of ensembles of ions.

For decades, mass spectrometry has used ions to measure the mass-to-charge ratio of molecules once lifted into the gas phase[1]. Denatured and native electrospray ionization of intact proteins and their complexes pose many complications due to sample heterogeneity and large charge state envelopes in the m/z domain. To simplify analysis, Charge Detection Mass Spectrometry (CDMS)[2-5] has enabled the generation of true mass spectra with the direct readout of an ion’s integer charge value (e.g. 35+, 34+, etc.). Here, we bring the commercially available Orbitrap mass analyzer into the CDMS universe to multiplex and regularize this approach for creation of mass spectra on extremely complex samples. The approach is shown to be general for measuring complex proteoform mixtures and their complexes, covering a range of masses from 8 kDa to 3.2 MDa, with Orbitrap-based Individual Ion Mass Spectrometry (I2MS) capable of precise mass determination where ensemble measurements fail or require prior separation.

Previously, single ion sensitivity has been demonstrated using an Orbitrap mass analyzer[6, 7], analogous to linear geometry trapping instruments more frequently used for CDMS. In 2018, we used an Orbitrap analyzer to acquire one-at-a-time measurements on single ions, and showed that >20-fold increases of resolving power over ensemble measurements were possible by centroiding m/z values of individual ions[8]. Here, we extend this methodology substantially with a complete workflow for CDMS that we call I2MS (Fig. 1) to directly measure the charge of individual ions. Further, we show this method to be fully established and functional for the simultaneous measurement of 200 individual ions.

Figure 1.

The five step workflow for the I2MS process.

Includes data acquisition and post-acquisition processing steps that culminate in the plotting of an I2MS spectrum (Step 5). A detailed explanation of each of these five steps is found in the accompanying Online Methods section.

In Step 1 (Fig. 1) of the I2MS approach, ~120 ions are observed per acquisition in a random-style trapping event. Parallelizing ion observation to >100 ions per acquisition decreases the data acquisition burden by ~100-fold over current CDMS techniques[9]. Due to the low ion count needed for I2MS, protein solutions can be diluted down by ~three orders of magnitude (into the high pM to low nM range). In Step 2 (Fig. 1), the frequency of each ion signal is determined and analyzed independently. At this processing stage, precise information for each ion including frequency, intensity, and m/z values is established. Step 3 (Fig. 1) determines the ion’s signal strength using a data plotting and analysis process that assesses the current induced by an ion on the Orbitrap detection electrodes as a function of acquisition time. For simplicity we call this signal strength determination the “STORI” process[10] with the slope of a STORI plot being proportional to the charge of the ion (Step 3, Fig. 1). In Step 4, the charge of the ion is determined by a slope-to-charge calibration function. The STORI slope of an ion with unknown charge is assigned the closest integer charge state on the calibration function, determined just once for each instrument (two different Q Exactive-style instruments were used in this study). Using the integer charge (z) and mass-to-charge ratio (m/z), it is possible to determine the mass of each ion and produce a spectrum in the true mass domain with quite different spectral properties and increased resolution via centroiding and binning single ion signals[8]. As a result of our charge assignment algorithm, further explained in the Online Methods, charge misassignments do occur at the 1–4% level of the total ion count in I2MS spectra; these artifacts are recognized readily as satellite peaks due to the ±1 error in charge state assignment. Further experimental optimization and software development is currently underway to completely eliminate charge misassignment. The I2MS process further uses validated raw data reduction and a reduced DC offset for the Orbitrap center electrode (−5 to −1 kilovolt) enabling greater than 70% of the ions to survive the ~2–4 second detection period (Supplemental Fig. 1). A more in-depth description of these steps is available in the Online Methods. We now focus on the application of the method to complex proteoform mixtures and detection of large complexes via native I2MS.

Initial experiments on a mixture of intact protein standards from 8–47 kDa highlighted an immediately apparent set of differences in I2MS readout from conventional ensemble mass spectrometry analysis. Basic figures of merit for the standards are shown in Supplemental Fig. 2c, illustrating high resolving power for intact proteins and the accurate determination of their monoisotopic mass. Note that I2MS produces simple spectra with an x-axis of mass, removing the normal requirement for the inference step to convert m/z to mass. The remarkably clean baselines of I2MS spectra proved highly insensitive to partially desolvated ions (Supplemental Fig. 2a vs. 2b), which extended the dynamic range for intact mass determination even when complexity was high (vide infra). The strict filtering criteria during STORI slope determination removed non-stable ion signals stemming from desolvation, ion fragmentation in the Orbitrap analyzer (such as neutral loss), and intermittent signals produced from electronics (electronic noise). In essence, I2MS transforms the instrument into an actual MASS spectrometer, removing the normal requirement for the inference step in conventional mass spectrometry.

We then used I2MS to analyze a mixture of heavy and light chains resulting from disulfide reduction of an IgG1 (Supplemental Fig. 3). Here, the removal of non-stable ion signals yielded highly sensitive I2MS spectra of species and adduct identification that was not possible in ensemble measurements. Ultimately, these mass domain readouts greatly simplify data analysis for non-experts.

To stress-test the I2MS process, it was compared to standard direct infusion with automatic gain control (AGC) set to 1×106 charges/spectrum for high complexity samples like the entire <30 kDa range of the human proteome from HEK293 T-Cells fractionated by GELFrEE technology[11]. A silver-stained SDS-PAGE gel (inset of Fig. 2a) shows the range of proteins that were analyzed by both standard MS and I2MS approaches (Fig. 2a vs. 2b). Although no protein masses could be determined using standard MS (Fig. 2a), I2MS (Fig. 2b) detected over 500 proteoform masses that were assigned using the intact mass tag approach[12, 13] and referencing prior top-down proteomics data where tandem MS (LC-MS/MS) was used to create a list of high quality proteoforms present in the Human Proteoform Atlas[14]. A comparable number of ions (~10 million) was utilized to generate the ensemble and I2MS spectra. While all identified proteoforms via I2MS analysis are listed in Supplemental Dataset 1, small insets in Fig. 2b illustrate 17 identified proteoforms that are listed in Supplemental Table 1. In addition, isotopic distributions for higher mass proteoforms between 20 – 25 kDa illustrated in the red box on Fig. 2b have not previously been identified and characterized by prior LC-MS/MS analysis. For conventional mass spectrometry, larger proteins (>20 kDa) are burdened with charge state distributions consisting of dozens of visible charge states. Appropriating a fixed number of ions between this many charge states lowers the S:N of each one. Consolidating the many low S:N charge states into one higher S:N mass distribution greatly increases the probability for their detection, annotation and characterization. I2MS, with the high intra-spectral dynamic range of ~2,000 clearly identifies these high mass features usually not detectable in lower dynamic range traditional m/z analysis. I2MS reduces charge state ambiguity and consolidates each charge state distribution into one mass channel, increasing detection capabilities for larger mass species in mixtures with extreme complexity.

Figure 2.

Deconvolution of Complex Proteoform Mixtures.

Comparison of the number of proteoform assignments possible using either conventional spectral acquisition with an ensemble of ions in a population (a) or the I2MS process (b). This comparative analysis involved directly infusing a mixture of 0 – 30 kDa proteins created by fractionation of whole extracts from HEK-293 human cells. The inset in panel (a) shows an analytical, silver-stained gel, with the fraction analyzed highlighted with a red rectangle. The 3 insets in panel (b) highlight 17 of the 550 total proteoforms assigned, which are listed in Supplemental Table 1 and include their proteoform record numbers (PFRs). The 20–25 kDa region highlighted with a red rectangle in panel (b) corresponds to proteoforms of higher mass proteins in the complex mixture previously unidentified using LC-MS/MS for top-down proteomics.

Another major advantage of I2MS is the ability to work with samples containing large protein complexes that are lifted into the gas phase using native electrospray ionization (so called native MS)[15]. Pressing I2MS further, we examined standard protein complexes[16, 17] and also applied it to engineered virus-like particles (VLP). As demonstrated by Supplemental Fig. 4, protein complexes including pyruvate kinase (Supplemental Fig. 4a, b, and c) and GroEL (Supplemental Fig. 4d, e, and f) with well resolved charge states validate I2MS implementation for native complexes. The charge state distributions from I2MS match those determined using standard MS, and the relative intensity of the charge assignments mirror their corresponding spectra in the m/z domain. Accurately assigned individual ions from resolved charge states produce well-defined peaks in the mass domain for pyruvate kinase (231.9 kDa) and GroEL (801.2 kDa). These higher mass assemblies are accompanied with spectral resolution values of 560 – 800. Further, minor species of the pyruvate kinase complex formed from loss of a C-terminal lysine in one or two of the subunits within the tetramer (Supplemental Fig. 4b) are easily recognized.

The single-ion analysis advanced here was conducted on two virus-like particles (VLPs) engineered from viral capsids carrying varying amounts of DNA and mRNA cargo loads after heterologous production in E. coli. Wild-type (WT) or mutant (MINI) versions of these VLPs were constructed from the MS2 wild-type coat protein or a Ser37Pro mutant, respectively. Recently, these particles have been investigated through electron microscopy and X-ray crystallography to determine capsid structure[18-20]. Briefly, WT and MINI VLPs are composed of exactly 180 and 60 copies of the 13.7 kDa coat protein (CP), generating theoretical masses of 2.47 and 0.82 MDa, and particles with 27 and 17 nm diameters. In turn, cavity volumes can be calculated around 10,300 and 2,600 nm3, assuming an approximate spherical nature of the VLPs. The WT and MINI monomer amino acid composition was verified through MS analysis in denatured mode (Supplemental Fig. 5). Although the structure of the VLPs is established, characterization of cargo such as therapeutic proteins and RNA would be highly valuable to quickly determine the stoichiometry and loading of drug-delivery within the particles.

To probe composition and cargo of the VLPs, I2MS data for the WT and MINI particles were compared to each other and contrasted with data from standard MS data produced in the m/z domain (Fig. 3c,d and Fig. 3e,f). No visible charge states were present in the m/z domain spectra of either the WT (Fig. 3c) or MINI (Fig. 3d) species. The lack of distinct charge states in the native MS data was attributed to the molecular heterogeneity of varying lengths of DNA and mRNA[19]. Given this molecular heterogeneity and no distinct charge state peaks present in the m/z-domain spectra, no information regarding mass of the VLP plus cargo could be discerned. This is an often-encountered situation in obtaining size distributions of highly complex samples in mass spectrometry and the reason other methods like gel electrophoresis and light scattering are used to provide broad molecular weight estimations.

Figure 3.

Mass Determination of Large Native Complexes.

Crystal structure renderings (a,b) and mass spectra of assembled WT (c,e,g) vs. MINI (d,f,h) MS2 virus-like particles, with approximate masses and capsid diameters of 3.2 vs. 1 MDa and 27 vs.17 nm, respectively. In panels (c) and (d), the standard MS spectra plotted in m/z space results in a loss of charge state resolution due to heterogeneity of cargo inside the capsid of the virus-like particles. However, the corresponding I2MS spectra in panels (e) and (f) allow mass assignment to the distribution of particles even without resolution of the individual charge states (a normal requirement of conventional data produced by electrospray MS). The individual ions assigned to their underlying charge distributions are shown in panels (g) and (h). Mass and resolution (Res.) at full-width of the half-maximum (FWHM) of the peak heights are labeled beside each peak.

In contrast to standard native MS, the I2MS spectra produced mass distributions of the WT (Fig. 3e) and MINI (Fig. 3f) species with masses of 3190 ±38 kDa and 990 ±16 kDa (mean ± s.d.). Assigned charge states for the set of single ions were centered around +137 and +67 for WT (Fig. 3g) and MINI species (Fig. 3h) with the expected Gaussian-like distributions. In addition, consolidating charge states in the creation of the VLP I2MS mass spectra increased resolution greater than 2-fold when compared to their m/z spectra counterparts. Determined I2MS masses minus proposed capsid masses reveal WT and MINI cargo masses of approximately 0.72 ± 0.038 MDa and 0.17 ± 0.016 MDa, respectively. Although the length of the DNA base pairs and mRNA nucleotides that form the VLP cargo loads vary across a wide distribution of lengths, MINI VLPs having a much smaller volume incorporate shorter DNA and mRNA strands than their WT counterparts. DNA and mRNA lengths averaged 900 (±400) vs. 375 (±150) base pairs (bp) and 500 (±250) vs. 100 (±50) nucleotides (nt) for WT vs. MINI species[19]. The mass distribution for the WT particles determined via I2MS (Fig. 3e) is consistent with the encapsulation of 1 average length DNA (900 bp) and mRNA (500 nt) strand together in the VLP (~0.74 MDa). The WT mass distribution with a FWHM of only 75 kDa means that as a larger DNA strand is added as cargo, it is partnered with a smaller length mRNA component during production and self-assembly within the cytosol of E. coli. The mass distribution for the MINI particles determined via I2MS (Fig. 3f) is consistent with VLP encapsulation of either 1 DNA strand or multiple smaller mRNA strands as cargo. The incorporation of the average length DNA strand (375 bp) matches the higher mass tailing feature around 1.06 MDa. In turn, the smaller VLP volume can contain multiple (~5) of the smaller ~100 nt mRNA strands, accounting for the majority of the MINI mass distribution. The ratio, 4.2, of the cargo masses () corresponds to the ratio of cavity volumes, (4.0) further supporting assignment of cargo content based on I2MS data for these spherical WT and MINI VLPs. The measurement of complex masses through I2MS, unobtainable via traditional mass spectrometry, reveals pertinent information to their composition and cargo. As a result, I2MS analysis is highly applicable to VLP engineering, can be extended to other promising systems including Qβ virions, adeno-associated virus (AAV), and hepatitis B, and can yield important information for more complex systems including bacterial microcompartments which contain an assortment of coat protein sizes.

With increased dynamic range, resolution, clean baselines and intuitive data interpretation, I2MS shows increased value for determination of mass values particularly as mixture complexity increases. Being an Orbitrap-based method, I2MS provides a readily deployable technique for single molecular MS of denatured proteins and native complexes with masses ranging from 5 kDa to 3.5 MDa. The ability to accurately determine mass values directly from dilute and complex samples regardless of charge state overlap will greatly inform and expand the utility of mass spectrometric analysis of proteins and their assemblies in improved characterization of molecular complexity for both endogenous and synthetic biology.

Online Methods

Cell growth, lysis, and protein fractionation

Human embryonic kidney cells HEK-293 (ATCC CRL-1573) as reported within the Life Sciences Reporting Summary were grown on 15 cm plates in DMEM media (Life Technologies) supplemented with 10% fetal bovine serum (Life Technologies), containing 1% Penicillin/Streptavidin (Life Technologies), and they were incubated at 37°C temperature under 5% CO2 atmosphere. Triplicate of confluent plates containing 1 × 107 cells were treated with 100 nM rapamycin or a DMSO control for 8 hours and harvested with 0.05% trypsin (Life Technologies) for 5 min. For all I2MS experiments reported here, cells were grown at control conditions. Cells were pelleted at 300 x g for 10 min. and washed twice with PBS (Life Technologies). Cells were incubated with lysis buffer (1% w/v N-lauroylsarcosine, 20 mM Tris Base pH 7.5, 100 mM NaCl, and HALT Inhibitor Cocktail (Thermo Fisher Scientific) for 10 min. on ice, 750 units of Benzonaze nuclease (Sigma) was added, the solution was incubated for 15 min. at 37°C, and the protein content was estimated using Pierce BCA assay (Thermo Fisher Scientific). Aliquots of 400 μg of protein were precipitated for 30 min. at −80°C with 4 volumes of cold acetone, resuspended in GELFrEE Tris-Acetate Sample Buffer (Expedeon), and proteins were fractionated on a GELFrEE 8100 (Expedeon) system using an 8% cartridge according to manufacture procedures. Only the first fraction containing proteins up to 30 kDa was collected.

Denatured mode top-down proteomics using LC-MS/MS

Fractions were methanol/chloroform/water precipitated[11], resuspended in buffer A (5% acetonitrile, 94.8% water, and 0.2% formic acid), and submitted to LC-MS/MS. Samples were fractionated on a Dionex UltiMate 3000 LC system (Thermo Fisher Scientific) using an in-house packed PLRP-S (5 μm particle, 1000 Å pore size) (Agilent) column (75 μm ID × 25 cm length) after loading on a PLRP-S trap column (150 μm ID × 2 cm). The loading pump was operated at 3 μL/min., and the nano pump gradient flow rate was set at 300 nL/min. Buffer A and buffer B (5% water, 94.8% acetonitrile, and 0.2% formic acid) were used for LC separation under the following gradient: 5% B at 0 min., 15% B at 5 min., 55% B at 55 min., 95% B from 58 to 61 min., 5% B from 64 to 80 min. Column was kept at constant temperature (55°C), the outlet was online with a 15 μm i.d. electrospray emitter (New Objective) packed with 3 mm of PLRP-S resin, and attached to a nanoelectrospray ionization source built in-house. LC was online with an Orbitrap Fusion Lumos (Thermo Fisher Scientific) mass spectrometer operating in “protein mode” with 2 mTorr of N2 pressure. Transfer capillary temperature was set at 320 °C, ion funnel RF was set at 30%, and a 15V of source CID was applied. MS[1] spectra were acquired at 120,000 resolving power (at 200 m/z), AGC target value of 500,000 charges/acquisition, 200 msec. of maximum injection time, and 4 μscans. Data-dependent top two MS[2] method using 23 NCE for higher-energy collisional dissociation (HCD) was used to generate MS[2] spectra that were acquired at 60,000 resolving power @ 200 m/z, with target AGC values of 500,000 charges/acquisition, 800 msec. of maximum injection time, and 4 μscans. Precursors were quadrupole isolated using a 3 m/z isolation window, dynamic exclusion of 60 sec. duration, and threshold of 2×104 intensity.

Protein Identification

Raw files from LC-MS/MS experiments were searched against a UniProt Homo sapiens (Taxon ID: 9606, Proteome ID: UP000005640) protein database using TDPortal (http://galaxy.kelleher.northwestern.edu/), as previously described[21, 22]. Proteoform search space was generated in silico allowing up to 11 PTMs or sequence variations per candidate proteoform. Two different search types were performed in parallel: Absolute Mass search with precursor tolerance of 2.2 Da and fragments tolerance of 10 ppm, and a Biomarker search with 10 ppm precursor and fragment tolerance. Estimation of false discovery rate (FDR) at the protein entry and proteoform level was performed[23]. Proteoforms above a 10% FDR threshold were reported and used for GELFrEE fraction I2MS spectrum annotation.

Proteoforms passing the FDR cutoff were transformed into theoretical +1 (M+H) isotopic distributions[23]. The I2MS mass spectrum was also converted into singly-charged data for comparison. With the identified proteoforms serving as theoretical distributions, the I2MS distributions were matched against the theoretical using an isotope fitting routine with a 25 ppm tolerance[24, 25]. Matches were loaded into TDValidator for visualization and manual validation (Proteinaceous, Inc.).

Virus-like Particle Expression and Purification

Glycerol stocks of E. coli DH10B cells containing plasmids encoding for either the MS2 wild-type coat protein or the mutant (mini) variant were grown in 10 mL 2xYT overnight at 37°C. These cells were then subcultured into 1 L 2xYT to a starting OD600 of 0.05 and grown with shaking at 37°C. When the OD600 reached 0.5, protein expression was induced with 0.1% arabinose. Expression continued overnight for ~ 16 h. After harvesting the cells by centrifugation, the cellular protein was extracted and VLPs were purified according to standard procedures using ammonium sulfate precipitation and fast protein liquid chromatography (FPLC)[20].

Additional Sample Preparation and MS Analysis

To produce the 4 protein mixture (Supplemental Fig. 2) ubiquitin (bovine; ID: P0CH28), myoglobin (equine heart; ID: P68082), carbonic anhydrase (bovine; ID: P00921), and enolase (baker’s yeast; ID: P00924) purchased from Sigma Aldrich were added together to produce a 2 mL solution with final protein concentrations of 20, 40, 50, and 90 nM, respectively. To produce the reduced version of the NIST antibody standard (Supplemental Fig. 3) 30 μL (300 μg) of stock solution (NIST) was mixed with 300 μL of 6 M GdHCl. 30 mM TCEP-HCl was added, the solution was incubated at 37°C for 1.5 h, and finally diluted to a 500 nM concentration. HEK293 GELFrEE fraction 1 (Fig. 1) was methanol/chloroform/water precipitated[11], resuspended in HPLC grade water and run through a HiPPR (Thermo Fisher Scientific) detergent removal column, following the manufacture’s protocol. The fraction protein concentration was determined with a BSA protocol (Thermo Fisher Scientfic) and diluted to a 5 nM concentration. Before dilution all samples were desalted six times at 10,000 x g for 5 min. in 3 kDa Amicon Ultra centrifugal filters (Merck Milipore) with 100 mM ammonium acetate. Once desalted, all samples were diluted and electrosprayed under denaturing conditions in a 40% acetonitrile and 0.2% acetic acid solution.

Pyruvate kinase and GroEL protein complexes (Supplemental Fig. 4) were purchased from Sigma Aldrich and diluted to final concentrations of 1 μM. GroEL was precipitated as previously described[26]. Both protein complexes were desalted six times at 10,000 x g for 4 min. in 100 kDa Amicon Ultra centrifugal filters (Merck Milipore) wtih 100 mM ammonium acetate. Once desalted, both samples were diluted and electrosprayed under native conditions in 100 mM ammonium acetate.

Denatured samples with molecular weights <100 kDa were introduced into a Q Exactive Plus (Thermo Fisher Scientific) mass spectrometer with a custom nano-electrospray source as described previously[27] using +0.8 to 1.6 kV spray voltage. Due to significant ion decay rates observed for denatured proteins with molecular weights between 8–150 kDa, the Q Exactive Plus mass spectrometer was modified to reduce the central electrode voltage of the Orbitrap to −1 kV during ion detection. Samples >200 kDa ionized with native-mode electrospray were introduced into a Q Exactive Ultra High Mass Range (UHMR) (Thermo Fisher Scientific) mass spectrometer with a Nanospray Flex Ion Source (Thermo Fisher Scientific) with spray voltages between +1.4 to 2 kV. As large native complexes/VLPs are large relative to the size of the background gas molecules, collisions are ineffective at either reducing the momentum of the species or transferring enough internal energy to cause dissociation. As a result, less decay for these types of ions were observed within the Orbitrap analyzer and thus the central electrode voltage for the Q Exactive UHMR was kept at −5 kV for ion detection. For both instruments, acquisition rates were reduced down to 0.5 spectra/sec, corresponding to a 2 second detection period per acquisition event and the HCD-pressure was set between 0 to 0.5 (UHV Pressure < 5×10−11 torr) to reduce ion and background-gas collisions from a more typical setting of 1 (UHV Pressure ~9×10−11 torr). Enhanced Fourier-Transform (eFT) was turned off, as it was not necessary to process single ion signals in this manner. A source temperature of 320°C, an in-source collision-induced dissociation (SID) value of 0 to 15 V, and an in-source trapping (UHMR) value of −50 to −150 V was optimized on a per species basis.

Ion Collection and Data Acquisition of the I2MS Process

The overall scheme describing I2MS is shown in Fig. 1. Below is an extended technical description of this new process.

Ion Collection (Step 1)

Similar to traditional CDMS techniques, I2MS evaluates single ion signals to produce a true mass spectrum. However, Orbitrap CDMS analysis differs from traditional linear ion trap CDMS in the efficiency of single ion collection[2-5]. With the implementation of image current detectors our single ion definition is not confined to one ion per acquisition event, but one ion at a defined frequency per acquisition, meaning that multiple ions are analyzed per acquisition as long as they do not correspond to multiple ions at the same m/z value. To lower the number of ions entering the Orbitrap analyzer, denatured samples were diluted to concentrations between 5–50 nM, while samples of large complexes were analyzed by native electrospray in the 1–3 μM regime. In addition, Automatic Gain Control (AGC) was disabled and ion injection times were set between 0.03 – 1 ms to manually attenuate signal to minimize multi-ion events and maximize the collection individual ions. As this is a random style trapping event, the exact number of individual ion signals collected varies slightly for each acquisition event. In some cases, the ion optics including the C-trap entrance lens were detuned to further reduce the number of ions entering the Orbitrap analyzer.

Individual Ion STORI Slope and Centroid m/z Determination (Fig. 1, Steps 2 and 3)

Within each acquisition event, every detected ion was analyzed individually. Both the mass-to-charge ratio (m/z) and charge (z) were necessary to determine the mass of the ion in question. In FT-MS Orbitrap analysis, m/z was determined from the frequency of ion rotation around the central electrode. To determine the centroid m/z value, the apex of the profile peak was determined by a quadratic fit to the three most intense m/z points. The z was measured from the rate of the induced charge on the Orbitrap outer electrode, otherwise known as Selective Temporal Overview of Resonant Ions (STORI) as previously described[10]. Briefly, STORI analysis was the integration of the ion induced charge over the course of the time domain acquisition period at the specified single ion frequency value. The slope of the STORI plot remained constant until either the end of the acquisition event or premature ion loss within the Orbitrap analyzer. The STORI plot slope of a particular ion was proportional to the charge of the observed ion. Supplementary Fig. 6 demonstrates the dependence depicting increasing STORI slopes for the +13 charge state of ubiquitin, the +25 charge state of myoglobin, and the +39 charge state of carbonic anhydrase. STORI slope is only dependent on the charge of the ion and not the m/z or mass of the species in question. The charge detection limit for I2MS analysis is such that ions only with a charge of +6 or greater can be measured. Ion signals with a charge state below that threshold cannot be distinguished from electronic noise. Electronic noise in FT decreases as a function of and longer acquisition time or hardware modifications will be necessary to further decrease this charge detection limit.

The variability in STORI slope decreases proportionately to the . To accurately determine the charge of an ion from its STORI slope, long acquisition times are necessary for which the ion signal must persist (ideally for the full 2 s detection period). Many factors contribute to how long an ion signal persists. A few include the amount of background gas that enters the Ultra-high vacuum (UHV) region of the mass spectrometer, velocity of the ions, length of the acquisition, and finally the collisional cross-section of the ion. Although we have successfully increased ion signal persistence through the reduction of ion kinetic energy (decreasing ion velocity) and limiting the amount of background gas entering the UHV region of the mass spectrometer, reference 8 gives a brief overview of how collisional cross-section can impact individual ion decay for denatured protein ions.

Charge Determination and Slope Averaging Technique (Fig. 1, Step 4)

In order to accurately assign the unknown z of an ion with a known STORI slope value, each ion slope was calibrated to a charge calibration function. Supplemental Fig. 7 illustrates 2 linear calibration functions created from isolated ions with known charge states analyzed with a Q Exactive Plus (red points/trace) and Q Exactive UHMR (blue points/trace). To produce the trace for the Q Exactive Plus various known charge states from denatured ubiquitin, myoglobin, carbonic anhydrase, and enolase were analyzed. To produce the trace for the Q Exactive UHMR various known charge states from natively electrosprayed carbonic anhydrase, monoclonal NIST antibody standard, pyruvate kinase, and groEL were analyzed. Thousands of single ions from each known charge state were analyzed and their median STORI slope value was calculated (represented by each black point). Each calibration file was manually curated to ensure the STORI slopes utilized to calculate the average value did not include slopes from two ion signals at the same m/z value or slopes from contaminating species that overlap in m/z space. The points on the calibration were fit with a linear regression yielding the equation below with a strong statistical relationship (R2 = 0.9997):

It is good to note that the charge calibration curve was specific to each instrument. Small differences in ion trajectories into the Orbitrap and differing imperfections of the Orbitrap analyzer for each instrument affect the slope of the calibration function. However, once the calibration was produced, the instrument required no recalibration as long as the injection energy of the ions from the C-Trap into the Orbitrap analyzer was not changed. To support this conclusion, testing on the Q Exactive UHMR over a 3-month period revealed, on average, only a 0.1% drift in the median charge state assigned using initial STORI plots. This negligible shift in the calibration points over time demonstrates the longitudinally robust nature of the detection process.

To determine the quantized z with a >96% rate of correct assignment, the ion STORI slope was scored based on proximity to the nearest median charge calibration value and assigned. Without any statistical processing methods, current instrumentation limits the correct charge assignment of ions to approximately 50%. The other half of ion signals are mis-assigned to a charge state that either high or low by one integer value. Supplemental Fig. 8 shows the charge assignment of approximately 800 single ion signals for the isolated myoglobin +20 charge state collected over a 2 second transient. Although all the ion signals collected were known to correspond to a +20 charge, slight variations in their STORI slopes resulted in the incorrect charge assignment of 47.5% of the ion signals collected, with the incorrect charge assignments having a Gaussian-like distribution to higher and lower charge states. A slope averaging mechanism, more commonly known as the central limit theorem, was implemented to correct the misassignments.

When STORI slope values from a known charge state were averaged, there was a reduction in charge misassignment proportional to , where n is equal to the number of ion STORI slope values being averaged. Supplemental Fig. 8 illustrates the implementation of this averaging mechanism, where the +20 charge assignment for the myoglobin species improved to 81.1% and 96% when averaging 4 and 16 ion STORI slope values, respectively, and subsequently assigning charge. This averaging of 16 STORI slopes was cemented into the algorithm for data acquisition utilized in the creation of all denatured I2MS spectra in this report to accurately assign charge to individual analyte ions from 6 – 100 kDa carrying +6 to +80 charge states in the 500 – 5,000 m/z range. Averaged STORI slope counts increased up to 100 to assign >96% of individual ions to the correct charge state for native species in the 10,000 – 23,000 m/z range. Before averaging STORI slope values, first all ions collected were ordered as a function of their m/z values. A slope tolerance for STORI slope values 30% above or below the median slope value of the points being averaged in addition to a maximum 5 m/z window for ordered STORI slopes to be averaged provided strict guidelines to remove outlier values. STORI slope values are averaged pre-charge calibration and the resulting calculated z was parsed back to all input single ions. Once an accurate quantized charge was calculated for each single ion the mass was determined:

Extensive testing revealed that for the Orbitrap detector, with a pre-amplifier lacking a high-pass filter, there was no deviation from the linearity of the linear function that converts STORI slope to charge for either the Q Exactive Plus or Q Exactive UHMR for ions with differing m/z or mass values. This coupled with the high correlation of the calibration functions R2=0.9997 and R2=0.9999, respectively, showed that extrapolation to higher non-surveyed charge states for VLP samples was reliably accomplished. To further validate this assertion, the collection of 2 GroEL +70 ion signals at the same frequency value within one acquisition, accurately yields a STORI slope that corresponds to a single +140 (+70 × 2) ion signal.

Plotting of I2MS Spectra (Fig. 1, Step 5)

Once the mass of each ion was calculated, the thousands of acquired individual ions were binned to create the I2MS spectrum. Histogram bin sizes varied based on the species analyzed. Masses under 200 kDa ions were binned in 0.2 Da increments. As masses increased beyond 200 and 1,000 kDa and resolution decreased, bin sizes were increased to 30 and 500 Da, respectively. In tandem, the calculated charge of each ion utilized to create the I2MS spectrum was binned in quantized domains to validate charge assignment. A simple protein example demonstrating these outputs is shown in Supplemental Fig. 9. Supplemental Fig. 9a,b,c represent the traditional m/z mass spectrum of myoglobin, corresponding to the I2MS spectrum in the mass domain, and single ion charge assignment histogram for ions utilized to create the I2MS spectrum, respectively. Supplemental Fig. 9a mirrors the charge state distribution in Supplemental Fig. 9c.

Statistical Analyses

Graph visualizations and calculations were preformed using Thermo Xcalibur Qual Browser software, Microsoft Excel 2016 software, or custom compiled C# code. Although experiments were generally performed twice, thousands of individual scan events were utilized to produce I2MS spectra with statistically meaningful mass distributions. These experimental mass distributions were then verified with previous datasets or matched to theoretical mass distributions. Furthermore, the method as a whole was anchored to a set of standards to validate I2MS as an accurate and reproducible platform.

Data Availability

Datasets utilized for the I2MS analyses can be found on the MassiVE repository, MSV000083840 doi:10.25345/C5XD01. Other data that support the findings of this study are available from the corresponding authors upon request.

Code Availability

Custom code associated with the I2MS creation process is available from the corresponding authors upon reasonable request.