Base Dynamics in the HhaI Protein Binding Site.

Protein-DNA interactions play an important role in numerous biological functions within the living cell. In many of these interactions, the DNA helix is significantly distorted upon protein-DNA complex formation. The HhaI restriction-modification system is one such system, where the methylation target is flipped out of the helix when bound to the methyltransferase. However, the base flipping mechanism is not well understood. The dynamics of the binding site of the HhaI methyltransferase and endonuclease (underlined) within the DNA oligomer [d(G1A2T3A4G5C6G7C8T9A10T11C12)]2 are studied using deuterium solid-state NMR (SSNMR). SSNMR spectra obtained from DNAs deuterated on the base of nucleotides within and flanking the [5'-GCGC-3']2 sequence indicate that all of these positions are structurally flexible. Previously, conformational flexibility within the phosphodiester backbone and furanose ring within the target sequence has been observed and hypothesized to play a role in the distortion mechanism. However, whether that distortion was occurring through an active or passive mechanism remained unclear. These NMR data demonstrate that although the [5'-GCGC-3']2 sequence is dynamic, the target cytosine is not passively flipping out of the double-helix on the millisecond-picosecond time scale. Additionally, although previous studies have shown that both the furanose ring and phosphodiester backbone experience a change in dynamics upon methylation, which may play a role in recognition and cleavage by the endonuclease, our observations here indicate that methylation has no effect on the dynamics of the base itself.

Kinetics and thermodynamics of BI-BII interconversion altered by T:G mismatches in DNA.

T:G mismatches in DNA result in humans primarily from deamination of methylated CpG sites. They are repaired by redundant systems, such as thymine DNA glycosylase (TDG) and methyl-binding domain enzyme (MBD4), and maintenance of these sites has been implicated in epigenetic processes. The process by which these enzymes identify a canonical DNA base in the incorrect basepairing context remains a mystery. However, the conserved contacts of the repair enzymes with the DNA backbone suggests a role for protein-phosphate interaction in the recognition and repair processes. We have used 31P NMR to investigate the energetics of DNA backbone BI-BII interconversion, and for this work have focused on alterations to the activation barriers to interconversion and the effect of a mismatch compared with canonical DNA. We have found that alterations to the ΔG of interconversion for T:G basepairs are remarkably similar to U:G basepairs in the form of stepwise differences in ΔG of 1-2 kcal/mol greater than equivalent steps in unmodified DNA, suggesting a universality of this result for TDG substrates. Likewise, we see perturbations to the free energy (∼1 kcal/mol) and enthalpy (2-5 kcal/mol) of activation for the BI-BII interconversion localized to the phosphates flanking the mismatch. Overall our results strongly suggest that the perturbed backbone energetics in T:G basepairs play a significant role in the recognition process of DNA repair enzymes.

Single-Base Lesions and Mismatches Alter the Backbone Conformational Dynamics in DNA.

DNA damage has been implicated in numerous human diseases, particularly cancer, and the aging process. Single-base lesions and mismatches in DNA can be cytotoxic or mutagenic and are recognized by a DNA glycosylase during the process of base excision repair. Altered local dynamics and conformational properties in damaged DNAs have previously been suggested to assist in recognition and specificity. Herein, we use solution nuclear magnetic resonance to quantify changes in BI-BII backbone conformational dynamics due to the presence of single-base lesions in DNA, including uracil, dihydrouracil, 1,N6-ethenoadenine, and T:G mismatches. Stepwise changes to the %BII and ΔG of the BI-BII dynamic equilibrium compared to those of unmodified sequences were observed. Additionally, the equilibrium skews toward endothermicity for the phosphates nearest the lesion/mismatched base pair. Finally, the phosphates with the greatest alterations correlate with those most relevant to the repair of enzyme binding. All of these results suggest local conformational rearrangement of the DNA backbone may play a role in lesion recognition by repair enzymes.

New Insights into the Interactions of a DNA Oligonucleotide with mPEGylated-PAMAM by Circular Dichroism and Solution NMR.

Dendrimers are well-defined, highly branched, synthetic three-dimensional molecules with a large number of reactive end groups. PAMAM dendrimers form stable complexes with DNA chemistries and constitute an important class of nonviral, cationic vectors in gene delivery. The aim of this study is to examine the interactions of a 12 bp DNA oligonucletide with PAMAM-G2 and mPEG- b-PAMAM-G3 having eight surface amine groups under physiological conditions, using constant DNA concentration but varying dendrimer concentration. 1D 31P NMR, 2D NOESY, and CD spectroscopic methods were employed to investigate the interactions between the dendrimer and the DNA. The CD experiments carried out with a constant DNA concentration of 25 µM and dendrimer concentrations from 0 to 100 µM indicated minimal change to the chirality of the DNA for both types of dendrimers. While the PAMAM-G2 dendrimer caused aggregation of the majority of the DNA, the 2D NMR data of the DNA with an mPEG- b-PAMAM-G3 dendrimer indicated general broadening of the 1D 31P peaks from the DNA phosphates, a small number of 1H chemical shift perturbations (CSPs), and reduction of specific 1H-1H NOE intensities. These data suggest there is minimal structural alteration of the DNA in the complex and indicate preferential binding of the dendrimer to the central AATT region of the DNA sequence. The results herein are the first such results demonstrating a soluble DNA complex with the mPEG- b-PAMAM-G3 dendrimer analyzed by multidimensional NMR.

