Next Article in Journal
Assay Development and Identification of the First Plasmodium falciparum 7,8-dihydro-6-hydroxymethylpterin-pyrophosphokinase Inhibitors
Next Article in Special Issue
Determination of Total Sennosides and Sennosides A, B, and A1 in Senna Leaflets, Pods, and Tablets by Two-Dimensional qNMR
Previous Article in Journal
Application of MOS Gas Sensors Coupled with Chemometrics Methods to Predict the Amount of Sugar and Carbohydrates in Potatoes
Previous Article in Special Issue
A Novel qNMR Application for the Quantification of Vegetable Oils Used as Adulterants in Essential Oils
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

NMR Spectroscopy Applied to the Metabolic Analysis of Natural Extracts of Cannabis sativa

by
Maria Francesca Colella
1,
Rosachiara Antonia Salvino
1,
Martina Gaglianò
1,
Federica Litrenta
2,
Cesare Oliviero Rossi
1,
Adolfo Le Pera
3 and
Giuseppina De Luca
1,*
1
Department of Chemistry and Chemical Technologies (CTC), University of Calabria—UNICAL, Via P. Bucci 14C, 87036 Arcavacata di Rende, Italy
2
Department of Biomedical, Dental and Morphological and Functional Imaging Sciences (Biomorf), University of Messina, Polo Universitario dell’Annunziata, 98168 Messina, Italy
3
Calabra Maceri e Servizi s.p.a., Via M. Polo 54, 87036 Rende, Italy
*
Author to whom correspondence should be addressed.
Molecules 2022, 27(11), 3509; https://doi.org/10.3390/molecules27113509
Submission received: 22 April 2022 / Revised: 23 May 2022 / Accepted: 26 May 2022 / Published: 30 May 2022

Abstract

:
Cannabis sativa is a herbaceous multiple-use species commonly employed to produce fiber, oil, and medicine. It is now becoming popular for the high nutritional properties of its seed oil and for the pharmacological activity of its cannabinoid fraction in inflorescences. The present study aims to apply nuclear magnetic resonance (NMR) spectroscopy to provide useful qualitative and quantitative information on the chemical composition of seed and flower Cannabis extracts obtained by ultra-sound-assisted extraction, and to evaluate NMR as an alternative to the official procedure for the quantification of cannabinoids. The estimation of the optimal ω-6/ω-3 ratio from the 1H NMR spectrum for the seed extracts of the Futura 75 variety and the quantitative results from the 1H and 13C NMR spectra for the inflorescence extracts of the Tiborszallasi and Kompolti varieties demonstrate that NMR technology represents a good alternative to classical chromatography, supplying sufficiently precise, sensitive, rapid, and informative data without any sample pre-treatment. In addition, different extraction procedures were tested and evaluated to compare the elaboration of spectral data with the principal component analysis (PCA) statistical method and the quantitative NMR results: the extracts obtained with higher polarity solvents (acetone or ethanol) were poor in psychotropic agents (THC < LOD) but had an appreciable percentage of both cannabinoids and triacylgliceroles (TAGs). These bioactive-rich extracts could be used in the food and pharmaceutical industries, opening new pathways for the production of functional foods and supplements.

Graphical Abstract

1. Introduction

Cannabis sativa is a fast-growing annual dioecious weed, probably native to Central Asia and the Indian subcontinent [1,2], belonging to the Cannabaceae family (order Urticales) [3]. Despite its critical taxonomic definition, because of its complex chemical composition and the presence of several spontaneous generations of hybrid species, nowadays, classifying Cannabis as a highly polymorphic and hybridized monotypic genus (Cannabis sativa L.) is the most accepted definition. Cannabis is one of the oldest and most versatile sources of intoxicating resin, textile fiber, and phytocannabinoids, which are extracted from different parts of the plant, especially from the inflorescence and for seed oil. Hemp seed oil, obtained from Cannabis sativa L. seeds, is highly appreciated for its nutritional, anti-inflammatory, antioxidant, and immune-stimulating properties [4]. It is practically free of cannabinoids [5], so it has no psychoactive action but, like other common vegetable oils, it is rich in essential fatty acids [6]. As reported in several works, this oil is a rich source of ω-3 and ω-6 polyunsaturated fatty acids (almost 80%), in particular, linoleic acid (LA) and α-linolenic acid (αLA), with a ω-6/ω-3 ratio approximately equal to 3:1 [7]. Although various factors, such as cultivation area, cultivar, seed origin, agronomic cultivation practices, etc., affect both the chemical composition and the ω-6/ω-3 ratio [4,7], this ratio is considered an optimal nutritional value in the prevention of the risk of coronary heart disease [8,9]. Due to this characteristic, cannabis seed oils are authorized and widely used in the food sector [10], such as in the production of functional foods. Despite the growing interest in this product, specific regulations to evaluate its analytical quality parameters are still lacking [7]. In this context, it would be desirable to find methodologies that can provide useful and rapid information both on the chemical composition and on the important ω-6/ω-3 ratio.
The female inflorescences of the Cannabis plant have been widely used in the traditional medicine of different populations thanks to the pharmacological activity of some phytocannabinoids present in large quantities in these parts of the plant [11]. Phytocannabinoids are a class of terpenophenolic compounds with a 21-carbon backbone: 120 of these molecules naturally present in the plant have been identified and isolated to date [3,12,13]. These natural molecules are different from synthetic cannabinoids, generally used as therapeutic agents, and from endocannabinoids, which are endogenous lipid-based retrograde neurotransmitters capable of interacting with cannabinoid receptors in the human body [14]. The renewed and recent interest in cannabis is due to the identification of these molecules, whose different pharmacological activities such as anti-inflammatory action, cell growth inhibition, and tumor regression seem to be supported by numerous experimental evidence [15,16,17]. The chemical structures of the most common cannabinoids present in the cannabis plant are shown in Figure 1. Among them, the most representative are the well-known psychotropic agent Δ9-trans-tetrahydrocannabinol (Δ9-THC), cannabidiol (CBD), and its precursor cannabidiolic acid (CBDA). Compared to THC, CBD shows therapeutic benefits without euphoric or dysphoric effects, which is an advantage for clinical applications [15,16,17]. CBD has become very popular over the years for its health benefits, and nowadays is commercially available as a dietary supplement, a lotion, and most importantly, as a CBD oil. Indeed, the interest from the scientific community in the therapeutic potential of CBD oil is growing every day [18]. The reason is simple: it has already been used in various scientific studies for the treatment of numerous health problems, and is now recognized as one of the main elements of the so-called “therapeutic Cannabis” [19,20,21]. Industrial hemp crops, with a low THC content, have always been exploited as food and as a source of textile fibers, but they disappeared in the 1970s due to their association with the type of plants rich in THC [22]. The reintroduction of the cultivation of some hemp cultivars to produce fibers and seeds with a THC content lower than 0.2% w/w took place only several years later, i.e., in 2009, by means of an appropriate regulation published by the European Union [23]. Nowadays, in many countries, Cannabis sativa cultivations and medicines have been legalized under certain conditions due to their immense prospects in various medicinal applications [24,25]. The Italian legislation on C. sativa cultivation is somewhat ambiguous regarding the legal and illegal uses and cultivation of the plant, and differs based on the concentration of psychoactive cannabinoids. The law 242/2016 “Dispositions for the promotion of cultivation and supply chain of agro-industrial hemp” [26] is the most recent regulation in that direction, and is the reference text governing industrial hemp production in Italy for fiber or other industrial uses different from pharmaceutics, with cultivation based around certified seeds [27]. This measure establishes that the THC level must not exceed 0.2 %. However, even more recently, on 4 November 2019, the Italian Ministry of Health approved and ratified the “Definition of maximum levels of tetrahydrocannabinol (THC) in food” (GU n.11, 15-1-2020) [28]. This document fixes the content of THC at a maximum of 2 mg per kilo (0.0002%) in hemp seeds, flour, and derived foods and at a maximum of 5 mg per kilo (0.0005%) for the oil obtained from hemp seeds. It should be noted that the list of regulated foods provided in the appendix includes only seeds, flours and oil, but it seems that it will soon be updated based on new scientific evidence. Currently, the “Union method for the quantitative determination of the Δ9-tetrahydrocannabinol content in hemp varieties”, described by the annex III of the Commission Delegated Regulation (EU) 2017/1155 (last updated on 15 February 2017), is the only official procedure that member states must use for the quantitative determination of THC by gas chromatography (GC) after extraction with a suitable solvent [29]. It describes in detail everything about sampling, including sample dimensions, drying and storage, and techniques and reagents for the extraction and determination of THC, and it provides an allowed tolerance equal to 0.03% in absolute value. However, this official method is quite laborious, and the scientific community is always looking for advanced methodologies that will allow us to rapidly analyze natural mixtures without requiring manipulations or separations. An effective alternative to classical analytical methodologies could be the use of nuclear magnetic resonance (NMR) spectroscopy. NMR spectroscopy is a powerful and versatile technique that has progressively become a well-established tool in different areas of scientific research such as medicine, biology, and chemistry. The importance of NMR in the structural investigation of chemical compounds in liquid or solid phases is widely known, and its applicative power in addition to mass spectrometry (MS) has brought satisfactory results, particularly in the metabolic characterization of complex mixtures such as foods or natural extracts [30]. Despite it yielding relatively low-sensitivity measurements compared to MS (10 to 100 times better) with lower limits of detection (LOD), typically with an order of magnitude around micromolar [31], high-resolution NMR is becoming increasingly popular for fingerprinting as well as profiling. In particular, compared to MS, NMR spectroscopy is non-destructive, non-biased, non-invasive, does not damage analytes, and allows the use of samples such as tissues obtained, e.g., from biopsies, for further experiments [30]. Moreover, this technique is often fast and with low operating costs, it is easily quantifiable, and requires little or no chromatographic separation, sample treatment, or chemical derivatization. NMR is also a multinuclear technique that permits the routine and contemporaneous identification of a wide range of metabolites (such as sugars, organic acids, alcohols, polyols, and other highly polar compounds) in a highly quantitative and reproducible way thanks to 1D and 2D experiments. In addition, the combination of this high-throughput technique with chemometric methods is extremely advantageous because it gives the possibility to visualize, maximize, and therefore analyze the useful information contained in the experimental NMR data. The research presented in this work is placed in this scenario and aims to apply NMR methodologies to the study of natural extracts from the seeds and inflorescences of different cultivars of C. sativa with a THC/CBD << 1 ratio [32,33,34,35].
Specifically, the work aims to: (a) characterize the chemical profiles of the inflorescences and seeds for different varieties of C. sativa grown in Calabria (South of Italy) via NMR spectroscopy and, in particular, by using 1D (1H NMR, 13C NMR) and 2D (1H COSY, 1H-13C HMQC, 1H J-Res) experiments; (b) evaluate the extraction efficiency of different common solvents such as hexane, acetone and ethanol; (c) perform a multivariate statistical analysis (principal component analysis—PCA) based on 1D NMR data to discriminate samples coming from different extractive processes and/or varieties by identifying correlations between the metabolites that influence each metabolic profile; (d) and perform a quantification via NMR of the main cannabinoids (CBD, CBDA, and eventually Δ9-THC) using different internal standards.

