Simulation of DNA bases in water: Comparison of the Monte Carlo algorithm with molecular mechanics force fields

The interaction between the nucleic acid bases and solvent molecules has an important effect in various biochemical processes. We have calculated total energy and free energy of the solvation of DNA bases in water by Monte Carlo simulation. Adenine, guanine, cytosine, and thymine were first optimized in the gas phase and then placed in a cubic box of water. We have used the TIP3 model for water and OPLS for the nucleic acid bases. The canonical (T, V, N) ensemble at 25 degrees C and Metropolis sampling technique have been used. Good agreement with other available computational data was obtained. Radial distribution functions of water around each site of adenine, guanine, cytosine, and thymine have been computed and the results have shown the ability of the sites for hydrogen bonding and other interactions. The computations have shown that guanine has the highest value of solvation free energy and N7 and N6 in adenine and guanine, N3 in cytosine, and N3 and O4 in thymine have the largest radial distribution function. Monte Carlo simulation has also been performed using the CHARMM program under the same conditions, and the results of two procedures are compared.     

Simulation of DNA bases in water: Comparison of the Monte Carlo algorithm with molecular mechanics force fields

The interaction between the nucleic acid bases and solvent molecules has an important effect in various biochemical processes. We have calculated total energy and free energy of the solvation of DNA bases in water by Monte Carlo simulation. Adenine, guanine, cytosine, and thymine were first optimized in the gas phase and then placed in a cubic box of water. We have used the TIP3 model for water and OPLS for the nucleic acid bases. The canonical (T, V, N) ensemble at 25°C and Metropolis sampling technique have been used. Good agreement with other available computational data was obtained. Radial distribution functions of water around each site of adenine, guanine, cytosine, and thymine have been computed and the results have shown the ability of the sites for hydrogen bonding and other interactions. The computations have shown that guanine has the highest value of solvation free energy and N7 and N6 in adenine and guanine, N3 in cytosine, and N3 and O4 in thymine have the largest radial distribution function. Monte Carlo simulation has also been performed using the CHARMM program under the same conditions, and the results of two procedures are compared.     

Selective effect of ionization of nitrogen bases on the thermostability of AT- and GC-pairs of DNA

The ionization specific influence of nitrogen bases on thermostability of AT- and GC-pairs of DNA has been investigated by the method of DNAs melting temperature analysis. It has been shown that the change of temperature interval of DNA helix-coil transition when changing pH environment is due to specific ionization of AT- and GC-pairs of nitrogen bases.     

Raman spectroscopic measurement of base stacking in solutions of adenosine, AMP, ATP, and oligoadenylates

Measurements of the colligative properties of nucleosides and their derivatives have shown that bases form transient aggregates in solution [Ts'o (1967) J. Am. Chem. Soc. 89, 3612-3622]. Aggregation of nucleotides cannot be measured by osmometry due to the presence of counterions. Sedimentation measurements are difficult to obtain and have been complicated by differences in pH [Ferguson et al. (1974) Biophys. Chem. 1, 325-337]. Raman studies of oligonucleotides have shown that the intensities due to base vibrational modes depend on the extent of base stacking, but this dependence has not been quantitated. We have measured this dependence by relating changes in the Raman spectra of nucleotides and nucleosides with previous measurements of colligative properties. Visible Raman spectra of ATP, AMP, and adenosine, taken over a range of concentrations from 1 to 1000 mM, show that the peak intensity ratio (I1305 + I1380)/I1340 varies linearly with the log of the concentration for all three bases. This concentration-dependent change correlates with published molal osmotic coefficient data for functionally similar bases with a correlation coefficient of 0.99. In contrast, UV resonance Raman spectra of the same bases show changes that vary linearly with concentration.     

Molecular modelling of the 3-D structure of RNA tetraloops with different nucleotide sequences.

One surprisingly common element of RNA secondary structure consists of a hairpin capped by a four-base loop (or the tetraloop). Recently the 3-D structures of two RNA-tetraloops have been determined by NMR-studies. Both structures have a similar architecture: the first and the last bases of the loop form a hydrogen bonded pair which is stacked on the stem base pair. We have analysed the ability of tetraloops, with the other combinations of the first and the fourth bases, to adopt such a 'diloop' conformation using computer modelling. The analysis has shown that the 'diloop' conformation has many covalent and steric constraints which give a possibility for reliable structural predictions. As a result, a set of the tetraloop 3-D structures in which hydrogen bonded pairing of the first and the last bases does not cause covalent and steric hindrances has been selected. In most cases several predicted 3-D structures corresponded to one tetraloop sequence. Taking into consideration the folding pathway of RNA hairpins we have resolved this ambiguity and predicted the most probable 3-D structure for every possible nucleotide sequence of the tetraloop. On the basis of these results a conclusion has been drawn on the possible reasons of the tetraloop phylogenetic preference.     

Accumulation of U14 small nuclear RNA in Saccharomyces cerevisiae requires box C, box D, and a 5'', 3'' terminal stem.

U14 is one of several nucleolar small nuclear RNAs required for normal processing of rRNA. Functional mapping of U14 from Saccharomyces cerevisiae has yielded a number of mutants defective in U14 accumulation or function. In this study, we have further defined three structural elements required for U14 accumulation. The essential elements include the U14-conserved box C and box D sequences and a 5', 3' terminal stem. The box elements are coconserved among several nucleolar small nuclear RNAs and have been implicated in binding of the protein fibrillarin. New mutational results show that the first GA bases of the box C sequence UGAUGA are essential, and two vital bases in box D have also been identified. An intragenic suppressor of a lethal box C mutant has been isolated and shown to contain a new box C-like PyGAUG sequence two bases upstream of normal box C. The importance of the terminal stem was confirmed from new compensatory base changes and the finding that accumulation defects in the box elements can be complemented by extending the terminal stem. The results suggest that the observed defects in accumulation reflect U14 instability and that protein binding to one or more of these elements is required for metabolic stability.     

Substrate specificity of Aspergillus clavatus intracellular acid RNAase]

