Skip to main content
  • Research Article
  • Open access
  • Published:

Molecular recognition of flunarizine dihydrochloride and β-cyclodextrin inclusion complex by NMR and computational approaches



Flunarizine dihydrochloride (FLN) is used in the prophylactic treatment of migraine, vertigo, occlusive peripheral vascular disease and epilepsy. Cyclodextrins (CDs) are chiral, truncated cone shaped macrocycles known for their inner hydrophobic and outer hydrophilic site. They form complexes with hydrophobic drug molecules and enhance the solubility and bioavailability of such compounds by enhancing drug permeability through mucosal tissues. NMR spectroscopy and computational docking have been recognized as an important tool for the interaction study of CDs-drug inclusion complexes in solution state.


The structural assignments of FLN and β-CD protons were determined by 1H NMR and 2D 1H-1H COSY NMR spectroscopy. 1H NMR spectroscopic studies of FLN, β-CD and their mixtures confirmed the formation of β-CD-FLN inclusion complex in solution. 1H NMR titration data for β-CD-FLN inclusion complex showed 1:1 stoichiometry, an association constant of K a  = 157 M−1 and change in Gibbs free energy of ∆G = − 12.65 kJ mol−1. The binding constant of the β-CD inclusion complex with two nearly similar structures, FLN and cetirizine dihydrochloride, were compared. Two-dimensional 1H-1H ROESY spectral data and molecular docking studies showed the modes of penetration of the aromatic rings from the wider rim side into the β-CD cavity. The possible geometrical structures of the β-CD-FLN inclusion complex have been proposed in which aromatic rings protrude close to the narrower rim of the β-CD truncated cone.


NMR spectroscopic studies of FLN, β-CD and FLN:β-CD mixtures confirmed the formation of 1:1 inclusion complex in solution at room temperature. Two-dimensional 1H-1H ROESY together with molecular docking study confirmed that the F-substituted aromatic ring of FLN penetrates into β-CD truncated cone and the tail of aromatic rings were proximal to narrower rim of β-CD. The splitting of aromatic signals of FLN in the presence of β-CD suggests chiral differentiation of the guest FLN by β-CD.


Migraine is a severe headache often unilateral, commonly accompanied by nausea, vomiting, and extreme sensitivity to sound and light. Flunarizine dihydrochloride (FLN) is a large hydrophobic fluorinated piperazine derivative, used in the prophylactic treatment of migraine, vertigo, occlusive peripheral vascular disease and epilepsy [1]. FLN (Fig. 1a) is a di-fluorinated derivative of cinnarizine and a poorly water-soluble drug. FLN is a selective calcium entry blocker with calmodulin binding properties and histamine H1 blocking activity. It is also known to prevent hepatitis C virus membrane fusion in a genotype-dependent manner [2] and to suppress endothelial angiopoietin-2 in a calcium-dependent fashion in sepsis [3]. FLN is reportedly effective against hepatitis C virus activity, preferably for the genotype 2 viruses [4].

Fig. 1
figure 1

Structural representation of a FLN (guest) ligand and b β-CD (host) macrocycle [Source Adopted from “NMR and molecular modelling studies on the interaction of fluconazole with β-cyclodextrin” by S.K. Upadhyay et al. (2009) Chemistry Central Journal 3:9]

Cyclodextrins (CDs) or cycloamyloses are truncated cone-shaped macrocycles (Fig. 1b) produced from starch through enzymatic degradation. CDs are a family of cyclic oligosaccharides and have been studied extensively as supramolecular hosts [5, 6]. The three common CDs are crystalline, homogeneous, nonhygroscopic substances, consisting of six (α-), seven (β-), and eight (γ-) d-glucose units, respectively, linked by α-d-(1→4) glycosidic bonds (Fig. 1b) [5, 6]. The glucose residue in CD has 4C1 (chair) conformation [5]. The primary hydroxyl groups (n) are located at the narrower rim whereas the wider rim is lined with secondary hydroxyl groups (2n). The outer surfaces of the CDs are highly hydrophilic due to the presence of a large number of hydroxyl groups but the central cavities are relatively hydrophobic (Fig. 1b). The outer dimension of these three common CDs are constant at 0.78 nm but their inner dimensions are variable, being 0.57 nm for α-, 0.78 nm for β-, and 0.95 nm for γ-CD respectively [6]. The H-3′ and H-5′ protons of these CDs are located in the hydrophobic central cavity whereas other protons (H-1′, H-2′, H-4′ and H-6′) are located at the outer surface (Fig. 1b), which is relatively hydrophilic. These properties facilitate their aqueous solubility and ability to encapsulate hydrophobic moieties within their central cavities through non-covalent interactions. CDs form host–guest inclusion complexes upon penetration of guest molecule in the central cavity of host CDs.