2. Materials and Methods

2.1. Plant Material and Extraction Procedure

In this work, two different varieties of hemp were considered—Tiborszallasi and Kompolti—for both the metabolic characterization of inflorescences and for the quantitative and statistical analyses. NMR assignments were made on both varieties. Hemp inflorescences from Tiborszallasi and Kompolti grown in Calabria were harvested in September 2020, i.e., in the ripening period for both cases, selecting a reasonable number of plants for each cultivar located a few meters away in the same crop. All the collected hemp inflorescences were naturally air-dried, manually separated from twigs, and finely chopped. After this procedure, the samples were stored in the dark at 4.0 °C until analysis. Moreover, the chemical composition of the seed extracts was also investigated via NMR. The seed samples were collected from the Futura 75 cultivar by selecting, as for the inflorescences, an appropriate number of plants representative of the entire crop. The dry seeds were ground into a powder and stored in the dark at 4.0 °C until analysis.
All inflorescence and seed samples were provided by “Calabria Maceri e Servizi S.p.A.” (Rende, CS, Italy), while the crops were produced by the farm “Le Querce S.r.l” (Montalto Uffugo, CS, Italy).
The storage, pre-treatment sampling, and extraction procedures for the inflorescences of the Tiborszallasi and Kompolti varieties were mostly in accordance with the official “Union method for the quantitative determination of the Δ9-tetrahydrocannabinol content in hemp varieties” (Annex III of the Commission Delegated Regulation (EU) No. 639/2014, 11 March 2014) [29]. In order to evaluate the efficiency of the extraction and thus highlight any differences between the various extracts, different extraction solvents commonly available in chemical laboratories and with increasing polarity were chosen. The solvents used were n-hexane, acetone, and ethanol. For each sample, 1.0 g of dried, chopped, and stored inflorescence was extracted with 45 mL of solvent at room temperature for 20 min using an ultrasonic bath (30 kHz frequency). The obtained extracts were centrifuged for 5 min at 3000 rpm, the solutions were paper filtered, and the residues were extracted once more using the same procedure with another 45 mL of the same solvent. Lastly, the solvents were completely removed under vacuum. Starting from the same dried inflorescence matrix, 24 extractions were carried out for each variety. For the Tiborszallasi variety, 9 extractions were performed using ethanol and acetone and 6 using hexane, while for the Kompolti variety, 9 extractions were performed using hexane and ethanol and 6 using acetone. For the quantitative analysis, for both varieties, each extraction was carried out in triplicate to calculate an average value for the extraction yield and estimate the relative error. Then, three samples were collected for each solvent for a total of nine extracts for each variety. In addition, a quantitative analysis with gas chromatography (GC) using the flame ionization detector (FID) method was conducted on samples of the Tiborszallasi variety prepared from the same dried inflorescence matrix as the NMR samples by following the protocol reported in the literature [29]. It should be emphasized that, given the chemical complexity of C. sativa, the extraction and collection of its various bioactive compounds are not simple and, for this reason, both solvents and different extraction methods are reported in the literature, ranging from microwave-assisted extraction to supercritical fluid extraction [36].
The extraction procedure for seeds of Futura 75 were based on dynamic-maceration ultrasound-assisted extraction (UAE; Hielcher UP 100Hz, 100 W pulse, 30 kHz frequency), using ethanol as the solvent. Then, 2.00 g of seeds—dried, chopped and stored—were extracted with 45 mL of ethanol at room temperature for 20 min under magnetic stirring. The solution was then paper filtered, evaporated under vacuum at 30 °C, and the residue was extracted with the same procedure one more time with another 45 mL of same solvent [37]. The schematic experimental steps are shown in Figure 2.

2.2. Chemicals and Solvents

Pure solvents, ethanol (Absolute, ≥99.8%—VWR Chemicals, Briare, France), n-hexane (Laboratory Reagent, ≥95%—Sigma Aldrich, Darmstadt, Germany) and acetone, (ACS Reagent, ≥99.5%—Sigma Aldrich, Darmstadt, Germany) were used for the cannabinoid extractions. Deuterated chloroform (CDCl3—99.95 atom % D) as the solvent for NMR sample preparation, and anthracene, benzoic acid and 3-(trimethylsilyl) propionic-2,2,3,3-d4 acid (TMSP-d4—98 atom % D) as internal analytical standards for the quantitative NMR analysis were purchased from Sigma-Aldrich (Milan, Italy).

2.3. NMR Sample Preparation, Experimentation, and Data Processing

To prepare the NMR sample, after the evaporation under vacuum, 30.0 mg of seed extract of Futura 75 was dissolved in CDCl3 (~1 mL) directly in a 5 mm o.d. NMR tube. In this solution, the 1H NMR spectrum (spectral width (SW) of 14.00 ppm, 128 free induction decays (FIDs) and a relaxation delay of 5.0 s) and the 2D 1H COSY experiment (SW of 14.00 ppm on both dimensions, 2K data points, 40 scans, and 256 increments) were recorded. Two other similar extraction procedures were repeated on the same starting matrix of the dried seeds. From these extracts, 1D 1H NMR spectra were recorded to be used for reproducibility and standard deviation in the calculation of the essential fatty acids ratio. 1H NMR spectra were manually phased, baseline-corrected, and the chemical shifts were reported with respect to the TMS signal used as reference. From the 1H NMR spectra of these extracts, the main fatty acids ω-6/ω-3 ratio can be determined by combining the integrals, obtained after applying the deconvolution procedure, of three different signals: (a) the methyl protons of all the acyl groups (LA), with the exception of those of α-linolenic acid; (b) the methyl protons of ω-3 fatty acid (α-linolenic acid (αLA)); (c) the methylene protons of the linoleic and α-linolenic acyl groups; and using the relations [38]:
α LA = ( b ) ( b ) + ( a )
LA = 3 · ( c ) 4 · ( b ) 3 · [ ( b ) + ( a ) ]
The extract residues of the two inflorescence varieties were dissolved in 1.20 mL of CDCl3 and 600 µL of this solution was transferred into a 5 mm o.d. NMR tube. For the quantitative analysis, samples of hemp in CDCl3 were prepared by carefully weighing all the components and by adding 0.3 mg of internal standard (anthracene, benzoic acid and TMSP-d4). No additional treatment was necessary for the preparation of the NMR samples [39,40].
All the NMR experiments were performed on a Bruker Avance 500 MHz spectrometer (Bruker, Fällanden, Switzerland) working at a field strength of 11.74 T (500 MHz 1H Larmor frequency), equipped with a 5 mm multinuclear probe TBO (triple-resonance broadband observe) and a standard variable-temperature control unit BVT-3000 (Bruker, Fällanden, Switzerland). All isotropic spectra were recorded at room temperature with CDCl3 used as a field-frequency lock signal.
Spectral assignments of metabolites were based on the one-dimensional (1D) 1H, 13C, and 13C-{1H} NMR spectra, the bi-dimensional (2D) homo- and heteronuclear correlation NMR experiments (1H COSY, 1H-13C HMQC) and by comparison with the data reports in the literature [39,40,41]. In addition, Jij couplings between some pairs of protons were measured thanks to the homonuclear 2D experiment of 1H J-Resolved spectroscopy (J-Res) [42,43]. For each 1H NMR experiment, 128 FIDs were acquired using a spectral width of 14.00 ppm and a relaxation delay of 5.0 s. The 1D 13C-{1H} NMR spectra were recorded with proton broad-band decoupling, collecting 8K FIDs using a SW of 250.00 ppm and a relaxation delay of 5.0 s. For an accurate quantitative analysis of the metabolites present in the complex mixture, it was necessary to calibrate both the 90° pulses on the monitored nuclei (1H and 13C) and the T1 spin-lattice relaxation time. The T1 relaxation time was estimated by using the conventional inversion recovery experiment (10 increments from 0.5 ms to 30.0 s for 1H, and 24 increments from 0.1 ms to 300.0 s for 13C) [44]. The 1H quantitative NMR (qNMR) spectra were recorded using the same acquisition parameters described before but with a relaxation delay of 20.0 s. Instead, for 13C qNMR, quantification experiments (zgig Bruker pulse sequence) were performed, collecting 4000 FIDs using a SW of 250.00 ppm, a relaxation delay of 160.0 s, and an acquisition time of 10.0 s.
The initial relative quantification was obtained using Equation (3), in which the molar ratio M X M Y between the metabolites to quantify (X) and the internal standard (Y) is reported as a function of the ratio between their integral (IX and IY) and the ratio of resonant nuclei that generates the considered signal (NX and NY) [45,46].
M X M Y = I X I Y · N Y N X
Then, on the basis of the mass of extract used to prepare the NMR sample and the relative extraction yield, the absolute quantification was obtained in terms of the percentage of the dry weight of the hemp flowers.
The 1H COSY experiments were acquired using a SW of 14.00 ppm in both dimensions, 2K data points, 40 scans, and 256 increments; the 1H-13C HMQC spectra were recorded using SWs of 14.00 ppm (1H) and 250.00 ppm (13C), 2K datapoints, 512 scans, and 40 experiments, and 1H J-Res spectra were acquired with a SW of 12.00 ppm, 2K datapoints, 256 scans, and 48 experiments. A sine and a qsine filter were applied in both dimensions, F1 and F2, for the COSY and HMQC experiments, respectively, before being Fourier-transformed. Then, 1D NMR spectra were Fourier-transformed and manually phased, baseline-corrected, and aligned using the TMS signal as a reference. 13C-{1H} NMR spectra were filtered with 1.0 Hz line broadening before Fourier transformation. For the multivariate statistical analysis, the 1H NMR spectra were segmented in a rectangular bucket fixed at 0.05 ppm. The integration region was defined in the first proton spectra and next it was reported, once saved, in the other spectra bucketed automatically. In particular, variables were manually selected by choosing the regions of NMR spectra with characteristic signals of metabolites and eliminating regions with a poor signal-to-noise ratio. Therefore, regions selected for the subsequent statistical treatment were 0.50–1.15 ppm, 1.3–2.6 ppm, and 3.8–7.0 ppm. Once normalized, the integrals were organized in a data matrix that was mean-centered and scaled. All the data processing steps were carried out using TopSpin 3.6 software (Bruker BioSpin, Rheinstetten, Germany) (TopSpin, 2018) [47].

