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.     

Study of DNA base-Li doped SiC nanotubes in aqueous solutions: a computer simulation study

Due to the importance of soluble nanotubes in biological systems, computational research on DNA base functionalized nanotubes is of interest. This study presents the quantitative results of Monte Carlo simulations of Li-doped silicon carbide nanotubes and its nucleic acid base complexes in water. Each species was first modeled by quantum mechanical calculations and then Monte Carlo simulations were applied to study their properties in aqueous solution. Solvation free energies were computed to indicate the solvation behavior of these compounds. The computations show that solvation free energies of the complexes of DNA bases with Li-doped SiC nanotubes are in the order: thymine > cytosine > adenine > guanine. The results of complexation free energies were also used to study the stability of related structures, which indicate that thymine-Li-doped SiC nanotubes produce the most stable compound among the four DNA base complexes.     

Hydration of nucleic acid bases studied using novel atom-atom potential functions

New simple atom-atom potential functions for simulating behavior of nucleic acids and their fragments in aqueous solutions are suggested. These functions contains 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 molecular interacts with two adjacent hydrophilic centres of the base. Studies of interaction of the bases with several water molecules showed that water-water interactions 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.     

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.     

Application of a universal solvation model to nucleic acid bases: comparison of semiempirical molecular orbital theory, ab initio Hartree-Fock theory, and density functional theory

The free energies of solvation of six nucleic acid bases (adenine, cytosine, hypoxanthine, guanine, thymine, and uracil) in water and chloroform are calculated using CM2 class IV charges and SM5.42R atomic surface tensions. Using any of three approximations to the electronic wave function (AM1, Hartree-Fock, or DFT), we obtain good agreement with experiment for five cases where the experimental results are known for the partition coefficients between the two solvents. Decomposition of the solvation effects into bulk electrostatic contributions and first-solvation-shell effects shows that the partitioning is dominated by the former, and this illustrates the importance of using accurate partial atomic charges for modeling these molecules in aqueous solution.     

Hydration of DNA bases: analysis of crystallographic data.

We present a systematic analysis of water structure around nucleic acid bases. We have examined 28 crystal structures of oligonucleotides, and have studied the patterns of water around the four bases, guanine, cytosine, adenine, and thymine. The geometries of water positions were calculated up to 4.00 A from base atoms. We have found conformation-dependent differences in both the geometry and extent of hydration of the bases.     

Aqueous hydration of nucleic acid constituents: Monte Carlo computer simulation studies

Monte Carlo computer simulations were performed on dilute aqueous solutions of thymine, cytosine, uracil, adenine, guanine, the dimethyl phosphate anion in the gauche-gauche conformation and a ribose and deoxyribose derivative. The aqueous hydration of each molecule was analysed in terms of quasi-component distribution functions based on the Proximity Criterion, and partitioned into hydrophobic, hydrophilic and ionic contributions. Color stereo views of selected hydration complexes are also presented. A preliminary discussion of the transferability of functional group coordination numbers is given. The results enable to comment on two current problems related to the hydration of nucleic acids: a) the theory of Dickerson and coworkers on the role of water in the relative stability of the A and B form of DNA and b) the idea of water bridges and filaments emerging from the computer simulation results on the hydration of DNA fragments by Clementi.     

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/).     

Aqueous Hydration of Nucleic Acid Constituents: Monte Carlo Computer Simulation Studies

Abstract

Monte Carlo computer simulations were performed on dilute aqueous solutions of thymine, cytosine, uracil, adenine, guanine, the dimethyl phosphate anion in the gauche-gauche conformation and a ribose and deoxyribose derivative. The aqueous hydration of each molecule was analysed in terms of quasi-component distribution functions based on the Proximity Criterion, and partitioned into hydrophobic, hydrophilic and ionic contributions. Color stereo views of selected hydration complexes are also presented. A preliminary discussion of the transferability of functional group coordination numbers is given. The results enable to comment on two current problems related to the hydration of nucleic acids: a) the theory of Dickerson and coworkers on the role of water in the relative stability of the A and B forms of DNA and b) the idea of water bridges and filaments emerging from the computer simulation results on the hydration of DNA fragments by Clementi.     

