Tautomeric Equilibria in Relation to Pi-Electron Delocalization

Intramolecular and Water Mediated Tautomerism of Solvated Glycine

Understanding tautomerism and characterizing solvent effects on the dynamic processes pose significant challenges. Using enhanced-sampling molecular dynamics based on state-of-the-art deep learning potentials, we investigated the tautomeric equilibria of glycine in water. We observed that the tautomerism between neutral and zwitterionic glycine can occur through both intramolecular and intermolecular proton transfers. The latter proceeds involving a contact anionic-glycine-hydronium ion pair or separate cationic-glycine-hydroxide ion pair. These pathways with comparable barriers contribute almost equally to the reaction flux.

A Torquoselective Thermal 6π-Electrocyclization Approach to 1,4-Cyclohexadienes via Solvent-Aided Proton Transfer: Experimental and Theoretical Studies

The thermal 6π-electrocyclization of hexatriene typically delivers 1,3-cyclohexadiene (1,3-CHD). However, there is only limited success in directly synthesizing 1,4-cyclohexadiene (1,4-CHD) using such an approach, probably due to the difficulty in realizing thermally-forbidden 1,3-hydride shift after electrocyclic ring closure. The present study shows that by heating (2E,4E,6E)-hexatrienes bearing ester or ketone substituents at the C1-position in a mixture of toluene/MeOH or EtOH (2 : 1) solvents at 90-100 °C, 1,4-CHDs can be selectively synthesized. This is achieved through a torquoselective disrotatory 6π-electrocyclic ring closure followed by a proton-transfer process. The success of this method depends on the polar protic solvent-assisted intramolecular proton transfer from 1,3-CHD to 1,4-CHD, which has been confirmed by deuterium-labeling experiments. There are no reports to date for such a solvent-assisted isomerization. Density functional theory (DFT) studies have suggested that forming 1,3-CHD and subsequent isomerization is a thermodynamically feasible process, regardless of the functional groups involved. Two possible successive polar solvent-assisted proton-transfer pathways have been identified for isomerization.

Azaindolizine proton cranes attached to 7-hydroxyquinoline and 3-hydroxypyridine: a comparative theoretical study

Theoretical design of several proton cranes, based on 7-hydroxyquinoline and 3-hydroxypyridine as proton-transfer frames, has been attempted using ground and excited-state density functional theory (DFT) calculations in various environments. Imidazo[1,2-a]pyridine, pyrazolo[1,5-a]pyridine and benzimidazole were considered as proton crane units. The proton crane action requires the existence of a single enol-like form in the ground state, which under excitation goes to the end keto-like one through a series of consecutive excited-state intramolecular proton transfers (ESIPT) and twisting steps with the participation of a crane unit, resulting in a long-range intramolecular proton transfer. The results suggest that 3-hydroxypyridine is not suitable for a proton-transfer frame and 8-(imidazo[1,2-a]pyridin-2-yl)quinolin-7-ol and 8-(pyrazolo[1,5-a]pyridin-2-yl)quinolin-7-ol behave as non-conjugated proton cranes, instead of tautomeric re-arrangement in the latter, which was thought to be possible.

Theoretical investigation of tautomerism of 2- and 4-pyridones: origin, substituent and solvent effects

Computational investigation at the BHandHLYP/6-311+G(d,p) level of theory of the gas-phase tautomerism of 2- and 4-pyridones confirmed the slight prevalence of lactim in the case of the former, but its dominance in the case of the latter, as shown previously. Examination of aromaticity by using HOMA, EDDB, NBOdel, NICS and AICD led to the conclusion that tautomerization of 4-pyridone results in greater aromaticity gain. It is also driven by the Pauli repulsion relief, which was revealed by the tautomerization energy decomposition analysis. By contrast, in the case of 2-pyridone, lactim is favoured by orbital and electrostatic interactions and disfavoured by the Pauli repulsion. Aromaticity gain in this case is smaller. The position of the tautomeric equilibrium can be modulated by substituent inductive effects (Cl and F), inductive and resonance effects (NH2 and NO2), hydrogen bonding (NO2), and medium polarity, the increase of which increases lactam population.

On Prototropy and Bond Length Alternation in Neutral and Ionized Pyrimidine Bases and Their Model Azines in Vacuo

In this review, the complete tautomeric equilibria are derived for disubstituted pyrimidine nucleic acid bases starting from phenol, aniline, and their model compounds-monosubstituted aromatic azines. The differences in tautomeric preferences for isolated (gaseous) neutral pyrimidine bases and their model compounds are discussed in light of different functional groups, their positions within the six-membered ring, electronic effects, and intramolecular interactions. For the discussion of tautomeric preferences and for the analysis of internal effects, recent quantum-chemical results are taken into account and compared to some experimental ones. For each possible tautomer-rotamer of the title compounds, the bond length alternation, measured by means of the harmonic oscillator model of electron delocalization (HOMED) index, is examined. Significant HOMED similarities exist for mono- and disubstituted derivatives. The lack of parallelism between the geometric (HOMED) and energetic (ΔG) parameters for all possible isomers clearly shows that aromaticity is not the main factor that dictates tautomeric preferences for pyrimidine bases, particularly for uracil and thymine. The effects of one-electron loss (positive ionization) and one-electron gain (negative ionization) on prototropy and bond length alternation are also reviewed for pyrimidine bases and their models.

