NMR solution structure of an oxaliplatin 1,2-d(GG) intrastrand cross-link in a DNA dodecamer duplex

We have determined, at high resolution, the NMR solution structure of an oxaliplatin-GG DNA dodecamer in the AGGC sequence context by 2D NMR studies. Homonuclear assignment strategies resulted in unambiguous assignment of 203 out of 249 protons, which corresponds to assignment of approximately 81% of the protons. Assignments of H5' and H5" protons were tentative due to resonance overlap. The structure of the oxaliplatin duplex was calculated using the program CNS with a simulated annealing protocol. A total of 510 experimental restraints were employed in the structure calculation. Of 20 calculated structures, the 15 with the lowest energy were accepted as a family. The RMSD of the 15 lowest energy structures was 0.68 A, indicating good structural convergence. The theoretical NOESY spectrum obtained by back-calculation from the final average structure showed excellent agreement with the experimental data, indicating that the final structure was in good agreement with the experimental NMR data. Significant conformational differences were observed between the oxaliplatin-GG 12-mer DNA we studied and all previous solution structures of cisplatin-GG DNA duplexes. For example, the oxaliplatin-GG adduct shows much less distortion at the AG base-pair step than the cisplatin-GG adducts. In addition, the oxaliplatin-GG structure also has a narrow minor groove and an overall axis bend of about 31 degrees, both of which are very different from the recent NMR structures for the cisplatin-GG adducts. These structural differences may explain some of the biological differences between oxaliplatin- and cisplatin-GG adducts.     

Prediction of atomic structure from sequence for double helical DNA oligomers

DNA can adopt different conformations depending on the base sequence, solvent, electrolyte composition and concentration, pH, temperature, and interaction with proteins. Here we present a model for calculating the three-dimensional atomic structure of double-stranded DNA oligomers. A theoretical energy function is used for calculating the interactions within the base steps and an empirical backbone function is used to restrict the conformational space accessible to the bases and to account for the conformational coupling of neighboring steps in a sequence. Conformational searching on large structures or a large number of structures is possible, because each base step can be described by just two primary degrees of freedom (slide and shift). A genetic algorithm is used to search for low-energy structures in slide-shift space, and this allows very rapid optimization of DNA oligomers. The other base step parameters have been previously optimized for all possible slide-shift sequence combinations, and a heuristic algorithm is used to add the atomic details of the backbone conformation in the final step of the calculation. The structures obtained by this method are very similar to the corresponding X-ray crystal structures observed experimentally. The average RMSD is 2.24 Angstroms for a set of 20 oligomer structures. For 15 of these sequences, the X-ray crystal structure is the global energy minimum. The other 5 are bistable sequences that have B-form global energy minima but crystallize as A-DNA.     

Metropolis Monte Carlo calculations of DNA structure using internal coordinates and NMR distance restraints: An alternative method for generating a high-resolution solution structure

