Amino−Imino Tautomerism in Derivatives of Cytosine: Effect on Hydrogen-Bonding and Stacking Properties

Understanding the structural basis of the binding specificity of c-di-AMP to M. smegmatis RecA using computational biology approach

Mycobacterium tuberculosis RecA (MtRecA), a protein involved in DNA repair, homologous recombination and SOS pathway, contributes to the development of multidrug resistance. ATP binding-site in RecA has been a drug target to disable RecA dependent DNA repair. For the first time, experiments have shown the existence and binding of c-di-AMP to a novel allosteric site in the C-terminal-Domain (CTD) of Mycobacterium smegmatis RecA (MsRecA), a close homolog of MtRecA. In addition, it was observed that the c-di-AMP was not binding to Escherichia coli RecA (EcRecA). This article analyses the possible interactions of the three RecA homologs with the various c-di-AMP conformations to gain insights into the structural basis of the natural preference of c-di-AMP to MsRecA and not to EcRecA, using the structural biology tools. The comparative analysis, based on amino acid composition, homology, motifs, residue types, docking, molecular dynamics simulations and binding free energy calculations, indeed, conclusively indicates strong binding of c-di-AMP to MsRecA. Having very similar results as MsRecA, it is highly plausible for c-di-AMP to strongly bind MtRecA as well. These insights from the in-silico studies adds a new therapeutic approach against TB through design and development of novel allosteric inhibitors for the first time against MtRecA.Communicated by Ramaswamy H. Sarma.

Characterization of tautomeric forms of anti-cancer drug gemcitabine and their interconversion upon mechano-chemical treatment, using ATR-FTIR spectroscopy and complementary methods

Gemcitabine is a widely used anti-cancer drug of pyrimidine structure, which can exist as a free base molecular form in crystals. Tautomers are structural isomers of molecules, which interconvert via proton transfer. Mechano-chemistry studies reactions of solids under mechanical impact. We investigated gemcitabine free base for the presence of specific molecular tautomers, using ATR-FTIR spectroscopic analysis, powder XRD, optical microscopy and HPLC. The amino-keto tautomer has the characteristic infrared (IR) peak of the amino group at 3390 cm-1. For the first time, the imino-keto tautomer of gemcitabine free base was detected. The imino-keto tautomer has the characteristic IR peak of the =N-H group, and its peak due to the CO group in pyrimidine ring is shifted vs. that of the amino-keto tautomer. This serves as the unique spectroscopic "fingerprints" of these tautomers. The ATR-FTIR spectroscopic analysis shows that gemcitabine free base can be enriched with the amino-keto or the imino-keto tautomer. Further, we studied the transformation of gemcitabine free base in crystals between its tautomers under conditions of liquid-assisted grinding (LAG). The imino-keto tautomer undergoes tautomerization to the amino-keto tautomer, while the amino-keto tautomer remains stable. No destruction of molecules of gemcitabine free base, when present as either tautomer, occurs during LAG as was verified by the HPLC-UV analysis. LAG is a new, straightforward, facile and fast method to interconvert tautomers in crystals, and ATR-FTIR spectroscopy is a method of choice to study tautomerization reactions of pharmaceuticals. The presented approach is promising for analysis of crystals of drugs containing one or more than one tautomer, and the knowledge-driven design of composite materials, which contain specific tautomeric molecular forms of pyrimidines, purines and other biologically active heterocyclic compounds.

Combinations of Tautomeric Forms and Neutral-Cationic Forms in the Cocrystals of Sulfamethazine with Carboxylic Acids