Elucidating the relationship between red fluorescence and structural dynamics of carbon dots dispersed in different solvents

The mechanism of intrinsic fluorescence of carbon dots (CDs), the latest nanomaterial from the carbon family, was supposedly deciphered through multiple theories. However, the much sought-after persistent red emission of CDs as a foreseeable consequence of experiments remains elusive prompting the question of whether tuning of the red emission of CDs is a predictable outcome or a serendipitous coincidence. Herein, we tried to decode the same by exploring Alizarin Red S (ARS)-based red emitting CDs in different solvents with wisely chosen analytical tools. The findings are aptly supported by molecular dynamics studies through an experimental intuition-driven model-building approach. Parallel interception of the CDs with powder X-ray diffraction (pXRD) and photophysical spectroscopic studies revealed an important relationship between the solvent and CDs. Tautomerism, a well-known phenomenon with chemical entities, was found to be operative for CDs that greatly influence the Stokes shift and ultimately the fluorescence outcome. Most importantly, pXRD studies established the turbostratism of the CDs where the well-ordered graphitic structure of CDs gets disrupted with solvent molecules. The extent of such disruption is a function of solvent and CD composition that plays a formidable role in obtaining red fluorescence. Thus, for the first time, we demonstrate that the red emission of CDs is related to its structural integrity and if taken care of could be sustained, a tremendously desirable outcome for relevant applications.

DFT Meets Wave-Function Composite Methods for Characterizing Cytosine Tautomers in the Gas Phase

A general strategy for the accurate computation of structural and spectroscopic properties of biomolecule building blocks in the gas phase has been further improved and validated with a special reference to tautomeric equilibria. The main improvements concern the use of the cc-pVTZ-F12 basis set in both DFT and CCSD(T)-F12 computations, the inclusion of core-valence correlation in geometry optimizations by double hybrid functionals, and the use of the cc-pVQZ-F12 basis set for complete basis set extrapolation at the MP2-F12 level. The resulting model chemistry is applied to the challenging problem of cytosine tautomers in the gas phase. The results are in remarkable agreement with experiment concerning both rotational and vibrational spectroscopic parameters and permit their unbiased interpretation in terms of structural and thermochemical features. Together with the intrinsic interest of the studied molecule, the accuracy of the results obtained at reasonable cost without any empirical parameter suggests that the proposed composite method can be profitably employed for accurate investigations of other molecular bricks of life.

On Integral INICS Aromaticity of Pyridodiazepine Constitutional Isomers and Tautomers

The structure, energetics, and aromaticity of c.a. 100 constitutional isomers and tautomers of pyrido[m,n]diazepines (m = 1, 2; n = 2, 3, 4, 5; m ≠ n) were studied at the B3LYP/cc-pVTZ level. The pyrido[1,3]diazepines appear the most, while pyrido[2,4]diazepines are the least stable (ca. 26 kcal/mol). In the pyrido[1,n]diazepine group (n = 2-5), the [1,5] isomers are higher in energy by ca. 4.5 kcal/mol and the [1,4] ones by ca. 7 kcal/mol, and the pyrido[1,2]diazepines are the least stable (ca. 20 kcal/mol). All the most stable pyrido[1,n]diazepines have N-atoms near the ring's junction bond but on opposite sites. The most stable [2,n]-forms are also those with the pyridine ring N6-atom near the junction bond. Surprisingly, for the [1,2]-, [1,3]-, and [1,4]-isomer condensation types of pyridine and diazepine rings, the same N9 > N7 > N6 > N8 stability pattern obeys. The stability remains similar in a water medium simulated with the Polarizable Continuum Model of the solvent and is conserved when calculated using the CAM-B3LYP or BHandHlyp functionals. The ring's aromaticity in the pyridine[m,n]diazepines was established based on the integral INICS index resulting from the NICSzz-scan curves' integration. The integral INICS index is physically justified through its relation to the ringcurrent as demonstrated by Berger, R.J.F., et al. Phys. Chem. Chem. Phys. 2022, 24, 624. The six-membered pyrido rings have negative INICSZZ indices and can be aromatic only if they are not protonated at the N-atom. All protonated pyrido and seven-membered rings exhibit meaningful positive INICSZZ values and can be assigned as antiaromatic. However, some non-protonated pyrido rings also have substantial positive INICSZZ indices and are antiaromatic. A weak linear correlation (R2 = 0.72) between the INICSZZ values of the pyridine I(6) and diazepine I(7) rings exists and is a consequence of the communication between the π-electron systems of the two rings. The juxtaposition of the INICS descriptor of the six- and seven-membered rings and diverse electron density parameters at the Ring Critical Points (RCP) revealed good correlations only with the Electrostatic Potentials from the electrons and nuclei (ESPe and ESPn). The relationships with other RCP parameters like electron density and its Laplacian, total energy, and the Hamiltonian form of kinetic energy density were split into two parts: one nearly constant for the six-membered rings and one linearly correlating for the seven-membered rings. Thus, most of the electron density parameters at the RCP of the six-membered rings of pyridodiazepines practically do not change with the diazepine type and the labile proton position. In contrast, those of the seven-membered rings display aromaticity changes in the antiaromatic diazepine with its ring structural modifications.

