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.     

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.     

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.     

Molecular dynamics simulations of ligand-induced backbone conformational changes in the binding site of the periplasmic lysine-, arginine-, ornithine-binding protein

The periplasmic lysine-, arginine-, ornithine-binding protein (LAOBP) traps its ligands by a large hinge bending movement between two globular domains. The overall geometry of the binding site remains largely unchanged between the open (unliganded) and closed (liganded) forms, with only a small number of residues exhibiting limited movement of their side chains. However, in the case of the ornithine-bound structure, the backbone peptide bond between Asp11 and Thr12 undergoes a large rotation. Molecular dynamics simulations have been used to investigate the origin and mechanism of this backbone movement. Simulations allowing flexibility of a limited region and of the whole binding site, with and without bound ligands, suggest that this conformational change is induced by the binding of ornithine, leading to the stabilisation of an energetically favourable alternative conformation.     

Influence of DNA structure on adjacent site cooperative binding

Previous NMR studies of Hoechst 33258 with the d(CTTTTGCAAAAG)2 sequence have shown very strong (K2 > K1) cooperativity between two adjacent binding sites (Searle, M. S.; Embrey, K. J. Nucleic Acids Res. 1990, 18 (13), 3753- 3762). In contrast, surface plasmon resonance (SPR) results with the hairpin analog of the same sequence show significantly reduced cooperativity. In an effort to explain the difference, two-dimensional (2-D) NMR experiments were done on both duplex and hairpin. Hoechst 33258 and an amidine analog, DB183, show very strong cooperativity with the duplex DNA but much weaker cooperativity with the hairpin. The significantly lower thermal melting temperature (Tm) of the duplex (34.8 degrees C) in comparison to its hairpin analog (62.3 degrees C) supports the idea of a dynamic difference between the two DNA structures that can influence cooperativity in binding. These results confirm the role of conformational entropy in positive cooperativity in some DNA interactions.     

Thermodynamic analysis of DNA binding by a Bacillus single stranded DNA binding protein

Background

Single-stranded DNA binding proteins (SSB) are essential for DNA replication, repair, and recombination in all organisms. SSB works in concert with a variety of DNA metabolizing enzymes such as DNA polymerase.

Results

We have cloned and purified SSB from Bacillus anthracis (SSB_BA). In the absence of DNA, at concentrations ??100 ??g/ml, SSB_BA did not form a stable tetramer and appeared to resemble bacteriophage T4 gene 32 protein. Fluorescence anisotropy studies demonstrated that SSB_BA bound ssDNA with high affinity comparable to other prokaryotic SSBs. Thermodynamic analysis indicated both hydrophobic and ionic contributions to ssDNA binding. FRET analysis of oligo(dT)₇₀ binding suggested that SSB_BA forms a tetrameric assembly upon ssDNA binding. This report provides evidence of a bacterial SSB that utilizes a novel mechanism for DNA binding through the formation of a transient tetrameric structure.

Conclusions

Unlike other prokaryotic SSB proteins, SSB_BA from Bacillus anthracis appeared to be monomeric at concentrations ??100 ??g/ml as determined by SE-HPLC. SSB_BA retained its ability to bind ssDNA with very high affinity, comparable to SSB proteins which are tetrameric. In the presence of a long ssDNA template, SSB_BA appears to form a transient tetrameric structure. Its unique structure appears to be due to the cumulative effect of multiple key amino acid changes in its sequence during evolution, leading to perturbation of stable dimer and tetramer formation. The structural features of SSB_BA could promote facile assembly and disassembly of the protein-DNA complex required in processes such as DNA replication.     

A noncanonical binding site of linezolid revealed via molecular dynamics simulations

Linezolid, an antibiotic of oxazolidinone family, is a translation inhibitor. The mechanism of its action that consists in preventing the binding of aminoacyl-tRNA to the A-site of the large subunit of a ribosome was embraced on the basis of the X-ray structural analysis of the linezolid complexes with vacant bacterial ribosomes. However, the known structures of the linezolid complexes with bacterial ribosomes poorly explain the linezolid selectivity in suppression of protein biosynthesis, depending on the amino acid sequence of the nascent peptide. In the present study the most probable structure of the linezolid complex with a E. coli ribosome in the A,A/P,P-state that is in line with the results of biochemical studies of linezolid action has been obtained by molecular dynamics simulation methods.     

Investigation of the copper binding site and the role of histidine as a ligand in riboflavin binding protein

