Comparison of Five Oral Cannabidiol Preparations in Adult Humans: Pharmacokinetics, Body Composition, and Heart Rate Variability

Data supporting the physiological effects of cannabidiol (CBD) ingestion in humans are conflicting. Differences between CBD preparations and bioavailability may contribute to these discrepancies. Further, an influence of body composition on CBD bioavailability is feasible, but currently undocumented. The aims of this study were to: (1) compare the pharmacokinetics of five oral CBD preparations over 4 h; (2) examine the relationship between body composition and CBD pharmacokinetics; and, (3) explore the influence of CBD on heart rate variability. In total, five preparations of CBD, standardized to 30 mg, were orally administered to 15 healthy men and women (21–62 years) in a randomized, crossover design. Prior to and 60 min following CBD ingestion, heart rate variability was determined. Body composition was assessed using dual energy X-ray absorptiometry. Peak circulating CBD concentration, time to peak concentration, and area under the curve was superior in a preparation comprising 5% CBD concentration liquid. Fat free mass was a significant predictor (R2 = 0.365, p = 0.017) of time to peak concentration for this preparation. Several heart rate variability parameters, including peak frequency of the high frequency band, were favorably, but modestly modified following CBD ingestion. These data confirm an influence of CBD preparation and body composition on CBD bioavailability, and suggest that acute CBD ingestion may have a modest influence on autonomic regulation of heart rate.


Introduction
Cannabidiol (CBD), a non-psychoactive component of Cannabis sativa L., has been purported to have a variety of beneficial physiological effects including but not limited to, pain relief, anti-anxiety, anti-seizure, anti-depressant, anti-oxidant, and anti-inflammatory [1]. Despite the enthusiasm of the CBD industry for these claims, the empirical evidence supporting favorable physiological responses is inconsistent [2,3]. These discrepancies may in part be explained by differences in CBD bioavailability when administered orally [4], that in turn may be influenced by non-standardized CBD formulation, and/or the body size and composition of the recipient. With regards to the former, differences in formulation might include water vs. lipid solubility, administration as a food vs. a beverage, co-administration with additional ingredients, and/or preparation as a ready-made product vs. a powder to be mixed with liquid by the consumer before ingestion. Once ingested, the resultant circulating CBD concentration will be influenced by rates of absorption from the gut, breakdown during first pass metabolism, and potentially by the body size and composition of the consumer. For example, lean mass is positively associated with total blood volume [5,6]. Thus, it is feasible that adults with a higher lean mass may demonstrate lower circulating CBD concentrations on account of a larger blood volume in which to dilute the CBD. Alternatively, fat mass may influence CBD absorption as CBD is lipid soluble and can therefore potentially accumulate in adipose tissue in a manner similar to that previously reported for ∆ 9 -tetrahydrocannabinol, a psychoactive component of Cannabis sativa L. [7,8]. Although intuitive, the influence of body size and composition on CBD pharmacokinetics is currently undocumented. The purpose of the current study was to compare the pharmacokinetics of five oral CBD preparations over 4 h, and to examine the relationship between body composition and CBD pharmacokinetics.
One of the inconsistent physiological responses to CBD administration is the cardiovascular response [9]. Some studies have reported changes in heart rate and/or blood pressure following acute CBD consumption [10,11] while others have reported no measurable effect [12,13]. Changes in heart rate (in any direction) are almost always mediated by changes in the sympathovagal balance; this balance is typically described via heart rate variability. Heart rate variability refers to the inconsistency within the periods of time separating consecutive cardiac cycles. It is considered to have clinical relevance as it is able to predict future cardiac events and mortality [14][15][16]. To our knowledge, only one other study has examined the influence of CBD (specifically, a hemp-oil extract containing CBD) on heart rate variability in adult humans [17], and this study reported on short-term CBD use (3 and 6 weeks) and not an acute response. Accordingly, an additional purpose of the current study was to explore the acute influence of CBD ingestion on heart rate variability.

Participants
The progress of all participants throughout the trial (from screening and enrollment through to completion) is presented in Figure 1. A total of 16 participants were enrolled in the study. Due to repeated scheduling difficulties, one participant was removed from the study and replaced with a new recruit; a total of 15 participants completed all procedures in a crossover design. Selected physiological characteristics are presented in Table 1. Consistent with inclusion and exclusion criteria, the physiological characteristics were unremarkable, and typical of a broad sample of adults, free from overt cardio-metabolic disease.