Prototropy, Intramolecular Interactions, Electron Delocalization, and Physicochemical Properties of 1,8-dihydroxy-9-anthrone-DFT-D3 Study of Substituent Effects

1,8-dihydroxy-9-anthrone are tricyclic compounds with a ketone group in the middle ring and two hydroxyl groups substituted in the side-aromatic rings what results in formation of two intramolecular hydrogen bonds in which the oxygen atom from the ketone group is the proton acceptor. 1,8-dihydroxy-9-anthrones in which intramolecular proton transfer between C10 and CO in the middle ring occurs, can exist in a tautomeric keto-enol equilibrium. For anthralin, the most important representative of this group, this equilibrium has been studied previously, but it has not been studied for its derivatives. Substituents in the middle ring change the geometry of 1,8-dihydroxy-9-anthrones so they are also expected to affect the keto-enol equilibrium. It is also important to study the effect of intramolecular hydrogen bonds on the structure of both tautomeric forms. It was found that the nature of the substituent in the middle ring could affect the antioxidant properties of the investigated compound.

Favipiravir-Tautomeric and Complexation Properties in Solution

The tautomeric properties of favipiravir were investigated experimentally for the first time by using molecular spectroscopy (UV-Vis absorption, fluorescence and NMR), as well as DFT quantum-chemical calculations. According to the obtained results, the enol tautomer is substantially more stable in most of the organic solvents. In the presence of water, a keto form appears to be favored due to the specific solute-solvent interactions. Upon the addition of alkaline-earth-metal ions, deprotonation and complexation occurred simultaneously, giving the formation of 2 : 1 ligand : metal complexes. According to the theoretical simulations, the metal ion is captured between the carbonyl groups as a result of the size-fit effect.

A synthetically benign one-pot construction of enamino-xanthene dyes

Polyhydroxylated phenols are components of biomass and precursors of pigments in plants. This paper reports a novel entry to xanthene dyes, involving the reaction of 2,4,6-trihydroxybenzaldehyde with primary aliphatic amines. This catalyst-free synthesis exhibits a high atom economy and can be conducted under eco-friendly conditions and operational simplicity.

α-Helix in Cystathionine β-Synthase Enzyme Acts as an Electron Reservoir

The modulation of electron density at the Pyridoxal 5'-phosphate (PLP) catalytic center, because of charge transfer across the α-helix/PLP interface, is the determining factor for the enzymatic activities in the human Cystathionine β-Synthase (hCBS) enzyme. Applying density functional theory calculations, in conjunction with the real space density analysis, we investigated the charge density delocalization across the entire heme-α-helix-PLP electron communication channels. The electron delocalization due to hydrogen bonds at the heme/α-helix and α-helix/PLP interfaces are found to be extended over a very long range, as a result of redistribution of electron densities of the cofactors. Moreover, the internal hydrogen bonds of α-helix that are crucial for its secondary structure also participate in the electron redistribution through the structured hydrogen-bond network. α-Helix is found to accumulate the electron density at the ground state from both of the cofactors and behaves as an electron reservoir for catalytic reaction at the electrophilic center of PLP.

Controlling the Isomerization of Photoresponsive Molecules through a Limiting Tautomerization Strategy

Controlling the multistage photoresponsivity remains a challenge, in part, due to the spontaneous tautomerization between isomers. Herein, we present a strategy to access three independent states (linear, cyclic keto, and cyclic enolate) of crown ether (CE)-substituted donor-acceptor Stenhouse adducts (DASAs) by limiting the tautomerization of the closed isomers. The linear-cyclic keto isomerization is reversibly triggered by treatment with metal ions (Na⁺ or K⁺) and CE, while the linear-cyclic enolate isomerization is induced by green light and heat. Density functional theory and molecular dynamics calculation results suggest that the steric effect and supramolecular interaction between the electron-donating and electron-withdrawing moieties play an important role in hindering the tautomerization between cyclic keto and cyclic enolate DASA-CE. The strategy to influence key steps in the photoswitching process inspires well-controlled multistage isomerization of photoresponsive molecules.

Targeting Individual Tautomers in Equilibrium by Resonant Inelastic X-ray Scattering

Tautomerism is one of the most important forms of isomerism, owing to the facile interconversion between species and the large differences in chemical properties introduced by the proton transfer connecting the tautomers. Spectroscopic techniques are often used for the characterization of tautomers. In this context, separating the overlapping spectral response of coexisting tautomers is a long-standing challenge in chemistry. Here, we demonstrate that by using resonant inelastic X-ray scattering tuned to the core excited states at the site of proton exchange between tautomers one is able to experimentally disentangle the manifold of valence excited states of each tautomer in a mixture. The technique is applied to the prototypical keto-enol equilibrium of 3-hydroxypyridine in aqueous solution. We detect transitions from the occupied orbitals into the LUMO for each tautomer in solution, which report on intrinsic and hydrogen-bond-induced orbital polarization within the π and σ manifolds at the proton-transfer site.

Catalytic Amination of Phenols with Amines

