Reaction of hypochlorous acid with imidazole: Formation of 2-chloro- and 2-oxoimidazoles

Reaction of hypochlorous acid (HOCl) with imidazole (Im) taken as a model for the 5-membered ring of guanine, leading to the products 2-chloro- and 2-oxo-imidazoles was investigated at the B3LYP/6-31+G* and B3LYP/AUG-cc-pVDZ levels of density functional theory. For all cases, single point energy calculations were performed at the MP2/AUG-cc-pVDZ level of theory using the geometries optimised at the B3LYP/AUG-cc-pVDZ level. Intrinsic reaction coordinate calculations were performed to ensure genuineness of all the calculated transition states. Effect of aqueous media was investigated by solvating all the species involved in the reactions using the polarizable continuum model. It is found that 2-chloroimidazole (2-ClIm) can be formed following three different reaction schemes while 2-oxoimidazole (2-oxoIm) can be formed following two different reaction schemes. The calculated barrier energies show that formation of 2-oxoIm would be less favored than that of 2-ClIm, which explains the experimental observations on relative yields of 8-chlorodeoxyguanosine and 8-oxodeoxyguanosine.

Electronic transitions of guanine tautomers, their stacked dimers, trimers and sodium complexes

Planar and nonplanar geometries of the keto-N9H and keto-N7H tautomers of the guanine base of DNA as well as the hydrogen bonded complexes of these species with three water molecules each were optimized using the density functional theory at the B3LYP/6-31G** level. Geometries of the isolated bases were also optimized using the ab initio approach at the MP2/6-31G** level. The isolated keto-N9H and keto-N7H tautomers as well as their hydrogen bonded complexes with three water molecules each were solvated in bulk water employing the polarized continuum model (PCM) of the self-consistent reaction field theory (SCRF). Stacked dimers and trimers of both the tautomers of guanine were generated by placing the planar forms of the species at a fixed distance of 3.5 A from the neighboring one and rotating one molecule with respect to the other by 110 degrees for the keto-N9H form and 90 degrees for the keto-N7H form which corresponded to total energy minima at the B3LYP/6-31G** level. Geometry optimization for the cation of the monomer of guanine was performed at the same level of theory, and its solvation in bulk water was treated using the PCM model of the SCRF theory. The geometries of complexes of the two tautomers of guanine with a Na+ ion each were optimized at the B3LYP/6-31G** level, and the Na+ ion is predicted to bind with the keto-N9H tautomer preferentially. While the complex of the keto-N7H form of guanine with three water molecules in gas phase is slightly more stable than the corresponding complex of the keto-N9H form of guanine, the reverse is true in bulk water. Stacking interactions enhance the relative stability of the keto-N9H tautomer over that of the keto-N7H tautomer, suggesting that in bulk solutions, the former would be dominant. Electronic spectra of the isolated tautomers of guanine, those of their complexes with three water molecules each, the (keto-n9h and keto-n7h) cation of guanine, the complexes of the tautomers with a Na+ ion each, the stacked dimers and trimers of the two tautomers were calculated using configuration interaction involving single electron excitations (CIS). The relative absorption intensities of the two tautomers of guanine near 275 and 248 nm in the monomer, dimer, and trimer are predicated to be in the opposite order. Thus the absorption intensity oscillation observed using a guanine aqueous solution can be explained in terms of oscillation of relative populations of the two tautomers of the molecule. The 248 nm absorption peak would be appreciably red-shifted on formation of the cation of guanine. Binding of the Na+ ion with the two tautomers of guanine reduces intensities of their transitions appreciably and also it causes large red-shifts in the same.

A geometry optimization and molecular electrostatic potential mapping study of structure-activity relationship for some anti-Alzheimer agents

Molecular geometries of some substituted (pyrroloamino)pyridines which possess anti-Alzheimer activity were optimized and potential-derived CHelpG point charges were computed using ab initio SCF molecular orbital approach employing the 3-21G basis set. AM1 molecular orbital calculations were performed using these optimized geometries and thus optimized Hybridization. Displacement Charges (HDC) combined with Löwdin charges continuously distributed in three dimension were obtained. Molecular electrostatic potential (MEP) maps of the molecules were obtained in two ways: (i) using the HDC-based model with the help of which MEP minima near the molecules were located, and (ii) using the CHelpG point charges, MEP values on the van der Waals surfaces of the molecules were computed. The MEP maps computed using both the methods have negative MEP regions near the pyridine nitrogen atom which appears to be the main binding site of the molecules with the appropriate receptor. Both electrostatic interaction and lipophilic association between these molecules and the receptor appear to contribute to biological activity.