Monte Carlo Simulation of Hydration of the Nucleic Acid Fragments

Abstract

Systems containing a base or a base pair and 25 water molecules, as well as a helical stack and 30 water molecules per base pair, have been simulated. Changes in the base hydration shell structure, after the bases have been included into the pair and then into the base pair stack are discussed. Hydration shells of several configurations of the base pair stacks are discussed. Probabilities of formation of the hydrogen-bonded bridges of 1, 2 and 3 water molecules between hydrophilic centres have been estimated. The hydration shell structure was shown to depend on the nature of the base pair and on the stack configuration, while dependence of the global hydration shell characteristics on the stack configuration has been proved to be rather slight. The most typical structural elements of hydration shells, in the glycosidic (minor in B-like conformation) and non-glycosidic (major) grooves, for different configurations of AU and GC stacks, have been found and discussed. The number of hydrogen bonds between water molecules and bases per water molecule was shown to change upon transformation of the stack from A to B configuration. This result is discussed in connection with the reasons for B to A conformational transition and the concept of “water economy”. Hydration shell patterns of NH₂-groups of AU and GC helical stacks differ significantly.     

Use of “Loss-of-Contact” Substitutions to Identify Residues Involved in an Amino Acid-Base Pair Contact: Effect of Substitution of Gln18 of Lac Repressor by Gly,Ser, and Leu

Abstract

A procedure to identify which base pair of lac operator (lacO) a suspected contacting amino acid of Lac repressor (LacR) interacts with is presented. The procedure is to eliminate the ability of the amino acid under study to contact DNA, and then to determine at which base pair—if any—specificity is eliminated. To implement this procedure, four sets of Escherichia coli K-12 strains have been constructed. These strains permit: (i) the substitution of a selected amino acid of LacR by, respectively, Gly, Ser, Leu, or Gln, and (ii) the analysis of the specificity of the resulting substituted LacR with respect to base pairs 5,6,7,8,9, and 10 of lacO. This procedure has been applied to Gln18 of LacR. The preliminary data indicate that LacR(Gln18?Gly) is unable to distinguish between the O⁺ base pair G:C and the O^c base pair T:A at position 7 of lacO (K_Doc/K_DO ⁺ = 0.93). In contrast, LacR(Glnl8?Gly) discriminates O⁺ from O^c by a factor of 13 to 23 at each other position. The same qualitative pattern of results was obtained with LacR(Glnl8?Ser) and LacR(Gln18?Leu). Therefore, I propose that Gln18 contacts base pair 7 of lacO. This proposal is consistent with the contact predicted in Ebright, R. in Protein Structure, Folding, and Design. D. Oxender ed., Alan R. Liss, New York (1985), in press.     

Zero erucic acid trait of rapeseed (<Emphasis Type="Italic">Brassica napus</Emphasis> L.) results from a deletion of four base pairs in the <Emphasis Type="Italic">fatty acid elongase 1</Emphasis> gene

The fatty acid elongase 1 (FAE1) gene is a key gene in the erucic acid biosynthesis in rapeseed. The complete coding sequences of the FAE1 gene were isolated separately from eight high and zero erucic acid rapeseed cultivars (Brassica napus L.). A four base pair deletion between T1366 and G1369 in the FAE1 gene was found in a number of the cultivars, which leads to a frameshift mutation and a premature stop of the translation after the 466th amino acid residue. This deletion was predominantly found in the C-genome and rarely in the A-genome of B. napus. Expression of the gene isoforms with the four base pair deletion in a yeast system generated truncated proteins with no enzymatic activity and could not produce very long chain fatty acids as the control with an intact FAE1 gene did in yeast cells. In the developing rape seeds the FAE1 gene isoforms with the four base pair deletion were transcribed normally but failed to translate proteins to form a functional complex. The four base pair deletion proved to be a mutation responsible for the low erucic acid trait in rapeseed and independent from the point mutation reported by Han et al. (Plant Mol Biol 46:229–239, 2001). Gang Wu, Yuhua Wu contribute equally to this article.     

