Structural Analysis of Phospholipid using Hydrogen Abstraction

last 20 years. .... (QIT) time-of-flight (TOF) mass spectrometer, following our previous study.5. Two ... The final 15-mm-long section of the quartz g...
2 downloads 0 Views 2MB Size
Subscriber access provided by Kaohsiung Medical University

Article

Structural Analysis of Phospholipid using Hydrogen Abstraction Dissociation and Oxygen Attachment Dissociation in Tandem Mass Spectrometry Hidenori Takahashi, Yuji Shimabukuro, Daiki Asakawa, Shosei Yamauchi, Sadanori Sekiya, Shinichi Iwamoto, Motoi Wada, and Koichi Tanaka Anal. Chem., Just Accepted Manuscript • DOI: 10.1021/acs.analchem.8b00322 • Publication Date (Web): 24 May 2018 Downloaded from http://pubs.acs.org on May 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 28 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

Structural Analysis of Phospholipid using Hydrogen Abstraction Dissociation and Oxygen Attachment Dissociation in Tandem Mass Spectrometry

Hidenori Takahashi,1* Yuji Shimabukuro,2 Daiki Asakawa,3 Shosei Yamauchi,1 Sadanori Sekiya,1 Shinichi Iwamoto,1 Motoi Wada,2 and Koichi Tanaka1

1. Koichi

Tanaka

Mass

Spectrometry

Research

Laboratory,

Shimadzu

Corporation,

1

Nishinokyo-Kuwabaracho Nakagyo-ku, Kyoto 604-8511, Japan; 2. Graduate School of Science and Engineering, Doshisha University, 1-3 Kyotanabe, Kyoto 610-0321, Japan; 3. National Institute of Advanced Industrial Science and Technology (AIST), National Metrology Institute of Japan (NMIJ), Tsukuba Central 2, 1-1-1 Umezono, Tsukuba, Ibaraki, 305-8568, Japan

Correspondence to: Hidenori Takahashi Koichi

Tanaka

Mass

Spectrometry

Research

Laboratory,

Shimadzu

Corporation,

1

Nishinokyo-Kuwabaracho Nakagyo-ku, Kyoto 604-8511, Japan TEL: +81-75-823-1482, E-mail: [email protected]

ORCID: Hidenori Takahashi (0000-0001-6887-1724) Yuji Shimabukuro (0000-0003-4090-3995) Daiki Asakawa (0000-0002-9357-8420)

Keywords; Phospholipid, Radical-induced dissociation, Hydrogen atom, Oxygen atom, Hydroxyl radical, Tandem mass spectrometry

ACS Paragon Plus Environment 1

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 2 of 28

ABSTRACT Gas-phase hydrogen radicals were introduced into a quadrupole ion trap containing singly charged phospholipids to obtain structural fragmentation patterns in tandem mass spectrometry (MS/MS). Saturated and unsaturated phosphatidylcholines were used as a model phospholipid, whose chain-length ranges between 16 and 24.

The MS/MS spectrum yielded a continuous series of

fragment ions with a mass difference of 14 Da, representing the saturated fatty acyl chains.

The

fragment ions corresponding to the double-bond position within a single fatty acyl chain showed a characteristic mass difference of 12 Da.

The detection of these diagnostic product ions enabled the

structural analysis of double-bond isomers of phospholipids.

To further investigate the potential of

radical-induced dissociation for the isomeric analysis of phospholipids, gas-phase hydroxyl radicals and triplet oxygen atoms were employed in tandem mass spectrometry.

The methylene bridges

adjacent to the double-bond positions were selectively dissociated, accompanied by oxidation of the double bonds.

Tandem mass spectrometry incorporating multiple radical species facilitates the

structural analysis of isomeric phospholipids.

Graphic Abstract Sequential Fragments within Acyl-Chains

H• or

Ion Trap MS/MS

OH• O Double-bond Specific Fragments

ACS Paragon Plus Environment 2

Page 3 of 28 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

INTRODUCTION Tandem mass spectrometry (MS/MS) is an effective analytical technique for the characterization of structurally diverse biopolymers, including proteins/peptides, nucleic acids, carbohydrates, and lipids.

However, distinguishing isomeric compounds using the conventional

low-energy collision-induced dissociation (LE-CID) MS/MS technique often remains ambiguities. This is because LE-CID predominantly cleaves only labile groups of target ions and consequently yields diagnostic fragment ions unsuitable for isomeric structure analysis.

Several radical-induced

dissociation (RID) MS/MS techniques complementary to LE-CID MS/MS have been devised in the last 20 years.

The commercially available RID MS/MS of electron capture dissociation (ECD)1 and

that of electron transfer dissociation (ETD)2 have been successfully applied to identify and characterize proteins and their post-translational modifications (PTMs).3,4

Regarding the

fragmentation process, ECD/ETD are initiated by fragile radical intermediate formation through the electron capture/transfer of a multiply charged peptide.

The radical ions then induce dissociation of

the peptide backbone (N-Cα) bonds, preserving the labile modifications, which enables the isomer analysis of PTM peptides.

However, ECD/ETD cannot be applied to singly charged molecules

because electron attachment to singly-charged molecules produces only neutral species, which cannot be detected by mass spectrometry.

Although peptides and proteins can be ionized into

multiply charged ions by electrospray ionization, small molecules are generally ionized to only singly charged ions.

Therefore, it is quite difficult to characterize small molecules using

