# Characterization of new psychoactive substances by integrating benchtop NMR to multi-technique databases

Thomas Castaing-Cordier, Alejandra Benavides Restrepo, Damien Dubois, Virginie Ladroue, Fabrice Besacier, Audrey Buleté, Céline Charvoz, Anais Goupille, Denis Jacquemin, Patrick Giraudeau, et al.

► To cite this version:

Thomas Castaing-Cordier, Alejandra Benavides Restrepo, Damien Dubois, Virginie Ladroue, Fabrice Besacier, et al. Characterization of new psychoactive substances by integrating benchtop NMR to multi-technique databases. Drug Testing and Analysis, 2022, 14 (9), pp.1629-1638. 10.1002/dta.3332. hal-03800479

HAL Id: hal-03800479

https://hal.science/hal-03800479

Submitted on 31 May 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Characterization of new psychoactive substances by integrating benchtop NMR to multi-technique databases

Thomas Castaing-Cordier, Alejandra Benavides Restrepo, Damien Dubois, Virginie Ladroue, Fabrice Besacier, Audrey Buleté, Céline Charvoz, Anais Goupille, Denis Jacquemin, Patrick Giraudeau, et al.

► To cite this version:

Thomas Castaing-Cordier, Alejandra Benavides Restrepo, Damien Dubois, Virginie Ladroue, Fabrice Besacier, et al. Characterization of new psychoactive substances by integrating benchtop NMR to multi-technique databases. Drug Testing and Analysis, 2022, 14 (9), pp.1629-1638. 10.1002/dta.3332. hal-03800479

HAL Id: hal-03800479

https://hal.science/hal-03800479v1

Submitted on 31 May 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
from different families.²,³ These NPS include synthetic cannabinoids derived from marijuana, ketamine-like compounds which are used as anesthetics, fentanyloids which is a painkiller more powerful than morphine, and benzodiazepines that are used for curing anxiety, cathinones, and phenethylamines. The effects of those substances generally remain unknown and their consumption can lead to intoxication, overdose, or even death. Due to the growing consumption of NPS on the market as well as the difficulty to constantly adapt the regulation of European members, there is an urgent need for accessible and reliable analytical tools allowing improved tracking of NPS seized by the Police and/or purchased on the internet. Indeed, it is essential to accurately characterize NPS, to recognize molecules that were previously identified, and to elucidate the structure of unknown ones. Identification of known substances consists in comparing experimental data of a seized compound with reference data, which is quite straightforward. In contrast, elucidation is a more challenging task, necessary when identification fails and consisting in determining the structure of the unknowns in a sample from a set of combined analytical data. An efficient NMR-based approach named CASE (computer-assisted spectral elucidation) has been developed for structure elucidation of small molecules.⁴ CASE uses software to generate all possible molecular structures that are consistent with a particular set of 1D and 2D NMR data. However, CASE relies on the determination of the exact mass obtained by high-resolution mass spectrometry (HRMS), a technique which is not widely available in forensic science services. This motivated us to develop a workflow in which the elucidation is assisted with databases and does not require prior knowledge of the exact mass.

NPS characterization can be achieved by a variety of complementary analytical techniques. The first one is gas chromatography coupled to mass spectrometry (GC-MS), currently the most used in forensics.⁵ MS combined with chemical ionization (CI) may be used to provide the molecular weight, but a relatively expensive HRMS instrumentation is required to access accurate molecular weight. MS combined with electron ionization (EI) mode provides information on different fragments of the molecule, and hyphenated MS/MS (or MS¹) delivers further insights into the fragmentation schemes of the molecule. Another important technique is infrared (IR) spectroscopy which is mainly used to determine chemical functions within a molecule, but can also be applied for spectral recognition.⁶⁻⁸ Unfortunately, these methods provide limited structural information. NMR spectroscopy is probably the most powerful structure identification and elucidation tool for small molecules such as NPS. Indeed, it offers highly accurate and specific information on the chemical environment of all atoms through chemical shifts, while providing crucial input on atomic connectivities through J-couplings. Moreover, highly accurate quantification is achievable with NMR with limited sample preparation and a short acquisition time for 1D experiments. However, high-field NMR (¹H frequency > 300 MHz) is rarely used in forensics due to the high purchase, maintenance, and running costs in addition to the need for dedicated staff and to the bulkiness of the NMR instrument. However, during the last decade, new compact NMR apparatus emerged, with magnetic fields between 1 and 2.1 T. These instruments can be

installed on a benchtop, are easily transportable, and do not require maintenance (no cryogenic fluids). Nevertheless, benchtop NMR remains a young and underexplored technique for NPS characterization as compared with conventional high-field NMR.⁹,¹⁰ Also, benchtop NMR has to face intrinsic sensitivity and resolution limits, which are illustrated in Figure 1 in the case of two common NPS.

Such limitations make it difficult to only rely on benchtop NMR for structure identification or elucidation. Therefore, combinations with other methods like GC-MS or IR and/or with databases become a prerequisite for a reliable, fast, and automated characterization. Two

<figure><img src="image_2.png" /><figcaption>**FIGURE 1** Spectrum of (a) 370mM AMB-FUBINACA in DMSO-d6 and (b) 306mM 2-BMMP in DMSO-d6 acquired at high-field (top: 400 MHz) and on a benchtop spectrometer (bottom: 60 MHz). Peak assignment after phasing and baseline correction is indicated. The two spectra are obtained with the same conditions: eight scans, a repetition time (TR) of 30 s, and an acquisition time of 1.6 s at 299.6 K. Most of the signals are overlapped at 60 MHz especially for the aromatic area of AMB-FUBINACA. This illustrates the resolution limitation of benchtop NMR. Indeed, since the resolution scales linearly with the magnetic field B0 in theory, the 400 MHz spectrum is expected to be roughly seven times more resolved than at 60 MHz. As regards sensitivity, it scales with B0³/²; hence, the 400 MHz spectrum is expected to be approximately 17 times more sensitive than a 60 MHz spectrum recorded in identical conditions [Colour figure can be viewed at wileyonlinelibrary.com]</figcaption></figure>
from different families.²,³ These NPS include synthetic cannabinoids derived from marijuana, ketamine-like compounds which are used as anesthetics, fentanyloids which is a painkiller more powerful than morphine, and benzodiazepines that are used for curing anxiety, cathinones, and phenethylamines. The effects of those substances generally remain unknown and their consumption can lead to intoxication, overdose, or even death. Due to the growing consumption of NPS on the market as well as the difficulty to constantly adapt the regulation of European members, there is an urgent need for accessible and reliable analytical tools allowing improved tracking of NPS seized by the Police and/or purchased on the internet. Indeed, it is essential to accurately characterize NPS, to recognize molecules that were previously identified, and to elucidate the structure of unknown ones. Identification of known substances consists in comparing experimental data of a seized compound with reference data, which is quite straightforward. In contrast, elucidation is a more challenging task, necessary when identification fails and consisting in determining the structure of the unknowns in a sample from a set of combined analytical data. An efficient NMR-based approach named CASE (computer-assisted spectral elucidation) has been developed for structure elucidation of small molecules.⁴ CASE uses software to generate all possible molecular structures that are consistent with a particular set of 1D and 2D NMR data. However, CASE relies on the determination of the exact mass obtained by high-resolution mass spectrometry (HRMS), a technique which is not widely available in forensic science services. This motivated us to develop a workflow in which the elucidation is assisted with databases and does not require prior knowledge of the exact mass.

NPS characterization can be achieved by a variety of complementary analytical techniques. The first one is gas chromatography coupled to mass spectrometry (GC-MS), currently the most used in forensics.⁵ MS combined with chemical ionization (CI) may be used to provide the molecular weight, but a relatively expensive HRMS instrumentation is required to access accurate molecular weight. MS combined with electron ionization (EI) mode provides information on different fragments of the molecule, and hyphenated MS/MS (or MS¹) delivers further insights into the fragmentation schemes of the molecule. Another important technique is infrared (IR) spectroscopy which is mainly used to determine chemical functions within a molecule, but can also be applied for spectral recognition.⁶⁻⁸ Unfortunately, these methods provide limited structural information. NMR spectroscopy is probably the most powerful structure identification and elucidation tool for small molecules such as NPS. Indeed, it offers highly accurate and specific information on the chemical environment of all atoms through chemical shifts, while providing crucial input on atomic connectivities through J-couplings. Moreover, highly accurate quantification is achievable with NMR with limited sample preparation and a short acquisition time for 1D experiments. However, high-field NMR (¹H frequency > 300 MHz) is rarely used in forensics due to the high purchase, maintenance, and running costs in addition to the need for dedicated staff and to the bulkiness of the NMR instrument. However, during the last decade, new compact NMR apparatus emerged, with magnetic fields between 1 and 2.1 T. These instruments can be

installed on a benchtop, are easily transportable, and do not require maintenance (no cryogenic fluids). Nevertheless, benchtop NMR remains a young and underexplored technique for NPS characterization as compared with conventional high-field NMR.⁹,¹⁰ Also, benchtop NMR has to face intrinsic sensitivity and resolution limits, which are illustrated in Figure 1 in the case of two common NPS.

Such limitations make it difficult to only rely on benchtop NMR for structure identification or elucidation. Therefore, combinations with other methods like GC-MS or IR and/or with databases become a prerequisite for a reliable, fast, and automated characterization. Two

<figure><img src="image_2.png" /><figcaption>**FIGURE 1** Spectrum of (a) 370mM AMB-FUBINACA in DMSO-d6 and (b) 306mM 2-BMMP in DMSO-d6 acquired at high-field (top: 400 MHz) and on a benchtop spectrometer (bottom: 60 MHz). Peak assignment after phasing and baseline correction is indicated. The two spectra are obtained with the same conditions: eight scans, a repetition time (TR) of 30 s, and an acquisition time of 1.6 s at 299.6 K. Most of the signals are overlapped at 60 MHz especially for the aromatic area of AMB-FUBINACA. This illustrates the resolution limitation of benchtop NMR. Indeed, since the resolution scales linearly with the magnetic field B0 in theory, the 400 MHz spectrum is expected to be roughly seven times more resolved than at 60 MHz. As regards sensitivity, it scales with B0³/²; hence, the 400 MHz spectrum is expected to be approximately 17 times more sensitive than a 60 MHz spectrum recorded in identical conditions [Colour figure can be viewed at wileyonlinelibrary.com]</figcaption></figure>
such multi-method approaches were already explored in a forensic
framework¹¹ In a first study merging NMR and GC-MS, ¹H NMR
spectra of reference compounds were collected using an 80 MHz
instrument to create a reference library of 302 spectra of different
NPS classes.¹² Next 432 seized samples were analyzed by NMR and
GC-MS for cross-validation. ¹H NMR analysis nicely matched the
GC-MS results with a 93% consistency rate. Another study reported
the identification of drugs in two case samples, using a limited library
of 12 spectra on an 80 MHz benchtop spectrometer, visually com-
pared with data collected on a 600 MHz spectrometer.¹³ In agree-
ment with GC-MS analysis, the seized samples were found to contain
morphine, acetyl codeine, MAM, and MDMA. In spite of these
encouraging preliminary results, there is no general workflow incorpo-
rating benchtop NMR for the identification and elucidation of NPS
structures. A first reason is that the above-mentioned examples were
limited to 1D ¹H NMR spectroscopy and did not include 2D NMR
experiments, which should provide a much higher degree of confi-
dence in terms of structure characterization. Moreover, the reported
elucidation procedures rely on the determination of the exact mass as
an initial step, but as explained above, HRMS instruments are uncommon in forensic laboratories. It is also important to note that the
SWDRUGS group recommends using at least two techniques for
unambiguous identification, which motivated us to combine benchtop
NMR with another analytical approach.¹⁴