Riboflavin Binding Protein (RBP) binds copper in a 1:1 molar ratio, forming a distinct well-ordered type II site. The nature of this site has been examined using X-ray absorption and pulsed electron paramagnetic resonance (EPR) spectroscopies, revealing a four coordinate oxygen/nitrogen rich environment. On the basis of analysis of the Cambridge Structural Database, the average protein bound copper-ligand bond length of 1.96 A, obtained by extended x-ray absorption fine structure (EXAFS), is consistent with four coordinate Cu(I) and Cu(II) models that utilize mixed oxygen and nitrogen ligand distributions. These data suggest a Cu-O 3N coordination state for copper bound to RBP. While pulsed EPR studies including hyperfine sublevel correlation spectroscopy and electron nuclear double resonance show clear spectroscopic evidence for a histidine bound to the copper, inclusion of a histidine in the EXAFS simulation did not lead to any significant improvement in the fit.     

Water PMF for predicting the properties of water molecules in protein binding site

Water is an important component in living systems and deserves better understanding in chemistry and biology. However, due to the difficulty of investigating the water functions in protein structures, it is usually ignored in computational modeling, especially in the field of computer‐aided drug design. Here, using the potential of mean forces (PMFs) approach, we constructed a water PMF (wPMF) based on 3946 non‐redundant high resolution crystal structures. The extracted wPMF potential was first used to investigate the structure pattern of water and analyze the residue hydrophilicity. Then, the relationship between wPMF score and the B factor value of crystal waters was studied. It was found that wPMF agrees well with some previously reported experimental observations. In addition, the wPMF score was also tested in parallel with 3D‐RISM to measure the ability of retrieving experimentally observed waters, and showed comparable performance but with much less computational cost. In the end, we proposed a grid‐based clustering scheme together with a distance weighted wPMF score to further extend wPMF to predict the potential hydration sites of protein structure. From the test, this approach can predict the hydration site at the accuracy about 80% when the calculated score lower than ?4.0. It also allows the assessment of whether or not a given water molecule should be targeted for displacement in ligand design. Overall, the wPMF presented here provides an optional solution to many water related computational modeling problems, some of which can be highly valuable as part of a rational drug design strategy. © 2012 Wiley Periodicals, Inc.     

A substituted triaza crown ether as a binding site in DNA conjugates

Synthesis of an asymmetrically substituted triaza crown ether, its incorporation into the 3'-end and 5'-end of ninemer oligonucleotides, and the influence of various alkanediamine ligands on duplex thermostabilities are reported.     

31P NMR investigation of backbone dynamics in DNA binding sites

The backbone conformation of DNA plays an important role in the indirect readout mechanisms for protein--DNA recognition events. Thus, investigating the backbone dynamics of each step in DNA binding sequences provides useful information necessary for the characterization of these interactions. Here, we use 31P dynamic NMR to characterize the backbone conformation and dynamics in the Dickerson dodecamer, a sequence containing the EcoRI binding site, and confirm solid-state 2H NMR results showing that the C3pG4 and C9pG10 steps experience unique dynamics and that these dynamics are quenched upon cytosine methylation. In addition, we show that cytosine methylation affects the conformation and dynamics of neighboring nucleotide steps, but this effect is localized to only near neighbors and base-pairing partners. Last, we have been able to characterize the percent BII in each backbone step and illustrate that the C3pG4 and C9pG10 favor the noncanonical BII conformation, even at low temperatures. Our results demonstrate that 31P dynamic NMR provides a robust and efficient method for characterizing the backbone dynamics in DNA. This allows simple, rapid determination of sequence-dependent dynamical information, providing a useful method for studying trends in protein-DNA recognition events.     

Backbone motions of free and pheromone-bound major urinary protein I studied by molecular dynamics simulation

Molecular motions of free and pheromone-bound mouse major urinary protein I, previously investigated by NMR relaxation, were simulated in 30 ns molecular dynamics (MD) runs. The backbone flexibility was described in terms of order parameters and correlation times, commonly used in the NMR relaxation analysis. Special attention was paid to the effect of conformational changes on the nanosecond time scale. Time-dependent order parameters were determined in order to separate motions occurring on different time scales. As an alternative approach, slow conformational changes were identified from the backbone torsion angle variances, and "conformationally filtered" order parameters were calculated for well-defined conformation states. A comparison of the data obtained for the free and pheromone-bound protein showed that some residues are more rigid in the bound form, but a larger portion of the protein becomes more flexible upon the pheromone binding. This finding is in general agreement with the NMR results. The higher flexibility observed on the fast (fs-ps) time scale was typically observed for the residues exhibiting higher conformational freedom on the ns time scale. An inspection of the hydrogen bond network provided a structural explanation for the flexibility differences between the free and pheromone-bound proteins in the simulations.     

Mapping the binding site topology of amyloid protein aggregates using multivalent ligands

