Methylation of a cytidine deaminase inhibitor, 1-(beta-D-ribofuranosyl)-2-pyrimidone (i.e., zebularine (zeb)), which produces 1-(beta-D-ribofuranosyl)-5-methyl-2-pyrimidinone (d5), has been investigated using density functional theory models. The optimized structures of zeb and d5 and the valence orbitals primarily responsible for the methylation in d5 are presented using state-of-the-art interactive (on a computer or online) three-dimensional (3D) graphics in a portable document format (pdf) file, 3D-PDF (http://www.web3d.org/x3d/vrml/ ). The facility to embed 3D molecular structures into pdf documents has been developed jointly at Swinburne University of Technology and the National Computational Infrastructure, the Australian National University. The methyl fragment in the base moiety shows little effect on the sugar puckering but apparently affects anisotropic properties, such as condensed Fukui functions. Binding energy spectra, both valence space and core space, are noticeably affected; in particular, in the outer-valence space (e.g., IP < 20 eV). The methyl fragment delocalizes and diffuses into almost all valence space, but orbitals 8 (57a, IP = 12.57 eV), 18 (47a, IP = 14.70 eV), and 37 (28a, IP = 22.15 eV) are identified as fingerprint for the methyl fragment. In the inner shell, however, the impact of the methyl can be localized and identified by chemical shift. A small, global, red shift is found for the O-K, N-K and sugar C-K spectra, whereas the base C-K spectrum exhibits apparent methyl-related changes.

A computational approach, which involves the combination of the OPLS force field and molecular orbital MNDO, AM^, and P M~ methods, has been developed to describe the effects of a large, molecular mechanically simulated environment on the Hamiltonian of a quantum chemical system. To test the validity of the combined quantum mechanical/molecular mechanical (QM/MM) potential, a systematic study of the structures and energies of neutral and charged hydrogen-bonded complexes has been carried out, including comparisons with pure semiempirical calculations and available experimental and ab initio data. It is shown that, in many cases, the hybrid QM/MM potential behaves better than do related MNDO/M, AM^, and P M~ methods. As a case in point, the drawback of AM^ favoring bifurcated H-bonded structures over single ones is not presented in the combined A M~/ O P L S scheme. Possible ways of improvement of the combined OM/MM potential are discussed. 0

Two very different quantum mechanically based energy decomposition analyses (EDA) schemes are employed to study the dominant energy differences between the eclipsed and staggered ferrocene conformers. One is the extended transition state (ETS) based on the Amsterdam Density Functional (ADF) package and the other is natural EDA (NEDA) based in the General Atomic and Molecular Electronic Structure System (GAMESS) package. It reveals that in addition to the model (theory and basis set), the fragmentation channels more significantly affect the interaction energy terms (ΔE) between the conformers. It is discovered that such an interaction energy can be absorbed into the pre-partitioned fragment channels so that to affect the interaction energies in a particular conformer of Fc. To avoid this, the present study employs a complete fragment channel—the fragments of ferrocene are individual neutral atoms. It therefore discovers that the major difference between the ferrocene conformers is due to the quantum mechanical Pauli repulsive energy and orbital attractive energy, leading to the eclipsed ferrocene the energy preferred structure. The NEDA scheme further indicates that the sum of attractive (negative) polarization (POL) and charge transfer (CL) energies prefers the eclipsed ferrocene. The repulsive (positive) deformation (DEF) energy, which is dominated by the cyclopentadienyle (Cp) rings, prefers the staggered ferrocene. Again, the cancellation results in a small energy residue in favour of the eclipsed ferrocene, in agreement with the ETS scheme. Further Natural Bond Orbital (NBO) analysis indicates that all NBO energies, total Lewis (no Fe) and lone pair (LP) deletion all prefer the eclipsed Fc conformer. The most significant energy preferring the eclipsed ferrocene without cancellation is the interactions between the donor lone pairs (LP) of the Fe atom and the acceptor antibond (BD*) NBOs of all C–C and C–H bonds in the ligand, LP(Fe)-BD*(C–C & C–H), which strongly stabilizes the eclipsed (D5h) conformation by −457.6 kcal·mol−1.