ECD/ETD. To address this limitation, the RID technique of hydrogen attachment/abstraction dissociation (HAD) was developed by our research group in 2016.5-7

Since HAD employs

interactions between the target precursor ions and a neutral hydrogen radical (H•) inside an ion trap, even singly charged molecules are dissociated in both positive- and negative-ion modes. HAD is potentially applicable to all types of ionization methods.

Therefore,

Previously, we demonstrated the

usefulness of this technique for peptide sequence analysis combined with matrix-assisted laser

ACS Paragon Plus Environment 3

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 28

desorption/ionization (MALDI),8-9 which preferentially generates singly charged molecules. Consequently, HAD was shown to be effective in characterizing not only multiply charged peptides and proteins, but also singly charged small molecules. In this study, HAD was applied to the structural characterization of small molecules of phospholipids for the assignment of double-bond positions in the fatty acyl chains and sn-1/sn-2 positions.

The chain length, degree of saturation, and position of the double bond of the fatty acids

determine their function of biomolecules.

The conventional LE-CID method provides limited

structural information for lipid analysis, because only the labile polar head-group is fragmented by this approach. phospholipids.10

Therefore, LE-CID cannot effectively discriminate the double-bond isomers of In contrast to LE-CID, RID provides detailed structural information for

biomolecules while preserving the labile groups, as reported in our previous study of HAD for PTM analysis5-7 and several studies of RID-based techniques.11-15

To expand the original concept of

HAD-MS/MS, gas-phase hydroxyl radicals (OH•) and triplet oxygen atoms (3O) were also employed for phospholipid dissociation.

Double-bond specific fragmentation within the fatty acyl chains can

be triggered by either OH• or 3O.

Similar double-bond specific fragmentation was reported under

ozone-induced dissociation (OzID) by Thomas and co-workers.16

OzID can identify double-bond

positions and sn- positions, even in complex biological mixtures.17-18

OzID also provides an easily

interpreted fragment spectrum and consequently a higher signal-to-noise ratio.

However, ozone

induces unwanted degradation of the mass spectrometer by oxidizing the electrodes and insulators in the vacuum chamber. Because the life time of ozone is on the order of several minutes,19 ozone interacts with both the target ions and apparatus materials inside the vacuum chamber.

Compared

to OzID, OH• and 3O are readily recombined to molecular form on the chamber wall,20,21 while the accumulated molecules are easily removed through baking the system. Therefore, the unexpected degradation of the mass spectrometer by reactive species is avoided.

Herein, we report how

MS/MS analysis incorporating H•, OH•, and 3O facilitates the structural analysis of isomeric compounds.

ACS Paragon Plus Environment 4

Page 5 of 28 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

EXPERIMENTS Materials Lipid

standards

(Lyso-PC

18:1(9Z),

PC

18:0/18:1(9Z),

PC

18:1/18:0(9Z),

PC

16:0/24:4(5Z,8Z,11Z,14Z)) and the MALDI matrix were purchased from Avanti Polar Lipids (Alabaster, AL) and Sigma-Aldrich (St. Louis, MO), respectively. dissolved in 100% methanol.

The sample solution was

2,5-Dihydroxybenzoic acid (DHB) matrix was dissolved at a

concentration of 10 mg/mL in 50% acetonitrile/0.1% aqueous tetrafluoroacetic acid (v/v).

The

sample solution (0.5 µL) and matrix solution (0.5 µL) were mixed well on a MALDI sample target. C60 fullerene powder was purchased from Sigma-Aldrich (St. Louis, MO).

C60 powder was

dissolved in 100% toluene, and the solution (1 µL) was spotted thoroughly on the MALDI sample target.

Tandem Mass Spectrometer All MS/MS experiments were performed using a prototype MALDI quadrupole ion trap (QIT) time-of-flight (TOF) mass spectrometer, following our previous study.5

Two

1.5-mm-diameter holes were arranged in the ring electrode of the ion trap for the entrance and exit of radical injection. electrode.

The radical sources were aligned with the axis of the two holes in the ring

A nitrogen laser (wavelength 337 nm) was employed for MALDI ion generation.

The

generated ions were subsequently transported to the ion trap through an aperture on an end-cap electrode.

Helium gas was used for cooling the trapped ions, and argon was used as the collision

gas for conventional CID. isolation (DAWI).22

Precursor ion isolation was achieved by digital asymmetric wave

The reaction time of the radical–ion interaction was set to 0.5–1 s.

The

pressure within the ion trap chamber was maintained below 5 × 10−4 Pa before the radical injection. Trapped ions were extracted with an accelerating voltage of 10 kV in reflectron mode.

Mass

spectra were acquired by averaging 50–100 single shots. H• was generated by passing hydrogen gas through a hot tungsten capillary (length of 50

ACS Paragon Plus Environment 5

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 28

mm, inner diameter of 1 mm) heated by thermal radiation from a surrounding hot tungsten filament. The temperature of the capillary was estimated to be approximately 2300 K.

The flow rate of the

input hydrogen gas through the capillary was controlled to below 5 sccm (1 sccm ≈ 4.5 × 1017 atoms/s) by a mass-flow controller (Model 3660, Kofloc, Kyoto, Japan). OH• and

3

O were generated by an in-house built radical source based on a

microwave-driven capacitively coupled plasma generator, as shown in Figure S1.23

A 2.45-GHz

microwave was applied to the quarter-inch (6.35 mm) outer diameter copper tube via an N-type radio frequency (RF) connector and a coaxial cable from a solid-state type microwave power generator (Tokyo Keiki Inc., TMEB101B00B).

