Tautomeric equilibrium, proton affinity and mass spectrometry fragmentation of flexible hydrogen-bonded precursors and rigid \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {N}\longrightarrow \hbox {BF}_2$$\end{document}N⟶BF2 fluorescent dyes

The stability of two groups of conformationally locked molecules, similar in topology, but differing only by the type of the bridge rigidifying their structure, is studied. The series of the less-rigid 2-phenacylheterocyclic compounds and their stiff difluoroboranyl derivatives are investigated for the determination of the effect of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {NCH}_3$$\end{document}NCH3/S/O replacement in a five-membered heterocyclic ring and the presence of a strong electron-donating group on the tautomeric equilibrium, protonation affinity, and fragmentation pattern observed in the structural elucidation by means of mass spectrometry technique. The results of the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\omega $$\end{document}ωB97X-D/6-311++G(d,p) calculations, the topological analysis of electron density as well as the experimental MS measurements show the importance of the number of heteroatoms, their properties, and location in the molecule for the rational design of the systems of desired stable tautomers or the favorable protonation sites. The obtained data allow for the understanding of the fundamentals of the novel highly fluorescent difluoroborates fragmentation behavior, vital for their structural elucidation with the application of high-resolution tandem mass spectrometry methods.

Tautomeric equilibrium, proton affinity and mass spectrometry fragmentation of flexible hydrogen-bonded precursors and rigid N −→ BF 2 fluorescent dyes Małgorzata A. Kaczorowska 1

