Mass spectrometry and molecular modeling studies on the inclusion complexes between alendronate and β-cyclodextrin

Complexation of alendronate sodium (AlnNa) with β-cyclodextrin (β-CD) was studied by means of ESI-mass spectrometry. The experimental results show that stable 1:1 inclusion complexes between selected bisphosphonates and β-CD were formed. In addition, complexes with different stoichiometry were observed. DFT/B3LYP calculations were performed to elucidate the different inclusion behavior between alendronate and β-CD. Molecular modeling showed that the inclusion complex of Aln-β-CD where the two phosphonate groups bound to the central carbon atom of bisphosphonate were inserted into the cavity of β-CD from its “top” side was thermodynamically more favorable than when they were inserted from its “bottom” side; the complexation energy was −74.05 versus −60.85 kcal/mol. The calculations indicated that the formation of conventional hydrogen bonds was the main factor for non-covalent β-CD:Aln complex formation and stabilization in the gas phase.


Introduction
Alendronate sodium (AlnNa, Fig. 1) is a drug from the bisphosphonate group, which is commonly used to treat diseases that result from increased bone resorption. It is used in postmenopausal, senile and steroid-induced osteoporosis as well as for treatment of the Paget's bone disease and bone metastases [1,2]. Bisphosphonates are characterized by two phosphonate groups bound to the central carbon atom. They are analogues of pyrophosphate in which the central oxygen atom is replaced by a carbon atom, resulting in a P-C-P group. Bisphosphonates are characterized by a very low bioavailability [3]. After oral intake, they are absorbed through the gastrointestinal tract at a low rate. Furthermore, orally administered bisphosphonates may cause irritation and damage to the esophageal mucosa [4,5]. One possibility to improve the bioavailability of the drug substance and to reduce its harmful effects is to produce complexes with macrocyclic oligosaccharides, for example b-cyclodextrin.
b-Cyclodextrin (b-CD, Fig. 1) is a macrocyclic oligosaccharide consisting of seven D-(?)glucopyranose units linked by a-1,4-glycosidic bonds. b-CD has an internal cavity with a truncated cone shape. It is about 8 Å deep and 6.0-6.5 Å wide. The 14 secondary hydroxyls on the C-2 and C-3 carbon atoms of the glucopyranose units are situated at the wider rim and the seven primary hydroxyls on the C-6 carbons are situated on the narrower rim of the cone. Most hydrogen atoms are directed to the middle part of the cone, resulting in a rather hydrophobic cavity along with edges of a more hydrophilic character [6,7].
The complexation of drugs with cyclodextrins has been widely used in the pharmaceutical industry. The first pharmaceutical product in the world, prostaglandin E2/b-CD (Prostarmon E TM sublingual tablets), was marketed in Japan in 1976 by Ono Pharmaceutical Co. Twelve years later, the first European cyclodextrin-containing pharmaceutical product, piroxicam/b-CD (Brexin Ò tablets), was marketed and in 1997, the first US-approved product, itraconazole/2-hydroxypropyl-b-CD oral solution (Sporanox Ò ) was introduced [8,9]. Nowadays 30-40 different drugs are marketed worldwide in different types of cyclodextrin complex formulations, including different types of tablets (i.e., conventional, chewing and sublingual tablets), oral capsules, parenteral solutions, suppositories, nasal sprays, eye drop solutions and dermal products [7]. The inclusion of drug molecules into CDs leads to important modifications to the physical, chemical and biological properties of the molecules. Complexation can potentially enhance the solubility of insoluble drugs, improve the efficacy, chemical stability and bioavailability of poorly soluble drugs, reduce their toxicity by making the drugs effective at lower doses, and control their rate of release [10,11]. It can also moderate drug odors and flavors [12].
The ability of b-CD to form inclusion complexes with drug substances is mainly due to hydrophobic interactions, van der Waals and electrostatic interactions, as well as the formation of hydrogen bonds [13,14]. Due to the lack of covalent bonds in these complexes, CDs are treated as excipients, and the resulting complexes are not considered as new active substances [15]. b-CD is mentioned in the US Pharmacopoeia (USP/NF), European Pharmacopoeia (Ph. Eur.) and also in the Japanese Pharmaceutical Codex (JPC). In 1995, the Joint FAO/WHO Expert Committee on Food Additives (JECFA) announced the ADI (Acceptable Daily Intake) value for b-CD at 0-5 mg/bw kg. The US Federal Register approved b-CD as GRAS (Generally Recognized as Safe) in 1997 [16,17].
Many physicochemical methods have been effectively used to characterize the stoichiometric ratio of the inclusion complexes between b-cyclodextrin and a guest, including UV spectroscopy, fluorescence measurements, mass spectrometry, nuclear magnetic resonance (NMR), etc. [18][19][20][21]. Electrospray ionization mass spectrometry (ESI-MS) is known for its high accuracy and sensitivity in mass assignment. Therefore, it is particularly suitable to obtain insight into the formation of complex species and their stoichiometry [22][23][24]. The main aim of our investigation was to examine whether alendronate sodium can form inclusion complexes with b-cyclodextrin, which could possibly be used to improve the bioavailability of alendronate sodium and diminish its side effects. The investigation was performed by way of ESI-MS experiments in order to see whether ions derived from a complex could be formed after introducing a b-CD and sodium alendronate mixture into the mass spectrometer, and by way of molecular modeling.

Chemicals and reagents
Alendronate sodium was obtained from Cipla Ltd. (India). b-CD was purchased from Cyclolab Ltd. (Hungary) and used without further treatment. All reagents were HPLC grade and were used as such. Acetonitrile and methanol were from Labscan and formic acid (FA) from Park Scientific. Double distilled water was used throughout.

