Molecular energy of ?G = -12.65 KJ

Molecular recognition of flunarizine dihydrochloride and ?-cyclodextrin inclusion complex by NMR and computational approachAbstractWe have carried out a comprehensive structural study on ?-cyclodextrin (?-CD) and flunarizine dihydrochloride (FLN) interaction using high resolution NMR and computational molecular docking approach. We report here the complete structural assignments of FLN using 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 Ka=157 M-1 and change in Gibbs free energy of ?G = -12.65 KJ mol-1. We have compared binding constant of ?-CD inclusion complex with two nearly similar structures, FLN and cetirizine. Two-dimensional 1H-1H ROESY spectral data and molecular docking studies equally showed the modes of penetration of the aromatic rings from the wider rim side into ?-CD cavity. The possible geometrical structures of the ?-CD-FLN inclusion complex have been proposed in which, aromatic rings protrude close to narrower rim of ?-CD truncated cone. The splitting of aromatic signals of FLN in the presence of ?-CD suggest chiral differentiation of the guest FLN by ?-CD.   Key wordsFlunarizine dihydrochloride, ?-cyclodextrin, inclusion complex, NMR spectroscopy, Molecular recognition, ROESY, molecular docking1. IntroductionThe human 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 (Figure 1) is a di-fluorinated derivative of cinnarizine and is 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. Cyclodextrins (CDs) or cycloamyloses are truncated cone shaped macrocycles produced from starch through enzymatic degradation. CDs is 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 (Figure 1) 5, 6. The glucose residue in CDs have 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 large number of hydroxyl groups but the central cavities are relatively hydrophobic (Figure 1). 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. 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 outer surface, 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-28, microencapsulation 29, nanotechnologies 30-33, remediation 34, decontamination 35 and biotechnology 36. The unique properties of CDs allow their various applications in many areas 37-39.  CDs are used to prepare inclusion complexes with pharmaceuticals for biomedical applications and biomedicine 22, 31, 36-38.  CDs are widely used in food industry as food additives, stabilizing flavors, 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-43. Upon complexation with pharmaceutical compounds, CDs form inclusion complexes with 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 are 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, proteins 44, nanoparticles 45, and add value addition of taste and colour of tooth paste 46.  Among various known spectroscopic methods such as Ultraviolet-visible (UV-Vis), Fourier-transform infrared spectroscopy (FTIR) for the studies of inclusion complexes of 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-49. NMR spectroscopic titration data can be used to determine the stoichiometry and association constant of the host-guest complexes 50-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. 2D 1H-1H Rotating-frame Overhauser Effect SpectroscopY (ROESY) has been found to be very 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-52. 2D 1H-1H ROESY provides useful information about the location and depth of inclusion of guest molecule into CD cavity 47, 50-52. The formation of an inclusion complex of a guest molecule with CDs results in the 1H chemical shift changes (??) in both the host and guest protons. 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 guests molecules 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 approach thus elucidating the structure of the ?-CD-FLN inclusion complex.   2. Materials and Methods2.1 MaterialsFlunarizine dihydrochloride (FLN) was a kind gift from by Geno Pharmaceutical Ltd. India. ?-cyclodextrin (?-CD) was obtained from Geertrui Haest, Cerestar Application Centre, Food & Pharma Specialities, France. These materials were used as obtained. 2.2 NMR Spectroscopy1H 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 300K during all NMR experiments. The mixing time (?mix) for 2D 1H-1H ROESY spectra was 500 ms under the spin lock 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: ?? = ?(complex) – ?(free)2.3 Molecular docking studiesMolecular docking studies were performed using Autodock Vina 1.1.2 53. Three dimensional coordinates of ?-CD (PDB Id: 1DMB) were sourced from http://www.rcsb.org/ while FLN were 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 keeping 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. 3. Results and discussion 3.1 2D 1H-1H ROESY spectrum of ?-CD-FLN mixture and structure of the ?-CD-FLN inclusion complex by NMR and molecular docking approachThe understanding of host-guest supramolecular structure is important for pharmaceutical industry for development of drug-cyclodextrin 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 mixture of ?-CD and FLN was analyzed. The NMR spectroscopic studies and assignments of ?-CD and FLN protons are discussed in section 3.2. The 2D 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 it should have been. 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 Figure 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 Figure S1 (Supplementary Information). On the basis of 1:1 stoichiometry of the ?-CD-FLN inclusion complex (See Section 3.4) and 2D 1H-1H ROESY spectral data, it can be inferred that F-containing aromatic ring preferentially enters the ?-CD cavity to form the inclusion complex. Also, the non-observance of cross peak between H-3′ and H-1 (Figure 2B) suggested 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 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 is 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 provides us not only the mode of inclusion but also depth of penetration inside the ?-CD cavity during complexation process. The best docked model of ?-CD: FLN complex is shown in Figure 3. It is evident that the mode of penetration of FLN guest in to 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. Cetirizine dihydrochloride (CTZ), an antihistamine drug used to treat allergies formed inclusion complexes in which the penetration of CTZ into the ?-CD cavity was from wider rim side 56. Our 2D 1H-1H ROESY and molecular docking approach together provide information about the penetration of FLN from the wider rim ride 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 into complexation (Figure 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 (Figure 4B, F, H). The docking binding affinity for best energy minimized ?-CD:FLN complex was obtained to be -5.4 kcal/mol, which is favourable for such type of complex 50. The ensemble of all possible computational docking conformations of ?-CD:FLN complex is shown in Figure S2 (Supplementary Information). Based on molecular docking studies performed, it is apparent that all aromatic rings dock into the ?-CD:FLN inclusion complexation but F-containing aromatic ring participates more favourably than the phenyl ring 50.3.2 1H NMR spectral assignments and chemical shift change data of ?-CDThe assignment of the ?-CD protons, in the spectra of mixture of ?-CD and FLN, 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 a FLN- ?-CD mixture showing ?-CD regions are shown in Figure S3 (Supplementary Information). On 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-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 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 the H-3′ and H-5′ of ?-CD in the 1H NMR spectra of mixtures of ?-CD-FLN thus confirm the formation of inclusion complex between ?-CD and FLN analogous to previous studies 47, 50-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 Figure 5 and their ?? data are listed in Table 1.3.3 1H NMR spectral assignments and chemical shift change data of FLNThe resonance assignment of guest FLN aromatic ring protons in the free as well as host ?-CD bound state was 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 Figure 6. 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 benzene ring. For instance, ortho-, meta-, and para-proton signals of fluorobenzene appear at 6.99 ppm, 7.24 ppm 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) is 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 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 Figure S4 (Supplementary Information). The 1H NMR and 2D 1H-1H COSY data of ?-CD and FLN in bound and unbound states provide us the information about their assignments, shielding and deshielding effects and splitting pattern. 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 are required. Two-dimensional 1H-1H ROESY and molecular docking studies further applied to understand ?-CD-FLN inclusion complex structure (See Section 3.1). 3.4 Stoichiometry and association constant of ?-CD-FLN complexNext, 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 was established with the help of the Scott’s method 57. In Scott’s equation,FLNt/??obs = FLNt/??c + 1/Ka??c (1)where FLNt is the molar concentration of the guest, ??obs the observed chemical shift change for a given FLNt concentration, ??c the chemical shift change between a pure sample of complex and the free component at the saturation. The plot of chemical shift changes (??) for the ?-CD protons against FLN in the form of FLN/??obs versus FLN appeared to be linear fits (Figure 7) suggesting 1:1 stoichiometry for the ?-CD-FLN inclusion complex. The slope of the plot (Figure 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. We were also interested to probe the differences between binding constants of two nearly similar structures and complex formation mode. FLN has some structural similarity with cetirizine dihydrochloride (CTZ), which is an antihistamine drug used to treat allergy symptoms. The binding constant of CTZ was reported earlier to be 70 M-1 56, which is nearly half of the binding constant calculated for FLN. This could be explained due to the presence of extra phenyl and F-substituted ring. The ?G associated during ?-CD and FLN inclusion complex was calculated using standard equation (2):?G = – RTlnKa (2)where R is the universal gas constant (J mol-1 K-1), T is temperature (Kelvin) and Ka is the binding constant. The ?G value was calculated to be -12.65 KJ mol-1. 4. Conclusions The 1H NMR spectral data of pure FLN, pure ?-CD and mixtures of ?-CD and FLN in D2O confirmed the complexation between the ?-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 studies. The 2D 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 pharmaceutical industry.AcknowledgementsSKU is currently funded (DST INSPIRE Faculty Award) by the Department of Science & Technology, Govt. of India.Competing financial interestsNoneFigure and Table LegendsFigure 1. Structural representation of (A) FLN (guest) ligand and (B) ?-CD (host) macrocycle 51. Figure 2. Partial 2D 1H-1H ROESY (500 MHz) spectra of mixture of ?-CD and FLN showing interactions between FLN and ?-CD cavity protons (?mix = 500 ms).Figure 3. Computational best molecular docked conformation model of ?-CD:FLN inclusion complexes 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. Also see Figure 4 and Supplementary Figure S1 for other docked conformations obtained during docking. ?-CD shown as ball and stick with surface while FLN 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 (http://www.cgl.ucsf.edu/chimera). Figure 4: Different binding conformations obtained during FLN docked into ?-CD cavity. The docking affinity. ?-CD shown as ball and stick with surface while FLN 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 (http://www.cgl.ucsf.edu/chimera). 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.Figure 6. Part of the 2D 1H-1H COSY spectrum (500 MHz) of a mixture of ?-CD and FLN, displaying through bond interaction of aromatic protons of FLN.Figure 7. Scott’s plot showing 1:1 stoichiometry for the ?-CD-FLN complex.Table 1. 1H NMR (500 MHz) chemical shift change (??) data for the ?-CD protons in the presence of FLN.FLN/?-CD H-1′ H-2′ H-3′ H-4′ H-5′ H-6’0.2 -0.056 -0.001 -0.105 -0.022 -0.117 -0.0090.4 -0.067 -0.003 -0.157 -0.033 -0.184 -0.0190.7 -0.078 -0.005 -0.203 -0.043 -0.239 -0.0451.0 -0.084 -0.008 -0.227 -0.058 -0.275 -0.0521.3 -0.086 -0.011 -0.244 -0.060 -0.294 -0.0531.8 -0.088 -0.007 -0.258 -0.062 -0.319 -0.057Negative values indicate upfield shift changes