The total length of the copper tube was 30.6 cm,

corresponding to 2.5 wavelengths of the 2.45-GHz microwave.

A quartz glass tube (inner diameter:

14 mm, outer diameter: 18 mm) enclosed a tapered copper electrode having a 1-mm-diameter gas outlet fixed to the tip of the copper tube.

The final 15-mm-long section of the quartz glass tube

tapered to the exit hole diameter of 4 mm.

An aluminum cylinder surrounded the tapered quartz

glass tube to maintain a plasma plume at the sharp tapered electrode tip by intensifying the local microwave electric field.

The flow rate of the discharge gas was controlled by a needle valve

(US-916P-P6.35, Fujikin, Japan).

We maintained a background gas pressure of the radical source

and the ion trap at 3 Pa and 0.05 Pa, respectively, which corresponds to a gas flow rate of ≈100 sccm. The higher background pressure (i.e., higher gas flow rate) consequently decreased the product ion intensity due to the decrease of the transport efficiency of the radicals into the ion trap.

The

presence of radical species inside the radical source was confirmed by optical emission spectroscopy acquired with a compact optical spectrometer (Ocean optics USB 2000+ and Ocean optics Flame-S). A Langmuir probe was used to measure the plasma parameters of the electron density and electron temperature.

The optical emission spectroscopy and plasma parameter measurements were

performed in a separate experimental setup which was disconnected from the mass spectrometer as shown in Figure S2.

ACS Paragon Plus Environment 6

Page 7 of 28 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

Calculations All electron structure calculations were performed with the Gaussian 16 program.24

The

geometries of radicals and molecules were obtained by gradient optimizations with density functional theory calculations using the M06-2X hybrid functional25 and the 6-31G(d) basis set. The structures were characterized by frequency calculations as local energy minima.

To establish

the energetics for fragmentation, the transition state geometries were also optimized at the M06-2X/6-31G(d) level and confirmed by vibrational frequency analysis, showing one imaginary frequency.

The connection of transition states to the reactants and intermediates was checked by

intrinsic reaction coordinate analysis26 of five forward and five reverse steps, starting from the transition state geometry.

The harmonic frequencies obtained by frequency analysis were used to

obtain zero-point energy corrections.

ACS Paragon Plus Environment 7

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 28

RESULTS AND DISCUSSION HAD-MS/MS of phospholipid In order to identify the HAD mechanism of phospholipids, we used phosphatidylcholine (PC 18:1(9Z)), which contains a single fatty acid chain, as the model compound.

Figure 1A shows

the HAD spectrum of singly protonated phosphatidylcholine (Lyso-PC 18:1(9Z), 522.36 Da, 10 pmol) formed through the interaction of H• with the isolated precursor ions inside the ion trap.

To

enhance the product ion intensity, supplemental collisional activation27 was applied to the precursor ions after H• injection (HAcaD).

The non-dissociative precursor ions were isolated prior to the

supplemental activation to eject the directly fragmented product ions, which undergo multiple reactions to form unfavorable products for the interpretation of the MS/MS spectrum. activation energy was optimized to obtain sufficient product ions for each analyte. was generated by a thermal cracking cell, following our previous study.5

The

The injected H•

The HAD spectrum of the

model phospholipid shows a continuous series of fragment ions with the mass difference of 14 Da, which represents a CH2 group, between m/z of 312 and 492.

In contrast, the conventional LE-CID

predominantly provides a head-group specific fragmentation, resulting in the undetectable sequential C–C bond cleavage as shown in Figure S3.

Since the HAD spectrum shows the fragmentation

produced by the C–C bond cleavage, HAD-MS/MS is a potentially useful method for the characterization of the double-bond position in fatty acids.

We first focused on the mechanism of

the C–C bond cleavage and then on the analysis of the double-bond position in phospholipids, as highlighted shortly. In order to determine the most probable fragmentation mechanism of HAD, the isotopic distribution of the precursor ion was compared before and after H• injection (shown in the inset of Figure 1A).

To avoid the complexities induced by the overlapping isotopic peaks of the product

ions, only the lightest isotopic peak was isolated prior to the H• injection.

In our previous study of

HAD for peptide analysis, both H• attachments and H• abstractions to and from the precursor ions showed significant results, indicating that H• was preferentially attached to the carbonyl oxygen

ACS Paragon Plus Environment 8

Page 9 of 28 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

atoms and abstracted from the amide hydrogen or β hydrogen atoms of the peptide backbone. contrast, only H• abstractions are observed after HAD of the phospholipid.

In

This result suggests that

the abundant reactive hydrogen atoms within the fatty acyl chain are dominantly abstracted by injected H•.

Additionally, the fragment ions observed between m/z of 312 and 492 are generated by

the C–C bond cleavage which accompanies double-bond formation.

Therefore, the fragmentation

of phospholipids is probably caused by induced hydrogen abstraction by H•. generated hydrogen-deficient radical underwent a C–C bond cleavage.

Subsequently, the

In this study, we focused on

the formation of the fragment ion at m/z 312, with the mechanism summarized in Scheme 1.

H•

abstracts a hydrogen atom in the acyl chain in Lyso-PC 18:1(9Z) to yield hydrogen-deficient Lyso-PC 18:1(9Z), R1, and H2.

The corresponding transition state barrier for the hydrogen

