We report the activities of 62 bisphosphonates as inhibitors of the Leishmania major mevalonate/isoprene biosynthesis pathway enzyme, farnesyl pyrophosphate synthase. The compounds investigated exhibit activities (IC(50) values) ranging from approximately 100 nM to approximately 80 microM (corresponding to K(i) values as low as 10 nM). The most active compounds were found to be zoledronate (whose single-crystal X-ray structure is reported), pyridinyl-ethane-1-hydroxy-1,1-bisphosphonates or picolyl aminomethylene bisphosphonates. However, N-alicyclic aminomethylene bisphosphonates, such as incadronate (N-cycloheptyl aminomethylene bisphosphonate), as well as aliphatic aminomethylene bisphosphonates containing short (n = 4, 5) alkyl chains, were also active, with IC(50) values in the 200-1700 nM range (corresponding to K(i) values of approximately 20-170 nM). Bisphosphonates containing longer or multiple (N,N-) alkyl substitutions were inactive, as were aromatic species lacking an o- or m-nitrogen atom in the ring, or possessing multiple halogen substitutions or a p-amino group. To put these observations on a more quantitative structural basis, we used three-dimensional quantitative structure-activity relationship techniques: comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA), to investigate which structural features correlated with high activity. Training set results (N = 62 compounds) yielded good correlations with each technique (R(2) = 0.87 and 0.88, respectively), and were further validated by using a training/test set approach. Test set results (N = 24 compounds) indicated that IC(50) values could be predicted within factors of 2.9 and 2.7 for the CoMFA and CoMSIA methods, respectively. The CoMSIA fields indicated that a positive charge in the bisphosphonate side chain and a hydrophobic feature contributed significantly to activity. Overall, these results are of general interest since they represent the first detailed quantitative structure-activity relationship study of the inhibition of an expressed farnesyl pyrophosphate synthase enzyme by bisphosphonate inhibitors and that the activity of these inhibitors can be predicted within about a factor of 3 by using 3D-QSAR techniques.

The effects of a series of 102 bisphosphonates on the inhibition of growth of Entamoeba histolytica and Plasmodium falciparum in vitro have been determined, and selected compounds were further investigated for their in vivo activity. Forty-seven compounds tested were active (IC 50 < 200 µM) versus E. histolytica growth in vitro. The most active compounds (IC 50 ∼ 4-9 µM) were nitrogen-containing bisphosphonates with relatively large aromatic side chains. Simple n-alkyl-1-hydroxy-1,1-bisphosphonates, known inhibitors of the enzyme farnesylpyrophosphate (FPP) synthase, were also active, with optimal activity being found with C9-C10 side chains. However, numerous other nitrogen-containing bisphosphonates known to be potent FPP synthase inhibitors, such as risedronate or pamidronate, had little or no activity. Several pyridine-derived bisphosphonates were quite active (IC 50 ∼ 10-20 µM), and this activity was shown to correlate with the basicity of the aromatic group, with activity decreasing with increasing pK a values. The activities of all compounds were tested versus a human nasopharyngeal carcinoma (KB) cell line to enable an estimate of the therapeutic index (TI). Five bisphosphonates were selected and then screened for their ability to delay the development of amebic liver abscess formation in an E. histolytica infected hamster model. Two compounds were found to decrease liver abscess formation at 10 mg/kg ip with little or no effect on normal liver mass. With P. falciparum, 35 compounds had IC 50 values <200 µM in an in vitro assay. The most active compounds were also simple n-alkyl-1-hydroxy-1,1-bisphosphonates, having IC 50 values around 1 µM. Five compounds were again selected for in vivo investigation in a Plasmodium berghei ANKA BALB/c mouse suppressive test. The most active compound, a C9 n-alkyl side chain containing bisphosphonate, caused an 80% reduction in parasitemia with no overt toxicity. Taken together, these results show that bisphosphonates appear to be useful lead compounds for the development of novel antiamebic and antimalarial drugs.

