1Department of Chemistry, Nanoscience Center, University of Jyväskylä, P.O. Box 35, University of Jyväskylä, Finland
2Institut für Chemie und Biochemie, Freie Universität Berlin, Takustrasse 3, 14195 Berlin, Germany
3Michael Barber Centre for Collaborative Mass Spectrometry, Manchester Institute of Biotechnology, School of Chemistry, The University of Manchester, Princess Street, Manchester, UK
Pyridinearenes have previously been considered as anion binding hosts due to the electron-poor nature of the pyridine ring. Herein, we demonstrate the encapsulation of Me4N+ cations inside a dimeric hydrogen-bonded pyridinearene capsule, which contradicts with earlier assumptions. The complexation of a cationic guest inside the pyridinearene dimer has been detected and studied by multiple gas-phase techniques, ESI-QTOF-MS, IRMPD, and DT-IMMS experiments, as well as DFT calculations. The comparison of classical resorcinarenes with pyridinearenes by MS and NMR experiments reveals clear differences in their host–guest chemistry and implies that cation encapsulation in pyridinearene is an anion-driven process.
Keywords: cation binding; DFT calculations; ion mobility mass spectrometry; macrocycles; pyridinearenes; resorcinarenes
Resorcinarenes and their derivatives are known for the molecular recognition properties of their self-assembled dimeric and hexameric capsules, which can encapsulate cationic and neutral guests [1-3]. Pyridinearenes  are analogous macrocycles to resorcinarenes. Whereas resorcinarenes are cyclic tetramers of resorcinol, pyridinearenes are formed from 2,6-dihydroxypyridine (see Scheme 1). Although the synthesis of pyridinearenes dates back to 2001 , their host–guest chemistry is still under-explored. Both macrocycles are concave and are known to form capsular assemblies via intermolecular hydrogen bonding [5,6]. Pyridine is significantly less electron-rich than benzene. Consequently, pyridinearene capsules were originally assumed to encapsulate anionic guests inside their cavity due to the π-acidic character of the aromatic walls [7,8]. Mattey et al. detected 1:1 complex formation with PF6− and BF4− by mass spectrometry, however, without ion mobility mass spectrometry, the location of the anion could not be verified and the anions were assumed to interact with the pyridinearene cavity. Inclusion complexes of anions within pyridinearene dimers were also theoretically studied by DFT calculations, but using a truncated pyridinearene dimer model . However, previous studies have also shown that the π-acidic character of pyridine rings, such as 2-oxo-6-oxypyridine, is rather weak [9,10]. Thus, pyridinearenes may be expected to show a dual binding behavior towards anions and cations . In addition to anion complexes also 2:1 complex formation with neutral carboxylic acids and amides have been previously detected by ESI-MS . Very recently, with the help of ion mobility mass spectrometry (IM-MS), we showed that pyridinearenes favor encapsulation of neutral molecules over anionic species and anions are in fact complexed in an exo-position (exclusion complexation) between the lower-rim alkyl chains . A PF6− anion was bound to the lower rim of the pyridinearene dimer via CH–anion and CH–F interactions, while the neutral guest was hosted inside the dimer. Calculated electrostatic potential (ESP) surfaces revealed that the cavity does not possess any significant partial positive potential on the surface of the cavity, except on the N–H hydrogen atoms. They are, however, on tangential positions along the capsules surface, and therefore do not significantly contribute to anion binding . More recently, we were also able to demonstrate that cationic transition metal complexes can template the formation of pyridinearene hexameric capsules in the gas phase .
Here, we report our novel findings on the ability of dimeric tetraisobutylpyridinearene (compound 1 in Scheme 1) to encapsulate cationic guests. Despite of the obvious structural similarities between the dimeric resorcinarene and pyridinearene capsules, we highlight here unique host–guest properties of pyridinearene capsules. In marked contrast to the corresponding resorcinarene capsules, cation binding is clearly feasible, when anions bind in an exo-site and support cation encapsulation by "through-wall" electrostatic interactions.
Results and Discussion
We have previously shown that tetraisobutylpyridinearene forms hydrogen-bonded dimers with eight intermolecular N–H···O(amide) hydrogen bonds in the solid state, in solution and in the gas phase . Resorcinarene capsules of similar size are well-known for their ability to encapsulate small alkylammonium cations inside the dimer, especially quaternary ammonium cations [12-14]. As the cavity sizes of both pyridinearene and resorcinarene dimers are comparable, alkylammonium cations were chosen as the guests for complexation studies with ESI-Q-TOF mass spectrometry. Complex formation was tested with the following series of cationic guests: MeNH3+, Me2NH2+, Me3NH+, Me4N+, EtNH3+, Et2NH2+, Et4N+ and Pr4N+, which were used as the corresponding Cl− or Br− salts. None of the cations MeNH3+, Me2NH2+, Me3NH+, EtNH3+, Et2NH2+ or Pr4N+ formed complexes with 1 and 12, neither using electrospray ionization (ESI) in the positive nor the negative mode. Pr4N+ is certainly too large to fit inside the dimeric capsule, whereas the other non-complexing cations are protic and – as rather strong hydrogen-bond donors – may interfere with intermolecular hydrogen bonds of the dimer. Smaller quaternary cations Me4N+ and Et4N+ were observed to form 1:1 and 2:1 host–guest complexes. With Me4N+, [12 + Me4N]+ was observed in the positive mode, and ions [1 − 2H + Me4N]− and [12 − 2H + Me4N]− in the negative mode (Figure 1). Complexation of Et4N+ was observed only in (−)ESIMS as [1 − 2H + Et4N]− and [12 − 2H + Et4N]− ions. With the chloride and bromide salts, the abundance of cation complexes was modest, but it significantly increased, in both the positive and negative modes, when the counterion was changed to PF6−, BF4− or I−, which have previously been observed to form exo complexes with pyridinearenes . In addition, ternary ion pair complexes such as [12 + Me4N + 2A]− and [12 – H + Me4N + A]− (A = anion, i.e., PF6−, BF4−, or I−) were detected in the negative mode.
The location (endo or exo) of the cation is of structural interest and of vital importance to understand the supramolecular chemistry of pyridinearenes. Interestingly, ternary complexes of 12 containing both cation and solvent are not observed in either of the ESI-MS modes, even though in earlier studies such complexes were observed with solvent and anion . This points to an endo location of the cation. Ion mobility mass spectrometry (IM-MS) is a powerful tool to study structural features such as the endo/exo complexation of supramolecular complexes [11,15-18]. Drift tube ion mobility mass spectrometry (DT-IMMS) enables the determination of structure-related collision cross section (CCS) values without complicated calibration. As depicted in Figure 1, DT-IM-MS revealed that drift times for the ions [12 + H]+, [12 + Me4N]+ and [12 + Na]+ are very similar. In fact, the obtained DTCCSN2 value for [12 + Me4N]+ (385.6 ± 0.40 Å2) is, within the error, the same as for [12 + H]+ (384.8 ± 0.41 Å2) and smaller than for the [12 + Na]+ complex, thus indicating inclusion complexation. This clearly indicates that the Me4N+ complex of 12 has the same rotation average diameter as the protonated dimer and thus the cation is encapsulated inside the cavity. Calculated [20,21] theoretical DTMCCSN2 values indicate that an exo-complexation of the cation should result in a ca. 10 Å2 larger CCS (415.6 Å2 vs 405.7 Å2) value as compared to the endo-complex. Ions [12 − H]− (DTCCSN2 = 387.0 ± 0.44 Å2) and [12 − 2H + Me4N]− (DTCCSN2 = 389.9 ± 0.41 Å2) also exhibit effectively the same DTCCSN2 values and show similar drift times (Figure 1 and Figure S2, Table S1, Supporting Information File 1) in their arrival time distributions (ATDs). Ions [12 + PF6]− and [12 – H + Me4N + PF6]− that presumably carry one exo-complexed anion result in ca. 9 Å2 larger DTCCSN2 values of 395.9 ± 0.50 Å2 and 399.6 ± 0.49 Å2. For [12 + Me4N + 2PF6]− with two attached counter ions, the DTCCSN2 value is even larger (5–8 Å2 larger compared to ions with one anion). This clearly demonstrates the anions to be located in the outer periphery, while the cation is encapsulated inside the dimer.
The difference in endo/exo complexation is clear in IM-MS, when the larger PF6− anion is used in the experiment. For complexes with smaller anions BF4− and I− the structural conclusions are more difficult to draw based only on CCS values, as the difference between the endo and exo complexes is smaller. To obtain complementary information, the ternary complex ions [12 + Me4N + 2A]− (A = PF6−, BF4−, I−) were further investigated by infrared multiphoton dissociation (IRMPD). IRMPD is an MS/MS technique, which can be utilized to study the fragmentation of supramolecular complexes . In these experiments, the [12 + Me4N + 2A]− ions were mass-selected and irradiated for 20 to 250 ms with a laser power of 95% of 25 W. At short irradiation times, the main dissociation product was [12 + Me4N − 2H]−, which is produced after elimination of two molecules of HA (HA = H+A−, Schematic presentation of dissociation in Figure S4, Supporting Information File 1). This ion further dissociates to a neutral monomer and to the [1 − 2H + Me4N]− ion, which can only result from two HA eliminations and dissociation of the H-bonded capsule. Elimination of two HA molecules and formation of [1 + Me4N − 2H]− would be unexpected from a complex with exo-complexed cation. Minor fragments resulting from a direct elimination of an ion pair are observed and [12 + A]− appears in spectra. This likely results from a relocation of the anion to a position close to the hydrogen-bonded seam of the dimer, which is then partially opened to release the ion pair (pathway 2 in Figure S4, Supporting Information File 1). As all studied ternary complexes behaved similarly in the experiments, it can be stated that Me4N+ is located inside the cavity, while the PF6−, BF4− and I− anions reside in the exo-binding site.
The binding efficiency of pyridinearene towards the Me4N+ cation was compared to that of resorcinarene, which is known to bind small cations with high affinity. 1, 2 and Me4NPF6 were mixed in 1:1:1 ratio and measured by ESI-QTOF-MS. As seen in Figure 2, the resorcinarene dimer [22 + Me4N]+ was detected as the base peak in the positive mode, while [12 + Me4N]+ was hardly observed. In addition to the pure dimers, the formation of the heterodimer [1·2 + Me4N]+ was also observed. In the negative mode, heterodimers [1·2 − H]− and [1·2 − 2H + Me4N]− were observed with higher intensity than the corresponding homodimeric capsules of 1 or 2 (Figure S1, Supporting Information File 1). This is surprising due to differences in H-bonding geometry between dimers of 1 or 2. Previously in case of hexameric capsules heterohexamers of pyridinearene and resorcinarene were not observed . The PF6− anion was observed to be bound to the pyridinearene dimer in [12 + PF6]− and the heterodimer [1·2 + PF6]−, but interestingly the abundance of [1·2 + PF6]− was half of [12 + PF6]−, and PF6− complexes with 2 or 22 had even lower abundances. This clearly shows that the anion has a higher affinity to pyridinearene than to resorcinarene. Inclusion complexation of solvent was also only detected with the pyridinearene homodimer (for example as the [12 + MeCN + PF6]− ion). This mixed host experiment shows the resorcinarene affinity towards Me4N+ to be higher than that of pyridinearene. Also, the lack of anion and/or solvent complex formation with 2 or 22 in the gas phase indicates that this binding feature is unique to pyridinearenes and it does not take place with resorcinarene. The high abundance of ternary ion pair complexes with pyridinearene, but not with resorcinarene indicates that the cation complexation process to pyridinearene is anion driven at least to some extent. The reduced electron density in the pyridine rings is compensated by favorable "through-wall" electrostatic interactions between the exo-anion(s) and the endo-cation.
In solution, by 1H NMR, complexation of Me4N+ cation was observed with 2, but not with 1 (Figures S5 and S6, Supporting Information File 1). In a sample of 10 mM 2 and Me4NPF6 (1:3 ratio in CDCl3/CD3CN 4:3 v:v mixture) a 0.35 ppm upfield shift was observed for the Me4N+ signal. In a sample of 1 and Me4NPF6 no shift is observed. This clearly indicates that also in solution 1 does not possess a similar affinity towards the Me4N+ cation as 2. Also, all our attempts to obtain a solid-state single crystal structure of 1 with a cationic guest were unsuccessful. It is possible that the observation of such complexes requires special conditions present in the ESI source.
To obtain more detailed insights into the observed gas-phase structures, DFT calculations were carried out for 1, 1OH, 12, [12 + Me4Nendo]+, [12 + Me4Nexo1]+, [12 + Me4Nexo2]+, [12 + Me4Nendo + Iendo], [12 + Me4Nendo + Iexo], [12 + Me4Nendo + 2Iexo]−, [12 + Iexo]−, and [12 + 2Iexo]2−(Figure S7, Supporting Information File 1) without (PBE0/def2-TZVP) and with (PBE0-D3/def2-TZVP) dispersion correction because in the gas phase dispersion plays an important role in the formation of supramolecular complexes [24-30]. Conformational analysis was first performed for both tautomers of the monomer to ensure that further calculations were carried out for the most stable tautomer. The conformational analysis revealed that the pyridone tautomer (1) is ≈150 kJ/mol lower in energy than the dihydroxy tautomer (1OH). The geometry optimizations showed that the geometry of [12 + Me4Nendo]+ is similar without and with dispersion correction (hydrogen bond lengths are listed in Table S2, Supporting Information File 1). This is expected, because the standard hybrid functionals are able to describe strong hydrogen bonds such as NH···O and OH···O in a reasonable manner [31,32]. Due to the encapsulated guest molecule, the [12 + Me4Nendo]+ complex has slightly longer NH···O and OH···O bonds than empty 12. The addition of I− anion(s) at the lower rim(s) of [12 + Me4Nendo]+ has hardly any influence on the length of NH···O and OH···O bonds as exemplified by the hydrogen bond lengths of [12 + Me4Nendo + Iexo] and [12 + Me4Nendo + 2Iexo]−. The same trend is also observed for cation-free complexes ([12 + Iexo]− and [12 + 2Iexo]2−) if their NH···O and OH···O bond lengths are compared to those in 12.
Calculations were also carried out for complexes, where exo binding of the cation was considered. The geometry optimization yields two different geometries, [12 + Me4Nexo1]+ and [12 + Me4Nexo2]+ (Figure S7, Supporting Information File 1). In [12 + Me4Nexo1]+, the Me4N+ cation interacts with the lower rim isobutyl groups only when the dispersion interaction is included in calculations, whereas at the PBE0/def2-TZVP level [12 + Me4Nexo1]+ is not a stable minimum on the potential energy surface and results in an additional exo complex, namely [12 + Me4Nexo2]+, in which the Me4N+ cation resides in proximity to the seam of hydrogen bonds. The NH···O and OH···O bond distances for [12 + Me4Nexo1]+ are almost identical with 12 indicating that the coordination of the Me4N+ cation at the lower rim has only a minor effect on the hydrogen-bond network of 12. To verify complexation of separate ion pairs, the geometry optimization was carried out also for [12 + Me4Nendo + Iendo] , which showed that the cavity of 12 is too small for the simultaneous complexation of anionic and cationic guests, resulting in a partial rupture of the hydrogen-bonding seam and a ≈100 kJ/mol weaker interaction energy compared to [12 + Me4Nendo + Iexo] (Figure S7 and Table S2, Supporting Information File 1).
To illustrate the unsuitability of the lower rim for binding cationic guests, we mapped the ESP surfaces for 1, 12, [12 + Me4Nendo]+, [12 + Me4Nexo1]+, and [12 + Me4Nexo2]+. Although the ESP surface of 1 show some π-acidic character for the cavity, as depicted in Figure 3, the electron-poor areas of 12, [12 + Me4Nendo]+, [12 + Me4Nexo1]+, and [12 + Me4Nexo2]+ (see also Figures S9 and S10, Supporting Information File 1) are more concentrated on the lower rim isobutyl groups than on the seam of the hydrogen bonds that actually contain electron-rich oxygen atoms able to interact with the Me4N+ cation through ion–dipole interactions. In the gas phase and in the absence of the dispersion force, the positively charged Me4N+ cation, therefore, prefers interactions with the electron-rich oxygen atoms of the hydrogen-bond network to the positively charged lower rim isobutyl groups.
The calculated interaction energies are listed in Table S2 (Supporting Information File 1). When the dispersion correction and/or the counterion(s) are taken into account, the energy analysis supports the findings obtained from ESI-MS studies, geometry optimizations and ESP surfaces: The Me4N+ cation is most likely bound inside the cavity of 12. It is also important to note that the calculated dispersion-corrected interaction energy of [12 + Me4Nexo1]+ (−37 kJ/mol) is roughly two fifth of the interaction energy of [12 + Me4Nexo2]+ (−98 kJ/mol) because there is no favorable ion–dipole interaction between 12 and Me4N+ in [12 + Me4Nexo1]+. Even though the encapsulation of the Me4N+ cation is a favorable process already even without anions, the calculated interaction energies show that the complexation of cationic guest can be enhanced by the additional exo-complexation of anions to the lower rims. The interaction energies of [12 + Me4Nendo]+, [12 + Me4Nendo + Iexo], and [12 + Me4Nendo + 2Iexo]−, increases from −152 kJ/mol to −342 kJ/mol, and −521 kJ/mol, respectively, when the amount of I− anions is increased in the complex. Moreover, the calculated interaction energies illustrate well that the dispersion interaction has considerable influence on the interaction energies of [12 + Me4Nendo]+ (−126.4 kJ/mol), [12 + Me4Nexo1]+ (−37.5 kJ/mol), [12 + Me4Nexo2]+ (−23.2 kJ/mol), [12 + Me4Nendo + Iexo] (−130.5 kJ/mol), and [12 + Me4Nendo + 2Iexo]− (−133.7 kJ/mol), although their geometries are similar at both levels of theory. This means that pyridinearene dimers and their complexes are already sufficiently large supramolecular systems in which the omnipresent dispersion interaction can add up to a substantial force due to multiple interaction sites.
Calculations were also carried out for 2, 22, [22 + Me4Nendo]+, [22 + Iexo]−, and [22 + 2Iexo]2−, at the PBE0-D3/def2-TZVP level of theory (Figure S8, Supporting Information File 1) to explain the observed differences between 1 and 2 in the experimental ESI-MS studies. If the interaction energies of [22 + Me4Nendo]+ (−220 kJ/mol) and [12 + Me4Nendo]+ (−152 kJ/mol) are compared it is evident that 22 is a better host for the Me4N+ cation than 12. On the other hand, when comparing the interaction energies of complexes [22 + Iexo]− (−149 kJ/mol), [22 + 2Iexo]2− (−170 kJ/mol), [12 + Iexo]− (−175 kJ/mol), and [12 + 2Iexo]2− (−223 kJ/mol), it can be stated that 12 has a stronger affinity towards anions than 22. These results are fully in line with the MS studies and calculated ESP surfaces. By comparing the calculated ESP surfaces of 1 and 2 in Figure 3, it is a clear that the π-basic character of the cavity is more pronounced in 2 than in 1. The investigation of ESP surfaces of 1, 12, 2, and 22 also reveal that the lower rims of pyridinearenes 1 and 12 are more electron deficient than 2 and 22, explaining the weaker affinity of 2 and 22 towards anions. This result is also underpinned by the calculated interaction energies of anionic complexes; −149 kJ/mol ([22 + Iexo]−) vs −175 kJ/mol ([12 + Iexo]−) and −170 kJ/mol ([22 + 2Iexo]2−) vs −223 kJ/mol ([12 + 2Iexo]2−) that show a stronger affinity of 12 towards anions compared to 22.
In conclusion, we have shown that the pyridinearene dimer, which until recently was considered to be an anion receptor, can bind also cationic guests. The Me4N+ cation was observed to bind to the pyridinearene dimer, and all mass spectrometric data and theoretical calculations show undoubtedly that the cation is located inside the cavity of the dimer in the gas phase. It is an interesting fact, that the binding properties of pyridinearene differ from earlier reports [7,8]. However, cation binding to pyridinearene is clearly not as strong as with resorcinarene, which is known for its excellent cation receptor properties. A comparison of these two macrocyclic hosts reveals significant differences in their binding properties. Pyridinearene appears to have a better affinity towards neutral guests and has also a higher affinity to complex anions at the exo-binding sites at the lower rim, whereas resorcinarene clearly has a higher affinity towards cations. However, the pyridinearenes' ability to form exo-complexes with anions can assist its ability to bind cations in endo-positions.
Compounds 1 and 2 have been prepared according to the reported procedures [11,33,34]. Salts were commercially available (From Aldrich, Fluka and TCI) and used as received. Mass spectrometric experiments have been performed with an ABSciex QSTAR Elite ESI-Q-TOF mass spectrometer. Ion mobility mass spectrometry (IM-MS) experiments were conducted with a Waters Synapt G2 equipped with a linear drift cell and Agilent 6560 Ion mobility Q-TOF mass spectrometer. All CCS values were obtained using nitrogen as a drift gas and stepped field methods. IRMPD experiments were performed with an Ionspec QFT-7 ESI-FT-ICR with a 7 T superconducting magnet. Samples for all mass spectrometric experiments were prepared with 10 or 20 µM concentration and 1:3 host–guest ratio in acetone. Theoretical CCS values were calculated using IMoS [20,21]. NMR experiments were performed with a Bruker Avance III HD 300 NMR spectrometer. Samples were prepared at 10 mM concentration and 1:3 host–guest ratio in CDCl3/CD3CN (4:3, v:v) mixture. DFT calculations were performed by Spartan’ 16 and Gaussian 09 (D01) software packages. More detailed information of the experiments and parametrization can be found in Supporting Information File 1.
Supporting Information File 1:
Experimental details and supplementary information.
We wish to thank The Academy of Finland (projects 284562, 278743 and 312514 to E.K. and 285855 to J.O.M.), the University of Jyväskylä and the Deutsche Forschungsgemeinschaft (Core Facility BioSupraMol) for financial support. J.O.M wishes to acknowledge Finnish Grid and Cloud Infrastructure (persistent identifier urn:nbn:fi:research-infras-2016072533), CSC - IT Center for Science, Finland, and Prof. Heikki M. Tuononen for computational resources. A. Kruve is grateful for a postdoctoral fellowship from the Alexander von Humboldt Foundation. Prof. Magdalena Zimnicka is gratefully acknowledged for her assistance and helpful discussions regarding the determination of theoretical CCS values.
Lee, T.-C.; Kalenius, E.; Lazar, A. I.; Assaf, K. I.; Kuhnert, N.; Grün, C. H.; Jänis, J.; Scherman, O. A.; Nau, W. M. Nat. Chem.2013,5, 376–382. doi:10.1038/nchem.1618 Return to citation in text: 
Dearden, D. V.; Ferrell, T. A.; Asplund, M. C.; Zilch, L. W.; Julian, R. R.; Jarrold, M. F. J. Phys. Chem. A2009,113, 989–997. doi:10.1021/jp808771n Return to citation in text: 
Jurček, O.; Bonakdarzadeh, P.; Kalenius, E.; Linnanto, J. M.; Groessl, M.; Knochenmuss, R.; Ihalainen, J. A.; Rissanen, K. Angew. Chem., Int. Ed.2015,54, 15462–15467. doi:10.1002/anie.201506539 Return to citation in text: 
Gabelica, V.; Shvartsburg, A. A.; Afonso, C.; Barran, P.; Benesch, J. L.; Bleiholder, C.; Bowers, M. T.; Bilbao, A.; Bush, M. F.; Campbell, J. L.; Campuzano, I. D.; Causon, T.; Clowers, B. H.; Creaser, C. S.; De Pauw, E.; Far, J.; FernandezLima, F.; Fjeldsted, J. C.; Giles, K.; Groessl, M.; Hogan, C. J.; Hann, S.; Kim, H. I.; Kurulugama, R. T.; May, J. C.; McLean, J. A.; Pagel, K.; Richardson, K.; Ridgeway, M. E.; Rosu, F.; Sobott, F.; Thalassinos, K.; Valentine, S. J.; Wyttenbach, T. Mass Spectrom. Rev.2019,38, 291–320. doi:10.1002/mas.21585 Return to citation in text: 
Tunstad, L. M.; Tucker, J. A.; Dalcanale, E.; Weiser, J.; Bryant, J. A.; Sherman, J. C.; Helgeson, R. C.; Knobler, C. B.; Cram, D. J. J. Org. Chem.1989,54, 1305–1312. doi:10.1021/jo00267a015 Return to citation in text: 
Gabelica, V.; Shvartsburg, A. A.; Afonso, C.; Barran, P.; Benesch, J. L.; Bleiholder, C.; Bowers, M. T.; Bilbao, A.; Bush, M. F.; Campbell, J. L.; Campuzano, I. D.; Causon, T.; Clowers, B. H.; Creaser, C. S.; De Pauw, E.; Far, J.; FernandezLima, F.; Fjeldsted, J. C.; Giles, K.; Groessl, M.; Hogan, C. J.; Hann, S.; Kim, H. I.; Kurulugama, R. T.; May, J. C.; McLean, J. A.; Pagel, K.; Richardson, K.; Ridgeway, M. E.; Rosu, F.; Sobott, F.; Thalassinos, K.; Valentine, S. J.; Wyttenbach, T. Mass Spectrom. Rev.2019,38, 291–320. doi:10.1002/mas.21585