abstraction (TS1-1) has the energy of 34 kJ/mol; the products of R1 and H2 are more stable than the reactants by 42 kJ/mol.

Because the temperature of H• in the HAD experiment was approximately

2300 K, approximately 50% of H• possesses a kinetic energy exceeding 34 kJ/mol and can react with the phospholipids.

Next, the hydrogen-deficient Lyso-PC 18:1(9Z), R1, undergoes radical-induced

homolytic C–C bond cleavage, leading to double-bond formation. homolytic C–C bond cleavage (TS1-2) is 139 kJ/mol.

The transition state barrier for

The corresponding bond cleavage produced

the intermediate IM1 and the formation of the fragment ion at m/z 312 requires 106 kJ/mol. Because the fragmentation of C–C bond cleavage by HAD proceeds as an endoergic reaction, supplemental activation is necessary to obtain sufficient HAD products in this experiment. Next, we focused on the characterization of the phospholipid by HAD.

Figure 1B shows

the HAD spectrum of singly protonated phosphatidylcholine (PC 18:0/18:1(9Z), 788.6 Da, 10 pmol), which contains two fatty acid chains.

The HAD spectrum of the model phospholipid yields a

continuous series of fragment ions with a mass difference of 14 Da, which represents a CH2 group between m/z 562 and 758.

In addition, fragment ions corresponding to the double-bond position

show the characteristic mass difference of 26 Da between m/z 648 and 674, indicating the HC=CH group.

The fragment ion from the cleavage of the C=C double bond is observed at m/z 660,

ACS Paragon Plus Environment 9

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 28

demonstrating that the double-bond position can be identified by HAD. These diagnostic fragment ions enable the discrimination of double-bond isomers of phospholipids.

The inset of Figure 1B

shows the expanded view of the region around the fragment ions at 522 and 524 Da of PC 18:0/18:1(9Z) and that of the regioisomer PC 18:1(9Z)/18:0.

The fragment ions at 524 and 522 Da

represent the product ion from the acyl chain loss at sn-2 for PC 18:0/18:1 and the regioisomer PC 18:1/18:0, respectively.

The fragment ion peaks from the sn-2 chain loss are more intense than

those of the sn-1 chain loss in both cases.

This result suggests that the acyl chain loss from the sn-2

position is preferentially induced by H• abstraction from reactive α-hydrogens of the acyl chain at sn-2.

This feature is effective in the diagnostic process to determine the sn- position of the fatty

acyl chain.

ACS Paragon Plus Environment 10

Page 11 of 28 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

Figure 1. (A) HAcaD spectrum of Lyso-PC 18:1(9Z) ([M+H]+=522.36 Da). isotopic distributions of precursor ion before and after H• injection. PC 18:0/18:1(9Z) ([M+H]+=788.6 Da).

Inset shows

(B) HAcaD spectrum of

Inset shows detailed product ion spectra between 520

and 530 Da for PC 18:0/18:1(9Z) and for the regioisomer PC 18:1(9Z)/18:0.

ACS Paragon Plus Environment 11

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 12 of 28

Scheme 1. Fragmentation pathway of PC 18:1(9Z) by H• abstraction. The relative energy (kJ/mol) was obtained by M06-2X/6-31G(d) level calculations, including zero-point vibrational energies.

Generation of gas-phase H•, OH•, and 3O As described above, HAD is initiated by H• abstraction from a phospholipid.

HAD can

cleave the C–C bond in a fatty acyl chain with similar probabilities, because the probabilities of H• abstraction in fatty acyl chain would not depend on the C-C bond position.

As a result, the HAD

MS/MS spectra provided detailed information on the molecular structure of the phospholipids. However, when the target lipids contain complex fatty acid mixtures, the HAD spectrum is

ACS Paragon Plus Environment 12

Page 13 of 28 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

complicated by the production of many different fragment ions, which impedes the interpretation of the HAD spectrum and deteriorates the signal-to-noise ratio. As described in the introduction, the properties of phospholipids are strongly affected by the double bond position in the fatty acyl chain. reaction with an electrophilic reagent.

The C=C double bond readily undergoes an addition

The reactive species containing oxygen atoms attaches to the

double bond and induces the C=C bond-specific fragmentation.

Herein, we used O2 and H2O as the

radical source, in order to determine the double bond position in the fatty acyl chain.

Firstly, we

analyzed the composition of neutral reactive species generated by the microwave discharge of H2O and that of O2.

The microwave discharge can dissociate a wide variety of reactive gases, yielding

reactive neutral species with lower operating gas pressure (< 10 Pa) and lower input power (< 100 W),28,29 which is suitable for incorporating into a conventional mass spectrometer. Figure 2 shows the typical optical emission spectra for the microwave discharge of H2O with an input power of 100 W and a gas flow rate of 100 sccm.

Under typical experimental

conditions, the estimated electron density and the electron temperature on the capillary axis at 3 cm downstream from the end of the capillary were 1×1010 cm-3 and 3 eV, respectively.

The brightest

peak is observed at 309 nm, corresponding to the ultraviolet-range OH emission band.

Other

intense peaks are found at 656 nm, 486 nm, 434 nm, and 410 nm, corresponding to H-atom emission lines (Balmer-series spectra) and 777 nm and 844 nm for O-atom emission lines.

This result shows

that OH•, H•, and 3O are successfully generated by the microwave discharge of H2O.

Figure S4

shows the optical emission spectrum for the microwave discharge of pure O2 gas, which indicates that 3O is predominantly produced in the discharge process.