The cocrystals of sulfamethazine with different acids, namely, 2-mercaptophenylcarboxylic acid, 2,6-pyridinedicarboxylic acid, 4-(4-hydroxyphenylazo)phenylcarboxylic acid, 3-(4-hydroxyphenyl)propanoic acid, and 4-(phenyl)phenylcarboxylic acid, are studied here. Each has distinct notable supramolecular features. The pyrimidin-2-amine unit of the sulfamethazine provided unique examples of cocrystals in which amidine and imidine forms or neutral and protonated forms of sulfamethazine are observed in 2:2 ratios. Hence, this study provides avenues to explore cocrystals with tautomeric forms together in a cocrystal and also neutral and protonated cocrystal partners as apparent multicomponents in cocrystals. Among the cocrystals, three of them have the amidine form of the sulfamethazine in respective self-assembly. The cocrystal of 2-mercapto-phenylcarboxylic acid with sulfamethazine has the amidine form and it has the distinction of having S-H···π interactions. The cocrystal of sulfamethazine with 2,6-pyridinecarboxylic acid is a rare example of a 1:1 cocrystal of sulfmethazine with dicarboxylic acid. It has methanol molecules as a solvent of crystallization. Sulfamethazine forms a hydrated cocrystal with 4-(4-hydroxyphenylazo)-phenylcarboxylic acid that has conventional R 2 2(8) synthons of amidine hydrogen-bonding with carboxylic acid. The phenolic part of the acid component is anchored to the water molecule and provides a robust self-assembly. The hydrated cocrystal of sulfamethazine with 3-(4-hydroxyphenyl)propanoic acid (2:2 cocrystal) has two independent molecules of sulfamethazine, one in amidine form and the other in imidine form. It has two neutral carboxylic acids anchored through complementary hydrogen bonds and also has two water molecules of crystallization. The cocrystal of sulfamethazine with 4-(phenyl)phenylcarboxylic acid is also a 2:2 cocrystal. It is a di-hydrate, which has a neutral and protonated form of sulfamethazine. The neutral form is hydrogen-bonded to a neutral carboxylic acid, whereas the protonated form is charge-assisted hydrogen-bonded to the corresponding carboxylate anion.

Quantum mechanical studies on the existence of AC mismatches through prototopic tautomerization pathway in adenine and cytosine recognition

The gas phase studies on the electronic stabilities and thermochemical properties of metastable AC mismatches have been performed and the results are analyzed to explain the features of experimentally available AC mismatches. The hydrogen bonding patterns observed in these mismatches are relevant to the formation of stable AC mismatches. In these AC pairs, the H migration mechanism to generate other tautomeric forms is not observed, which shows the compatibility of H bonding capacity of the sites involved in H bonds. The presence of hydrogen bond type –N[Formula: see text]H–N– may contribute to better AC pairing, hence cisA1-1cisC2, transA2-cisC1 and transA1-1cisC2 are found to be more favorable pairs compared to other AC pairs.

Detailed theoretical investigation of excited-state intramolecular proton transfer mechanism of a new chromophore II

In the present work, TDDFT has been used to investigate the excited state intramolecular proton transfer (ESIPT) mechanism of a new chromophore II [Sensors and Actuators B: Chemical. 202 (2014) 1190]. The calculated absorption and fluorescence spectra agree well with experimental results. In addition, two types of II configurations are found in the first excited state (S1), which can be ascribed to the ESIPT reaction. Based on analysis of the calculated infrared (IR) spectra of O-H stretching vibration as well as the hydrogen bonding energies, the strengthening of the hydrogen bond in the S1 state has been confirmed. The frontier molecular orbitals (MOs), Hirshfeld charge distribution and the Natural bond orbital (NBO) have also been analyzed, which displays the tendency of the ESIPT process. Finally, potential energy curves of the S0 and S1 states were constructed, demonstrating that the ESIPT reaction can be facilitated based on the photo-excitation.

Understanding the amino ↔ imino tautomeric preference in (imidazole)imidazolidine-N-aryl(alkyl) systems: a case study of moxonidine drug and insights from the Cambridge structural database (CSD)

Moxonidine drug crystallizes exclusively as its imino tautomer. Conjugation plays a key role in the tautomer stability.

Nature and magnitude of aromatic stacking of nucleic acid bases