Substrate specificity of intracellular acid RNAse from Aspergillus clavatus, has been studied using different RNAs, synthetic polynucleotides and diribonucleoside monophosphates as substrates. The enzyme was shown to be a RNAse, non-specific to the chemical nature of bases adjacent to the disrupted phosphodiesther bonds in the molecules of RNA. It has been demonstrated that the order of nucleotide release from RNA coincides with the order of weakening of the enzyme binding to substrates XpY, depending on the base X. Purine bases increase substrates XpY binding with the enzyme and hamper their splitting. The effect of pyrimidine bases on adsorption and catalytic functions of the enzyme is contrary to that of purine bases cited above.     

Radiation-induced DNA damage and its repair

Application of modern methods of organic chemistry and recombinant DNA technologies has provided new insights in the field of DNA radiation damage and its repair. An overview of the chemical nature of the lesions inflicted on DNA by ionizing radiation is presented. The structures of 29 different DNA modified base or sugar residues are shown in comprehensive formation schemes. A fraction of radiation-induced modified bases is spontaneously released from the DNA chain during irradiation. Another part remains attached to the DNA chain backbone and for its characterization mild formic acid or enzymatic hydrolysis have been used. Starting from the chemical formulae of the altered base residues, the specific repair enzymes and their modes of action are discussed. Various glycosylases and endonucleases have been purified to homogeneity, and in some cases the gene which encodes the protein cloned. Using methods derived from Maxam and Gilbert sequencing procedures and DNA fragment 32P-labelled at one end, it has been shown that the alkali-labile sites in DNA induced by radiation are strongly dependent on the DNA base sequence. Enzymatic methods have been used to analyse the DNA base defects produced by gamma-irradiation of cells under in vivo conditions. Structures of modified bases were the same as those observed when DNA was irradiated in aqueous solution.     

Temperature dependence of the electronic absorption spectra of polynucleotides and their components]

The temperature dependence of UV-absorption spectra of solutions nucleic bases, nucleosides, nucleotides and metilated bases, uncapable of tautomerization has been studied. The nature of such dependence, its connection with hypochromic effect is discussed. It is shown that for some methods of investigating polynucleotides it is necessary to take into account the temperature changes spectra of monomers.     

On the mechanism of action of free purine bases on DNA synthesis in serum-starved L-cells treated with platelet extract

It has previously been found that there is a synergistic effect of free purine bases and low concentrations of dialyzed platelet extract on net synthesis of DNA in serum-starved fibroblast-like mouse L-cells. Experiments with a mutant line of L-cells that was deficient in hypoxanthine phosphoribosyl transferase (EC 2.4.2.8) indicated that purine bases had a stimulatory effect only if they were incorporated into cellular ribonucleotides. In the present paper it was shown that platelet extract induced the incorporation of hypoxanthine or adenine into both ATP and DNA. The induced net synthesis of DNA appears to take place in the nuclei and it requires that platelet extract is present in the medium only initially while free purine bases have to be present only later in the period of the experiment when DNA is being synthesized. The induction of both incorporation of free purine bases into DNA and ATP and of net DNA synthesis is dependent on heat-labile components in platelet extract. The extract cannot be substituted for by platelet derived growth factor.     