Summary A new method, a restrained Monte Carlo (rMC) calculation, is demonstrated for generating high-resolution structures of DNA oligonucleotides in solution from interproton distance restraints and bounds derived from complete relaxation matrix analysis of two-dimensional nuclear Overhauser effect (NOE) spectral peak intensities. As in the case of restrained molecular dynamics (rMD) refinement of structures, the experimental distance restraints and bounds are incorporated as a pseudo-energy term (or penalty function) into the mathematical expression for the molecular energy. However, the use of generalized helical parameters, rather than Cartesian coordinates, to define DNA conformation increases efficiency by decreasing by an order of magnitude the number of parameters needed to describe a conformation and by simplifying the potential energy profile. The Metropolis Monte Carlo method is employed to simulate an annealing process. The rMC method was applied to experimental 2D NOE data from the octamer duplex d(GTA-TAATG)·d(CATTATAC). Using starting structures from different locations in conformational space (e.g. A-DNA and B-DNA), the rMC calculations readily converged, with a root-mean-square deviation (RMSD) of <0.3 Å between structures generated using different protocols and starting structures. Theoretical 2D NOE peak intensities were calculated for the rMC-generated structures using the complete relaxation matrix program CORMA, enabling a comparison with experimental intensities via residual indices. Simulation of the vicinal proton coupling constants was carried out for the structures generated, enabling a comparison with the experimental deoxyribose ring coupling constants, which were not utilized in the structure determination in the case of the rMC simulations. Agreement with experimental 2D NOE and scalar coupling data was good in all cases. The rMC structures are quite similar to that refined by a traditional restrained MD approach (RMSD<0.5 Å) despite the different force fields used and despite the fact that MD refinement was conducted with additional restraints imposed on the endocyclic torsion angles of deoxyriboses. The computational time required for the rMC and rMD calculations is about the same. A comparison of structural parameters is made and some limitations of both methods are discussed with regard to the average nature of the experimental restraints used in the refinement.Abbreviations MC Monte Carlo - rMC restrained Monte Carlo - MD molecular dynamics - rMD restrained molecular dynamics - DG distance geometry - EM energy minimization - 2D NOE two-dimensional nuclear Overhauser effect - DQF-COSY double-quantum-filtered correlation spectroscopy - RMSD root-mean-square deviation To whom correspondence should be addressed.     

Spatial distribution of GC- and AT-rich DNA sequences within human chromosome territories

Previous topological analyses of DNA sequence organization in the interphase chromosome mainly focused on the spatial distribution of individual gene copies within chromosome territories. In order to achieve a more comprehensive view into the subchromosomal arrangement of DNA, we isolated the GC-richest/gene-richest fraction (H3 isochores) as well as AT-richest/gene-poorest fraction of human genomic DNA (L1+L2 isochores) and visualized the respective DNA within individual chromosome territories by means of dual-color FISH. Application of confocal laser scanning microscopy and dedicated 3D image analysis software, which differentiated territory subvolumes by peeling shells one voxel in width, revealed a significant difference in the intraterritorial distribution of these two DNA sequence classes. While the H3 isochores were found localized in all subvolumes of the territories at similar frequency, simultaneously detected L1+L2 isochores were observed more to the interior of the same chromosome territories. Thus the GC-rich sequences display a much higher variability in their intraterritorial localization than AT-rich DNA fragments.     

NMR structure of an intracellular loop peptide derived from prostaglandin EP3alpha receptor

We found that a peptide (EP3a: TIKALVSRCRAKAAV) corresponding to the N-terminal site of the intracellular third loop of human prostaglandin EP3α receptor could activate G protein α-subunit directly. The activity was almost same as Mastoparan-X, a G protein activating peptide from wasp venom. The three-dimensional molecular structure of the peptide in SDS-d₂₅ micelles was determined by 2D ¹H NMR spectroscopy. The structure of EP3a consists of a positive charge cluster on the C-terminal helical site. The cluster was also found in several corresponding receptor peptides. Therefore, the positive charge cluster on the helical structure might play a crucial role in activation of G protein.     

Solution structure of SWI1 AT-rich interaction domain from Saccharomyces cerevisiae and its nonspecific binding to DNA

SWI1 is a subunit of the SWI/SNF complex involved in chromatin remodeling. It contains an AT-rich interaction domain (ARID) which has the potential DNA binding activity. In this study, we determined the solution structure of the SWI1 ARID domain from Saccharomyces cerevisiae by nuclear magnetic resonance spectroscopy. Yeast SWI1 ARID domain is composed of seven alpha helices, six of which are conserved among the ARID family. In addition, the DNA-binding activity of the SWI1 ARID domain was confirmed by chemical shift perturbation assay. Similar to its human homolog, the yeast SWI1 ARID domain binds DNA nonspecifically.     

NMR solution structure of the theta subunit of DNA polymerase III from Escherichia coli