Databases can be used with other analytical techniques such as IR. For instance, in a study by Jones et al., 221 samples—most of them in mixtures—were screened using IR and Raman.¹⁵ Authors compared the spectra to the database, and if a reference matched, they subtracted the reference spectra to the unknown one and made another comparison with the database. Only 41% of samples were unambiguously identified. Other samples were sent to MS and NMR, for a full elucidation procedure. With this approach, 33 samples were identified and added to the IR and Raman library.

About NMR quantification, some previous works were carried
out on benchtop NMR spectra. Naqi et al. quantified seized MDMA
by using ¹H-qNMR, UHPLC, and UHPLC-MS.¹⁶ Maleic acid and
MDMA-d₅ were used, respectively, as an internal reference for NMR
and UHPLC measurements. The MDMA concentrations determined
by UHPLC and NMR were found comparable with no significant sta-
tistical difference as revealed by ANOVA single factor analysis. Quan-
titative NMR with benchtop NMR was applied to drugs with similar
molecular structures. For example, Hussain et al. were able to quantify
the mass weight of MDMA contained in a tablet.¹⁷ The values
obtained on eight tablets were between 209 and 212 mg and are
comparable with results given by GC-MS.

In this context, this work aims to develop and evaluate a general
 workflow based on optimized 1D and 2D NMR measurements per-
formed on a benchtop spectrometer, combined with IR spectroscopy,
for creating a multi-technique database allowing the characterization
of pure and mixed NPS. Both identification and elucidation are consid-
ered in this workflow, which uses accessible analytical techniques
while keeping in mind the demand of end-users in terms of reliability,
robustness, and experiment time. After describing its implementation

and the choice of the key parameters, the workflow was blind-tested
on six real-seized samples, and the quantification capabilities of
benchtop NMR were also assessed.

2 | MATERIALS AND METHODS

2.1 | Sample preparation

Fifty-seven samples were used in this study. Five of them were pur-
chased from Lipomed, seven samples were seized by the Finnish
police, while all the others came from seizures by the French police.
Also, six unknown seizures were obtained from the French police to
allow validating the database. DMSO-d₆ purchased from Eurisotop
with a purity of 99.8% was used as an NMR solvent. TMS from Acros
Organics with a purity of 99.9% and TSP from Eurisotop with a purity
of 99.8% were used as chemical shift and concentration references,
respectively. All seizures were dissolved in DMSO-d₆ to obtain a con-
centration as close as possible to 300 mM (see Table S1), and 10 µl of
TMS was added to calibrate ¹H and ¹³C chemical shifts.

## 2.2 | NMR

All the spectra were recorded at 26.5°C using a ¹H, ¹⁹F, ¹³C 60 MHz
Spinsolve spectrometer from Magritek equipped with a 20-position
autosampler and a gradient coil in the direction of the static magnetic
field B₀. The duration of 90° hard pulses was 14 µs at 0 dBW for ¹H,
60 µs at 0 dBW for ¹⁹F, and 65.9 µs at -14 dBW for ¹³C. All acquisi-
tions were made with the SpinsolveExpert software (v.1.41.07). More-
over, an Avance III HD 700 MHz spectrometer with a ¹H/¹³C/¹⁵N
cryoprobe from Bruker was used to confirm the identification/
elucidation of the unknown NPS in samples.

Three NMR pulse sequences were chosen for our study, 1D ¹H, 1D ¹⁹F, and 2D ¹H-¹³C heteronuclear single quantum coherence (HSQC). The HSQC sequence was improved by optimization with the help of a test molecule, named 2-FDCK (see Table S1). Three variants of the HSQC sequence were tested, namely, Echo-Antiecho (EA) HSQC with WALTZ decoupling, phase sensitive (PH) HSQC with WALTZ decoupling, PH HSQC with Multiplicity Edition (ME), and WALTZ decoupling. PH HSQC was found in average twice more sensitive than EA HSQC (see S2). The PH HSQC ME (see S3) was found to be the most informative since it allowed determining the multiplicity *n* of ¹³C(¹H)*n* but this pulse sequence was 20% less sensitive than PH HSQC. The most important acquisitions parameters are gathered in Table 1.

The NMR processing was performed with MestreNova™ software
(version 14.1). The 1D free induction decays were first multiplied by
an exponential apodization function (0.3 Hz line broadening) and
zero-filled to a factor 4 then manually phased and baseline corrected
with a Whittaker smoother algorithm. The 2D free induction decays
were first multiplied by a cosine apodization function and zero-filled
to a factor 4 in both dimensions, then manually phased and baseline
such multi-method approaches were already explored in a forensic framework¹¹ In a first study merging NMR and GC-MS, ¹H NMR spectra of reference compounds were collected using an 80 MHz instrument to create a reference library of 302 spectra of different NPS classes.¹² Next 432 seized samples were analyzed by NMR and GC-MS for cross-validation. ¹H NMR analysis nicely matched the GC-MS results with a 93% consistency rate. Another study reported the identification of drugs in two case samples, using a limited library of 12 spectra on an 80 MHz benchtop spectrometer, visually compared with data collected on a 600 MHz spectrometer.¹³ In agreement with GC-MS analysis, the seized samples were found to contain morphine, acetyl codeine, MAM, and MDMA. In spite of these encouraging preliminary results, there is no general workflow incorporating benchtop NMR for the identification and elucidation of NPS structures. A first reason is that the above-mentioned examples were limited to 1D ¹H NMR spectroscopy and did not include 2D NMR experiments, which should provide a much higher degree of confidence in terms of structure characterization. Moreover, the reported elucidation procedures rely on the determination of the exact mass as an initial step, but as explained above, HRMS instruments are uncommon in forensic laboratories. It is also important to note that the SWDRUGS group recommends using at least two techniques for unambiguous identification, which motivated us to combine benchtop NMR with another analytical approach.¹⁴

Databases can be used with other analytical techniques such as IR. For instance, in a study by Jones et al., 221 samples—most of them in mixtures—were screened using IR and Raman.¹⁵ Authors compared the spectra to the database, and if a reference matched, they subtracted the reference spectra to the unknown one and made another comparison with the database. Only 41% of samples were unambiguously identified. Other samples were sent to MS and NMR, for a full elucidation procedure. With this approach, 33 samples were identified and added to the IR and Raman library.

About NMR quantification, some previous works were carried out on benchtop NMR spectra. Naqi et al. quantified seized MDMA by using ¹H-qNMR, UHPLC, and UHPLC-MS.¹⁶ Maleic acid and MDMA-d₅ were used, respectively, as an internal reference for NMR and UHPLC measurements. The MDMA concentrations determined by UHPLC and NMR were found comparable with no significant statistical difference as revealed by ANOVA single factor analysis. Quantitative NMR with benchtop NMR was applied to drugs with similar molecular structures. For example, Hussain et al. were able to quantify the mass weight of MDMA contained in a tablet.¹⁷ The values obtained on eight tablets were between 209 and 212 mg and are comparable with results given by GC-MS.

In this context, this work aims to develop and evaluate a general workflow based on optimized 1D and 2D NMR measurements performed on a benchtop spectrometer, combined with IR spectroscopy, for creating a multi-technique database allowing the characterization of pure and mixed NPS. Both identification and elucidation are considered in this workflow, which uses accessible analytical techniques while keeping in mind the demand of end-users in terms of reliability, robustness, and experiment time. After describing its implementation

and the choice of the key parameters, the workflow was blind-tested on six real-seized samples, and the quantification capabilities of benchtop NMR were also assessed.

## 2 | MATERIALS AND METHODS

### 2.1 | Sample preparation

Fifty-seven samples were used in this study. Five of them were purchased from Lipomed, seven samples were seized by the Finnish police, while all the others came from seizures by the French police. Also, six unknown seizures were obtained from the French police to allow validating the database. DMSO-d₆ purchased from Eurisotop with a purity of 99.8% was used as an NMR solvent. TMS from Acros Organics with a purity of 99.9% and TSP from Eurisotop with a purity of 99.8% were used as chemical shift and concentration references, respectively. All seizures were dissolved in DMSO-d₆ to obtain a concentration as close as possible to 300 mM (see Table S1), and 10 µl of TMS was added to calibrate ¹H and ¹³C chemical shifts.

### 2.2 | NMR

All the spectra were recorded at 26.5°C using a ¹H, ¹⁹F, ¹³C 60 MHz Spinsolve spectrometer from Magritek equipped with a 20-position autosampler and a gradient coil in the direction of the static magnetic field B₀. The duration of 90° hard pulses was 14 µs at 0 dBW for ¹H, 60 µs at 0 dBW for ¹⁹F, and 65.9 µs at -14 dBW for ¹³C. All acquisitions were made with the SpinsolveExpert software (v.1.41.07). Moreover, an Avance III HD 700 MHz spectrometer with a ¹H/¹³C/¹⁵N cryoprobe from Bruker was used to confirm the identification/ elucidation of the unknown NPS in samples.

Three NMR pulse sequences were chosen for our study, 1D ¹H, 1D ¹⁹F, and 2D ¹H-¹³C heteronuclear single quantum coherence (HSQC). The HSQC sequence was improved by optimization with the help of a test molecule, named 2-FDCK (see Table S1). Three variants of the HSQC sequence were tested, namely, Echo-Antiecho (EA) HSQC with WALTZ decoupling, phase sensitive (PH) HSQC with WALTZ decoupling, PH HSQC with Multiplicity Edition (ME), and WALTZ decoupling. PH HSQC was found in average twice more sensitive than EA HSQC (see S2). The PH HSQC ME (see S3) was found to be the most informative since it allowed determining the multiplicity *n* of ¹³C(¹H)*n* but this pulse sequence was 20% less sensitive than PH HSQC. The most important acquisitions parameters are gathered in Table 1.

The NMR processing was performed with MestreNova™ software (version 14.1). The 1D free induction decays were first multiplied by an exponential apodization function (0.3 Hz line broadening) and zero-filled to a factor 4 then manually phased and baseline corrected with a Whittaker smoother algorithm. The 2D free induction decays were first multiplied by a cosine apodization function and zero-filled to a factor 4 in both dimensions, then manually phased and baseline
<table><caption>TABLE 1 NMR acquisitions parameters</caption><thead><tr><th>Parameters</th><th>1D <sup>1</sup>H</th><th>1D <sup>19</sup>F</th><th>2D <sup>1</sup>H-<sup>13</sup>C HSQC</th></tr></thead><tbody><tr><td>NS</td><td>8</td><td>360</td><td>24</td></tr><tr><td>DW (μs)</td><td>200</td><td>50</td><td>100</td></tr><tr><td>AQ (ms)</td><td>1638.4</td><td>204.8</td><td>409.6</td></tr><tr><td>TR (s)</td><td>30</td><td>1</td><td>2</td></tr><tr><td>NI</td><td>—</td><td>—</td><td>128</td></tr><tr><td>T<sub>exp</sub> (min)</td><td>4</td><td>6</td><td>103</td></tr></tbody></table>