The crystal structure of the DNA-binding drug berenil: molecular modelling studies of berenil-DNA complexes.

The crystal structure of the DNA minor-groove DNA-binding drug berenil has been determined. Molecular-modelling techniques have been used to establish plausible binding modes of the structure to A-T sequences. These have shown that specific hydrogen bonds are possible between the amidine groups of the drug molecule and 02 atoms of thymine, although global energy minimisations tended to emphasise electrostatic interactions with phosphate groups rather than these hydrogen bonds with bases.     

Hydration of nucleic bases in dilute aqueous solutions. Apparent molar adiabatic and isothermal compressibilities, apparent molar volumes and their temperature slopes at 25 degrees C

The concentration increment of the ultrasound velocity has been measured with an accuracy of +/- 0.03 cm/s in dilute aqueous solutions of a variety of nucleic bases and their derivatives in the concentration range 0.5-1.5 mg/g H2O at temperatures of 15-35 degrees C. A new method for the precise measurement of ultrasound velocity in small volumes of liquids has been used. The values of the apparent molar adiabatic compressibilities plus the corresponding temperature slopes, apparent molar volumes with their temperature slopes, and apparent molar isothermal compressibilities at infinite dilution have been obtained. The regularities describing the signs of these values and their dependence on the chemical structure of the solute have been revealed. It is shown that these regularities can be described as a consequence of partial 'normalization' of some of the properties of water around the bases, namely, weaker structural contribution to compressibility, less negative temperature slope of compressibility and less negative structural contribution to the coefficient of thermal expansion of water.     

Incorporation of a complete set of deoxyadenosine and thymidine analogues suitable for the study of protein nucleic acid interactions into oligodeoxynucleotides. Application to the EcoRV restriction endonuclease and modification methylase

A complete set of dA and T analogues designed for the study of protein DNA interactions has been prepared. These modified bases have been designed by considering the groups on the dA and T bases that are accessible to proteins when these bases are incorporated into double-helical B-DNA [Seeman, N. C., Rosenberg, J. M., & Rich, A. (1976) Proc. Natl. Acad. Sci. U.S.A. 73, 804-808]. Each of the positions on the two bases, having the potential to interact with proteins, have been subject to nondisruptive, conservative change. Typically a particular group (e.g., the 6-NH2 of dA or the 5-CH3 of T) has been replaced with a hydrogen atom. Occasionally keto groups (the 2- and 4-keto oxygen atoms of T) have been replaced with sulfur. The base set has been incorporated into the self-complementary dodecamer d(GACGATATCGTC) at the central d(ATAT) sequence. Melting temperature determination shows that the modified bases do not destabilize the double helix. Additionally, circular dichroism spectroscopy shows that almost all the altered bases have very little effect on overall oligodeoxynucleotide conformation and that most of the modified oligomers have a B-DNA type structure. d(GATATC) is the recognition sequence for the EcoRV restriction modification system. Initial rate measurements (at a single oligodeoxynucleotide concentration of 20 microM) have been carried out with both the EcoRV restriction endonuclease and modification methylase. This has enabled a preliminary identification of the groups of the dA and T bases within the d(GATATC) sequence that make important contacts to both proteins.     

The Z-form of bacteriophage lambda DNA, modified in situ

Using the methods of chemical modification, restriction analysis and immune-electron microscopy it has been shown that the definite regions of the bacteriophage lambda DNA contain unpaired bases in situ. The distribution map of such sites along the genome has been constructed. The correlation of the in situ modification and the reaction with anti-Z-DNA antibodies is shown for the 44972 bp site of bacteriophage DNA. The possibility of the existence of Z-form DNA in situ is discussed.     

Reinterpretation of fluorescence of terbium ion-DNA complexes

Terbium (Tb3+) fluorescence was used to investigate local non-denaturation perturbations of double-helical DNA structure induced in this nucleic acid by various physical and chemical agents. It has been shown that the interaction of Tb3+ with DNA into which single-strand or double-strand breaks have been introduced by DNase I or by low doses of ionizing radiation does not influence the fluorescence of the lanthanide cation. On the other hand, interaction of terbium with DNA modified by the antitumour drug cis-diamminedichloroplatinum(II) at low levels of binding and by low doses of ultraviolet radiation (wavelength 254 nm) has been shown to result in substantial enhancement of the fluorescence of this cation. It has been proposed that the terbium fluorescent probe can also be exploited successfully for the purpose of analysing the guanine bases present in distorted double-stranded regions of DNA, in which only the vertical stacking of the base-pairs is altered.     