This review summarises recent advances in quantum chemical calculations of base-stacking forces in nucleic acids. We explain in detail the very complex relationship between the gas-phase base-stacking energies, as revealed by quantum chemical (QM) calculations, and the highly variable roles of these interactions in nucleic acids. This issue is rarely discussed in quantum chemical and physical chemistry literature. We further extensively discuss methods that are available for base-stacking studies, complexity of comparison of stacking calculations with gas phase experiments, balance of forces in stacked complexes of nucleic acid bases, and the relation between QM and force field descriptions. We also review all recent calculations on base-stacking systems, including details analysis of the B-DNA stacking. Specific attention is paid to the highest accuracy QM calculations, to the decomposition of the interactions, and development of dispersion-balanced DFT methods. Future prospects of computational studies of base stacking are discussed.

Infrared spectra of protonated uracil, thymine and cytosine

The gas-phase structures of protonated uracil, thymine, and cytosine are probed by using mid-infrared multiple-photon dissociation (IRMPD) spectroscopy performed at the Free Electron Laser facility of the Centre Laser Infrarouge d'Orsay (CLIO), France. Experimental infrared (IR) spectra are recorded for ions that were generated by electrospray ionization, isolated, and then irradiated in a quadrupole ion trap; the results are compared to the calculated infrared absorption spectra of the different low-lying isomers (computed at the B3LYP/6-31++G(d,p) level). For each protonated base, the global energy minimum corresponds to an enolic tautomer, whose infrared absorption spectrum matched very well with the experimental IRMPD spectrum, with the exception of a very weak IRMPD signal observed at about 1800 cm(-1) in the case of the three protonated bases. This signal is likely to be the signature of the second-energy-lying oxo tautomer. We thus conclude that within our experimental conditions, two tautomeric ions are formed which coexist in the quadrupole ion trap.

Comprehensive study of the effects of methylation on tautomeric equilibria of nucleic acid bases

Minor tautomers of nucleic acid bases can result by intramolecular proton transfer. These rare tautomers could be stabilized through the addition of methyl groups to DNA bases. A comprehensive theoretical study of tautomers of methylated derivatives of guanine, adenine, cytosine, thymine, and uracil was performed. Molecular geometries of all tautomers were obtained at the density functional theory and MP2 levels with the 6-31G(d,p) basis set, and single-point calculations were performed at the CCSD(T)/6-311G(d,p) level. Tautomers obtained by protonation at the preferred protonation site for methylated isolated bases were compared to their nonmethylated counterparts. The effects of methylation on the relative stabilities of nucleic acid base tautomers are also studied and discussed in this work. The results suggest that some sites on the bases may not be mutagenic and may even stabilize the canonical Watson-Crick form. The results also indicate that a number of methylation sites can stabilize the tautomers, suggesting possible mechanisms for mutagenic changes.

Adenine tautomers: relative stabilities, ionization energies, and mismatch with cytosine

In this study, we have investigated 12 tautomers of the DNA base adenine at the BP86/TZ2P and BP86/QZ4P levels of density functional theory. The vertical and adiabatic ionization energies of all tautomers were determined as the difference in energy between the radical cation and the corresponding neutral system. Furthermore, an evaluation is made for the eigenvalue spectra calculated with the SAOP functional, which is shown to lead to substantial improvements for orbital energies compared to BP86. We have also explored the correlations between the Kohn-Sham orbitals of the different tautomers at the BP86/QZ4P and SAOP/QZ4P levels. Finally, we discuss implications of the existence of the tautomeric forms of adenine for the DNA replication.

Dimers of and Tautomerism between 2-Pyrimidinethiol and 2(1H)-Pyrimidinethione: A Density Functional Theory (DFT) Study