Note: NS is the number of scans, DW is the dwell time, AQ is the Acquisition time, TR is the repetition time (including the acquisition time and the recovery delay), and NI is the number of increments in the indirect dimension of 2D HSQC spectra. T<sub>exp</sub> indicates the experiment time.

corrected with Whittaker smoother algorithm. <sup>1</sup>H and <sup>13</sup>C chemical shifts were referenced to TMS, and <sup>19</sup>F chemical shifts were referenced to the internal lock substance of the spectrometer.

## 2.3 | IR

Analyses were carried out with a Bruker tensor 27 spectrometer in Attenuated Total Reflectance (ATR) mode. Each spectrum was acquired with 32 scans with a 400–4000 cm<sup>-1</sup> spectral range and a 2 cm<sup>-1</sup> spectral resolution. The processing of the spectra was achieved with the Opus™ (v8.5) software. CO<sub>2</sub>, H<sub>2</sub>O compensation, baseline, and smoothing were applied to all IR spectra to obtain more reliable NPS IR fingerprints.

## 2.4 | ACD/Labs software

All analytical data were gathered in the ACD/Labs software (version 2020.1.0 of July 15, 2020) from Advanced Chemistry Development, Inc (ACD/Labs). First, a project with all NMR experiments and the dedicated structure was created for each molecule and added to the database after peak picking and assignment of all signals and correlations. Then the corresponding IR spectrum was added to the database after an automatic peak picking.

## 2.5 | Database parameters

Once the database has been created, two search methods are available in the ACD/Labs software to compare an unknown spectrum with one from the database. The first one, denoted “similarity search”, compares the shape of the spectra only, disregarding peak picking and the multiplicity of the signals. The second one, called “peak searching”, compares the peak picking of the NMR signal, an approach which is, obviously, significantly operator-dependent. With both search methods, the comparison between the experimental

spectrum and the database yields a result called the Hit Quality Index (HQI) which quantifies the agreement between the database spectrums the most closely matching and the spectrum of the unknown compound.

In the “similarity search” method, HQI can be calculated with two algorithms detailed below. In both cases, the query spectrum is first indexed and split into N small regions. Next, all regions are integrated and, depending on the calculated integral, an index in the [-127,+127] range is assigned to each region.

The first algorithm relies on the calculation of the Euclidian difference and treats each spectrum as a vector in a N-dimensional space. Thus, the resulting HQI is a measurement of the angle between the vectors giving a value between 0 (best match) and the square root of 2 (worst match). In the ACD/Labs software, this value is next scaled to deliver a value between 0 (worst match) and 100 (best match) for the HQI display:

$$ \text{HQI} = \frac{1}{N \times 127} \sqrt{\sum_{i=1}^{L} \text{abs}(p_i^{\text{EXP}} - p_i^{\text{DB}})^2}, \quad (1) $$

with L the compared regions, N the number of used indexes, and p<sub>i</sub> the index from the experimental (EXP) spectrum and the database spectrum (DB).

The second algorithm relies on the calculation of the absolute distance based on the difference for each data point with respect to the total area under the curve (i.e., the sum of all the intensities). In other words, each data point represents a percentage of the total area. Thus, the larger the intensity value of a data point, the larger its weight:

$$ \text{HQI} = \frac{1}{N \times 127} \sum_{i=1}^{L} | (p_i^{\text{EXP}} - p_i^{\text{DB}}) |. \quad (2) $$

Only one HQI algorithm is available for “peak searching” calculated as follows:

$$ \text{HQI} = \frac{1}{L} \sum_{k=1}^{L} \frac{N_k}{M_k} \times 100\%, \quad (3) $$

with L the total number of peak regions in the experimental spectrum, N<sub>k</sub> the number of peaks in the k-peak region of the experimental spectrum, and M<sub>k</sub> the number of peaks in the k-peak region of the DB spectrum.

Finally, for improving the database search, two parameters can be chosen by the operator. The first one is the HQI threshold value under which the results are filtered out. The second one, only available for “peak searching”, is the looseness factor (LF) that allows considering the variability in the chemical shifts of the different peaks obtained in NMR. LF corresponds to the highest difference in ppm which can be accepted to provide a match between the experimental spectra and the database (both in 1D and 2D).
<table><caption>TABLE 1 NMR acquisitions parameters</caption><thead><tr><th>Parameters</th><th>1D <sup>1</sup>H</th><th>1D <sup>19</sup>F</th><th>2D <sup>1</sup>H-<sup>13</sup>C HSQC</th></tr></thead><tbody><tr><td>NS</td><td>8</td><td>360</td><td>24</td></tr><tr><td>DW (μs)</td><td>200</td><td>50</td><td>100</td></tr><tr><td>AQ (ms)</td><td>1638.4</td><td>204.8</td><td>409.6</td></tr><tr><td>TR (s)</td><td>30</td><td>1</td><td>2</td></tr><tr><td>NI</td><td>—</td><td>—</td><td>128</td></tr><tr><td>T<sub>exp</sub> (min)</td><td>4</td><td>6</td><td>103</td></tr></tbody></table>

Note: NS is the number of scans, DW is the dwell time, AQ is the Acquisition time, TR is the repetition time (including the acquisition time and the recovery delay), and NI is the number of increments in the indirect dimension of 2D HSQC spectra. T<sub>exp</sub> indicates the experiment time.

corrected with Whittaker smoother algorithm. <sup>1</sup>H and <sup>13</sup>C chemical shifts were referenced to TMS, and <sup>19</sup>F chemical shifts were referenced to the internal lock substance of the spectrometer.

## 2.3 | IR

Analyses were carried out with a Bruker tensor 27 spectrometer in Attenuated Total Reflectance (ATR) mode. Each spectrum was acquired with 32 scans with a 400–4000 cm<sup>-1</sup> spectral range and a 2 cm<sup>-1</sup> spectral resolution. The processing of the spectra was achieved with the Opus™ (v8.5) software. CO<sub>2</sub>, H<sub>2</sub>O compensation, baseline, and smoothing were applied to all IR spectra to obtain more reliable NPS IR fingerprints.

## 2.4 | ACD/Labs software

All analytical data were gathered in the ACD/Labs software (version 2020.1.0 of July 15, 2020) from Advanced Chemistry Development, Inc (ACD/Labs). First, a project with all NMR experiments and the dedicated structure was created for each molecule and added to the database after peak picking and assignment of all signals and correlations. Then the corresponding IR spectrum was added to the database after an automatic peak picking.

## 2.5 | Database parameters

Once the database has been created, two search methods are available in the ACD/Labs software to compare an unknown spectrum with one from the database. The first one, denoted “similarity search”, compares the shape of the spectra only, disregarding peak picking and the multiplicity of the signals. The second one, called “peak searching”, compares the peak picking of the NMR signal, an approach which is, obviously, significantly operator-dependent. With both search methods, the comparison between the experimental

spectrum and the database yields a result called the Hit Quality Index (HQI) which quantifies the agreement between the database spectrums the most closely matching and the spectrum of the unknown compound.

In the “similarity search” method, HQI can be calculated with two algorithms detailed below. In both cases, the query spectrum is first indexed and split into N small regions. Next, all regions are integrated and, depending on the calculated integral, an index in the [-127;+127] range is assigned to each region.

The first algorithm relies on the calculation of the Euclidian difference and treats each spectrum as a vector in a N-dimensional space. Thus, the resulting HQI is a measurement of the angle between the vectors giving a value between 0 (best match) and the square root of 2 (worst match). In the ACD/Labs software, this value is next scaled to deliver a value between 0 (worst match) and 100 (best match) for the HQI display:

$$ \text{HQI} = \frac{1}{N \times 127} \sqrt{\sum_{i=1}^{L} \text{abs}(p_i^{\text{EXP}} - p_i^{\text{DB}})^2}, \quad (1) $$

with L the compared regions, N the number of used indexes, and p<sub>i</sub> the index from the experimental (EXP) spectrum and the database spectrum (DB).

The second algorithm relies on the calculation of the absolute distance based on the difference for each data point with respect to the total area under the curve (i.e., the sum of all the intensities). In other words, each data point represents a percentage of the total area. Thus, the larger the intensity value of a data point, the larger its weight:

$$ \text{HQI} = \frac{1}{N \times 127} \sum_{i=1}^{L} | (p_i^{\text{EXP}} - p_i^{\text{DB}}) |. \quad (2) $$

Only one HQI algorithm is available for “peak searching” calculated as follows:

$$ \text{HQI} = \frac{1}{L} \sum_{k=1}^{L} \frac{N_k}{M_k} \times 100\%, \quad (3) $$

with L the total number of peak regions in the experimental spectrum, N<sub>k</sub> the number of peaks in the k-peak region of the experimental spectrum, and M<sub>k</sub> the number of peaks in the k-peak region of the DB spectrum.

Finally, for improving the database search, two parameters can be chosen by the operator. The first one is the HQI threshold value under which the results are filtered out. The second one, only available for “peak searching”, is the looseness factor (LF) that allows considering the variability in the chemical shifts of the different peaks obtained in NMR. LF corresponds to the highest difference in ppm which can be accepted to provide a match between the experimental spectra and the database (both in 1D and 2D).
## 2.6 | Analytical time

The overall data acquisition and analysis times should be kept as short as possible, to remain compatible with the typical duration of a police custody (48 to 72 h in France, for instance). In this study, optimized NMR experiments last a total time of 113 min in total for each sample (see Table 1) and IR experiments take less than 1 min per acquisition. Finally, the data processing steps following acquisition are not very time-consuming (a few minutes per experiment at most). Thus, for identification, the largest share of the processing time originates in the acquisition of analytical data. In contrast, for the elucidation of unknown NPS, the analysis of the spectral data by a trained operator is likely the rate-limiting step, which is also dependent on the complexity of the problem to solve.

## 2.7 | SWDRUGS recommendations

The choice of creating a database containing IR and NMR analytical data is perfectly justified by the recommendations of the SWDRUGS group.¹⁴ All analytical techniques can be separated into three different categories. Category A techniques provide the highest level of selectivity through structural information, category B techniques deliver an intermediate level of selectivity through physical or chemical characteristics but lack structural information, and category C techniques give rise to a low level of selectivity but provide general/chemical class informations. Unambiguous structure identification requires combining one category A technique to another from any category. According to SWDRUG, NMR (category

A) and IR (category A) convey enough information to validate the proposed structure.

# 3 | RESULTS AND DISCUSSIONS

## 3.1 | NMR methods

Among the broad variety of NMR experiments employed for structural analysis, a reduced set of three experiments was chosen for our identification and elucidation workflow: 1D ¹H, 1D ¹⁹F, and 2D ¹H-¹³C HSQC experiments. A systematic order was set up for the acquisition of these NMR experiments. First, 1D ¹⁹F spectroscopy was used because ¹⁹F nuclei were present in a minority of NPS.¹⁸,¹⁹ This experiment allows separating unknown molecules into two categories, those with ¹⁹F and those without (Figure 2). Second, ¹H-¹³C HSQC was used as it is one of the most sensitive and informative heteronuclear 2D pulse sequences. 2D HSQC spectra were acquired with spectral edition which allows to label the ¹³C-¹H multiplicity for each correlation. CH and CH₃ are represented in red (positive) and CH₂ in blue (negative). This second experiment allows creating a double entry table with ¹³C chemical shifts in rows and ¹H chemical shifts in columns encoded with C-Hₐ parity. Therefore, provided that the LF parameter is correctly adjusted, the spectrum of the unknown NPS should either yield a unique answer or deliver no match (i.e., HQI 100 or 0). In some cases, however, the HSQC spectra recorded on a benchtop spectrometer could be limited by sensitivity due to the low natural abundance of ¹³C nuclei. This is why 1D ¹H NMR spectra were also incorporated in the workflow, offering a more sensitive, yet