Template-directed primer extension catalyzed by the Tetrahymena ribozyme.

The Tetrahymena ribozyme has been shown to catalyze an RNA polymerase-like reaction in which an RNA primer is extended by the sequential addition of pN nucleotides derived from GpN dinucleotides, where N = A, C, or U. Here, we show that this reaction is influenced by the presence of a template; bases that can form Watson-Crick base pairs with a template add as much as 25-fold more efficiently than mismatched bases. A mutant enzyme with an altered guanosine binding site can catalyze template-directed primer extension with all four bases when supplied with dinucleotides of the form 2-aminopurine-pN.     

105 Activity of DNA ligase on substrates containing non-canonical structures

The expansion of trinucleotide repeat tracts (e.g. (CAG)_n tracts) has been shown to contribute to genomic instability and has been implicated in the pathogenesis of several neurodegenerative diseases, including Huntington’s Disease and Fragile X syndrome (Kovtun et al., 2008). While the molecular mechanism of this expansion is unknown, the ability of trinucleotide repeat sequences to form non-canonical secondary structures, such as hairpins, has been implicated as a multifaceted source of error (Gacy et al., 1995). Non-canonical DNA secondary structures have been shown to impact the action of enzymes in the base excision repair (BER) pathway, by which oxidatively damaged bases are removed. More specifically, there is evidence that trinucleotide repeat-containing DNA mistakenly enters long-patch BER, which can potentially lead to the incorporation of extra nucleobases by DNA polymerase (Jarem et al., 2011). The final enzyme in the BER pathway is DNA Ligase, which catalyses the formation of a phosphodiester bond to seal a nick site (Taylor et al., 2011). When extra nucleotides have been added during an erroneous long-patch BER process, the action of DNA ligase may expand the repeat tract by incorporating these additional bases into duplex DNA. In this study, DNA constructs containing (CAG)_n hairpins at various distances from a nick site are used to investigate the ability of DNA Ligase to ligate substrates containing non-canonical secondary structure back into duplex DNA.     

Characterisation of mRNAs encoding the precursor for human apolipoprotein CI.

cDNA clones encoding human apolipoprotein CI have been isolated from an adult liver cDNA library. Apo CI mRNA was shown to have two species of approximately 580 and 560 bases by RNA blot hybridisation. The intracellular precursor of apo CI was inferred from the cDNA sequence to be an 83 amino acid polypeptide consisting of the 57 residue mature protein and an additional 26 residue amino terminal signal peptide. The 5' untranslated regions of the messages are 63 and 40 bases as determined by primer extension and the 3' untranslated region 111 bases. A polyadenylation signal is situated 10 bases 3' of the poly(A) tall. The mRNA level of apo CI in human liver was significantly greater than that of apo All and apo E.     

An approach for searching insertions in bacterial genes leading to the phase shift of triplet periodicity

The concept of the phase shift of triplet periodicity (TP) was used for searching potential DNA insertions in genes from 17 bacterial genomes. A mathematical algorithm for detection of these insertions has been developed. This approach can detect potential insertions and deletions with lengths that are not multiples of three bases, especially insertions of relatively large DNA fragments (>100 bases). New similarity measure between triplet matrixes was employed to improve the sensitivity for detecting the TP phase shift. Sequences of 17,220 bacterial genes with each consisting of more than 1,200 bases were analyzed, and the presence of a TP phase shift has been shown in ～16% of analysed genes (2,809 genes), which is about 4 times more than that detected in our previous work. We propose that shifts of the TP phase may indicate the shifts of reading frame in genes after insertions of the DNA fragments with lengths that are not multiples of three bases. A relationship between the phase shifts of TP and the frame shifts in genes is discussed.     

Role of DNA damage in the inactivation of Saccharomyces cerevisiae yeasts by 313-nm ultraviolet light

The relative contribution of respiration photoinhibition and DNA damage in the lethal effect induced by 313 nm ultraviolet light (UV) has been investigated in some strains of the yeast Saccharomyces cerevisiae. It has been shown that cells inactivation is essentially due to photo-induced damage to DNA. By photoreactivation experiments it has been found that dimers of the pyrimidine bases are the main lethal photoproducts induced in the DNA by 313 nm ultraviolet light.