Given the wide prevalence and ready availability of both phenols and amines, aniline synthesis through direct coupling between these starting materials would be extremely attractive. Herein, we describe a rhodium-catalyzed amination of phenols, which provides concise access to diverse anilines, with water as the sole byproduct. The arenophilic rhodium catalyst facilitates the inherently difficult keto-enol tautomerization of phenols by means of π-coordination, allowing for the subsequent dehydrative condensation with amines. We demonstrate the generality of this redox-neutral catalysis by carrying out reactions of a large array of phenols with various electronic properties and a wide variety of primary and secondary amines. Several examples of late-stage functionalization of structurally complex bioactive molecules, including pharmaceuticals, further illustrate the potential broad utility of the method.

Theoretical Study of the Geometry of Dibenzoazepine Analogues

The geometry of dibenzoazepine analogues—typical multifunctional drugs—was investigated to find the geometrical parameters sensitive to the substitution of the central seven-membered ring. Exploration of the crystal structure database (CSD) shows that the geometrical parameter sensitive to the substitution of the carbon atom distance of the central ring not included in the aromatic rings to the plane through the carbon atoms common for the central ring and the aromatic side rings. Presence of the double bond in the central ring was reflected in its partial aromaticity expressed by the HOMED parameter. Some derivatives of 5H-dibenzo[b,f]azepine with flat conformation of the central ring are characterized by mobility of the electron density comparable to the mobility in the aromatic side rings. Influence of the surrounding on the investigated compounds was confirmed by comparison of the optimized molecules and the molecules in the crystal state where the packing forces can influence the molecular geometry.

Push–Pull Effect on the Gas-Phase Basicity of Nitriles: Transmission of the Resonance Effects by Methylenecyclopropene and Cyclopropenimine π-Systems Substituted by Two Identical Strong Electron Donors

The gas-phase basicity of nitriles can be enhanced by a push–pull effect. The role of the intercalated scaffold between the pushing group (electron-donor) and the pulling (electron-acceptor) nitrile group is crucial in the basicity enhancement, simultaneously having a transmission function and an intrinsic contribution to the basicity. In this study, we examine the methylenecyclopropene and the N-analog, cyclopropenimine, as the smallest cyclic π systems that can be considered for resonance propagation in a push–pull system, as well as their derivatives possessing two strong pushing groups (X) attached symmetrically to the cyclopropene scaffold. For basicity and push–pull effect investigations, we apply theoretical methods (DFT and G2). The effects of geometrical and rotational isomerism on the basicity are explored. We establish that the protonation of the cyano group is always favored. The push–pull effect of strong electron donor X substituents is very similar and the two π-systems appear to be good relays for this effect. The effects of groups in the two cyclopropene series are found to be proportional to the effects in the directly substituted nitrile series X–C≡N. In parallel to the basicity, changes in electron delocalization caused by protonation are also assessed on the basis of aromaticity indices. The calculated proton affinities of the nitrile series reported in this study enrich the gas-phase basicity scale of nitriles to around 1000 kJ mol−1.

Tautomeric equilibrium, proton affinity and mass spectrometry fragmentation of flexible hydrogen-bonded precursors and rigid [Formula: see text] fluorescent dyes

The stability of two groups of conformationally locked molecules, similar in topology, but differing only by the type of the bridge rigidifying their structure, is studied. The series of the less-rigid 2-phenacylheterocyclic compounds and their stiff difluoroboranyl derivatives are investigated for the determination of the effect of [Formula: see text]/S/O replacement in a five-membered heterocyclic ring and the presence of a strong electron-donating group on the tautomeric equilibrium, protonation affinity, and fragmentation pattern observed in the structural elucidation by means of mass spectrometry technique. The results of the [Formula: see text]B97X-D/6-311++G(d,p) calculations, the topological analysis of electron density as well as the experimental MS measurements show the importance of the number of heteroatoms, their properties, and location in the molecule for the rational design of the systems of desired stable tautomers or the favorable protonation sites. The obtained data allow for the understanding of the fundamentals of the novel highly fluorescent difluoroborates fragmentation behavior, vital for their structural elucidation with the application of high-resolution tandem mass spectrometry methods.

Evidence of IR-Induced Chemistry in a Neat Solid: Tautomerization of Thiotropolone by Thermal, Electronic, and Vibrational Excitations

Thiotropolone isolated in argon and xenon matrices (as monomers) or in a neat solid (as the crystalline or amorphous state) at low temperature was found to exist only in the thione-enol form. Visible light irradiation (λ ≥ 400 nm) leads to thione-enol → thiol-keto tautomerization in matrices and under neat solid conditions at 15 K. The assignment of the IR spectra of the two thiotropolone tautomers (thione-enol and thiol-keto) was carried out with the support of B3LYP/6-311+G(2d,p) computations. The thiol-keto form generated in situ in a neat solid was found to tautomerize back to the thione-enol upon annealing up to 100 K. Gaussian-4 (G4) computations estimate that such a tautomerization process has an energy barrier of ∼25 kJ mol^-1, which is consistent with the observations. Moreover, it was found that narrowband IR irradiation of the thiol-keto form in a neat solid, at the frequency of its CH stretching overtones/combination modes, also induces tautomerization to the thione-enol form. Such a result constitutes an important demonstration of vibrationally induced chemistry under neat solid conditions.

