A non-empirical molecular orbital study on the relative stabilities of adenine and guanine tautomers

A density functional theory study of the molecular structure, reactivity, and spectroscopic properties of 2-(2-mercaptophenyl)-1-azaazulene tautomers and rotamers

Five stable tautomer and rotamers of the 2-(2-Mercaptophenyl)-1-azaazulene (thiol, thione, R1, R2, and R3) molecules were studied using density functional theory (DFT). The geometries of the studied tautomer and rotamers were fully optimized at the B3LYP/6-31G(d,p) level. Thermodynamic calculations were performed at M06-2X/6-311G++(2d,2p) and ωB97XD/6-311G++(2d,2p) in the gas phase and ethanol solution conditions modeled by the solvation model based on density (SMD). The kinetic constant of tautomer and rotamers conversion was calculated in the temperature range 270-320 K using variational transition state theory (VTST) accompanied by one-dimensional wigner tunneling correction. Energy refinement at CCSD(T)/6-311++G(2d,2p) in the gas phase has been calculated. All the studied DFT methods qualitatively give similar tautomer stability orders in the gas phase. The ethanol solvent causes some reordering of the relative stability of 2-(2-Mercaptophenyl)-1-azaazulene conformers. The transition states for the 2-(2-Mercaptophenyl)-1-azaazulene tautomerization and rotamerization processes were also determined. The reactivity, electric dipole moment, and spectroscopic properties of the studied tautomer and rotamers were computed. The hyper-Rayleigh scattering (βHRS), and depolarization ratio (DR) exhibited promising optical properties when nonlinear optical properties were calculated.

Investigation of tautomerism of 1,3,5-triazine derivative, stability, and acidity of its tautomers from density functional theory

Recent studies have identified N2,N4-bis(4-fluorophenethyl)-N6-(3-(dimethylamino)propyl)-1,3,5-triazine-2,4,6-triamine (1TZ(7,8,9)) as a potent, pure antagonist that inhibits thermosensory transient receptor potential vanilloid 1 channel (TRPV1) channel activity. This study provides theoretical data on the stability and acidity of the tautomers of this molecule. We show that this triazine can exist as three predominant tautomers (2TZ(5,7,8), 4TZ(3,7,9), 7TZ(1,8,9)). In the aqueous phase, equilibrium constants calculations show that only the tautomeric equilibria between 1TZ(7,8,9) and the three most stable triazines can be present which suggests that these three tautomeric equilibria would be the basis of 1TZ(7,8,9)'s biological activity.

Automated exploration of the low-energy chemical space with fast quantum chemical methods

We propose and discuss an efficient scheme for the in silico sampling for parts of the molecular chemical space by semiempirical tight-binding methods combined with a meta-dynamics driven search algorithm. The focus of this work is set on the generation of proper thermodynamic ensembles at a quantum chemical level for conformers, but similar procedures for protonation states, tautomerism and non-covalent complex geometries are also discussed. The conformational ensembles consisting of all significantly populated minimum energy structures normally form the basis of further, mostly DFT computational work, such as the calculation of spectra or macroscopic properties. By using basic quantum chemical methods, electronic effects or possible bond breaking/formation are accounted for and a very reasonable initial energetic ranking of the candidate structures is obtained. Due to the huge computational speedup gained by the fast low-cost quantum chemical methods, overall short computation times even for systems with hundreds of atoms (typically drug-sized molecules) are achieved. Furthermore, specialized applications, such as sampling with implicit solvation models or constrained conformational sampling for transition-states, metal-, surface-, or noncovalently bound complexes are discussed, opening many possible applications in modern computational chemistry and drug discovery. The procedures have been implemented in a freely available computer code called CREST, that makes use of the fast and reliable GFNn-xTB methods.

Tautomer preference in PDB complexes and its impact on structure-based drug discovery

Tautomer enrichment is a key step of ligand preparation prior to virtual screening. In this paper, we have investigated how tautomer preference in various media (water, gas phase, and crystal) compares to tautomer preference at the active site of the protein by analyzing the different possible H-bonding contacts for a set of 13 tautomeric structures. In addition, we have explored the impact of four different protocols for the enumeration of tautomers in virtual screening by using Flap, Glide, and Gold as docking tools on seven targets of the DUD data set. Excluding targets in which the binding does not involve tautomeric atoms (HSP90, p38, and VEGFR2), we found that the average receiver operating characteristic curve enrichment at 10% was 0.25 (Gold), 0.24 (Glide), and 0.50 (Flap) by considering only tautomers predicted to be unstable in water versus 0.41 (Gold), 0.56 (Glide), 0.51 (Flap) by limiting the enumeration process only to the predicted most stable tautomer. The inclusion of all tautomers (stable and unstable) yielded slightly poorer results than considering only the most stable form in water.

Ultraviolet photolysis of adenine: Dissociation via the π1σ* state

High resolution total kinetic energy release (TKER) spectra of the H atom fragments resulting from photodissociation of jet-cooled adenine molecules at 17 wavelengths in the range 280>lambda(phot)>214 nm are reported. TKER spectra obtained at lambda(phot)>233 nm display broad, isotropic profiles that peak at low TKER ( approximately 1800 cm(-1)) and are largely insensitive to the choice of excitation wavelength. The bulk of these products is attributed to unintended multiphoton dissociation processes. TKER spectra recorded at lambda(phot)<or=233 nm display additional fast structure, which is attributed to N(9)-H bond fission on the 1pisigma* potential energy surface (PES). Analysis of the kinetic energies and recoil anisotropies of the H atoms responsible for the fast structure suggests excitation to two 1pisigma* excited states (the (1)L(a) and (1)B(b) states) at lambda(phot) approximately 230 nm, both of which dissociate to yield H atoms together with ground state adeninyl fragments by radiationless transfer through conical intersections with the 1pisigma* PES. Parallels with the photochemistry exhibited by other, smaller heteroaromatics (pyrrole, imidazole, phenol, etc.) are highlighted, as are inconsistencies between the present conclusions and those reached in two other recent studies of excited state adenine molecules.