, Anna Kaczmarek-Kędziera 2 & Borys Ośmiałowski 2*
The stability of two groups of conformationally locked molecules, similar in topology, but differing only by the type of the bridge rigidifying their structure, is studied. The series of the less-rigid 2-phenacylheterocyclic compounds and their stiff difluoroboranyl derivatives are investigated for the determination of the effect of NCH 3 /S/O replacement in a five-membered heterocyclic ring and the presence of a strong electron-donating group on the tautomeric equilibrium, protonation affinity, and fragmentation pattern observed in the structural elucidation by means of mass spectrometry technique. The results of the ωB97X-D/6-311++G(d,p) calculations, the topological analysis of electron density as well as the experimental MS measurements show the importance of the number of heteroatoms, their properties, and location in the molecule for the rational design of the systems of desired stable tautomers or the favorable protonation sites. The obtained data allow for the understanding of the fundamentals of the novel highly fluorescent difluoroborates fragmentation behavior, vital for their structural elucidation with the application of high-resolution tandem mass spectrometry methods.
Dyes containing a disubstituted boron center constitute a well-known class of fluorophores, applied in various branches of science and technology, such as biological labeling 1,2 , construction of new drug carriers, anticancer or antimicrobial photosensitizing agents 3,4 or dye-sensitized solar cells 5,6 . They owe the wide popularity to their extraordinary properties, namely large fluorescence quantum yield and easy chemical modifications. The rational design, and further, the precise tuning of these luminophore characteristics is vital and can be achieved as a consequence of the good understanding of the relationships between the molecular structure and their optical features [7][8][9][10] . However, even small structural modification can significantly affect the functional properties and, consequently, the applicability of dyes 11 . Moreover, the advanced process of the molecular design is usually multidimensional and influence of several factors at the same time can undergo the synergistic or anti-cooperative combination. This fact rationalizes the requirement of application of computational chemistry tools in the initial step of the procedure.
In order to verify the properties of the theoretically designed fluorophores, the synthesis of the selected species can be performed. Thus, in the context of the present study, it is also crucial to control precisely the difluoroboranyl derivatives synthesis processes, which often requires confirmation of the structure, properties and-what follows-the reactivity of the precursors and intermediate products. The precursors for these system synthesis are known to exhibit relatively high flexibility due to the presence of several single carbon-carbon bonds, enabling free rotation. For a group of the investigated precursors, such a rotation can strongly affect their rigidity by shifting the hydrogen bond donor towards the hydrogen bond acceptor. This transformation can lead to the formation of the intramolecular hydrogen bond which considerably stabilizes the molecule. In the Experimental and theoretical methods Computational methodology. Full geometry optimization has been performed for all of the investigated species with a long-range corrected hybrid functional containing an empirical Grimme D2 dispersion correction 30 , ωB97X-D 31 , together with a conventional Pople basis set augmented with diffuse functions, 6-311++G(d,p) [32][33][34] , in vacuum. This has been proven to provide a safe general-purpose approach for prediction of molecular geometries, tautomeric equilibria and reaction barrier heights and delivers the data for comparison with the MS experiment [35][36][37] . The tautomeric equilibria were also investigated with two other multipurpose hybrid functionals, namely the most popular three-parameter B3LYP 38 one and the Truhlar M06-2X 39 Minnesota functional containing 54% of exact exchange, recommended for the main group thermochemistry and non-covalent interactions 39,40 . As the wave-function formalism reference, also the Hartree-Fock 41 and MP2 [42][43][44][45]   www.nature.com/scientificreports/ results are provided in SI. They are compared with the double hybrid functional B2PLYP containing the combination of the exact HF exchange and MP2-like correlation contribution, additionally supplemented with the empirical Grimme dispersion D3 correction. Such functionals are perceived as the significant improvement upon the conventional hybrid functionals in context of dispersion intramolecular contribution description for relative energies or weak intermolecular interactions [46][47][48] . Additionally, as the reference, the single-point energy  calculations for the investigated tautomers has been performed within the high-level local correlation approach  in the domain-based local pair natural orbitals coupled cluster version, namely DLPNO-CCSD(T), providing  linear scalling vital for high accuracy calculations for large systems [49][50][51][52][53] . The B2PLYP-D/6-311++G(d,p) and DLPNO-CCSD(T)/aug-cc-pVTZ results are presented in Fig. S4. The character of the stationary points has been confirmed with harmonic vibrational analysis. The reaction paths for proton transfer have been verified with the intrinsic reaction coordinate calculations.
Proton affinities (PA) and gas phase basicities (GB) have been estimated for the model protonation reaction, acid [X+H] + −→ base [X] + H + , as the enthalpy and Gibbs free energy components, respectively: where X denotes one of the analyzed systems, H is the enthalpy change, calculated as the sum of electronic energy and thermal correction to enthalpy at 298.15 K and G is a Gibbs free energy change. The enthalpy of the proton at 298.15 K was assumed to be equal to 5/2RT = 0.002368 Hartree, and its Gibbs free energy − 0.0100 Hartree 54 . It has been previously shown, that the time-efficient DFT approaches provide the reliable values of proton affinity for a wide spectrum of organic moieties [54][55][56][57][58] , therefore the ωB97X-D/6-311++G(d,p) approach has been applied for PAs and GBs as well. Corresponding Wiberg bond indexes are presented. The calculations have been performed in Gaussian16 and Gaussian09 packages 59,60 and the DLPNO-CCSD(T) calculations have been carried out in ORCA 4.2.1 61,62 .
In order to shed more light on the protonation and fragmentation processes in the investigated systems, the Atom-In-Molecules analysis has been performed within AIMAll program [63][64][65] . The obtained electron density ρ , laplacian of electron density 2 ρ , total energy density H b 65 and the estimation of the energy of the intermolecular hydrogen bonds calculated according to Espinosa et al. 66 as E HB = 0.5· V, where V is the potential energy density in bond critical point, are summarized in Supplementary Information.

Chemicals.
All the reagents and solvents necessary for the synthesis were purchased from Sigma-Aldrich and were used without further purification. The methods applied for the synthesis of 2-phenacylheterocycles and their difluoroboranyl derivatives (shown at Fig. 1) have been described in detail elsewhere 10 .
Mass spectrometry experiments. The Q-Exactive Orbitrap mass spectrometer (Thermo Fisher Scientific, Bremen, Germany) has been used to perform all of the mass spectrometry experiments. The examined compounds were dissolved in methanol and formic acid solution (98:2, vol/vol) before being introduced into the mass spectrometer using an electrospray source (Advion BioSciences ltd., Ithaca, NY, USA). Data acquisition and analysis were performed using the Thermo Xcalibur software (ver. 4.1.31.9). All of the mass spectra were obtained in positive ion mode within the m/z range of 50-750, at the resolution of 140,000 (m/z 200). MS/MS experiments were carried out with the use of higher energy collisional dissociation mode (HCD) with normalized collision energy set individually for each compound (usually in a range of 15-45 eV). The resulting mass spectra were analyzed manually.