Solid state 2H NMR analysis of furanose ring dynamics in DNA containing uracil.

DNA damage has been implicated in numerous human diseases, particularly cancer, and the aging process. Single-base lesions, such as uracil, in DNA can be cytotoxic or mutagenic and are recognized by a DNA glycosylase during the process of base excision repair. Increased dynamic properties in lesion-containing DNAs have been suggested to assist recognition and specificity. Deuterium solid-state nuclear magnetic resonance (SSNMR) has been used to directly observe local dynamics of the furanose ring within a uracil:adenine (U:A) base pair and compared to a normal thymine:adenine (T:A) base pair. Quadrupole echo lineshapes, , and relaxation data were collected, and computer modeling was performed. The results indicate that the relaxation times are identical within the experimental error, the solid lineshapes are essentially indistinguishable above the noise level, and our lineshapes are best fit with a model that does not have significant local motions. Therefore, U:A base pair furanose rings appear to have essentially identical dynamic properties as a normal T:A base pair, and the local dynamics of the furanose ring are unlikely to be the sole arbiter for uracil recognition and specificity in U:A base pairs.

Spectroscopic analysis of interactions between alkylated silanes and alumina nanoporous membranes.

Transport across alumina nanoporous membranes can be altered via surface attachment of alkylated trimethoxysilane compounds. The mechanism of attachment has been previously assumed to be monolayer silane coverage through full chemisorption regardless of reaction conditions. This chemisorption arises via covalent Si-O-Al bond formation resulting from condensation between the three putative silanols (due to hydrolysis of the three Si-OCH(3) bonds) and hydroxides present on the alumina surface. If this model was correct, methanol would be produced in large quantities in the reaction solution, and the methoxy moieties would no longer be present on the silane molecule. The results presented in this paper utilized FT-IR and both solution and solid-state NMR to examine the chemical nature of octadecyltrimethoxysilane (ODTMS) present on the alumina surface. The FT-IR results confirm the presence of the silane on the membrane. The (1)H solution NMR results indicate small but detectable methanol production during attachment. The solid-state NMR results demonstrate that the methoxy proton NMR integrated peak intensities remain in nearly the same ratios present in the free silane, concluding that the majority of methoxy groups are intact while the silane is attached to the membrane surface. These three results suggest that monolayer surface coverage and chemisorption through full covalent bonding is not the primary means of attachment for ODTMS on the surface of alumina nanomembranes under these reaction conditions.

Furanose dynamics in the HhaI methyltransferase target DNA studied by solution and solid-state NMR relaxation.

Both solid-state and solution NMR relaxation measurements are routinely used to quantify the internal dynamics of biomolecules, but in very few cases have these two techniques been applied to the same system, and even fewer attempts have been made so far to describe the results obtained through these two methods through a common theoretical framework. We have previously collected both solution 13C and solid-state 2H relaxation measurements for multiple nuclei within the furanose rings of several nucleotides of the DNA sequence recognized by HhaI methyltransferase. The data demonstrated that the furanose rings within the GCGC recognition sequence are very flexible, with the furanose rings of the cytidine, which is the methylation target, experiencing the most extensive motions. To interpret these experimental results quantitatively, we have developed a dynamic model of furanose rings based on the analysis of solid-state 2H line shapes. The motions are modeled by treating bond reorientations as Brownian excursions within a restoring potential. By applying this model, we are able to reproduce the rates of 2H spin-lattice relaxation in the solid and 13C spin-lattice relaxation in solution using comparable restoring force constants and internal diffusion coefficients. As expected, the 13C relaxation rates in solution are less sensitive to motions that are slower than overall molecular tumbling than to the details of global molecular reorientation, but are somewhat more sensitive to motions in the immediate region of the Larmor frequency. Thus, we conclude that the local internal motions of this DNA oligomer in solution and in the hydrated solid state are virtually the same, and we validate an approach to the conjoint analysis of solution and solid-state NMR relaxation and line shapes data, with wide applicability to many biophysical problems.

Backbone dynamics in the DNA HhaI protein binding site.