In the ultraviolet region of the O2

discharge, the OH emission band does not show sufficient intensity compared to the background molecular spectra.

ACS Paragon Plus Environment 13

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 14 of 28

Figure 2. Optical emission spectra of the microwave discharge of H2O. (A) Ultraviolet region, (B) visible-light region.

To confirm whether or not the generated radicals are efficiently transported into the ion trap from the radical source, we observed the attachment of H• and O to the fullerene ion C60+•, which is highly reactive with radical species.30

Figure 3 shows the isotopic distribution of C60+• after the

irradiation of the microwave discharge products from (A) H2O and (B) pure O2 at the reaction time of 1 s.

The attachment of more than 4H• and 2O to C60+• are clearly observed, indicating that the

generated radicals are successfully introduced to the ion trap. O attachments are observed.

In the case of pure O2 discharge, only

Assuming that the H• temperature does not alter the H• attachment

coefficient, the number density of H• inside the ion trap is estimated to be 4×1010 cm-3 from the observed H• attachment reaction rate (4 atom attachments per 1 s) and the reported absolute value of the H• attachment coefficient to C60+ (≈1×10-10 cm-3 s).31 Because H• attachment to C60+• is the radical recombination reaction, the threshold energy of the H• attachment to C60+• should be zero or quite low, resulting in the temperature independence of the H• attachment coefficient.

Whereas, the

number density of O and that of OH• cannot be estimated from Figure 3, because we cannot find the absolute value of their attachment coefficients in any literature.

Further study is required to

estimate the number density of radicals inside the ion trap using laser-induced fluorescence (LIF) measurements.32

ACS Paragon Plus Environment 14

Page 15 of 28 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

Figure 3. Mass spectra of H• and 3O attachments to fullerene ion C60+ induced by (A) H2O discharge and (B) pure O2 discharge at the reaction time of 1 s.

Oxygen Attachment Dissociation of Phospholipid Next, we applied the microwave-driven radical source to the fragmentation of phospholipids. Figure 4A shows the MS/MS spectra of the protonated Lyso-PC 18:1(9Z) obtained by the microwave discharge of H2O vapor.

The irradiation of the microwave discharge products from H2O produced

the fragment ions at m/z 424.21 and 426.22, which were identified as C18H34NO8P and C18H36NO8P, respectively. The proposed structure of the fragment ions is shown in Scheme 2.

The fragment ions

would be formed by oxidation of the C=C double bond accompanied by C–C bond cleavage adjacent to the oxidation site. reference.

To confirm the fragmentation pathway, we used saturated phospholipid as the

As shown in Figure S5, product ions corresponding to the cleavage of fatty acyl chains

were absent from the saturated phospholipid, which suggests that the C–C bond dissociation is triggered by the specific oxidation of the double bond.

The fragmentation of the fatty acyl chain is

induced by the interaction between the C=C double bond and either OH• or 3O, because the

ACS Paragon Plus Environment 15

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 16 of 28

microwave discharge of the H2O vapor produced H•, OH•, and 3O, as described in the previous section.

Notably, microwave discharge of H2O produced H•, which can induce C–C bond cleavage,

as shown in Figure 1. However, the fragment ions due to the C–C bond cleavage were absent in Figure 4A.

The efficiency of hydrogen abstraction by H• appeared to be lower than that of

oxidation of the C=C double bond by OH• or 3O.

In contrast, the product ions observed below m/z

350 are attributable to H• interaction, because these product ions are also observed in the HAD spectrum of Figure 1. To elucidate the fragmentation process of OH• and 3O individually, we performed an MS/MS experiment on Lyso-PC(18:1(9Z)) using the microwave discharge of pure oxygen gas.

As

shown in Figure 4B, several characteristic fragment ions are observed at m/z 398.19 (C16H32NO8P), 412.22 (C17H34NO8P), 415.21 (C17H37NO8P), 424.21 (C18H34NO8P), and 440.21 (C18H34NO9P). The proposed structure of the fragment ions is shown in Scheme 2. A variety of product ions should be generated by the subsequent reaction following the double-bond oxidation because these product ions are not observed for the saturated phospholipid with no double bonds (see Figure S6). In addition, background oxygen molecules (0.05 Pa ≈ 1.3×1013 cm-3) are more abundant than 3O inside the ion trap because the degree of dissociation of the oxygen molecules is less than several tens of percent in the microwave discharge.

The oxygen molecules are reported to react with alkyl

radicals to yield acetyl peroxide radicals.33

This highly reactive species generates a variety of

fragments.

Therefore, these characteristic fragment ions in the mass spectrum likely arose from

secondary interactions between non-dissociative oxidized precursor ions and the background oxygen gas.

To validate this proposed secondary reaction with background oxygen, oxygen gas was

introduced directly into the ion trap after the injection of the gas-phase H•, OH•, and 3O mixture from the microwave discharge of H2O vapor (Figure 5C). similar to those observed in Figure 5B.

The observed product ions are quite

The secondary interactions between non-dissociative

oxidized ions and background H2O gas are avoided, yielding straightforward double-bond specific fragmentation using the microwave discharge of H2O vapor because the reactivity of H2O vapor is

ACS Paragon Plus Environment 16

Page 17 of 28 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

quite low compared to that of oxygen.

DFT Calculations for C=C Bond Dissociation As previously described, the double-bond specific fragmentation is initiated by the attachment of 3O and/or OH• to the double bond.

In order to study the most probable fragmentation