<figure><img src="image_3.png" /><figcaption>**FIGURE 2** Workflow for structural identification and elucidation. For the heteronuclear single quantum coherence (HSQC) comparison in the elucidation section, the looseness factor (LF) is incremented by x (a nonzero natural number) [Colour figure can be viewed at wileyonlinelibrary.com]</figcaption></figure>
## 2.6 | Analytical time

The overall data acquisition and analysis times should be kept as short as possible, to remain compatible with the typical duration of a police custody (48 to 72 h in France, for instance). In this study, optimized NMR experiments last a total time of 113 min in total for each sample (see Table 1) and IR experiments take less than 1 min per acquisition. Finally, the data processing steps following acquisition are not very time-consuming (a few minutes per experiment at most). Thus, for identification, the largest share of the processing time originates in the acquisition of analytical data. In contrast, for the elucidation of unknown NPS, the analysis of the spectral data by a trained operator is likely the rate-limiting step, which is also dependent on the complexity of the problem to solve.

## 2.7 | SWDRUGS recommendations

The choice of creating a database containing IR and NMR analytical data is perfectly justified by the recommendations of the SWDRUGS group.¹⁴ All analytical techniques can be separated into three different categories. Category A techniques provide the highest level of selectivity through structural information, category B techniques deliver an intermediate level of selectivity through physical or chemical characteristics but lack structural information, and category C techniques give rise to a low level of selectivity but provide general/chemical class informations. Unambiguous structure identification requires combining one category A technique to another from any category. According to SWDRUG, NMR (category

A) and IR (category A) convey enough information to validate the proposed structure.

# 3 | RESULTS AND DISCUSSIONS

## 3.1 | NMR methods

Among the broad variety of NMR experiments employed for structural analysis, a reduced set of three experiments was chosen for our identification and elucidation workflow: 1D ¹H, 1D ¹⁹F, and 2D ¹H-¹³C HSQC experiments. A systematic order was set up for the acquisition of these NMR experiments. First, 1D ¹⁹F spectroscopy was used because ¹⁹F nuclei were present in a minority of NPS.¹⁸,¹⁹ This experiment allows separating unknown molecules into two categories, those with ¹⁹F and those without (Figure 2). Second, ¹H-¹³C HSQC was used as it is one of the most sensitive and informative heteronuclear 2D pulse sequences. 2D HSQC spectra were acquired with spectral edition which allows to label the ¹³C-¹H multiplicity for each correlation. CH and CH₃ are represented in red (positive) and CH₂ in blue (negative). This second experiment allows creating a double entry table with ¹³C chemical shifts in rows and ¹H chemical shifts in columns encoded with C-Hₐ parity. Therefore, provided that the LF parameter is correctly adjusted, the spectrum of the unknown NPS should either yield a unique answer or deliver no match (i.e., HQI 100 or 0). In some cases, however, the HSQC spectra recorded on a benchtop spectrometer could be limited by sensitivity due to the low natural abundance of ¹³C nuclei. This is why 1D ¹H NMR spectra were also incorporated in the workflow, offering a more sensitive, yet

<figure><img src="image_3.png" /><figcaption>**FIGURE 2** Workflow for structural identification and elucidation. For the heteronuclear single quantum coherence (HSQC) comparison in the elucidation section, the looseness factor (LF) is incremented by x (a nonzero natural number) [Colour figure can be viewed at wileyonlinelibrary.com]</figcaption></figure>
less resolved technique to compare the spectral patterns of the unknown molecule with those included in the database.

## 3.2 | Multi-technique workflow

The choice of the analytical workflow is essential regarding the complexity of the sample embedding NPS. The strategy to build for identifying and elucidating NPS in seized samples has to face difficulties met to detect and quantify NPS even in mixture.

For a more comprehensive and integrated drug investigation, different multi-methods workflows have been reported and allow for elucidating NPS¹¹. They all combine multi-analytical techniques like Raman, DSC, IR, NMR, MS, and coupled chromatographic techniques with chemoinformatic tools and commercial, in-house databases for assisting spectral assignment, identification, and elucidation.

An applicative case where seized blotter paper²⁰ was analyzed by GC-MS, UHPLC-HRMS as well as 1D and 2D NMR in combination with multi-methods chemometric tools like ACD/Labs, Mass Hunter, and XCalibur. This workflow led to the identification of a LSD homologous product. In the case of seized NPS in mixture with herbs, the same workflow was utilized after NPS extraction and successfully revealed a new synthetic cannabinoid.¹⁸

Since high-field NMR and HRMS are expensive instruments and not owned by a majority of forensic laboratories, there is a need for cheaper and more accessible alternative workflows. We propose 1D and 2D optimized experiments implemented on liquid-state benchtop instruments in combination with IR and dedicated associated databases to identify and also elucidate NPS. Note that our approach aims at being quite general and may not be adapted to NPS hidden or sprayed in materials that would require a more specific strategy.

For our approach, an experimental database was created with 57 entries, mainly coming from cathinones, cannabinoids, amphetamines, arylcyclohexylamines, and fentanyloids (see Table S1). Before elaborating the multi-technique workflow, three model molecules were used to determine the best search method, the optimum algorithm for HQI calculation, the most adapted HQI threshold, as well as the optimum LF value. Table 2 summarizes the optimized parameters,

<table><caption>TABLE 2 Selected parameters for spectral comparison with ACD/Labs software</caption><thead><tr><th>Experiment</th><th>1D <sup>1</sup>H</th><th>1D <sup>19</sup>F</th><th><sup>1</sup>H-<sup>13</sup>C HSQC</th><th>IR</th></tr></thead><tbody><tr><td>Search method</td><td>Similarity</td><td>Peak</td><td>Peak</td><td>Similarity</td></tr><tr><td>HQI algorithm</td><td>Absolute</td><td>Peak</td><td>Peak</td><td>Euclidian</td></tr><tr><td>HQI limit (%)</td><td>70</td><td>70</td><td>a</td><td>50</td></tr><tr><td>LF (ppm)</td><td>b</td><td>0.3</td><td>0.1/1</td><td>b</td></tr></tbody></table>

Note: a: No HQI threshold for HSQC comparison because the two possibilities are only 0 or 100. b: No LF for similarity search, LF is defined at 0.1 ppm in the <sup>1</sup>H dimension and 1 ppm in the <sup>13</sup>C dimension for HSQC comparison.

Abbreviations: HQI, Hit Quality Index; HSQC, heteronuclear single quantum coherence; IR, infrared; LF, looseness factor.

details of the optimization process being provided in the SI (see S4 and S5).

As described above, the workflow starting point is the 1D <sup>19</sup>F analysis to determine if the molecule possesses one or several <sup>19</sup>F (Figure 2). Then, comparison between the measured HSQC spectrum and the database is carried out. The result of this comparison dictates the choice between the identification and elucidation parts of the workflow. When the HSQC pulse sequence gives a single match for the unknown structure, the leftmost part of the workflow in Figure 2 is used, corresponding to identification. Nevertheless, according to the SWGDRUG guidelines, one needs to ascertain the structure of the NPS with another technique. This is done by recording the IR spectrum. There is also a small chance that HSQC provides several matches or that not all signals are assigned. In this case, both 1D <sup>1</sup>H NMR and IR can be used to determine the most probable structure. The rightmost side of the workflow displayed in Figure 2 is dedicated to cases where the HSQC comparison returns no match. In such scenario, one would try to pinpoint the closest structure present inside the database, in order to provide structural insights into the molecular structure of the unknown NPS. To this end, the HSQC comparison is performed again with less restrictive LF values which are incremented until one match is obtained. This can help in determining the family of the unknown structure or part of its skeleton. Finally, IR and 1D <sup>1</sup>H NMR can be used to refine the proposed structure which can be validated by comparison between the experimental HSQC spectrum and a 2D HSQC spectrum predicted by the software. Of course, full elucidation with a database can remain difficult, especially when considering the limited performance of benchtop NMR. There will be cases where conventional structure elucidation at high field, relying on a complete set of 2D experiments (including correlation spectroscopy and heteronuclear multiple bond correlation), will be required to determine the exact structure. However, the present approach is directly accessible and advantageously requires a limited level of expertise compared with a full traditional structure elucidation. Therefore, it may prove useful as an initial structure elucidation tool, able of sorting out a large number of cases in a time compatible with the requirements of the police services.

## 3.3 | Database validation on real cases

A blind validation of the workflow was carried out on a set of six unknown seizures (numbered from 1 to 6) to evaluate if the database would be able to cope with concrete cases and could be routinely used in practice. All relevant data associated to these cases are given in supplementary information.

Sample no. 1 gives one signal on the 1D <sup>19</sup>F spectrum. AMB-FUBINACA (HQI = 100), AB-FUBINACA (100), and MDMB-FUBINACA (100) were returned as matches by comparison with the fluorine spectra present in the database. It indicates that the structure is part of the cannabinoid family and more specifically of the FUBINACA subfamily. Then, HSQC comparison provided only one hit, namely, AB-FUBINACA (100). The structure was then easily validated
less resolved technique to compare the spectral patterns of the unknown molecule with those included in the database.

## 3.2 | Multi-technique workflow

The choice of the analytical workflow is essential regarding the complexity of the sample embedding NPS. The strategy to build for identifying and elucidating NPS in seized samples has to face difficulties met to detect and quantify NPS even in mixture.

For a more comprehensive and integrated drug investigation, different multi-methods workflows have been reported and allow for elucidating NPS¹¹. They all combine multi-analytical techniques like Raman, DSC, IR, NMR, MS, and coupled chromatographic techniques with chemoinformatic tools and commercial, in-house databases for assisting spectral assignment, identification, and elucidation.

An applicative case where seized blotter paper²⁰ was analyzed by GC-MS, UHPLC-HRMS as well as 1D and 2D NMR in combination with multi-methods chemometric tools like ACD/Labs, Mass Hunter, and XCalibur. This workflow led to the identification of a LSD homologous product. In the case of seized NPS in mixture with herbs, the same workflow was utilized after NPS extraction and successfully revealed a new synthetic cannabinoid.¹⁸

Since high-field NMR and HRMS are expensive instruments and not owned by a majority of forensic laboratories, there is a need for cheaper and more accessible alternative workflows. We propose 1D and 2D optimized experiments implemented on liquid-state benchtop instruments in combination with IR and dedicated associated databases to identify and also elucidate NPS. Note that our approach aims at being quite general and may not be adapted to NPS hidden or sprayed in materials that would require a more specific strategy.

For our approach, an experimental database was created with 57 entries, mainly coming from cathinones, cannabinoids, amphetamines, arylcyclohexylamines, and fentanyloids (see Table S1). Before elaborating the multi-technique workflow, three model molecules were used to determine the best search method, the optimum algorithm for HQI calculation, the most adapted HQI threshold, as well as the optimum LF value. Table 2 summarizes the optimized parameters,

