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

New Insights into the Structure and Reactivity of Uracil Derivatives in Different Solvents-A Computational Study

Ab initio calculations were carried out to understand the reactivity and stability of some uracil derivatives, cytosine, 1-methyl cytosine, and cytidine in solvents, water, dimethyl sulfoxide (DMSO), n-octanol, and chloroform. Geometries were fully optimized at MP2 and B3LYP using the 6-31+G(d,p) basis set by applying the Solvation Model on Density (SMD) in solvent systems. The syn conformer of cytidine (cytidine II) is the most stable conformer in the gas phase, while the anticonformer (cytidine IV) is most stable in all of the solvents. Solvation free energy and polarizability values in different solvents decrease in the order water > DMSO > n-octanol > chloroform, while dipole moment, first-order hyperpolarizability, and HOMO-LUMO energy gap values follow the order of polar protic solvent (water and n-octanol) > polar aprotic solvent (DMSO) > nonpolar solvent (chloroform). The solvation free energy, dipole moment, polarizability, and first-order hyperpolarizability values also follow the order of cytosine > 1-methyl cytosine > cytidine. To illustrate that the molecular properties correlate well with the reactivity of the molecules, ab initio calculations were carried out for the reaction of uracil derivatives with Br₂ in the gas phase, water, DMSO, n-octanol, and chloroform. All ground and transition state geometries were fully optimized at B3LYP/6-31+G(d,p), and energies were also calculated at G3MP2 for cytosine and 1-methyl cytosine. For cytosine and 1-methyl cytosine, Gibbs energies of activation decrease with the polarity of the solvent that is chloroform > n-octanol > DMSO > water, while the Gibbs energies of activation for the reaction with cytidine decrease in the order of water > DMSO > n-octanol > chloroform. These results suggest that solvent polarity is very important for the stability and reactivity of uracil derivatives. Hydrogen bonding may also play an important role mainly for cytidine. Free energies of activation decrease with the size of the molecule, i.e., cytosine > 1-methyl cytosine > cytidine.

Electromechanical Signatures for DNA Sequencing through a Mechanosensitive Nanopore

Biological nanopores have been extensively used for DNA base detection since these pores are widely available and tunable through mutations. Distinguishing bases of nucleic acids by passing them through nanopores has so far primarily relied on electrical signals-specifically, ionic currents through the nanopores. However, the low signal-to-noise ratio makes detection of ionic currents difficult. In this study, we show that the initially closed mechanosensitive channel of large conductance (MscL) protein pore opens for single-stranded DNA (ssDNA) translocation under an applied electric field. As each nucleotide translocates through the pore, a unique mechanical signal is observed-specifically, the tension in the membrane containing the MscL pore is different for each nucleotide. In addition to the membrane tension, we found that the ionic current is also different for the four nucleotide types. The initially closed MscL adapts its opening for nucleotide translocation due to the flexibility of the pore. This unique operation of MscL provides single nucleotide resolution in both electrical and mechanical signals. Finally, we also show that the speed of DNA translocation is roughly 1 order of magnitude slower in MscL compared to Mycobacterium smegmatis porin A (MspA), suggesting MscL to be an attractive protein pore for DNA sequencing.

Predicting methylation status of human DNA sequences by pseudo-trinucleotide composition

DNA methylation plays a key role in the regulation of gene expression. The most common type of DNA modification consists of the methylation of cytosine in the CpG dinucleotide. The detections of DNA methylation have been determined mostly by experimental methods; however, these methods were time-consuming, expensive, and difficult to meet the requirements of modern large-scale sequencing technology. Accordingly, it is necessary to develop automatic and reliable prediction methods for DNA methylation. In this study, the pseudo-trinucleotide composition was proposed, and a novel method was developed by support vector machine (SVM) with the pseudo-trinucleotide composition as input parameter to represent DNA sequence for DNA methylation prediction. The model was evaluated on two datasets, including a dataset of Rollins (dataset_1) and a dataset collected healthy human records from the MethDB database (dataset_2). For dataset_1, the Matthews correlation coefficient (MCC) and accuracy (ACC) by jackknife validation were 0.8051 and 0.6098, respectively. For dataset_2, the MCC and ACC were 0.8500 and 0.7203, respectively. The good prediction results reveal that the pseudo-trinucleotide composition is an effective representation method for DNA sequence and plays a very important role in the prediction of DNA function.

Solvation free energies of glutamate and its metal complexes: a computer simulation study

Most biological ion channels demonstrate a high degree of selectivity for one type of ion more than others, and in many cases, how they control attaining this is still not clear. So we have studied on some metal ion compounds of glutamate. The Glutamate and its meal ion compounds (Ca(2+), Na(+), K(+) and Li(+)) were first modeled by ab initio calculations and then Monte Carlo simulation was used to calculate solvation free energies and also the complexes free energies for the related structures. The results indicated that Glutamate-Ca(2+) have more stability in water than other metal ion. Also, it was found out that the more movement in ions; less stability of the structure would result. This trend can be seen both in gas and liquid phase.

"False" thymine-1H-Enol guanine base pair. low misinsertion rate by DNA polymerase explained by computational chemistry consideration

Formation of correct TA and GC and "false" thymine-1H-enol guanine (TGenol) base pairs is here considered to control nucleotide insertion into DNA via low substrate concentration Michaelis-Menten controlled kinetics. Contributions of base pairing to formation of Gibbs free energies in water solution, DeltaDeltaG, are calculated for the correct and false base pairs with the semi-empiric MNDO/PM3 method for base pairing energies in vacuum and the BEM method for hydration effects. The results for DeltaDeltaG indicate equal insertion rates for correct base pairing and a 10(-3)-10(-4) error probability for false insertion controlled by the TGenol false pair.

Hydration and stability of nucleic acid bases and base pairs

Empirical, quantum chemical calculations and molecular dynamics simulations of the role of a solvent on tautomerism of nucleic acid bases and structure and properties of nucleic acid base pairs are summarized. Attention was paid to microhydrated (by one and two water molecules) complexes, for which structures found by scanning of empirical potential surfaces were recalculated at a correlated ab initio level. Additionally, isolated as well as mono- and dihydrated H-bonded, T-shaped and stacked structures of all possible nucleic acid base pairs were studied at the same theoretical levels. We demonstrate the strong influence of a solvent on the tautomeric equilibrium between the tautomers of bases and on the spatial arrangement of the bases in a base pair. The results provide clear evidence that the prevalence of either the stacked or hydrogen-bonded structures of the base pairs in the solvent is not determined only by its bulk properties, but rather by specific hydrophilic interactions of the base pair with a small number of solvent molecules.