Circulating CBD and Pharmacokinetics
Each of the CBD preparations contained a standardized CBD dose of 30 mg. Th characteristics of each of the preparations are described in Table 2. Circulating concentra tions of CBD and the calculated pharmacokinetics are presented in Figure 2 and Table  respectively. Baseline (Time 0) concentrations for all CBD preparations, for all partic pants, were below the limit of quantitation. By the final blood collection (4 h), averag concentration was above the limit of quantitation for each preparation, but not for eac participant (1 participant below for 178, 2 participants below for 203, 0 participants below for 340, 1 participant below for 472, and 2 participants below for 707). Informal visu inspection of Figure 2 suggested preparation 707 evoked the highest circulating CBD con centration, the shortest time to maximal concentration, and the greatest area under th curve during the blood collection period. Formal statistical analysis, detailed in Table  supported several of these interpretations, although not all comparisons between prepa ration 707 and the other preparations attained statistical significance (p < 0.05). It was no possible to calculate some pharmacokinetic parameters (i.e., AUCinf, t½, Ke and Vd) for a of the CBD preparations on account of insufficient values above the limit of quantitatio during the initial 1-2 h of blood collection.

Circulating CBD and Pharmacokinetics
Each of the CBD preparations contained a standardized CBD dose of 30 mg. The characteristics of each of the preparations are described in Table 2. Circulating concentrations of CBD and the calculated pharmacokinetics are presented in Figure 2 and Table 3, respectively. Baseline (Time 0) concentrations for all CBD preparations, for all participants, were below the limit of quantitation. By the final blood collection (4 h), average concentration was above the limit of quantitation for each preparation, but not for each participant (1 participant below for 178, 2 participants below for 203, 0 participants below for 340, 1 participant below for 472, and 2 participants below for 707). Informal visual inspection of Figure 2 suggested preparation 707 evoked the highest circulating CBD concentration, the shortest time to maximal concentration, and the greatest area under the curve during the blood collection period. Formal statistical analysis, detailed in Table 3, supported several of these interpretations, although not all comparisons between preparation 707 and the other preparations attained statistical significance (p < 0.05). It was not possible to calculate some pharmacokinetic parameters (i.e., AUC inf , t 1 2 , K e and V d ) for all of the CBD preparations on account of insufficient values above the limit of quantitation during the initial 1-2 h of blood collection.   20,178 ± 9989 a n = 8 8024 ± 6630 a n = 15 Data are mean and SD. Limit of quantitation: 0.25 ng/mL. Values below limit of quantitation were classed as "missing". n: number of observations used to calculate parameter. T max : the time to maximum concentration. C max : the maximum concentration. AUC 0-4 : the area under the curve representing total cannabidiol exposure between 0 and 4 h. AUC 0-inf : an estimate of the total exposure to cannabidiol over time. t 1 2 : the amount of time it takes to decrease the circulating concentration to half of its initial value. K a : the rate at which the cannabidiol is absorbed into the body. K e : the rate at which the cannabidiol is removed from the body. V d : the volume of distribution, an estimate of the degree to which cannabidiol is distributed in the body tissue vs. the plasma. Values sharing the same superscript letter are different (p < 0.05).

Body Composition
Within each CBD preparation, there was considerable variability in the calculated pharmacokinetic parameters. Correlation analysis revealed several associations between some of these parameters and selected physiological characteristics.

Heart Rate Variability
Heart rate variability was assessed immediately prior to and 60 min following CBD ingestion. Circulating CBD concentrations for preparations 178 and 203 were below the limit of quantitation at 60 min, thus heart rate variability data associated with these preparations were excluded from statistical analysis. The heart rate variability data are presented in Table 4. There were main effects of time (all p < 0.05) for heart rate (decreased), R-to-R interval (increased), peak frequency of the high frequency band (decreased), Poincaré plot standard deviation perpendicular the line of identity (increased), and Poincaré plot standard deviation along the line of identity (decreased). In addition, there were several parameters with main effects of time that did not attain statistical significance (0.05 < p < 0.08) including, baseline width of the R-to-R interval histogram (increased), and the ratio of Poincaré plots standard deviation perpendicular to along the line of identity (decreased). There were no time x CBD preparation interactions (all p > 0.10). There were no appreciable changes in heart rate and blood pressure across the 4 h of data collection (Table 5).  Data are mean ± SD. Data were collected over 10 min during paced breathing. R-to-R int : R-to-R interval. SDNN: Standard deviation of N-to-N intervals. RMSSD: Root mean square of successive R-to-R interval differences. R-to-R Triangular Index: Integral of the density of the R-to-R interval histogram divided by its height. TINN: Baseline width of the R-to-R interval histogram. b Denotes non-significant effect of time (0.05 < p < 0.08). There were no time x CBD preparation interactions (all p > 0.10). Circulating CBD concentrations for preparations 178 and 203 were below the limit of quantitation at 60 min, thus heart rate variability data associated with these preparations were excluded from statistical analysis.