A key process in the development of neurodegenerative diseases such as Alzheimer''s and Parkinson''s diseases is the aggregation of proteins to produce fibrillary aggregates with a cross β-sheet structure, amyloid. The development of reagents that can bind these aggregates with high affinity and selectivity has potential for early disease diagnosis. By linking two benzothiazole aniline (BTA) head groups with different length polyethylene glycol (PEG) spacers, fluorescent probes that bind amyloid fibrils with low nanomolar affinity have been obtained. Dissociation constants measured for interaction with Aβ, α-synuclein and tau fibrils show that the length of the linker determines binding affinity and selectivity. These compounds were successfully used to image α-synuclein aggregates in vitro and in the post-mortem brain tissue of patients with Parkinson''s disease. The results demonstrate that multivalent ligands offer a powerful approach to obtain high affinity, selective reagents to bind the fibrillary aggregates that form in neurodegenerative disease.

Multivalent ligands offer a powerful approach to obtain high affinity reagents to bind the aggregates that form in neurodegenerative disease. Selectivity for different proteins was achieved by using different linkers to connect the head groups.     

Fluorescence turn-on detection of a protein through the displaced single-stranded DNA binding protein binding to a molecular beacon

A new approach has been developed for the highly sensitive and selective sensing of a protein. Lysozyme binding to its aptamer prevents SSB protein binding, and the subsequent binding of the free SSB protein to a molecular beacon results in a turn-on fluorescence signal, which can be used for lysozyme quantification.     

Backbone and side-chain ordering in a small protein

We investigate the relation between backbone and side-chain ordering in a small protein. For this purpose, we have performed multicanonical simulations of the villin headpiece subdomain HP-36, an often used toy model in protein studies. Concepts of circular statistics are introduced to analyze side-chain fluctuations. In contrast to earlier studies on homopolypeptides [Wei et al., J. Phys. Chem. B 111, 4244 (2007)], we do not find collective effects leading to a separate transition. Rather, side-chain ordering is spread over a wide temperature range. Our results indicate a thermal hierarchy of ordering events, with side-chain ordering appearing at temperatures below the helix-coil transition but above the folding transition. We conjecture that this thermal hierarchy reflects an underlying temporal order, and that side-chain ordering facilitates the search for the correct backbone topology.     

Parametrization of calcium binding site in proteins and molecular dynamics simulation on phospholipase A2

The empirical energy parameters for a calcium ion and its ligands in proteins were determined within a pairwise additive framework. The interaction energies of Ca²⁺-water, Ca²⁺-peptide group and Ca²⁺-carboxyl group systems were calculated using the ab initio molecular orbital method with basis sets of double zeta quality including polarization or diffuse functions. The resulting potential energy surfaces served as references for the determination of the nonbonded parameters in the empirical energy function. The nonadditive corrections for the Ca²⁺-ligand pair potentials are incorporated implicitly in the nonbonded paremeters by treating three-body (1:2 complex) or seven-body (1:6 complex) systems in reference calculations. Ligand polarizations induced by Ca²⁺ are estimated from the partial atomic charges of two-body (1:1 complex) systems. The charge sets were determined by scaling so as to reproduce the reference potential energy surfaces. The newly determined parameter set was used in a stochastic boundary molecular dynamics simulation of phospholipase A₂. The solvated structure of the Ca²⁺-binding site obtained from an X-ray crystallographic study is well reproduced by the parameter set.     

Correlation of the side-chain hubs with the functional residues in DNA binding protein structures

The structure networks of DNA-binding proteins have been constructed and analyzed. The detailed analysis of the networks indicates a strong relation between the positions of the residues interacting with DNA and those that form extensive interactions within the protein structure (called hubs). This study shows that the functional residues in these proteins are held in place by efficient scaffolding of the structure using side-chain interactions, thus highlighting the role of these side-chain hubs with respect to the functional residues in the protein structure.     

Lipid binding to the carotenoid binding site in photosynthetic reaction centers

Lipid binding to the carotenoid binding site near the inactive bacteriochlorophyll monomer was probed in the reaction centers of carotenoid-less mutant, R-26 from Rhodobacter sphaeroides. Recently, a marked light-induced change of the local dielectric constant in the vicinity of the inactive bacteriochlorophyll monomer was reported in wild type that was attributed to structural changes that ultimately lengthened the lifetime of the charge-separated state by 3 orders of magnitude (Deshmukh, S. S.; Williams, J. C.; Allen, J. P.; Kalman, L. Biochemistry 2011, 50, 340). Here in the R-26 reaction centers, the combination of light-induced structural changes and lipid binding resulted in a 5 orders of magnitude increase in the lifetime of the charge-separated state involving the oxidized dimer and the reduced primary quinone in proteoliposomes. Only saturated phospholipids with fatty acid chains of 12 and 14 carbon atoms long were bound successfully at 8 °C by cooling the reaction center protein slowly from room temperature. In addition to reporting a dramatic increase of the lifetime of the charge-separated state at physiologically relevant temperatures, this study reveals a novel lipid binding site in photosynthetic reaction center. These results shed light on a new potential application of the reaction center in energy storage as a light-driven biocapacitor since the charges separated by ～30 ? in a low-dielectric medium can be prevented from recombination for hours.