pathway, we next investigated the reaction of Lyso-PC(18:1(9Z)) with 3O and OH• by performing DFT calculation. Firstly, we discussed the reaction between 3O and Lyso-PC 18:1(9Z) and the calculated reaction pathway is shown in Scheme 3.

The 3O interacts with the C=C double bond in

Lyso-PC 18:1(9Z) to produce the intermediate complex IM3-1, and the corresponding energy was 22 kJ/mol.

Subsequently, the 3O bonds to the C=C double bond to yield R3 through the transition state

TS3-1.

The transition state barrier for 3O addition is 17 kJ/mol.

R3 is a triplet state diradical,

containing oxyl and alkyl radicals, as shown in Scheme 3. According to the investigation on oxidation of hydrocarbon, the alkyl radical immediately reacts with O2 and we assumed that the alkyl radical in R3 was converted into the peroxide radical.

R3 and O2 adduct on R3 undergo a variety

of radical-induced dissociation, leading to fragment ions at m/z 398.19, 412.22, 415.21, 424.21, and 440.21.

As previously described, the secondary reaction between the analyte radical and O2 was

observed in Figure 5C.

Therefore, fragment ions at m/z 398.19, 412.22, 415.21, 424.21, 426.22 and

440.21 in Figure 5B and Figure 5C are attributable to the interactions between R3 and the abundant background O2 gas in Scheme 3.

The generated radical immediately reacts with O2 and a variety of

fragment ions were generated by subsequent radical-induced dissociation, when O2 was present in the ion trap.

The reaction between the analyte radical and O2 hampers the identification of the

double bond position in the fatty acyl chain. In contrast to O2, H2O does not react with the analyte radical.

As a result, the use of

discharge products from H2O for MS/MS provides the information for identifying the double bond position in the fatty acyl chain.

A comparison of Figure 4A and 4B reveals that the fragment ion at

m/z 426.22 is only generated by H2O discharge containing H• and OH•.

This result suggests that

ACS Paragon Plus Environment 17

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 18 of 28

the fragment ion at 426.22 is attributable to either OH• attachment to [PC(18:1(9Z))+H]+ or H• attachment to R3.

As previously estimated, the number density of H• inside the ion trap is about 3

orders of magnitude lower than that of background O2 or H2O.

In addition, the peak intensity of

[M+H+O]+ is only 3% of the precursor ion intensity in Figure 4A, and the number density of the unstable radical R3 should be much lower than that of the stable oxidized product [M+H+O]+ with an epoxide group.

Thus, it is supposed that the sequential radical reaction of H• attachment to R3 is

a negligible process.

Herein, we focused on the reaction between OH• and the protonated Lyso-PC

18:1(9Z), in order to understand the formation mechanism of fragment ion at m/z 426.22.

Scheme

4 shows the reaction pathway of OH• and the protonated Lyso-PC 18:1(9Z) obtained by DFT calculation.

The OH• interacts with the C=C double bond in Lyso-PC 18:1(9Z) to produce the

intermediate complex IM4-1. was 24 kJ/mol.

The interaction energy of the protonated Lyso-PC 18:1(9Z) and OH•

Subsequently, the OH• bonds to the C=C double bond to yield the radical R4

through the transition state TS4-1.

The transition state barrier for OH• addition is 19 kJ/mol.

As

described in Scheme 1, the transition state barrier for hydrogen abstraction by H• is 34 kJ/mol, which is much higher than that for OH• addition. Consequently, the Lyso-PC 18:1(9Z) preferentially reacts with OH• and the formation of the hydrogen-deficient Lyso-PC 18:1(9Z) radical is suppressed in the presence of OH•.

The radical R4 then undergoes radical-induced C–C bond cleavage,

yielding to enol and alkyl radicals. by 128 kJ/mol.

The corresponding transition state TS4-2 is less stable than R4

Finally, the fragment ion at m/z 426 is formed through the intermediate complex

IM4-2 in exoergic fragmentation at 38 kJ/mol.

The formation process of the fragment ion at m/z

426.22 is reasonably explained by DFT calculation. Regarding the HAD-MS/MS of the phospholipid shown in Figure 1, supplemental activation is necessary to obtain sufficient fragment ion signals.

In contrast, sufficient fragment ions are

obtained without supplemental activation of non-dissociative precursor ions when using the microwave discharge of H2O vapor and that of pure O2 gas.

The irradiation of the phospholipids by

OH• and 3O of the phospholipid yields stable oxidized products and subsequent fragmentations

ACS Paragon Plus Environment 18

Page 19 of 28 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

proceed simultaneously as exoergic reactions. As a result, supplemental activation is unnecessary. Fragmentation by HAD proceeds via endoergic reactions, which requires supplemental activation to generate fragment ions.

The attachment of 3O and OH• to the double bond induced fragmentation.

In contrast, the triplet state diradical formed by 3O attachment provide more stable singlet state species by intersystem crossing.

Therefore, R3 would result in a stable form of [M+H+O]+ with an

epoxide group, which is more stable than the reactants by 388 kJ/mol (free form of 3O and Lyso-PC 18:1(9Z)), as shown in Scheme 3. radical-driven dissociation.

The epoxide is singlet molecule and does not induce further

In order to investigate the properties of [M+H+O]+, we performed

CID-MS3 experiments for [M+H+O]+ (Figure 5).

The collisional activation for [M+H+O]+ yields

only a head-group fragment at m/z 184.07 with the absence of double-bond fragmentation.

This