The formation of catalytically competent enzyme–substrate complex is not a bottleneck in lesion excision by human alkyladenine DNA glycosylase

Human alkyladenine DNA glycosylase (AAG) protects DNA from alkylated and deaminated purine lesions. AAG flips out the damaged nucleotide from the double helix of DNA and catalyzes the hydrolysis of the N-glycosidic bond to release the damaged base. To understand better, how the step of nucleotide eversion influences the overall catalytic process, we performed a pre-steady-state kinetic analysis of AAG interaction with specific DNA-substrates, 13-base pair duplexes containing in the 7th position 1-N6-ethenoadenine (εA), hypoxanthine (Hx), and the stable product analogue tetrahydrofuran (F). The combination of the fluorescence of tryptophan, 2-aminopurine, and 1-N6-ethenoadenine was used to record conformational changes of the enzyme and DNA during the processes of DNA lesion recognition, damaged base eversion, excision of the N-glycosidic bond, and product release. The thermal stability of the duplexes characterized by the temperature of melting, T_m, and the rates of spontaneous opening of individual nucleotide base pairs were determined by NMR spectroscopy. The data show that the relative thermal stability of duplexes containing a particular base pair in position 7, (T_m(F/T)?T_m(εA/T)?T_m(Hx/T)?T_m(A/T)) correlates with the rate of reversible spontaneous opening of the base pair. However, in contrast to that, the catalytic lesion excision rate is two orders of magnitude higher for Hx-containing substrates than for substrates containing εA, proving that catalytic activity is not correlated with the stability of the damaged base pair. Our study reveals that the formation of the catalytically competent enzyme–substrate complex is not the bottleneck controlling the catalytic activity of AAG.     

Variable Temperature Infrared Spectroscopy of Cytosine-Guanine Base Pairs: Tautomerism Versus Polarization

Abstract

This report describes an infrared (IR) spectroscopic study of a model cytosine - guanine base pair. This base pair is part of a self-consistent experimental system based on lipophilic ribose derivatives of cytidine (C), guanosine (G) and O⁶-methylguanosine (O⁶MeG) that are soluble in non-aqueous, low dielectric solvents at appreciable concentrations. Previous experiments on this system have revealed different rotation dynamics for the amino bonds within the CG base pair, an observation that could be explained by the presence of rare tautomers (P.O. Lowdin, Reviews of Modern Physics 35, 724 (1963)), or by mutual polarization of the base pairs (L.D. Williams, N.G. Williams and B.R. Shaw, J. Am. Chem. Soc. 112, 829 (1990)). The IR spectra in the OH and NH stretching region indicate formation of hydrogen-bonded CG base pairs and self associates in 1,2-dichlorobenzene over a temperature range from 10 to 290K. Changes in the lineshapes and intensities of the IR bands with temperature correlate with phase transitions of the solvent but no evidence is seen for an OH stretching band that would indicate the formation of hydroxyl tautomers within base pairs. Similarly, the relative intensities of the C=O stretching bands of CG in cyclohexane solution remain constant over this same temperature range, confirming that within the base pair, the tautomeric states of the bases remain essentially unperturbed in the 2-amino/6-keto form of G and the 2-keto/4-amino form of C. The spectra of O⁶-MeG aid in the band assignments, since this molecule is frozen in an equivalent of the 2-amino/6-hydroxyl tautomer, but without the OH group and its associated stretching band. We conclude that the probability of tautomerism does not appear to be sufficient to explain the different rotation dynamics for the two amino bonds of the CG base pair. Rather it is argued that mutual polarization within the base pair, which would increase the bond order of the amino bond of C within the base pair, can explain the results without the formation of unconventional tautomers.     

The Role of ions in the acid melting of DNA

The use of a simple chemical equilibrium and of the law of mass action leads to correct prediction when applied to the study of the role of the counterions in the process of acid melting of the DNA molecule. This approach allows an estimate to be given of the enthalpy variation per base pair associated with this process. Experiments have been carried out to test (a) the linearity of the dependence of pH_m vs. 1/T; (b) the stabilizing effect of Na⁺ concentration. The enthalphy variation per base pair, deduced from the slope of the pH_m vs. 1/T plot, and from the number of H⁺ bound per base pair, is in agreement with direct measurement.^6,8     