2.4. Statistical Analysis: Principal Component Analysis (PCA)

Multivariate analysis was performed on the 48 1D 1H NMR experiments recorded from the samples prepared from the extracts obtained with different extractive solvents. There were 24 samples for both varieties: Tiborszallasi and Kompolti. It should be noted that each sample prepared came from a different extract, i.e., for example, the 9 samples of Kompolti in hexane were obtained from 9 extractive procedures of the corresponding inflorescences. Among the possible statistical approaches, in this work, principal component analysis (PCA) was used as an explorative method. PCA is a technique able to reduce the dimensionality of a multivariate problem without losing information. This mathematical treatment allows us to clearly visualize samples in a two- or three-dimensional space and reveals trends in the data or groupings of samples (clusters) based on their similarity [48,49,50]. This methodology, applied to the data matrix of the bucketing 1H NMR spectra of all extract samples, obtained as described in the previous subsection, allowed us to obtain two datasets of 24 samples and 130 variables for each variety that were used as starting points to carry out PCA analysis using R software (R software (R Core Team, 2019)) (Vienna, Austria) [51].

2.5. Chromatographic Experiments

GC-FID analysis was intended to verify the qNMR results for the Tiborszallasi variety. It was carried out by using the protocol reported in [29] and a gas chromatograph (GC) equipped with a split/splitless injector and a flame ionization detector (FID) (Dani Master GC1000, Dani instrument, Milan, Italy).

3. Results and Discussion

3.1. NMR Characterization of Seeds Extracts

Figure 3 shows the 1D 1H NMR spectrum of a sample of hemp seed extracts prepared as described in Section 2.1.
From the proton spectrum, it is possible to recognize the classic profile of the esters of polyunsaturated fatty acids (PUFAs) [52,53] and, from the 2D 1H COSY spectrum, it is possible to trace all the correlations of the triacylglycerols (TAGs) containing both saturated and unsaturated fatty acids, as shown in Figure 4.
Let us consider, for example, the strong signal generated by the α-methylenic protons of all acyl chains (HA, Figure 4) detected between 2.27 ppm and 2.37 ppm. Clearly, the resonance frequency of the group of nuclei (labelled in Figure 4) for all fatty acids was not exactly the same, and it depended on the nature of the chain to which they belonged. Consequently, it was not possible to determine a well-resolved multiplicity to the final complex signals appearing in the spectrum due to the overlap of all these signals with slightly different chemical shifts. For this reason, it was indicated as a multiplet (m) and the corresponding range of chemical shift is reported in Table 1, in which the correlations obtained from the COSY spectrum for these protons are also reported. The same reasoning is valid for all the other signals reported as multiplets in Table 1, which summarizes all assignments for hemp seed extracts obtained from the 1H and 1H COSY spectra.
Since the integrals of the 1H NMR signals were proportional to the number of hydrogen atoms present in each functional group and, overall, to the number of functional groups present in the sample, from the combination of the integrals of different signals it is possible to calculate the concentration of fatty acids in general, and the ω-6/ω-3 ratio in particular. To this end, three different signals in the protonic spectra were considered: (a) the multiplet at 0.88 ppm due to the overlapping triplet signals of the methyl protons of all the acyl groups (LA), with the exception of those of α-linolenic acid; (b) the triplet at 0.97 ppm generated by the methyl protons of ω-3 fatty acid (α-linolenic acid; (αLA)); (c) the multiplet at 2.72–2.86 ppm generated by the diallylic protons of the linoleic and α-linolenic acyl groups. By combining the area of these signals, using the relations (1) and (2) that take into account the number of equivalent nuclei in each group, the concentrations of αLA and LA were calculated, from which the ω-6/ω-3 ratio was obtained [38].
Calculations were made considering three different samples for reproducibility and to give a main value and a standard deviation. The value obtained for the ω-6/ω-3 ratio using Equations (1) and (2) was 2.93 ± 0.07. As can be seen, this value obtained by 1H NMR is very close to the value of 3:1 for ω-6/ω-3 considered optimal for human dietary purposes, since it is able to prevent various diseases such as diabetes, cardiovascular disease, cancer, and other chronic diseases [4,9,53,54]. On the other hand, the growing interest in hemp seed oil in other fields such as pharmaceuticals and cosmetics [53] has resulted in a constant search for methods that allow fast and systematic quality control: the ω-6/ω-3 ratio is one of these quality parameters [55]. Our result is doubly important because, on the one hand, it indicates the quality of the oil extracts from the seeds of the Futura 75 cultivar grown in Calabria and, on the other hand, it confirms that NMR is a reliable quantitative platform for the fast screening of hemp oil quality. Indeed, the measurement of this ratio is based only on the recording and analysis of the 1H NMR spectra obtained directly from the seed extracts without further derivatization, as is required by the gas chromatography (GC) method, the common and validated method used to determine the composition of oil in terms of fatty acids [56]. Moreover, our result agrees with that reported in the paper of Siudem et al. [38] in which the authors analyzed six different samples of hemp seed oils and calculated the ω-6/ω-3 ratio by 1H NMR for each of them using the relationships (1) and (2). The authors compared these data with those from CG method applied to the same hemp seed oils and found a substantial agreement between them, which proves the effectiveness of the method. It is worth noting that, in a recent paper [53], the 1H NMR methodology—as well as being used in the evaluation of the ω-6/ω-3 ratio—is successfully combined with chemometric methods in order to observe the differences in several oil samples due to the different time and storage conditions of the oils. Hence, once again, this demonstrates how NMR can be used for the rapid and reliable analysis of hemp seed oil quality as an alternative to the more common classical analytical methods.

3.2. NMR Characterization of Flower Extracts

Figure 5 shows the 1H and 13C NMR spectra of inflorescence ethanolic extracts for the Tiborszallasi variety. As is evident from the figure, the two spectra exhibit a complex distribution of resonances due to strongly overlapped signals of cannabinoids with similar molecular structures (Figure 1), from which it is very difficult to recognize single metabolites through simple inspection.
To overcome these limitations that are characteristic of the 1D spectra of complex mixtures, the identification of the various cannabinoids present in the extracts (CBD, CBDA, CBG, THC) was carried out using the 2D correlation experiments 1H COSY, 1H-13C HMQC, and 1H J-Res and the data reported in the literature regarding the NMR characterizations of many single isolated cannabinoids [40,41,42]. It is worth underlining that, even in recent papers [35,36,57], the NMR methodology is mostly used to characterize fractions or single components of C. sativa obtained by separation techniques. In Figure 6, 1H COSY and 1H-13C HMQC spectra are shown together with the corresponding enlargement on the CBD and CBDA correlations taken as an example of the metabolic identification procedure adopted. In order to fully characterize the inflorescence extract, the 2D 1H J-Res experiments were performed on the same sample [42,43]. This type of experiment is capable of separating coupling and chemical shift information into two orthogonal dimensions, allowing for multiplet analysis. The one-dimensional spectrum (1H or 13C) is shown in the F2 dimension, while the couplings related to each chemical shift can be read in the indirect dimension, F1, at the chemical shift value of the multiplet signal to be studied. In Figure 7, the 2D 1H J-Res spectrum recorded on the sample of Tiborszallasi is reported, and, as an example, the case of H-9cis signal of CBD is reported in the enlargement. The signal of this proton in the 1D 1H NMR spectrum is a multiplet from which it is impossible to extract any useful information. On the contrary, the extrapolation of the column in correspondence to the chemical shift of H-9cis in the F2 dimension of the 1H J-Res spectrum gives the 1D profile of the considered signal in which it is possible to measure all the coupling constants. Indeed, as can be seen from the enlargement of Figure 7, this signal is a doublet of a quadruplet due to the interaction of the proton H-9cis with the proton H-9trans (doublet: 2J9cis-trans = 2.6 Hz) and with the methyl proton H-10 (quadruplet: 4J9cis-10 = 0.9 Hz).
The assignments for the various cannabinoids present in the ethanolic extract, CBD, CBDA, and CBG of the hemp inflorescences of the Tiborszallasi variety and their relative experimental information obtained from the NMR spectra (1D and 2D) are summarized in Table 2. These results agree with many other literature data [35,39,40,41] and demonstrate how the non-separative NMR technique, combining 1D and 2D experiments, can be successfully applied directly to a complex mixture, such as a Cannabis extract, for chemical characterization. An interesting result of the present study is the measurement of a large number of JH-H for the main cannabinoids, carried out again directly on the Cannabis extract through the analysis of the 2D J-Res spectra. It must be highlighted that, in the ethanolic extract, the signals related to THC were not detected in the recorded spectra; this means that, in this solvent, the quantity of THC extracted was below the NMR detection limit [31]. This also occurred for the acetone extracts while the spectra of the hexane extracts showed a small broad peak isolated at 6.40 ppm corresponding to the proton H-10 of Δ9-THC. This means that hexane, which is also the solvent indicated in the official method for determining the amount of Δ9-THC in hemp, was the most effective at extracting cannabinoids compared to the other solvents used. Once the assignment of cannabinoids in the Tiborszallasi variety was determined, a comparison between the NMR spectra of the samples relating to the Kompolti variety was possible. Except for the proton spectrum of the hexane extract of the Kompolti variety, in which Δ9-THC signals were absent, no further differences in terms of profile were distinguished in the spectra of any other samples of either variety. All these spectra are reported in the Supplementary Materials, where it is also possible to find the comparison between 1H NMR spectra from both varieties for the different solvents.
In addition to the cannabinoids, by comparison of the protonic spectra of the hemp seed oil and hemp inflorescence extracts, it was possible to recognize in this last spectrum some signals that referred to the triacylclycerols constituent. In particular, this profile was easily recognizable in the samples obtained by using ethanol as the extracting solvent. All the signal assignments are reported in Figure 8. Regarding the solvents, it must be underlined that many are used for cannabinoid extraction but, to date, there is no specific protocol. However, it is not surprising that we were able to better recognize the fatty acid profiles in the ethanol extracts since, from the literature data [58], it seems that ethanol has a greater extraction power than acetone and hexane. On the other hand, hexane, while showing the worst performance in terms of total yield, leads to cleaner extracts with fewer contaminants and extracts richer in cannabinoids [59]. For this reason, hexane is used for the selective analysis of cannabinoids in the official method established by the European Commission [29]. It should be noted that this was also evident in our spectra recorded on the hexane extracts. Indeed, these extracts have a cleaner profile than the other two solvents and allowed us, for the Tiborszallasi variety, to quantify the THC content together with the other cannabinoids.

3.3. Multivariate Analysis