result supports the idea that the observed O attached precursor ions are predominantly in the epoxide form which is more stable than the labile head-group.

The CID-MS3 spectrum of [M+H+O]+

generated by H2O discharge also shows an identical result, as indicated in Figure S7. In summary, fragmentation using the microwave discharge of H2O vapor is a potentially useful method for the characterization of the double-bond position in phospholipids.

To confirm

the usefulness of the method, we analyzed PC 20:4/14:0, which contains four double bonds, as a model phospholipid.

Figure 6 shows the fragment ions at m/z 608.39, 648.42, 688.47, and 728.48,

which originate from the cleavage of the C–C bond adjacent to the C=C double bond accompanied by the addition of an oxygen atom.

Therefore, the position of double bonds can be successfully

identified by MS/MS experiments using the microwave discharge of H2O vapor.

ACS Paragon Plus Environment 19

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 20 of 28

Figure 4. MS/MS spectra of phosphatidylcholine (PC 18:1(9Z)) with dissociation induced by (a) microwave discharge of H2O, (b) that of O2, and (c) that of H2O followed by O2 gas injection.

ACS Paragon Plus Environment 20

Page 21 of 28 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

Scheme 2. Proposed structure of fragment ions observed in Figure 4.

ACS Paragon Plus Environment 21

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 22 of 28

Scheme 3. Fragmentation pathway of PC 18:1(9Z) by 3O attachment. The relative energy (kJ/mol) was obtained by M06-2X/6-31G(d) level calculations, including zero-point vibrational energies.

ACS Paragon Plus Environment 22

Page 23 of 28 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

Scheme 4. Fragmentation pathway of PC 18:1(9Z) by OH• attachment. The relative energy (kJ/mol) was obtained by M06-2X/6-31G(d) level calculations, including zero-point vibrational energies.

ACS Paragon Plus Environment 23

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

Figure 5.

Page 24 of 28

CID-MS3 spectrum of isolated oxygen attached to phosphatidylcholine (Lyso-PC

18:1(9Z), [M+H+O]+ = 538.34 Da).

Oxygen attachment was induced by the injection of 3O from

the microwave discharge of pure O2 gas.

ACS Paragon Plus Environment 24

Page 25 of 28 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

Figure 6. MS/MS spectrum of phosphatidylcholine (PC 16:0/24:4) induced by H•, 3O, and OH• mixture (microwave discharge of H2O).

ACS Paragon Plus Environment 25

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 26 of 28

CONCLUSIONS The HAD-MS/MS of phospholipids was demonstrated using MALDI-QIT-TOF MS to identify the double-bond positions within the acyl chains and assign the sn- positions of the acyl chains.

The HAD-MS/MS spectra showed a continuous series of fragment ions with a mass

difference of 14 Da for the saturated acyl chains and characteristic ion peaks with a mass difference of 12 Da for unsaturated acyl chains.

These diagnostic fragment ions enabled the assignment of

double-bond positions within the acyl chains. To expand the original concept of HAD-MS/MS, OH• and 3O were employed for phospholipid dissociation to obtain the double-bond specific fragmentation analogous to that by OzID.

As a result, the methylene bridges adjacent to the double

bonds were selectively fragmented by either OH• or 3O, following the oxidation of the double bonds. These fragment ions are useful for diagnostics to realize detailed structural analyses of lipids in complex biological mixtures using MS/MS.

ACKNOWLEDGMENTS All electron structure calculations were performed at the Research Center for Computational Science, Okazaki, Japan.

ACS Paragon Plus Environment 26

Page 27 of 28 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

REFERENCES 1.

Zubarev, R. A.; Kelleher, N. L.; McLafferty, F. W. A nonergodic process. J. Am. Chem. Soc. 1998, 120(13), 3265-3266.

2.

Syka, J. E.; Coon, J. J.; Schroeder, M. J.; Shabanowitz, J.; Hunt, D. F. Proc. Natl. Acad. Sci.

U.S.A. 2004, 2004 101(26), 9528-9533. 3.

Stensballe, A.; Jensen, O. N.; Olsen, J. V.; Haselmann, K. F.; Zubarev, R. A. Rapid Commun.

Mass Spectrom. 2000, 2000 14(19), 1793-1800. 4.

Mikesh, L. M.; Ueberheide, B.; Chi, A.; Coon, J. J.; Syka, J. E.; Shabanowitz, J.; Hunt, D. F.

Biochim. Biophys. Acta. 2006, 2006 1764(12), 1811-1822. 5.

Takahashi, H.; Sekiya, S.; Nishikaze, T.; Kodera, K.; Iwamoto, S.; Wada, M.; Tanaka, K. Anal.

Chem. 2016, 88 (7), 3810–3816. 6.

Asakawa, D.; Takahashi, H.; Iwamoto, S.; Tanaka, K. Anal. Chem. 2018, 2018 90, 2701−2707.

7.

Asakawa, D.; Takahashi, H.; Iwamoto, S.; Tanaka, K. Phys. Chem. Chem. Phys. in press, DOI: 10.1039/C8CP00733K.

8.

Tanaka, K.; Waki, H.; Ido, Y.; Akita, S.; Yoshida, Y.; Yoshida, T.; Matsuo, T. Rapid Commun. Mass

Spectrom. 1988, 1988 2(8), 151-153. 9.

Karas, M.; Hillenkamp, F. Anal. Chem. 1988, 1988 60(20), 2299-2301.