CDs are extensively studied in various areas of chemistry including macrocyclic [7], supramolecular [8, 9], agro [10], click [11], analytical [12], chromatography [13, 14], sugar-based surfactants [15], foods [16], catalysis [17, 18], membranes [19], textiles [20], cosmetics [21, 22], fragrance and aromas [23, 24], enzyme technology [25], pharmacy and medicine [26,27,28], microencapsulation [29], nanotechnologies [30,31,32,33], remediation [34], decontamination [35] and biotechnology [36]. The unique properties of CDs allow their various applications in many areas [37,38,39,40]. CDs are used to prepare inclusion complexes with pharmaceuticals for biomedical applications and biomedicine [22, 31, 36,37,38]. CDs are widely used in food industry as food additives, stabilizing flavours, to remove undesirable compounds such as cholesterol, and also as agents to avoid microbiological contaminations in the food [16]. CDs can be used to enhance solubility, bioavailability and stability of pharmaceuticals [41,42,43]. Upon complexation with pharmaceutical compounds, CDs form inclusion complexes with the ability to alter the physiochemical properties of the complexed drug. Various drugs such as nimesulide, omeprazole, piroxicam, mitomycin, diclofenac sodium, indomethacin and others complexed with CDs are approved and available in the market [42]. Inclusion complexes with dimethyl-β-CD are used in the preparation of vaccine Deptacel (Sanofi Group, Pasteur) for protection against diphtheria, tetanus and pertussis. CDs are also used to stabilize sensitive substances to light or oxygen [44], proteins [45], nanoparticles [46], and add value addition of taste and colour of toothpaste [44].

Among various known spectroscopic methods such as Ultraviolet–visible (UV–Vis), Fourier-transform infrared (FTIR) spectroscopy for the studies of inclusion complexes between host CDs and guest molecules, Nuclear Magnetic Resonance (NMR) spectroscopy is considered as one of the most significant analytical tool for understanding the interaction between host and guest molecules [47]. This technique provides not only the structural assignments of host and guest molecules but also data on the inclusion complex formation. Further NMR spectroscopy could also offer valuable information on chiral recognition or chiral discrimination or both [47,48,49]. NMR titration data can be used to determine the stoichiometry and association constant of the host–guest complexes [50,51,52].

Two-dimensional (2D) NMR method such as 1H-1H COSY (COrrelation SpectroscopY) is a useful technique, which provides information on the 1H signals arising from neighbouring protons connected through bonds and protons signals emerging from up to 4 bonds can be captured. Two-dimensional 1H-1H Rotating-frame Overhauser Effect SpectroscopY (ROESY) has been found to be useful for the investigation of the interaction between CD and guest molecule as the Nuclear Overhausser Effect (NOE) cross-peaks are observed between the protons that are close in space even if they are not bonded [47, 50,51,52]. Two-dimensional 1H-1H ROESY provides useful information about the location and depth of inclusion of guest molecule into CD cavity [47, 50,51,52].

The formation of inclusion complex of a guest molecule with CDs results in the 1H chemical shift changes (∆δ) in both the host and guest protons. The inclusion of a molecule inside the hydrophobic cavity of CD is mainly characterized by the chemical shift variation of the CD protons located inside the central cavity (H-3′ and H-5′), whereas other CD protons (H-1′, H-2′, H-4′ and H-6′) are less affected. During host–guest inclusion complex formation the guest molecule protons generally show downfield chemical shift changes but sometimes upfield chemical shift changes are also observed [47].

These analytical procedures revealing the structural details of complexes are used in pharmaceutical industries for characterization. In order to understand correct inclusion architecture of interaction between guest FLN and host β-CD, we report here a high-resolution NMR spectroscopic and computer-based molecular docking study. We describe our results based on the 1H NMR spectral data with chemical shift changes, 2D 1H-1H COSY spectrum for assignment of protons and 1H-1H ROESY spectrum together with molecular docking approaches thus elucidating the structure of the β-CD-FLN inclusion complex.

Materials and methods


Flunarizine dihydrochloride (FLN) was a kind gift from Geno Pharmaceutical Ltd. India. β-cyclodextrin (β-CD) was obtained from Geertrui Haest, Cerestar Application Centre, Food & Pharma Specialities, France. These materials were used as obtained.

NMR spectroscopy

1H NMR and 2D 1H-1H NMR (COSY, ROESY) spectra were recorded on a JEOL α-500 MHz instrument in D2O. The sample temperature was maintained at 300 K during all NMR experiments. The mixing time (τmix) for 2D 1H-1H ROESY spectra was 500 ms under the spinlock condition using standard 1H-1H ROESY pulse sequences. The chemical shift values (δ) are reported in ppm. No external indicator was used and HDO peak at 4.80 ppm was considered as an internal reference throughout this work. 1H NMR spectra of six samples of mixtures of β-CD and FLN with FLN/β-CD molar ratios ranging from 0.2 to 1.8 were recorded. The FLN/β-CD molar ratios were calculated by direct NMR integration of their appropriate signals. The concentration of β-CD was kept constant at 10 mM while that of FLN was varied from 2.0 to 18.0 mM. Chemical shifts changes (Δδ) were calculated according to the formula: \( \Delta\delta \, = \,\delta_{{({\text{complex}})}} {-} \,\delta_{{({\text{free}})}} \)