Results and discussion
Tautomeric equilibria. The initial part of the computational study is devoted to the tautomeric equilibria in neutral 1-6 species and is focused on the determination of the stability order of the tautomers and thus presume about the composition of the reaction mixture upon the standard conditions. In the present study, 1-3 series (R=H) allow for the analysis of the effect of the nature of the heterocyclic five-membered ring on the tautomeric equilibrium, while 4-6 (R=N(CH 3 ) 2 ) gives the information about the influence of the electron-donating dimethylamino group. On the other hand, the protonation of the terminal dimethylamino group, which occurs spontaneously for instance in the process of ESI-HRMS analysis, provides the strong electron-withdrawing effect on the molecule under investigation. Therefore, Fig. 2 presents the relative energies of keto (K), enol (O) and enamine (E) tautomers of 1-6 and [4+H] + , [5+H] + and [6+H] + systems with protonated dimethylamino group. The full set of relative energies and relative Gibbs free energies for the three tautomers of nine systems obtained in different methods are provided in Supplementary Information as histograms in Fig. S4 and S5. One should also notice that for the comparison for all of the investigated structures the alike keto tautomers have been selected-all of them are significantly bend and possess the carbonyl group pointing up to the direction of the N(CH 3 )/S/O atom in heterocyclic ring. Consequently, for all of the analyzed phenacylheterocycles and their N(CH 3 ) 2 -derivatives this is the lowest energy keto tautomer. Additionally, at least two other rotamers of the keto form can be determined: with carbonyl group pointing downwards in the direction of the heterocyclic amine/ imine nitrogen and in carbonyl placed in the plane perpendicular to the heterocyclic ring. All of these keto rotamers for 1-6 lie in the 5 kcal/mol energy interval. Enol and enamine tautomers for all of the 2-phenacylheterocycles and their dimethylaminated derivatives are considered as planar forms stabilized with the intramolecular www.nature.com/scientificreports/ Several general trends can be noticed. First, for R=H (systems 1-3) the observed energy differences are smaller than 2.3 kcal/mol and for imidazole and thiazole remain below the accuracy of the applied computational methodology. Therefore, one can assume the easy N7-to-O13 proton transfer and significant abundance of all tautomers in the reaction mixture. The substitution of the investigated species with an electron-donating terminal N(CH 3 ) 2 group leads to the additional stabilization of the enol form in the case of 4 and 6. The protonation of the terminal N(CH 3 ) 2 group significantly destabilizes the keto tautomer in all three investigated systems, thus causing its smallest abundance. Nevertheless, the obtained energy differences among tautomers are in most cases small enough to remain inconclusive. The qualitative tendencies of the ωB97X-D/6-311++G(d,p) calculations are confirmed by the B2PLYP-D/6-311++G(d,p) and DLPNO-CCSD(T)/aug-cc-pVTZ data (Fig. S4). www.nature.com/scientificreports/ The general remark that can be drawn from these results is that the reaction mixtures in standard conditions contain multitude of low energy isomers for all analyzed systems. However, it should be underlined that due to the MS protocol applied, ab initio calculations presented here are performed in vacuum in order to achieve the closest correspondence between the experimental MS setup and computational study. This, however, can significantly alter the results with respect to other experimental measurements carried out in solution. Therefore, the direct comparison of the experimental data with the present results is in principle not straightforward taking into account comparisons between MS and NMR techniques. It should be mentioned that on the basis of the 1 H and 13 C NMR measurements in CDCl 3 and computational approach 67 of a series of similar 1-methyl-2-phenacylbenzimidazole systems, enamine tautomer is the least stable and therefore not present in the reaction mixture. On the other hand keto has been found to remain the most stable tautomer independently on the substituent present in the phenyl ring, and enol form have been shown to undergo an additional stabilizing effect from the intramolecular N . . . H-O interaction. The K and O energetic order is reversed in the case of corresponding oxazole derivatives 68 .