The catalytic core of Escherichia coli DNA polymerase III contains three tightly associated subunits (alpha, epsilon, and theta). The theta subunit is the smallest, but the least understood of the three. As a first step in a program aimed at understanding its function, the structure of the theta subunit has been determined by triple-resonance multidimensional NMR spectroscopy. Although only a small protein, theta was difficult to assign fully because approximately one-third of the protein is unstructured, and some sections of the remaining structured parts undergo intermediate intramolecular exchange. The secondary structure was deduced from the characteristic nuclear Overhauser effect patterns, the 3J(HN alpha) coupling constants and the consensus chemical shift index. The C-terminal third of the protein, which has many charged and hydrophilic amino acid residues, has no well-defined secondary structure and exists in a highly dynamic state. The N-terminal two-thirds has three helical segments (Gln10-Asp19, Glu38-Glu43, and His47-Glu54), one short extended segment (Pro34-Ala37), and a long loop (Ala20-Glu29), of which part may undergo intermediate conformational exchange. Solution of the three-dimensional structure by NMR techniques revealed that the helices fold in such a way that the surface of theta is bipolar, with one face of the protein containing most of the acidic residues and the other face containing most of the long chain basic residues. Preliminary chemical shift mapping experiments with a domain of the epsilon subunit have identified a loop region (Ala20-Glu29) in theta as the site of association with epsilon.     

NMR structure of duplex DNA containing the α‐OH‐PdG·dA base pair: A mutagenic intermediate of acrolein

Acrolein, a cell metabolic product and main component of cigarette smoke, reacts with DNA generating α‐OH‐PdG lesions, which have the ability to pair with dATP during replication thereby causing G to T transversions. We describe the solution structure of an 11‐mer DNA duplex containing the mutagenic α‐OH‐PdG·dA base pair intermediate, as determined by solution nuclear magnetic resonance (NMR) spectroscopy and retrained molecular dynamics (MD) simulations. The NMR data support a mostly regular right‐handed helix that is only perturbed at its center by the presence of the lesion. Undamaged residues of the duplex are in anti orientation, forming standard Watson‐Crick base pairs alignments. Duplication of proton signals at and near the damaged base pair reveals the presence of two enantiomeric duplexes, thus establishing the exocyclic nature of the lesion. The α‐OH‐PdG adduct assumes a syn conformation pairing to its partner dA base that is protonated at pH 6.6. The three‐dimensional structure obtained by restrained molecular dynamics simulations show hydrogen bond interactions that stabilize α‐OH‐PdG in a syn conformation and across the lesion containing base pair. We discuss the implications of the structures for the mutagenic bypass of acrolein lesions. © 2010 Wiley Periodicals, Inc. Biopolymers 93: 391–401, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

HMG-D complexed to a bulge DNA: an NMR model

An NMR model is presented for the structure of HMG-D, one of the DROSOPHILA: counterparts of mammalian HMG1/2 proteins, bound to a particular distorted DNA structure, a dA(2) DNA bulge. The complex is in fast to intermediate exchange on the NMR chemical shift time scale and suffers substantial linebroadening for the majority of interfacial resonances. This essentially precludes determination of a high-resolution structure for the interface based on NMR data alone. However, by introducing a small number of additional constraints based on chemical shift and linewidth footprinting combined with analogies to known structures, an ensemble of model structures was generated using a computational strategy equivalent to that for a conventional NMR structure determination. We find that the base pair adjacent to the dA(2) bulge is not formed and that the protein recognizes this feature in forming the complex; intermolecular NOE enhancements are observed from the sidechain of Thr 33 to all four nucleotides of the DNA sequence step adjacent to the bulge. Our results form the first experimental demonstration that when binding to deformed DNA, non-sequence-specific HMG proteins recognize the junction between duplex and nonduplex DNA. Similarities and differences of the present structural model relative to other HMG-DNA complex structures are discussed.     

