Effect of the sulfur atom on S2 and S4 positions of the uracil ring in different DNA:RNA hybrid microhelixes with three nucleotide base pairs

The effect of the sulphur atom on the uracil ring was analyzed in different DNA:RNA microhelixes with three nucleotide base-pairs, including uridine, 2-thiouridine, 4-thiouridine, 2,4-dithiouridine, cytidine, adenosine and guanosine. Distinct backbone and helical parameters were optimized at different density functional (DFT) levels. The Watson-Crick pair with 2-thiouridine appears weaker than with uridine, but its interaction with water molecules appears easier. Two types of microhelixes were found, depending on the H-bond of H2' hydroxyl atom: A-type appears with the ribose ring in ³ E-envelope C_3' -endo, and B-type in ² E-envelope C_2' -endo. B-type is less common but it is more stable and with higher dipole-moment. The sulphur atoms significantly increase the dipole-moment of the microhelix, as well as the rise and propeller twist parameters. Simulations with four Na atoms H-bonded to the phosphate groups, and further hydration with explicit water molecules were carried out. A re-definition of the numerical value calculation of several base-pair and base-stacking parameters is suggested.

Atomic charges for conformationally rich molecules obtained through a modified principal component regression

A modification of the principal component regression model is proposed for obtaining a fixed set of atomic charges (referred to as dipole-derived charges) optimized for reproducing the dipole moment of a conformationally rich molecule, i.e., a molecule with multiple local minima on the potential energy surface. The method does not require any adjustable parameters and requires the geometries of conformers, their dipole moments and atomic polar tensor (APT) charges as the only input data. The fixed atomic charges generated by the method not only reproduce the molecular dipole moment in all the conformers accurately, but are also numerically close to the APT charges, thereby ensuring accurate reproduction of the dipole moment variations caused by small geometrical distortions (e.g., by vibrations) of the conformers. The proposed method has been applied to canonical 2'-deoxyribonucleotides, the model DNA monomers, and the dipole-derived charges have been shown to outperform both the averaged APT and RESP charges in reproducing the dipole moments of large sets of conformers, thus demonstrating a potential usefulness of the dipole-derived charges as a 'reference point' for modeling polarization effects in conformationally rich molecules.

Can we treat ab initio atomic charges and bond orders as conformation-independent electronic structure descriptors?

It is shown that atomic charges and bond orders of 2′-deoxycytidine depend on the molecule conformation.

Nucleotides containing variously modified sugars: energetics, structure, and mechanical properties

The influence of various sugar residue modifications on intrinsic energetic, conformational, and mechanical properties of 2'-deoxyribonucleotide-5'-monophosphates (dNs) was comprehensively investigated using modern quantum chemical approaches. In total, fourteen sugar modifications, including double bonds and heteroatoms (S and N) inside the sugar ring, as well as fluorination in various positions, were analyzed. Among hundreds of possible conformational states of dNs, only two - AI and BI, corresponding to the most biologically significant forms of a double-helical DNA, were considered for each dN. It was established that the most of the studied modifications tend to strongly stabilize either AI or BI conformation of dNs both in the gas phase and in aqueous solution (modelled by implicit solvent models). Therefore, some of these modifications can be used as a tool for reducing structural polymorphism of nucleic acids in solution as well as for designing oligonucleotides with specific structural features. The evaluation of relaxed force constants (RFC) for glycosidic bonds suggests that the majority of the studied modifications of the sugar residue yield increased strengths of glycosidic bonds in dNs, and can therefore be used for designing modified nucleic acids with an increased resistance to abasic lesions. The most significant reinforcement of the glycosidic bond occurs in dNs containing the CF2 group instead of the O4' oxygen and the fluorine atom at the 2'-α-position. The calculation of the RFC and vibrational root-mean-square (VRMS) deviations for conformational degrees of freedom revealed a strong dependence between mechanical properties of dNs and their energetic characteristics. In particular, electronic energies of AI and BI conformers of dNs calculated in vacuo are closely connected with the values of relaxed force constants (RFC) for the δ angle: the higher RFC(δ) values correspond to more energetically favorable conformers.

THE EFFECT OF CHLORINATION OF NUCLEOTIDE BASES ON THE CONFORMATIONAL PROPERTIES OF THYMIDINE MONOPHOSPHATE

Recent studies on Escherichia coli bacteria cultivation, in which DNA thymine was replaced with 5-chlorouracil have refreshed the problem of understanding the changes to physical properties of DNA monomers resultant from chemical modifications. These studies have shown that the replacement did not affect the normal activities and division of the bacteria, but has significantly reduced its life span. In this paper a comparative analysis was carried out by the methods of computational experiment of a set of 687 possible conformers of natural monomeric DNA unit (2'-deoxyribonucleotide thymidine monophosphate) and 660 conformers of 5-chloro-2'-deoxyuridine monophosphate - a similar molecules in which the natural nitrogenous base thymine is substituted with 5-chlorouracil. Structures of stable conformers of the modified deoxyribonucleotide have been obtained and physical factors, which determine their variation from the conformers of the unmodified molecule have been analyzed. A comparative analysis of the elastic properties of conformers of investigated molecules and non-covalent interactions present in them was conducted. The results can be usedfor planning experiments on synthesis of artficial DNA suitable for incorporation into living organisms.

DFT studies of the formation of furanoid derivatives of ammonium chlorides

B3LYP/6-31+G** level computations were performed on the formation of ammonium salts during the reaction of (S)-1,4-anhydro-5-chloro-2,3,5-trideoxypentitol (1) (2S,5S)-2,5-anhydro-6-chloro-1,3,4,6-tetradeoxyhexitol (2) and methyl 5-chloro-2,3,5-trideoxy-β-D-pentofuranoside (3) with ammonia in order to describe the reaction pathway in detail. All the structures were fully optimized in the gas phase, in chloroform and water. In addition, the gas phase activation barrier heights were estimated at B3LYP/6-311++G**, MPWIK/6-31+G**, MPWIK/6-311++G** and MP2/6-311++G(2d,2p)//MPWIK/6-31+G** levels of theory. All the calculations in solvents were performed the using polarizable continuum model (PCM) and the B3LYP functional with the 6-31+G** basis set. A detailed description of all the stationary points is presented, and the conformational behavior of the five-membered ring is discussed in the gas phase and in the solvents. The conversion of the reactant complexes into ion pairs is accompanied by a strong energy decrease in the gas phase and in all the solvents. The overall process is strongly unfavorable in the gas phase, but takes place readily in high-polarity solvents.

The conformational behavior, geometry and energy parameters of Menshutkin-like reaction of O-isopropylidene-protected glycofuranoid mesylates in view of DFT calculations

The formation of pyridinium salts in the transformation of three O-isopropylidene-protected mesylates of furanoid sugar derivatives under pyridine action is considered at the B3LYP/6-31+G** computation level. All the structures were optimized in the gas phase, in chloroform and water. Activation barrier heights in the gas phase were also estimated at the B3LYP/6-311++G**, MPW1K/6-31+G** and MPW1K/6-311++G** levels. The conducted calculations, both in the gas phase (regardless of the computation level) and in solvents, revealed the barrier height increasing order as follows: 1>2>3 for the three reactions studied. The conformational behavior of the five-membered ring is discussed in the gas phase and in solvents. The fused dioxolane ring makes the furanoid ring less likely to undergo conformational changes. In the case of reaction 3, the furanoid ring shape does not change either in the gas phase or in solvents. All conformers are close to E0 or (0)E.