Asymmetric Enamide-Imine Tautomerism in the Kinetic Resolution of Tertiary Alcohols

An efficient protocol for kinetic resolution of tertiary alcohols has been developed through an unprecedented asymmetric enamide-imine tautomerism process enabled by chiral phosphoric acid catalysis. A broad range of racemic 2-arylsulfonamido tertiary allyl alcohols could be kinetically resolved with excellent kinetic resolution performances (with s-factor up to >200). This method is particularly effective for a series of 1,1-dialkyl substituted allyl alcohols, which produced chiral tertiary alcohols that would be difficult to access via other asymmetric methods. Facile and versatile transformations of the chiral α-hydroxy imine and enamide products, especially the efficient stereodivergent synthesis of all four stereoisomers of β-amino tertiary alcohols using one enantiomer of the catalyst, demonstrated the value of this kinetic resolution method.

Impact of the Characteristics of Quantum Chemical Databases on Machine Learning Prediction of Tautomerization Energies

An essential aspect for adequate predictions of chemical properties by machine learning models is the database used for training them. However, studies that analyze how the content and structure of the databases used for training impact the prediction quality are scarce. In this work, we analyze and quantify the relationships learned by a machine learning model (Neural Network) trained on five different reference databases (QM9, PC9, ANI-1E, ANI-1, and ANI-1x) to predict tautomerization energies from molecules in Tautobase. For this, characteristics such as the number of heavy atoms in a molecule, number of atoms of a given element, bond composition, or initial geometry on the quality of the predictions are considered. The results indicate that training on a chemically diverse database is crucial for obtaining good results and also that conformational sampling can partly compensate for limited coverage of chemical diversity. The overall best-performing reference database (ANI-1x) performs on average by 1 kcal/mol better than PC9, which, however, contains about 2 orders of magnitude fewer reference structures. On the other hand, PC9 is chemically more diverse by a factor of ∼5 as quantified by the number of atom-in-molecule-based fragments (amons) it contains compared with the ANI family of databases. A quantitative measure for deficiencies is the Kullback-Leibler divergence between reference and target distributions. It is explicitly demonstrated that when certain types of bonds need to be covered in the target database (Tautobase) but are undersampled in the reference databases, the resulting predictions are poor. Examples of this include the poor performance of all databases analyzed to predict C(sp²)-C(sp²) double bonds close to heteroatoms and azoles containing N-N and N-O bonds. Analysis of the results with a Tree MAP algorithm provides deeper understanding of specific deficiencies in predicting tautomerization energies by the reference datasets due to inadequate coverage of chemical space. Capitalizing on this information can be used to either improve existing databases or generate new databases of sufficient diversity for a range of machine learning (ML) applications in chemistry.

Pyrrole-bridged quinones: π-electronic systems that modulate electronic structures by tautomerism and deprotonation

A new series of π-extended quinone derivatives containing a pyrrole bridge exhibited NH/OH-type tautomerization and anion binding along with deprotonation that induced near-infrared absorption and ion-pairing assemblies.

Facile room-temperature self-assembly of extended cation-free guanine-quartet network on Mo-doped Au(111) surface

Guanine-quadruplex, consisting of several stacked guanine-quartets (GQs), has emerged as an important category of novel molecular targets with applications from nanoelectronic devices to anticancer drugs. Incorporation of metal cations into a GQ structure is utilized to form stable G-quadruplexes, while formation of a cation-free GQ network has been challenging. Here we report the room temperature (RT) molecular self-assembly of extended pristine GQ networks on an Au(111) surface. An implanted molybdenum atom within the Au(111) surface is used to nucleate and stabilize the cation-free GQ network. Additionally, decoration of the Au(111) surface with 7-armchair graphene nanoribbons (7-AGNRs) enhances the GQ domain size by suppressing the influence of the disordered phase nucleated from Au step edges. Scanning tunneling microscopy/spectroscopy (STM/STS) and density functional theory (DFT) calculations confirm the formation of GQ networks and unravel the nucleation and growth mechanism. Our work, utilizing a hetero-atom doped substrate, provides a facile approach to enhance the stability and domain size of the GQ self-assembly, which would be applicable for other molecular structures.

The thermochemical, structural, and spectroscopic analyses of the tautomers of sulfur and selenium modified emissive nucleobases

In this study, density functional theory (DFT) and time dependent density functional theory (TD-DFT) are used to investigate the stabilities and spectral properties [IR, UV–vis, and two-photon absorption (2PA)] of two sets of modified RNA nucleobase tautomers. The modifications introduce either a sulfur or selenium atom to form an isothiazolo[4,3-d]pyrimidine or isoselenazolo[4,3-d]pyrimidine heterocyclic core, respectively. The relative stabilities of both sets of modified tautomers determined with B3LYP/6-31++G(d,p) reveal that in water (with a polarizable continuum model), the 6-keto-2-amino tautomer of guanine and the rare 4-imino-2-keto tautomer of cytosine may be present at significant populations, whereas the 6-enol-2-amino tautomer of guanine is more common in the gas phase. The identification of these modified tautomers is possible due to the natural differences in their vibrational modes and hence IR spectra. Furthermore, the photophysical properties of both these sets of modified tautomers indicate that excitation and emission energies are shifted relative to their more abundant form in both one photon absorption and emission and 2PA spectra, as determined at the B3LYP/6-31++G(d,p) and CAM-B3LYP/aug-cc-pVDZ levels of theory, respectively. Even though the 2PA cross sections in water for all of the species are small (0.3–2.3 GM), the modified cytosine tautomer shows promise, as its cross section is larger than the more dominant form. The spectral separation between the dominant form and the tautomers of isoselenazole and isothiazole modified cytosine and guanine are relatively similar, suggesting both modifications could be useful in elucidating the tautomers from their more abundant counterparts.