Pyrazolo[3,4-d]pyrimidine Nucleic Acids: Adjustment of the dA-dT to the dG-dC Base Pair Stability

Abstract

The pyrazolo[3,4-d]pyrimidine-4,6-diamine nucleosides 2b-d stabilize the dA-dT base pair significantly when the dA-residue is replaced. Oligonucleotide duplexes incorporating 2b-d show a 4–6°C T _m increase per modification. The 7-bromo compound 2b harmonizes the stability of the dA-dT vs. the dG-dC pair. According to this the stability of such duplexes depends no longer on the base pair composition of a DNA molecule.     

Why the tautomerization of the G·C Watson–Crick base pair via the DPT does not cause point mutations during DNA replication? QM and QTAIM comprehensive analysis

The ground-state tautomerization of the G·C Watson–Crick base pair by the double proton transfer (DPT) was comprehensively studied in vacuo and in the continuum with a low dielectric constant (??=?4), corresponding to a hydrophobic interface of protein–nucleic acid interactions, using DFT and MP2 levels of quantum-mechanical (QM) theory and quantum theory “Atoms in molecules” (QTAIM). Based on the sweeps of the electron-topological, geometric, polar, and energetic parameters, which describe the course of the G·C???G^?·C^? tautomerization (mutagenic tautomers of the G and C bases are marked with an asterisk) through the DPT along the intrinsic reaction coordinate (IRC), it was proved that it is, strictly speaking, a concerted asynchronous process both at the DFT and MP2 levels of theory, in which protons move with a small time gap in vacuum, while this time delay noticeably increases in the continuum with ??=?4. It was demonstrated using the conductor-like polarizable continuum model (CPCM) that the continuum with ??=?4 does not qualitatively affect the course of the tautomerization reaction. The DPT in the G·C Watson–Crick base pair occurs without any intermediates both in vacuum and in the continuum with ??=?4 at the DFT/MP2 levels of theory. The nine key points along the IRC of the G·C base pair tautomerization, which could be considered as electron-topological “fingerprints” of a concerted asynchronous process of the tautomerization via the DPT, have been identified and fully characterized. These key points have been used to define the reactant, transition state, and product regions of the DPT reaction in the G·C base pair. Analysis of the energetic characteristics of the H-bonds allows us to arrive at a definite conclusion that the middle N1H?N3/N3H?N1 and the lower N2H?O2/N2H?O2 parallel H-bonds in the G·C/G^?·C^? base pairs, respectively, are anticooperative, that is, the strengthening of the middle H-bond is accompanied by the weakening of the lower H-bond. At that point, the upper N4H?O6 and O6H?N4 H-bonds in the G·C and G^?·C^? base pairs, respectively, remain constant at the changes of the middle and the lower H-bonds at the beginning and at the ending of the G·C???G^?·C^? tautomerization. Aiming to answer the question posed in the title of the article, we established that the G^?·C^? Löwdin’s base pair satisfies all the requirements necessary to cause point mutations in DNA except its lifetime, which is much less than the period of time required for the replication machinery to forcibly dissociate a base pair into the monomers (several ns) during DNA replication. So, from the physicochemical point of view, the G^?·C^? Löwdin’s base pair cannot be considered as a source of point mutations arising during DNA replication.     

Effects of osmolytes on stable UUCG tetraloops and their preference for a CG closing base pair