Density functional theory (BLYP, B3LYP, B3P86, B3PW91) with the 6-31+G(d,p), 6-311+G(d,p), and cc-pVTZ basis sets has been used to calculate structural parameters, relative energies, and vibrational spectra of 2-pyrimidinethiol (1) and 2(1H)-pyrimidinethione (2) and their hydrogen-bonded homodimers (C(2) 3, C(2h) [4](double dagger), C(2h) 5), monohydrates, and dihydrates and a heterodimer (6). Several transition state structures proposed for the tautomerization process have also been examined. At the B3PW91/6-311+G(d,p)//B3PW91/6-31+G(d,p) level of theory 2-pyrimidinethiol (1) is predicted to be 3.41 kcal/mol more stable (E(rel)) than 2(1H)-pyrimidinethione (2) in the gas phase and 2 is predicted to be 6.47 kcal/mol more stable than 1 in aqueous medium. An unfavorable planar intramolecular strained four center transition state (TS1) for the tautomerization of 1 and 2 in the gas-phase lies 29.07 kcal/mol higher in energy than 2-pyrimidinethiol (1). The C(2) 2-pyrimidinethiol dimer (3) is 6.84 kcal/mol lower in energy than the C(2) homodimer transition state structure ([11](double dagger)) that connects dimers 3 and 4. Transition state [11](double dagger) provides a facile pathway for tautomerization between 1 and 2 in the gas phase (monomer-dimer promoted tautomerization). The hydrogen bonded 2-pyrimidinethiol- - -H(2)O and 2-pyrimidinethiol- - -2H(2)O structures are predicted to be 1.27 and 1.55 kcal/mol, respectively, higher in energy than 2(1H)-pyrimidinethione- - -H(2)O and 2(1H)-pyrimidinethione- - -2H(2)O. Water promoted tautomerization via cyclic transition states involving one water molecule (TS- - -H(2)O, [12](double dagger)) and two water molecules (TS- - -2H(2)O, [13](double dagger)) lie 11.42 and 11.44 kcal/mol, respectively, higher in energy than 2-pyrimidinethiol- - -H(2)O and 2-pyrimidinethiol- - -2H(2)O. Thus, the hydrated transition states [12](double dagger) and [13](double dagger) are involved in the tautomerism between 1 and 2 in aqueous medium.

Photoinduced Amino−Imino Tautomerism: An Infrared Study of 2-Amino-5-methylpyridine in a Low-Temperature Argon Matrix

The photoreaction of 2-amino-5-methylpyridine was investigated by matrix-isolation infrared spectroscopy and DFT calculation. Photoinduced reversible amino (N=C-NH(2))-imino (NH-C=NH) tautomerism was found between 2-amino-5-methylpyridine and 5-methyl-2(1H)-pyridinimine; the amino tautomer changes to the imino tautomer by UV irradiation (340>lambda>or= 300 nm) and the reverse change occurs by longer-wavelength light irradiation (420>lambda>or= 340 nm). The results of the CASSCF calculation revealed that the amino-imino tautomerism proceeds in vibrational relaxation process from electronic excited state to the ground state. The IR spectra of 2-amino-5-methylpyridine in the T(1) state and 5-methyl-2-pyridinamino radical were also obtained by UV irradiation (lambda>or= 300 nm).

Benzoderivatives of Nucleic Acid Bases as Modified DNA Building Blocks

The tautomeric properties of benzoderivatives of the canonical nucleic acid bases have been studied by using different computational approaches. Attention has been paid to the impact of the benzene group in altering the tautomeric preferences of the canonical bases both in the gas phase and in aqueous solution. To this end, relative solvation free energies of the tautomers determined from Self-Consistent Reaction Field continuum calculations and Monte Carlo-Free Energy Perturbation are combined with gas-phase tautomerization free energies determined from quantum mechanical calculations. The results provide a detailed picture of the tautomeric preferences of the benzoderivatives of nucleic acid bases. This information is used to examine the recognition properties of the preferred tautomers of the benzo-fused derivatives, paying particular attention to the ability to form Watson-Crick hydrogen-bonding and stacking interactions as well as to the hydrophobic nature of the modified bases. The implications of present results on the potential use of benzo-fused bases as potential building blocks in modified DNA duplexes are examined.

Singlet excited-state lifetimes of cytosine derivatives measured by femtosecond transient absorption

