cgDNA: a software package for the prediction of sequence-dependent coarse-grain free energies of B-form DNA

cgDNA is a package for the prediction of sequence-dependent configuration-space free energies for B-form DNA at the coarse-grain level of rigid bases. For a fragment of any given length and sequence, cgDNA calculates the configuration of the associated free energy minimizer, i.e. the relative positions and orientations of each base, along with a stiffness matrix, which together govern differences in free energies. The model predicts non-local (i.e. beyond base-pair step) sequence dependence of the free energy minimizer. Configurations can be input or output in either the Curves+ definition of the usual helical DNA structural variables, or as a PDB file of coordinates of base atoms. We illustrate the cgDNA package by comparing predictions of free energy minimizers from (a) the cgDNA model, (b) time-averaged atomistic molecular dynamics (or MD) simulations, and (c) NMR or X-ray experimental observation, for (i) the Dickerson–Drew dodecamer and (ii) three oligomers containing A-tracts. The cgDNA predictions are rather close to those of the MD simulations, but many orders of magnitude faster to compute. Both the cgDNA and MD predictions are in reasonable agreement with the available experimental data. Our conclusion is that cgDNA can serve as a highly efficient tool for studying structural variations in B-form DNA over a wide range of sequences.     

The analysis of stress-induced duplex destabilization in long genomic DNA sequences.

We present a method for calculating predicted locations and extents of stress-induced DNA duplex destabilization (SIDD) as functions of base sequence and stress level in long DNA molecules. The base pair denaturation energies are assigned individually, so the influences of near neighbors, methylated bases, adducts, or lesions can be included. Sample calculations indicate that copolymeric energetics give results that are close to those derived when full near-neighbor energetics are used; small but potentially informative differences occur only in the calculated SIDD properties of moderately destabilized regions. The method presented here for analyzing long sequences calculates the destabilization properties within windows of fixed length N, with successive windows displaced by an offset distance d(o). The final values of the relevant destabilization parameters for each base pair are calculated as weighted averages of the values computed for each window in which that base pair appears. This approach implicitly assumes that the strength of the direct coupling between remote base pairs that is induced by the imposed stress attenuates with their separation distance. This strategy enables calculations of the destabilization properties of DNA sequences of any length, up to and including complete chromosomes. We illustrate its utility by calculating the destabilization properties of the entire E. coli genomic DNA sequence. A preliminary analysis of the results shows that promoters are associated with SIDD regions in a highly statistically significant manner, suggesting that SIDD attributes may prove useful in the computational prediction of promoter locations in prokaryotes.     

Thermodynamic parameters for loop formation in RNA and DNA hairpin tetraloops.

We determined the melting temperatures (Tm) and thermodynamic parameters of 15 RNA and 19 DNA hairpins at 1 M NaCl, 0.01 M sodium phosphate, 0.1 mM EDTA, at pH 7. All these hairpins have loops of four bases, the most common loop size in 16S and 23S ribosomal RNAs. The RNA hairpins varied in loop sequence, loop-closing base pair (A.U, C.G, or G.C), base sequence of the stem, and stem size (four or five base pairs). The DNA hairpins varied in loop sequence, loop-closing base pair (C.G, or G.C), and base sequence of the four base-pair stem. Thermodynamic properties of a hairpin may be represented by nearest-neighbor interactions of the stem plus contributions from the loop. Thus, we obtained thermodynamic parameters for the formation of RNA and DNA tetraloops. For the tetraloops we studied, a free energy of loop formation (at 37 degrees C) of about +3 kcal/mol is most common for either RNA or DNA. There are extra stable loops with delta G degrees 37 near +1 kcal/mol, but the sequences are not necessarily the same for RNA and DNA. The closing base pair is also important; changing from C.G to G.C lowered the stability of several tetraloops in both RNA and DNA. These values will be useful in predicting RNA and DNA secondary structures.     

Nonplanar DNA Base Pairs

Abstract

Three-dimensional structures of a representative set of more than 30 hydrogen-bonded nucleic acids base pairs have been studied by reliable ab initio quantum mechanical methods. We show that many hydrogen-bonded nucleic acid base pairs are intrinsically nonplanar, mainly due to the partial sp³ hybridization of nitrogen atoms of their amino groups and secondary electrostatic interactions. This finding extends the variability of intermolecular interactions of DNA bases in that i) flexibility of the base pairs is larger than has been assumed before, and ii) attractive proton-proton acceptor interactions oriented out of the base pair plane are allowed. For example, all four G…A mismatch base pairs are propeller twisted, and the energy preferences for the nonplanar structures range from less than 0.1 kcal/mol to 1.8 kcal/mol. We predict that nonplanarity of the amino group of guanine in the G(anti)…A(anti) pair of the ApG step of the d(CCAAGATTGG)₂ crystal structure is an important stabilizing factor that improves the energy of this structure by almost 3 kcal/mol. Currently used empirical potentials are not accurate enough to properly cover the interactions associated with amino-group and base-pair nonplanarity.     

SSCP primer design based on single-strand DNA structure predicted by a DNA folding program.

To predict alterations in single-strand DNA mobility in non-denaturing electrophoretic gels, Zuker's RNA folding program was modified. Energy files utilized by the LRNA RNA folding algorithm were modified to emulate folding of single-strand DNA. Energy files were modified to disallow G-T base pairing. Stacking energies were corrected for DNA thermodynamics. Constraints on loop nucleotide sequences were removed. The LRNA RNA folding algorithm using the DNA fold energy files was applied to predict folding of PCR generated single-strand DNA molecules from polymorphic human ALDH2 and TPH alleles. The DNA-Fold version 1.0 program was used to design primers to create and abolish SSCP mobility shifts. Primers were made that add a 5' tag sequence or alter complementarity to an internal sequence. Differences in DNA secondary structure were assessed by SSCP analysis and compared to single-strand DNA secondary structure predictions. Results demonstrate that alterations in single-strand DNA conformation may be predicted using DNA-Fold 1.0.     

Theoretical studies on interaction of anticancer drugs (dacarbazine,procarbazine and triethylenemelamine) with normal (AT and GC) and mismatch (GG,CC, AA and TT) base pairs

This study is an attempt to gain a better understanding of the physicochemical interaction between novel anticancer drugs and DNA bases. We have employed quantum chemical tools to explore the interaction of a few anticancer drugs [namely procarbazine (PR), dacarbazine (DC) and triethylenemelamine (TR)] with isolated normal (GC and AT) and mismatch (AA, CC, GG and TT) base pairs. The molecular geometries, electronic structural stability, vibrational energies, chemical reactivity and other electronic properties were studied using MP2/6-311+G^**, B3LYP/6-311+G^** and M05-2X/6-311+G^** methods. The optimised geometries of the usual and mismatch base pairs are almost planar whereas the geometries of drug-interacting complexes deviate from planarity. The presence of steric hindrance and π-bond overlaps between C–C bonds in the complexes has distorted the planarity of the four- and five-member rings in the base pairs. Among the three drugs chosen, DC and PR bond well with normal and mismatch base pairs with large interaction energy. The electron density (ED) difference maps of the most stable GG–DC, GG–PR and GG–TR drug-interacting complexes show the information about sharing of ED and gain or loss of ED within the interacting molecules. The stabilisation energy of the charge transfer interaction between the relevant donor–acceptor orbital of GG–DC and GC–DC complexes has been found to be around 16 kcal/mol and GG–PR and GC–PR complexes has been found to be around 12 kcal/mol. But, for the GG–TR and GC–TR complexes, the stabilisation energy is found to be less than 6 kcal/mol. Moreover, the topological analysis of hydrogen bond network of DC and PR drug-interacting complexes have high electron and Laplacian density with structural stability at the bond critical points (BCPs), while compared TR drug-interacting complexes by atoms in molecules and natural bond orbital analysis. Finally, we may conclude that the drugs DC and PR are highly efficient drugs to target normal and mismatch base pair for control and inhibition of DNA replication.     

DNA analysis servers: plot.it,bend.it,model.it and IS

The WWW servers at http://www.icgeb.trieste.it/dna/ are dedicated to the analysis of user-submitted DNA sequences; plot.it creates parametric plots of 45 physicochemical, as well as statistical, parameters; bend.it calculates DNA curvature according to various methods. Both programs provide 1D as well as 2D plots that allow localisation of peculiar segments within the query. The server model.it creates 3D models of canonical or bent DNA starting from sequence data and presents the results in the form of a standard PDB file, directly viewable on the user's PC using any molecule manipulation program. The recently established introns server allows statistical evaluation of introns in various taxonomic groups and the comparison of taxonomic groups in terms of length, base composition, intron type etc. The options include the analysis of splice sites and a probability test for exon-shuffling.     

Selective Preference of Parallel DNA Triplexes Is Due to the Disruption of Hoogsteen Hydrogen Bonds Caused by the Severe Nonisostericity between the G*GC and T*AT Triplets

Implications of DNA, RNA and RNA.DNA hybrid triplexes in diverse biological functions, diseases and therapeutic applications call for a thorough understanding of their structure-function relationships. Despite exhaustive studies mechanistic rationale for the discriminatory preference of parallel DNA triplexes with G_*GC & T_*AT triplets still remains elusive. Here, we show that the highest nonisostericity between the G_*GC & T_*AT triplets imposes extensive stereochemical rearrangements contributing to context dependent triplex destabilisation through selective disruption of Hoogsteen scheme of hydrogen bonds. MD simulations of nineteen DNA triplexes with an assortment of sequence milieu reveal for the first time fresh insights into the nature and extent of destabilization from a single (non-overlapping), double (overlapping) and multiple pairs of nonisosteric base triplets (NIBTs). It is found that a solitary pair of NIBTs, feasible either at a G_*GC/T_*AT or T_*AT/G_*GC triplex junction, does not impinge significantly on triplex stability. But two overlapping pairs of NIBTs resulting from either a T_*AT or a G_*GC interruption disrupt Hoogsteen pair to a noncanonical mismatch destabilizing the triplex by ~10 to 14 kcal/mol, implying that their frequent incidence in multiples, especially, in short sequences could even hinder triplex formation. The results provide (i) an unambiguous and generalised mechanistic rationale for the discriminatory trait of parallel triplexes, including those studied experimentally (ii) clarity for the prevalence of antiparallel triplexes and (iii) comprehensive perspectives on the sequence dependent influence of nonisosteric base triplets useful in the rational design of TFO’s against potential triplex target sites.     

VADAR: a web server for quantitative evaluation of protein structure quality

VADAR (Volume Area Dihedral Angle Reporter) is a comprehensive web server for quantitative protein structure evaluation. It accepts Protein Data Bank (PDB) formatted files or PDB accession numbers as input and calculates, identifies, graphs, reports and/or evaluates a large number (>30) of key structural parameters both for individual residues and for the entire protein. These include excluded volume, accessible surface area, backbone and side chain dihedral angles, secondary structure, hydrogen bonding partners, hydrogen bond energies, steric quality, solvation free energy as well as local and overall fold quality. These derived parameters can be used to rapidly identify both general and residue-specific problems within newly determined protein structures. The VADAR web server is freely accessible at http://redpoll.pharmacy.ualberta.ca/vadar.     

Conformations of DNA Duplexes Containing 8-Oxoguanine

Abstract

As a step towards elucidating the mechanisms of mutagenesis induced by irradiation and oxidation, we study the incorporation of 8-oxoguanine (OG) into duplex DNA. Molecular modelling is used to reveal changes in DNA conformational parameters due to mispairs within the sequences d(A₅XA₅) · d(T₅YT₅) and d(G₅XG₅) · d(C₅YC₅) where one of the bases of the bases of the central X:Y pair is OG and the other A T, G or C. The G:C to OG:C replacements in DNA duplexes produce only minor conformational changes, similar to normal base sequence effects. The calculations suggest that both OG(syn): G and OG(syn): A mispairs can also be introduced without drastic distortion of sugar-phosphate backbone. The distortions produced by OG-containing mispairs are also found to be sequence dependent Overall these calculations suggest that the G→OG conversion could be an important factor in the irradiative or oxidative damage of DNA.     

Conformational changes of the phenyl and naphthyl isocyanate-DNA adducts during DNA replication and by minor groove binding molecules

DNA lesions produced by aromatic isocyanates have an extra bulky group on the nucleotide bases, with the capability of forming stacking interaction within a DNA helix. In this work, we investigated the conformation of the 2′-deoxyadenosine and 2′-deoxycytidine derivatives tethering a phenyl or naphthyl group, introduced in a DNA duplex. The chemical modification experiments using KMnO₄ and 1-cyclohexyl-3 -(2-morpholinoethyl) carbodiimide metho-p-toluenesulfonate have shown that the 2′-deoxycytidine lesions form the base pair with guanine while the 2′-deoxyadenosine lesions have less ability of forming the base pair with thymine in solution. Nevertheless, the kinetic analysis shows that these DNA lesions are compatible with DNA ligase and DNA polymerase reactions, as much as natural DNA bases. We suggest that the adduct lesions have a capability of adopting dual conformations, depending on the difference in their interaction energies between stacking of the attached aromatic group and base pairing through hydrogen bonds. It is also presented that the attached aromatic groups change their orientation by interacting with the minor groove binding netropsin, distamycin and synthetic polyamide. The nucleotide derivatives would be useful for enhancing the phenotypic diversity of DNA molecules and for exploring new non-natural nucleotides.     

An unnatural base pair system for efficient PCR amplification and functionalization of DNA molecules

Toward the expansion of the genetic alphabet, we present an unnatural base pair system for efficient PCR amplification, enabling the site-specific incorporation of extra functional components into DNA. This system can be applied to conventional PCR protocols employing DNA templates containing unnatural bases, natural and unnatural base triphosphates, and a 3′→5′ exonuclease-proficient DNA polymerase. For highly faithful and efficient PCR amplification involving the unnatural base pairing, we identified the natural-base sequences surrounding the unnatural bases in DNA templates by an in vitro selection technique, using a DNA library containing the unnatural base. The system facilitates the site-specific incorporation of a variety of modified unnatural bases, linked with functional groups of interest, into amplified DNA. DNA fragments (0.15 amol) containing the unnatural base pair can be amplified 10⁷-fold by 30 cycles of PCR, with <1% total mutation rate of the unnatural base pair site. Using the system, we demonstrated efficient PCR amplification and functionalization of DNA fragments for the extremely sensitive detection of zeptomol-scale target DNA molecules from mixtures with excess amounts (pmol scale) of foreign DNA species. This unnatural base pair system will be applicable to a wide range of DNA/RNA-based technologies.     

MutL: conducting the cell's response to mismatched and misaligned DNA

Base pair mismatches in DNA arise from errors in DNA replication, recombination, and biochemical modification of bases. Mismatches are inherently transient. They are resolved passively by DNA replication, or actively by enzymatic removal and resynthesis of one of the bases. The first step in removal is recognition of strand discontinuity by one of the MutS proteins. Mismatches arising from errors in DNA replication are repaired in favor of the base on the template strand, but other mismatches trigger base excision or nucleotide excision repair (NER), or non‐repair pathways such as hypermutation, cell cycle arrest, or apoptosis. We argue that MutL homologues play a key role in determining biologic outcome by recruiting and/or activating effector proteins in response to lesion recognition by MutS. We suggest that the process is regulated by conformational changes in MutL caused by cycles of ATP binding and hydrolysis, and by physiologic changes which influence effector availability.     

The P1 phage replication protein RepA contacts an otherwise inaccessible thymine N3 proton by DNA distortion or base flipping

The RepA protein from bacteriophage P1 binds DNA to initiate replication. RepA covers one face of the DNA and the binding site has a completely conserved T that directly faces RepA from the minor groove at position +7. Although all four bases can be distinguished through contacts in the major groove of B-form DNA, contacts in the minor groove cannot easily distinguish between A and T bases. Therefore the 100% conservation at this position cannot be accounted for by direct contacts approaching into the minor groove of B-form DNA. RepA binding sites with modified base pairs at position +7 were used to investigate contacts with RepA. The data show that RepA contacts the N3 proton of T at position +7 and that the T=A hydrogen bonds are already broken in the DNA before RepA binds. To accommodate the N3 proton contact the T₊₇ /A_+7^′ base pair must be distorted. One possibility is that T₊₇ is flipped out of the helix. The energetics of the contact allows RepA to distinguish between all four bases, accounting for the observed high sequence conservation. After protein binding, base pair distortion or base flipping could initiate DNA melting as the second step in DNA replication.     

Probing the role of hydrogen bonds in the stability of base pairs in double-helical DNA

Aromatic stacking and hydrogen bonding between nucleobases are two of the key interactions responsible for stabilization of DNA double-helical structures. The present work aims at defining the specific contributions of these interactions to the stability of individual base pairs in DNA. The two DNA double helices investigated are formed, respectively, by the palindromic base sequences 5'-dCCAACGTTGG-3' and 5'-dCGCAGATCTGCG-3'. The strength of the N==H...N inter-base hydrogen bond in each base pair is characterized from the measurement of the protium-deuterium fractionation factor of the corresponding imino proton using NMR spectroscopy. The structural stability of each base pair is evaluated from the exchange rate of the imino proton, measured by NMR. The results reveal that the fractionation factors of the imino protons in the two DNA double helices investigated fall within a narrow range of values, between 0.92 and 1.0. In contrast, the free energies of structural stabilization for individual base pairs span 3.5 kcal/mol, from 5.2 to 8.7 kcal/mol (at 15 degrees C). These findings indicate that, in the two DNA double helices investigated, the strength of N==H...N inter-base hydrogen bonds does not change significantly depending on the nature or the sequence context of the base pair. Hence, the variations in structural stability detected by proton exchange do not involve changes in the strength of inter-base hydrogen bonds. Instead, the results suggest that the energetic identity of a base pair is determined by the number of inter-base hydrogen bonds, and by the stacking interactions with neighboring base pairs.     

Thermodynamic properties of the specific binding between Ag+ ions and C:C mismatched base pairs in duplex DNA

Metal-mediated base pairs formed by the interaction between metal ions and artificial bases in oligonucleotides have been developed for potential applications in nanotechnology. We recently found that a natural C:C mismatched base pair bound to an Ag(+) ion to generate a novel metal-mediated base pair in duplex DNA. Preparation of the novel C-Ag-C base pair involving natural bases is more convenient than that of metal-mediated base pairs involving artificial bases because time-consuming base synthesis is not required. Here, we examined the thermodynamic properties of the binding between the Ag(+) ion and each of single and double C:C mismatched base pair in duplex DNA by isothermal titration calorimetry. The Ag(+) ion specifically bound to the C:C mismatched base pair at a 1:1 molar ratio with 10(6) M(-1) binding constant, which was significantly larger than those for nonspecific metal ion-DNA interactions. The specific binding between the Ag(+) ion and the single C:C mismatched base pair was mainly driven by the positive dehydration entropy change and the negative binding enthalpy change. In the interaction between the Ag(+) ion and each of the consecutive and interrupted double C:C mismatched base pairs, stoichiometric binding at a 1:1 molar ratio was achieved in each step of the first and second Ag(+) binding. The binding affinity for the second Ag(+) binding was similar to that for the first Ag(+) binding. Stoichiometric binding without interference and negative cooperativity may be favorable for aligning multiple Ag(+) ions in duplex DNA for applications of the metal-mediated base pairs in nanotechnology.     

Cognate base‐pair selectivity of hydrophobic unnatural bases in DNA ligation by T4 DNA ligase

We present cognate base pair selectivity in template‐dependent ligation by T4 DNA ligase using a hydrophobic unnatural base pair (UBP), Ds‐Pa. T4 DNA ligase efficiently recognizes the Ds‐Pa pairing at the conjugation position, and Ds excludes the noncognate pairings with the natural bases. Our results indicate that the hydrophobic base pairing is allowed in enzymatic ligation with higher cognate base‐pair selectivity, relative to the hydrogen‐bond interactions between pairing bases. The efficient ligation using Ds‐Pa can be employed in recombinant DNA technology using genetic alphabet expansion, toward the creation of semi‐synthetic organisms containing UBPs.     

Thermodynamic, spectroscopic, and equilibrium binding studies of DNA sequence context effects in six 22-base pair deoxyoligonucleotides.

Effects of different end sequences on stability, circular dichroism spectra (CD), and enzyme binding properties were investigated for six 22-base pair, non-self-complementary duplex DNA oligomers. The center sequences of these deoxyoligonucleotides have 8-14 base pairs in common and are flanked on both sides by sequences differing in context and A-T content. Temperature-induced melting transitions monitored by differential scanning calorimetry (DSC) and ultraviolet absorbance were measured for the six duplexes in buffered 115 mM Na(+) solutions. Values of the melting transition enthalpy, DeltaH(cal), and entropy, DeltaS(cal), were obtained directly from DSC experiments. Melting transition parameters, DeltaH(vH) and DeltaS(vH), were also estimated from van't Hoff analysis of optical melting curves collected as a function of DNA concentration, assuming a two-state melting transition. Melting free energies (20 degrees C) of the six DNAs evaluated from DSC experiments ranged from -18.7 to -32.7 kcal/mol. van't Hoff estimates of the free energies ranged from -18.5 to -48.0 kcal/mol. With either method, the trends in free energy as a function of sequence were identical. Equilibrium binding by BamHI restriction endonuclease to the 22-base pair DNAs was also investigated. The central eight base pairs of all six molecules, 5'-A-GGATCC-A-3', contained a BamHI recognition sequence bounded by A-T base pairs. Magnesium free binding assays were performed by titering BamHI against a constant concentration of each of the deoxyoligonucleotide substrates and analyzing reaction products by gel retardation. Binding isotherms of the total amount of bound DNA versus protein concentration were constructed which provided semiquantitative estimates of the equilibrium dissociation constants for dissociation of BamHI from the six DNA oligomers. Dissociation constants ranged from 0.5 x 10(-)(9) to 12.0 x 10(-)(9) M with corresponding binding free energies of -12.5 to -10.6 (+/-0. 1) kcal/mol. An inverse relationship is found when binding and stability are compared.