Osmolytes have the potential to affect the stability of secondary structure motifs and alter preferences for conserved nucleic acid sequences in the cell. To contribute to the understanding of the in vivo function of RNA we observed the effects of different classes of osmolytes on the UNCG tetraloop motif. UNCG tetraloops are the most common and stable of the RNA tetraloops and are nucleation sites for RNA folding. They also have a significant thermodynamic preference for a CG closing base pair. The thermal denaturation of model hairpins containing UUCG loops was monitored using UV-Vis spectroscopy in the presence of osmolytes with different chemical properties. Interestingly, all of the osmolytes tested destabilized the hairpins, but all had little effect on the thermodynamic preference for a CG base pair, except for polyethylene glycol (PEG) 200. PEG 200 destabilized the loop with the CG closing base pair relative to the loop with a GC closing base pair. The destabilization was linear with increasing concentrations of PEG 200, and the slope of this relationship was not perturbed by changes in the hairpin stem outside of the closing pair. This result suggests that in the presence of PEG 200, the UUCG loop with a GC closing base pair may retain some preferential interactions with the cosolute that are lost in the presence of the CG closing base pair. These results reveal that relatively small structural changes may influence how osmolytes tune the stability, and thus the function of a secondary structure motif in vivo.     

Molecular-Mechanical studies of the mismatched base analogs of d(CGCGAATTCGCG)2:d(CGTGAATTCGCG)2, d(CGAGAATTCGCG)2, d(CGCGAATTCACG)2, d(CGCGAATTCTCG)2, and d(CGCAGAATTCGCG)·d(CGCGAATTCGCG)

Molecular-mechanical simulations have been carried out on “mismatched base” analogs of the DNA double-helical structure d(CGCGAATTCGCG)₂, in which the base pairs CG at the 3 and 10 positions have been replaced by CA, AG, TC, and TG base pairs, as well as an insertion analog in which an extra adenine has been incorporated into one strand of the above structure between bases 3 and 4. The results of these simulations (calculated relative stabilities, structures, and nmr ring-current shifts) have been compared with calorimetric and nmr data. The calculated relative stabilities of the double-helical parent dodecamer and the various “wobble” base pairs qualitatively correlate with the experimental melting temperatures. The base-pairing structure for the GT wobble pair is in agreement with that previously determined from nmr experiments. For the GA base pair, the structure with both bases anti has a slightly more favorable energy from base pairing and stacking than a structure with non-Watson-Crick H-bonding with adenine syn, in agreement with nmr experiments. The CA wobble base is calculated to favor an adenine 6NH₂ …? cytosine N3 H-bond over cytosine 4NH₂ …? adenine N1, again, in agreement with nmr experiments. There is no definitive experimental data on the TC base pair, but the existence of (somewhat long and weak) H-bonds involving cytosine 4NH₂ …? thymine 4CO and cytosine N3 …? thymine HN3 seems reasonable. We find a structure in which the extra adenine base of the insertion analogs sits “inside” the double helix.     

Structure of an 11-mer DNA Duplex Containing the Carbocyclic Nucleotide Analog: 2′-Deoxyaristeromycin

Abstract

2′-deoxyaristeromycin (dAr) is a nucleoside analogue that is resistant to the action of DNA glycosylases. High-resolution NMR spectroscopy and molecular dynamics simulations were used to determine the three-dimensional structure of an 11-mer DNA containing a single dAr?T base pair at its center. Analysis of the spectra revealed the existence of a right-handed duplex in solution, stabilized by Watson-Crick hydrogen bonding and base-stacking interactions. The carbocyclic sugar adopted a C1′-exo conformation and sugars of the 3′-flanking base pair had puckers in the O4′—endo range. The dAr?T base pair was mildly propeller twisted, and the dAr analogue showed a positive roll with the 3′-flanking base. Our findings indicate that the observed resistance of dAr-containing oligodeoxynucleotides to the catalytic action of DNA glycosylases relates to its electronic properties rather than structure, and validate the use of dAr and related carbocyclic nucleoside analogues for biological and structure/function relationship studies.     

Investigation of the base-pairing structure of the anticodon hairpin from E. coli initiator tRNA by high-resolution nmr