<table><caption>TABLE 2 Selected parameters for spectral comparison with ACD/Labs software</caption><thead><tr><th>Experiment</th><th>1D <sup>1</sup>H</th><th>1D <sup>19</sup>F</th><th><sup>1</sup>H-<sup>13</sup>C HSQC</th><th>IR</th></tr></thead><tbody><tr><td>Search method</td><td>Similarity</td><td>Peak</td><td>Peak</td><td>Similarity</td></tr><tr><td>HQI algorithm</td><td>Absolute</td><td>Peak</td><td>Peak</td><td>Euclidian</td></tr><tr><td>HQI limit (%)</td><td>70</td><td>70</td><td>a</td><td>50</td></tr><tr><td>LF (ppm)</td><td>b</td><td>0.3</td><td>0.1/1</td><td>b</td></tr></tbody></table>

Note: a: No HQI threshold for HSQC comparison because the two possibilities are only 0 or 100. b: No LF for similarity search, LF is defined at 0.1 ppm in the <sup>1</sup>H dimension and 1 ppm in the <sup>13</sup>C dimension for HSQC comparison.

Abbreviations: HQI, Hit Quality Index; HSQC, heteronuclear single quantum coherence; IR, infrared; LF, looseness factor.

details of the optimization process being provided in the SI (see S4 and S5).

As described above, the workflow starting point is the 1D <sup>19</sup>F analysis to determine if the molecule possesses one or several <sup>19</sup>F (Figure 2). Then, comparison between the measured HSQC spectrum and the database is carried out. The result of this comparison dictates the choice between the identification and elucidation parts of the workflow. When the HSQC pulse sequence gives a single match for the unknown structure, the leftmost part of the workflow in Figure 2 is used, corresponding to identification. Nevertheless, according to the SWGDRUG guidelines, one needs to ascertain the structure of the NPS with another technique. This is done by recording the IR spectrum. There is also a small chance that HSQC provides several matches or that not all signals are assigned. In this case, both 1D <sup>1</sup>H NMR and IR can be used to determine the most probable structure. The rightmost side of the workflow displayed in Figure 2 is dedicated to cases where the HSQC comparison returns no match. In such scenario, one would try to pinpoint the closest structure present inside the database, in order to provide structural insights into the molecular structure of the unknown NPS. To this end, the HSQC comparison is performed again with less restrictive LF values which are incremented until one match is obtained. This can help in determining the family of the unknown structure or part of its skeleton. Finally, IR and 1D <sup>1</sup>H NMR can be used to refine the proposed structure which can be validated by comparison between the experimental HSQC spectrum and a 2D HSQC spectrum predicted by the software. Of course, full elucidation with a database can remain difficult, especially when considering the limited performance of benchtop NMR. There will be cases where conventional structure elucidation at high field, relying on a complete set of 2D experiments (including correlation spectroscopy and heteronuclear multiple bond correlation), will be required to determine the exact structure. However, the present approach is directly accessible and advantageously requires a limited level of expertise compared with a full traditional structure elucidation. Therefore, it may prove useful as an initial structure elucidation tool, able of sorting out a large number of cases in a time compatible with the requirements of the police services.

## 3.3 | Database validation on real cases

A blind validation of the workflow was carried out on a set of six unknown seizures (numbered from 1 to 6) to evaluate if the database would be able to cope with concrete cases and could be routinely used in practice. All relevant data associated to these cases are given in supplementary information.

Sample no. 1 gives one signal on the 1D <sup>19</sup>F spectrum. AMB-FUBINACA (HQI = 100), AB-FUBINACA (100), and MDMB-FUBINACA (100) were returned as matches by comparison with the fluorine spectra present in the database. It indicates that the structure is part of the cannabinoid family and more specifically of the FUBINACA subfamily. Then, HSQC comparison provided only one hit, namely, AB-FUBINACA (100). The structure was then easily validated
by its IR spectrum. Indeed, for IR, the four top matches all belong to the cannabinoid family: AB-FUBINACA (89.43), AB-CHMINACA (66.56), ADB-CHMINACA (58.88), and APINACA (53.52), which unambiguously validates AB-FUBINACA as the actual sample no. 1 (see S6).

Sample no. 2 followed a similar workflow even if no signal appeared in the 1D ¹⁹F spectrum, indicating the absence of fluorine. Only one match was found with HSQC, namely, acetyl fentanyl (AF, HQI = 100). The structure was next confirmed by IR since AF (86.54) is the top match, its HQI largely exceeding the one of the second hit, PHENF with HQI of 58.85 (see S7).

Likewise, sample no. 3, that returned no match based on the 1D ¹⁹F spectrum, provided a single hit after the HSQC comparison, 2C-E (100). However, when the ¹H spectrum in the 2D projection was manually checked, it appeared that two signals in the aromatic region (6.9 and 7.3 ppm) could not be assigned to this structure. A comparison with 1D ¹H experiments confirmed that those two signals belong to another NPS, 2C-I (HQI 93.70), that was also present in the same sample and whose structure was similar to 2C-E (HQI 95.85). In short, 2C-I could be identified by the direct HSQC comparison whereas the presence of 2C-E was found by the ¹H NMR comparison and confirmed by visual inspection of the operator (see S8).

These first three examples highlight the efficiency of the identification component of the workflow since the database was able to systematically provide the structure of NPS obtained from seizures. In addition, the identification of a (simple) mixture could be achieved thanks to our workflow.

The three next samples correspond to cases where the elucidation part of the workflow was applied, because the identification failed to provide a clear hit. As illustrated in Figure 3, sample no. 4 provided no signal in 1D ¹⁹F spectrum, and we moved directly to the HSQC comparison. No match could be obtained using a LF of 0.1 ppm for ¹H and 1 ppm for ¹³C. Next, we tried to find structures that would be close to the unknown one (elucidation by homology) by comparing the experimental HSQC spectrum with the database while concomitantly incrementing LF by steps of 0.1 ppm for ¹H and 1 ppm

for ¹³C. With such procedure, we obtained O-PCE as a first match, with a LF of 0.6 ppm for ¹H and 6 ppm for ¹³C. Then, a 1D ¹H comparison was performed and the first two structures were O-PCE (91.04) and MXE (87.33). All ¹H signals had been integrated, and we found that the aromatic region integrated to four protons in agreement with a di-substituted phenyl ring, with no major modification of the spectral pattern between the unknown structure and O-PCE. Therefore, we attempted to better clarify the nature of substituting group on the aromatic region with the IR experiment. Again, O-PCE came as a first hit (57.80) and a band at 710 cm⁻¹ was found as a possible indicator of the presence of chlorine atom(s) in the structure (see S9). Due to the high number of IR bands and their overlap, the first hit for IR was lower than the first NMR hit for O-PCE.

Three different positional isomers were possible in that case, with the CI center in *ortho* (Figure 3c), *meta* (Figure 3b), or *para* (Figure 3a), and prediction was therefore needed to identify the structure. In the HSQC spectrum three correlations could be peak picked 8.02/132.4, 7.58/127.9, and 7.60/131.7 ppm, the latter apparently corresponding to two different protons (the signal intensity was much higher for this one than for the two others correlations). With those information at hand, *para* substitution could already be discarded as only two correlations would be expected for this isomer. Next, the Predictor ACD/Labs module was used to predict the chemical shift of the three isomers (Figure 4). The closest match with the actual measured spectrum of sample no. 4 was obtained for the *ortho* isomer; this molecule is called NEK (Figure 3c). For further validation, the structure was later confirmed by high-field experiments performed at 700 MHz.

