Study to Probe Subsistence of Host-Guest Inclusion Complexes of α and β-Cyclodextrins with Biologically Potent Drugs for Safety Regulatory Dischargement

Host-guest interaction of two significant drugs, phenylephrine hydrochloride and synephrine with α and β-cyclodextrins were studied systematically. Initially two simple but reliable physicochemical techniques namely conductance and surface tension were employed to find out saturation concentration for the inclusion and its stoichiometry. The obtained 1:1 stoichiometry was further confirmed by two spectrometric methods, UV-Vis study and spectrofluorimetry. Significant shifts in IR stretching frequency also support the inclusion process. Relative stabilities of the inclusion complexes were established by the association constants obtained from UV-Vis spectroscopic measurements, program based mathematical calculation of conductivity data. Calculations of the thermodynamic parameters dictates thermodynamic feasibility of the inclusion process. Spectrofluorometric measurement scaffolds the UV-Vis spectroscopic measurement validating stability of the ICs once again. Mass spectroscopic measurement gives the molecular ion peaks corresponding to the inclusion complex of 1:1 molar ratio of host and guest molecules. The mechanism of inclusion was drawn by 1H-NMR and 2D ROESY spectroscopic analysis. Surface texture of the inclusion complexes was studied by SEM. Finally, the cytotoxic activities of the inclusion complexes were analyzed and found, Cell viability also balances for non-toxic behavior of the ICs. Moreover, all the studies reveal the formation of inclusion complexes of two ephedra free, alternatively emerging drugs (after their banned product having ephedra) SNP, PEH with α and β-CD which enriches the drug delivery system with their regulatory release without any chemical modification.

In supramolecular chemistry of cyclodextrins various guest molecules having hydrophobic part, influenced by non-covalent interaction, can be inserted into the hydrophobic cavity of cyclodextrin molecules. Cyclodextrins (α-CD, β-CD, γ-CD) having six, seven and eight numbers of glucopyranose units respectively, (Fig. 1) produced from starch by the enzymatic conversion, have different cavity sizes. Inclusion complexes (ICs) with structures of higher complexity in the solid state and solution phase can increase the aqueous solubility of various drugs cum bio-active molecules of merely water solubility which leads to the development of drug delivery systems 1,2 . Chiral Separation of molecules using cyclodextrins as chiral additives are also possible by applying capillary electrophoresis (CE) and electrochemical detection (ED) method [3][4][5][6] . Structural characterization of Host-Guest inclusion complexes of α-CD and β-CD with two bio-active molecules, PEH and SNP were done over here in terms of geometry and structural preferences by means of a variety of physical and spectroscopic methods in solid state and solution phase.
Phenylephrine hydrochloride (PEH) (Fig. 1) is a selective α 1 -adrenergic receptor agonist of the phenethylamine class used primarily in cold and flu conditions as an antipyretic, analgesic drug to relief pain 7 . In the  Cyclodextrins, mostly α-cyclodextrin, are found to form complexes with the dietary fat which are stable enough to undergo enzymatic hydrolysis by lipase. This restrains accumulation of fat in human body 21,22 . Hence, inclusion complex of SNP and cyclodextrins can be of a great deal for the weight loss/weight management dietary food supplement for sportsman or obese person 23 .