Mono- and Diamination of 4,6-Dichloropyrimidine, 2,6-Dichloropyrazine and 1,3-Dichloroisoquinoline with Adamantane-Containing Amines

N-heteroaryl substituted adamantane-containing amines are of substantial interest for their perspective antiviral and psychotherapeutic activities. Chlorine atom at alpha-position of N-heterocycles has been substituted by the amino group using convenient nucleophilic substitution reactions with a series of adamantylalkylamines. The prototropic equilibrium in these compounds was studied using NMR spectroscopy. The introduction of the second amino substituent in 4-amino-6-chloropyrimidine, 2-amino-chloropyrazine, and 1-amino-3-chloroisoquinoline was achieved using Pd(0) catalysis.

Finding the direct energy-structure correlations in intramolecular aromaticity assisted hydrogen bonding (AAHB)

A predictive model for intramolecular hydrogen bond energy (E_HB) calculation of polyaromatic ortho-hydroxyaldehydes based on a set of small, functionalized hydrocarbons is developed. The complete data set of 18 compounds was used for this study. The model is based on one of four optional categories of molecular descriptors: geometric, spectroscopic, bond order and topological indices. The model of Wiberg bond indices (WBIs) as descriptors of the CC involved bond based on stepwise regression has acceptable prediction abilities for 14 structures of ortho-hydroxyformylobenzo[a]pyrene derivatives already at the semiempirical level. The presented correlation enables a significantly more rapid and quantitative description of the hydrogen bonding strength than the much more time-consuming MTA method. Thus, WBIs are shown to provide a reliable means for fast prescreening of the energy of chelate hydrogen bonds potentially for any polyaromatic derivatives.

Biomimetic hydrogen-bonding cascade for chemical activation: telling a nucleophile from a base

Hydrogen bonding-assisted polarization is an effective strategy to promote bond-making and bond-breaking chemical reactions. Taking inspiration from the catalytic triad of serine protease active sites, we have devised a conformationally well-defined benzimidazole platform that can be systematically functionalized to install multiple hydrogen bonding donor (HBD) and acceptor (HBA) pairs in a serial fashion. We found that an increasing number of interdependent and mutually reinforcing HBD–HBA contacts facilitate the bond-forming reaction of a fluorescence-quenching aldehyde group with the cyanide ion, while suppressing the undesired Brønsted acid–base reaction. The most advanced system, evolved through iterative rule-finding studies, reacts rapidly and selectively with CN⁻ to produce a large (>180-fold) enhancement in the fluorescence intensity at λ_max = 450 nm.

Biomimetic cascade hydrogen bonds promote covalent capture of a nucleophile by polarizing the electrophilic reaction site, while suppressing non-productive acid–base chemistry as the competing reaction pathway.

Synthesis of Porphyrinquinone and Doubly-Fused Diporphyrin Quinone Through Oxidation of Diarylporphyrins Using a Hypervalent Iodine Compound

Treatment of a meso-diarylporphyrin with PhI(OAc)₂ in the presence of BF₃  ⋅ OEt₂ and propionic acid affords the corresponding porphyrinquinone in a high yield (91%). A novel quinone derived from meso-meso β-β doubly-fused diporphyrin was obtained as the sole byproduct (16% yield), which exhibits strong panchromatic absorption between 300 and 1000 nm. It has a low HOMO-LUMO gap owing to expanded and low-symmetry π-planes.

Imine or Enamine? Insights and Predictive Guidelines from the Electronic Effect of Substituents in H-Bonded Salicylimines

Imine and enamine bonds decorate the skeleton of numerous reagents, catalysts, and organic materials. However, it is difficult to isolate at will a single tautomer, as dynamic equilibria occur easily, even in the solid state, and are sensitive to electronic and steric effect, including π-conjugation and H-bonding. Here, using as model Schiff bases generated from salicylaldehydes and TRIS in a set of linear free energy relationships (LFER), we disclose how the formation of either imines or enamines can be controlled and provide a comprehensive framework that captures the structural underpinning of this prediction. This work highlights the potentiality of tailor-made designs en route to compounds with desirable functionality.

Purine tautomeric preferences and bond-length alternation in relation with protonation-deprotonation and alkali metal cationization