The workflow for sample no. 5 started directly with HSQC since no signal was observed in the 1D ¹⁹F spectrum. No match could be found with the identification parameters given in Table 2. As above, we incremented LF until 0.5 ppm for ¹H and 5 ppm for ¹³C, values leading to PV9 as the closest structure. The HSQC spectra of the unknown compound and of PV9 are reported in S10. In the aromatic region PV9 presents three correlations representative of a monosubstituted species, whereas substance no. 5 shows two correlations, a pattern consistent with a *para* bisubstitution (CH₃ correlation at 0.73

<figure><img src="image_4.png" /><figcaption>**FIGURE 3** Workflow for the determination of unknown sample no. 4, finally yielding NEK as a structure (a) *para* substitution, (b) *meta* substitution, and (c) *ortho* substitution [Colour figure can be viewed at wileyonlinelibrary.com]</figcaption></figure>
by its IR spectrum. Indeed, for IR, the four top matches all belong to the cannabinoid family: AB-FUBINACA (89.43), AB-CHMINACA (66.56), ADB-CHMINACA (58.88), and APINACA (53.52), which unambiguously validates AB-FUBINACA as the actual sample no. 1 (see S6).

Sample no. 2 followed a similar workflow even if no signal appeared in the 1D ¹⁹F spectrum, indicating the absence of fluorine. Only one match was found with HSQC, namely, acetyl fentanyl (AF, HQI = 100). The structure was next confirmed by IR since AF (86.54) is the top match, its HQI largely exceeding the one of the second hit, PHENF with HQI of 58.85 (see S7).

Likewise, sample no. 3, that returned no match based on the 1D ¹⁹F spectrum, provided a single hit after the HSQC comparison, 2C-E (100). However, when the ¹H spectrum in the 2D projection was manually checked, it appeared that two signals in the aromatic region (6.9 and 7.3 ppm) could not be assigned to this structure. A comparison with 1D ¹H experiments confirmed that those two signals belong to another NPS, 2C-I (HQI 93.70), that was also present in the same sample and whose structure was similar to 2C-E (HQI 95.85). In short, 2C-I could be identified by the direct HSQC comparison whereas the presence of 2C-E was found by the ¹H NMR comparison and confirmed by visual inspection of the operator (see S8).

These first three examples highlight the efficiency of the identification component of the workflow since the database was able to systematically provide the structure of NPS obtained from seizures. In addition, the identification of a (simple) mixture could be achieved thanks to our workflow.

The three next samples correspond to cases where the elucidation part of the workflow was applied, because the identification failed to provide a clear hit. As illustrated in Figure 3, sample no. 4 provided no signal in 1D ¹⁹F spectrum, and we moved directly to the HSQC comparison. No match could be obtained using a LF of 0.1 ppm for ¹H and 1 ppm for ¹³C. Next, we tried to find structures that would be close to the unknown one (elucidation by homology) by comparing the experimental HSQC spectrum with the database while concomitantly incrementing LF by steps of 0.1 ppm for ¹H and 1 ppm

for ¹³C. With such procedure, we obtained O-PCE as a first match, with a LF of 0.6 ppm for ¹H and 6 ppm for ¹³C. Then, a 1D ¹H comparison was performed and the first two structures were O-PCE (91.04) and MXE (87.33). All ¹H signals had been integrated, and we found that the aromatic region integrated to four protons in agreement with a di-substituted phenyl ring, with no major modification of the spectral pattern between the unknown structure and O-PCE. Therefore, we attempted to better clarify the nature of substituting group on the aromatic region with the IR experiment. Again, O-PCE came as a first hit (57.80) and a band at 710 cm⁻¹ was found as a possible indicator of the presence of chlorine atom(s) in the structure (see S9). Due to the high number of IR bands and their overlap, the first hit for IR was lower than the first NMR hit for O-PCE.

Three different positional isomers were possible in that case, with the CI center in *ortho* (Figure 3c), *meta* (Figure 3b), or *para* (Figure 3a), and prediction was therefore needed to identify the structure. In the HSQC spectrum three correlations could be peak picked 8.02/132.4, 7.58/127.9, and 7.60/131.7 ppm, the latter apparently corresponding to two different protons (the signal intensity was much higher for this one than for the two others correlations). With those information at hand, *para* substitution could already be discarded as only two correlations would be expected for this isomer. Next, the Predictor ACD/Labs module was used to predict the chemical shift of the three isomers (Figure 4). The closest match with the actual measured spectrum of sample no. 4 was obtained for the *ortho* isomer; this molecule is called NEK (Figure 3c). For further validation, the structure was later confirmed by high-field experiments performed at 700 MHz.

The workflow for sample no. 5 started directly with HSQC since no signal was observed in the 1D ¹⁹F spectrum. No match could be found with the identification parameters given in Table 2. As above, we incremented LF until 0.5 ppm for ¹H and 5 ppm for ¹³C, values leading to PV9 as the closest structure. The HSQC spectra of the unknown compound and of PV9 are reported in S10. In the aromatic region PV9 presents three correlations representative of a monosubstituted species, whereas substance no. 5 shows two correlations, a pattern consistent with a *para* bisubstitution (CH₃ correlation at 0.73

<figure><img src="image_4.png" /><figcaption>**FIGURE 3** Workflow for the determination of unknown sample no. 4, finally yielding NEK as a structure (a) *para* substitution, (b) *meta* substitution, and (c) *ortho* substitution [Colour figure can be viewed at wileyonlinelibrary.com]</figcaption></figure>
<figure><img src="image_5.png" /><figcaption>FIGURE 4 Comparison between predicted (a) para, (b) meta, (c) ortho, and (d) experimental heteronuclear single quantum coherence (HSQC) spectra for unknown sample no. 4 [Colour figure can be viewed at wileyonlinelibrary.com]</figcaption></figure>

and 13.47 ppm). Correlations for protons number 9, 12, 13, 14, and 15 are identical in the two spectra. A correlation at 1.98 ppm ($^1$H) and 29 ppm ($^{13}$C) is present for both structures, this CH$_2$ being the beginning of the side chain. Then, two CH$_2$ and one CH$_3$ remained unasigned with this HSQC comparison. Then, a 1D $^1$H NMR comparison was achieved to obtain more information on the nature of the lateral chain. Two different chain lengths were present in the three first hit: propyl (PVP, MDPV) and hexyl (PV9). However, visually, those spectra strongly differed from the one of our unknown structure. In addition, unknown compound no. 5 could not contain only methyl or ethyl groups since the NMR spectra would respectively show a doublet or a triplet coupled to a quadruplet in these two cases (see S10 Figure S1). IR comparison was also made, and the same structures as before were identified as the closest to our structure. Finally, the HSQC prediction was done considering a butyl side chain (3 CH$_2$ and 1 CH$_3$) and it returned a map closely resembling one of the unknown samples (Figure S10a). Thus, the structure was identified as MPHP with a butyl side chain, and this was confirmed by 700 MHz experiments (see S10).

The two previous cases demonstrate that elucidation with the help of a database can be carried out relying on benchtop NMR combined with IR spectroscopy. The two new structures were eventually added to the database. Indeed, the database can be improved while previously unknown NPS are discovered, allowing a continuous improvement of its identification and elucidation abilities.

For the last case, sample no. 6, the usual workflow was followed. The 1D $^{19}$F comparison was not performed since no fluorine signal could be observed on the spectrum. Then, the comparison with HSQC provided no match with the LF factors given in Table 2. Therefore, LF was incremented until a match was obtained with a LF of 0.5 for $^1$H and 5 ppm for $^{13}$C. The proposed structure with this comparison is

NEK, the sample that was added to the database just after processing sample no. 4. To confirm that structure, 1D $^1$H comparison was used and provided NEK (92.50), O-PCE (89.72), and HEXEN (86.56) as the three best matches. However, when spectral patterns of those two spectra were compared, it could be seen that additional signals were present in sample no. 6 (see S11). The presence of a mixture explains why LF needed to be increased in order to obtain a match with HSQC. The two components of this mixture were determined and validated at 700 MHz (see S11): It turns out that the sample is a mixture of NEK and MDMC. This shows that although benchtop NMR failed to identify all the analytes in the mixture, the major component (NEK) could be identified relying on our workflow.

To summarize, six cases were examined using the previously constructed workflow. In the first three cases, AF, AB-FUBINACA, and a mixture of 2C-I and 2C-E, were straightforwardly identified with the help of their reference spectra in the database. The two next samples were elucidated with the help of the database to determine the structures of NEK and MPHP (see S10 and S11). Finally, the last sample could not be fully characterized by our workflow which nevertheless identified NEK as the major component of the mixture. The strategy based on structure similarity was also reported with LC/ESI/HRMS in the field of a non-targeted screening.²¹ This approach also made it possible to access semiquantitative data of main interest for tracking drugs.

## 3.4 Quantification

To complete our study, the purity of the fully characterized NPS was measured by quantitative benchtop NMR. Purity determination is an additional data to trace the manufacturers and to further dismantle
<figure><img src="image_5.png" /><figcaption>FIGURE 4 Comparison between predicted (a) para, (b) meta, (c) ortho, and (d) experimental heteronuclear single quantum coherence (HSQC) spectra for unknown sample no. 4 [Colour figure can be viewed at wileyonlinelibrary.com]</figcaption></figure>

and 13.47 ppm). Correlations for protons number 9, 12, 13, 14, and 15 are identical in the two spectra. A correlation at 1.98 ppm ($^{1}$H) and 29 ppm ($^{13}$C) is present for both structures, this CH$_{2}$ being the beginning of the side chain. Then, two CH$_{2}$ and one CH$_{3}$ remained unasigned with this HSQC comparison. Then, a 1D $^{1}$H NMR comparison was achieved to obtain more information on the nature of the lateral chain. Two different chain lengths were present in the three first hit: propyl (PVP, MDPV) and hexyl (PV9). However, visually, those spectra strongly differed from the one of our unknown structure. In addition, unknown compound no. 5 could not contain only methyl or ethyl groups since the NMR spectra would respectively show a doublet or a triplet coupled to a quadruplet in these two cases (see S10 Figure S1). IR comparison was also made, and the same structures as before were identified as the closest to our structure. Finally, the HSQC prediction was done considering a butyl side chain (3 CH$_{2}$ and 1 CH$_{3}$) and it returned a map closely resembling one of the unknown samples (Figure S10a). Thus, the structure was identified as MPHP with a butyl side chain, and this was confirmed by 700 MHz experiments (see S10).

The two previous cases demonstrate that elucidation with the help of a database can be carried out relying on benchtop NMR combined with IR spectroscopy. The two new structures were eventually added to the database. Indeed, the database can be improved while previously unknown NPS are discovered, allowing a continuous improvement of its identification and elucidation abilities.

For the last case, sample no. 6, the usual workflow was followed. The 1D $^{19}$F comparison was not performed since no fluorine signal could be observed on the spectrum. Then, the comparison with HSQC provided no match with the LF factors given in Table 2. Therefore, LF was incremented until a match was obtained with a LF of 0.5 for $^{1}$H and 5 ppm for $^{13}$C. The proposed structure with this comparison is

NEK, the sample that was added to the database just after processing sample no. 4. To confirm that structure, 1D $^{1}$H comparison was used and provided NEK (92.50), O-PCE (89.72), and HEXEN (86.56) as the three best matches. However, when spectral patterns of those two spectra were compared, it could be seen that additional signals were present in sample no. 6 (see S11). The presence of a mixture explains why LF needed to be increased in order to obtain a match with HSQC. The two components of this mixture were determined and validated at 700 MHz (see S11): It turns out that the sample is a mixture of NEK and MDMC. This shows that although benchtop NMR failed to identify all the analytes in the mixture, the major component (NEK) could be identified relying on our workflow.

To summarize, six cases were examined using the previously constructed workflow. In the first three cases, AF, AB-FUBINACA, and a mixture of 2C-I and 2C-E, were straightforwardly identified with the help of their reference spectra in the database. The two next samples were elucidated with the help of the database to determine the structures of NEK and MPHP (see S10 and S11). Finally, the last sample could not be fully characterized by our workflow which nevertheless identified NEK as the major component of the mixture. The strategy based on structure similarity was also reported with LC/ESI/HRMS in the field of a non-targeted screening.²¹ This approach also made it possible to access semiquantitative data of main interest for tracking drugs.

## 3.4 Quantification

To complete our study, the purity of the fully characterized NPS was measured by quantitative benchtop NMR. Purity determination is an additional data to trace the manufacturers and to further dismantle
drug networks. It is generally performed by GC-FID or LC-DAD analysis, but these techniques require finding the corresponding reference material, which is not systematically accessible and can also be prohibitively expensive. Therefore, NMR appears as a valuable alternative because, on the one hand, quantification can be achieved without physical separation of the sample components, and on the other hand, multiple compounds can be precisely quantified with a single internal or external reference, which does not necessarily have to be close to the analyte (in contrast to chromatography). An internal reference was not used here since avoiding overlap with NMR signals was beyond reach for the benchtop apparatus. Therefore, the choice was made to determine the purity using TSP as an external reference. Purity determination was carried out relying on the following equation:

$$P_{\text{NPS}} = P_{\text{TSP}} \times \left(\frac{A_x}{A_{\text{TSP}}}\right) \times \left(\frac{N_{\text{TSP}}}{N_x}\right) \times \left(\frac{M_x}{M_{\text{TSP}}}\right) \times \left(\frac{m_{\text{TSP}}}{m_x}\right), \quad (4)$$

where $P_{\text{TSP}}$ is the purity of the external standard TSP, $A_x$ the integral of the analyte signal, $A_{\text{TSP}}$ the integral of the reference signal, $N_{\text{TSP}}$ the number of reference protons, $N_x$ the number of analyte protons, $M_x$ the molecular weight of the analyte (mg/mol), $M_{\text{TSP}}$ the molecular weight of the reference (mg/mol), $m_{\text{TSP}}$ the mass of reference compound (g), and $m_x$ the mass of sample (g).

Then, the precision of the purity was determined on five successive spectra for each sample, and it was calculated with the following formula.

$$P(\%) = CV(A_x) + CV(A_{\text{TSP}}) + \left(\frac{\Delta m_x}{m_x}\right) \times 100 + \left(\frac{\Delta m_{\text{TSP}}}{m_{\text{TSP}}}\right) \times 100, \quad (5)$$

where $CV(A_x)$ is the standard deviation of the reference integral determined with five repetitions, $CV(A_{\text{TSP}})$ is the standard deviation of the analyte integral determined with five repetitions, $\Delta m_x$ and $\Delta m_{\text{TSP}}$ equal to 0.01 mg, and $m_x$ and $m_{\text{TSP}}$ have been defined above.

NMR quantification with 1D ¹H benchtop NMR can be difficult, because it is necessary to find at least one signal per compound which is not overlapped with others and where the baseline reaches zero on both sides to reach an accurate (true and precise) quantification. Thus, purity of the substance, nature of the substance, and sample/matrix interferences could impact the accuracy of the quantification. For each of the previous samples, the signal which was