Discussion
The goals of the current study were to compare the pharmacokinetics of five oral CBD preparations over 4 h, to examine the relationship between body composition and CBD pharmacokinetics, and to explore the influence of acute CBD ingestion on heart rate variability. Our primary findings were: (1) compared with most of the other preparations, the preparation comprising 5% CBD concentration liquid (preparation 707; Table 2) evoked the fastest T max , and the greatest C max ; (2) within each CBD preparation, there was considerable variability in the calculated pharmacokinetic parameters; some of the variability for some of the preparations could be explained by body size and composition; and, (3) CBD had only a modest effect on some of the parameters used to describe heart rate variability.
The high degree of variability in CBD absorption and the resulting pharmacokinetic parameters following CBD ingestion have been previously described [4,12,[18][19][20]. Potential explanations for variable, and often poor, CBD absorption have pertained to properties of the CBD itself (e.g., poor water solubility) and also variability brought about by the significant first pass metabolism it must endure following oral ingestion. To counter these potential absorption hurdles, CBD preparations have been modified, with mixed success, in a variety of ways, including generation of synthetic CBD [20], creation of water-soluble CBD powders [12], development of self-emulsifying delivery systems [18,20], and encapsulation of CBD within gelatin matrix pellets [21] and liposomes [22]. One important consideration when evaluating CBD formulations is the pharmacokinetic goal and intended use of the CBD. For example, if the indication for the CBD is to treat acute pain, then a faster T max and higher C max may be desirable, and also may help to decrease the risk of overdose due to premature repeat self-administration. Alternatively, as a chronic treatment for anxiety, a larger AUC may be preferrable if a user follows a regular dosing schedule. The CBD preparations under investigation in the current study generated a wide range of values for pharmacokinetic parameters (Table 3). Preparation 707 appeared to be the superior preparation as, compared with most of the others, it generally evoked the fastest T max , greatest C max , and largest AUC. Perhaps more important than comparisons within the current study are comparisons of the preparations to others described within the literature. In this regard, a very thorough review of CBD pharmacokinetics was published in 2018 [19]. We have supplemented this published information to include data from the current study, and other data published since the release of the original review (Table 6). An important caveat to these comparisons is that the dose of CBD administered across the reviewed studies ranged between 5 and 6000 mg; the influence of CBD dose on pharmacokinetics, especially AUC, has been previously described [19,23], thus direct comparisons to appreciably smaller/larger doses are problematic. Further, these are informal and not definitive statistical/quantifiable comparisons. With these caveats in mind, preparation 707 from the current study appears to have the fastest T max , but a smaller C max when compared with the Self-Emulsifying Drug Delivery System, despite the latter comprising a slightly smaller dose (25 mg) [18]. The favorable CBD bioavailability of the Self-Emulsifying Drug Delivery System was attributed to augmented solubility within the gastrointestinal tract [18]. A similar rationale was employed during the development of preparation 707; some of the components of the preparation, such as gum arabic and medium-chain triglycerides (MCTs) have been shown to improve bioavailability of dietary supplements [24,25] and may have contributed to the superior performance in the current study, and its favorable comparison to the other published formulations. : the amount of time it takes to decrease the circulating concentration to half of its initial value. K a : the rate at which the cannabidiol is absorbed into the body. K e : the rate at which the cannabidiol is removed from the body. V d : the volume of distribution, an estimate of the degree to which cannabidiol is distributed in the body tissue vs. the plasma. The information presented in this table represents an update to a previous review [19]. * PTL-gelatin matrix pellets technology-based formulation; ** MCT-CBD-medium-chain triglycerides; *** SEEDS-CBD-self-emulsifying drug delivery system. **** Caliper CBD Water Soluble is identical to Preparation 472. The background color highlights the preparations investigated in the current study.
Another important consideration for CBD use is terminal elimination half-life. In the context of pharmacokinetic research, undetectable baseline circulating CBD concentrations can be critical for valid derivation of pharmacokinetic parameters. The terminal elimination half-life of CBD administered to fasted humans varies greatly (between 6 and 32 h) [37,38] depending on mode of administration and CBD preparation. This variability is further illustrated by the K e values (i.e., the rate of CBD removal from the body) presented in Tables 3 and 6. In the current study, baseline (Time 0) concentrations for all CBD preparations, for all participants, were below the limit of quantitation. This provides support for our rationale for a minimum 72 h "washout" period between consecutive laboratory visits, although it should be noted that other studies have employed considerably longer duration (14 days) CBD abstention periods between repeated study visits [18].
An influence of body size and composition on CBD pharmacokinetics has been previously speculated [18] but to our knowledge is currently undocumented. In the current study we used DEXA to quantify body composition. Some of the inter-personal variability, for some of the preparations, could be explained by body size and composition as reflected by correlations between several of the DEXA-derived variables and CBD pharmacokinetic parameters. It is not uncommon for DEXA-derived variables to be interrelated (for example, high fat free mass is often associated with high body mass), thus we employed stepwise regression to further explore potential relations between body composition and CBD pharmacokinetics. With respect to our superior CBD preparation, preparation 707, fat free mass was a significant predictor of T max (r = −0.60, p = 0.017).
Somewhat counter-intuitively, the association was negative, implying that a larger fat free mass was associated with a faster T max . In light of the positive association between fat free mass and total blood volume [5,6], it may seem unlikely that peak concentration is reached earlier in adults with greater blood volumes. Circulating concentrations of all substances, including endogenous hormones and active ingredients in dietary supplements, are determined by a combination of rate of entry and rate of clearance from the blood. The shorter time to peak concentration in adults with greater fat free mass may reflect a greater rate of clearance by metabolically active and relatively well-perfused tissues (such as skeletal muscle). Noteworthy, fat free mass did not predict C max or AUC for preparation 707. Alternatively, body mass index was positively associated with T max for preparation 178 (r = 0.63, p = 0.016). The preparation 178 T max was among the slowest of the compared preparations (approximately 3.5 hours); it is feasible that for preparations absorbed more slowly (i.e., slower rate of entry into the blood), body size (and presumably blood volume) may contribute to time to peak circulating concentration. From a broader perspective, CBD bioavailability may be determined not only by preparation/formulation, but perhaps also by body size and composition. The implication is that future dosing guidelines may better serve users of CBD when additional consideration is given to expected rates of CBD absorption based on preparation/formulation characteristics (e.g., presence/absence of ingredients used to promote solubility), requirements of intended use (e.g., acute needs from a single dose vs. chronic repeated dosing), and the body size and composition of the user. Dosing guidelines should be a focus of future research as we are unable to make definitive recommendations based on our current data.
An additional goal of the current study was to explore the influence of acute CBD ingestion on heart rate variability, the inconsistencies in the periods of time separating consecutive cardiac cycles. It is an important predictor of future cardiac events [14,15] and, in some instances, longevity [16]. While there are many methods and different parameters used to quantify heart rate variability, generally speaking, greater variability is considered desirable. To our knowledge, this is the first exploration of a potential influence of acute CBD on heart rate variability in adult humans. Only one previous study has examined the short-term influence of CBD (administered as a hemp oil extract) on heart rate variability [17] but the experimental design did not include an acute measurement. Other studies have compared heart rate variability in users and non-users of Cannabis sativa L. (i.e., CBD plus ∆ 9 -Tetrahydrocannabinol and other potentially active ingredients), thus making assertions as to the independent influence of CBD difficult [39]. In the current study, acute CBD ingestion decreased heart rate, peak frequency of the high frequency band, and Poincaré plot standard deviation along the line of identity, and increased R-to-R interval and Poincaré plot standard deviation perpendicular to the line of identity (Table 4). Poincaré plots are a useful tool in the heart rate variability analysis arsenal; they are often able to recognize patterns and rhythms in the R-to-R interval data that are sometimes overlooked using spectral analysis [40,41]. From an interpretation perspective, the standard deviation perpendicular to the line of identity reflects short-term heart rate variability and is thought to be regulated by parasympathetic (vagal) input, while the standard deviation along the line of identity reflects total heart rate variability and is determined by a combination of sympathetic and parasympathetic input [40]. When considered together, the ratio of these two variables (i.e., the SD2/SD1 ratio; Table 4) represents sympathovagal balance. In light of the opposing direction of change in these Poincaré plot data, in addition to the small magnitude of change in the other heart rate variability parameters, it appears likely that acute CBD evokes only a modest, and physiologically irrelevant effect on heart rate variability, at least when studied at rest. In the short-term study [17], neither 3 nor 6 weeks of daily CBD ingestion influenced heart rate variability measured at rest or during mental stress activities. One important point worthy of emphasis in the current study-our primary focus was the pharmacokinetics of the different CBD preparations and the potential interaction with body size and composition. Accordingly, we did not include a placebo-control within our experimental design. Thus, any conclusions with respect to our exploration of the acute influence of CBD on heart rate variability must be considered with caution. For example, it is plausible that CBD does not exert any acute influence on heart rate variability and we are simply reporting on the influence of 60 min of semi-recumbent rest.
CBD may or may not have appreciable, and potentially favorable effects on physiological function. Conflicting empirical evidence describing the response to CBD can be partially explained by considerable inconsistencies in its bioavailability. In the current study we identified a CBD preparation that compared well with other preparations described in the literature [19], and with the other preparations incorporated within our design. Further, we demonstrate that some of the pharmacokinetic parameters of this superior preparation were influenced by body size and composition. These findings highlight the need to optimize CBD preparations and personalized dosing strategies in order to confirm/refute the physiological relevance of CBD as a potential therapeutic. Finally, we document, for the first time, a modest effect of acute CBD on selected parameters of resting heart rate variability in healthy adults, alleviating any potential concerns that CBD may unfavorably impact autonomic regulation of heart rate.

Materials and Methods
A randomized, double-blind, repeated measures cross-over study design was employed. The institutional review board at Colorado State University reviewed and approved all procedures in accordance with the principles established in the Declaration of Helsinki. Written informed consent was provided by all participants prior to commencement of any study activity.

Participants
Healthy adult men and women were invited to participate. Inclusion criteria included aged 18 years and over, body mass greater than 50 kg, absence of any known gastrointestinal or metabolic disease, and willing to refrain from all products containing CBD for 3 days prior to each study visit. Exclusion criteria included pregnancy, breastfeeding, known food allergies, autoimmune disorders, celiac disease, inflammatory bowel disease, gastrointestinal cancers, and history of diabetes. In addition, anyone who reported experiencing a previous adverse reaction to ingesting products containing Cannabis sativa L. were excluded from participation.

Protocol Overview
Following a screening visit consisting of a review of medical history and assessment of body composition, participants reported to the laboratory on five separate mornings in a crossover design. In light of the considerable variability in CBD absorption and the resulting pharmacokinetic parameters following CBD ingestion [4,12,[18][19][20], a crossover design was employed in an attempt to minimize the influence of inter-individual variability. Each morning began with a baseline blood sample and determination of blood pressure, and heart rate variability. Participants then ingested one of five CBD preparations (described below). Venous blood was sampled over the next 4 hours. Heart rate variability was reassessed after one hour. Heart rate and blood pressure were also recorded over the next four hours.

Procedures
Participant screening comprised completion of a detailed medical history questionnaire and assessment of body size and composition using a physician's digital scale and dual energy X-ray absorptiometry (Hologic, Discovery W, QDR Series, Bedford, MA, USA) as previously described [42].
The remaining five laboratory visits were identical in all aspects except for the CBD preparation that was administered. Participants reported to the laboratory on five separate mornings, each preceded by a 12 h fast. The terminal elimination half-life of CBD administered to fasted humans varies between 6 and 32 h [37,38]; to facilitate negligible baseline circulating CBD concentrations in each participant at the start of each study visit, every visit was preceded by a minimum 72 h abstention from all products containing CBD. The time of arrival was kept constant for each participant. On arrival, participants were instrumented for measurement of heart rate (3-lead electrocardiogram (ECG)) and blood pressure (auscultation) using a physiological monitor (IntelliVue MP5 Patient Monitor, Philips Healthcare, Andover, MA, USA); a venous catheter was introduced to an antecubital vein.
Heart rate variability was determined immediately prior to, and 60 min following CBD ingestion. During 11 min of paced-breathing (metronome: 6 breaths per minute) raw ECG signals were recorded using a personal computer and analogue-to-digital convertor (WinDaq, Dataq Instruments Inc., Akron, OH, USA) [43] and analyzed offline using commercially available software (Kubios HRV), as previously described [44,45]. The first minute was considered habituation and was discarded from the analysis. Additional measurements of heart rate and blood pressure were made prior to, and 30, 60, 120, 180 and 240 min after CBD ingestion. All measurements and recordings were completed in a temperature-controlled (20-22 • C), dimly lit room. During the measurements and recordings, participants were situated on a medical bed, in a comfortable, semi-recumbent posture.
Then, 90 min following CBD ingestion, participants were provided with a standardized breakfast consisting of bagel sandwich and choice of a non-alcoholic beverage. The breakfast was different between participants but was constant within laboratory visits (i.e., whatever a participant ate during the first visit, they ate for all of the subsequent visits).

Cannabidiol Preparations
In total, five different CBD preparations were provided by Caliper Foods (Commerce City, CO, USA). Each of the preparations contained a standardized CBD dose of 30 mg. The characteristics of each of the preparations are described in Table 2. Broadly speaking, the CBD preparations differed in their solubility (i.e., water vs. lipid), the concentration of the CBD powder, and additional ingredients used in each specific formulation. All preparations were administered in 227 mL (8 oz) of water and were consumed within~30 s of administration.
Plasma samples were prepared for LC-MS/MS analysis by protein precipitation and dispersive solid phase extraction. In total, 200 µL of sample was mixed with 400 µL of ice-cold acetonitrile containing 20 ng/mL CBD-D3 in a microcentrifuge tube and vortexed for 30 s to precipitate proteins. Next, 200 mg of dSPE sorbent was added, vortexed for 30 s, and centrifuged for at 14,000 rpm for 5 min. Sample supernatants were then transferred to an autosampler vial for LC-MS/MS analysis.
Samples were analyzed with an Agilent 1290 UHPLC coupled to an Agilent 6460 triple quadruple mass spectrometer equipped with an Agilent Jet Stream electrospray ionization source (Agilent, Santa Clara, CA, USA). Cannabinoids were first chromagraphically separated on a Phenomenex Phenyl Hexyl column (3.0 mm × 50 mm, 2.6 µm) held at 40 • C. A sample volume of 10 µL was injected and a mixture of water with 0.1% formic acid (A) and acetonitrile with 0.1% formic acid (B) at a flow rate of 0.4 mL/min. The gradient elution used was 40% B for 0.5 min, increasing to 100% B at 2 min, and held at 100% B for 1.5 min. The ionization source conditions used were as follows: positive polarity, nebulizer 45 psi; gas flow of 10 L/min at 300 • C; sheath gas flow of 12 L/min at 375 • C. The ion transitions monitored for CBD were 315.2 → 193.1/123 m/z and 318.2 → 196.1/123.1 m/z for CBD-D3. CBD was confirmed by retention time and the product ion ratio correlation between the sample peaks and corresponding standards (±20%). The limit of detection and limit of quantitation for CBD in this analysis was 0.1 ng/mL and 0.25 ng/mL, respectively. The data collection and processing were performed by using Agilent MassHunter Quantitative software (v.B.08.01). Quantitative analysis was performed with linear regression using a 6-point calibration curves from 0.25 ng/mL to 50 ng/mL.

Pharmacokinetic Analysis
Pharmacokinetic analysis of the circulating concentrations of CBD for each of the preparations was completed using dedicated software (Phoenix WinNonlin v8.2, Certara, NJ, USA). Values below the limit of quantitation were classified as "missing". Areas under the CBD concentration curves were calculated using the trapezoidal method.

Statistical Analysis
All data, unless otherwise stated, are expressed as mean and standard deviation. Statistical calculations were performed using dedicated software (Prism v8.4.3, GraphPad Software, San Diego, CA, USA). Differences in the pharmacokinetic properties between the CBD preparations were examined using 1-way analysis of variance mixed-effect models with Tukey tests employed to further explore identified main effects. Similarly, differences in the characteristics describing heart rate variability before/after CBD ingestion were also examined using 1-way analysis of variance mixed-effect models, and Tukey test when appropriate. Relations between CBD pharmacokinetic parameters and body size and composition values were explored using Pearson correlations, and further examined with forward stepwise regression when pharmacokinetics parameters were correlated with multiple body composition variables (SigmaStat 3.0, Systat Software Inc., San Jose, CA, USA). The level of statistical significance was set at p < 0.05.  Informed Consent Statement: Informed consent was obtained from all subjects involved in the study.

Data Availability Statement:
The data presented in this study are available as supplementary material.