Subscriber access provided by YORK UNIV
Article
Graphical Ranking of Divisors to Get the Most out of a Resolution-Enhanced Kendrick Mass Defect Plot Sayaka Nakamura, Robert Bernard Cody, Hiroaki Sato, and Thierry Nicolas Jean Fouquet Anal. Chem., Just Accepted Manuscript • DOI: 10.1021/acs.analchem.8b04371 • Publication Date (Web): 24 Dec 2018 Downloaded from http://pubs.acs.org on December 24, 2018
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 18 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Analytical Chemistry
1
Graphical Ranking of Divisors to Get the Most out of a Resolution-Enhanced
2
Kendrick Mass Defect Plot
3
Sayaka Nakamura,1 Robert B. Cody,2* Hiroaki Sato,1 Thierry Fouquet1*
4
1
5
Technology (AIST), Tsukuba, Ibaraki, Japan. 2 JEOL USA Inc., Peabody, MA, USA.
Research Institute for Sustainable Chemistry, National Institute of Advanced Industrial Science and
6 7 8
Abstract: Resolution-enhanced Kendrick mass defect (KMD) analysis using the new concept of
9
fractional base units (repeating unit R divided by an integer X, R/X as a mathematical moiety) is now
10
a powerful data processing tool to unravel complex mass spectra of polymers. It enhances a regular
11
KMD analysis using the chemical moiety R to compute mass defects with unprecedented separation of
12
ion series differing by their isotopic or co-monomeric content, end-groups or charge states in highly
13
visual KMD plots. The value of the divisor X dictates the gain of separating power from the regular to
14
the resolution-enhanced KMD plot and its choice strongly affects the ease and speed of data
15
interpretation. A simple tool to help selecting the best values of X depending on the users’ needs is
16
mandatory to rationalize the analysis and avoid a time-consuming trial-and-error methodology. We
17
propose two graphical representations intuitively ranking the well-suited divisors for the appropriate
18
separation of isotopes and/or co-oligomers for copolymeric mass spectral data. Rankings are extended
19
to any type of dataset from homopolymeric blends to terpolymers by generalizing the formulas with
20
three variables beyond the specific separation of isotopes. The “RANK” functions are now available in
21
commercial or home-made spreadsheet (available upon request) to interactively select divisors and
22
compute the associated KMD plots.
23 24 25
Table of Contents
26 27
ACS Paragon Plus Environment
Analytical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
28
Page 2 of 18
Introduction
29
It has been a long journey from the change of basis proposed by Kendrick fifty years ago1 to the
30
compositional mapping of complex mass spectra from hydrocarbons2-4 and polymers5-10 using Kendrick
31
mass defect (KMD) plots. A “regular” KMD analysis first consists in peak-picking a high-resolution /
32
high-accuracy mass spectrum displaying peaks spaced by one or several repetitive moieties. Mass-to-
33
charge ratios (m/z) are then turned into Kendrick masses by an affine transformation
34
KM(R) = 𝑚/𝑧 ∙
round(R) R
(1)
35
with R the exact mass of a repeating unit used as “base unit”. The associated “defects” are evaluated
36
using
37
(2)
KMD = round(KM) ―KM
38
Plotting KMDs as a function of m/z produces visual and informative “maps”. Horizontally aligned
39
points correspond to ions of identical elemental composition modulo R while oblique alignments reveal
40
variations of elemental compositions in any moiety but R.3-8 All these steps may be done using
41
commercial programs dedicated to the KMD analysis or spreadsheets. It makes KMD analysis a very
42
simple data processing tool available to everyone equipped with a computer.
43
Analysts have taken a renewed interest in KMD analysis11 thanks to the recent introduction of
44
“resolution-enhanced” KMDs computed by replacing the “chemical” moiety R with a “mathematical”
45
moiety R/X, a fractional base unit with X being a positive integer.12 R
46
KM(R,X) = m/z ∙
round(X) R X
(3)
47
The associated KMD plots variably separate ion series differing by their end-groups,13 co-
48
monomeric content14 or charge states15 depending on the value of X, taking full advantage of the spectral
49
width of a KMD plot. Resolution-enhanced KMDs are also compatible with low-accuracy data as the
50
decrease of accuracy is counterbalanced by the improved separating power of the plot.12 In spite of its
51
undeniable capabilities for polymer and carbonaceous ions,16 the question remains: how to choose the
52
optimal divisor X of the fractional base unit R/X?
53
In a previous letter to the editor, we demonstrated that a “recommended range” of divisors X falls
54
between round(2/3R) and round(2R).17 For the basic methylene moiety R = CH2, the range of
55
recommended divisors is X = {9, …, 28} (= 20 divisors) which means twenty resolution-enhanced
56
KMD plots should be drawn to find the most favorable map. The range becomes larger with 60 divisors
57
and plots for a poly(ethylene glycol) backbone (R = 44.0262, X = {29, ..., 88}) up to 534 divisors and
58
plots for a poly(phenylsulfone) sample (R = 400.0768, X = {267, …, 800}) or more. A systematic
59
computation of all the plots would cancel the benefit of KMD analysis as a high-speed data processing
60
tool.
61
To reduce the trial-and-error to its minimum, a user would advantageously select a divisor from a
62
short list to instantly compute the KMD plots depending on his/her expectations. A resolution-enhanced
ACS Paragon Plus Environment
Page 3 of 18 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Analytical Chemistry
63
KMD analysis of homo-, co- or terpolymers (one, two or three repeating units within the same
64
backbone) may indeed separate ion series depending on their isotopic composition and/or end-groups
65
and/or co-monomeric content. A first tool to choose the well-fitted divisors X has been proposed
66
elsewhere17 based on the separation of isotopes but valid for homopolymers only. A more convenient
67
and visual approach should be developed that can extend beyond homopolymers. This article describes
68
the concepts of two “ranking function” named RANK1 and RANK2 intuitively positioning the optimal
69
divisors in graphs for a resolution-enhanced KMD analysis of copolymers separating the ion series by
70
one co-monomer and/or the isotopic content. These functions are generalized for the selection of the
71
best divisors for a resolution-enhanced KMD analysis separating ions by any moiety such as end-groups
72
for homopolymers or three repeating units for the notoriously complex case of terpolymers.
73 74
Experimental
75
Polymer samples
76
(H, OH)-ended poly(ethylene oxide-block-propylene oxide-block-ethylene oxide) triblock
77
copolymer8,18 (P(EO-b-PO-b-EO), repeating units: EO = 44.0262 and PO = 58.0419, ~10 wt% EO,
78
Mn~1100 g mol-1) was from Sigma Aldrich (St Louis, US). PEO “monostearate” was purchased from
79
Wako Pure Chemicals Industries (Osaka, Japan).13 A terpolyester formed upon the polymerization of
80
adipic acid condensed on 1,4-butanediol, 2,2-dimethyl-1,3-propanediol and 3-methyl-1,5-pentanediol
81
(repeating units: C10H16O4, C11H18O4 and C12H20O4, W-2300) was kindly supplied by a Japanese
82
company (name upon request).
83 84
Mass spectrometry
85
High-resolution mass spectra were recorded using a JMS-S3000 MALDI SpiralTOFTM mass
86
spectrometer (JEOL ,Tokyo, Japan)19 with a solvent-free sample preparation20 using {(2E)-2-methyl-3-
87
[4-(2-methyl-2-propanyl)phenyl]-2-propen-1-ylidene}malononitrile aka DCTB as matrix (Tokyo
88
Chemical Industry, Japan) and sodium trifluoroacetate NaTFA as cationizing agent (Sigma Aldrich).
89
MSTornado control (JEOL) was used for data recording, mMass 5.521 for data treatment and Illustrator
90
CC (Adobe Inc., San Jose, CA, USA) for artworks.
91 92
Kendrick mass defect analysis
93
Points were peak-picked from the mass spectra with or without deisotoping using mMass (maximum
94
charge state: 1, isotope mass tolerance: 0.01 m/z) and a relative intensity threshold set at 0.5%. KMD
95
plots were independently computed using either msRepeatFinder v.2 (JEOL), Mass Mountaineer™
96
(massmountaineer.com) and a home-made “Kendo” Excel sheet (available upon request). A KMD plot
97
is a bubble chart displaying KMDs (y-axis) as a function of m/z (x-axis) and intensities (via the size of
98
disks). DP plots and ranking functions RANK1 and RANK2 were generated using either Mass
99
Mountaineer or “Kendo”. A DP plot (DP: degree of polymerization) is a bubble chart displaying the
ACS Paragon Plus Environment
Analytical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 4 of 18
100
referenced KMDs strictly identical to the discrete content in one or the other repeating unit using the
101
sum of the masses of the end-groups and adducted ion as reference.22
102 103
Results and discussion
104
Invalidity of the default divisors for copolymers: RANK1
105
The variation of KMD replacing one
12C
carbon atom by its first
13C
isotope in the elemental
106
composition of any ion (∆m = ―12 + 13.003355 = +1.003355) using R/X as the base unit is noted
107
KMD ― 12C + 13C(R,X). Within the “recommended range” of divisors, KMD ― 12C + 13C(R,X) varies linearly
108
with the minimal step per iteration of X.17 The minimum of its absolute value is reached for X =
109
round(R) so the separating power for the 13C isotopes of the associated KMD plot is minimal. It strictly
110
corresponds to the regular KMD plot computed with R as the base unit (Equations S1-S4) and accounts
111
for the condensed appearance of all the plots in the literature using a chemical moiety as the base unit.
112
Decreasing or increasing X iteratively from X = round(R), the separating power for the 13C isotopes of
113
the associated resolution-enhanced KMD plots is increasing stepwise. From these results compiled with
114
a PEO homopolymer standard,17 a default divisor X = round(R)-1 was proposed for the resolution-
115
enhanced KMD analysis with the fractional base unit R/(round(R)-1) extensively used in practice.11-14
116
The regular KMD plot from a deisotoped mass spectrum of a P(EO-b-PO-b-EO) triblock copolymer
117
(EO or EO/44 as the base unit) is depicted in Figure 1A. It displays a typical condensed and unresolved
118
cloud of points covering a narrow 0.1 range offering the worst possible visualization of data. The
119
resolution-enhanced KMD plot using X = 43 as the default divisor (round(EO) = 44, X = round(EO)-1
120
= 43) is depicted in Figure 1B. If EO/PO co-oligomers are greatly separated depending on their content
121
in PO (one line per fixed composition in PO and a varying content in EO, e.g. EOxPO14-16 highlighted
122
with grey rectangles), the gain of separation is too great. EOxPO16 co-oligomers indeed appear between
123
EOxPO14 and EOxPO15 series due to a severe aliasing22 (blue, red, green and yellow points being aliased
124
once, twice, three and four times) making the plot poorly intuitive. A manual anti-aliasing can be done
125
once all the ion series have been assigned by adding an integer N = {0,1,2…} to aliased KMDs. The
126
anti-aliased resolution-enhanced KMD plot using X = 43 as the default divisor is depicted is Figure 1C
127
with a spectral width of 4 instead of 122 (i.e. N = 0-4,) but displaying EO/PO co-oligomers in the correct
128
order. The whole procedure remains impractical and requires a complete assignment which ruins the
129
benefit of a KMD analysis. In other words, the default divisor mildly expanding the 13C isotopes of a
130
homopolymer is no longer valid when considering the separating power of a resolution-enhanced KMD
131
plot for a generic R2 repeating unit from a P(R1-co-R2) copolymer.
ACS Paragon Plus Environment
Page 5 of 18 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Analytical Chemistry
132 133
Figure 1. (A) Regular KMD plot (deisotoped sodiated (H,OH)-ended P(EO-b-PO-b-EO), base unit:
134
EO or EO/44). (B) Aliased resolution-enhanced KMD plot using EO/43 as base unit. (C) Anti-
135
aliased plot with a spectral width at ~4. (D) Variation of KMDs with PO using EO/X as base unit.
136
(E) RANK1(EO,PO). Resolution-enhanced KMD plot using (F) EO/47 and (G) EO/50 as the base
137
units with no aliasing.
138 139
Instead of KMD ― 12C + 13C(R,X), KMD +R2(R1,X) is a more relevant function dictating the variation
140
of KMDs while adding R2 to any generic elemental composition using R1/X as the base unit (Figure
141
1D with R1 = EO and R2 = PO). As for KMD ― 12C + 13C(EO,X), KMD+PO(EO,X) is either erratic or null
142
for X below round(2/3EO) or above round(2EO) while it varies linearly with the smallest step at each
143
iteration of X comprised between round(2/3EO) and round(2EO). It validates the generic recommended
144
range of divisors regardless of R and the variation of elemental composition. With the exact mass of
145
one repeating unit as the only input required for the data processing, a list of all the possible divisors
ACS Paragon Plus Environment
Analytical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 6 of 18
146
and their associated variation of KMD is predetermined as a very first step towards the selection of a
147
few well-suited divisors.
148
Identifying the optimal divisors to avoid aliasing for the resolution-enhanced KMD plot is not
149
obvious from the graph of KMD+PO(EO,X). Instead, a first ranking function RANK1 is proposed being
150
simply the absolute values of KMD+PO(EO,X):
151
RANK1(R1,R2) = |KMD +R2(R1,X)|
(4)
152
RANK1 is a function of the exact mass of the repeating units R1 and R2 only while X falls in the
153
recommended range itself a function of R1 (X ∈ round 3R1 ,…,round(2R1) ). RANK1 tends towards
154
0.5 for the highest value of KMD+R2(R1,X) and tends to 0 for the lowest value of KMD+R2(R1,X).
155
Minima are thus associated with the divisors X that result in minimal separation for R2 in the resolution-
156
enhanced KMD plot computed with R1/X as the base unit. A user can then easily infer how a KMD
157
plot would change depending on the value of the divisor X simply by looking at the position of the
158
points in the RANK1 graph.
{
(2 )
}
159
RANK1(EO,PO) is depicted in Figure 1E with a logarithmic scale for the y-axis to emphasize the
160
global and local minima (Figure S1 for the regular scale). The global minimum is reached for X = 44
161
corresponding to the regular KMD analysis. Choosing a value of X associated with a local minimum of
162
RANK1 improves the separating power of the plot as reflected by the position of the points, the higher
163
the minimum the higher the separation until aliasing occurs. The resolution-enhanced KMD plots
164
computed with EO/47 (Figure 1F) and EO/50 (Figure 1G) separate the EO/PO co-oligomers to a higher
165
and higher extent with the clusters of points covering a ~0.5 and ~1 range on the verge of aliasing but
166
maximizing the advantage of the resolution-enhanced analysis. For X = 50, it resembles the anti-aliased
167
KMD plot computed with EO/43 (Figure 1C) with EO/PO co-oligomers separated in an intuitive order
168
PO14 < PO15 < PO16 but no anti-aliasing is needed. Of note, aliasing may occur with local minima but it
169
is not possible to foresee to what extent as it depends on the number of ion series (e.g. end-groups) and
170
the number of lines per ion series (e.g. isotopes) which are initially unknown. The graphical ranking of
171
the divisors using RANK1 and their subsequent selection represents nevertheless a great savings in time
172
as compared to the systematic computation of all the 60 resolution-enhanced KMD plots (Figure S2)
173
for a selection of the best map a posteriori.
174 175
Application 1: resolution-enhanced DP plot
176
A direct application consists in the easy computation of “DP plots” (DP: degree of polymerization)
177
using the “referenced KMD” concept.22 Knowing the nature of the end-groups and the adduct ion of a
178
given copolymeric distribution, a DP plot displays the discrete composition in the co-monomers as y-
179
and x-axes in lieu of the abstruse KMD and m/z. The regular DP plot computed from the mass spectrum
180
of the sodiated (H,OH)-ended P(EO-b-PO-b-EO) copolymer after automatic deisotoping using EO (or
181
EO/44) and PO (or PO/58) as base units and (H2O, Na) as reference is depicted in Figure 2A. The
ACS Paragon Plus Environment
Page 7 of 18 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Analytical Chemistry
182
content in EO (from 0 to ~14) and PO (from ~10 to 24) is barely evaluated but the plot remains fuzzy
183
with poor point alignment due to the low separation power of the regular KMD plots and the limited
184
accuracy of the mass measurements.
185 186
Figure 2. (A) Regular DP plot (deisotoped sodiated (H,OH)-ended P(EO-b-PO-b-EO), base units:
187
EO/44 and PO/58). (B) RANK1(PO,EO). (C) Resolution-enhanced DP plot using EO/47 and PO/54
188
as the base units.
189 190
A resolution-enhanced DP plot using the default divisors X = 43 for EO (round(EO) = 44, X =
191
round(R)-1 = 43) and X = 57 for PO (round(PO) = 58, X = round(R)-1 = 57) exhibits a strong aliasing
192
due to an over-expansion of the KMD dimension highlighted by non-integer DPs (Figure S3).
193
Automatic anti-aliasing is possible, but it complicates the data treatment with additional parameters
194
(minimal and maximal DPs, tolerance).22 Instead, well-suited divisors are readily selected from RANK1
195
computed with EO (Figure 1C) and PO as base units (Figure 2B). Local minima associated with X =
196
47 for EO and X = 54 for PO lead to a limited expansion of the KMD dimension. The resulting DP plot
197
using EO/47 and PO/54 is depicted in Figure 2C with no aliasing but a striking re-alignment of points
198
along integer values as the discrete content in EO and PO.
199 200
RANK1 with three variables: R1, R2 and 13C
ACS Paragon Plus Environment
Analytical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 8 of 18
201
In the first section, only one variation of KMDs has been considered – either KMD ― 12C + 13C(R,X) or
202
KMD+R2(R1,X) – but their combination would be of interest for a controlled gain of separating power
203
depending on the content in 13C and R2. Figure 3A displays KMD ― 12C + 13C(EO,X) (red dots) together with
204
KMD+PO(EO,X) (black dots) with the range of divisors X = {round(2/3EO),…,round(2EO)} being
205
common to both functions. It becomes less and less obvious to spot the well-suited divisors for the
206
resolution-enhanced KMD analysis.
207
Based on Equation 4, a generalized definition of RANK1 using three variables becomes: RANK1(R1,R2,13C) = |KMD +R2(R1,X)| + |KMD ― 12C + 13C(R,X)|
208
(5)
209
with R1 the repeating unit used as base unit for the regular KMD plot, R2 the other repeating unit and
210
X the divisor for the resolution-enhanced KMD plot using R1/X as base unit. From its definition, it is
211
obvious that if KMD +R2(R1,X) and KMD ― 12C + 13C(R,X) are small (low separating power in the KMD
212
plot for a change of elemental composition in R2 and 13C), then RANK1 is small. If KMD +R2(R1,X)
213
and/or KMD ― 12C + 13C(R,X) tends to 0.5 then RANK1 tends to 1. It thus ranks the divisors in terms of
214
variation of KMD for both R2 and 13C, with the lowest separating power for the global minimum and
215
higher separating powers for local minima (the higher the minimum the higher the separation by both
216
R2 and 13C until aliasing occurs). The associated plot using EO as R1, PO as R2 and 13C is depicted in
217
Figure 3B and displays a series of global (X=44, regular KMD analysis) and local minima (e.g. X = 88,
218
X = 41, X = 47…).
219
The regular KMD plot computed from the raw mass spectrum of the P(EO-b-PO-b-EO) triblock
220
copolymer (no deisotoping; generic composition EOxPOy and 12C / 13Cz isotopes forming the isotopic
221
distribution, 12C being the first peak of the pattern) using X = 44 (global minimum of RANK1) displays
222
a condensed and unresolved cluster of points (Figure S4). Zooming in is mandatory but errors of mass
223
measurement lead to misaligned points departing from the horizontal alignment.
224
From the third lowest minimum of RANK1 in Figure 3B associated with X = 41, the resolution-
225
enhanced KMD plot computed with EO/41 as the base unit now displays an expanded cluster spreading
226
over ~0.7 (Figure 3C). Clear horizontal alignments of points are assigned to a fixed composition in
227
12C/13C
228
rectangles, EOxPO15 and EOxPO16). It is an indisputable enhancement of data visualization, while the
229
speed of processing is considerably increased as compared to a systematic analysis of the 60 plots in
230
the hope of coming across a favourable situation (Figure S4).
isotopes (red rectangles in inset in Figure 3C, EOxPO16, 13C1, 13C2) or a fixed PO content (grey
ACS Paragon Plus Environment
Page 9 of 18 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Analytical Chemistry
231 232
Figure 3. (A) Variation of KMDs adding a PO repeating unit (black dots) or replacing 12C by 13C
233
(red dots) with EO/X as the base unit. (B) RANK1(EO,PO,13C). (C) Resolution-enhanced KMD plot
234
(raw sodiated (H,OH)-ended P(EO-b-PO-b-EO), base unit: EO/41).
235 236
Inverted separation by R2 and 13C: RANK2
237
Ranking the divisors associated with a high separating power for both R2 and 13C is not favored as
238
it will lead to extreme expansion of the KMD dimension for both changes of composition and severe
239
aliasing. The last possibility is a high separating power for a change of composition in R2 and a low
240
separating power of KMDs for the isotopic composition in
241
resolution-enhanced KMD plot would display a set of condensed groups of points, every group differing
242
from the other in its content in R2 while isotope patterns would be clustered within one group. An
243
intuitive function to rank the divisors X associated with a high value of KMD +R2(R1,X) and a low value
244
of KMD ― 12C + 13C(R1,X) or vice-versa may be |KMD +R2(R1,X)| ― |KMD ― 12C + 13C(R1,X)|. It remains
245
too simple as it is not sensitive to the relative gain of separation between R2 and
246
|KMD +R2(R1,X)|/|KMD ― 12C + 13C(R1,X)| should indeed be either maximal or minimal to ensure the
247
best separation of R2 relatively to 13C.
248
A refined function RANK2 is thus:
13C
ACS Paragon Plus Environment
– or vice versa. For example, a
13C.
The ratio
Analytical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 10 of 18
|KMD +R2(R1,X)| ― |KMD ― 12C + 13C(R1,X)|
RANK2(R1,R2,13C ) = |KMD
249
+R2(R1,X)|
(6)
+ |KMD ― 12C + 13C(R1,X)|
250
If |KMD +R2(R1,X)| and |KMD ― 12C + 13C(R1,X)| are low or high, then RANK2 tends towards 0. If
251
|KMD +R2(R1,X)| or |KMD ― 12C + 13C(R1,X)| is high then RANK2 tends towards +1 for
252
|KMD +R2(R1,X)|
253
|KMD ― 12C + 13C(R1,X)|. RANK2 thus reveals the best divisors for a high separating power by R2 and a
254
low separating power by 13C isotopes as global and local maxima, and vice versa as global and local
255
minima.
>
|KMD ― 12C + 13C(R1,X)|)
or
towards
-1
for
|KMD +R2(R1,X)|
250 > 164 > 193 > 136 > 236 > 178, the broader the ellipses).
379 380
Conclusion
381
Two functions have been developed to graphically rank the divisors depending on the user’s need
382
for a rationalized resolution-enhanced KMD analysis of polymer ions. Based on very simple formulas
383
and requiring no inputs other than the mass of the repeating units (i.e. the same input as for any KMD
384
analysis), they are now available in a commercial program and in our home-made “Kendo” Excel sheet
385
available upon request (Figure S10 for screenshots). For an ultimate user-friendly analysis, divisors X
386
are directly selected from RANK1 / RANK2 and the associated resolution-enhanced KMD plots
387
instantly computed. Other user-friendly selection tools are conceivable such as a scroll bar to adjust the
388
enhancement of separation in the KMD plot (i.e. choosing an expected KMD +R2(R1) for the program
389
to find the divisor X associated with the closest actual KMD +R2(R1,X)). It is nevertheless a matter of
390
programming and comfort of use only as any selection of the well-suited divisors will be based on either
391
KMD +R2(R1,X) or the RANKs functions.
392
The application to terpolymeric mass spectral data replacing the variation of KMDs for 13C isotopes
393
by a generic variation for a third repeating unit R3 will be invaluable for analysts facing complex
ACS Paragon Plus Environment
Analytical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
394
industrial formulations, turning the uninformative regular KMD plot of a terpolymer into a set of
395
copolymeric or homopolymeric clusters in the resolution-enhanced plots via a sequential analysis.
396 397
ASSOCIATED CONTENT
398
Supporting Information
399
Identity of KM(R) and KM(R, X = round(R)); RANK1 using a regular and a logarithmic scale;
400
Resolution-enhanced KMD plots from the deisotoped mass spectrum of P(EO-b-PO-b-EO) using EO/X
401
(X = 29,…,88); Regular and resolution-enhanced DP plots for P(EO-b-PO-b-EO) using EO and PO or
402
EO/43 and PO/57. Resolution-enhanced KMD plots from the raw mass spectrum of P(EO-b-PO-b-EO)
403
using EO/X (X = 29,…,88); Variation of KMDs and RANKs functions for R1 = EO, R2 = C18H34O, R3
404
= C2H4. Resolution-enhanced KMD plots from the deisotoped mass spectrum of PEO monostearate
405
using EO/X (X = 29,…,88); Resolution-enhanced KMD plots from the mass spectrum of a terpolyester
406
using C10H16O4/X (X = 133,…,257); Variation of KMDs and RANKs functions for R1= C10H16O4, R2
407
= C11H18O4, R3 = C12H20O4.; Resolution-enhanced KMD plots of the selected ion series using
408
C10H16O4/X (X = 136, 150, 164, 178, 193, 207, 221, 236 and 250, local maxima of RANK2);
409
Screenshots from our Excel sheet with RANK functions.
410 411
AUTHOR INFORMATION
412
Corresponding Authors
413
* Email:
[email protected] 414
ORCID: Thierry Fouquet: 0000-0002-9473-9425
415
* Email:
[email protected] 416
ORCID: Robert Cody: 0000-0002-6624-8530
417 418
REFERENCES
419
(1) Kendrick, E. A Mass Scale Based on CH2 = 14.0000 for High Resolution Mass Spectrometry of
420
Organic Compounds. Anal. Chem. 1963, 35, 2146-2154.
421
(2) Hsu, C. S.; Qian, K.; Chen, Y. C. An innovative approach to data analysis in hydrocarbon
422
characterization by on-line liquid chromatography-mass spectrometry. Anal. Chim. Acta 1992, 264, 79-
423
89.
424
(3) Krajewski, L. C.; Rodgers, R. P.; Marshall, A. G. 126 264 Assigned Chemical Formulas from an
425
Atmospheric Pressure Photoionization 9.4 T Fourier Transform Positive Ion Cyclotron Resonance
426
Mass Spectrum. Anal. Chem. 2017, 89, 11318–11324.
427
(4) Hughey, C. A.; Hendrickson, C. L.; Rodgers, R. P.; Marshall, A. G. Kendrick Mass Defect
428
Spectrum: A Compact Visual Analysis for Ultrahigh-Resolution Broadband Mass Spectra. Anal. Chem.
429
2001, 73, 4676–4681.
ACS Paragon Plus Environment
Page 16 of 18
Page 17 of 18 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Analytical Chemistry
430
(5) Sato, S.; Nakamura, S.; Teramoto, K.; Sato, T. Structural Characterization of Polymers by MALDI
431
Spiral-TOF Mass Spectrometry Combined with Kendrick Mass Defect Analysis. J. Am. Soc. Mass
432
Spectrom. 2014, 25, 1346-1355.
433
(6) Fouquet, T.; Aizawa, H.; Sato, H. Taking MALDI SpiralTOF high-resolution mass spectrometry
434
and mass defect analysis to the next level with ethylene vinyl acetate vinyl alcohol terpolymers. Rapid
435
Commun. Mass Spectrom. 2016, 30, 1818-1822.
436
(7) Fouquet, T.; Mertz, G.; Delmee, M.; Becker, C.; Bardon, J.; Sato, H. The Definitive Evidence of a
437
Plasma Copolymerization of Alkyl and Perfluorinated Acrylates Using High Resolution Mass
438
Spectrometry and Mass Defect Analysis. Plasma Process. Polym. 2016, 13, 862-868.
439
(8) Cody, R. B; Fouquet, T. Paper spray and Kendrick mass defect analysis of block and random
440
ethylene oxide/propylene oxide copolymers. Anal. Chim. Acta 2017, 989, 38-44.
441
(9) Qi, Y.; Hempelmann, R.; Volmer, D. A. Two-dimensional mass defect matrix plots for mapping
442
genealogical links in mixtures of lignin depolymerisation products. Anal. Bioanal. Chem. 2016, 408,
443
4835-4843.
444
(10) Dier, T. K. F.; Egele, K.; Fossog, V.; Hempelmann, R.; Volmer, D. A. Enhanced Mass Defect
445
Filtering To Simplify and Classify Complex Mixtures of Lignin Degradation Products. Anal. Chem.
446
2016, 88, 1328-1335.
447
(11) Morgan, T. E.; Ellacott, S. H.; Wootton, C. A.; Barrow, M. P.; Bristow, A W. T.; Perrier, S.;
448
O'Connor P. B. Coupling Electron Capture Dissociation and the Modified Kendrick Mass Defect for
449
Sequencing of a Poly(2-ethyl-2-oxazoline) Polymer. Anal. Chem. 2018, 90, 11710-11715.
450
(12) Fouquet, T.; Sato, H. Extension of the Kendrick Mass Defect Analysis of Homopolymers to Low
451
Resolution and High Mass Range Mass Spectra Using Fractional Base Units. Anal. Chem. 2017, 89,
452
2682-2686.
453
(13) Fouquet, T.; Shimada, H.; Maeno, K.; Ito, K.; Ozeki, Y.; Kitagawa, S.; Ohtani, H.; Sato, H. High-
454
resolution Kendrick Mass Defect Analysis of Poly(ethylene oxide)-based Non-ionic Surfactants and
455
Their Degradation Products. J. Oleo Sci. 2017, 9, 1061-1072.
456
(14) Fouquet, T.; Sato, H. Improving the Resolution of Kendrick Mass Defect Analysis for Polymer
457
Ions with Fractional Base Units. Mass Spectrom. (Tokyo) 2017, 6, A0055.
458
(15) Fouquet, T.; Cody, R. B.; Ozeki, Y.; Kitagawa, S.; Ohtani, H.; Sato, H. On the Kendrick Mass
459
Defect Plots of Multiply Charged Polymer Ions: Splits, Misalignments, and How to Correct Them. J.
460
Am. Soc. Mass Spectrom. 2018, 29, 1611-1626.
461
(16) Zheng, G.; Morimoto, M.; Sato, H.; Fouquet, T. Resolution-enhanced Kendrick mass defect plots
462
for the data processing of mass spectra from wood and coal hydrothermal extracts. Fuel, 2019, 235,
463
944-953.
464
(17) Fouquet, T.; Sato, H. How to choose the best fractional base unit for a high‐resolution Kendrick
465
mass defect analysis of polymer ions. Rapid. Commun. Mass Spectrom. 2017, 31, 1067-1072.
ACS Paragon Plus Environment
Analytical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
466
(18) Sato, T.; Teramoto, K.; Satoh, T.; Ueda, Y. Ethylene Oxide-Propylene Oxide block Copolymer
467
Analysis by MALDI Spiral-TOF/TOFKobunshi Ronbunshu 2012, 69, 406-415.
468
(19) Satoh, T.; Tsuno, H.; Iwanaga, M.; Kammei, Y. The design and characteristic features of a new
469
time-of-flight mass spectrometer with a spiral ion trajectory. J. Am. Soc. Mass Spectrom. 2005, 16,
470
1969-1975.
471
(20) Trimpin, S.; Keune, S.; Räder, H.J.; Müllen, K. Solvent-free MALDI-MS: Developmental
472
improvements in the reliability and the potential of MALDI in the analysis of synthetic polymers and
473
giant organic molecules. J. Am. Soc. Mass Spectrom. 2006, 17, 661-671.
474
(21) Strohalm, M.; Kavan, D.; Novák, P.; Volný, M.; Havlícek, V. mMass 3: A Cross-Platform
475
Software Environment for Precise Analysis of Mass Spectrometric Data. Anal. Chem. 2010, 82, 4648-
476
4651.
477
(22) Fouquet, T.; Cody, R.B.; Sato, H. Capabilities of the remainders of nominal Kendrick masses and
478
the referenced Kendrick mass defects for copolymer ions. J. Mass Spectrom. 2017, 52, 618-624.
ACS Paragon Plus Environment
Page 18 of 18