Lifetimes of the lowest excited singlet (S1) electronic states of various derivatives of the pyrimidine nucleobase cytosine (Cyt) were measured by the femtosecond transient absorption technique. The bases were excited in room-temperature aqueous solution at 265 nm using approximately 200 fs pump pulses from a titanium-sapphire laser system. The decay of excited-state absorption (ESA) at visible probe wavelengths was used to determine the S1 lifetimes of a variety of modified Cyt compounds at different pH values by global fitting. Identical lifetimes were observed for Cyt and cytidine (Cyd) within experimental uncertainty, but ESA by the ribonucleoside was considerably stronger, suggesting that the ribose group increases the oscillator strength of the S1 --> SN transition. The S1 lifetime of the important minor base 5-methylcytosine (m5Cyt) is 7.2 +/- 0.4 ps at pH 6.8. The same lifetime was measured for the ribonucleoside 5-methylcytidine, but sugar substitution again increased the strength of the ESA signal. Protonation of Cyd and m5Cyt at low pH led to a modest decrease in their S1 lifetimes. On the other hand, deprotonation of Cyt and m5Cyt significantly increased the lifetime of their respective S1 states. These trends support the intermediacy of the n,pi* state localized on the carbonyl oxygen in the nonradiative decay mechanism of Cyt. Longer S1 lifetimes were observed for 5-fluorocytosine and N4-acetylcytosine. Collectively, these results illustrate the great potential of femtosecond laser spectroscopy for investigating excited-state dynamics in DNA and DNA components.

2-Pyridinethiol/2-pyridinethione tautomeric equilibrium. A comparative experimental and computational study

The gas phase and solvent dependent preference of the tautomerization between 2-pyridinethiol (2SH) and 2-pyridinethione (2S) has been assessed using variable temperature Fourier transform infrared (FTIR) experiments, as well as ab initio and density functional theory computations. No spectroscopic evidence (nu(S)(-)(H) stretch) for 2SH was observed in toluene, C(6)D(6), heptane, or methylene chloride solutions. Although, C(s)() 2SH is 2.61 kcal/mol more stable than C(s)() 2S (CCSD(T)/cc-pVTZ//B3LYP/6-311+G(3df,2p)+ZPE), cyclohexane solvent-field relative energies (IPCM-MP2/6-311+G(3df,2p)) favor 2S by 1.96 kcal/mol. This is in accord with the FTIR observations and in quantitative agreement with the -2.6 kcal/mol solution (toluene or C(6)D(6)) calorimetric enthalpy for the 2S/2SH tautomerization favoring the thione. As the intramolecular transition state for the 2S, 2SH tautomerization (2TS) lies 25 (CBS-Q) to 30 kcal/mol (CCSD/cc-pVTZ) higher in energy than either tautomer, tautomerization probably occurs in the hydrogen bonded dimer. The B3LYP/6-311+G(3df,2p) optimized C(2) 2SH dimer is 10.23 kcal/mol + ZPE higher in energy than the C(2)(h)() 2S dimer and is only 2.95 kcal/mol + ZPE lower in energy than the C(2) 2TS dimer transition state. Dimerization equilibrium measurements (FTIR, C(6)D(6)) over the temperature range 22-63 degrees C agree: K(eq)(298) = 165 +/- 40 M(-)(1), DeltaH = -7.0 +/- 0.7 kcal/mol, and DeltaS = -13.4 +/- 3.0 cal/(mol deg). The difference between experimental and B3LYP/6-311+G(3df,2p) [-34.62 cal/(mol deg)] entropy changes is due to solvent effects. The B3LYP/6-311+G(3df,2p) nucleus independent chemical shifts (NICS) are -8.8 and -3.5 ppm 1 A above the 2SH and 2S ring centers, respectively, and the thiol is aromatic. Although the thione is not aromatic, it is stabilized by the thioamide resonance. In solvent, the large 2S dipole, 2-3 times greater than 2SH, favors the thione tautomer and, in conclusion, 2S is thermodynamically more stable than 2SH in solution.