Principal component analysis (PCA) was applied to discriminate hemp flowers samples based on their 1H NMR spectra [48,49,50]. The main aim was to observe how the cannabinoid profiles changed for extracts in relation to the nature of solvent extraction and the efficiency of the extraction procedure. PCA was carried out on two different data matrices, one for each hemp variety. For Tiborszallasi, the data matrix consisted of 24 samples (9 extracts for ethanol and acetone, respectively, 6 extracts for hexane) and 130 variables. For the Kompolti variety, the data matrix consisted of 24 sample (9 extracts for hexane and ethanol, respectively, 6 extracts for acetone) and 130 variables.
Figure 9 and Figure 10 show the 3D score plot for the Tiborszallasi and Kompolti variety samples, respectively, with a cumulative percentage of explained data variance for the three first PCs equal to: (a) 88.1 % (45.5% for PC1, 30.2% for PC2, 12.4% for PC3) for the Tiborszallasi variety and (b) 84.5% (51.8% for PC1, 17.5% for PC2, 15.2% for PC3) for the Kompolti variety. This means that data loss was negligible in both cases. The score plot effectively summarizes the relationship between the samples and highlights what was not discovered by the simple comparison among the protonic spectra: clear and well-defined separation between samples was observed and every type of extract was clustered into one defined region. Indeed, PC1 allows us to discriminate among extracts with non-polar solvent and extracts with higher polarity: hexane extracts had, in fact, positive values of PC1, while the other ones had negative values of PC1. Instead, PC2 seems to best discriminate between the acetone and ethanol extracts. The same result was obtained for both varieties.
However, in order to determine which metabolites have greater influence on the discrimination and evaluation of the quality and efficiency of the extraction procedure, an accurate discussion about loadings is necessary.
In Figure 11, a 2D biplot of the first two PCs (PC1 and PC2) for the Tiborszallasi variety is reported; this shows the sample differentiation and the changes in the metabolite concentrations from one extract to another and, consequently, the variables responsible for the sample clustering can be observed in the score diagram.
THC is the marker that gives the first clear distinction between hexane and all the other extracts: the THC loadings had positive values of PC1 and negative values of PC2, and these were very close to the scores of the hexane extracts. This means, as also demonstrated by the loadings of CBD and CBDA, that hexane is more efficient in the extraction of cannabinoids and of THC in particular. Indeed, the latter was below the NMR detection limit in the acetone and ethanol extracts. These results are not surprising since it is well known from the literature [34,60] that the polarity of the solvent affects the chemical composition of the cannabinoids present in the extracts.
Focusing on the TAGs previously detected during the signal assignment step, negative values of PC1 were obtained for these loadings. This indicates that extraction with a higher polarity solvent (acetone or ethanol) obtains samples with a lower percentage of cannabinoids but that are richer in fatty acids than hexane extracts. This peculiarity could be exploited to obtain extracts rich in both bioactive compounds—cannabinoids and TAGs—that could potentially be used in the food and pharmaceutical industries to produce functional foods and supplements.
The same conclusions were also reached for the Kompolti variety, whose 2D biplot of the first two PCs is reported in the Supplementary Materials.
The present study has shown that the exploratory PCA method, although simple, is able to differentiate samples coming from different extraction solvents and, in addition, it highlights which cannabinoids form the basis of this differentiation. Chemometrics-aided NMR techniques performed on C. sativa components (inflorescences, seed oils, leaves, and other less valuable parts of the plant), or on fractions of them, are widely used to determine different properties of C. sativa, as evidenced by the plethora of studies reported in the literature [34,41,61,62,63,64,65] and the references reported therein. In many of these works, the combination of 1H-NMR spectra with chemometric tools was used for the metabolomic differentiation of inflorescence extracts from different cultivars and to identify particular markers in order to discriminate different plant chemotypes [41,61]. In other studies [62,66], targeted and non-targeted NMR methodologies were used to identify and quantify compounds of different classes present in the inflorescence extracts of several C. sativa cultivars and to monitor their variations in three different harvested stages. The cultivars analyzed in these studies had a THC content always below the legal limit, while the quantities of the other cannabinoids in the extracts were affected by the harvest time and by the solvent. Although it is not easy to make comparisons given the quantity of factors affecting the composition of cannabinoids, the results obtained in this study substantially agree with the studies described above, and therefore seem to demonstrate the reliability of this method and its possible application in the routine analysis of cannabinoids.

3.4. Quantitative Analysis of Inflorescences

As a support to the results obtained from the PCA analysis and to evaluate if the NMR methodology could be a valid tool in the quantitative determination of the cannabinoids extracted as a function of the solvent used, an NMR quantification was also carried out using 1H and 13C NMR spectra and different internal standards.
Tiborszallasi- and Kompolti-type hemp flower samples used for the quantification were obtained, as mentioned before, by ultrasound-assisted extraction using three different common solvents of increasing polarity: hexane, acetone and ethanol. The complete procedure, from extraction to NMR measurement, was performed in triplicate to evaluate repeatability, to calculate an average value for the extraction yield, and to estimate the relative error. The procedure below is described for the Tiborszallasi samples, but exactly the same was performed for the Kompolti variety.
The average experimental extraction yields for each procedure were: (a) 17.9 ± 0.3% for extraction with hexane; (b) 9.9 ± 0.8% for extraction with acetone; (c) 19.6 ± 0.9% for extraction with ethanol.
For the quantification of the cannabinoids present in the extracts, both 1D 1H NMR and 13C NMR spectra were used. The quantification of cannabinoids by protonic spectra included the use of three different internal standards: anthracene, in accordance with literature [20], benzoic acid, and 3-(trimethylsilyl)propionic-2,2,3,3-d4 acid (TMSP-d4). These compounds are highly pure (≥ 99.9%), have low volatility, are chemically inert, and are not similar in structure to the cannabinoids to be quantified. Indeed, in the case of NMR spectroscopy, it is necessary that they generate well-isolated signals in the spectrum that do not overlap with the peaks assigned to other metabolites in the mixture. In this case, signals of aromatic protons of anthracene and benzoic acid ranged between 7.4 ppm and 8.5 ppm in the protonic spectrum, and the singlet of the methyl group of TMSP-d4 was set to 0.00 ppm.
The signals considered for the quantification, using Equation (3), were the following: for CBD, the H-9trans olefinic proton signal was set at δ = 4.65 ppm; for CBDA, the H-9cis olefinic proton signal was set at δ = 4.38 ppm; and for Δ9-THC, the H-10 proton signal was set at δ = 6.40 ppm. The results obtained for the three standards used are reported in Table 2. For further confirmation of the 1H qNMR results, quantitative 13C NMR spectra were also recorded on the hemp flower hexane extracts of Tiborszallasi, optimizing the acquisition parameters (see Section 2.3) and the operative conditions [40]. However, in this case, only TMSP-d4 was used as an internal standard for 13C qNMR because the aromatic carbons of anthracene and benzoic acid generated signals in the region between 125 ppm and 134 ppm, thereby overlapping with the cannabinoid signals. Moreover, given the low sensitivity of the 13C nucleus, the 13C qNMR was useful for the quantitative determination of CBDA, the only cannabinoid whose signals had an acceptable signal-to-noise ratio. The CBDA content (% on dry weight) obtained via 13C qNMR for Tiborszallasi in hexane extract was equal to 6.2 ± 0.9, which was substantially in agreement with that obtained from protonic spectra using the same internal standard (Table 2).
Observing the results reported in Table 2, the value obtained for the three cannabinoids CBD, CBD and Δ9-THC, for each type of extract using different internal standards, were fairly reproducible. This means that benzoic acid and TMSP-d4 seem to be valid alternatives to anthracene as internal standards for quantification. Indeed, anthracene is the standard that is usually used for the quantification of cannabinoids via NMR, but it has many drawbacks related to its toxicity, its poor solubility in common solvents, and long recording times due to long T1 [67]. Furthermore, the data show that hexane extracts are richer in cannabinoids than other solvents, as had already emerged in the biplot analysis, and, most interestingly, the quantities obtained for these extracts via NMR matched very well with the experimental data acquired following the official European procedure for the determination of cannabinoid Δ9-THC, as reported in Table 3. Indeed, the table shows the comparison, for the same sample of Tiborszallasi inflorescences, between the qNMR data of Δ9-THC via NMR and those obtained with the official GC-FID technique [29]. In the table, the GC-FID value concerning the quantity of CBD/CBDA obtained from the same sample is also reported. It should be noted that the GC-FID procedure does not allow discrimination between CBDA and CBD: before extraction, it includes a pretreatment of inflorescences at a high temperature that leads to the complete decarbossilation of the acidic form of this cannabinoid (CBDA) to the neutral (CBD). So, this unique value should be considered as the sum of the quantities of CBDA and CBD initially present in flowers. As can be seen from the table, also in this case, the data via NMR are in perfect agreement with those obtained via GC-FID. It is worth noting that the NMR technique has been proven to be a reliable and powerful tool for the quantification of different natural products and, especially in recent years, has also been successfully applied to the quantification of CBD and other cannabinoids directly on hemp extracts coming from different cultivars using both the 1H and 13C qNMR methodologies. [40,44,68,69,70,71]. Quantitative 1H-NMR is the most widely used method in the quantification of natural extracts and has been shown to have a good level of accuracy and reproducibility. However, its application to hemp extracts presents several drawbacks, due both to the presence of contaminants and to the overlapping of different signals that require the extensive use of deconvolution processes [72]. The use of 13C q-NMR, introduced by Marchetti et al. [40], partially removes these weaknesses and at the same time offers sufficiently precise and sensitive results. As expected, our results substantially agree with these previous reports on the 1H and 13C qNMR investigations, even if it is necessary to point out that a direct comparison on the quantities of cannabinoids found in the different extracts of the cultivar analyzed is practically impossible given, as we have already highlighted, the large quantity of variables, i.e., cultivar, geographical origin, harvesting period, agronomic practices, extraction methodologies, etc., that affect the composition of cannabinoids. However, the quantitative results of the present study highlight, once again, the remarkable potentialities of the NMR technique which was able to quantify the main metabolites present in the hemp inflorescence extracts we analyzed as they were, without the further treatment or derivatization required by the official technique [73,74]. Moreover, as reported recently by Dadiotis et al. [75] concerning the quantitative analysis of cannabinoids in hemp extracts using, in a complementary way, the 1H-NMR and 1H-1H COSY NMR spectra, these quantitative data via NMR are comparable with those acquired with other more consolidated techniques applied to the same extracts. This evidence shows good correspondence between the various quantification techniques, as also confirmed by the data we obtained on the Timborszallase cultivar given the satisfactory agreement between the NMR and GC-FID data of the hexane solvent.
In addition, the quantitative data obtained for the different solvents confirm the purely qualitative indications given by the PCA analysis, which proved to be very informative and fast.
These results were obtained for the Tiborszallasi variety but, as previously mentioned, the quantification of the main cannabinoids by 1H and 13C qNMR was also performed for the Kompolti variety and the results are reported in the Supplementary Materials.