Quantum chemical calculations were carried out for deprotonated (P-) and protonated purine (PH+) and for adducts with one alkali metal cation (P-M+ and PM+, where M+ is Li+ or Na+) in the gas phase {B3LYP/6-311+G(d,p)}, a model of perfectly apolar environment, and for selected structures in aqueous solution {PCM(water)//B3LYP/6-311+G(d,p)}, a reference polar medium for biological studies. All potential isomers of purine derivatives were considered, the favored structures indicated, and the preferred sites for protonation/deprotonation and cationization reactions determined. Proton and metal cation basicities of purine in the gas phase were discussed and compared with those of imidazole and pyrimidine. Bond-length alternations in the P, PH+, P-M+, and PM+ forms were quantitatively measured using the harmonic oscillator model of electron delocalization (HOMED) indices and compared with those for P. Variations of the HOMED values when proceeding from the purine structural building blocks, pyrimidine and imidazole, to the bicyclic purine system were also examined. Generally, the isolated NH isomers exhibit a strongly delocalized π-system (HOMED > 0.8). Deprotonation slightly increases the HOMED values, whereas protonation and cationization change the HOMED indices in different way. For bidentate M+-adducts, the HOMED values are larger than 0.9 like for the largely delocalized P-. The HOMED values correlate well in a comprehensive relationship with the relative Gibbs energies (ΔG) calculated for individual isomers whatever the purine form is, neutral, protonated, or cationized. When PCM-DFT model was utilized for P-, PH+, PM+, and P-M+ (M+ = Li+) both electron delocalization and relative stability are different from those for the molecules in vacuo. The solvation effects cause a slight increase in HOMEDs, whereas the ΔEs decrease, but in different ways. Hence, contribution of particular isomers in the isomeric mixtures of PH+, PM+, and P-M+ also varies. HOMED variations for the favored neutral, deprotonated, protonated, and lithiated forms of purine in the gas phase and aqueous solution.

Tautomeric Equilibria of Nucleobases in the Hachimoji Expanded Genetic Alphabet

Evolution has yielded biopolymers that are constructed from exactly four building blocks and are able to support Darwinian evolution. Synthetic biology aims to extend this alphabet, and we recently showed that 8-letter (hachimoji) DNA can support rule-based information encoding. One source of replicative error in non-natural DNA-like systems, however, is the occurrence of alternative tautomeric forms, which pair differently. Unfortunately, little is known about how structural modifications impact free-energy differences between tautomers of the non-natural nucleobases used in the hachimoji expanded genetic alphabet. Determining experimental tautomer ratios is technically difficult, and so, strategies for improving hachimoji DNA replication efficiency will benefit from accurate computational predictions of equilibrium tautomeric ratios. We now report that high-level quantum-chemical calculations in aqueous solution by the embedded cluster reference interaction site model, benchmarked against free-energy molecular simulations for solvation thermodynamics, provide useful quantitative information on the tautomer ratios of both Watson-Crick and hachimoji nucleobases. In agreement with previous computational studies, all four Watson-Crick nucleobases adopt essentially only one tautomer in water. This is not the case, however, for non-natural nucleobases and their analogues. For example, although the enols of isoguanine and a series of related purines are not populated in water, these heterocycles possess N₁-H and N₃-H keto tautomers that are similar in energy, thereby adversely impacting accurate nucleobase pairing. These robust computational strategies offer a firm basis for improving experimental measurements of tautomeric ratios, which are currently limited to studying molecules that exist only as two tautomers in solution.

A Hidden Side of the Conformational Mobility of the Quercetin Molecule Caused by the Rotations of the O3H, O5H and O7H Hydroxyl Groups: In Silico Scrupulous Study

In this study at the MP2/6-311++G(2df,pd)//B3LYP/6-311++G(d,p) level of quantum-mechanical theory it was explored conformational variety of the isolated quercetin molecule due to the mirror-symmetrical hindered turnings of the O3H, O5H and O7H hydroxyl groups, belonging to the A and C rings, around the exocyclic C–O bonds. These dipole active conformational transformations proceed through the 72 transition states (TSs; C1 point symmetry) with non-orthogonal orientation of the hydroxyl groups relatively the plane of the A or C rings of the molecule (HO7C7C8/HO7C7C6 = ±(89.9–93.3), HO5C5C10 = ±(108.9–114.4) and HO3C3C4 = ±(113.6–118.8 degrees) (here and below signs ‘±’ corresponds to the enantiomers)) with Gibbs free energy barrier of activation ΔΔGTS in the range 3.51–16.17 kcal·mol−1 under the standard conditions (T = 298.1 K and pressure 1 atm): ΔΔGTSO7H (3.51–4.27) < ΔΔGTSO3H (9.04–11.26) < ΔΔGTSO5H (12.34–16.17 kcal mol−1). Conformational dynamics of the O3H and O5H groups is partially controlled by the intramolecular specific interactions O3H…O4, C2′/C6′H…O3, O3H…C2′/C6′, O5H…O4 and O4…O5, which are flexible and cooperative. Dipole-active interconversions of the enantiomers of the non-planar conformers of the quercetin molecule (C1 point symmetry) is realized via the 24 TSs with C1 point symmetry (HO3C3C2C1 = ±(11.0–19.1), HC2′/C6′C1′C2 = ±(0.6–2.9) and C3C2C1′C2′/C3C2C1′C6′ = ±(1.7–9.1) degree; ΔΔGTS = 1.65–5.59 kcal·mol−1), which are stabilized by the participation of the intramolecular C2′/C6′H…O1 and O3H…HC2′/C6′ H-bonds. Investigated conformational rearrangements are rather quick processes, since the time, which is necessary to acquire thermal equilibrium does not exceed 6.5 ns.

A dielectric and spectrophotometric study of the tautomerization of 2-hydroxypyridine and 2-mercaptopyridine in water

The basic ionization (pk₁) and acidic ionization (pk₂) constants and equilibrium constant (K_T) of 2HPy and 2MPy were determined. The pk₁(s) of their N- and X-methyl derivatives (X = O, S) were also determined. The equilibrium constant of 2MPy is approximately 60 times greater than its oxygen analog, 2HPy. The micro-ionization constants of the functional groups, –NH (pk_A and pk_C) and –XH (pk_B and pk_D), were determined to provide further insights into the ionization equilibria of these N-heteroaromatic XH compounds (2HPy and 2MPy). The relaxation time of water (τ) in aqueous solutions of 2HPy and 2MPy are collectively used with the K_T values to determine the forward (k_f) and backward (k_b) rate constants of tautomerization. Subsequently, the k_f and k_b are used to provide the rationale for the K_T and τ values.

The basic ionization (pk₁) and acidic ionization (pk₂) constants and equilibrium constant (K_T) of 2HPy and 2MPy were determined.

Using Chou’s 5-Step Rule to Evaluate the Stability of Tautomers: Susceptibility of 2-[(Phenylimino)-methyl]-cyclohexane-1,3-diones to Tautomerization Based on the Calculated Gibbs Free Energies

Gibbs free energies, based on DFT (Density Functional Theory) calculations, prove that enaminone (2-(anilinemethylidene)cyclohexane-1,3-dione) and ketamine (2-[(phenylimino)-methyl]cyclohexane-1,3-dione) are the most and least stable tautomeric forms of the studied systems, respectively. 1H and 13C NMR spectra prove that 2-(anilinemethylidene)cyclohexane-1,3-diones are the only tautomeric species present in dimethylsulfoxide solution (a very weak signal can be seen only for the p-methoxy derivatives). The zwitterionic character of these enaminones is strengthened by naphthoannulation and by the insertion of the electron-withdrawing substituent into the benzene ring (the latter weakens the intramolecular hydrogen bond in the compound). Substituent and naphtoannulation have no effect on the stability of the studied tautomers. Slight twisting of the benzene ring, with respect to the CArNC plane (seen in the crystalline state), was proven to also take place in vacuum and in solution.

Electronic effects in tautomeric equilibria: the case of chiral imines from d-glucamine and 2-hydroxyacetophenones

A one-pot procedure for preparing a series of chiral imines by direct condensation of d-glucamine with 2-hydroxyacetophenones is described. Under conventional acetylation an unexpected mixture of two different peracetylated molecules is obtained, one with an open enamine structure, and the other incorporating an N-acetyl-1,3-oxazolidine into the acyclic skeleton. Surprisingly, both molecules coexist within the crystal's unit cell, as inferred from single-crystal X-ray analysis of a 5-bromo-substituted aryl derivative. Moreover, the 1,3-oxazolidine ring exists as rotational conformers (E,Z) owing to the restricted rotation around the N-acetyl bond. The equilibrium involving imine and enamine structures has been assessed in detail, providing in addition linear free-energy relationships between the tautomerization constants (K_T) and the electronic effect of the substituents.

Revisiting the Structure and Chemistry of 3(5)-Substituted Pyrazoles

Pyrazoles are known as versatile scaffolds in organic synthesis and medicinal chemistry, often used as starting materials for the preparation of more complex heterocyclic systems with relevance in the pharmaceutical field. Pyrazoles are also interesting compounds from a structural viewpoint, mainly because they exhibit tautomerism. This phenomenon may influence their reactivity, with possible impact on the synthetic strategies where pyrazoles take part, as well as on the biological activities of targets bearing a pyrazole moiety, since a change in structure translates into changes in properties. Investigations of the structure of pyrazoles that unravel the tautomeric and conformational preferences are therefore of upmost relevance. 3(5)-Aminopyrazoles are largely explored as precursors in the synthesis of condensed heterocyclic systems, namely pyrazolo[1,5-a]pyrimidines. However, the information available in the literature concerning the structure and chemistry of 3(5)-aminopyrazoles is scarce and disperse. We provide a revision of data on the present subject, based on investigations using theoretical and experimental methods, together with the applications of the compounds in synthesis. It is expected that the combined information will contribute to a deeper understanding of structure/reactivity relationships in this class of heterocycles, with a positive impact in the design of synthetic methods, where they take part.

Acid-catalyzed synthesis of condensed polycyclic diaryl ethers from arenols

Diaryl ethers containing condensed polycyclic aryl groups were synthesized from arenols in the presence of a catalytic amount of p-chlorobenzenesulfonic acid. Symmetrical binaphthyl, biphenanthryl, and bipyrenyl ethers were obtained in high yields. Unsymmetrical derivatives were also synthesized from 9-phenanthrol and arenols, using combinations of reactive and unreactive substrates.