Compounds similar to acyclovir. VII. Attempts to construct an anti-herpetic agent (6-hydroxy-2-oxa-4-hexenyl derivatives of nucleic bases

Based on the available data on the acyclovir's mechanism of action we attempted to predict the antiherpetic activity of 6-hydroxy-2-oxahexen-4-yl derivatives of nucleic bases. In terms of this model 9-(6-hydroxy-2-oxahexen-4-yl) guanine might be active. 6-Hydroxy-2-oxahexen-4-yl derivatives of adenine, guanine, cytosine, thymine, uracil, 1,2,4-triazole-3 and 1,2,4-triazole-5-carboxamide have been synthesized and their activity against herpes virus I investigated. The guanine derivative proved to possess rather high activity (chemotherapeutical index 8).     

Expanding the Genetic Code: Unnatural Base Pairs in Biological Systems

Molecular Biology - The genetic code is considered to use five nucleic bases (adenine, guanine, cytosine, thymine and uracil), which form two pairs for encoding information in DNA and two pairs for...     

Compounds similar to acyclovir. II. Synthesis of acyclic analogs of 2'-deoxynucleosides

Acyclic analogues of nucleosides, viz.9- and 3-(1,6-dihydroxy-4-oxahex-3-yl)adenine, 9-(1,6-dihydroxy-4-oxahex-3-yl)guanine, 1-(1,6-dihydroxy-4-oxahex-3-yl)cytosine and thymine, with C3'-C4' bond of the furanose ring cleaved, have been prepared by condensation of trimethylsilyl derivatives of nucleic acid bases with 1,4,6-triacetoxy-3-oxahexane in the presence of Lewis acids followed by treatment with metanolic ammonia.     

Purine and pyrimidine contents of some desoxypentose nucleic acids

The distribution of purines and pyrimidines in desoxypentose nucleic acids prepared from a variety of animal and plant sources has been studied. 1. The nucleic acids were prepared from calf thymus, calf kidney, sheep spleen, horse spleen, chicken erythrocyte, turtle erythrocyte, trout sperm, shad testes, sea urchin sperm, wheat germ, and Pneumococcus Type III. 2. Separate hydrolyses were carried out for the determination of purines and pyrimidines. These procedures permitted nearly quantitative recovery of nucleic acid phosphorus in many of the preparations examined. 3. In the case of those preparations where a quantitative recovery was obtained it can be concluded that no bases other than adenine, guanine, thymine, and cytosine were present in appreciable amounts. 4. The distribution of purines and pyrimidines in all the nucleic acids studied renders the tetranucleotide hypothesis untenable. 5. The results of the analyses have indicated no great differences in the composition of these nucleic acids with respect to purines and pyrimidines.     

The nucleic acids of some insect viruses

Purine and pyrimidine bases have been estimated from the desoxyribonucleic acids of eleven insect viruses. Their proportions vary in the different species in a balanced way so that the molar ratios adenine:thymine and guanine:cytosine are constant and close to unity, whereas adenine + thymine:guanine + cytosine ranges from 0.71 to 1.87. This ratio is identical for some biologically dissimilar viruses, and no general parallelism is evident between DNA composition and biological relationship. Two different viruses from one host have distinct DNA's.     

Fluorometric determination of adenine and its derivatives by reaction with glyoxal hydrate trimer

A fluorometric determination of adenine has been devised. It was found that adenine gives strong fluorescence when reacting with glyoxal hydrate trimer in an acidic medium. The maximum wavelengths of excitation and emission spectra were 328 and 382 mμ, respectively. This reaction was successfully applied to determination of 0.2–1.0 or 2.0–10.0 mμmoles of adenine and its derivatives, in which the hydrogen atom at position 9 is substituted with ribose or substituted ribose. Other nucleic acid bases, guanine, uracil, thymine, and cytosine, did not interfere with this fluorometric determination of adenine.     

Polymethylene derivatives of nucleic bases with ω-functional groups: VII. Cytotoxicity in the series of <Emphasis Type="Italic">N</Emphasis>-(2-oxocyclohexyl)-ω-oxoalkyl-substituted purines and pyrimidines

New polymethylene derivatives of the nucleic bases with β-diketo function in the ω-position have been synthesized by alkylation of uracil, thymine, cytosine, hypoxanthine, adenine, and N ²-isobutyryl guanine with 2-(ω-chloroalkanoyl)cyclohexanones. The physicochemical characteristics of compounds synthesized and their effect on tumor K562 and HCT116 cell lines have been studied.     

How flexible are DNA constituents? The quantum-mechanical study

Relaxed force constants (RFCs) and vibrational root-mean-square deviations have been evaluated by the original calculation method for conformational parameters of the DNA structural units and their constituents: nucleic acid bases (uracile, thymine, cytosine, adenine and guanine) and their 'building blocks' (benzene, pyrimidine, imidazole and purine molecules), as well as the DNA backbone structural units: tetrahydrofuran, 1,2-dideoxyribose, methanol and orthophosphoric acid. It has been found that the RFCs for nomenclature torsions beta, gamma, epsilon; and sugar pseudorotation angle P in 1,2-dideoxyribose are sensible to the molecule conformation and their values are in the range of 1-25 kcal/(mole·rad2) obeying the inequality K(γ)> K(ε) > K(ρ) > K(β). The RFCs values for endocyclic torsions of nucleic acid bases six-member rings lie within 15-45?kcal/(mole·rad2) in pyrimidines and within 20-60?kcal/(mole·rad2) in purines. It is shown that the quantum zero-point motion effectively neglects the amino group non-planarity in cytosine, adenine and partially in guanine.     

Determination of trace amounts of 5-methylcytosine in DNA by reverse-phase high-performance liquid chromatography

A high-performance liquid chromatographic method to separate five major bases (cytosine, thymine, guanine, adenine, and uracil) and three minor methylated bases (5-methylcytosine, N6-methyladenine, and 7-methylguanine) has been developed using a volatile mobile phase under isocratic conditions. It is extended to quantitate 5-methylcytosine in trace amounts (1 in 20,000 cytosine residues). The suitability of the method has been verified by estimating 5-methylcytosine in DNAs of phi X174 and pBR322. The method has been applied to quantitate the extent of cytosine methylation in DNA of larval silk glands of Bombyx mori. Our results confirm that the pupal DNA of Drosophila melanogaster does not contain detectable amounts of 5-methylcytosine.     

Comparative study of permanganate oxidation reactions of nucleotide bases by spectroscopy

The permanganate oxidation of free nucleotide bases was successfully studied in aqueous solution of tetraethylammonium chloride using spectroscopic techniques. The reaction was highly selective toward thymine and uracil, less with cytosine, very little reaction on guanine, and no reaction on adenine.     

Quantitative determination of 5-methylcytosine in DNA by reverse-phase high-performance liquid chromatography

A method to separate the four major bases (cytosine, guanine, thymine and adenine) and the two minor modified bases (5-methylcytosine and 6N-methyladenine) in DNA has been developed. For optimal separation, several different buffer systems are available for isocratic elution. The 12 5-methylcytosine (5-mC) residues in the plasmid pBR322 can be determined with a deviation of less than 3% of the expected value and have been used for internal standardization. Formic acid hydrolysis of bases and probably of DNA does not lead to the deamination of cytosine or 5-mC and thus can be used routinely for DNA hydrolysis. Adenovirus or baculovirus DNA does not contain detectable amounts of 5-mC. The distribution of 5-mC in hamster cell DNA appears to be nonrandom in that different 5'-CpG-3'-containing restriction sites are methylated to different extents.