Analytical and in silico study of the inclusion complexes between tropane alkaloids atropine and scopolamine with cyclodextrins

Datura stramonium L. (Solanaceae) possesses a rich tropane alkaloids (TAs) spectrum. The plant contains, in particular, the allelopathic compounds scopolamine and atropine, which are poorly soluble in water, thus limiting their use in agrochemical formulations as biocidal and deterrent agents against herbivore insects. The efficacy of the hydrophobic TAs extracts could be increased with the improvement of their dissolution/leaching properties. This is important for improving screening and test performance and for elucidating the activity of environmentally friendly agricultural approaches, with new perspectives for the production and use of those biodegradable insecticidal products. The present study explores the aspects of atropine and scopolamine complexation with cyclodextrin (CDs) through FT-IR and UV–Vis spectroscopies. In addition, the structures of the inclusion complex of atropine, scopolamine and β-CD have been investigated by molecular modeling techniques. The results obtained indicate that β-CDs are a promising carriers for improving the properties of TAs, therefore increasing their application potential in agrochemical formulations.


Introduction
The allelopathy, a phenomenon of biochemical interactions among plants, is considered one of the possible alternatives for reaching sustainable weed management (Cheng and Cheng 2015;Singh et al. 2001Singh et al. , 2003. Allelopathic compounds (ACs) are bioactive secondary metabolites produced by plants that exert ecological functions to counterbalance abiotic and biotic stressors (Cheng and Cheng 2015;Gross 2009;Mithofer and Boland 2012). They play a key role in the agroecosystems leading to a wide array of interactions between plants such as crop-crop, crop-weed, and tree-crops (Farooq et al. 2011). Allelopathy is consequently considered a valid option for controlling weeds, either through directly utilizing natural allelopathic interactions between plants or by using ACs as natural herbicides (Jabran et al. 2015;Singh et al. 2003). ACs play a major role in plant diversity, Datura stramonium L., known also with the common names of thorn apple, jimsonweed, or devil's snare, is a plant of the Solanaceae family which possesses a rich alkaloid spectrum (Berkov et al. 2006). The main TAs components are atropine and scopolamine (Fig. 1), and their proportions differ considerably between plant parts (leafs, roots, stems, capsules, petioles and seeds) and developmental stages, but also for environment, geographic location, temperature and season (Jakabová et al. 2012). Iranbakhsh et al. identified the highest concentration of atropine and scopolamine in leafs in a vegetative phase (0.127% dried-weight, with atropine:scopolamine ratio of 0.4:1), while their content has been showed to decrease rapidly in the generative phase (0.050% with a ratio of 1.5:1) (Iranbakhsh et al. 2006). In the petioles, the trend from the vegetative to the generative phase showed a decrease in the concentration of the two alkaloids from 0.122 to 0.082% and an increase in their ratio from 2:1 to 3:1 (Iranbakhsh et al. 2006). In addition to scopolamine and atropine, which account for about 80% of total ATs, D. stramonium L. also contains other alkaloids, including hyoscyamine, cuscohygrine, and toluidine (Steenkamp et al. 2004). All these components are poorly soluble in water, thus limiting their use in agrochemical formulations.
The cyclic oligosaccharides cyclodextrins (CDs) can be used to increase solubility, stability, and bioavailability of several classes of molecules, comprising ACs. CDs are cylinder-shaped nanostructure composed of six or more α-1,4linked glucose units. In particular, CDs of six, seven, and eight α-d-glucose residues are known as α-, β-, and γ-CD, respectively ( Fig. s1) (Crini 2014).
CDs possess a hydrophilic exterior, making them soluble in aqueous media, and a hydrophobic cavity that can host guest molecules of different sizes. Non-polar chemicals, such as ACs, can be encapsulated through hydrophobic interactions inside the cavity, forming nanoscale molecular inclusion complexes in solution, and in the solid-state (Bayomi et al. 2002;Kfoury et al. 2016;Nguyen et al. 2013). The main characteristics of CDs are bioadaptability, absorbability, interaction with bio-membranes, soil resistance, and biodegradability (Dodziuk 2006;Sharma and Baldi 2016). CD complexation leads to alterations in the physicochemical properties of guest molecules, such as solubility, chemical stability, dispersibility, flowability, wettability and dissolution rate. All these factors give to these complexes a considerable potential in sustainable agriculture (Mohapatra et al. 2018). For instance, host molecules entrap hydrophobic pesticides inside their cavity and enhance their solubility, protect them from hydrolysis, reduce their volatility, and increase their bioavailability (Villaverde et al. 2004(Villaverde et al. ,2005Yanez et al. 2012). AC/CDs conjugations can be useful for accurate, ground-level administration, time-controlled-release, and time extended-release of ACs in combination with other substances/tools. Moreover, multifunctional CDs have the advantages of protecting sensitive compounds from the effect of heat-or light-induced degradation (Bayomi et al. 2002;Garcia et al. 2014).
Hydroxypropyl-β-cyclodextrin (HPβCD) is widely used as a solubilizing agent and carrier with rapid dissolution for much poorly water-soluble ACs (Perchyonok et al. 2014). Several inclusion complexes with CDs have been reported in the promising field of environmental friendly pesticides (Benfeito et al. 2013;Gao et al. 2019b;Geng et al. 2018;Kim et al. 2019;Villaverde 2007;Villaverde et al. 2004Villaverde et al. , 2005Yanez et al. 2012).
In order to explore rationally new pesticide formulations with enhanced insecticidal or herbicidal activity, it is necessary to assess the influence of CDs on the properties of the guest molecules.
The present study explores the aspects of TAs complexation with CDs obtained by the "wet mixing method" described below. TAs have been obtained from D. stramonium L. extract. The increased TAs hydrosolubility, following the complexation with different types of CDs, have been followed through FT-IR and UV-Vis spectroscopies. Moreover, the structural characteristic of the inclusion complex has been also evaluated with the molecular docking method. The information presented in this research can prove useful in future trends in regard to biocides.

Plant material extraction
Datura Stramonium L. plants were harvested from Banat (Western Romania) in September 2019, at its physiological maturity (the fruits had ripened, and dispersal of the seeds has occurred).

Chemicals
Scopolamine, atropine, cyclodextrins, acetonitrile, methanol HPLC grade, ≥ 99.9% were purchased from Sigma-Aldrich (Germany). The water used in the study was double distilled and deionized with a Milli-Q RG (Millipore, Burlington, MA, USA) water purification system. The standard solution of scopolamine and atropine (10 mg) was dissolved in methanol (100 mL).

Tropane alkaloids extraction procedure
Plant material, in powder (500 mg), obtained from D. stramonium L. species was subjected to extraction three times (30 min each), using chloroform (15 mL), methanol, ammonia (25%) (15: 15: 0.1 (v/v/v) and an ultrasound device (Sonicor Inc., West Babylon, NY, USA). The extract was kept at room temperature for one hour, filtered through filter paper, and washed twice with 1 mL CHCl 3 . The solvent was evaporated to dryness. The dry residue was dissolved in 5 mL of CHCl 3 and 2 mL of H 2 SO 4 (0.5 M). The CHCl 3 organic fraction was removed, whereas the aqueous solution was adjusted to pH 10 by the addition of ammonium hydroxide (25%) in an ice-bath. From the solution, the alkaloids were extracted once with 2 ml and twice with 1 ml of CHCl 3 . Following the addition of anhydrous Na 2 SO 4, the solution was filtered, and the filtrate washed again with CHCl 3 (1-2 mL). The CHCl 3 was further evaporated to dryness under vacuum at 40 °C (Kamada et al. 1986).

Preparation of the inclusion complex
The inclusion complexes of TAs in β-CD were prepared by the wet mixing method according to the following procedure. From the extract of D. stramonium L., one system was made by using the raw alkaloids tropane (mainly composed by atropine and scopolamine), whereas the purchased purified forms have been considered for preparing an additional system.
For each system, five mixtures were prepared, each containing approximately 500 mg TAs (brought to a final volume of 25 mL, far beyond its hydrosolubility) and different β-CD quantities weighed with the accuracy of ± 0.1 mg.
Weight was performed with a semi-micro analytical balance (resolution ± 0.1 mg) (Sartorius). Mixtures were placed in an agate mortar, wetted, and kneaded every 10 min. Subsequently, mixtures were further transferred in a 25 mL volumetric flask and filled to the mark with distilled water and mixed for one hour.
In addition, a saturated solution of TAs in water (revealing some solubility in the absence of CD) was prepared. All solutions were filtered through disposable filters (0.2 µm pore size).

FT-IR and UV-Vis spectroscopy
The infrared absorption spectra were recorded with a Fourier transform infrared spectroscopy (FT-IR), model "460 Plus", Jasco Products Co. Detector: DLATGS with Peltier element (KRS-5), MCT-N (-750 cm −1 ), resolution of 0.9 cm −1 and S/N 15.000:1. Apodization function: Happ-Genzel. Samples were included in a potassium bromide tablet (potassium bromide of spectral purity "Uvasol," Merck, was used). Data processing was done with "Jascow" software (Jasco Company). The FT-IR spectrum of the inclusion complex was compared with those of TAs and hydroxypropyl β-CD (HPβCD). Since water traces from samples interfere with FT-IR spectra, after the mixing operation, the product has been kept for 10 h at 120 °C in an oven. Solubility changes were followed by tracking the absorption spectrum in the ultraviolet domain with a double beam spectrophotometer, model PG Instruments, stray light 0.015% T (220 nm and 340 nm), wavelength accuracy 0.3 nm (automatic wavelength correction), photometric range 0.2-3 Abs, using a UV-Vis, UV WIN 5.05 software. Samples were placed in quartz cuvettes with a 1 cm length path.

Compound design and optimization
The three-dimensional structures of atropine, scopolamine and β-cyclodextrin (β-CD) were designed using USCF Chimera 1.11.2 software (Pettersen et al. 2004). The optimization process for all three compounds was performed with the GAMESS software (Barca 2020;Gordon and Schmidt 2005). For detailed information, see the Supporting Information.

Molecular dynamics (MD) and docking simulations
The molecule of β-CD was submitted to MD simulation with AMBER14 (Case 2014). The optimized conformation of β-CD was employed for the derivation of atomic charges with the RESP methodology (Bayly et al. 1993). Cyclodextrin was solvated in water and periodic boundary conditions were applied with a cutoff of 10.0 Å. The system was heated to 300 K, followed by pressure equilibration at 1 atm. A 150 ns production run was performed using the pmemd implementation in AMBER. All docking calculations were carried out with the AutoDock Vina (Trott and Olson 2010). The structure of β-CD was considered rigid while the docked molecules of atropine and scopolamine were defined as flexible. All molecular graphics were generated using the USCF Chimera 1.11.2 (Pettersen et al. 2004). For further details, see the Supporting Information.

FT-IR characterization of HPβCD-tropane alkaloids inclusion complex
The comparison of FT-IR spectra of the three systems: TAs, extract of D. stramonium L., HPβCD and the inclusion complex of the two components obtained with the wet mixing method is depicted in Fig. 2 and S2. The changes in shift, shape, and intensity of the FT-IR absorption peaks providing clear evidence for the occurrence of the inclusion (Ge et al. 2012).  (Darekar et al. 2016;Gao, 2019a;Su et al. 2012).
Spectral comparison (Fig. 2 and S2) indicated that the two components are not physically mixed but have undergone a significant alteration of their vibration properties (mainly atropine and scopolamine) following their inclusion in the CD cavity with the subsequent restricted dynamic motion. This observation is in accordance with other studies that have been performed in CD inclusion complexes and  . 3 The solubility of pure tropane alkaloids in the presence of β-CD: UV-Vis spectra of tropane alkaloids with different β-CD concentration (a) and the correlation between the absorbance (at 278 nm) and β-CD concentration (b). The solubility of tropane alkaloids raw mixture in the presence of β-CD: UV-Vis spectra of tropane alkaloids mixture with different β -CD concentration (c) and the correlation between the absorbance (at 254 nm) and β -CD concentration (d) by a MD and docking simulations performed in this study (Asztemborska et al. 2019;Mansouri et al. 2013;Wszelaka-Rylik and Gierycz 2015).
In the inclusion complex, the characteristic FT-IR absorption peaks of alkaloid tropane and HPβCD changed. The FT-IR spectrum of the inclusion complex showed no features similar to pure tropane alkaloids. In particular, the absorption peak of alkaloid tropane at 1738 cm −1 (C = O ester) strongly reduce its intensity whereas the peaks at 1650 cm −1 (N-H bending), 1496 cm −1 (C = C aromatic) together with the signals around 1400 cm −1 (C-H bending alkane) almost disappeared, while other bands were obscured by very intense and broad HPβCD bands. Other apparent changes in terms of shift and shapes in several absorption peaks evidence that alkaloids tropane have been included in the HPβCD to form an inclusion complex. Some changes observed in the O-H stretching vibration (absorption peak at 3460 cm −1 ), could allow probing the alteration of the H-bonded environments upon complexation. The O-H bending vibrations of CH 2 OH and C-H stretching vibrations of HPβCD got distorted due to the inclusion of guests to the host cavity, as has been reported in other researches (Darekar et al. 2016).

UV-Vis evaluation of tropane alkaloids water solubility following CDs complexation
To outline the effect of β-CD to enhance atropine and scopolamine water solubility, the absorption spectrum of TAs saturated aqueous solution (sample n° 1, Fig. s3) was compared to some solutions in which, in addition to TAs (atropine and scopolamine from purified extract of D. stramonium L.), increasing quantities of β-CD have been introduced (sample no. 2-6, Fig. s3). Table 1 collects the β-CD solution concentrations (samples n o 1-6) and the related optical absorbance values recorded for purified atropine and scopolamine extract of D. stramonium L and for the raw TAs. Figure 3a shows the relative experimental UV-Vis spectra and the relationship between the concentrations of β-CD (m/v %) and the absorbance of the saturated solution of purified atropine and scopolamine from the extract of D. stramonium L.. The mixture presents low solubility even in the absence of CD (spectrum no. 1) according to absorbance results (Fig. 3a). Increasing concentration of β-CD, in turn, enhances the absorbance of TAs, hence contributing to their water solubility.
The absorption plot (Fig. 3b) shows a second-degree polynomial correlation (R 2 = 0.9996) with increasing β-CD concentration. This correlation points in the direction of saturation phenomena, which is probably because the complex purified TAs-CD has a certain upper limit of solubility.
Experimental UV-Vis spectra and the relationship between the concentrations of β-CD and the absorbance of the saturated solution of raw TAs recorded at 254 nm are shown in Fig. 3c, d (Table 1).
UV-Vis absorption spectra were recorded on filtered solutions in 220-300 nm range. This system reflects better the natural tropane alkaloid components of the extract of D. stramonium L.. As atropine and scopolamine are insoluble in water, their spectra practically coincide with the baseline (spectrum no. 1) (Fig. S3), while in the presence of β-CD, an increase in absorbance was observed (no. 2-6). As can be seen, β-CD increases tropane alkaloid solubility Fig. 4 The effect of various CDs on tropane alkaloids solubility. UV-Vis spectra of tropane alkaloids with different types of CDs: α-cyclodextrin (α-CD), β-cyclodextrin (β-CD), γ-cyclodextrin (γ-CD), hydroxypropyl β-cyclodextrin (HPβCD) and randomly methylated βcyclodextrin (RAMEB) mixture with β-CD (a) and the enlarged scale (b) 1 3 significantly. Unlike pure compounds, an excellent linear correlation (R 2 = 0.9995) was observed between β-CD concentration and TAs mixture solubility in water (in the range tested). This difference can be attributed to the higher solubility of TAs-β-CD complexes compared to allelochemical (pure mixed atropine and scopolamine)-β-CD complexes (at the reported concentrations). This interpretation is fairly in agreement with the β-CD solubilization effect, being higher with raw natural TAs mixture than with pure atropine and scopolamine mixture (Castagne et al. 2010). Figure 4 represents a comparison of the solubilization effects for TAs in complex with different types of CDs: α-cyclodextrin (α-CD), β-cyclodextrin (β-CD), γ-cyclodextrin (γ-CD), hydroxypropyl β-cyclodextrin (HPβCD) and randomly methylated β-cyclodextrin (RAMEB).
Because of the difference in the number of glucopyranose units, each CD has a distinct molecular weight, cavity size, mobility and water solubility. Since α-CD has the fewest glucopyranose units (six), they have the smallest cavity size of 0.57nnm, followed by β-CD which has seven glucopyranose units and a cavity size of 0.78 nm, while, γ-CD, which has eight glucopyranose units, has a cavity size of 0.95 nm. Water solubility depends mainly from the hydroxyl groups interacting much less with each other and much more with water. Since γ-CD ring is more strained compared to α-CD and β-CD, its hydroxyl groups are more distant and more able to interact with water, providing greater solubility. In summary due to the larger size and hence more hydroxyl groups available, the aqueous solubility of γ-CD is greater than α-CD, while β-CD possess strong intramolecular bonding between its hydroxyl groups and thus possess the lowest water solubility. However, the modification of β-CD in HPβCD leads to a large increase in its water solubility (> 60% at 25 °C).
The comparison between the different CDs (Fig. 4) shows a clear superiority effect of RAMEB, respect to all other CDs, in the enhancement of the aqueous solubility of TAs. However, RAMEB possesses poor permeability and high cytotoxicity (Kiss 2010) and its use should be avoided. Therefore, compared to the remaining CDs, the dimensions of the inner cavity of the safer β-CD and HPβCD appear the best to include ATs molecules which lead to their improved hydrosolubility.
To exclude the hypothesis that the solubilization effect is due to a solvotropic effect generated by the CDs and not to the inclusion of TAs in the hydrophobic cavity, a solubilization test was carried out with an equivalent quantity of glucose. The glucose employed contained the same amount of glucopyranose units as the CDs without having a cyclic structure, therefore being deprived of the signature characteristics of CDs. The comparison indicated (Fig. 4b) that the presence of glucose virtually does not have any effect in the solubility of TAs. The outcome is that the solubilization effect can be attributed on the inclusion properties of the cyclic structure of CDs.

Molecular modeling study
The initial conformations of atropine, scopolamine and β-CD were optimized in order to obtain the conformation with the lowest potential energy. The results are presented in Fig. S4 and S5 in Supporting information. The optimized structures do not present any substantial difference compared to the initial conformation (Fig. S5). For β-CD, this observation is more pronounced due to the circular nature of the molecule (Fig. S5c) that does not allow great flexibility.
The analysis of the MD simulation showed that the circular conformation of the molecule impedes the movement of the atoms. The RMS changes observed during the production run are more pronounced in the first 60 ns of the simulation whereas between 70 and 150 ns the conformation of cyclodextrin does not present high RMS fluctuations (Fig.  S6a). The mean RMS over the whole simulation process is 2.25 ± 0.424 Å, a measure that suggests small conformational changes from its starting conformation. The small changes of β-CD show that the solvent molecules may interact favorably with the solute due to the presence of ester and -OH groups in β-CD. The number of -OH groups in cyclodextrin may also lead to increased interactions with other organic molecules via formation of hydrogen bonds. Thus, β-CD proves to be an excellent drug delivery system. The clustering analysis of the MD simulation revealed the presence of three representative conformations for β-CD (Fig. S6b). The superimposition of the structures revealed only small conformational changes ( Fig. S6b and Table S1).
The dominant conformation of β-CD (Fig. S6b, black) was employed in the molecular docking experiments. The results of the docking simulations showed that the atropine molecule presents higher binding affinity to β-CD in contrast to scopolamine. The calculated docking score for atropine and scopolamine is − 4.5 and − 4.0 kcal mol −1 , respectively. The docking scores show that atropine presents a better binding affinity than scopolamine. The analysis of the docking conformations presents the hydrogen bonds (HBs) forming between atropine/ scopolamine and β-CD (Fig. 5a-magenta and 5b-cyan). As expected, the -OH groups of β-CD are involved directly with the formation of HBs. Additionally, atropine and scopolamine are positioned in the cavity of β-CD ( Fig. 5c and d).
The two molecules are anchored in the cavity of cyclodextrin via hydrogen bonds forming with the respective sides of β-CD. The positioning of the molecules inside the cavity of cyclodextrin may improve their water solubility as described in previous sections and thus make β-CD an ideal carrier molecule. Moreover, the HBs formed between the alkaloids and β-CD may improve the bioavailability of the alkaloids and further explain the stability of both the carrier molecule and the substances attached, as reported in the experimental sections above.

Conclusions
A number of complexes of CDs, CD polymers, and conjugations of CDs have been designed and evaluated for practical use in the field of sustainable agriculture. Because of their multifunctional characteristics and bioadaptability, CDs can mitigate the undesirable properties of AC molecules through the formation of inclusion complexes or as AC/CDs conjugates. One of the characteristics of CDs is their tendency to include in their hydrophobic cavity active molecules with hydrophobic properties (partial to complete hydrophobicity). The experimental intermolecular complex presented in this study shows an increased water solubility of tropane alkaloids (TAs) obtained from D. stramonium L. extract due to the hydrophilicity of the external surface of CDs. The increased TAs hydrosolubility has been evaluated, through UV-Vis spectroscopy, following the complexation with different types of CDs. The degree of solubilization effect is dependent on the nature of TAs mixture, the size and characteristics of the cyclic molecule (α-, ß-, HPβ-or -CD) and the CDs concentration. The dimensions of the inner cavity of the safer β-CD and its hydroxyled derivative HPβCD, appear the best to include ATs molecules which lead to their improved hydrosolubility respect to other CDs. FT-IR analysis demonstrated that TAs can be effectively Docking conformation highlighting HB interactions, as yellow dotted lines, for: a atropine (magenta) and β-cyclodextrin and b scopolamine (cyan) and β-cyclodextrin. Positioning of atropine (magenta) and scopolamine (cyan) in the cavity of β-CD, shown as surface: c side view and d top view encapsulated in HPβCD forming inclusion complexes with enhanced hydrosolubility. In addition, molecular docking gives evidences of the better binding affinity of atropine than scopolamine with β-CD.
The results obtained indicate that β-CD and HPβCD are the most efficient carriers for improving the properties of tropane alkaloids. The information gained can prove useful in future trends concerning biocides for the application of tropane alkaloids in agrochemical formulations as natural and biodegradable pesticide products.
Funding: Open access funding provided by Università degli Studi di Sassari within the CRUI-CARE Agreement. This research received no external funding.

Conflicts of interest There is no conflict of interest.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.