The high-resolution nmr spectrum of the anticodon hairpin from E. coli tRNA^fMet has been obtained at a number of different temperatures. The positions of the resonances from interior Watson-Crick base pairs are well accounted for (within 0.1 ppm) by a semi-empirical ring current shift theory, but the terminal base pairs are susceptible to the exact orientation of adjacent bases in single-stranded regions. From a careful examination of the exact way in which resonances disappear at elevated temperatures, we conclude that melting in the nmr experiments occurs when the lifetime of a base pair is reduced to several milliseconds. On the basis of these experiments we are able to assign an nmr T_m to each individual base pair and these should be useful in interpreting the melting behavior of the intact molecule. An “extra” resonance is observed at ～11.3 ppm and, on the basis of its position and temperature sensitivity, it is tentatively assigned to the ring nitrogen proton of a “protected” U residue in the anticodon loop. A strong preference for stacking of a nonbase-paired A residue on an adjacent GC base pair is observed even at temperatures in excess of 52°C.     

Investigation of the thermal unfolding of secondary and tertiary structure in E. coli tRNAfMet by high-resolution nmr

The high-resolution (300 MHz) proton nmr spectrum of E. coli tRNA^fMet has been examined in 0.17M NaCl, with and without Mg²⁺, and at various temperatures. In light of recent studies of other E. coli tRNA and fragments of tRNA^fMet, some low field (11–15 ppm) resonances previously assigned to secondary structure base pairs are reassigned to a tertiary structure A₁₄–S⁴U₈ base pair and a protected uridine residue in the anticodon loop. These two resonances and other low field resonances which are assigned to secondary structure base pairs are used to monitor the thermal unfolding of the molecule. In the absence of Mg²⁺ the tertiary structure base pair is present only to ～45°C, but in the presence of Mg²⁺ it remains until at least 70°C. Analysis of the temperature dependence of other low field resonances indicates that the melting of the dihydrouridine stem occurs more or less simultaneously with the loss of tertiary structure. The observation of the resonance from the A₁₄–S⁴U₈ base pair proves that tertiary structure is present in this molecule below 40°C, even in the absence of Mg²⁺.     

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     

Solution Structure of a GAAG Tetraloop in Helix 6 of SRP RNA from Pyrococcus furiosus

The NMR structure of a 12-mer RNA derived from the helix 6 of SRP RNA from Pyrococcus furiosus, whose loop-closing base pair is U:G, was determined, and the structural and thermodynamic properties of the RNA were compared with those of a mutant RNA with the C:G closing base pair. Although the structures of the two RNAs are similar to each other and adopt the GNRR motif, the conformational stabilities are significantly different to each other. It was suggested that weaker stacking interaction of the GAAG loop with the U:G closing base pair in 12-mer RNA causes the lower conformational stability.     

Silver(I)-mediated Hoogsteen-type base pairs

Metal-mediated Hoogsteen-type base pairs are useful for the construction of DNA duplexes containing contiguous stretches of metal ions along the helical axis. To fine-tune the stability of such base pairs and the selectivity toward different metal ions, the availability of a selection of artificial nucleobases is highly desirable. In this study, we follow a theoretical approach utilizing dispersion-corrected density functional methods to evaluate a variety of artificial nucleobases as candidates for metal-mediated Hoogsteen-type base pairs. We focus on silver(I)-mediated Hoogsteen- and reverse Hoogsteen-type base pairs formed between 1-deaza- and 1,3-dideazapurine-derived nucleobases, respectively, and cytosine. Apart from two coordinative bonds, these base pairs are stabilized by a hydrogen bond. We elucidate the impact of different substituents at the C6 position and the presence or absence of an endocyclic N3 nitrogen atom on the overall stability of a base pair and concomitantly on the strength of the hydrogen and coordinative bonds. All artificial base pairs investigated in this study are less stable than the experimentally established benchmark base pair C–Ag⁺–G. The base pair formed from 1,3-dideaza-6-methoxypurine is isoenergetic to the experimentally observed C–Ag⁺–C base pair. This makes 1,3-dideaza-6-methoxypurine a promising candidate for the use as an artificial nucleobase in DNA.     

Structural basis for a novel mechanism of DNA bridging and alignment in eukaryotic DSB DNA repair

