Determinations of atomic partial charges for nucleic acid constituents from x-ray diffraction data. I. 2'-Deoxycytidine-5'-monophosphate

Metal Ion Modeling Using Classical Mechanics

Metal ions play significant roles in numerous fields including chemistry, geochemistry, biochemistry, and materials science. With computational tools increasingly becoming important in chemical research, methods have emerged to effectively face the challenge of modeling metal ions in the gas, aqueous, and solid phases. Herein, we review both quantum and classical modeling strategies for metal ion-containing systems that have been developed over the past few decades. This Review focuses on classical metal ion modeling based on unpolarized models (including the nonbonded, bonded, cationic dummy atom, and combined models), polarizable models (e.g., the fluctuating charge, Drude oscillator, and the induced dipole models), the angular overlap model, and valence bond-based models. Quantum mechanical studies of metal ion-containing systems at the semiempirical, ab initio, and density functional levels of theory are reviewed as well with a particular focus on how these methods inform classical modeling efforts. Finally, conclusions and future prospects and directions are offered that will further enhance the classical modeling of metal ion-containing systems.

Four Ways to Oligonucleotides Without Phosphoimidazolides

Emergence of the very first RNA or RNA-like oligomers from simple nucleotide precursors is one of the most intriguing questions of the origin of life research. In the current paper, we analyse the mechanism of four non-enzymatic template-free scenarios suggested for the oligomerization of chemically non-modified cyclic and acyclic nucleotides in the literature. We show that amines may have a twofold role in these syntheses: due to their high affinity to bind protons they may activate the phosphorus of the phosphate group via proton transfer reactions, or indirectly they may serve as charge compensating species and influence the self-assembling of nucleotides to supramolecular architectures compatible with the oligomerization reactions. Effect of cations and pH on the reactions is also discussed.

Solvated ensemble averaging in the calculation of partial atomic charges

In the calculation of partial atomic charges, for use in molecular mechanics or dynamics simulations, it is common practice to select only a single conformation for the molecule of interest. For molecules that contain rotatable bonds, it is preferable to compute the charges from several relevant conformations. We present here results from a charge derivation protocol that determines the partial charges by averaging charges computed for conformations selected from explicitly solvated MD simulations, performed under periodic boundary conditions. This approach leads to partial charges that are weighted by a realistic population of conformations and that are suitable for condensed phase simulations. This protocol can, in principle, be applied to any class of molecule and to nonaqueous solvation. Carbohydrates contain numerous hydroxyl groups that exist in an ensemble of orientations in solution, and in this report we apply ensemble averaging to a series of methyl glycosides. We report the extent to which ensemble averaging leads to charge convergence among the various monosaccharides and among the constituent atoms within a given monosaccharide. Due to the large number of conformations (200) in our ensembles, we are able to compute statistically relevant standard deviations for the partial charges. An analysis of the standard deviations allows us to assess the extent to which equivalent atom types may, nevertheless, require unique partial charges. The configurations of the hydroxyl groups exert considerable influence on internal energies, and the limits of ensemble averaged charges are discussed in terms of these properties.

Modelling the effect of structure and base sequence on DNA molecular electronics

DNA is a material that has the potential to be used in nanoelectronic devices as an active component. However, the electronic properties of DNA responsible for its conducting behaviour remain controversial. Here we use a self-consistent quantum molecular dynamics method to study the effect of DNA structure and base sequence on the energy involved when electrons are added or removed from isolated molecules and the transfer of the injected charge along the molecular axis when an electric field is applied. Our results show that the addition or removal of an electron from DNA molecules is most exothermic for poly(dC)-poly(dG) in its B-form and poly(dA)-poly(dT) in its A-form, and least exothermic in its Z-form. Additionally, when an electric field is applied to a charged DNA molecule along its axis, there is electron transfer through the molecule, regardless of the number and sign of the injected charge, the molecular structure and the base sequence. Results from these simulations provide useful information that is hard to obtain from experiments and needs to be considered for further modelling aiming to improve charge transport efficiency in nanoelectronic devices based on DNA.

Conformation and proton configuration of pyrimidine deoxynucleoside oxidation damage products in water

Emerging data strongly suggest that the oxidation of DNA bases can contribute to genomic instability. Structural changes to DNA, induced by base oxidation, may reduce the fidelity of DNA replication and interfere with sequence-specific DNA-protein interactions. We have examined the structures of a series of pyrimidine deoxynucleoside oxidation damage products in aqueous solution. The modified nucleosides studied include the deoxynucleoside derivatives of 5-hydroxyuracil, 5-hydroxycytosine, 5-(hydroxymethyl)uracil, 5-(hydroxymethyl)cytosine, 5-formyluracil, and 5-formylcytosine. The influence of base oxidation on ionization constants, sugar conformation, and tautomeric configuration has been determined on the basis of UV, proton, and nitrogen NMR spectra of the (15)N-enriched derivatives. The potential biological consequences of the structural perturbations resulting from base oxidation are discussed.

Atomic charges for DNA constituents derived from single-crystal X-ray diffraction data

We have derived a complete set of atomic charges for DNA from very high resolution, low temperature, single-crystal X-ray diffraction data, collected for a variety of nucleosides and nucleotides: cytidine; deoxycytidine 5'-monophosphate; deoxythymidine; guanosine 5'-monophosphate; deoxyadenosine; adenosine. This set of charges represents the first experimentally based parameterization of an important term in the energy function used in most modeling of DNA. The resulting charges are in good agreement with chemical intuition and experimental observations. They also agree qualitatively with the theoretically derived values now commonly used, but numerous and significant quantitative differences are observed. Possible reasons for the quantitative disagreement are discussed. An averaged set of charges (derived from the experimental results), which can be used in DNA modeling calculations, is presented.

Theoretical analysis of deoxycytidine-5'-phosphate protonation

The electron density distribution in deoxycytidine-5'-monophosphate (5'-dCMP) molecule and dianion has been studied by the method of CNDO/2. The comparison between the results of calculation for the neutral molecule and the data obtained by Pearlman and Kim shows that there is a linear correlation between the atomic charges calculated using quantum chemistry and those derived from X-ray results. However, partial charges for the deoxyribose fragment are correlated in a nonlinear manner. The influence of the protons added to the cytosine and phosphate residues on the atomic charges and bond orders of deoxy-cytidine-5'-monophosphate has been analyzed here. The conclusion has been drawn that the semiempirical quantum-chemical CNDO/2 technique is applicable to the mononucleotide studies.