2-phenacylheterocycles.
Proton affinity of 1-3. Proton affinity values for the 2-phenacylheterocyclic compounds as well as their -N(CH 3 ) 2 -substituted derivatives are presented in Fig. 3. Additionally, the gas-phase basicity for all these systems is provided in Table S10. Both of these notions, as well as the simple energetic criteria, provide the alike informations for the preferred protonation sites in 1-3. The generally observed tendencies point out that in the neutral systems, the sites most prone to protonation are the N7 atom in heterocyclic ring and the C10 atom. The only system which seem to favour the O13 protonation is the enamine form of imidazole. However, one should also notice, that the obtained PA and GB values can be very similar for several protonation sites (for instance PA equal to 239.75 and 242.20 kcal/mol for enolic form of 1 at N7 and C10 sites, respectively), thus not allowing for the definite conclusions about the single prevailing protonated form and-what followsindicating the probable presence of several protonated forms of similar stability in the experimental conditions. The more detailed discussion of the numerical values of the PA and GB for 1-3, the Wiberg bond indexes and feasible fragmentation mechanisms is provided in SI.
In order to shed the light on the fragmentation pattern, the computational analysis of the mechanism of the possible proton transfer pathways in [1+H] + has been performed and depicted in Fig. 4. The estimated energy barriers demonstrate that starting from the lowest-energy 1.3 keto tautomer of 1 protonated at the C10 site, the process in two directions can proceed: the lower energy barrier lead to the intramolecular rotation and later on, the C10-O13 proton transfer to the enamine form 1.1, with a barrier of 50.64 kcal/mol. The other possibility is the generation of the protonated enolic form 1.5, which however requires crossing of the high (over 100 kcal/ mol) overall energy barriers. This would be surely impossible for the spontaneous process in standard conditions, nevertheless taking into account the conditions under which HCD experiments are carried out, it remains viable during the fragmentation process.
The Atom-In-Molecules analysis allows the estimation of the intramolecular hydrogen bond energy, E HB , according to Espinosa et al. 66 The corresponding data for enol and enamine forms of 1-6 and their protonated counterparts are summarized in Table S13 and on histograms in Fig. S9. One can clearly notice the substantial differences in the strength of this interaction depending on the system and its tautomer. Enol form is significantly favoured, exhibiting the strongest stabilizing effect of H-bond, however its strength diminishes in the sequence of N(CH 3 )/S/O in heterocyclic ring in both series (namely pristine 1-3 and dimethylaminated 4-6). Similar dependence on the heterocyclic ring composition is presented by the intramolecular hydrogen bond in the enolic form protonated at the terminal dimethylamine group.
Fragmentation process of 1-3. The results of the DFT calculations show that in the case of compounds 1-3, heterocyclic nitrogen atoms (-N=) are the privileged protonation sites (compare Fig. 5, panel (a)). The weakest bonds according to the calculated Wiberg bond indexes (see Table S11), most susceptible for breaking upon the MS collisions, include the detachment of the methyl group from the imidazole ring (N9-C30), and (depending on the protonation site) the C10-C11 or C11-O13 bonds. The ESI HCD MS/MS mass spectra of the three ions [1+H] + , [2+H] + and [3+H] + are dominated by signals corresponding to singly charged parent ions and to the [ C 7 H 5 O] + species (Fig. 5, panel (b)).
In -CH=O] + ions would involve the complex intramolecular rearrangements, the detachment of hydrogen atoms from the benzene rings, cleavage of several bonds and is unlikely, whereas the mechanism leading to the formation of ArCO + products requires the transfer of only one proton and simple intramolecular rearrangement (in case of protonated keto tautomeric form of singly charged ions of 1 this mechanism is even simpler). Furthermore, also the analysis of previous results of gas-phase fragmentation behavior of similar structure compounds leads to the conclusion that formation of [(C 6 H 5 )-C=O] + is most likely 69 .
The influence of the properties of the -N(CH 3 ) group and single elements S and O, located in the fivemembered ring on the higher energy collisional dissociation behavior of singly charged ions of 1-3 is observed in the case of minor HCD fragmentation processes. They are discussed in details in Supporting Information.