DNA sequence organization and transcription of the chicken genome

A new approach has been used to examine DNA sequence organization in the chicken genome. The interspersion pattern was determined by studying the fraction of labelled DNA fragments of different lengths that hybridized to an excess of short chicken repeated DNA sequences. The results indicate that chicken DNA has a pattern of sequence organization quite different than the standard ‘Xenopus’ or ‘Drosophila’ patterns. Two classes of unique sequences are found. One, 34% of the genome, consists of unique sequences approx. 4 kb long interspersed with repeated sequences. The second, non-interspersed fraction, 38% of the genome, consists of unique sequences found in long tracts, a minimum of approx. 22 kb in length. In an attempt to determine whether a relationship exists between DNA sequence organization and the distribution of structural genes we have isolated chicken DNA sequences belonging to different interspersion classes and tested each for the presence of structural genes by hybridization to excess poly(A)⁺ mRNA. Sequences complementary to poly(A)⁺ mRNA can be found with approximately the same frequency in both the non-interspersed fraction of the genome and a repeat-contiguous fraction enriched for interspersed sequences.     

DNA sequence,structure, and phylogenetic relationship of the small subunit rRNA coding region of mitochondrial DNA fromPodospora anserina

Summary DNA sequence analysis and the localization of the 5 and 3 termini by S1 mapping have shown that the mitochondrial (mt) small subunit rRNA coding region fromPodospora anserina is 1980 bp in length. The analogous coding region for mt rRNA is 1962 bp in maize, 1686 bp inSaccharomyces cerevisiae, and 956 bp in mammals, whereas its counterpart inEscherichia coli is 1542 bp. TheP. anserina mt 16S-like rRNA is 400 bases longer than that fromE. coli, but can be folded into a similar secondary structure. The additional bases appear to be clustered at specific locations, including extensions at the 5 and 3 termini. Comparison with secondary structure diagrams of 16S-like RNAs from several organisms allowed us to specify highly conserved and variable regions of this gene. Phylogenetic tree construction indicated that this gene is grouped with other mitochondrial genes, but most closely, as expected, with the fungal mitochondrial genes.     

Insertion sequence IS5 contains a sharply curved DNA structure at its terminus

Summary It was demonstrated that insertion sequence IS5 contains a sequence-directed bent (sharply curved) DNA structure at its terminus, close to one of its 16 bp terminal repeats. The minimal number of copies of IS5 related sequences and the locations of the latter on the Escherichia coli K12 W3110 chromosome were determined. Evidence is presented of the occurrence of IS5 mediated translocation and duplication of a large DNA segment on the E. coli chromosome.     

Distinctive features in the structure and dynamics of the DNA repeat sequence GGCGGG

G-rich DNA has been known to form a variety of folded and multistranded structures, with even single base modifications causing important structural changes. But, very little is known about the dynamic characteristics of the structures, which may play crucial roles in facilitating the structural transitions. In this background, we report here NMR investigations on the structure and dynamics of a DNA repeat sequence GGCGGG in aqueous solution containing Na+ ions at neutral pH. The chosen sequence d-TGGCGGGT forms a parallel quadruplex with a C-tetrad in the middle, formed by symmetrical pairing of four Cs in a plane via NH2-O2 H-bonds. 13C relaxation measurements at natural abundance for C' sugar carbons provided valuable insight into the sequence specific dynamism of G and C-tetrads in the quadruplex. The C4 tetrad seems to introduce high conformational dynamism at milli- to micro-second time scale in the quadruplex. Concomitantly, there is a decrease in the pico-second time scale dynamics. Interestingly, these effects are seen more prominently at the G-tetrads on the 3' end of C-tetrad than on its 5' end. These observations would have important implications for the roles the tetrads may play in many biological functions.     

Molecular Dynamics Simulation in Solvent of the Estrogen Receptor Protein DNA Binding Domain in Complex with a Non-Consensus Estrogen Response Element DNA Sequence