Experimental Section
Materials. Phenylephrine hydrochloride, Synephrine, α and β-cyclodextrin of puris grade of purity ≥98.0% were purchased from Sigma-Aldrich and were kept in a refrigerator as received and used right away.  Apparatus. Utilizing JASCO V-530 UV−Vis spectrophotometer, UV−visible spectra were recorded with a wavelength accuracy of ±0.5 nm. Cell temperature during the experiment was controlled from 298.15 K to 308.15 K with a digital thermostat.
Studies on surface tension at the experimental temperatures with the accuracy of ±0.1 mN m −1 were done by employing K9 digital TENSIOMETER (Krüss GmbH, Hamburg, Germany) which uses the platinum ring detachment technique. The temperature of the experimental solutions was kept constant at 298.15 K by circulating thermostat water through a double-walled glass vessel holding solution.
Proper instrumentation of METTLER-TOLEDO Seven Multi conductivity meter provides specific conductivity values with an uncertainty of ±1.0 μS m −1 . Constancy in temperature at the specific value of the solutions under experiment was maintained with an auto-thermostatic water bath. HPLC-grade water with a specific conductance of 6.0 μS m −1 was utilized for conductivity measurement. Calibration of the Systronics Type CD -30 conductivity cell was done using 0.01 M freshly prepared aqueous solution of KCl.
Fluorescence spectra were noted via JASCO V-530 UV/VIS Spectrophotometer, at 25 °C in a Hellma quartz cuvette (250-400 nm spectral range, 2.0 mL volume, 10 mm path length) equipped with a magnetic stir bar. To a solution of [SNP/PEH] (100 µM, 1 mL) in deionized water (Millipore, 18.2 MΏ. Cm) was prepared with α-CD and β-CD. (200 µM) in the stock fluorescence spectra were recorded after 1 hr of mixing time. The output range of the machine was nearly about 2 analogs (+/−10 volts).  HRMS spectra of the solid ICs were recorded on a quadrupole time-of-flight (Q-TOF) high-resolution instrument with positive-mode electrospray ionization taking the methanol solution of the solid ICs.
FTIR spectral analysis was performed on a Perkin-Elmer FTIR spectrometer in the scanning range of 4000−400 cm −1 . According to the KBr disk method the disks were made in 1:100 ratios of sample and KBr. Studies were carried out at room temperature and at a humidity of 45%.
SEM: Scanning Electron Microscope (JSM-6360) was aided to perform the analysis and obtain the data's. It also discusses about the morphological patterns and particle size of the Inclusion Complex.
Antimicrobial activity assay: In this experiment (gram negative E. coli), (gram positive B. subtilis) were considered as model organism. This test was done according to the Agar cup method. In brief, spread plate technique was applied to inoculate the organisms in Muller-Hinton agar and the compounds were applied in agar cup at 1 mg/ml concentration in separate plates and incubated at 37 °C for 24 hrs. Double distilled water was used as the control. Antimicrobial activity was determined by means of the zone of inhibition surroundings agar cup. Each of the experiments was done in triplicate.  Cell viability assay: In this experiment pure sample as well as ICs, SNP + α-CD, SNP + β-CD, PEH + α-CD and PEH + β-CD were added in the nutrient agar broth and E. coli and B. subtilis were inoculated. After 24 hrs of incubation at 37 °C cells were plated, and colony count was completed. Growth in nutrient broth without the ICs was taken as the control. All the experiment was done in triplicate. Level of significance (p) for all experiment was set to 0.05.
Procedure. All the solutions under experiment were prepared after checking the solubility of the PEH, SNP and CDs in triply distilled, deionized and degassed water. METTLER TOLEDO AG-285 analytical balance with an uncertainty of ±0.1 mg at 298.15 K was used to weigh all the experimental materials. Loss of materials caused by evaporation during mixing and working with the solutions was minimized by taking sufficient precautions. For the preparation of the solid inclusion complexes, 20 mL 1.0 (mM) solutions of α and β-CD were prepared separately with triply distilled, deionized and degassed water which, allowed to stir for 6 hours on a magnetic stirrer. Then, 20 mL 1.0 (mM) aqueous solutions of SNP/PEH were added drop wise to the previously prepared aqueous solution of α-CD or β-CD making the ultimate equimolar mixture and were continued to stir for 48 hours at 55-60 °C. The suspensions obtained after cooling the mixture to 5 °C were filtered to obtain white crystalline powder, which were then dried in air and preserved in vacuum desiccators for further use.