4. Conclusions

Cannabis sativa is a fast-growing plant currently grown all over the world that is gaining popularity in various fields of research for its biological and pharmaceutical properties. Actually, C. sativa is widely recognized and appreciated for the high nutritional and health-promoting properties of the oil obtained from its seeds, together with the pharmacological activity mainly associated with psychoactive and non-psychoactive cannabinoids and the chemical components mainly extracted from the inflorescences. In this work, NMR spectroscopy was applied to analyze extracts from the seeds and inflorescences of different varieties of Cannabis sativa grown in Calabria in order to explore the potentialities of this technique for the qualitative and quantitative analysis of the extracts, and to evaluate the possibility of using it as an alternative to the most common methods in the quantification of cannabinoids present in inflorescence extracts. The quantitative NMR results obtained from two varieties of hemp inflorescence extracts, using different internal standards and solvents, demonstrated the high potentiality of the proposed technique in this field of application. Indeed, the NMR technique was able to quantify the main cannabinoids present in the extracts, the quantitative data were reproducible, and—most importantly—the data from the hexane solvent were congruent with the data obtained by the GC-FID method. Moreover, while this last methodology is not able to distinguish CBD and CBDA, using the NMR method, it was possible to separate the two contributions and quantify them. This proves, once again, the analytical power of the NMR technique which is not only able to offer the same results obtained from the official method, including the evaluation of THC, but can indicate more informative data without performing particular treatments on the sample.
In addition to the characterization and the quantitative study, different extraction procedures were tested and evaluated by NMR spectroscopy with the aim of obtaining inflorescence extracts poor in psychotropic agents and rich in medical cannabinoids and triacylglicerols (TAGs), which have an ω-6/ω-3 ratio that has been found to be excellent from a nutritional point of view. Specifically, extracts of inflorescences obtained by ultrasound-assisted solute–solvent extraction using hexane, acetone and ethanol as solvents were studied. By elaborating the spectral data with a statistical method (PCA) together with the qNMR approach, it was possible to conclude that hexane was more efficient in the extraction of cannabinoids (THC included) than the TAG constituents, while extraction with a higher polarity solvent (acetone or ethanol) obtained samples free from THC (THC content < LOD), rich in TAGs, and with a lower percentage of cannabinoids. This evidence can be exploited to obtain extracts rich in bioactive compounds (both cannabinoids and TAGs) that could potentially be used in the food and pharmaceutical industries, opening new paths for the production of functional foods and supplements.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules27113509/s1, Figure S1: Comparison between 1H NMR spectra ethanol (blue), acetone (red) and hexane (green) extracts for the Tiborszallasi variety; Figure S2: Comparison between 1H NMR spectra of ethanol (blue), hexane (red) and acetone (green) extracts for the Kompolti variety; Figure S3: Comparison between the enlarged region [6.25 ppm–6.5ppm] of the 1H NMR spectra from the hexane extract of the Tiborszallasi (blue) and Kompolti (purple) varieties. A broad peak isolated at 6.40 ppm corresponding to the proton H-10 of Δ9-THC appears in the proton spectra of Tiborszallasi, while this signal is undetectable in the 1H NMR spectrum acquired for Kompolti; Figure S4: Biplot of PCA carried out on NMR spectra of acetone (red dots), hexane (blue dots) and ethanol (orange dots) extracts of the Kompolti variety of hemp. The score plot shows the first two PCs (PC1 and PC2) with their respective variations. R2X(PC1) = 51.8%, R2X(PC2) = 17.5%; Table S1: The 1H NMR data of main cannabinoids in Kompolti inflorescences.

Author Contributions

Conceptualization, M.F.C. and G.D.L.; methodology, G.D.L.; validation, C.O.R., A.L.P.; formal analysis, M.G.; investigation, M.F.C., R.A.S., F.L.; data curation, M.F.C. and R.A.S.; writing—original draft preparation, M.F.C., R.A.S., G.D.L.; writing—review and editing, G.D.L. and C.O.R.; supervision, G.D.L. and C.O.R.; All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the University of Calabria and POR Calabria—FSE/FESR 2014–2020.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data generated or analyzed during this study are included in this published article.

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

Samples of the compounds are not available from the authors.