the least overlapped with others was used (see S12). Calculated purity values are displayed in Table 3. The two pure identified structures have a purity higher than 90% with a precision of approximately 2%. For the identified mixture, the purities of 2C-I and 2C-E were determined as 66% and 21%, respectively, with again a precision close to near 2%. MPHP was determined to be pure at more than 90%. Finally, the purity of NEK was determined to be much higher than 100%, and this awkward result is related to the absence of signal free of overlapping in the 1D spectrum. Signals integrated for the determination of purity belonged to the aromatic region which were presumably overlapped with signals of a minor component that could not be detected with a magnetic field under 100 MHz (some signals could be detected but not identified at 700 MHz). With a benchtop NMR spectrometer, quantitative 1D NMR could be limited when seized samples are complex, for example, embedding different types of compounds and impurities. Thus, in these cases 2D quantitative NMR could be relevant, as well as the combination with other chromatographic-based separations in case of highly complex mixtures.

## 4 | CONCLUSION

An integrative workflow was created to help identifying NPS. It includes benchtop 1D and 2D NMR, IR, and prediction together with an experimental database to characterize NPS according to the directives of the SWGDRUGS workgroup. The database includes data for 57 compounds but will be enriched with more samples in the near future. The efficiency of the workflow was evaluated on six seizures. Most structures were unambiguously identified or elucidated, even in the case of mixtures. This evaluation also highlighted the limitation of the workflow when significant overlap between peaks occurs. In some cases, an expert or habituated eye will most likely be needed, but we are confident that the workflow presented in this study would facilitate the identification process.

Additionally, purity was also determined for each sample thanks to 1D ¹H NMR. Precise purities were obtained except for one case where the NMR signal was heavily overlapped with the one of another molecule. To the best of our knowledge, this work stands as the first example of using benchtop NMR in an integrated workflow for identifying, elucidating and quantifying NPS. This workflow can serve as an alternative to the use of GC-MS: It could allow forensic laboratories to efficiently characterize and quantify a wider panel of drugs and especially NPS.

The present workflow and database work well for simple elucidation cases, for more difficult ones, the creation of a bigger database including more reference structures is needed. In addition, the proposed workflow could be combined with the current workflows already used by the police, by including GC-MS in the database to maximize its potential. Another appealing perspective would be to combine those developments with machine learning, to create an even more automated and reliable elucidation procedure.

<table><caption>TABLE 3 Purity obtained for NPS samples 1 to 6 and associated precision</caption><thead><tr><td>Sample no.</td><td>Molecule</td><td>Purity (%)</td><td>Precision (%)</td></tr></thead><tbody><tr><td>1</td><td>AB-FUBINACA</td><td>90.2</td><td>1.8</td></tr><tr><td>2</td><td>AF</td><td>99.0</td><td>2.2</td></tr><tr><td>3</td><td>2C-E</td><td>66.4</td><td>2.4</td></tr><tr><td>3</td><td>2C-I</td><td>21.3</td><td>3.1</td></tr><tr><td>4</td><td>NEK</td><td>113.5</td><td>2.6</td></tr><tr><td>5</td><td>MPHP</td><td>94.8</td><td>4.5</td></tr></tbody></table>
drug networks. It is generally performed by GC-FID or LC-DAD analysis, but these techniques require finding the corresponding reference material, which is not systematically accessible and can also be prohibitively expensive. Therefore, NMR appears as a valuable alternative because, on the one hand, quantification can be achieved without physical separation of the sample components, and on the other hand, multiple compounds can be precisely quantified with a single internal or external reference, which does not necessarily have to be close to the analyte (in contrast to chromatography). An internal reference was not used here since avoiding overlap with NMR signals was beyond reach for the benchtop apparatus. Therefore, the choice was made to determine the purity using TSP as an external reference. Purity determination was carried out relying on the following equation:

$$P_{\text{NPS}} = P_{\text{TSP}} \times \left(\frac{A_x}{A_{\text{TSP}}}\right) \times \left(\frac{N_{\text{TSP}}}{N_x}\right) \times \left(\frac{M_x}{M_{\text{TSP}}}\right) \times \left(\frac{m_{\text{TSP}}}{m_x}\right), \quad (4)$$

where $P_{\text{TSP}}$ is the purity of the external standard TSP, $A_x$ the integral of the analyte signal, $A_{\text{TSP}}$ the integral of the reference signal, $N_{\text{TSP}}$ the number of reference protons, $N_x$ the number of analyte protons, $M_x$ the molecular weight of the analyte (mg/mol), $M_{\text{TSP}}$ the molecular weight of the reference (mg/mol), $m_{\text{TSP}}$ the mass of reference compound (g), and $m_x$ the mass of sample (g).

Then, the precision of the purity was determined on five successive spectra for each sample, and it was calculated with the following formula.

$$P(\%) = CV(A_x) + CV(A_{\text{TSP}}) + \left(\frac{\Delta m_x}{m_x}\right) \times 100 + \left(\frac{\Delta m_{\text{TSP}}}{m_{\text{TSP}}}\right) \times 100, \quad (5)$$

where $CV(A_x)$ is the standard deviation of the reference integral determined with five repetitions, $CV(A_{\text{TSP}})$ is the standard deviation of the analyte integral determined with five repetitions, $\Delta m_x$ and $\Delta m_{\text{TSP}}$ equal to 0.01 mg, and $m_x$ and $m_{\text{TSP}}$ have been defined above.

NMR quantification with 1D ¹H benchtop NMR can be difficult, because it is necessary to find at least one signal per compound which is not overlapped with others and where the baseline reaches zero on both sides to reach an accurate (true and precise) quantification. Thus, purity of the substance, nature of the substance, and sample/matrix interferences could impact the accuracy of the quantification. For each of the previous samples, the signal which was

the least overlapped with others was used (see S12). Calculated purity values are displayed in Table 3. The two pure identified structures have a purity higher than 90% with a precision of approximately 2%. For the identified mixture, the purities of 2C-I and 2C-E were determined as 66% and 21%, respectively, with again a precision close to near 2%. MPHP was determined to be pure at more than 90%. Finally, the purity of NEK was determined to be much higher than 100%, and this awkward result is related to the absence of signal free of overlapping in the 1D spectrum. Signals integrated for the determination of purity belonged to the aromatic region which were presumably overlapped with signals of a minor component that could not be detected with a magnetic field under 100 MHz (some signals could be detected but not identified at 700 MHz). With a benchtop NMR spectrometer, quantitative 1D NMR could be limited when seized samples are complex, for example, embedding different types of compounds and impurities. Thus, in these cases 2D quantitative NMR could be relevant, as well as the combination with other chromatographic-based separations in case of highly complex mixtures.

## 4 | CONCLUSION

An integrative workflow was created to help identifying NPS. It includes benchtop 1D and 2D NMR, IR, and prediction together with an experimental database to characterize NPS according to the directives of the SWGDRUGS workgroup. The database includes data for 57 compounds but will be enriched with more samples in the near future. The efficiency of the workflow was evaluated on six seizures. Most structures were unambiguously identified or elucidated, even in the case of mixtures. This evaluation also highlighted the limitation of the workflow when significant overlap between peaks occurs. In some cases, an expert or habituated eye will most likely be needed, but we are confident that the workflow presented in this study would facilitate the identification process.

Additionally, purity was also determined for each sample thanks to 1D ¹H NMR. Precise purities were obtained except for one case where the NMR signal was heavily overlapped with the one of another molecule. To the best of our knowledge, this work stands as the first example of using benchtop NMR in an integrated workflow for identifying, elucidating and quantifying NPS. This workflow can serve as an alternative to the use of GC-MS: It could allow forensic laboratories to efficiently characterize and quantify a wider panel of drugs and especially NPS.

The present workflow and database work well for simple elucidation cases, for more difficult ones, the creation of a bigger database including more reference structures is needed. In addition, the proposed workflow could be combined with the current workflows already used by the police, by including GC-MS in the database to maximize its potential. Another appealing perspective would be to combine those developments with machine learning, to create an even more automated and reliable elucidation procedure.

<table><caption>TABLE 3 Purity obtained for NPS samples 1 to 6 and associated precision</caption><thead><tr><td>Sample no.</td><td>Molecule</td><td>Purity (%)</td><td>Precision (%)</td></tr></thead><tbody><tr><td>1</td><td>AB-FUBINACA</td><td>90.2</td><td>1.8</td></tr><tr><td>2</td><td>AF</td><td>99.0</td><td>2.2</td></tr><tr><td>3</td><td>2C-E</td><td>66.4</td><td>2.4</td></tr><tr><td>3</td><td>2C-I</td><td>21.3</td><td>3.1</td></tr><tr><td>4</td><td>NEK</td><td>113.5</td><td>2.6</td></tr><tr><td>5</td><td>MPHP</td><td>94.8</td><td>4.5</td></tr></tbody></table>
## ACKNOWLEDGEMENTS

The authors acknowledge Magritek, and especially Hélène Freichels, for her help with the autosampler and Craig Eccles for his help with the SpinsolveExpert software. The authors are also indebted to Benoît Charrier for his help in the maintenance of the spectrometers. Shahriar Jahanbakht (ACD/Labs) is warmly acknowledged for his clear explanations about the various features of the ACD/Labs software. We thank the Finnish National Bureau of Investigation for giving us some samples to complete our database. The authors from the CEISAM lab acknowledge the French National Infrastructure for Metabolomics and Fluxomics MetaboHUB-ANR-11-INBS-0010 (www.metabohub.fr) and the Corsaire metabolomics core facility (Biogenouest). This work includes NMR experiments carried out on the CEISAM NMR platform. J.F. thanks Sandrine Bouchet for an unfailing assistance.

## CONFLICT OF INTEREST

None.

## DATA AVAILABILITY STATEMENT

Data available on request from the authors.

## ORCID

Fabrice Besacier https://orcid.org/0000-0003-4990-8901

<figure><img src="image_6.png" /></figure>

Jonathan Farjon https://orcid.org/0000-0003-0707-2026

<figure><img src="image_7.png" /></figure>

## REFERENCES

1. Nouveaux produits de synthèse - Synthèse des connaissances - OFDT. https://www.ofdt.fr/produits-et-addictions/de-z/nouveaux-produits-de-synthèse/ Accessed July 16, 2021.

2. Actualités européennes - OFDT. https://www.ofdt.fr/europe-et-international/actualites-europeennes/ Accessed July 15, 2021.

3. Lukić V, Micić R, Arsić B, Nedović B, Radosavljević Ž. Overview of the major classes of new psychoactive substances, psychoactive effects, analytical determination and conformational analysis of selected ille- gal drugs. Open Chemistry. 2021;19(1):60-106. doi:10.1515/chem- 2021-0196

4. Elyashberg M. Identification and structure elucidation by NMR spec- troscopy. TrAC. 2015;69:88-97. doi:10.1016/j.trac.2015.02.014

5. Elie MP, Elie LE, Baron MG. Keeping pace with NPS releases: fast GC- MS screening of legal high products. Drug Test Anal. 2013;5(5): 281-290. doi:10.1002/dta.1434

6. Angerer V, Bisel P, Moosmann B, Westphal F, Auwärter V. Separation and structural characterization of the new synthetic cannabinoid JWH-018 cyclohexyl methyl derivative “NE-CHMIMO” using flash chromatography, GC-MS, IR and NMR spectroscopy. Forensic Sci Int. 2016;266:93-98. doi:10.1016/j.forsciint.2016.05.031