10. Shimma, S.; Sugiura, Y.; Hayasaka, T.; Zaima, N.; Matsumoto, M.; Setou, M. Anal. Chem. 2008 2008, 08 80(3), 878-885. 11. Pham, H. T.; Ly, T.; Trevitt, A. J.; Mitchell, T. W.; Blanksby, S. J. Anal. Chem. 2012, 2012 84 (17), 7525-7532. 12. Li, P.; Hoffmann, W. D.; Jackson, G. P. Int. J. Mass Spectrom. 2016, 2016, 403, 1-7. 13. Li, P.; Jackson, G. P. J. Mass Spectrom. 2017, 2017 52, 271-282. 14. Pham, H. T.; Julian, R. R. Analyst, 2016, 2016 141(4), 1273-1278. 15. Rustam, Y. H.; Reid, G. E.; Anal. Chem. 2017, 017 90(1), 374-397. 16. Thomas, M. C.; Mitchell, T. W.; Harman, D. G.; Deeley, J. M.;Nealon, J. R.; Murphy, R. C.; Blanksby, S. J. Anal. Chem. 2007, 2007 79, 5013−5022. 17. Poad, B. L.; Pham, H. T.; Thomas, M. C.; Nealon, J. R.; Campbell, J. L.; Mitchell, T. W.; Blanksby, S. J. J. Am. Soc. Mass Spectrom. 2010, 2010 21(12), 1989-1999. 18. Poad, B. L.; Green, M. R.; Kirk, J. M.; Tomczyk, N.; Mitchell, T. W.; Blanksby, S. J. Anal. Chem. 2017, 2017 89(7), 4223-4229. 19. Zhukov, V., Popova, I.; Yates Jr, J. T. J. Vac. Sci. Technol. 2000, 2000 18(3), 992-994. 20. Smith, W. V. J. Chem. Phys., 1943, 1943 11(3), 110-125. 21. Magne, L.; Coitout, H.; Cernogora, G.; Gousset, G. J. de Phys. III. 1993, 1993 3(9), 1871-1889. 22. Brancia, F. L.; McCullough, B.; Entwistle, A.; Grossmann, J. G.; Ding, L. J. Am. Soc. Mass

Spectrom. 2010, 2010 21(9), 1530-1533. 23. Shimabukuro, Y.; Takahashi, H.; Wada, M. Jpn. J. Appl. Phys. 2017, 2017 57(1S), 01AA02. 24. Gaussian 16; Revision A.03; Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb,

ACS Paragon Plus Environment 27

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 28 of 28

M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Petersson, G. A.; Nakatsuji, H.; Li, X.; Caricato, M.; Marenich, A. V.; Bloino, J.; Janesko, B. G.; Gomperts, R.; Mennucci, B.; Hratchian, H. P.; Ortiz, J. V.; Izmaylov, A. F.; Sonnenberg, J. L.; Williams-Young, D.; Ding, F.; Lipparini, F.; Egidi, F.; Goings, J.; Peng, B.; Petrone, A.; Henderson, T.; Ranasinghe, D.; Zakrzewski, V. G.; Gao, J.; Rega, N.; Zheng, G.; Liang, W.; Hada, M.; Ehara, M.; Toyota, K.; Fukuda, R.; Hasegawa, J.; Ishida, M.; Nakajima, T.; Honda, Y.; Kitao, O.; Nakai, H.; Vreven, T.; Throssell, K.; Montgomery, J. J. A.; Peralta, J. E.; Ogliaro, F.; Bearpark, M. J.; Heyd, J. J.; Brothers, E. N.; Kudin, K. N.; Staroverov, V. N.; Keith, T. A.; Kobayashi, R.; Normand, J.; Raghavachari, K.; Rendell, A. P.; Burant, J. C.; Iyengar, S. S.; Tomasi, J.; Cossi, M.; Millam, J. M.; Klene, M.; Adamo, C.; Cammi, R.; Ochterski, J. W.; Martin, R. L.; Morokuma, K.; Farkas, O.; Foresman, J. B.; Fox, D. J., Gaussian, Inc., Wallingford CT 2016. 2016 25. Zhao, Y.; Truhlar, D. G., Theor. Chem. Acc. 2008, 2008 120, 215-241. 26. Fukui, K., Acc. Chem. Res. 1981, 1981 14, 363-368. 27. Swaney, D. L.; McAlister, G. C.; Wirtala, M.; Schwartz, J. C.; Syka, J. E.; Coon, J. J. Anal. Chem. 2007, 2007 79(2), 477-485. 28. Hopwood, J. Plasma Sources Sci. Technol. 1992, 1992 1, 109-116. 29. Ohachi, T.; Yamabe, N.; Shimomura, H.; Shimamura, T.; Ariyada,O.; Wada, M. J. Cryst. Growth 2009, 2009 311 2987-2991. 30. Demirev, P. A. Rapid Commun. Mass Spectrom. 2000, 2000 14(9), 777-781. 31. Petrie,S; Becker, H; Baranov, V. I.: Bohme, D. K. Int. J. Mass Spectrom. Ion Proc. 1995, 1995 145, 79-88. 32. Amorim, J.; Baravian, G.; Touzeau, M.; Jolly, J; J. Appl. Phys. 1994, 1994 76(3), 1487-1493.

33. Asakawa, D.; Kuramochi, A.; Takahashi, E.; Saito, N. Phys. Chem. Chem. Phys., 2018, 20(2), 1082-1090.

ACS Paragon Plus Environment 28