References

  1. Small, E. Classification of Cannabis sativa L. In Relation to Agricultural, Biotechnological, Medical and Recreational Utilization. In Cannabis sativa L.—Botany and Biotechnology, 1st ed.; Chandra, S., Hemant, L., ElSohly, M.A., Eds.; Springer International Publishing: Cham, Switzerland, 2017; pp. 1–62. [Google Scholar]
  2. Li, H.L. An archaeological and historical account of Cannabis in China. Econ. Bot. 1974, 28, 437–447. [Google Scholar] [CrossRef]
  3. ElSolhy, M.A.; Tadwan, M.M.; Gul, W.; Chandra, S.; Galal, A. Phytochemistry of Cannabis sativa L. Prog. Chem. Org. Nat. Prod. 2017, 103, 1–36. [Google Scholar]
  4. Crescente, G.; Piccolella, S.; Esposito, A.; Scognamiglio, M.; Fiorentino, A.; Pacifico, S. Chemical composition and nutraceutical properties of hemp seeds: An ancient food with actual functional value. Phytochem. Rev. 2018, 17, 733–749. [Google Scholar] [CrossRef]
  5. Tura, M.; Ansorena, D.; Astiasarán, I.; Mandrioli, M.; Toschi, T.G. Evaluation of Hemp Seed Oils Stability under Accelerated Storage Test. Antioxidants 2022, 11, 490. [Google Scholar] [CrossRef]
  6. Banskota, A.H.; Jones, A.; Hui, J.P.M.; Stefanova, R. Triacylglycerols and Other Lipids Profiling of Hemp By-Products. Molecules 2022, 27, 2339. [Google Scholar] [CrossRef]
  7. Spano, M.; Di Matteo, G.; Rapa, M.; Ciano, S.; Ingallina, C.; Cesa, S.; Menghini, L.; Carradori, S.; Giusti, A.M.; Di Sotto, A.; et al. Commercial Hemp Seed Oils: A Multimethodological Characterization. Appl. Sci. 2020, 10, 6933. [Google Scholar] [CrossRef]
  8. Teleszko, M.; Zając, A.; Rusak, T. Hemp Seeds of the Polish ‘Bialobrzeskie’ and ‘Henola’ Varieties (Cannabis sativa L. var. sativa) as Prospective Plant Sources for Food Production. Molecules 2022, 27, 1448. [Google Scholar] [CrossRef]
  9. Simopoulos, A.P. The importance of the omega-6/omega-3 fatty acid ratio in cardiovascular disease and other chronic disease. Exp. Biol. Med. 2008, 233, 674–688. [Google Scholar] [CrossRef]
  10. Ministero della Salute. Produzione e Commercializzazione di Prodotti a Base di Semi di Canapa Per L’utilizzo nci Settori Dell’alimentazione Umana; Ministero della Salute: Roma, Italy, 2009; pp. 1–4.
  11. Mercuri, A.M.; Accorsi, C.A.; Bandini Mazzanti, M. The long history of Cannabis and its cultivation by Romans in central Italy, shown by pollen records from Lago Albano and Lago di Nemi. Veget. Hist. Archaeobot. 2002, 11, 263–276. [Google Scholar] [CrossRef]
  12. ElSolhy, M.A.; Slade, D. Chemical Constituents of marijuana: The complex mixture of natural cannabinoids. Life Sci. 2005, 78, 539–548. [Google Scholar] [CrossRef]
  13. Radwan, M.M.; Chandra, S.; Gul, S.; Elsohly, M.A. Cannabinoids, phenolics, terpenes and alkaloids of cannabis. Molecules 2021, 26, 2774. [Google Scholar] [CrossRef] [PubMed]
  14. Bautista, J.L.; Yu, S.; Tian, L. Flavonoids in Cannabis sativa: Biosynthesis, bioactivities, and biotechnology. ACS Omega 2021, 6, 5119–5123. [Google Scholar] [CrossRef] [PubMed]
  15. Baker, D.; Pryce, G.; Giovannoni, G.; Thompson, A.J. The therapeutic potential of Cannabis. Lancet. Neurol. 2003, 2, 191–298. [Google Scholar] [CrossRef]
  16. Izzo, A.A.; Borrelli, F.; Capasso, R.; Di Marzo, V.; Mechoulam, R. Non-psychotropic plant cannabinoids: New therapeutic opportunities from an ancient herb. Trends Pharmacol. Sci. 2009, 30, 515–527. [Google Scholar] [CrossRef] [PubMed]
  17. Odieka, A.E.; Obuzor, G.U.; Oyedeji, O.O.; Gondwe, M.; Hosu, Y.S.; Oyedeji, A.O. The Medicinal Natural Products of Cannabis sativa Linn.: A Review. Molecules 2022, 27, 1689. [Google Scholar] [CrossRef]
  18. Iftikhar, A.; Zafar, U.; Ahmed, W.; Shabbir, M.A.; Sameen, A.; Sahar, A.; Bhat, Z.F.; Kowalczewski, P.Ł.; Jarzębski, M.; Aadil, R.M. Applications of Cannabis sativa L. in Food and Its Therapeutic Potential: From a Prohibited Drug to a Nutritional Supplement. Molecules 2021, 26, 7699. [Google Scholar] [CrossRef]
  19. Russo, E.B. Cannabis Therapeutics and the Future of Neurology. Front. Integr. Neurosci. 2018, 12, 51. [Google Scholar] [CrossRef] [PubMed]
  20. Stasiłowicz, A.; Tomala, A.; Podolak, I.; Cielecka-Piontek, J. Cannabis sativa L. as a Natural Drug Meeting the Criteria of a Multitarget Approach to Treatment. Int. J. Mol. Sci. 2021, 22, 778. [Google Scholar] [CrossRef]
  21. Kopustinskiene, D.M.; Masteikova, R.; Lazauskas, R.; Bernatoniene, J. Cannabis sativa L. Bioactive Compounds and Their Protective Role in Oxidative Stress and Inflammation. Antioxidants 2022, 11, 660. [Google Scholar] [CrossRef]
  22. Pisanti, S.; Bifulco, M. Medical Cannabis: A plurimillennial history of an evergreen. J. Cell. Physiol. 2019, 234, 8342–8351. [Google Scholar] [CrossRef]
  23. European Commission. Regulation (EC) No 1107/2009 of the European Parliament and of the Council of 21 October 2009 concerning the placing of plant protection products on the market and repealing Council Directives 79/117/EEC and 91/414/EEC. Off. J. Eur. Union 2009, 309, 1–50. Available online: https://eur-lex.europa.eu/legal-content/EN/TXT/?uri=celex%3A32009R1107 (accessed on 14 April 2022).
  24. Schachtsiek, J.; Warzecha, H.; Kayser, O.; Stehle, F. Current Perspectives on Biotechnological Cannabinoid Production in Plants. Planta Med. 2018, 84, 214–220. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Brunetti, P.; Pichini, S.; Pacifici, R.; Busardò, F.P.; del Rio, A. Herbal preparations of medical cannabis: A vademecum for prescribing doctors. Medicina 2020, 56, 237. [Google Scholar] [CrossRef]
  26. Gazzetta Ufficiale. Legge 2 Dicembre 2016, n.242. Disposizioni per la Promozione della Coltivazione e della Filiera Agroindustriale della Canapa (16G00258), GU Serie Generale n. 304 del 30-12-2016. Available online: https://www.gazzettaufficiale.it/eli/id/2016/12/30/16G00258/sg (accessed on 20 April 2022).
  27. Circolare—Camera dei Deputati. Circolare 31 Luglio 2018, Prot. 2018/43586. Aspetti Giuridico-Operativi Connessi al Fenomeno della Commercializzazione delle Infiorescenze della Canapa Tessile a Basso Tenore di THC e Relazioni con la Normativa Sugli Stupefacenti. Available online: https://www.camera.it/temiap/2018/12/14/OCD177-3851.pdf (accessed on 2 April 2022).
  28. Gazzetta Ufficiale. Decreto 4 Novembre 2019 Definizione di Livelli Massimi di Tetraidrocannabinolo (THC) Negli Alimenti. (20A00016) (GU Serie Generale n.11 del 15-01-2020). Available online: https://www.gazzettaufficiale.it/eli/id/2020/01/15/20A00016/sg (accessed on 20 April 2022).
  29. European Union. Annex III, Commission Delegated Regulation (EU) 2017/1155 of 15 February 2017 Amending Delegated Regulation (EU) No 639/2014 as Regards the Control Measures Relating to the Cultivation of Hemp, Certain Provisions on the Greening Payment, the Payment for Young Farmers in Control of a Legal Person, the Calculation of the per Unit Amount in the Framework of Voluntary Coupled Support, the Fractions of Payment Entitlements and Certain Notification Requirements Relating to the Single Area Payment Scheme and the Voluntary Coupled Support, and Amending Annex X to Regulation (EU) No 1307/2013 of the European Parliament and of the Council. Available online: https://eur-lex.europa.eu/eli/reg_del/2017/1155/oj (accessed on 2 April 2022).
  30. Vignoli, A.; Ghini, V.; Meoni, G.; Licari, C.; Takis, P.G.; Tenori, L.; Turano, P.; Luchinat, C. High-Throughput Metabolomics by 1D NMR. Angew. Chem. Int. 2019, 58, 968–994. [Google Scholar] [CrossRef] [PubMed]
  31. Eisenmann, P.; Ehlers, M.; Weinert, C.H.; Tzvetkova, P.; Silber, M.; Rist, M.J.; Luy, B.; Muhle-Goll, C. Untargeted NMR Spectroscopic Analysis of the Metabolic Variety of New Apple Cultivars. Metabolites 2016, 6, 29. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Chandra, S.; Lata, H.; Khan, I.A.; ElSohly, M.A. Cannabis sativa L.—Botany and Horticulture. In Cannabis sativa L.—Botany and Biotechnology, 1st ed.; Chandra, S., Hemant, L., ElSohly, M.A., Eds.; Springer International Publishing: Cham, Switzerland, 2017; pp. 79–100. [Google Scholar]
  33. Aizpurua-Olaizola, O.; Soydaner, U.; Öztürk, E.; Schibano, D.; Simsir, Y.; Navarro, P.; Etxebarria, N.; Usobiaga, A. Evolution of the Cannabinoid and Terpene Content during the Growth of Cannabis sativa Plants from Different Chemotypes. J. Nat. Prod. 2016, 79, 324–331. [Google Scholar] [CrossRef]
  34. Brighenti, V.; Protti, M.; Anceschi, L.; Zanardi, C.; Mercolini, L.; Pellati, F. Emerging challenges in the extraction, analysis and bioanalysis of cannabidiol and related compounds. J. Pharm. Biomed. Anal. 2021, 192, 113633. [Google Scholar] [CrossRef]
  35. Ohtsuki, T.; Friesen, J.B.; Chen, S.N.; McAlpine, J.B.; Pauli, G.F. Selective Preparation and High Dynamic-Range Analysis of Cannabinoids in “CBD Oil” and Other Cannabis sativa Preparations. J. Nat. Prod. 2022, 85, 634–646. [Google Scholar] [CrossRef]
  36. Liu, Y.; Liu, H.-Y.; Li, S.-H.; Ma, W.; Wu, D.-T.; Li, H.-B.; Xiao, A.-P.; Liu, L.-L.; Zhu, F.; Gan, R.-Y. Cannabis sativa bioactive compounds and their extraction, separation, purification, and identification technologies: An updated review. TrAC 2022, 149, 116554. [Google Scholar] [CrossRef]
  37. Rezvankhah, A.; Emam-Djomeh, Z.; Safari, M.; Askari, G.; Salami, M. Investigation on the extraction yield, quality, and thermal properties of hemp seed oil during ultrasound-assisted extraction: A comparative study. J. Food Process. Preserv. 2018, 42, e13766. [Google Scholar] [CrossRef]
  38. Siudem, P.; Wawer, I.; Paradowska, K. Rapid evaluation of edible hemp oil quality using NMR and FT-IR spectroscopy. J. Mol. Struct. 2019, 1177, 204–208. [Google Scholar] [CrossRef]
  39. Hazekamp, A.; Choi, Y.H.; Verpoorte, R. Quantitative Analysis of Cannabinoids from Cannabis sativa using 1H-NMR. Chem. Pharm. Bull. 2004, 52, 718–721. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Marchetti, L.; Brighenti, V.; Rossi, M.C.; Sperlea, J.; Pellati, F.; Bertelli, D. Use of 13C-qNMR Spectroscopy for the Analysis of Non-Psychoactive Cannabinoids in Fibre-Type Cannabis sativa L. (Hemp). Molecules 2019, 24, 1138. [Google Scholar] [CrossRef] [Green Version]
  41. Choi, Y.H.; Kim, H.K.; Hazekamp, A.; Erkelens, C.; Lefeber, A.W.M.; Verpoorte, R. Metabolomic differentiation of Cannabis sativa cultivars using 1H NMR spectroscopy and principal component analysis. J. Nat. Prod. 2004, 67, 953–957. [Google Scholar] [CrossRef] [PubMed]
  42. Ludwig, C.; Viant, M.R. Two-dimensional J-resolved NMR spectroscopy: Review of a key methodology in the metabolomics toolbox. Phytochem. Anal. 2010, 21, 22–32. [Google Scholar] [CrossRef] [PubMed]
  43. Huang, Y.; Zhang, Z.; Chen, H.; Feng, J.; Cai, S.; Chen, Z. A high-resolution 2D J-resolved NMR detection technique for metabolite analyses of biological samples. Sci. Rep. 2015, 5, 8390. [Google Scholar] [CrossRef] [Green Version]
  44. Brighenti, V.; Marchetti, L.; Anceschi, L.; Protti, M.; Verri, P.; Pollastro, F.; Mercolini, L.; Bertelli, D.; Zanardi, C.; Pellati, F. Separation and non-separation methods for the analysis of cannabinoids in Cannabis sativa L. J. Pharm. Biomed. 2021, 206, 114346. [Google Scholar] [CrossRef]
  45. Kumar Bharti, S.; Roy, R. Quantitative 1H NMR Spectroscopy. TrAC 2012, 35, 5–26. [Google Scholar]
  46. Araneda, J.F.; Chu, T.; Leclerc, M.C.; Riegel, S.D.; Spingarn, N. Quantitative analysis of cannabinoids using benchtop NMR instruments. Anal. Methods 2020, 12, 4853–4857. [Google Scholar] [CrossRef]
  47. TopSpin. 2018. Available online: http://www.bruker-biospin.com/topspin.html (accessed on 15 March 2022).
  48. Varmuz, K.; Filzmoser, P. Introduction to Multivariate Statistical Analysis in Chemometrics, 1st ed.; CRC Press Taylor & Francis Group: Boca Raton, FL, USA, 2009; pp. 59–101. [Google Scholar]
  49. Ebrahimi, P.; Viereck, N.; Bro, R.; Engelsen, S.B. Chemometric Analysis of NMR Spectra. In Modern Magnetic Resonance, 2nd ed.; Webb, A.G., Ed.; Springer International Publishing: Cham, Switzerland, 2018; pp. 1649–1668. [Google Scholar]
  50. Ren, S.; Hinzman, A.A.; Kang, E.L.; Szczesniak, R.; Lu, L.J. Computational and Statistical Analysis of Metabolomics Data. Metabolomics 2015, 11, 1492–1513. [Google Scholar] [CrossRef]
  51. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2019; Available online: https://www.R-project.org/ (accessed on 21 March 2022).
  52. Popescu, R.; Costinel, D.; Dinca, O.R.; Marinescu, A.; Stefanescu, I.; Ionete, R.E. Discrimination of vegetable oils using NMR spectroscopy and chemometrics. Food Control 2015, 48, 84–90. [Google Scholar] [CrossRef]
  53. Siudem, P.; Wawer, I.; Kowalska, V.; Paradowska, K. Rapid 1H NMR and chemometric methods in verification of hemp-seed oil quality. J. Pharm. Biomed. 2022, 212, 114650. [Google Scholar] [CrossRef] [PubMed]
  54. Mikulcová, V.; Kašpárková, V.; Humpolíček, P.; Bunková, L. Formulation, Characterization and Properties of Hemp Seed Oil and Its Emulsions. Molecules 2017, 22, 700. [Google Scholar] [CrossRef] [PubMed]
  55. Farinon, B.; Costantini, L.; Molinari, R.; Di Matteo, G.; Garzoli, S.; Ferri, S.; Ceccantoni, B.; Mannina, L.; Merendino, N. Effect of malting on nutritional and antioxidant properties of the seeds of two industrial hemp (Cannabis sativa L.) cultivars. Food Chem. 2022, 370, 131348. [Google Scholar] [CrossRef]
  56. Borges, G.R.; Birk, L.; Scheid, C.; Morés, L.; Carasek, E.; Kitamura, R.O.S.; Roveri, F.; Eller, S.; de Oliveira Merib, J.; de Oliveria, T.F. Simple and straightforward analysis of cannabinoids in medicinal products by fast-GC–FID. Forensic Toxicol. 2020, 38, 531–535. [Google Scholar] [CrossRef]
  57. Mazzara, E.; Torresi, J.; Fico, G.; Papini, A.; Kulbaka, N.; Dall’Acqua, S.; Sut, S.; Garzoli, S.; Mustafa, A.M.; Cappellacci, L.; et al. A Comprehensive Phytochemical Analysis of Terpenes, Polyphenols and Cannabinoids, and Micromorphological Characterization of 9 Commercial Varieties of Cannabis sativa L. Plants 2022, 11, 891. [Google Scholar] [CrossRef]
  58. Brighenti, V.; Pellati, F.; Steinbach, M.; Maran, D.; Benvenuti, S. Development of a new extraction technique and HPLC method for the analysis of non-psychoactive cannabinoids in Cannabis sativa L. fibosa (hemp). J. Pharm. Biomed. Anal. 2017, 143, 228–236. [Google Scholar] [CrossRef]
  59. Pellati, F.; Brighenti, V.; Sperlea, J.; Marchetti, L.; Bertelli, D.; Benvenuti, S. New methods for the comprehensive analysis of bioactive compounds in Cannabis sativa L. (hemp). Molecules 2018, 23, 2639. [Google Scholar] [CrossRef] [Green Version]
  60. Valizadehderakhshan, M.; Shahbazi, A.; Kazem-Rostami, M.; Todd, M.S.; Bhowmik, A.; Wang, L. Extraction of Cannabinoids from Cannabis sativa L. (Hemp)—Review. Agriculture 2021, 11, 384. [Google Scholar] [CrossRef]
  61. Peschel, W.; Politi, M. 1H NMR and HPLC/DAD for Cannabis sativa L. Chemotype distinction, extract profiling and specification. Talanta 2015, 140, 150–165. [Google Scholar] [CrossRef]
  62. Spano, M.; Di Matteo, G.; Ingallina, C.; Botta, B.; Quaglio, D.; Ghirga, F.; Balducci, S.; Cammarone, S.; Campiglia, E.; Giusti, A.M.; et al. Emerging challenges in the extraction, analysis and bioanalysis of cannabidiol and related compounds. Molecules 2021, 26, 2912. [Google Scholar] [CrossRef] [PubMed]
  63. Monti, M.C.; Frei, P.; Weber, S.; Scheurer, E.; · Mercer-Chalmers-Bender, K. Beyond Δ9-tetrahydrocannabinol and cannabidiol: Chemical differentiation of cannabis varieties applying targeted and untargeted analysis. Anal. Bioanal. Chem. 2022, 414, 3847–3862. [Google Scholar] [CrossRef] [PubMed]
  64. Palmieri, S.; Mascini, M.; Ricci, A.; Fanti, F.; Ottaviani, C.; Lo Sterzo, C.; Sergi, M. Identification of Cannabis sativa L. (hemp) Retailers by Means of Multivariate Analysis of Cannabinoids. Molecules 2019, 24, 3602. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Mudge, E.M.; Murch, S.J.; Brown, P.N. Chemometric Analysis of Cannabinoids: Chemotaxonomy and Domestication Syndrome. Sci. Rep. 2018, 8, 13090. [Google Scholar] [CrossRef] [PubMed]
  66. Ingallina, C.; Sobolev, A.P.; Circi, S.; Spano, M.; Fraschetti, C.; Filippi, A.; Di Sotto, A.; Di Giacomo, S.; Mazzoccanti, G.; Gasparrini, F.; et al. Cannabis sativa L. Inflorescences from Monoecious Cultivars Grown in Central Italy: An Untargeted Chemical Characterization from Early Flowering to Ripening. Molecules 2020, 25, 1908. [Google Scholar] [CrossRef] [PubMed]
  67. Olasehinde, T.A.; Olaniran, A.O. Neurotoxicity of anthracene and benz[a]anthracene involves oxidative stress-induced neuronal damage, cholinergic dysfunction and disruption of monoaminergic and purinergic enzymes. Toxicol. Res. 2022, 80, 105312. [Google Scholar] [CrossRef]
  68. Nagy, D.U.; Cianfaglione, K.; Maggi, F.; Sut, S.; Dall’Acqua, S. Chemical Characterization of Leaves, Male and Female Flowers from Spontaneous Cannabis (Cannabis sativa L.) Growing in Hungary. Chem. Biodivers. 2019, 16, e1800562. [Google Scholar] [CrossRef]
  69. Kostas Ioannidis, K.; Dadiotis, E.; Mitsis, V.; Melliou, E.; Magiatis, P. Biotechnological Approaches on Two High CBD and CBG Cannabis sativa L. (Cannabaceae) Varieties: In Vitro Regeneration and Phytochemical Consistency Evaluation of Micropropagated Plants Using Quantitative 1H-NMR. Molecules 2020, 25, 5928. [Google Scholar] [CrossRef]
  70. Siciliano, C.; Bartella, L.; Mazzotti, F.; Aiello, D.; Napoli, A.; De Luca, P.; Temperini, A. 1H NMR quantification of cannabidiol (CBD) in industrial products derived from Cannabis sativa L. (hemp) seeds. IOP Conf. Ser. Mater. Sci. Eng. 2019, 572, 012010. [Google Scholar] [CrossRef]
  71. Barthlott, I.; Scharinger, A.; Golombek, P.; Kuballa, T.; Lachenmeier, D.W. A Quantitative 1H NMR Method for Screening Cannabinoids in CBD Oils. Toxics 2021, 9, 136. [Google Scholar] [CrossRef]
  72. Risoluti, R.; Gullifa, G.; Battistini, A.; Materazzi, S. Monitoring of cannabinoids in hemp flours by MicroNIR/Chemometrics. Talanta 2020, 211, 120672. [Google Scholar] [CrossRef] [PubMed]
  73. Sgrò, S.; Lavezzi, B.; Caprari, C.; Polito, M.; D’Elia, M.; Lago, G.; Furlan, G.; Girotti, S.; Ferri, E.N. Delta9-THC determination by the EU official method: Evaluation of measurement uncertainty and compliance assessment of hemp samples. Anal. Bioanal. Chem. 2021, 413, 3399–3410. [Google Scholar] [CrossRef] [PubMed]
  74. Nahar, L.; Guo, M.; Sarker, S.D. Gas chromatographic analysis of naturally occurring cannabinoids: A review of literature published during the past decade. Phytochem. Anal. 2020, 31, 135–146. [Google Scholar] [CrossRef] [PubMed]
  75. Dadiotis, E.; Mitsis, V.; Melliou, E.; Magiatis, P. Direct Quantitation of Phytocannabinoids by One-Dimensional 1H qNMR and Two-Dimensional 1H-1H COSY qNMR in Complex Natural Mixtures. Molecules 2022, 27, 2965. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Chemical structure and nuclei numbering of molecular fragments in hemp principal cannabinoids.