Dimethylamino derivatives of 2-phenacylheterocycles. Proton affinitty of 4-6. The substitution
of the BF 2 -containing fluorophores with the terminal N(CH 3 ) 2 group has been shown to drastically affect their photophysical properties due to the more pronounced charge transfer character of the S 0 −→ S 1 transition 10 .
In the case of their precursors however, although modifying the electron distribution within the molecule by www.nature.com/scientificreports/ dimethyloamine group can be observed, it seems to not influence much the protonation affinity and thus the ESI HCD MS/MS fragmentation pattern in the gas phase. Nevertheless, in most cases, the PA as well as the GB values seem to be slightly higher for the N(CH 3 ) 2 -substituted 2-phenacylheterocycles than for their pristine analogs, what remains in agreement with chemical intuition (compare Fig. 3 and numerical values of GB in Table S10.). Taking into account the distance between the heteroatoms from the molecular core (N7, N/S/O9 and O13) and in the terminal dimethylamino group, the intramolecular proton transfer between N26 and any other protonation site in the system arises to be improbable. Yet, the MS technique applied in the present study, provides the source of the external protons, which can be attached at any point. Therefore, two cases of the proton transfer process are considered here and they are presented as panel (a) and panel (b) of Fig. 6, for a mechanism not involving the terminal N(CH 3 ) 2 group and the other one, in which originally dimethylamino substituent is www.nature.com/scientificreports/ protonated. In the first path, as it was observed for 1, the diprotonated hydroxyl group in 4.1 is the least stable and requires the highest energy barrier crossing. On the other hand, the proton transfer from the C10 of initial, most stable 4.3 keto form to O13, preceeded by the intramolcular rotation, delivers another relatively stable 4.5 isomer, which could possibly be observed in noticeable amounts even in standard conditions. Similar processes involving the proton transfer in the central part of the molecule, between the heterocyclic ring and enol group, for the system originally protonated at the terminal dimethylamino substituent, demand crossing of the significantly higher energy barrier (of the order of 60 kcal/mol) than in the case when both protons are placed close to each other. Still, taking into account the collision energy of MS/MS measurements, these processes cannot be definitely excluded from consideration of fragmentation pattern.
Fragmentation process of 4-6. The introduction of the terminal dimethylamine group to the 2-phenacylheterocycle adds another possible protonation site, which occurs to be competitive with respect the C10, O13 or N/S/O9 atoms (Fig. 3, Table S10). The calculated Wiberg bond indexes (Tab. S14) indicate that in addition to the weakest points in the systems observed for 1-3 series, the detachment of CH 3 groups from the terminal dimethylamino substituent could be considered. In the case of the protonation of the nitrogen from the terminal group, also the detachment of the whole -N(CH 3 ) 2 substituent becomes plausible since the corresponding Wiberg bond indexes for C16-N26 bond remain smaller than 1.0. The results of the performed tandem mass spectrometry experiments (Fig. 7) clearly demonstrate that the main dissociation processes of dimethylamino derivatives are similar to that found in the fragmentation of unmodified 2-phenacylheterocycles, and are related to protonation of heterocyclic ring nitrogen atoms. The HCD MS/MS spectra of singly charged ions of compounds 4-6 are dominated by signals which can be assigned to [ C 9 H 10 NO] + ([(CH 3 ) 2 N(C 6 H 4 )-C=O] + ) species, corresponding to [ C 7 H 5 O] + ions formed in the HCD of protonated molecules 1-3.
The higher energy collisional dissociation of all ions containing -N(CH 3 ) 2 group is also associated with protonation of nitrogen atoms from these moieties and as a consequence elimination of the entire or part of a dimethylamino substituent occurs (i.e., [4- 3 +H] +· ). Furthermore, the presence of the dimethylamino group affects other minor fragmentation processes, because the terminal -N(CH 3 ) 2 modifies the electron distribution within the molecules. This effect strongly depends on the elemental composition of the heterocyclic moiety or, more precisely, on the properties of atoms/groups that are present in the ring (N(CH 3 ), S, O) and is discussed in details in Supplementary Information.
Difluoroborates. Proton affinity of 7-9. Due to their rigidity, difluoroborate systems analyzed in the present study appear computationally less demanding. Therefore only one neutral form for 7, 8 and 9 has been determined and consequently protonated by the ESI proton in all possible protonation sites in order to establish the energetic preferences and favorable proton affinities for fragmentation process justification. It has been undoubtedly shown that the highest proton affinity is found for a C10 protonation in all systems (234.32 kcal/ mol for 7, 229.49 kcal/mol for 8 and 227.26 kcal/mol for 9). In the case of imidazole (7) this leads to the lowest energy protonated isomer 7.1 (see Fig. 8). Thereafter, the proton can be transferred to the O13 site (7.2) crossing the barrier of 69.65 kcal/mol and generating the structure of the relative energy equal to 21.58 kcal/mol. The further proton transfer can occur to the N7 site via the smaller barrier of 36.08 kcal/mol to the 7.3 isomer. The proton relocation to N9 site requires the energy cost of 67.57 kcal/mol and the relatively high energy 7.4 form (35.54 kcal/mol) is obtained. Although not feasible by the straightforward intramolecular process, the external www.nature.com/scientificreports/ ESI proton can also attack the N26 site generating the 7.5 structure of the relative energy 8.91 kcal/mol. Again, whereas some of the above-mentioned barrier are prohibitively, they can be easily crossed assuming the high collision energy in ESI HCD MS/MS. Thus, one can expect a numerous proton transfer processes occurring in the investigated system, weakening instantaneously different bonds in the molecule and allowing for the complex fragmentation processes. The careful analysis of the Wiberg bond indexes for 7 favors the destruction of the part of the structure containing BF 2 unit, since the O13-B29 (Wiberg bond index for the C10-protonated structure 0.5574), N7-B29 (0.6050) and B-F bonds (0.7799 and 0.7612) are the weakest in the whole system. The modification of the heterocyclic ring composition from imidazole in 7 to thiazole in 8 does affect much neither the structure of the neutral form of the system nor the sequence of the protonation affinities, still for C10 protonation affinity being the largest (229.49 kcal/mol), followed by the N(CH 3 ) 2 PA (222.55 kcal/mol), fluorine site protonation with the HF detachment (211.29 kcal/mol), O13 protonation (PA equal to 207.13 kcal/ mol), N7 protonation (PA 206.30 kcal/mol) and finally-the least readily protonated is the S9 site (184.21 kcal/ mol). This PA's ordering agrees with the energetic sequence of the protonated forms (compare Supplementary  Information). The protonation affinity order for 9 directly follows the sequence observed both for 7 and 8. The detailed discussion of the protonation affinity at different sites together with the Wiberg bond indexes analysis and energy barriers for protion transfer for 7-9 is provided in SI.  (Fig. 7). It should be emphasized that, although the process of elimination of a methyl radical was observed in the case of gas-phase fragmentation of various BODIPYs molecules 17 and dimethylamino derivatives of 2-phenacylheterocycles, it was not a privileged/dominant process, which is consistent with the even electron rule. In general, exceptions to this rule are rare but have been reported for conjugated structures and compounds containing aromatic rings 70 . The formation of [X-CH · 3 +H] +· products (where X=7, 8,9), regardless of the structure of the five-membered heterocyclic ring of the studied molecules, indicates that the methyl radicals originate from the dimethylamino groups. Furthermore, because the studied difluoroborates 7-9 differ from the corresponding 4-6 compounds only by the presence of the BF 2 group, the difference in HCD fragmentation behavior must be related to the properties of boron and/or fluorine atoms. The thesis can be supported by the results of a theoretical study of the influence of the boron atom and -BF 2 group on the acidity of alcohols (i.e., CH 3 OH , BH 2 CH 2 OH , BF 2 CH 2 OH ) in the gas phase, performed by Kheirjou et al. 71 . They reported that alcohol boronation considerably enhances its acidity, which can be related to the stabilization of the alkoxy ion due to the overlap of the unoccupied orbital of the boron atom with the electron pairs of the negative oxygen. The results of Natural Bond Orbital (NBO) analysis have shown that an increase in the acidity of the boron-containing alcohols is due to the charge transfer from the negative oxygen (in deprotonated structure) to the empty orbital of -BH 2 and -BF 2 71 . It could be assumed that in the case of the analyzed difluoroborates a similar phenomenon occurs. However, the results of DFT calculations performed for all compounds containing the BF 2 group (7-9) have shown, that in the most stable forms of ions the boron atoms are coordinated by heterocyclic nitrogen   www.nature.com/scientificreports/ atoms. Consequently, these nitrogen atoms become a less favorable sites for the ionizing protons compared to (-N=) atoms of 2-phenacylheterocycles. Fragmentation processes related to location of ionizing protons at (-N=) atoms of difluoroborates, leading to formation of [ C 9 H 10 NO] + ions, occurs, but, are not the dominant dissociation reactions, as was the case with the dimethylamino derivatives of 2-phenacylheterocycles ([C 9 H 10 NO] + -major products in HCD of singly charged ions of 4-6, analogous to dominant [ C 7 H 5 O] + ions formed in HCD of protonated 1-3 molecules). Instead, the dominant dissociation processes of compounds 7-9 are associated with the protonation (ESI) of the nitrogen atoms of the dimethylamino groups and lead to the privileged loss of methyl radicals. This process is rationalized by the large values of the corresponding proton affinity and gas-phase basicity at the N26 centre (Tab. S17) and the Wiberg bond indexes for N26-C32 and N26-C36 bonds lower than 1.0 (Tab. S18).
More details on the secondary fragmentation processes of difluoroborates including the elimination of the HF, HOBF 2 or CH 3 SBF 2 moieties, is provided in SI.