Molecular docking studies

Molecular docking studies were performed using Autodock Vina 1.1.2 [53]. Three-dimensional coordinates of β-CD (PDB Id: 1DMB) were sourced from while FLN was sourced from the UCSF ZINC database (ZINC19360739) [54]. Molecular docking of FLN into β-CD cavity was carried out following the methods as reported previously [50]. The grid centre of docking coordinates were x = − 6.89 Å, y = − 7.65 Å and z = 4.34 Å. The grid dimensions were 54 Å , 56 Å and 44 Å in x, y and z-axes respectively.

Results and discussion

2D 1H-1H ROESY spectrum of β-CD-FLN mixture and structure of the β-CD-FLN inclusion complex by NMR and molecular docking approaches

The understanding of host–guest supramolecular structure is important for the pharmaceutical industry for development of drug-CD based new formulations. In order to clearly establish the identity of the aromatic ring involved in complexation between β-CD and FLN, a 2D 1H-1H ROESY spectrum of the mixture of β-CD and FLN was analyzed. The NMR spectroscopic studies and assignments of β-CD and FLN protons are discussed in “1H NMR spectral assignments and chemical shift change data of β-CD” and “1H NMR spectral assignments and chemical shift change data of FLN” sections. Two-dimensional 1H-1H ROESY spectrum exhibited strong cross-correlation peak between the cavity protons of β-CD and the protons of the F-substituted aromatic ring of FLN thereby confirming the penetration of F-substituted aromatic rings into the β-CD cavity. The cross peaks between phenyl ring protons and β-CD cavity protons were also observed but these were relatively weak. It is apparent from the 2D 1H-1H ROESY spectrum that H-1 exhibited cross peak with only H-5′ while H-2 displayed cross peaks with both the H-3′ and H-5′ protons. The quality of 2D 1H-1H ROESY spectrum is not as good as required. Expansions of the parts of the 2D 1H-1H ROESY spectrum showing cross peaks between protons of aromatic rings of FLN and β-CD cavity protons are shown in Fig. 2. The full 2D 1H-1H ROESY spectrum displaying protons of β-CD and FLN and their NOE cross-correlation peaks close to proposed interaction site are presented in Additional file 1: Figure S1.

Fig. 2
figure 2

a and b Partial 2D 1H-1H ROESY (500 MHz) spectra of the mixture of β-CD and FLN showing interactions between FLN and β-CD cavity protons (τmix = 500 ms)

On the basis of 1:1 stoichiometry of the β-CD-FLN inclusion complex (see “Stoichiometry and association constant of β-CD-FLN complex” section) and 2D 1H-1H ROESY spectral data, it can be inferred that F-containing aromatic ring preferentially enters into the β-CD cavity to form the inclusion complex. Also, the non-observance of the cross peak between H-3′ and H-1 (Fig. 2b) suggested the position of H-1 towards narrower rim side. The penetration from wider rim side would have brought H-1 closer to H-3′ also. It appears that there are interactions between the phenyl ring and β-CD but the amount is lower compared to complex formed involving F-containing aromatic ring. The penetration of FLN into β-CD cavity was reported to be from wider rim side based on 2D 1H-1H ROESY results [55] without clear inclusion architecture. The plausible mode of inclusion and structure of the β-CD-FLN inclusion complex cannot be achieved only from 2D 1H-1H ROESY spectral data and therefore, another approach was required. In order to understand the β-CD-FLN inclusion complex structure, computer-based molecular docking was performed using Autodock Vina 1.1.2 [53]. Molecular docking studies provide us not only the mode of inclusion but also the depth of penetration inside the β-CD cavity during complexation process. The best-docked model of β-CD:FLN complex is shown in Fig. 3. It is evident that the mode of penetration of FLN guest into the β-CD cavity was from the wider rim side and similar to 2D 1H-1H ROESY results the F-containing aromatic ring participates more favourably than phenyl ring. We compared our result with cetirizine dihydrochloride (CTZ) which has some structural similarity with FLN. CTZ, an antihistamine drug used to treat allergies, formed 1:1 inclusion complexes in which the penetration of CTZ into the β-CD cavity was from wider rim side [56]. Similarly our 2D 1H-1H ROESY and molecular docking approach together provide information about the penetration of FLN from the wider rim side of the β-CD cavity [56]. Moreover, F-containing aromatic ring of FLN positioned towards the narrower rim of the β-CD truncated cone, which is also observed in 2D 1H-1H ROESY spectrum containing the cross peak between H-1, H-2 and H-5′. In the other dockings conformations models, it is apparent that the phenyl ring also participates in complexation (Fig. 4). Interestingly, similar to 2D 1H-1H ROESY results, the phenyl ring protrudes on the opening of the narrower rim side of the β-CD cavity (Fig. 4b, f, h). The docking binding affinity for the best energy minimized β-CD:FLN complex was obtained to be − 5.4 kcal mol−1, which is favourable for such type of complex [50]. The ensemble of all possible computational docking conformations of β-CD:FLN complex is shown in Additional file 1: Figure S2. Based on molecular docking studies performed, it is apparent that all aromatic rings dock into the β-CD cavity but F-containing aromatic ring participates more favourably than the phenyl ring.