Figure 1. Chemical structure and nuclei numbering of molecular fragments in hemp principal cannabinoids.
Molecules 27 03509 g001
Figure 2. Schematic experimental steps involved in the extraction, sample preparation, and NMR characterization of C. sativa seeds and inflorescences.
Figure 2. Schematic experimental steps involved in the extraction, sample preparation, and NMR characterization of C. sativa seeds and inflorescences.
Molecules 27 03509 g002
Figure 3. 1H NMR spectrum (500 MHz) of hemp seed oil dissolved in CDCl3 recorder at 298K, obtained with ultrasound-assisted extraction (UAE) procedure.
Figure 3. 1H NMR spectrum (500 MHz) of hemp seed oil dissolved in CDCl3 recorder at 298K, obtained with ultrasound-assisted extraction (UAE) procedure.
Molecules 27 03509 g003
Figure 4. (a) General structure and nuclei labelling of molecular fragments in triacylglicerols (TAGs); (b) 2D map 1H COSY of the ethanolic extract with the UAE of C. sativa seeds. In the figure, all cross peaks corresponding to the homuncular correlations have been highlighted.
Figure 4. (a) General structure and nuclei labelling of molecular fragments in triacylglicerols (TAGs); (b) 2D map 1H COSY of the ethanolic extract with the UAE of C. sativa seeds. In the figure, all cross peaks corresponding to the homuncular correlations have been highlighted.
Molecules 27 03509 g004
Figure 5. NMR spectra (500 MHz) of inflorescence extracts dissolved in CDCl3 recorded at 298K. (a) 1H NMR spectrum recorded using zg30 Bruker standard pulses sequence; for each experiment, 128 FIDs were accumulated using a spectral width of 14.00 ppm and a relaxation delay of 5 s. (b) 13C-{1H} NMR spectrum (zgig Bruker pulse sequence) performed with proton broad-band decoupling, collecting 8K free induction decays (FIDs) and using a spectral width of 250.00 ppm and a relaxation delay of 5s. One-dimensional NMR FIDs were Fourier-transformed, phased, baseline-corrected, and aligned using the TMS signal as a reference. 13C-{1H} NMR spectra were filtered with 1 Hz line broadening before Fourier transformation.
Figure 5. NMR spectra (500 MHz) of inflorescence extracts dissolved in CDCl3 recorded at 298K. (a) 1H NMR spectrum recorded using zg30 Bruker standard pulses sequence; for each experiment, 128 FIDs were accumulated using a spectral width of 14.00 ppm and a relaxation delay of 5 s. (b) 13C-{1H} NMR spectrum (zgig Bruker pulse sequence) performed with proton broad-band decoupling, collecting 8K free induction decays (FIDs) and using a spectral width of 250.00 ppm and a relaxation delay of 5s. One-dimensional NMR FIDs were Fourier-transformed, phased, baseline-corrected, and aligned using the TMS signal as a reference. 13C-{1H} NMR spectra were filtered with 1 Hz line broadening before Fourier transformation.
Molecules 27 03509 g005
Figure 6. (a) Structure of CBD and CBDA reported together with atom numbering adopted; (b) 500 MHz 1H COSY spectrum; and (c) 1H-13C HMQC NMR spectrum of inflorescence ethanolic extract sample of Tiborszallasi, a variety dissolved in CDCl3 and recorded at 298K.
Figure 6. (a) Structure of CBD and CBDA reported together with atom numbering adopted; (b) 500 MHz 1H COSY spectrum; and (c) 1H-13C HMQC NMR spectrum of inflorescence ethanolic extract sample of Tiborszallasi, a variety dissolved in CDCl3 and recorded at 298K.
Molecules 27 03509 g006aMolecules 27 03509 g006b
Figure 7. 1H J-Res NMR spectrum of inflorescence ethanolic extract samples of Tiborszallasi dissolved in CDCl3 and recorded at 298K. The projection of the H-9cis signal (doublet of quadruplets) of CBD is reported in the enlargement.
Figure 7. 1H J-Res NMR spectrum of inflorescence ethanolic extract samples of Tiborszallasi dissolved in CDCl3 and recorded at 298K. The projection of the H-9cis signal (doublet of quadruplets) of CBD is reported in the enlargement.
Molecules 27 03509 g007
Figure 8. 1H NMR spectrum of an ethanolic extract of C. sativa inflorescences (Tiborszallasi variety). For the enlarged regions of (a) [0.8 ppm–2.2 ppm], (b) [2.2 ppm–3.00 ppm], (c) [3.7 ppm–4.8 ppm], and (d) [4.8 ppm–6.5 ppm] the signal attributions are shown. Signal assignments are reported with different colors: orange for TAGs; grey for the glycerol backbone; blue for CBD/CBDA; green for CBG/CBGA.
Figure 8. 1H NMR spectrum of an ethanolic extract of C. sativa inflorescences (Tiborszallasi variety). For the enlarged regions of (a) [0.8 ppm–2.2 ppm], (b) [2.2 ppm–3.00 ppm], (c) [3.7 ppm–4.8 ppm], and (d) [4.8 ppm–6.5 ppm] the signal attributions are shown. Signal assignments are reported with different colors: orange for TAGs; grey for the glycerol backbone; blue for CBD/CBDA; green for CBG/CBGA.
Molecules 27 03509 g008
Figure 9. Principal component analysis (PCA) of hexane (brown dots), acetone (purple dots) and ethanol (green dots) extracts of the Tiborszallasi variety of hemp. The score plot shows the first three PCs (PC1, PC2 and PC3) with their respective variations. R2X(PC1) = 45.5%, R2X(PC2) = 30.2%, R2X(PC3) = 12.4%.
Figure 9. Principal component analysis (PCA) of hexane (brown dots), acetone (purple dots) and ethanol (green dots) extracts of the Tiborszallasi variety of hemp. The score plot shows the first three PCs (PC1, PC2 and PC3) with their respective variations. R2X(PC1) = 45.5%, R2X(PC2) = 30.2%, R2X(PC3) = 12.4%.
Molecules 27 03509 g009
Figure 10. Principal component analysis (PCA) of hexane (blue dots), acetone (red dots), and ethanol (orange dots) extracts for the Kompolti variety of hemp. The score plot shows the first three PCs (PC1, PC2 and PC3) with their respective variation. R2X(PC1) = 51.8%, R2X(PC2) = 17.5%, R2X(PC3) = 15.2%.
Figure 10. Principal component analysis (PCA) of hexane (blue dots), acetone (red dots), and ethanol (orange dots) extracts for the Kompolti variety of hemp. The score plot shows the first three PCs (PC1, PC2 and PC3) with their respective variation. R2X(PC1) = 51.8%, R2X(PC2) = 17.5%, R2X(PC3) = 15.2%.
Molecules 27 03509 g010
Figure 11. Biplot of PCA carried out on NMR spectra of hexane (brown dots), acetone (purple dots) and ethanol (green dots) extracts of the Tiborszallasi variety of hemp. The score plot shows the first two PCs (PC1 and PC2) with their respective variations. R2X(PC1) = 45.5%, R2X(PC2) = 30.2%.
Figure 11. Biplot of PCA carried out on NMR spectra of hexane (brown dots), acetone (purple dots) and ethanol (green dots) extracts of the Tiborszallasi variety of hemp. The score plot shows the first two PCs (PC1 and PC2) with their respective variations. R2X(PC1) = 45.5%, R2X(PC2) = 30.2%.
Molecules 27 03509 g011
Table 1. 1H NMR chemical shifts and 1H/1H correlations of fatty acid protons in triacylglicerols (TAGs) in CDCl3 for hemp seed extracts.
Table 1. 1H NMR chemical shifts and 1H/1H correlations of fatty acid protons in triacylglicerols (TAGs) in CDCl3 for hemp seed extracts.
PositionδH, Multiplicity a (J in Hz)COSY
A2.27–2.37, mE
B1.55–1.67, mC, A
C1.23–1.39, mGω6 B, D
D1.98–2.11, mGω3, C, F, E
E5.28–5.42, mD, F
F2.72–2.86, mD, E
Gω30.97, tD
Gω60.88, mC
H, L (Gly a)4.14, dd (11.88, 5.93)H’, L’, I
I (Gly a)5.26, mH, H’, L, L’
H’, L’ (Gly a)4.29, dd (11.88, 4.31)H, L, I
a Abbreviations: d—doublet; dd—doublet of doublet; t—triplet; m—multiplet; Gly—Glycerol.
Table 2. 1H and 13C chemical shifts of the main cannabinoids in the flower extracts of Cannabis sativa (Tiborszallasi variety) in CDCl3.
Table 2. 1H and 13C chemical shifts of the main cannabinoids in the flower extracts of Cannabis sativa (Tiborszallasi variety) in CDCl3.
Compoundδ 1H ppm (Multiplicity *, 1H-1H J-Coupling—Hz)δ 13C ppm
CBDH33.86 (ddt; JH3-H4 = 13.00 Hz (d), JH3-H2 = 3.51 Hz (d), JH3-H5 = 2.51 Hz (t))C337.01
H25.57C2124.14
H6a2.05–2.09C630.36
H6b2.22
H51.78–1.84 (ddd; JH5-H4 = 5.30 Hz (d), JH5-H6a = 1.30 Hz (d), JH5-H6b = 0.60 Hz (d))C528.35
H42.40 (dd; JH4-H3 = 13.00 Hz (d), JH4-H5 = 5.00 Hz (d))C446.16
H71.79 (d;3JH7-H2 = 0.50 Hz)C723.69
H9trans4.64 (dq; J9trans-9cis = 2.65 Hz (d), 3J9trans-10 = 1.50 Hz (q))C9110.81
H9cis4.53 (dq; J9cis-9trans = 2.65 Hz (d), 3J9cis-10 = 0.92 Hz (q))
H101.66 (dd; 3J10-9cis = 0.92 Hz (d), 3J10-9trans = 1.50 Hz (d))C1020.30
H2′6.26C2′109.56
H4′6.16C4′107.92
H1″2.43 (t)C1″ 35.46
H2″1.52–1.61C2″30.65
H3″, H4″1.27–1.32C3″31.48
C4″22.54
H5″0.86–0.88C5″14.04
CBDAH34.08C335.38
H25.55C2124.14
H6a2.05–2.09C630.36
H6b2.22
H51.79 (ddd; JH5-H4 = 5.30 Hz (d), JH5-H6a = 1.30 Hz (d), JH5-H6b = 0.60 Hz (d))C528.35
H42.40 (dd; JH4-H3 = 13.00 Hz (d), JH4-H5 = 5.00 Hz (d))C446.45
H71.79 (d; 3JH7-H2 = 0.50 Hz)C723.69
H9trans4.51 (dq; 3J9cis-9trans = 3.00 Hz (d); 3J9trans-10 = 1.76 Hz (q))C9111.21–111.25
H9cis4.39 (dm; 3J9cis-9trans = 3.00 Hz (d))
H101.70C1018.91
H4′6.21C4′111.21–111.25
H1″a2.81C1″36.68
H1″b2.92
H2″1.52–1.61C2″31.24
H3″, H4″1.27–1.32C3″31.94
C4″22.54
H5″0.86–0.88C5″14.04
CBGH26.24C2108.25
H5′, H4′2.04C4′32.28
C5″26.51
H6′5.12C6′125.08
H8′, H10′1.68C8′20.51
C10″23.44
* Abbreviations: d—doublet; t—triplet; q—quadruplet; m—multiplet; dd—doublet of doublet; ddd—doublet of doublet of doublet; ddt—doublet of doublet of triplet; dq—doublet of quadruplet; dm—doublet of multiplet.
Table 3. 1H NMR data of the main cannabinoids in Timborzallasi inflorescences compared with the GC-FID method.
Table 3. 1H NMR data of the main cannabinoids in Timborzallasi inflorescences compared with the GC-FID method.
CompoundqNMR on Flowers UAE ExtractsGC-FID
qNMR ISHexaneAcetoneEthanol
CBDA content *Anthracene6.3 ± 0.80.40 ± 0.10.31 ± 0.04Validated Laboratory Method
6.9 ± 0.2
Referred to CBD
after decarbossilation
Benzoic acid6.5 ± 0.80.40 ± 0.10.39 ± 0.06
TMSP-d46.4 ± 0.60.41 ± 0.10.41 ± 0.06
CBD content *Anthracene0.4 ± 0.14.60 ± 0.80.30 ± 0.1
Benzoic acid0.30 ± 0.064.54 ± 0.62.2 ± 0.1
TMSP-d40.4 ± 0.14.59 ± 0.62.9 ± 0.1
Δ9-THC content *Anthracene0.11 ± 0.44<LOD<LODRegulation (EU) N° 639/2014 [15]
0.09 ± 0.01
Benzoic acid0.07 ± 0.02
TMSP-d40.10 ± 0.02
* % of dry weight.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Colella, M.F.; Salvino, R.A.; Gaglianò, M.; Litrenta, F.; Oliviero Rossi, C.; Le Pera, A.; De Luca, G. NMR Spectroscopy Applied to the Metabolic Analysis of Natural Extracts of Cannabis sativa. Molecules 2022, 27, 3509. https://doi.org/10.3390/molecules27113509

AMA Style

Colella MF, Salvino RA, Gaglianò M, Litrenta F, Oliviero Rossi C, Le Pera A, De Luca G. NMR Spectroscopy Applied to the Metabolic Analysis of Natural Extracts of Cannabis sativa. Molecules. 2022; 27(11):3509. https://doi.org/10.3390/molecules27113509

Chicago/Turabian Style

Colella, Maria Francesca, Rosachiara Antonia Salvino, Martina Gaglianò, Federica Litrenta, Cesare Oliviero Rossi, Adolfo Le Pera, and Giuseppina De Luca. 2022. "NMR Spectroscopy Applied to the Metabolic Analysis of Natural Extracts of Cannabis sativa" Molecules 27, no. 11: 3509. https://doi.org/10.3390/molecules27113509

APA Style

Colella, M. F., Salvino, R. A., Gaglianò, M., Litrenta, F., Oliviero Rossi, C., Le Pera, A., & De Luca, G. (2022). NMR Spectroscopy Applied to the Metabolic Analysis of Natural Extracts of Cannabis sativa. Molecules, 27(11), 3509. https://doi.org/10.3390/molecules27113509

Article Metrics

Back to TopTop