7. Westphal F, Junge T, Girreser U, Greibl W, Doering C. Mass, NMR and IR spectroscopic characterization of pentedrone and pentylone and identification of their isocathinone by-products. Forensic Sci Int. 2012;217(1):157-167. doi:10.1016/j.forsciint.2011.10.045

8. Böselt L, Dötzer R, Steiner S, Stritzinger M, Salzmann S, Riniker S. Determining the regiochemistry and relative stereochemistry of small and druglike molecules using an alignment algorithm for infrared spectra. Anal Chem. 2020;92(13):9124-9131. doi:10.1021/acs. analchem.0c01399

9. Blümich B. Low-field and benchtop NMR. J Magn Reson. 2019;306: 27-35. doi:10.1016/j.jmr.2019.07.030

10. Castaing-Cordier T, Bouillaud D, Farjon J, Giraudeau P. Chapter Four - Recent advances in benchtop NMR spectroscopy and its applica- tions. Ann Rep NMR Spectros; Webb, G. A., Ed.; Academic Press. 2021;103:191-258. PMID. doi:10.1016/bs.arnmr.2021.02.003

11. Tsochatzis ED, Giannopoulos G, Lopes JA, Guillou C. Drugs analysis and citizens safety. Workflows for chemical analysis and structural elucidation of new psychoactive substances in the EU. Saf Sci. 2021; 137:105126. doi:10.1016/j.ssci.2020.105126

12. Antonides LH, Brignall RM, Costello A, et al. Rapid identification of novel psychoactive and other controlled substances using low-Field ¹H NMR spectroscopy. ACS Omega. 2019;4(4):7103-7112. doi:10. 1021/acsomega.9b00302

13. Zhong Y, Huang K, Luo Q et al. The application of a desktop NMR spectrometer in drug analysis. https://www.hindawi.com/journals/ ijac/2018/3104569/ Accessed January 31, 2022.

14. SWGDRUG Approved Recommendations. http://www.swgdrug.org/ approved.htm Accessed September 30, 2021.

15. Jones LE, Stewart A, Peters KL, et al. Infrared and raman screening of seized novel psychoactive substances: a large scale study of >200 Samples. Analyst. 2016;141(3):902-909. doi:10.1039/C5AN02326B

16. Naqi HA, Husbands SM, Blagbrough IS. ¹H quantitative NMR and UHPLC-MS analysis of seized MDMA/NPS mixtures and tablets from night-club venues. Anal Methods. 2019;11(37):4795-4807. doi:10. 1039/C9AY01403A

17. Hussain JH, Gilbert N, Costello A, et al. Quantification of MDMA in seized tablets using benchtop ¹H NMR spectroscopy in the absence of internal standards. Forensic Chem. 2020;20:100263. doi:10.1016/j. forc.2020.100263

18. Tsochatzis E, Lopes JA, Reniero F, Holland M, Åberg J, Guillou C. Identification of 1-Butyl-Lysergic acid diethylamide (1B-LSD) in seized blotter paper using an integrated workflow of analytical techniques and chemo-informatics. Molecules. 2020;25(3):1-10. doi:10. 3390/molecules25030712

19. Brandt SD, Kavanagh PV, Westphal F, et al. Synthetic cannabinoid receptor agonists: Analytical profiles and development of QMPSB, QMMSB, QMPCB, 2F-QMPSB, QMiPSB, and SGT-233. Drug Test Anal. 2021;13(1):175-196. doi:10.1002/dta.2913

20. Tsochatzis ED, Alberto Lopes J, Holland MV, Reniero F, Palmieri G, Guillou C. Identification and analytical characterization of a novel syn- thetic cannabinoid-type substance in herbal material in Europe. Mole- cules. 2021;26(4):1-16. doi:10.3390/molecules26040793

21. Malm L, Palm E, Souihi A, Plassmann M, Liigand J, Kruve A. Guide to semi-quantitative non-targeted screening using LC/ESI/HRMS. Mole- cules. 2021;26(12):3524. doi:10.3390/molecules26123524

## SUPPORTING INFORMATION

Additional supporting information can be found online in the Support- ing Information section at the end of this article.

**How to cite this article:** Castaing-Cordier T, Benavides Restrepo A, Dubois D, et al. Characterization of new psychoactive substances by integrating benchtop NMR to multi-technique databases. Drug Test Anal. 2022;14(9): 1629-1638. doi:10.1002/dta.3332
## ACKNOWLEDGEMENTS

The authors acknowledge Magritek, and especially Hélène Freichels, for her help with the autosampler and Craig Eccles for his help with the SpinsolveExpert software. The authors are also indebted to Benoît Charrier for his help in the maintenance of the spectrometers. Shahriar Jahanbakht (ACD/Labs) is warmly acknowledged for his clear explanations about the various features of the ACD/Labs software. We thank the Finnish National Bureau of Investigation for giving us some samples to complete our database. The authors from the CEISAM lab acknowledge the French National Infrastructure for Metabolomics and Fluxomics MetaboHUB-ANR-11-INBS-0010 (www.metabohub.fr) and the Corsaire metabolomics core facility (Biogenouest). This work includes NMR experiments carried out on the CEISAM NMR platform. J.F. thanks Sandrine Bouchet for an unfailing assistance.

## CONFLICT OF INTEREST

None.

## DATA AVAILABILITY STATEMENT

Data available on request from the authors.

## ORCID

Fabrice Besacier https://orcid.org/0000-0003-4990-8901

<figure><img src="image_6.png" /></figure>

Jonathan Farjon https://orcid.org/0000-0003-0707-2026

<figure><img src="image_7.png" /></figure>

## REFERENCES

1. Nouveaux produits de synthèse - Synthèse des connaissances - OFDT. https://www.ofdt.fr/produits-et-addictions/de-z/nouveaux-produits-de-synthèse/ Accessed July 16, 2021.

2. Actualités européennes - OFDT. https://www.ofdt.fr/europe-et-international/actualites-europeennes/ Accessed July 15, 2021.

3. Lukić V, Micić R, Arsić B, Nedović B, Radosavljević Ž. Overview of the major classes of new psychoactive substances, psychoactive effects, analytical determination and conformational analysis of selected ille- gal drugs. *Open Chemistry*. 2021;19(1):60-106. doi:10.1515/chem- 2021-0196

4. Elyashberg M. Identification and structure elucidation by NMR spectroscopy. *TrAC*. 2015;69:88-97. doi:10.1016/j.trac.2015.02.014

5. Elie MP, Elie LE, Baron MG. Keeping pace with NPS releases: fast GC-MS screening of legal high products. *Drug Test Anal*. 2013;5(5): 281-290. doi:10.1002/dta.1434

6. Angerer V, Bisel P, Moosmann B, Westphal F, Auwärter V. Separation and structural characterization of the new synthetic cannabinoid JWH-018 cyclohexyl methyl derivative “NE-CHMIMO” using flash chromatography, GC-MS, IR and NMR spectroscopy. *Forensic Sci Int*. 2016;266:93-98. doi:10.1016/j.forsciint.2016.05.031

7. Westphal F, Junge T, Girreser U, Greibl W, Doering C. Mass, NMR and IR spectroscopic characterization of pentedrone and pentylone and identification of their isocathinone by-products. *Forensic Sci Int*. 2012;217(1):157-167. doi:10.1016/j.forsciint.2011.10.045

8. Böselt L, Dötzer R, Steiner S, Stritzinger M, Salzmann S, Riniker S. Determining the regiochemistry and relative stereochemistry of small and druglike molecules using an alignment algorithm for infrared spectra. *Anal Chem*. 2020;92(13):9124-9131. doi:10.1021/acs. analchem.0c01399

9. Blümich B. Low-field and benchtop NMR. *J Magn Reson*. 2019;306: 27-35. doi:10.1016/j.jmr.2019.07.030

10. Castaing-Cordier T, Bouillaud D, Farjon J, Giraudeau P. Chapter Four - Recent advances in benchtop NMR spectroscopy and its applications. *Ann Rep NMR Spectros*; Webb, G. A., Ed.; Academic Press. 2021;103:191-258. PMID. doi:10.1016/bs.arnmr.2021.02.003

11. Tsochatzis ED, Giannopoulos G, Lopes JA, Guillou C. Drugs analysis and citizens safety. Workflows for chemical analysis and structural elucidation of new psychoactive substances in the EU. *Saf Sci*. 2021; 137:105126. doi:10.1016/j.ssci.2020.105126

12. Antonides LH, Brignall RM, Costello A, et al. Rapid identification of novel psychoactive and other controlled substances using low-Field ¹H NMR spectroscopy. ACS Omega. 2019;4(4):7103-7112. doi:10.1021/acsomega.9b00302

13. Zhong Y, Huang K, Luo Q et al. The application of a desktop NMR spectrometer in drug analysis. https://www.hindawi.com/journals/ ijac/2018/3104569/ Accessed January 31, 2022.

14. SWGDRUG Approved Recommendations. http://www.swgdrug.org/ approved.htm Accessed September 30, 2021.

15. Jones LE, Stewart A, Peters KL, et al. Infrared and raman screening of seized novel psychoactive substances: a large scale study of >200 Samples. *Analyst*. 2016;141(3):902-909. doi:10.1039/C5AN02326B

16. Naqi HA, Husbands SM, Blagbrough IS. ¹H quantitative NMR and UHPLC-MS analysis of seized MDMA/NPS mixtures and tablets from night-club venues. *Anal Methods*. 2019;11(37):4795-4807. doi:10.1039/C9AY01403A

17. Hussain JH, Gilbert N, Costello A, et al. Quantification of MDMA in seized tablets using benchtop ¹H NMR spectroscopy in the absence of internal standards. *Forensic Chem*. 2020;20:100263. doi:10.1016/j. forc.2020.100263

18. Tsochatzis E, Lopes JA, Reniero F, Holland M, Åberg J, Guillou C. Identification of 1-Butyl-Lysergic acid diethylamide (1B-LSD) in seized blotter paper using an integrated workflow of analytical techniques and chemo-informatics. *Molecules*. 2020;25(3):1-10. doi:10.3390/molecules25030712

19. Brandt SD, Kavanagh PV, Westphal F, et al. Synthetic cannabinoid receptor agonists: Analytical profiles and development of QMPSB, QMMSB, QMPCB, 2F-QMPSB, QMiPSB, and SGT-233. *Drug Test Anal*. 2021;13(1):175-196. doi:10.1002/dta.2913

20. Tsochatzis ED, Alberto Lopes J, Holland MV, Reniero F, Palmieri G, Guillou C. Identification and analytical characterization of a novel synthetic cannabinoid-type substance in herbal material in Europe. *Molecules*. 2021;26(4):1-16. doi:10.3390/molecules26040793

21. Malm L, Palm E, Souihi A, Plassmann M, Liigand J, Kruve A. Guide to semi-quantitative non-targeted screening using LC/ESI/HRMS. *Molecules*. 2021;26(12):3524. doi:10.3390/molecules26123524

## SUPPORTING INFORMATION

Additional supporting information can be found online in the Supporting Information section at the end of this article.

**How to cite this article:** Castaing-Cordier T, Benavides Restrepo A, Dubois D, et al. Characterization of new psychoactive substances by integrating benchtop NMR to multi-technique databases. *Drug Test Anal*. 2022;14(9): 1629-1638. doi:10.1002/dta.3332