Eukaryotic DNA polymerase mu of the PolX family can promote the association of the two 3′‐protruding ends of a DNA double‐strand break (DSB) being repaired (DNA synapsis) even in the absence of the core non‐homologous end‐joining (NHEJ) machinery. Here, we show that terminal deoxynucleotidyltransferase (TdT), a closely related PolX involved in V(D)J recombination, has the same property. We solved its crystal structure with an annealed DNA synapsis containing one micro‐homology (MH) base pair and one nascent base pair. This structure reveals how the N‐terminal domain and Loop 1 of Tdt cooperate for bridging the two DNA ends, providing a templating base in trans and limiting the MH search region to only two base pairs. A network of ordered water molecules is proposed to assist the incorporation of any nucleotide independently of the in trans templating base. These data are consistent with a recent model that explains the statistics of sequences synthesized in vivo by Tdt based solely on this dinucleotide step. Site‐directed mutagenesis and functional tests suggest that this structural model is also valid for Pol mu during NHEJ.     

Nicotiana chloroplast genome: X. Correlation between the DNA sequences and the isoelectric focusing patterns of the LS of Rubisco

Summary Comparison of the DNA sequences of the rbcL gene from three Nicotiana species reveals a high degree of homology among the 1431 bp in the coding region. Only eight base pair differences are observed between N. otophora and N. tabacum, and between N. otophora and N. acuminata. Four base pair differences are observed between N. acuminata and N. tabacum. Most changes are in the third position of the codon resulting in only two amino acid alterations when N. otophora and N. acuminata are compared with N. tabacum. Evidence is presented demonstrating that the amino acid compositions of the LS derived from the DNA sequence are related to the IEF cluster pattern. A single charged residue is responsible for the difference in cluster pattern.     

The paramutagenic line niv-44 has a 5 kb insert, Tam 2, in the chalcone synthase gene of Antirrhinum majus

Summary Several alleles of the nivea locus of Antirrhinum majus, both stable and unstable, have been characterised genetically (Harrison and Carpenter 1973 a, b). In this work the niv-44 allele is characterised at the molecular level. It contains a 5kb insertion element, Tam 2, which has 14 base pair inverted repeats. There is a three base pair duplication at the target site, which is at the first intron-exon boundary of the chalcone synthase gene. Tam 2 homologous sequences are present in multiple copies in several A. majus lines, including niv-53, and most have at least a 2.9 kb sequence in common with the copy at the chalcone synthase gene. Possible reasons for the apparent stability of the niv-44 allele and molecular explanations for the role of this allele in paramutation in A. majus are discussed.Dedicated to Professor Georg Melchers to celebrate his 50-year association with the journal     

Theoretical Investigation of the Molecular Structure of the πκ DNA Base Pair

Abstract

The structure of the nonclassical πκ base pair (7–methyl-oxoformycin … 2,4-diaminopyrimidine) was studied at the ab initio Hartree-Fock (HF) and MP2 levels using the 6–31G* and 6–31G** basis sets. The πκ base pair is bound by three parallel hydrogen bonds with the donor-acceptor-donor recognition pattern. Recently, these bases were proposed as an extension of the genetic alphabet from four to six letters (Piccirilli et al. Nature 343, 33(1990)). By the HF/6- 31G* method with full geometry optimization we calculated the 12 degree propeller twist for the minimum energy structure of this complex. The linearity of hydrogen bonds is preserved in the twisted structure by virtue of the pyramidal arrangement of the κ-base amino groups. The rings of both the π and κ molecules remain nearly planar. This nonplanar structure of the πκ base pair is only 0.1 kcal/mol more stable than the planar (Cs) conformation. The HF/6- 31G* level gas-phase interaction energy of πκ (—13.5 kcal/mol) calculated by us turned out to be nearly the same as the interaction energy obtained previously for the adenine-thymine base pair (—13.4 kcal/mol) at the same computational level. The inclusion of p-polarization functions on hydrogens, electron correlation effects (MP2/6–31G** level), and the correction for the basis set superposition error (BSSE) increase this energy to -14.0 kcal/mol.