Summary and conclusions
The knowledge of the stability of fluorophores in the light of their structure is important not only in their excited state where all photoreactions may proceed, but also in the ground state. The most probable reaction in ground state of BF 2 -carrying molecules is the loss of the BF 2 moiety leading to their precursors. Since those precursors exhibit tautomerism, their properties may, again, induce side reactions. Their knowledge and control can be vital for further utilization, and therefore a detailed information on the structural features of these systems with the emphasis on relative energies and stability of bonds is crucial. Based on those, one can even more precisely design compounds for special applications. In the light of the above, the structural changes introduced in the present study (H vs. NMe 2 and N vs. O vs. S) show how important are subtle modification for fragmentation.
All the theoretical calculations consistently support the mechanism of the fragmentation of 2-phenacylheterocycles as well as their difluoroborate counterparts, suggested on the basis of the experimental MS measurements. Wiberg bond indexes correlate well with the AIM electron density and together with the energy for the investigated reaction paths they rationalize the detachment of the benzoyl ion.
The structural elucidation of a set of 2-phenacylheterocycles which differ by the elements located in heterocyclic moieties ( NCH 3 /S/O), and the substituents attached to the benzene rings (presence/absence of -N(CH 3 ) 2 ), as well as for their difluoroboranyl derivatives has been performed by means of ESI HCD MS/MS technique. For all of the analyzed ions, gas-phase charge-directed fragmentation processes strongly depend on the number and properties of heteroatoms. The formation of dominant HCD products of all 2-phenacylheterocycles with and without a dimethylamino groups is related to the location of ionizing proton(s) at the heterocyclic ring nitrogen atoms (-N=) (or at the carbon atoms from -CH-moieties, followed by intramolecular proton transfers). The presence of additional -N(CH 3 ) 2 group bonded to the benzene ring causes the formation of minor fragmentation products associated with elimination of the entire or part of the dimethylamino substituent (i.e. 3 +H] +· ). It also has an effect on other, minor fragmentation processes (i.e. reduction of intensity of the signals corresponding to the ions formed as a result of protonation of heteroatoms located in the heterocyclic ring). In the case of difluoroboranyl derivatives of 2-phenacylheterocycles, HCD leads to the formation of dominant [X-CH · 3 +H] +· products, related to the protonation of nitrogen of -N(CH 3 ) 2 . Additionally, the results confirm that the loss of HOBF 2 , which has been observed in ESI MS/MS experiments performed for all difluoroboranyl derivatives (presented here and reported previously 28 ), can serve as a signature tag for the presence of the -O-BF 2 moiety. As the elimination of neutral HSBF 2 was observed for both, the compounds containing -SBF 2 28 group and those with -OBF 2 unit and sulfur atom located in the heterocyclic ring, this fragmentation process does not provide information about the location of the difluoroboranyl moiety and it does not play such analytical role as HOBF 2 elimination.
The NCH 3 /S/O replacement in the five-membered ring has a strong impact on minor HCD fragmentation reactions. For example, in the case of compounds containing a sulfur atom, the loss of H 2 S is observed for those with a heterocyclic oxygen atom, the elimination of CO occurs, and the presence of a N(CH 3 ) group in the ring leads to detachment of the methyl radical.
The results clearly show that ESI HCD MS/MS provides detailed information on all structural motifs of the analyzed 2-phenacylheterocycles and their difluoroboranyl derivatives and can be used as a fast and simple tool for the structural analysis of novel difluoroborates. Furthermore, samples for MS and MS/MS analysis do not require any special, time-consuming preparation, a small amount (ppm) of compounds is required for the experiments, and HCD provides simple mass spectra that are easy to interpret. However, the elucidation of the details of the HCD fragmentation mechanisms of the analyzed compounds requires the involvement of theoretical chemistry methods, such as DFT, which provides valuable complementary information.
It should be underlined that the detailed theoretical study of the secondary fragmentation processes is not provided here. One need also note that the energetic, enthalpic and Gibbs free energy data should be interpreted with caution, taking into account the particular experimental HCD conditions which do not allow for reaching of the thermal equilibrium.