The dynamics of the phosphodiester backbone in the [5'-GCGC-3'] 2 moiety of the DNA oligomer [d(G 1A 2T 3A 4 G 5 C 6 G 7 C 8T 9A 10T 11C 12)] 2 are studied using deuterium solid-state NMR (SSNMR). SSNMR spectra obtained from DNAs nonstereospecifically deuterated on the 5' methylene group of nucleotides within the [5'-GCGC-3'] 2 moiety indicated that all of these positions are structurally flexible. Previous work has shown that methylation reduces the amplitude of motion in the phosphodiester backbone and furanose ring of the same DNA, and our observations indicate that methylation perturbs backbone dynamics through not only a loss of mobility but also a change of direction of motion. These NMR data indicate that the [5'-GCGC-3'] 2 moiety is dynamic, with the largest amplitude motions occurring nearest the methylation site. The change of orientation of this moiety in DNA upon methylation may make the molecule less amenable to binding to the HhaI endonuclease.

Solid-state nuclear magnetic resonance spectroscopy studies of furanose ring dynamics in the DNA HhaI binding site.

The dynamics of the furanose rings in the GCGC moiety of the DNA oligomer [d(G 1A 2T 3A 4 G 5 C 6 G 7 C 8T 9A 10T 11C 12)] 2 are studied by using deuterium solid-state NMR (SSNMR). SSNMR spectra obtained from DNAs selectively deuterated on the furanose rings of nucleotides within the 5'-GCGC-3' moiety indicated that all of these positions are structurally flexible. The furanose ring within the deoxycytidine that is the methylation target displays the largest-amplitude structural changes according to the observed deuterium NMR line shapes, whereas the furanose rings of nucleotides more remote from the methylation site have less-mobile furanose rings (i.e., with puckering amplitudes < 0.3 A). Previous work has shown that methylation reduces the amplitude of motion in the phosphodiester backbone of the same DNA, and our observations indicate that methylation perturbs backbone dynamics through the furanose ring. These NMR data indicate that the 5'-GCGC-3' is dynamic, with the largest-amplitude motions occurring nearest the methylation site. The inherent flexibility of this moiety in DNA makes the molecule more amenable to the large-amplitude structural rearrangements that must occur when the DNA binds to the HhaI methyltransferase.

Contrasting views of the internal dynamics of the HhaI methyltransferase target DNA reported by solution and solid-state NMR spectroscopy.

Solution and solid-state NMR have been used conjointly to probe the internal motions of a DNA dodecamer containing the recognition site for the HhaI methyltransferase. The results strongly suggest that ns-mus motions contribute to the functionally relevant dynamic properties of nucleic acids during DNA methylation.

Solid-state NMR, crystallographic, and computational investigation of bisphosphonates and farnesyl diphosphate synthase-bisphosphonate complexes.

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.

Pyridinium-1-yl bisphosphonates are potent inhibitors of farnesyl diphosphate synthase and bone resorption.

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.

Selective in vitro effects of the farnesyl pyrophosphate synthase inhibitor risedronate on Trypanosoma cruzi.

We present the results of the first detailed study of the molecular and cellular basis of the antiproliferative effects of the bisphosphonate risedronate (Ris) on Trypanosoma cruzi, the causative agent of Chagas' disease. Ris and related compounds, which block poly-isoprenoid biosynthesis at the level of farnesyl pyrophosphate synthase, are currently used for the treatment of bone resorption disorders, but also display selective activity against trypanosomatid and apicomplexan parasites. Ris induced a dose-dependent effect on growth of the extracellular epimastigote form of T. cruzi; complete growth arrest and cell lysis ensued at 150 microM. Growth inhibition was associated with depletion of the parasite's endogenous sterols, but complete growth arrest and loss of cell viability took place before full depletion of these compounds, suggesting that disappearance of other essential poly-isoprenoids is involved in its anti-parasitic action. Ris had a variety of effects on cellular ultrastructure, including mitochondrial swelling, disorganisation of other organelles, such as reservosomes and the kinetoplast, together with the appearance of autophagic vesicles and progressive vacuolization of the cytoplasm. Ris had selective antiproliferative effects against the clinically relevant amastigote form of T. cruzi, and at 100 microM, was able to prevent completely the development of T. cruzi infection of murine muscle heart or Vero cells, and to cure cultures which were already infected. Ris induced drastic ultrastructural alterations in the intracellular parasites and blocked amastigote to trypomastigote differentiation, with no biochemical or ultrastructural effects on the host cells, which fully recovered their normal structure and activity after treatment. Ris is, therefore, a promising lead compound for the development of new drugs against T. cruzi.

Quantitative structure-activity relationships for gammadelta T cell activation by bisphosphonates.

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.

Effects of bisphosphonates on the growth of Entamoeba histolytica and Plasmodium species in vitro and in vivo.

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 microM) versus E. histolytica growth in vitro. The most active compounds (IC(50) approximately 4-9 microM) 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) approximately 10-20 microM), 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 microM 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 microM. 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.

3-D QSAR investigations of the inhibition of Leishmania major farnesyl pyrophosphate synthase by bisphosphonates.

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.

Activity of bisphosphonates against Trypanosoma brucei rhodesiense.

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.