Mass spectrometry measurements
A complex of alendronate with b-CD was prepared by the addition of an equimolar portion of the sodium alendronate to an aliquot of b-CD in aqueous solution (5 mmol/l). Then the whole solution was shaken for 3 min to mix it effectively before being introduced to the mass spectrometer. All ESI-MS experiments were performed on a MicrOTOF-QII mass spectrometer equipped with an electrospray ion source (Bruker-Daltonics, Bremen, Germany). The following settings were used: electrospray ionization (ESI) in positive and negative ion modes. The dry gas (nitrogen) flow rate was set at 8.0 l/min and the dry heater operated at 180°C. The capillary voltage was set at 4,500 V and the end plate offset at -500 V. Collision energy was varied in the range of 25-30 eV. MS data were recorded in full scan mode (from 50 to 3,000 m/z). Data processing was carried out with Chromeleon 6.8. Samples were introduced directly to the electrospray source of the MS using an LC pump and the mobile phase at a flow rate 0.

Molecular modeling
Calculations were carried out by Spartan'08 software [25] using the DFT/B3LYP method and the 6-31 ? G* basis set. The initial structure of alendronate was optimized with the use of the software. The geometry of b-CD was obtained from its crystallographic structure as determined by X-ray diffraction, which is archived in the Cambridge databank and then optimized at the B3LYP/6-31 ? G* level.
Modeling of the complex was performed by docking the optimized structure of the alendronate molecule into the b-CD cavity and allowing for full geometry optimization. Two different inclusion orientations were considered: one with two phosphonate groups bound to the central carbon atom of bisphosphonate inserted into the cavity of bCD from its ''top'' side (orientation 1) and, another, from its ''bottom'' side (orientation 2).
DE of the complex formation was calculated according to Eq. (1): In this equation E complex , E b-CD and E Aln are the total energies of the system and of the two separated subsystems, respectively.

Mass spectrometry measurements
The positive ESI mass spectrum of b-cyclodextrin and alendronate sodium mixture is reported in Fig. 2  The negative ESI mass spectrum of b-cyclodextrin and alendronate sodium mixture is reported in Fig. 3 (very low but with the correct isotopic profile). The singly charged ion at m/z 1,179.340 was not identified. The total ion current of positive ions was about 20 times higher than that of negative ions. Tandem mass spectrometry (MS/MS) was used to confirm the structures. The negative ion tandem mass spectrum of the intense doubly charged ion at m/z 690.682, attributed to [bCD ? AlnH-2H] 2 These results show that the Aln/b-CD complex is clearly indicated in the negative ion spectrum where it is registered as [bCD ? AlnH-2H] 2with a relatively strong peak at m/ z 690.682. Negative ions are more likely to occur because a dissociated anion of the sodium salt is present in the solution.
The results indicate not only the 1:1 Aln/b-CD complex, but also complexes with different stoichiometry that had formed in gas phase and could possibly be formed in solution; however, their corresponding peaks were less intense.

Molecular modeling studies
Calculations using the DFT/B3LYP method were performed in order to obtain some global information about the geometry of the host-guest complexes and to find the intermolecular interaction leading to b-CD and alendronate inclusion complexation. The formed hydrogen bonds, the calculated energies of alendronate and b-CD, the energy of the complex and the complex formation energy (DE) are shown in Table 1   After optimization, the alendronate molecule was not completely entrapped in b-CD. The complexation energies were negative, which demonstrates that the inclusion processes of the drug into b-CD were thermodynamically favorable. The lowest values of the complexation energy correspond to the most stable complexes. It is assumed that when the distance between the H atom to a donor (O) and an acceptor (O) is in the range of 0.8-2.8 Å , and the O-H-O angle is above 120.0°, a hydrogen bond is formed [26]. According to this criterion, we noted that the optimized structure in orientation 1 (Fig. 4a) reveals five intermolecular hydrogen bonds formed between the hydrogen atoms of the hydroxyl groups of the drug and the oxygen atoms of the hydroxyl groups of the glucopyranose units of b-CD, between the -OH groups of b-CD and the oxygen atoms of the phosphonate groups, and between the hydroxyl group of b-CD and the nitrogen atom of the drug. The optimized structure with the minimum energy of inclusion complex, i.e. orientation 2 (Fig. 4b), shows that a single intermolecular H bond was formed between a hydrogen atom of the hydroxyl group of the drug and an oxygen atom of the hydroxyl group of b-CD. These results confirm that hydrogen bonding interactions play an important role in the inclusion complexation process.

Conclusions
This experimental and theoretical study constitutes a physicochemical characterization of the interaction between b-CD and alendronate, a drug from the bisphosphonate group. The specific host-guest complexations between alendronate and b-CD were easily detected by ESI-MS in the gas phase, which suggests that hydrogen bonds were maintained upon vaporization and ionization. Furthermore, molecular modeling showed that the inclusion complex of Aln-b-CD where the two phosphonate groups bound to the central carbon atom of bisphosphonate were inserted into the cavity of b-CD from its ''top'' side was thermodynamically more favorable than when they Fig. 4 Energy-minimum structures obtained by DFT/ B3LYP calculations for the alendronate/b-CD complex in two different orientations: a orientation 1, b orientation 2. Dotted lines represent hydrogen bonds were inserted from its ''bottom'' side. The difference between the two complex formation energies was 13.2 kcal/ mol. It can be expected that the present findings will promote extensive future studies into the encapsulation of bisphosphonates to reduce toxicity to the esophageal mucosa.