Abstract

We investigated protein/DNA interactions, using molecular dynamics simulations computed between a 10 Angstom water layer model of the estrogen receptor (ER) protein DNA binding domain (DBD) amino acids and DNA of a non-consensus estrogen response element (ERE) consisting of 29 nucleotide base pairs. This ERE nucleotide sequence occurs naturally upstream of the Xenopus laevis Vitelligenin AI gene. The ER DBD is encoded by three exons. Namely, exons 2 and 3 which encode the two zinc binding motifs and a sequence of exon 4 which encodes a predicted alpha helix. We generated a computer model of the ER DBD using atomic coordinates derived from the average of 30 nuclear magnetic resonance (NMR) spectroscopy coordinate sets. Amino acids on the carboxyl end of the ER DBD were disordered in both X-ray crystallography and NMR determinations and no coordinates were reported. This disordered region includes 10 amino acids of a predicted alpha helix encoded in exon 4 at the exon 3/4 splice junction. These amino acids are known to be important in DNA binding and are also believed to function as a nuclear translocation signal sequence for the ER protein. We generated a computer model of the predicted alpha helix consisting of the 10 amino acids encoded in exon 4 and attached this helix to the carboxyl end of the ER DBD at the exon 3/4 splice junction site. We docked the ER DBD model within the DNA major groove halfsites of the 29 base pair non-consensus ERE and flanking nucleotides. We constructed a solvated model with the ER DBD/ERE complex surrounded by a ten Angstrom water layer and conducted molecular dynamics simulations. Hydrogen bonding interactions were monitored. In addition, van der Waals and electrostatic interaction energies were calculated. Amino acids of the ER DBD DNA recognition helix formed both direct and water mediated hydrogen bonds at cognate codon-anticodon nucleotide base and backbone sites within the ERE DNA right major groove halfsite. Amino acids of the ER DBD exon 4 encoded predicted alpha helix formed direct and water mediated H-bonds with base and backbone sites of their cognate codon-anticodon nucleotides within the minor grooves flanking the ERE DNA major groove halfsites. These interactions together induced bending of the DNA into the protein.     

NMR structure of a 4 x 4 nucleotide RNA internal loop from an R2 retrotransposon: identification of a three purine-purine sheared pair motif and comparison to MC-SYM predictions

The NMR solution structure is reported of a duplex, 5'GUGAAGCCCGU/3'UCACAGGAGGC, containing a 4 × 4 nucleotide internal loop from an R2 retrotransposon RNA. The loop contains three sheared purine-purine pairs and reveals a structural element found in other RNAs, which we refer to as the 3RRs motif. Optical melting measurements of the thermodynamics of the duplex indicate that the internal loop is 1.6 kcal/mol more stable at 37°C than predicted. The results identify the 3RRs motif as a common structural element that can facilitate prediction of 3D structure. Known examples include internal loops having the pairings: 5'GAA/3'AGG, 5'GAG/3'AGG, 5'GAA/3'AAG, and 5'AAG/3'AGG. The structural information is compared with predictions made with the MC-Sym program.     

NMR solution structure of the acylphosphatase from Escherichia coli

The solution structure of Escherichia coli acylphosphatase (E. coli AcP), a small enzyme catalyzing the hydrolysis of acylphosphates, was determined by (1)H and (15)N NMR and restrained modelling calculation. In analogy with the other members of AcP family, E. coli AcP shows an alpha/beta sandwich domain composed of four antiparallel and one parallel beta-strand, assembled in a five-stranded beta-sheet facing two antiparallel alpha-helices. The pairwise RMSD values calculated for the backbone atoms of E. coli and Sulfolobus solfataricus AcP, Bovine common type AcP and Horse muscle AcP are 2.18, 5.31 and 5.12 A, respectively. No significant differences are present in the active site region and the catalytic residue side chains are consistently positioned in the structures.