Result and Discussion
Job plot: Stoichiometry of inter molecular association between guest and host. The stoichiometry of the host-guest inclusion complexes was determined by employing the well-established Jobs method 24 . UV-Vis spectroscopic data were used in this technique to determine the stoichiometry of inclusion complexation. Absorption spectra of a set of solutions, prepared by mixing aqueous SNP/PEH solution with the aqueous   α-CD/β-CD solution in the sort of 0-1 mole fraction, were recorded at 298.15 K of temperature. Absorbance of the prepared set of the solutions were taken at λ max = 209 nm for SNP and λ max = 219 nm for PEH (Fig. 2). Jobs plots of (SNP + α-CD, SNP + β-CD, PEH + α-CD, PEH + β-CD) were obtained by plotting a graph, ΔA × R vs R. Where, ΔA is the deference in absorbance between the pure SNP/PEH and the solutions of the set, prepared with CDs (Tables S1-S4  Application of (K a ) to Van't Hoff equation Application of ( θ K a ) to Van't Hoff equation Application of (K a C ) to Van't Hoff equation Table 4. Thermodynamic parameters (ΔH°, ΔS°, ΔG°) calculated, using the association constants (K a , θ K a , K a C ) obtained from Benesi-Hildebrand method, nonlinear Program, program based mathematical calculation of non-linear changes in the conductivity data respectively. Mean errors in variables are as follows: than the triply distilled pure water, suggests SNP and PEH to have surface activity. This may be due to the simultaneous presence of phenyl ring as well as the -CH(OH)CH 2 NHCH 3 group to the opposite terminals of the same molecule. Study on the surface tension of diverse surface-active guest molecules with cyclodextrins strongly supports the inclusion phenomenon and the stoichiometry of the ICs [29][30][31][32] . In this work, surface tension of a fixed quantity of aqueous SNP/PEH solutions was studied at 298.15 K with the step wise addition of CDs solutions in same quantity (Tables S5, S6, Fig. 5). Whereas, according to literature as well as practically it is found that, there is an extremely slight change in the surface tension (γ) of CDs over a wide range of concentration in aqueous medium at 298. 15 K 33,34 . It signifies all the changes in the value of surface tension (γ) are associated with the SNP/ PEH. Being a surface phenomenon, more the number of surface active molecules in the surface of a solution, decreases more the surface tension of that solution. But permanent migration of surface active molecules from the surface to the bulk of the solution by means of solvation or many other stabilizing factors leads to the increase in surface tension (γ) of that solution. This is exactly the trend, what we observed in our experiment during step wise addition of CDs in the aqueous SNP/PEH solution (Tables S5, S6, Fig. 5). This is obviously; there is migration of surface active SNP/PEH molecules from the surface to the bulk of the solution by means of encapsulation of the SNP/PEH into the hydrophobic cavity of the CDs forming host-guest ICs 30,31,35 (Fig. 6). After a certain   (Table 1, Fig. 5). Accordingly, surface chemistry confers the inclusion phenomenon and appearance of sharp, single break point at the 1:1 molar concentration ratio of host and guest molecules for all the cases (SNP + α-CD, SNP + β-CD, PEH + α-CD, PEH + β-CD) establishes the 1:1 stoichiometry of host-guest ICs 36-39 .

Conductance: Molecular recognition of guest into host molecules and their stoichiometric ratio in ICs.
Conductimetric study is also another approach, which makes us able to conclude about the supramolecular host-guest interaction between the SNP/PEH and CDs and their stoichiometric ratio in the ICs 37,40 . Though, both the SNP and PEH are organic molecules, 10 (mM) aqueous solution of SNP and PEH shows appreciable conductivity. Being a hydrochloride salt, PEH shows higher conductivity than SNP for the same concentration. In the present work, conductivity of SNP and PEH were measured with the step wise increasing concentration of CDs, at three different temperatures from 298.15 K to 308.15 K with the interval of 5 K of temperature (Tables S5,  S6 and Fig. 7). It was found that, gradual increase in concentration of CDs leads to the decrease in conductivity (κ), of the aqueous SNP/PEH solutions (Fig. 7). The fruit full explanation for this observation comes through the decrease in the mobility of the conducting species in the solution due to molecular encapsulation of SNP/ PEH into the hydrophobic cavity of the CDs 36,41 (Fig. 6). Generation of a single break point in the conductivity curves after reaching a certain concentration of CDs, suggests, the molecular encapsulation of SNP/PEH into the cavity of CDs is 1:1 31,41 (Fig. 7). Corresponding concentration of SNP/PEH and CDs at the break points of the conductivity curve are listed in (Table 2). The near about equimolar concentration of SNP/PEH and CDs at the break points of the conductivity curve suggests the 1:1 stoichiometric ratio of the SNP/PEH into CDs of the ICs 36,37 (Fig. 3).
Ultraviolet Spectroscopy: The association constants (K a ) and Stability of the of the ICs. The binding ability of the guest into the host molecule and the stability of the inclusion complexes formed were explored by measuring the association constants (K a ) of the ICs. The UV-vis spectroscopic study enables us to determine the association constant (K a ) of the ICs in the solution phase 31 . Molar extinction coefficient (Δε) of SNP/PEH, depending upon the solvent polarity, should change while going from polar aqueous media to the apolar hydrophobic cavity of the CDs to form ICs 41,42 . To determine association constant (K a ), the changes in absorbances (ΔA) of SNP/PEH were measured with increasing concentration of CDs at the temperature range 298.15 K to 308.15 K. (Tables S7-S10) The λ max = 209 nm for SNP and λ max = 219 nm for PEH were considered to determine the association constant (K a ) in this case (Table 3). According to the Benesi-Hildebrand method to determine the association constant for the 1:1 host-guest inclusion complex, the double reciprocal plot was obtained by using the following equation 31,[43][44][45] .
a where, ΔA represents the difference in absorbances of PEH or SNP without CDs to the absorbances of the same with the CDs.
[DGs] refers to the concentration of the PEH and SNP. The association constants(K a ) of the inclusion complexes, listed in the (Table 3) were obtained by dividing the slope by the intercept of the plot given in the Figs S1, S2. UV-vis spectroscopic data were also used in a nonlinear program that practices the changes in absorbance of SNP/PEH due to its molecular recognition into the apolar cavities of CDs and the association constants ( θ K a ) were obtained 31,46 . There should be an equilibrium between host and the guest molecules to the formation of 1:1 ICs 47,48 .
The expression for the association constant θ K ( ) a can be obtained from the above equation as follows- where, A 0 is the absorbance of SNP/PEH molecules in the initial state, A obs denotes the absorbances of the same during the gradual addition of CDs and A refers to the final concentration of SNP/PEH molecules.
[CD] x and [DGs] x is the concentration of cyclodextrins added and SNP/PEH molecules respectively. The association constants ( θ K a ), obtained from the binding isotherm with the application non-linear program are listed in the Table 3.

Conductance: Program based mathematical calculation of non-linear changes in the conductivity data and association constants K
( ) a C . Non-linear changes in the conductivity data at the temperature ranging from 298.15 K to 308.15 K were utilized in the mathematical program and the association constants K ( ) a c for1:1 DGs-CDs ICs, listed in the Table 3 are frequently obtained [49][50][51] . The complexation reaction between DGs and CDs to produce ICs is supposed to proceed via the following chemical equilibrium

Fluorescence: Modified Benesi-Hildebrand equation and association constants. The association constants (K )
a F of the ICs in the solution phase were also determined using the spectrofluorometric data and the association constants, determined are found in good agreement with the data obtained from all the previously described methodology [52][53][54] (Table 3, Tables S11-S14). An enhancement of the intensities of the spectral lines  (Figs S3, S4). The observations enriched us with the knowledge that, a change in the molar extinction coefficient i.e. polarity of the environment, surrounding the chromophore, leaded by the encapsulation of the chromophore of the guest molecules (SNP and PEH) from the polar aqueous environment to the apolar hydrophobic cavity of the cyclodextrins. Sometimes, the enhancement in the intensities of the spectral lines are experienced due to the shielding of the excited singlet species of the chromophores from quenching and non-radiative decay with the protective microenvironment created by the hollow-circular, apolar cavity of the CDs 55,56 . The spectrofluorometric data were analyzed and run with the modified Benesi-Hildebrand equation to generate the double reciprocal plots, (Figs S5, S6) and the association constants (K a F ) of the ICs were obtained as the ratio of the intercept to slope of the plots.

UV-vis spectroscopy, Non-linear program based mathematical calculation, Fluorescence and
Conductance: Association constants and the thermodynamic parameters. According to the Van't Hoff Eq. (11) the various thermodynamic parameters for the formation of the inclusion complexes were derived from the ready available association constants (K a , θ K a , K a C ) obtained from Benesi-Hildebrand equation, nonlinear methods and non-linear changes in the conductivity data 31,37,41 . (Tables 4, S15-S20, Figs S7-S12) Calculation on the thermodynamic parameters of the formation of ICs, it is found that, both the changes in entropy and enthalpy of formation appears to be negative, suggesting an exothermic and entropy restricted rather Figure 14. (a,b,c) FTIR spectra of (a) SNP + β-CD, (b) SNP, (c) β-CD. than entropy driven process ( Table 4). The explanation on the decrease in entropy during inclusion complexation comes from the molecular association of the host and guest molecules to form inclusion complexes in the solution. Though the process is entropy restricted, the negative value of the enthalpy (Δ Δ θ H H , 0 0 , ΔH C0 ) makes the overall energy negative i.e. negative ΔG 0 and finds its spontaneity in the formation of ICs (Table 4).

H NMR and 2D ROESY NMR spectra analysis.
Inclusion of a molecule inside into the cavity of cyclodextrin consequences in the chemical shift of the interacting protons of both the guest and cyclodextrin in 1 H NMR spectra, due to their mutual shielding through space 57 . Encapsulation of aromatic guest molecule results the diamagnetic shielding of the interacting protons of cyclodextrin by the aromatic moiety 58 (Table S21). Cyclodextrin molecule has H3 and H5 hydrogens at inside of the conical cavity, specially, the H3 are located near the wider rim while H5 are positioned near the narrower rim and the other H1, H2 and H4 hydrogens are situated at the exterior of the cyclodextrin molecule (Fig. 1) 59 . In this work the molecular inclusions have been studied with the help of 1 H NMR spectra. The 1 H NMR spectra of the pure α-CD, β-CD, SNP + α-CD, SNP + β-CD, PEH + α-CD and PEH + β-CD systems are shown in Figs S13-S18 respectively, where the aromatic as well as signals of H3 and H5 protons of cyclodextrins may be observed with corresponding chemical shift (δ) values. The 1 H NMR spectra of the complexes reveal that the signals of interior H3 and H5 of α and β-CD plus that of the interacting aromatic protons of SNP/PEH showed substantial upfield shift confirming the formation of inclusion complexes 60 (Figs S13-S18).
2D ROESY NMR spectroscopy provides decisive evidence about the spatial closeness of the interacting atoms of the host and the guest by observing the intermolecular dipolar cross-correlations 61,62 . The protons which are situated within 0.4 nm in space may produce a rotating-frame NOE spectroscopy (ROESY) 63 . According to structure of α and β-CD, inclusion phenomenon inside into cyclodextrin cavity can be shown by the appearance of NOE cross-peaks between the protons of cyclodextrin and the protons of the aromatic guest identifying their spatial proximity 64,65 . To prove this, 2D ROESY spectra of the complexes of SNP and PEH with α and β-CD in D 2 O, were recorded, which shows significant correlation of aromatic protons of SNP and PEH with the H3 and H5 protons of α and β-CD, establishing the aromatic ring was encapsulated inside both the cyclodextrin cavities 66 (Figs 8-11). It may be detected that the H-6 protons of cyclodextrins were not influenced by the inclusion processes, suggesting that the SNP/PEH molecule was incorporated into the cyclodextrin cavity via the wider rim, not through the narrower rim as otherwise cross-peaks between the H6 and the guest would have been observed in the ROESY spectra 67 (Fig. 6).

HRMS Analysis of Inclusion complexes.
Mass spectroscopic study of the solid inclusion complexes of SNP/PEH with α and β-CD were done after the dissolution of the ICs in methanol. The spectra are shown in the Fig. 12 and Table S22. enlists the m/z values for the corresponding fragmentations added to the molecular ion peak. The appearance of peaks at the m/z 1140. 42  FTIR spectroscopy. Interpretation of the Infra-red spectroscopic data of the ICs as well as the pure host and guest molecules also reveals the veracity about the way by which the ICs are formed and supports the same circumstances of host -guest interaction as obtained from the 2D ROESY NMR spectroscopic study 70,71 . All the FTIR spectra of the solid inclusion complexes and the pure host and guest molecules were recorded by preparing KBr disk. The changes in the significant peak values in the IR spectra on going from the pure host and guest molecules to the inclusion complexes which are shown in the Figs 13-16, suggests the formation of ICs exploring the binding mode of the guests to the host molecules 60,72 . The IR stretching frequencies (cm −1 ) of noteworthy responsible for the corresponding chemical bonds are listed in the Table S23.    19 and 690.20 cm −1 respectively, shifting of these signals leads to the ready explanation that, the close proximity of the -C-H of the α-CD with the aromatic -C-H of the PEH as obtainable from the 2D ROESY spectra. (Fig. 15).
The respectively. Thus, FTIR spectral analysis also indorses the same as obtained from the 2D ROESY spectra (Fig. 16).
There is no sign of chemical reaction taking place while the formation of all the inclusion complexes, exemplified by the no appearance of additional signal in the IR spectra suggesting, all the shifting in signals appearing are responsible for the formation of ICs. well-known technique for analyzing the surface texture and particle size of solid materials. The surface morphological structures of α-CD, (SNP + α-CD) physical mixture, (SNP + α-CD) inclusion complex, β-CD, (SNP + β-CD) physical mixture, (SNP + β-CD) inclusion complex are shown in Fig. 17 respectively. From Fig. 17 it is clear, that the morphological structures that they are totally different from each other. Similarly, the surface morphological structures of α-CD, (PEH + α-CD) physical mixture, (PEH + α-CD) inclusion complex, β-CD, (PEH + β-CD) physical mixture, (PEH + β-CD) inclusion complex are shown in Fig. 18 respectively. From Fig. 18 it is clear, that the morphological structures that they are totally different from each other. Moreover, as the complexation by α-CD and β-CD can be viewed distinctly. This is an additional evidence about the formation of inclusion complexes of SNP/PEH with α and β-CD, may support the same evident from 2D ROESY NMR analysis.
Cytotoxic activity of the Inclusion complexes. No zone of inhibition is observed in case of both the gram-positive and gram-negative organisms 73,74 . There is no growth inhibition compared to control. These results suggest that ICs doesn't have any antimicrobial activity. So, it can be said that it is non-toxic for the cells. After that cell viability assay is finished. Here, we have found that cell viability of E. coli is 4.6% and 9% increase in presence of [SNP + β-CD, SNP + α-CD] and [PEH + β-CD, PEH + α-CD] respectively whereas the cell viability of B. subtilis is 3.2% and 8% increase in the presence of [SNP + β-CD, SNP + α-CD] and [PEH + β-CD, PEH + α-CD] correspondingly (Figs 19, 20). These consequences indicate that cell viability is positively regulated in occurrence of these ICs (Figs 19, 20). But there is very significant increase in growth when the samples are treated with (SNP + β-CD). So, this (SNP + β-CD) is more suitable for pharmaceutically active compounds. The outcome shows that both inclusion complexes have increased the capability of SNP inhibiting cell growth rather than PEH. Particularly, SNP, complexed with beta-cyclodextrin (β-CD) show the highest cytotoxic activity resting on E. coli and B. subtilis; with alpha-cyclodextrin (α-CD) the cytotoxic activity is found to be rather low.

Conclusions
The suggestion obtained from surface tension and conductometric study for the formation of 1:1 host-guest inclusion complexes of SNP and PEH with α and β-cyclodextrins are established by UV-Vis spectroscopy, spectrofluorimetry, 2D ROESY NMR spectrometry and SEM technique by analyzing surface texture of the solid inclusion complexes. The association constants obtained from all the well-established techniques dictates the stability of inclusion complexes formed and the thermodynamic parameters reveals truth about the feasibility of their formation. Removal of water molecules from the cavity of the CDs to make room for the guest molecule for accommodation while formation of inclusion complex, increases entropy of the process. Dimensional suitability being, one of the major stabilizing factor, the larger cavity size of β-CD (0.70 nm, diameter) compared to α-CD (0.56 nm, diameter), explains for the greater value of association constants and stability of the inclusion complexes formed with β-CD. The association constants, hence stability of the inclusion complexes of SNP with CDs were found more than that of the PEH. Because, -O-H group of SNP, being oriented to the para position may exert H-bonding interaction with CDs to some greater extent than that of the PEH, in which -OH group, being oriented at the meta -position can't travel the minimum distance for the formation of H-bond with the CDs. The hydrophobic and H-bonding interactions thus stabilizes the ICs. The Cytotoxicity and Cell viability also balances for non-toxic behavior of the ICs. Thus, inclusion complexes of the recently emerging two drugs, SNP and PEH (after their banned alternatives) stabilizes SNP and PEH from their chemical modification and conveys a new approach for regulatory release to the targeted site reducing overdoses.