gammadelta T cells are the first line of defense against many infectious organisms and are also involved in tumor cell surveillance and killing. They are stimulated by a broad range of small, phosphorus-containing antigens (phosphoantigens) as well as by the bisphosphonates commonly used in bone resorption therapy, such as pamidronate and risedronate. Here, we report the activation of gammadelta T cells by a broad range of bisphosphonates and develop a pharmacophore model for gammadelta T cell activation, in addition to using a comparative molecular similarity index analysis (CoMSIA) approach to make quantitative relationships between gammadelta T cell activation by bisphosphonates and their three-dimensional structures. The CoMSIA analyses yielded R(2) values of approximately 0.8-0.9 and q(2) values of approximately 0.5-0.6 for a training set of 45 compounds. Using an external test set, the activities (IC(50) values) of 16 compounds were predicted within a factor of 4.5, on average. The CoMSIA fields consisted of approximately 40% hydrophobic, approximately 40% electrostatic, and approximately 20% steric interactions. Since bisphosphonates are known to be potent, nanomolar inhibitors of the mevalonate/isoprene pathway enzyme farnesyl pyrophosphate synthase (FPPS), we also compared the pharmacophores for gammadelta T cell activation with those for FPPS inhibition, using the Catalyst program. The pharmacophores for gammadelta T cell activation and FPPS inhibition both consisted of two negative ionizable groups, a positive charge feature and an endocyclic carbon feature, all having very similar spatial dispositions. In addition, the CoMSIA fields were quite similar to those found for FPPS inhibition by bisphosphonates. The activities of the bisphosphonates in gammadelta T cell activation were highly correlated with their activities in FPPS inhibition: R = 0.88, p = 0.002, versus a human recombinant FPPS (N = 9 compounds); R = 0.82, p < 0.0001, for an expressed Leishmania major FPPS (N = 45 compounds). The bisphosphonate gammadelta T cell activation pharmacophore differs considerably, however, from that reported previously for gammadelta T cell activation by phosphoantigens (Gossman, W.; Oldfield, E. J. Med. Chem. 2002, 45, 4868-4874), suggesting different primary targets for the two classes of compounds. The ability to quite accurately predict the activity of bisphosphonates as gammadelta T cell activators by using 3D QSAR techniques can be expected to help facilitate the design of additional bisphosphonates for potential use in immunotherapy.

We report the results of a comparative molecular field analysis (CoMFA) investigation of the growth inhibition of the bloodstream form of Trypanosoma brucei rhodesiense trypomastigotes by bisphosphonates. A quantitative three-dimensional structure-activity relationship CoMFA model for a set of 26 bisphosphonates having a range of activity spanning approximately 3 orders of magnitude (minimum IC(50) = 220 nM; maximum IC(50) = 102 microM) yielded an R(2) value of 0.87 with a cross-validated R(2) value of 0.79. The predictive utility of this approach was tested for three sets of three compounds: the average pIC(50) error was 0.23. For the nitrogen-containing bisphosphonates, in general, the activity was aromatic- >> aliphatic-containing side chains. The activity of aromatic species lacking an alkyl ring substitution decreased from ortho to meta to para substitution; halogen substitutions also reduced activity. For the aliphatic bisphosphonates, the IC(50) values decreased nearly monotonically with increasing chain length (down to IC(50) = 2.0 microM for the n-C(11) alkyl side chain species). We also show, using a "rescue" experiment, that the molecular target of the nitrogen-containing bisphosphonate, risedronate, in T. b. rhodesiense is the enzyme farnesyl pyrophosphate synthase. In addition, we report the LD(50) values of bisphosphonates in a mammalian cell general toxicity screen and present a comparison between the therapeutic indices and the IC(50) values in the T. b. rhodesiense growth inhibition assay. Several bisphosphonates were found to have large therapeutic indices (> or =200:1) as well as low IC(50) values, suggesting their further investigation as antiparasitic agents against T. b. rhodesiense.