Fig. 3
figure 3

Computational best molecular docked conformation model of β-CD:FLN inclusion complex performed by Autodoc Vina 1.1.2 [53] showing penetration of F-substituted aromatic ring into β-CD cavity from wider rim side. The docking affinity was obtained to be − 5.4 kcal mol−1. Also, see Fig. 4 and Additional file 1: Figure S2 for other docked conformations obtained during docking. β-CD is shown as ball and stick with the surface while FLN is shown as stick bond. All atoms are shown in their elemental colour. Non-polar hydrogens are not shown for the sake of clarity. The figure was prepared using Chimera (

Fig. 4
figure 4

ai Different binding conformations obtained during FLN docked into β-CD cavity. The docking affinity is shown under the model. β-CD is shown as ball and stick with the surface while FLN is shown as stick bond. All atoms are shown in their elemental colour. Non-polar hydrogens are not shown for the sake of clarity. The figure was prepared using Chimera (

1H NMR spectral assignments and chemical shift change data of β-CD

The assignment of the β-CD protons, in the spectra of β-CD and FLN mixture, was made with the help of their 1H signals and 2D 1H-1H COSY spectral data [47]. Expansion of 2D 1H-1H COSY spectrum of an FLN:β-CD mixture showing β-CD regions are shown in Additional file 1: Figure S3. On the investigation of 1H NMR spectra of mixtures of β-CD and FLN, an upfield shift in H-3′ and H-5′ (located inside the central cavity) signals of β-CD was observed [47, 50,51,52, 56]. Other β-CD signals (H-1′, 2′, 4′, 6′) also exhibited shift changes but these were negligible compared to H-3′ and H-5′. In the presence of FLN, ∆δ for H-5′ were more pronounced than those of H-3′ signal of β-CD.

The upfield shift of 1H signals located inside the cavity, namely H-3′ and H-5′, have been attributed to the magnetic anisotropy effect in the β-CD cavity due to the inclusion of groups rich in π-electrons [51]. The continuous upfield shift changes of 1H signals observed in H-3′ and H-5′ of β-CD in the 1H NMR spectra of β-CD-FLN mixtures thus confirm the formation of the inclusion complex between β-CD and FLN [47, 50,51,52, 56]. Expansions of part of 1H NMR spectra of pure β-CD and mixture of β-CD and FLN in varying amounts of FLN are displayed in Fig. 5 and their ∆δ data are listed in Table 1.

Fig. 5
figure 5

A part of 1H NMR spectra (500 MHz) showing protons of β-CD in the absence, as well as in the presence, of varying amount of FLN

Table 1 1H NMR (500 MHz) chemical shift change (Δδ) data for the β-CD protons in the presence of FLN

1H NMR spectral assignments and chemical shift change data of FLN

The resonance assignment of guest FLN aromatic ring protons in the free as well as host β-CD bound state were achieved using 1H NMR as well as 2D 1H-1H COSY spectral data. Part of the 2D 1H-1H COSY spectrum of the mixture of β-CD and FLN displaying through bond cross connection peaks between aromatic protons of FLN is shown in Fig. 6.

Fig. 6
figure 6

Part of the 2D 1H-1H COSY spectrum (500 MHz) of a mixture of β-CD and FLN, displaying through the bond interaction of aromatic protons of FLN

The aromatic protons were observed as three signals, a triplet at 7.20 ppm integrating for four protons, a multiplet at 7.45 ppm for three protons and a multiplet at 7.58 ppm for six protons. Fluorine has a slight ‘donor substituent’ effect in the benzene ring. For instance, ortho-, meta-, and para-proton signals of fluorobenzene appear at 6.99, 7.24 and 7.08 ppm, respectively. In styrene, ortho, meta and para protons are increasingly shielded. In order to elucidate the question on H-1 and H-2 assignment, we examined for intramolecular NOE cross peaks between H-2 and H-3 in 1H-1H ROESY spectrum. The triplet at 7.20 ppm (J = 8.6 Hz) was assigned to H-1 protons and it showed the 1H-1H COSY interaction with the multiplet at 7.58 ppm, which was ascribed to H-9 and H-2 protons. The observed shape of H-1 and H-2 (like triplets) are undoubtedly from 1H-19F cross coupling interactions. It is well known that, for fluorobenzene derivatives, the coupling constants 3J(H, F) = 6.2–10.1 Hz and 4J(H, F) = 6.2–8.3 Hz. The multiplet at 7.45 ppm was due to H-10, 11 protons. In FLN:β-CD mixtures, the signal for H-2 and H-9 separated and the nature of H-2 resembles a triplet. A doublet at 6.97 ppm (J = 16.0 Hz), which appeared in the aromatic region was ascribed to H-8, while the H-7 was found resonating as a merged doublet of the triplet at 6.32 ppm.

The aromatic protons of FLN were deshielded and pattern of their 1H NMR peaks splitting in presence of β-CD suggests some chiral differentiation of guest FLN by host β-CD [48, 49]. The 1H NMR signal for H-9, 2 which appeared as a merged signal in the spectrum of unbound FLN, separated in the spectra of some FLN:β-CD mixtures. The 1H NMR spectra of expanded aromatic regions of guest FLN in the bound as well as unbound with host β-CD is shown in Additional file 1: Figure S4. The 1H NMR shielding and deshielding pattern of β-CD and FLN protons in the bound state indicate the involvement of aromatic ring group in complexation [47] but the identity of the aromatic rings penetrating into β-CD cavity could not be achieved and therefore further studies were required. Two-dimensional 1H-1H ROESY and molecular docking studies further applied to understand β-CD-FLN inclusion complex structure (see “2D 1H-1H ROESY spectrum of β-CD-FLN mixture and structure of the β-CD-FLN inclusion complex by NMR and molecular docking approach” section).

Stoichiometry and association constant of β-CD-FLN complex

Next, we wanted to determine the stoichiometry, association constant (Ka) and the Gibb’s free energy (∆G) of the β-CD-FLN inclusion complex. The stoichiometry and Ka of the β-CD-FLN complex were established with the help of the Scott’s method [57]. In Scott’s equation,

$$ \left[ {\text{FLN}} \right]_{\text{t}} /\Delta \delta_{\text{obs}} = \, \left[ {\text{FLN}} \right]_{\text{t}} /\Delta \delta_{\text{c}} + { 1}/K_{\text{a}} \Delta \delta_{\text{c}} $$

where [FLN]t is the molar concentration of the guest, Δδobs the observed chemical shift change for a given [FLN]t concentration, Δδc the chemical shift change between a pure sample of complex and the free component at the saturation. The plot of Δδ for the β-CD protons (H-3′ and H-5′) against [FLN] in the form of [FLN]/Δδobs versus [FLN] appeared to be linear fits (Fig. 7) suggesting 1:1 stoichiometry for the β-CD-FLN inclusion complex. The slope of the plot (Fig. 7) is thus equal to 1/Δδc and the intercept with the vertical axis to 1/KaΔδc allowing the estimation of Ka to be 157 M−1, which is the average of two Ka.

Fig. 7
figure 7

Scott’s plot showing 1:1 stoichiometry for the β-CD-FLN complex

We were also interested to probe the differences between binding constants of two nearly similar structures. The binding constant of β-CD-CTZ complex was reported earlier to be 70 M−1 [56], which is nearly half of the binding constant calculated for β-CD-FLN complex. This could be due to the structural differences between CTZ and FLN. The ∆G associated during β-CD and FLN inclusion complex was calculated using standard Eq. (2):

$$ \Delta {\text{G }} = \, {-}{\text{ RTln}}K_{a} $$

where R is the universal gas constant (J mol−1 K−1), T is temperature (Kelvin) and K a is the binding constant. The ∆G value was calculated to be − 12.65 kJ mol−1.


The 1H NMR spectral data of pure FLN, pure β-CD and mixtures of β-CD and FLN in D2O confirmed the complexation between β-CD and FLN. The 1H NMR together with 2D 1H-1H COSY spectral data provided the resonance assignment of host and guest molecules. The stoichiometry, association constant and the Gibb’s free energy were determined using 1H NMR titration data. Two-dimensional 1H-1H ROESY spectral data together with computational molecular docking simulation studies confirmed that F-substituted aromatic ring of guest penetrates into the host β-CD cavity from the wider rim side. The tail end aromatic rings of guest FLN were proximal near to narrower rim side of truncated host β-CD cone. The splitting of the most of the aromatic ring protons of the FLN, in the presence of β-CD, suggests some chiral differentiation of guest FLN by host β-CD. The structural studies of FLN-β-CD inclusion complex may open new avenues for new drug formulation in the pharmaceutical industry.


  1. Holmes B, Brogden RN, Heel RC, Speight TM, Avery GS (1984) Flunarizine. A review of its pharmacodynamic and pharmacokinetic properties and therapeutic use. Drugs 27:6–44.

    Article  CAS  Google Scholar 

  2. Perin PM, Haid S, Brown RJP, Doerrbecker J, Schulze K, Zeilinger C, von Schaewen M, Heller B, Vercauteren K, Luxenburger E, Baktash YM, Vondran FWR, Speerstra S, Awadh A, Mukhtarov F, Schang LM, Kirschning A, Müller R, Guzman CA, Kaderali L, Randall G, Meuleman P, Ploss A, Pietschmann T (2016) Flunarizine prevents hepatitis C virus membrane fusion in a genotype-dependent manner by targeting the potential fusion peptide within E1. Hepatology 63:49–62.

    Article  CAS  Google Scholar 

  3. Retzlaff J, Thamm K, Ghosh CC, Ziegler W, Haller H, Parikh SM, David S (2017) Flunarizine suppresses endothelial angiopoietin-2 in a calcium-dependent fashion in sepsis. Sci Rep 7:44113.

    Article  Google Scholar 

  4. Prigozhin DM, Modis Y (2016) Flunarizine arrests hepatitis C virus membrane fusion. Hepatology 63:14–16.

    Article  Google Scholar 

  5. Szejtli J (1998) Introduction and general overview of cyclodextrin chemistry. Chem Rev 98:1743–1754.

    Article  CAS  Google Scholar 

  6. Crini G (2014) Review: a history of cyclodextrins. Chem Rev 114:10940–10975.

    Article  CAS  Google Scholar 

  7. Steed JW, Atwood JL (2009) Molecular guests in solution. In: Supramolecular chemistry, 2nd ed.; Wiley, London, Chapter 6, p 307

  8. Harada A, Takashima Y, Yamaguchi H (2009) Cyclodextrin-based supramolecular polymers. Chem Soc Rev 38:875–882.

    Article  CAS  Google Scholar 

  9. Hapiot F, Tilloy S, Monflier E (2006) Cyclodextrins as supramolecular hosts for organometallic complexes. Chem Rev 106:767–781.

    Article  CAS  Google Scholar 

  10. Morillo E (2006) Application of cyclodextrins in agrochemistry. In: Dodziuk H (ed) Cyclodextrins and their complexes. Chemistry. Analytical methods, applications, Chapter 16. Wiley-VCH, Verlag GmbH & Co. KGaA, Weinheim, p 459

    Google Scholar 

  11. Faugeras P-A, Boëns B, Elchinger P-H, Brouillette F, Montplaisir D, Zerrouki R, Lucas R (2012) When cyclodextrins meet click chemistry. Eur J Org Chem 2012:4087–4105.

    Article  CAS  Google Scholar 

  12. Li S, Purdy WC (1992) Cyclodextrins and their applications in analytical chemistry. Chem Rev 92:1457–1470.

    Article  CAS  Google Scholar 

  13. Hinze WL (1981) Applications of cyclodextrins in chromatographic separations and purification methods. Sep Purif Methods 10:159–237.

    Article  CAS  Google Scholar 

  14. Xiao Y, Ng S-C, Tan TTY, Wang Y (2012) Recent development of cyclodextrin chiral stationary phases and their applications in chromatography. J Chromatogr A 1269:52–68.

    Article  CAS  Google Scholar 

  15. Dodziuk H (2006) Cyclodextrins and their complexes. Chemistry. Analytical methods, applications. Wiley-VCH, Verlag GmbH & Co. KGaA, Weinheim

    Book  Google Scholar 

  16. Astray G, Gonzalez-Barreiro C, Mejuto JC, Rial-Otero R, Simal-Gándara J (2009) A review on the use of cyclodextrins in foods. Food Hydrocolloids 23:1631–1640.

    Article  CAS  Google Scholar 

  17. Komiyama M, Monflier E (2006) Cyclodextrin catalysis. In: Dodziuk H (ed) Cyclodextrins and their complexes. Chemistry, analytical methods, applications, Chapter 4. Wiley-VCH, Verlag GmbH & Co. KGaA, Weinheim, p 93

    Chapter  Google Scholar 

  18. Bai CC, Tian BR, Zhao T, Huang Q, Wang ZZ (2017) Cyclodextrin-catalyzed organic synthesis: reactions, mechanisms, and applications. Molecules 22:1475.

    Article  Google Scholar 

  19. Trotta F (2016) Cyclodextrin in membranes. In: Drioli E, Giorno L (eds) Encyclopedia of membranes. Springer, Berlin.

    Google Scholar 

  20. Buschmann H-J, Knittel D, Schollmeyer E (2001) New textile applications of cyclodextrins. J Incl Phenom Macrocycl Chem 40:169–172.

    Article  CAS  Google Scholar 

  21. Buschmann HJ, Schollmeyer E (2002) Applications of cyclodextrins in cosmetic products: a review. J Cosmet Sci 53:185–191

    CAS  Google Scholar 

  22. Bilensoy E (ed) (2011) Cyclodextrins in pharmaceutics, cosmetics, and biomedicine: current and future industrial applications. Wiley, Hoboken, USA

    Google Scholar 

  23. Marques HMC (2011) A review on cyclodextrin encapsulation of essential oils and volatiles. Flavour Fragr J 25:313–326.

    Article  Google Scholar 

  24. Kant A, Linforth RS, Hort J, Taylor AJ (2004) Effect of β-cyclodextrin on aroma release and flavor perception. J Agric Food Chem 52:2028–2035.

    Article  CAS  Google Scholar 

  25. Villalonga R, Cao R, Fragoso A (2007) Supramolecular chemistry of cyclodextrins in enzyme technology. Chem Rev 107:3088–3116.

    Article  CAS  Google Scholar 

  26. Zhang J, Ma PX (2013) Cyclodextrin-based supramolecular systems for drug delivery: recent progress and future perspective. Adv Drug Deliv Rev 65:1215–1233.

    Article  CAS  Google Scholar 

  27. Alvarez-Lorenzo C, García-González CA, Concheiro A (2017) Cyclodextrins as versatile building blocks for regenerative medicine. J Control Release 268:269–281.

    Article  CAS  Google Scholar 

  28. van de Manakker F, Vermonden T, van Nostrum CF, Hennink WE (2009) Cyclodextrin-based polymeric materials: synthesis, properties, and pharmaceutical/biomedical applications. Biomacromolecules 10:3157–3175.

    Article  Google Scholar 

  29. Guo R, Wilson LD (2013) Cyclodextrin-based microcapsule materials—their preparation and physiochemical properties. Curr Org Chem 17:14–21.

    Article  CAS  Google Scholar 

  30. Varan G, Varan C, Erdoğar N, Hıncal AA, Bilensoy E (2017) Amphiphilic cyclodextrin nanoparticles. Int J Pharm 531:457–469.

    Article  CAS  Google Scholar 

  31. Mejia-Ariza R, Graña-Suárez L, Verboom W, Huskens J (2017) Cyclodextrin-based supramolecular nanoparticles for biomedical applications. J Mater Chem B 5:36–52.

    Article  CAS  Google Scholar 

  32. Adeoye O, Cabral-Marques H (2017) Cyclodextrin nanosystems in oral drug delivery: a mini review. Int J Pharm 531:521–531.

    Article  CAS  Google Scholar 

  33. Trotta F, Zanetti M, Cavalli R (2012) Cyclodextrin-based nanosponges as drug carriers. Beilstein J Org Chem 8:2091–2099.

    Article  CAS  Google Scholar 

  34. Landy D, Mallard I, Ponchel A, Monflier E, Fourmentin S (2012) Cyclodextrins for remediation technologies. In: Lichtfouse E, Schwarzbauer J, Robert D (eds) Environmental chemistry for a sustainable world. Springer, Dordrecht.

    Google Scholar 

  35. Morin-Crini N, Crini G (2013) Environmental applications of water-insoluble β-cyclodextrin–epichlorohydrin polymers. Prog Polym Sci 38:344–368.

    Article  CAS  Google Scholar 

  36. Szejtli J (1990) The cyclodextrins and their applications in biotechnology. Carbohydr Polym 12:375–392.

    Article  CAS  Google Scholar 

  37. Amiri S, Amiri S (2017) Cyclodextrins: properties and industrial applications. Wiley, Chichester

    Book  Google Scholar 

  38. Sliwa W, Girek T (eds) (2017) Cyclodextrins: properties and applications. Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim

    Google Scholar 

  39. Hedges AR (1998) Industrial applications of cyclodextrins. Chem Rev 98:2035–2044.

    Article  CAS  Google Scholar 

  40. Heath JR (2015) Nanotechnologies for biomedical science and translational medicine. Proc Natl Acad Sci USA 112:14436–14443.

    Article  CAS  Google Scholar 

  41. Uekama K, Hirayama F, Irie T (1998) Cyclodextrin drug carrier systems. Chem Rev 98:2045–2076.

    Article  CAS  Google Scholar 

  42. Davis ME, Brewster ME (2004) Cyclodextrin-based pharmaceutics: past, present and future. Nat Rev Drug Discov 3:1023–1035.

    Article  CAS  Google Scholar 

  43. Loftsson T, Duchêne D (2007) Cyclodextrins and their pharmaceutical applications. Int J Pharm 329:1–11.

    Article  CAS  Google Scholar 

  44. Del Valle EEM (2004) Cyclodextrins and their uses: a review. Process Biochem 39:1033–1046.

    Article  Google Scholar 

  45. Serno T, Geidobler R, Winter G (2011) Protein stabilization by cyclodextrins in the liquid and dried state. Adv Drug Deliv Rev 3:1086–1106.

    Article  Google Scholar 

  46. Noël S, Léger B, Ponchel A, Philippot K, Denicourt-Nowicki A, Roucoux A, Monflier E (2014) Cyclodextrin-based systems for the stabilization of metallic(0) nanoparticles and their versatile applications in catalysis. Catal Today 235:20–32.

    Article  Google Scholar 

  47. Schneider H-J, Hacket F, Rüdiger V, Ikeda H (1998) NMR studies of cyclodextrins and cyclodextrin complexes. Chem Rev 98:1755–1786.

    Article  CAS  Google Scholar 

  48. Kitae T, Nakayama T, Kano K (1998) Chiral recognition of α-amino acids by charged cyclodextrins through cooperative effects of Coulomb interaction and inclusion. J Chem Soc Perkin Trans 2:207–212. doi:

    Article  Google Scholar 

  49. Dodziuk H, Koźmiński W, Ejchart A (2004) NMR studies of chiral recognition by cyclodextrins. Chirality 16:90–105.

    Article  CAS  Google Scholar 

  50. Maheshwari A, Saraswat H, Upadhyay SK (2017) Structural insights into the inclusion complexes between clomiphene citrate and β-cyclodextrin: the mechanism of preferential isomeric selection. Chirality 29:451–457.

    Article  CAS  Google Scholar 

  51. Upadhyay SK, Kumar G (2009) NMR and molecular modelling studies on the interaction of fluconazole with β-cyclodextrin. Chem Cent J 3:9.

    Article  Google Scholar 

  52. Ali SM, Upadhyay SK (2008) Complexation study of midazolam hydrochloride with β-cyclodextrin: NMR spectroscopic study in solution. Magn Reson Chem 46:676–679.

    Article  CAS  Google Scholar 

  53. Trott O, Olson AJ (2010) AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization and multithreading. J Comput Chem 31:455–461.

    CAS  Google Scholar 

  54. Irwin JJ, Sterling T, Mysinger MM, Bolstad ES, Coleman RG (2012) ZINC: a free tool to discover chemistry for biology. J Chem Inf Model 52:1757–1768.

    Article  CAS  Google Scholar 

  55. Maheshwari A, Sharma D (2010) A comparative study of inclusion complexes of flunarizine with alpha (α-CD) and beta-cyclodextrin (β-CD). J Incl Phenom Macrocycl Chem 68:453–459.

    Article  CAS  Google Scholar 

  56. Ali SM, Upadhyay SK, Maheshwari A (2007) NMR spectroscopic study of inclusion complexes of cetirizine dihydrochloride and β-cyclodextrin in solution. Spectroscopy 21:177–182.

    Article  CAS  Google Scholar 

  57. Scott RL (1956) Some comments on the Benesi–Hildebrand equation. Recl Trav Chim Pays-Bas 75:787–789.

    Article  CAS  Google Scholar 

Download references

Authors’ Contributions

SKU and SMA conceived the study and designed the experiments. SKU analysed the NMR and computational docking data. SKU and SMA wrote the manuscript. Both authors read and approved the final manuscript.


SKU is currently funded (DST INSPIRE Faculty Award) by the Department of Science & Technology, Govt. of India. SKU thanks to Dr. Souvik Maiti for continuous support and Dr. S. Ramachandran for critical reading of the manuscript.

Competing interests

The authors declare that they have no competing interests.

Ethics approval and consent to participate

Not applicable.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Authors and Affiliations


Corresponding authors

Correspondence to Santosh Kumar Upadhyay or Syed Mashhood Ali.

Additional file

Additional file 1: Figure S1:

Full 2D 1H-1H ROESY spectrum (500 MHz) showing through space cross correlation peaks between β-CD protons and aromatic rings proton of FLN. Figure S2. Ensemble of different conformations obtained during FLN docked into β-CD cavity. The best docking conformation model shown as bond while other docking conformation modes are shown as wire frame. All atoms are shown in their elemental colour (β-CD and best docking conformer). The inclusion of guest FLN was from the wide rime side during all docking simulations performed. The phenyl ring of FLN are close to the narrow rim in some docking conformations. β-CD shown as ball and stick with surface while FLN shown as stick/wire bond. Non-polar hydrogens are not shown for sake of clarity. The figure was prepared using Chimera ( Figure S3. Expanded region of 2D 1H-1H COSY spectrum (500 MHz) of FLN:β-CD mixture showing β-CD region. The assignments of β-CD protons namely H-1′, H-2′, H-3′, H-4′, H-5′ and H-6′ was made with the help of 1H NMR and 2D 1H-1H COSY cross correlation peaks. Figure S4. Expansion of part of 1H NMR spectra (500 MHz) of FLN protons in the presence, as well as in the absence, of β-CD.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Upadhyay, S.K., Ali, S.M. Molecular recognition of flunarizine dihydrochloride and β-cyclodextrin inclusion complex by NMR and computational approaches. Chemistry Central Journal 12, 33 (2018).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: