Nature of the stacking of nucleic acid bases in water: a Monte Carlo simulation

The results of a Monte Carlo simulation of the hydration of uracil and thymine molecules, their stacked dimers and hydrogen-bonded base pairs are presented. Simulations have been performed in a cluster approximation. The semiempirical atom-atom potential functions have been used (cluster consisting of 200 water molecules). It has been shown that the stacking interactions of uracil and thymine molecules in water arise mainly due to the increase in the water-water interaction during the transition from monomers to dimer. It has been found out that stacked base associates are more preferable than base pairs in water. This preference is mainly due to the energetically more favourable structure of water around the stack.     

Crystal Structure of A Z-DNA Fragment Containing Thymine/2-Aminoadenine Base Pairs

Abstract

The Z-DNA structure has been shown to form in two crystals made from self-complementary DNA hexamers d(CGTDCG) and d(CDCGTG) which contain thymine/2-ammoadenine (TD) base pairs. The latter structure has been solved and refined to 1.3 Å resolution and it shows only small conformational changes due to the introduction of the TD base pairs in comparison with the structure of d(CG)3. Spectroscopic studies with these compounds demonstrate that DNA molecules containing 2-aminoadenine residues form Z-DNA slightly more easily than do those containing adenine nucleotides, but not as readily as the parent sequence containing only guanine-cytosine base pairs.     

Hydration of Uracil and Thymine Methylderivatives: a Monte Carlo Simulation

Abstract

The simulation performed shows that under methylation of uracil and thymine NH-groups the interaction energy between a base and water (U_wb) is increased. It is also detected that the increase in this energy was observed in the 1st and the 3rd sectors. These conclusions do not confirm the assumption made in the literature on the character of an interaction between methylated bases and water. According to this assumption, when the NH-groups are methylated, the energy of U_wb in these sectors decreases as a result of the van der Waals interactions between a methyl group and water, whose energy compensates the increase in the U_wb energy due to the breaking of an H-bond. Regularity of water molecules near a hydrophobic group under the hydration of polar molecules is detected for the first time.     

96 Study of the stacking of DNA homoassociates by the extended cluster approach

Metropolis Monte Carlo method based on the extended cluster approach (Danilov, Dailidonis, van Mourik, & Fruchtl, 2011a, 2011b; Dailidonis, Danilov, Früchtl, & van Mourik, 2011) is used to investigate adenine–adenine (AA), guanine–guanine (GG), thymine–thymine (TT), and cytosine–cytosine (CC) homoassociates in a cluster consisting of 400 water molecules. The starting structures taken were AA N(7) amino symmetric, TT N(3) – O(4) symmetric, GG N(1) – O(6) symmetric, and CC N(3) amino symmetric base pairs. A water spherical cluster with the density of water at room temperature and a radius sphere equal 13.9 Å was used, which corresponds to the most difficult conditions for the formation of stacks (see Abraham, 1982). In spite of such initial conditions, it is shown that during the simulation, each base pair is transformed into a more favorable stacked configuration. The results obtained allow to observe the whole process of convergence for the first time (for more information, visit the Website http://biophys.in.ua/).     

Hydration of Nucleic Acid Bases Studied Using Novel Atom-Atom Potential Functions

Abstract

New simple atom-atom potential functions for simulating behavior of nucleic acids and their fragments in aqueous solutions are suggested. These functions contain terms which are inversely proportional to the first (electrostatics), sixth (or tenth for the atoms, forming hydrogen bonds) and twelfth (repulsion of all the atoms) powers of interatomic distance. For the refinement of the potential function parameters calculations of ice lattice energy, potential energy and configuration of small clusters consisting of water and nucleic acid base molecules as well as Monte Carlo simulation of liquid water were performed. Calculations using new potential functions give rise to more linear hydrogen bonds between water and base molecules than using other potentials. Sites of preferential hydration of five nucleic bases—uracil, thymine, cytosine, guanine and adenine as well as of 6,6,9-trimethyladenine were found. In the most energetically favourable sites water molecule interacts with two adjacent hydrophilic centres of the base. Studies of interaction of the bases with several water molecules showed that water-water interaction play an important role in the arrangement of the nearest to the base water molecules. Hydrophilic centres are connected by “bridges” formed by hydrogen bonded water molecules. The results obtained are consistent with crystallographic and mass-spectrometric data.     

Perturbation of Hydrogen Bonds in the Adenine…Thymine Base Pair by Na+, Mg2+, Ca2+ and NH4 + Cations

Abstract

Perturbation of the hydrogen bonds in the adenine…thymine base pair by Na⁺, Mg²⁺, Ca²⁺ and NH₄ ⁺ cations has been investigated by means of ab initio SCF calculations with the STO-3G basis set. The geometry of adenine…thymine, as well as those of the perturbed pairs were optimized. Approach of any cation to thymine at 06 leads to destabilization of the adenine…thy mine pair; divalent cations (Mg²⁺, Ca²⁺) have a profound effect on the structure of the base pair. The approach of a cation to other available sites (thymine: O2, adenine N1 and N3) leads, on the other hand, to stabilization of the base pair. If a water molecule is placed between the cation and the base pair, the structure and stability of the base pair are changed only negligibly.     

Photochemical demonstration of stacked C.C+ base pairs in a novel DNA secondary structure

The secondary structure of the alternating polydeoxynucleotide sequence poly[d(C-T)] was studied as a function of pH by ultraviolet absorbance and circular dichroism spectroscopy and by the analysis of UV-induced photoproducts. As the pH was lowered, poly[d(C-T)] underwent a conformational transition that was characterized by changes in the long-wavelength region (280-320 nm) of the CD spectrum. These changes have previously been interpreted as evidence for the formation of a core of stacked, protonated C X C+ base pairs in a double-helical complex of poly[d(C-T)], with the thymidyl residues being looped out into the solvent [Gray, D. M., Vaughan, M., Ratliff, R. L., & Hayes, F. N. (1980) Nucleic Acids Res. 8, 3695-3707]. In the present work, poly[d(C-T)] was labeled with [U-14C]cytosine and [methyl-3H]thymine and irradiated at pH values both above and below the conformational transition point (monitored by CD spectroscopy). The distribution of radioactivity in uracil means value of uracil dimers, uracil means value of thymine dimers (the deamination products of cytosine means value of cytosine and cytosine means value of thymine dimers, respectively), and thymine-means value of thymine dimers was then determined. As the pH was decreased, we found an increase in the yield of uracil means value of uracil dimers and a decrease in the yield of uracil means value of thymine dimers, which occurred concomitantly with the change in the CD spectrum.(ABSTRACT TRUNCATED AT 250 WORDS)     

ANALYSIS OF RELATIVE POSITIONS OF RIBONUCLEOTIDE BASES IN A CRYSTAL STRUCTURE OF RIBOSOME

ABSTRACT

Relative positions of bases to bases in a crystal structure of ribosome were analyzed extensively. It was found that there is no clear relation between bases apart more than 15 Å and, thus, the relative location of bases can be analyzed within 15 Å of the reference bases. As for base pairing, major positioning was found to be due to the Watson-Crick type base pairs. Some other positions corresponding to non-Watson-Crick type base pairs were also found in some extents. As for base–base stacking, it was observed that the bases stacked to adenine base are dispersive. It was found that less non-Watson-Crick base pairs was found close to the protein binding site, suggesting that the protein components have a tendency to bind to the regular stem structures. The database of relative location of bases must be useful for improvement of structural determination and structural modeling systems.     

Identification of Escherichia coli Mismatch-specific Uracil DNA Glycosylase as a Robust Xanthine DNA Glycosylase

The gene for the mismatch-specific uracil DNA glycosylase (MUG) was identified in the Escherichia coli genome as a sequence homolog of the human thymine DNA glycosylase with activity against mismatched uracil base pairs. Examination of cell extracts led us to detect a previously unknown xanthine DNA glycosylase (XDG) activity in E. coli. DNA glycosylase assays with purified enzymes indicated the novel XDG activity is attributable to MUG. Here, we report a biochemical characterization of xanthine DNA glycosylase activity in MUG. The wild type MUG possesses more robust activity against xanthine than uracil and is active against all xanthine-containing DNA (C/X, T/X, G/X, A/X and single-stranded X). Analysis of potentials of mean force indicates that the double-stranded xanthine base pairs have a relatively narrow energetic difference in base flipping, whereas the tendency for uracil base flipping follows the order of C/U > G/U > T/U > A/U. Site-directed mutagenesis performed on conserved motifs revealed that Asn-140 and Ser-23 are important determinants for XDG activity in E. coli MUG. Molecular modeling and molecular dynamics simulations reveal distinct hydrogen-bonding patterns in the active site of E. coli MUG that account for the specificity differences between E. coli MUG and human thymine DNA glycosylase as well as that between the wild type MUG and the Asn-140 and Ser-23 mutants. This study underscores the role of the favorable binding interactions in modulating the specificity of DNA glycosylases.     

Alternation of DNA and Solvent Layers in the A Form of d(GGCGCC) Obtainhted by Ethanol Crystallization

Abstract

We have determined by X-ray crystallography the structure of the hexamer duplex d(GGCGCC)₂ in the A-form using ethanol as a precipitant. The same sequence had previously been crystallized in the B-form, but with 2-methyl-2, 4-pentanediol as a precipitant. It appears that ethanol precipitation is a useful method to induce the formation of A-form crystals of DNA. Packing of the molecules in the crystal has unique features: the known interaction of A-DNA duplexes between terminal base-pairs and the minor groove of neighbor molecules is combined with a superstructure consisting in an alternation of DNA layers and solvent layers (water/ions). This organization in layers has been observed before, also with hexamers in the A conformation which crystallize in the same space group (C222 ₁). The solvent layer has a precise thickness, although very few ordered water molecules can be detected. Another feature of this crystal is its large unit cell, which gives rise to an asymmetric unit with three hexamer duplexes. One of the three duplexes is quite different from the other two in several aspects: the number of base pairs per turn, the twist pattern, the mean value of the twist angle and the fact that one terminal base-pair is not stacked as part of the duplex and appears to be disordered. So the variability in conformation of this sequence is remarkable.     

Sequence and conformational effects on imino proton exchange in A.T- and A.U-containing DNA and RNA duplexes

¹H-nmr relaxation has been used to study the effect of sequence and conformation on imino proton exchange in adenine–thymine (A · T) and adenine–uracil (A · U) containing DNA and RNA duplexes. At low temperature, relaxation is caused by dipolar interactions between the imino and the adenine amino and AH2 protons, and at higher temperature, by exchange with the solvent protons. Although room temperature exchange rates vary between 3 and 12s^?1, the exchange activation energies (E_α) are insensitive to changes in the duplex sequence (alternating vs homopolymer duplexes), the conformation (B-form DNA vs A-form RNA), and the identity of the pyrimidine base (thymine vs uracil). The average value of the activation energy for the five duplexes studied, poly[d(A-T)], poly[d(A) · d(T)], poly[d(A-U)], Poly[d(A) · d(U)], and poly[r(A) · r(U)], was 16.8 ± 1.3 kcal/mol. In addition, we find that the average E_α for the A.T base pairs in a 43-base-pair restriction fragment is 16.4 ± 1.0 kcal/mol. This result is to be contrasted with the observation that the E_α of cytosine-containing duplexes depends on the sequence, conformation, and substituent groups on the purine and pyrimidine bases. Taken together, the data indicate that there is a common low-energy pathway for the escape of the thymine (uracil) imino protons from the double helix. The absolute values of the exchange rates in the simple sequence polymers are typically 3–10 times faster than in DNAs containing both A · T and G · C base pairs.     

Crystal structure of a Z-DNA fragment containing thymine/2-aminoadenine base pairs

The Z-DNA structure has been shown to form in two crystals made from self-complementary DNA hexamers d(CGTDCG) and d(CDCGTG) which contain thymine/2-aminoadenine (TD) base pairs. The latter structure has been solved and refined to 1.3 A resolution and it shows only small conformational changes due to the introduction of the TD base pairs in comparison with the structure of d(CG)3. Spectroscopic studies with these compounds demonstrate that DNA molecules containing 2-aminoadenine residues form Z-DNA slightly more easily than do those containing adenine nucleotides, but not as readily as the parent sequence containing only guanine-cytosine base pairs.     

The effects of substituted pyrimidines in DNAs on cleavage by sequence-specific endonucleases.

The rates of cleavage of DNAs containing substituents at position 5 of thymine or cytosine have been measured for a variety of sequence-specific endonucleases, so as to determine which features in the DNA sequence are being probed. Phage phi e DNA fully substituted with 5-hydroxymethyluracil is cleaved more slowly by enzymes whose recognition sequences contain A-T base pairs than are DNAs containing thymine, but both types of DNA are cleaved at similar rates by enzymes recognizing sequences composed only of G-C base pairs. Phage PBS2 DNA with uracil completely substituted for thymine is cleaved slowly by several enzymes which recognize sequences containing A-T base pairs (endonucleases Hpa I, HindII, and HindIII), while the rates of cleavage by other enzymes (endonucleases EcoRI and BamHI) are not affected. Phage lambda- and P22 DNAs containing 5-bromouracil are cleaved more slowly by several enzymes (endonucleases HindIII, Hpa I, BamHI) than are thymine-containing DNAs. Enzymes that recognize sequence isomers with the composition G:C:2A:2T (endonucleases EcoRI, Hpa I, HindIII) are not equally affected by substitution at position 5 of thymine, suggesting that they differ in their contacts with A-T base pairs. DNA containing glucosylated 5-hydroxymethylcytosine in place of cytosine is resistant to cleavage by all the endonucleases examined.     

Base Sequence Effects in Double Helical DNA. I. Potential Energy Estimates of Local Base Morphology

Abstract

A series of potential energy calculations have been carried out to estimate base sequence dependent structural differences in B-DNA. Attention has been focused on the simplest dimeric fragments that can be used to build long chains, computing the energy as a function of the orientation and displacement of the 16 possible base pair combinations within the double helix. Calculations have been performed, for simplicity, on free base pairs rather than complete nucleotide units. Conformational preferences and relative flexibilities are reported for various combinations of the roll, tilt, twist, lateral displacement, and propeller twist of individual residues. The predictions are compared with relevant experimental measures of conformation and flexibility, where available. The energy surfaces are found to fit into two distinct categories, some dimer duplexes preferring to bend in a symmetric fashion and others in a skewed manner. The effects of common chemical substitutions (uracil for thymine, 5-methyl cytosine for cytosine, and hypoxanthine for guanine) on the preferred arrangements of neighboring residues are also examined, and the interactions of the sugar-phosphate backbone are included in selected cases. As a first approximation, long range interactions between more distant neighbors, which may affect the local chain configuration, are ignored. A rotational isomeric state scheme is developed to describe the average configurations of individual dimers and is used to develop a static picture of overall double helical structure. The ability of the energetic scheme to account for documented examples of intrinsic B-DNA curvature is presented, and some new predictions of sequence directed chain bending are offered.     

G.T wobble base-pairing in Z-DNA at 1.0 A atomic resolution: the crystal structure of d(CGCGTG).

The DNA oligomer d(CGCGTG) crystallizes as a Z-DNA double helix containing two guanine-thymine base pair mismatches of the wobble type. The crystal diffracts to 1 A resolution and the structure has been solved and refined. At this resolution, a large amount of information is revealed about the organization of the water molecules in the lattice generally and more specifically around the wobble base pairs. By comparing this structure with the analogous high resolution structure of d(CGCGCG) we can visualize the structural changes as well as the reorganization of the solvent molecules associated with wobble base pairing. There is only a small distortion of the Z-DNA backbone resulting from introduction of the GT mismatched base pairs. The water molecules cluster around the wobble base pair taking up all of the hydrogen bonding capabilities of the bases due to wobble pairing. These bridging water molecules serve to stabilize the base-base interaction and, thus, may be generally important for base mispairing either in DNA or in RNA molecules.     

Theoretical study on the binding mechanism between N6-methyladenine and natural DNA bases

N6-methyladenine (m⁶A) is a rare base naturally occurring in DNA. It is different from the base adenine due to its N-CH₃. Therefore, the base not only pairs with thymine, but also with other DNA bases (cytosine, adenine and guanine). In this work, Møller-Plesset second-order (MP2) method has been used to investigate the binding mechanism between m⁶A and natural DNA bases in gas phase and in aqueous solution. The results show that N-CH₃ changed the way of N6-methyladenine binding to natural DNA bases. The binding style significantly influences the stability of base pairs. The trans-m⁶A:G and trans-m⁶A:C conformers are the most stable among all the base pairs. The existence of solvent can remarkably reduce the stability of the base pairs, and the DNA bases prefer pairing with trans-m⁶A to cis-m⁶A. Besides, the properties of these hydrogen bonds have been analyzed by atom in molecules (AIM) theory, natural bond orbital (NBO) analysis and Wiberg bond indexes (WBI). In addition, pairing with m⁶A decreases the binding energies compared to the normal Watson-Crick base pairs, it may explain the instability of the N6 site methylated DNA in theory.

Designing of MIP based QCM sensor having thymine recognition sites based on biomimicking DNA approach

Quartz crystal microbalance (QCM) sensors coated with molecular imprinted polymers (MIP) have been developed for the determination of thymine. In this method, methacryloylamidoadenine (MA-Ade) have used as a new monomer and thymine template for inspiration of DNA nucleobases interaction. The thymine can be simultaneously hydrogen binding to MA-Ade and fit into the shape-selective cavities. Thus, the interaction between nucleobases has an effect on the binding ability of the QCM sensors. The binding affinity of the thymine imprinted sensors has investigated by using the Langmuir isotherm. The thymine imprinted QCM electrodes have shown homogeneous binding sites for thymine (K_a: 1.0 × 10⁵ M⁻¹) while heterogeneous binding sites for uracil. On the other hand, recognition selectivity of the QCM sensor based on thymine imprinted polymer toward to uracil, ssDNA and ssRNA has been reported in this work.     

A possible prebiotic synthesis of thymine: uracil-formaldehyde-formic acid reaction.

When uracil is reacted with formaldehyde and formic acid in dilute aqueous solutions at 100-140 degrees C, 5-hydroxymethyluracil (5-HMU), methylenebiuracil (MBU) and thymine are formed. It has been shown that 5-HMU is an intermediate in the formation of MBU and thymine. In the presence of formic acid, 5-HMU gives MBU, thymine and in some cases uracil. The formation of thymine is generally favoured under acidic conditions, although small amounts of this base could also be obtained when the reactions were carried out under mildly basic conditions. A hydride ion transfer mechanism is suggested for some of these reactions. These results have relevance to the formation of thymine under prebiotic conditions.     

Hydration of uracil and thymine methylderivatives: a Monte Carlo simulation

The simulation performed shows that under methylation of uracil and thymine NH-groups the interaction energy between a base and water (Uwb) is increased. It is also detected that the increase in this energy was observed in the 1st and the 3rd sectors. These conclusions do not confirm the assumption made in the literature on the character of an interaction between methylated bases and water. According to this assumption, when the NH-groups are methylated, the energy of Uwb in these sectors decreases as a result of the van der Waals interactions between a methyl group and water, whose energy compensates the increase in the Uwb energy due to the breaking of an H-bond. Regularity of water molecules near a hydrophobic group under the hydration of polar molecules is detected for the first time.