A solid-state deuterium NMR study of localized mobility at the C9pG10 step in the DNA dodecamer [d(CGCGAATTCGCG)]2 is described. In contrast to the results of earlier deuterium NMR studies of furanose ring and backbone dynamics within the d(AATT) moiety, the furanose ring and helix backbone of dC9 display large amplitudes of motion on the 0.1 ms time scale at hydration levels characteristic of the B form structure. Solid-state deuterium NMR line shape data obtained from labeled dC9 DNA are interpreted using a composite motion model, in which the DNA oligomer is treated as rotating as a whole about the helix axis, while the base, furanose ring, and phosphodiester backbone execute localized motions. Consistent with past solid-state NMR studies, the amplitude and rate of the uniform rotation of the dC9-labeled oligomer are found to be sensitive to hydration level. Amplitudes of localized reorientational motions of C−D bonds in the furanose ring and backbone of dC9 are found to be larger than the librational amplitudes for the C−D bonds in the base of dC9, indicating that the pyrimidine base sugar does not move as a rigid entity and intersects a locally flexible region of the phosphodiester backbone. At hydration levels corresponding to 10−12 waters per nucleotide, Zeeman relaxation times for the furanose ring and backbone deuterons of dC9 in B form DNA equal 0.025 and 0.03 ms, respectively, and are the shortest relaxation times observed thus far for any deuteron in the DNA dodecamer at comparable hydration levels. The results of this solid-state NMR study suggest the existence of a significant dynamic component of sequence-specific recognition in this system.

Bisphosphonates are a class of molecules in widespread use in treating bone resorption diseases and are also of interest as immunomodulators and anti-infectives. They function by inhibiting the enzyme farnesyl diphosphate synthase (FPPS), but the details of how these molecules bind are not fully understood. Here, we report the results of a solid-state (13)C, (15)N, and (31)P magic-angle sample spinning (MAS) NMR and quantum chemical investigation of several bisphosphonates, both as pure compounds and when bound to FPPS, to provide information about side-chain and phosphonate backbone protonation states when bound to the enzyme. We then used computational docking methods (with the charges assigned by NMR) to predict how several bisphosphonates bind to FPPS. Finally, we used X-ray crystallography to determine the structures of two potent bisphosphonate inhibitors, finding good agreement with the computational results, opening up the possibility of using the combination of NMR, quantum chemistry and molecular docking to facilitate the design of other, novel prenytransferase inhibitors.

We report the design, synthesis and testing of a series of novel bisphosphonates, pyridinium-1-yl-hydroxy-bisphosphonates, based on the results of comparative molecular similarity indices analysis and pharmacophore modeling studies of farnesyl diphosphate synthase (FPPS) inhibition, human Vgamma2Vdelta2 T cell activation and bone resorption inhibition. The most potent molecules have high activity against an expressed FPPS from Leishmania major, in Dictyostelium discoideum growth inhibition, in gammadelta T cell activation and in an in vitro bone resorption assay. As such, they represent useful new leads for the discovery of new bone resorption, antiinfective and anticancer drugs.

Solid-state deuterium NMR is used to investigate perturbations of the local, internal dynamics in the EcoRI restriction binding site, -GAATTC- induced by cytidine methylation. Methylation of the cytidine base in this sequence is known to suppress hydrolysis by the EcoRI restriction enzyme. Previous solid-state deuterium NMR studies have detected large amplitude motions of the phosphate-sugar backbone at the AT-CG junction of the unmethylated DNA sequence. This study shows that methylation of the cytidine base in a CpG dinucleotide reduces the amplitudes of motions of the phosphate-sugar backbone. These observations suggest a direct link between suppression of the amplitudes of localized, internal motions of the sugar-phosphate backbone of the DNA and inhibition of restriction enzyme cleavage.