Special Issue in Honor of Professor Nick Hadjiliadis’ RetirementView this Special Issue
Synthesis, Solution, and Structural Characterization of Tetrahydrofuranyl-2,2-Bisphosphonic Acid Disodium Salt
Bisphosphonates are biologically relevant therapeutics for bone disorders and cancer. Reaction of -chlorobutyric acid, phosphorus acid, and phosphorus trichloride without the use of solvent gave the tetrahydrofuranyl-2,2-bisphosphonate sodium salt (). The was isolated, characterized in solution by , , and NMR spectroscopy and in solid state by single X-Ray crystallography. The crystal structure showed that the forms in the crystal infinite two-dimensional sheets stacked one parallel to the other. A comparison of the chelating properties of with similar hydroxyl bisphosphonate ligands shows that the strength of the Na–O(furanyl/hydroxyl) bond is directly related to the total charge of the ligand anion.
Bisphosphonates (BPs) represent a very important class of compounds known for their medicinal and other properties [1, 2]. In contrast to the naturally occurring pyrophosphate P-O-P group, these compounds contain a characteristic P-C-P bridge, which is chemically and enzymatically nonhydrolyzable. BPs exhibit specific affinity towards bone, which makes them an excellent therapeutic for bone resorption diseases (especially osteoporosis, Paget’s disease, tumour induced osteolysis, hypercalcemia originated from malignancy) by inhibition of farnesyl diphosphate synthase (FPPS) and for bone tumour caused by metastatic breast tumours [3–5]. The biological activity of bisphosphonates has been found to be dependent on the structure , lipophilicity , and bone binding affinity  of the compounds. For example, the structure of the side chain of bisphosphonates is important in determining the potency of individual bisphosphonates in biological models, and this includes a potential role for the side chain in modulating bone binding. The presence and position of an oxygen or nitrogen atom within the side chain of bisphosphonates are directly related to their relative potency and bone binding affinity. In general, the oxygen/nitrogen containing bisphosphonates exhibits high bone binding affinities, which might be expected to reduce efficacy. However, some of these molecules are very potent and the difference in potency has been attributed to the inhibition of FPPS. Therefore, the structural characterization of bisphosphonates is very important tool for the elucidation of the mechanisms of the bone binding and biological activity and for the design of new, more potent compounds.
Herein, we report the synthesis, solution characterization, and crystal structure of the sodium salt of tetrahydrofuranyl-2,2-bisphosphonate (). In the crystal the compound forms, the unusual for bisphosphonates, two-dimensional parallel sheets arrangement. A comparison of the chelating properties of with similar hydroxyl bisphosphonate ligands shows that the Na–O(furanyl/hydroxyl) bond length increases as the charge of the ligand decreases. Thus, the weak Na–O(furanyl) interaction is attributed to the 2 charge of the ligand.
2. Materials and Methods
The chemicals obtained from commercial sources were reagents grade and used as received. Solvents were freshly distilled and water was doubly distilled. -Chlorobutyric acid, phosphorus acid, and were purchased from Sigma-Aldrich. Sodium hydroxide and phosphorus trichloride were purchased from Merck. Ethyl alcohol and methyl alcohol were bought from BDH. The synthesis was performed in inert Argon atmosphere.
2.2. Synthesis and Crystallization of Tetrahydrofuranyl-2,2-Bisphosphonic Acid Disodium Salt ()
A mixture of -chlorobutyric acid (7.12 g, 58.1 mmol), phosphorus acid (3.38 g, 41.2 mmol), and phosphorus trichloride (10.9 g, 79.4 mmol) was heated slowly to reach the temperature C over an hour under constant flow of argon. The two-phase mixture was refluxed for 4 hours at C and the viscous yellow mass was stirred vigorously for 10 hours at C. Precooled water (125 mL) was added slowly to the mixture under stirring as evolution of gaseous HCl was observed. The water (110 mL) was constantly removed by distillation from the orange reaction mass. Another 150 mL of water was added to the mixture and allowed to cool down slowly to room temperature (C). The pH of the reaction mixture was adjusted to 3 by careful addition of 4 M NaOH solution. Absolute ethyl alcohol (150 mL) was poured slowly to the resulting yellow solution in order to induce crystallisation. White solid was formed after 12 hours stirring at C and collected by filtration. The yellow filtrate was again treated with ethyl alcohol (100 mL) to obtain another batch of white solid. Total three batches were collected and for the crude product the yield was calculated to be 9.75 g, (65% based on -chlorobutyric acid). Crystallisation was performed by layering concentrated aqueous solution (pH = 4.2) of with methyl alcohol. Samples for NMR measurements were prepared by dissolution of the crystals in . The atoms have been numbered according to the numbering in Scheme 1. NMR (, ppm): 3.90 (t, = 6.7 Hz, two protons H(4)); 2.30 (m, = 6.7 Hz, = 16.0 Hz, two protons, H(2)); 1.99 (p, = 6.7 Hz, two protons H(3)). NMR (, ppm): 81.3 (t, = 144 Hz), C(1); 70.3, C(4); 30.2, C(3); 26.8 (t, = 2.26 Hz), C(2). NMR (, ppm): 21.9 (t, = 16.0 Hz). Anal. Calcd. for (f.w.: 276.02): C, 17.41; H, 2.92. Found: C, 17.38; H, 2.89.
2.3. NMR Studies
The one-dimensional H, C, and P spectra were recorded in a 300 MHz Avance Bruker spectrophotometer at 300.13 MHz for , 121.46 MHz for , and 75.47 MHz for with a 5 mm multinucleus probe at ambient temperature (C) in solution. A -pulse width and 1-second relaxation delay were applied for , , and NMR spectra.
2.4. X-Ray Crystallography
Crystal data, details of data collection, and refinement of crystal structure of are provided in Table 1. The intensity data for the compound were collected at 100 K on an XCalibur III 4-cycle diffractometer, equipped with a CCD camera detector. The experimental data were collected using Mo K radiation ( = 0.7107 Å) at a crystal-to-detector distance 60 mm. The structure was solved by direct methods using the programs SHELXS-86 and refined anisotropically (nonhydrogen atoms) by full-matrix least squares on [9, 10]. The H atoms were calculated geometrically and refined with riding model with isotropic displacement parameters. The programs ORTEP-3  and Diamond-3 were used for diagrams and WINGX  was used to prepare material for publication.
3. Results and Discussion
Compound has been synthesized by heating at C the mixture of -chlorobutyric acid, phosphorus acid, and phosphorus trichloride (Scheme 1). This molecule has been previously synthesized by Kieczykowski et al.  following a similar procedure by heating the mixture at C in methanesulfonic acid. However, the product isolated from that reaction had always been contaminated with the methanesulfonate salt and the yield of the reaction was less than 30%. In contrast, the reaction without the use of solvent gave pure product and yield more than 60%. Upon completion of the reaction, it is quenched with water. The product was crystallized by addition of ethanol in an aqueous solution of at pH 3.0. It is worth to notice that the isolated compound is the disodium salt of the bisphosphonate anion although it has been isolated at pH lower than that of the monosodium salt (pH 4.3) by Kieczykowski et al. This difference is attributed to the different methodology used for the crystallization of the product.
3.2. Solution Characterization
The molecule was characterized in aqueous solution by , , and NMR spectroscopy. The J coupling constants were used to confirm the assignments made by the NMR spectra. The NMR spectrum gave one triplet due to the coupling ( = 16.0 Hz) of phosphorus nuclei with the two H(2) protons at 21.9 ppm. This chemical shift is close to the chemical shifts of other hydroxyl bisphosphonates such as etidronate (L1, Scheme 2) (22.9 ppm) . The spectrum of showed a triplet assigned to C(1). The splitting of this peak is assigned to the strong coupling of C(1) with the two phosphorus atoms ( = 144.0 Hz). The weak coupling of C(2) with the two phosphorus nuclei ( = 2.26 Hz) gives a triplet feature for this peak that can be easily distinguished from the peak assigned to C(3). The NMR spectrum exhibits a triplet for H(4) protons (due to the coupling with the two H(3) protons), a quintuplet for H(3) (due to the coupling with the two H(2) and two H(4) protons), and a multiplet for H(2) (due to the coupling with the two H(3) protons and the two P atoms).
3.3. X-Ray Crystallographic Structure
The molecular structure of is shown in Figures 1 and 2(a). Crystallographic data are provided in Table 1. Tables 3 and 4 contain the interatomic bond lengths and angles of . The bisphosphonate anion in the structure of possesses an overall charge, thus two of the phosphonate oxygen atoms are deprotonated (P–), two have protons attached (P–OH), and two form double bond with phosphorous (P=O). Both phosphorous atoms exhibit distorted tetrahedral geometries. The P–O bond lengths are close to 1.50 Å for the P-O(deprotonated) and 1.57 Å for the P–OH(protonated). Although the P–OH bonds can be located very easily by inspection of the bond distances, the P– and P=O cannot be distinguished because of the charge delocalization over the –P=O groups. The bond angles around P(1), P(2) range from 104.0(4) to 116.4(3). The O–P–O angles involving the two unprotonated oxygen atoms are the largest ones in both phosphonate groups. The two phosphorus centers are fairly close (3.087(3) Å), due to their bisectional positions relative to C(1) on the furanyl ring (P(1)–C(1)–P(2) = 114.7(5)). The P…P distance is a little larger than the distance (3.035(3) Å) observed in the crystal structure of nitrogen containing respective molecule, the pyrrolidine-2,2-diylbisphosphonic acid .
The coordination sphere is different for each of the two sodium atoms (Figures 3(a) and 3(b)). The Na(1) exhibits a distorted octahedral environment defined by O(5), O(4), O(2), O(6), O(3), and O(). The Na(2) exhibits a capped octahedral geometry with the O(1), O(2), O(3), O(5), O(6), and O() to define an octahedron and the O(7) added to the triangular face defined by O(1), O(3), and O(6). The Na–O bond distances range between 2.327(7) and 2.872(7) Å with the Na(2)–O(7) to exhibit the longest bond distance. The sodium-centered octahedrons and capped octahedrons form an edge-sharing linear array directed along the axis (1,1,0) of the crystal (Figure 3(a)). These arrays are connected to each other with the diphosphonate anions (Figure 3(b)) forming infinite 2D sheets. The two-dimensional networks are stacked one over the other along axis (Figure 3(c)). The distance between the mean planes defined from the atoms of two successive sheets is 11.668(5) Å. This distance is actually the length of the c axis. The furanyl rings are located above and below the 2D structures insulating the sodium ions in each layer. It is possible that this unusual for diphosphonates two-dimensional structure is induced by the hydrophobic furanyl rings which prevent any electrostatic interlayer interactions.
A recent study  has shown that bisphosphonates binding strength can be calculated by summarizing the interactions of the phosphonate groups, the side chain groups, the hydroxyl group, and the hydrophobic group of the ligand with the bone surface. Most importantly, it has been found that although the group is necessary for the strong binding of phosphonates on the bone, the interaction of the with the bone surface is weak. Selected bond lengths of –O/–O of similar tetrahydrofuranyl-2,2-bisphosphonate compounds extracted from their crystal structure are shown in Table 4. The comparison shows that the chelating properties of the ligands are in good agreement with their binding properties on bone surface. Thus, the bond distances between the donor oxygen atoms of bisphosphonate ligand and the cations might give an estimation of the bisphosphonates bone binding strength. The (/)–O(phosphoryl) bond distances are similar, within the range 2.261 to 2.873 Å, suggesting that the interaction of all phosphonate groups with the bone is of similar strength for all ligands. In all crystal structures of the calcium salts of the bisphoshonates, except the structure of etidronate, the metal ion is ligated only to the phosphonate oxygen atoms. The charge of the ligand in all calcium structures is . In contrast, the ions have been found to be ligated both to the phosphonate and the hydroxyl oxygen in the structures of the sodium salts of the respective ligands. However, the –O bond length is much shorter in the structures of monoanion bisphosphonate structures (~2.5 Å) than in the respective dianions (~ 2.8 Å) including . Apparently, this bond length comparison indicates that the –O(hydroxyl/furanyl) bond strength is directly dependent on the total charge of the bisphosphonate ligand; the bond length increases by decreasing the charge of the ligand. In addition, the weak interaction of hydroxyl group with the bone  also suggests that the total charge of the ligand that binds the bone surface is .
The –O bond distances found in the crystal structure of are close to the respective bond distance of 2 (L2 = pamidronate, Scheme 2) showing that the contribution of these bonds for bone binding is similar . However, the tetrahydrofuranyl side chain does not contain any group that will contribute to the bone binding, thus it is expected that will be a weaker bone binder than , , and (Scheme 2).
The tetrahydrofuranyl-2,2-bisphosphonic acid has been prepared with an efficient method without the use of solvent, producing pure product in high yield over 60%. The crystal structure of the complex showed that the disodium salt of bisphosphonate crystallizes forming two-dimensional sheets stacked parallel one over the other.
The results of this study show that the crystallographic characterization of a widely used class of drug molecules, the bisphosphonate salts, provide important information on their biological activity (bone binding, FPPS inhibition), and can be used for the design of new more active molecules. The bisphosphonates binding strength with bone can be calculated by summarizing the interactions of the phosphonate groups, the side chain groups, the hydroxyl group, and the hydrophobic group of the ligand with the bone surface. The crystallographic data of the / salts of bisphosphonate show that the interaction of the –OH group with the metal ions is weak which is in agreement with the results from a recent study  on the interaction of the bisphosphonates with bone surface. Furthermore a comparison of / bisphosphonate structures shows that the strength of the Na–O(hydroxyl/furanyl) bond reduces with the decrease of the total bisphosphonate anion charge.
The authors thank the Research Promotion Foundation of Cyprus for the financial support of this work with the proposal YEIA/BIO/0308(BIE)/03.
W. M. Abdou and A. A. Shaddy, “The development of bisphosphonates for therapeutic uses, and bisphosphonate structure-activity consideration,” Arkivoc, vol. 2009, no. 9, pp. 143–182, 2009.View at: Google Scholar
H. Fleisch, “Bisphosphonates: mechanisms of action,” Endocrine Reviews, vol. 19, no. 1, pp. 80–100, 1998.View at: Google Scholar
J. J. Body, R. Bartl, P. Burckhardt et al., “Current use of bisphosphonates in oncology,” Journal of Clinical Oncology, vol. 16, no. 12, pp. 3890–3899, 1998.View at: Google Scholar
G. M. Sheldrick, “SHELX-86: program for the solution of crystal structures,” 1990.View at: Google Scholar
G. M. Sheldrick, “SHELXL-97, program for crystal structure refinement,” 1997.View at: Google Scholar
L. J. Farrugia, “ORTEP-3 for windows—a version of ORTEP-III with a graphical user interface (GUI),” Journal of Applied Crystallography, vol. 30, no. 5, p. 565, 1997.View at: Google Scholar
L. J. Farrugia, “WinGX suite for small-molecule single-crystal crystallography,” Journal of Applied Crystallography, vol. 32, no. 4, pp. 837–838, 1999.View at: Google Scholar
G. R. Kieczykowski, R. B. Jobson, D. G. Melillo, D. F. Reinhold, V. J. Grenda, and I. Shinkai, “Preparation of (4-amino-1-hydroxybutylidene)bisphosphonic acid sodium salt, MK-217 (alendronate sodium). An improved procedure for the preparation of 1-hydroxy-1,1-bisphosphonic acids,” Journal of Organic Chemistry, vol. 60, no. 25, pp. 8310–8312, 1995.View at: Publisher Site | Google Scholar
I. Cukrowski, L. Popoviç, W. Barnard, S. O. Paul, P. H. van Rooyen, and D. C. Liles, “Modeling and spectroscopic studies of bisphosphonate-bone interactions. The Raman, NMR and crystallographic investigations of Ca-HEDP complexes,” Bone, vol. 41, no. 4, pp. 668–678, 2007.View at: Publisher Site | Google Scholar
V. A. Uchtman, “Structural investigations of calcium binding molecules. II. The crystal and molecular structures of calcium dihydrogen ethane-1-hydroxy-1,1-diphosphonate dihydrate, ; implications for polynuclear complex formation,” The Journal of Physical Chemistry, vol. 76, no. 9, pp. 1304–1310, 1972.View at: Google Scholar
L. M. Shkol'nikova, A. A. Masyuk, and E. G. Afonin, “Structural chemistry of 1-hydroxyethylidenediphosphonic acid complexes,” Koordinatsionnaya Khimiya, vol. 16, p. 902, 1990 (Russian).View at: Google Scholar