A review on point mutations via proton transfer in DNA base pairs in the absence and presence of electric fields.

This comprehensive review focuses on spontaneous mutations that may occur during DNA replication, the fundamental process responsible for transferring genetic information. In 1963, Löwdin postulated that these mutations are primarily a result of proton transfer reactions within the hydrogen-bonded DNA base pairs. The single and double proton transfer reactions within the base pairs in DNA result in zwitterions and rare tautomers, respectively. For persistent mutations, these products must be generated at high rates and should be thermodynamically stable. This review covers the proton transfer reactions studied experimentally and computationally. The review also examines the influence of externally applied electric fields on the thermodynamics and kinetics of proton transfer reactions within DNA base pairs, and their biological implications.

When Tautomers Matter: UV-Vis Absorption Spectra of Hypoxanthine in Aqueous Solution from Fully Atomistic Simulations.

The UV-Vis spectrum of the solvated purine derivative Hypoxanthine (HYX) is investigated using the Quantum Mechanics/Fluctuating Charges (QM/FQ) multiscale approach combined with a sampling of configurations through atomistic Molecular Dynamics (MD) simulations. Keto 1H7H and 1H9H tautomeric forms of HYX are the most stable in aqueous solution and form different stable complexes with the surrounding water molecules, ultimately affecting the electronic absorption spectra. The final simulated spectrum resulting from the combination of the individual spectra of tautomers agrees very well with most of the characteristics in the measured spectrum. The importance of considering the effect of the solute tautomers and, in parallel, the contribution of the different solvent arrangements around the solute when modeling spectral properties, is highlighted. In addition, the high quality of the computed spectra leads to suggesting an alternative way for acquiring tautomeric populations from combined computational/experimental spectra.

Significance of Triazole in Medicinal Chemistry: Advancement in Drug Design, Reward and Biological Activity.

One of the triazole tautomers, 1,2,4-triazole derivatives, has a wide range of biological activities that suggest its potential therapeutic utility in medicinal chemistry. These actions include anti-inflammatory, anti-cancer, anti-bacterial, anti-tuberculosis, and anti-diabetic effects. Using computational simulations and models, we investigate the structure-activity relationships of 1,2,4-triazoles, showing how various modifications to the triazole core yield a variety of clinical therapeutic benefits. The review highlights the anti-inflammatory effect of 1,2,4-triazoles in relation to their ability to disrupt significant inflammatory mediators and pathways. We present in-silico data that illuminate the triazoles' capacity to inhibit cell division, encourage apoptosis, and stop metastasis in a range of cancer models. This review looks at the bactericidal and bacteriostatic properties of 1,2,4-triazole derivatives, with a focus on their potential efficacy against multi-drug resistant bacterial infections and their usage in tuberculosis therapy. In order to better understand these substances' potential anti-diabetic benefits, this review also looks at how they affect glucose metabolism regulation and insulin responsiveness. Coordinated efforts are required to translate the efficacy of 1,2,4-triazole compounds in preclinical models into practical therapeutic benefits. Based on the information provided, it can be concluded that 1,2,4-triazole derivatives are a promising class of diverse therapeutic agents with potential utility in a range of disorders. Their development and improvement might herald a new era of medical care that will be immensely advantageous to both patients and the medical community as a whole. This comprehensive research, which is further reinforced by in-silico investigations, highlights the great medicinal potential of 1,2,4-triazoles. Additionally, this study encourages more research into these substances and their enhancement for use in pharmaceutical development.

Synthesis, In Silico and Kinetics Evaluation of N-(ß-d-glucopyranosyl)-2-arylimidazole-4(5)-carboxamides and N-(ß-d-glucopyranosyl)-4(5)-arylimidazole-2-carboxamides as Glycogen Phosphorylase Inhibitors.

Recently studied N-(ß-d-glucopyranosyl)-3-aryl-1,2,4-triazole-5-carboxamides have proven to be low micromolar inhibitors of glycogen phosphorylase (GP), a validated target for the treatment of type 2 diabetes mellitus. Since in other settings, the bioisosteric replacement of the 1,2,4-triazole moiety with imidazole resulted in significantly more efficient GP inhibitors, in silico calculations using Glide molecular docking along with unbound state DFT calculations were performed on N-(ß-d-glucopyranosyl)-arylimidazole-carboxamides, revealing their potential for strong GP inhibition. The syntheses of the target compounds involved the formation of an amide bond between per-O-acetylated ß-d-glucopyranosylamine and the corresponding arylimidazole-carboxylic acids. Kinetics experiments on rabbit muscle GPb revealed low micromolar inhibitors, with the best inhibition constants (Kis) of ~3-4 µM obtained for 1- and 2-naphthyl-substituted N-(ß-d-glucopyranosyl)-imidazolecarboxamides, 2b-c. The predicted protein-ligand interactions responsible for the observed potencies are discussed and will facilitate the structure-based design of other inhibitors targeting this important therapeutic target. Meanwhile, the importance of the careful consideration of ligand tautomeric states in binding calculations is highlighted, with the usefulness of DFT calculations in this regard proposed.

Overcoming NMR line broadening of nitrogen containing compounds: A simple solution.

This study presents a straightforward solution to the challenge of elucidating the structures of nitrogen containing compounds undergoing isomerization. When spectral line broadening occurs related to isomerization, be it prototropic tautomerism or bond rotations, this poses a significant obstacle to structural elucidation. By adding acids, we demonstrate a simple approach to overcome this issue and effectively sharpen NMR signals for acid stable prototropic tautomers as well as the conformational isomers containing a morpholine or piperazine ring.

On the Potential Role of the (Pseudo-) Jahn-Teller Effect in the Membrane Transport Processes: Enniatin B and Beauvericin.

The molecular structure of mycotoxins enniatin B and beauvericin, which are used as ionophores, was studied using density functional theory in various symmetry groups and singly charged states. We have shown that the charge addition or removal causes significant structural changes. Unlike the neutral C3 molecules, the stability of the charged C1 structures was explained by the Jahn-Teller or Pseudo-Jahn-Teller effect. This finding agrees with the available experimental X-ray structures of their metal complexes where electron density transfer from the metal can be expected. Hence, the membrane permeability of metal sandwich-structure complexes possessing antimicrobial activities is modulated by the conformational changes.

A DFT study on non-enzymatic degradations of anti-tuberculosis drug isoniazid.

CONTEXT: Isoniazid (INH) is one of the medications most used for tuberculosis (TB) treatment. However, long-term continuous therapy can cause hepatotoxicity and peripheral neuritis. The degradation of INH is an important aspect of the research in the field of drug stability as well as drug formulation for controlling release. It is thought that tautomerization, hydrolysis as well as nucleophilic substitutions can cause decrease in INH as non-enzymatic degradation. Therefore, it is crucial to understand the mechanisms and energies of the major reactions in order to provide reference for future drug formulation and application. This study is an effort to understand the kinetic and thermodynamic properties of the non-enzymatic degradation reactions. The chemical reaction phenomena are investigated using the density functional theory (DFT) method. This study shows that major degradation of INH can be done via tautomerization followed by hydrolysis. The general trends in nucleophilic degradation presented here are consistent with experimental pKa of nucleophiles. METHODS: All DFT calculations were performed using the Gaussian Software Packages (Gaussian 09 revision B.01 and GaussView 5.0.8). MOLEKEL 4.3 software was utilized to visualize the molecular graphics of all relevant species. The optimized molecular geometries were calculated using B3LYP/6-311 + G(d,p) level in the gas phase. The IEF-PCM/B3LYP/6-311 + G(d,p) level was selected for single-point and frequency calculations in aqueous media.

Computational Analysis of Histamine Protonation Effects on H1R Binding.

Despite numerous studies investigating histamine and its receptors, the impact of histamine protonation states on binding to the histamine H1-receptor (H1R) has remained elusive. Therefore, we assessed the influence of different histamine tautomers (τ-tautomer, π-tautomer) and charge states (mono- vs. dicationic) on the interaction with the ternary histamine-H1R-Gq complex. In atomistic molecular dynamics simulations, the τ-tautomer formed stable interactions with the receptor, while the π-tautomer induced a rotation of the histamine ring by 180° and formed only weaker hydrogen bonding interactions. This suggests that the τ-tautomer is more relevant for stabilization of the active ternary histamine-H1R-Gq complex. In addition to the two monocationic tautomers, the binding of dicationic histamine was investigated, whose interaction with the H1R had been observed in a previous experimental study. Our simulations showed that the dication is less compatible with the ternary histamine-H1R-Gq complex and rather induces an inactive conformation in the absence of the Gq protein. Our data thus indicate that the charge state of histamine critically affects its interactions with the H1R. Ultimately these findings might have implications for the future development of new ligands that stabilize distinct H1R activation states.

How Does Pseudo-Jahn-Teller Effect Induce the Photoprotective Potential of Curcumin?

In this paper, the molecular and electronic structure of curcumin is studied. High-symmetric gas-phase tautomers and their deprotonated forms in various symmetry groups are identified. The stability of lower-symmetry structures was explained by using the Pseudo-Jahn-Teller (PJT) effect. This effect leads to stable structures of different symmetries for the neutral enol and keto forms. The presented analysis demonstrated the potential significance of the PJT effect, which may modulate the setting of electronic and vibrational (vibronic) energy levels upon photodynamic processes. The PJT effect may rationalize the photoprotection action and activity of naturally occurring symmetric dyes.

Design and Synthesis of 3-(ß-d-Glucopyranosyl)-4-amino/4-guanidino Pyrazole Derivatives and Analysis of Their Glycogen Phosphorylase Inhibitory Potential.

Glycogen phosphorylase (GP) is a key regulator of glucose levels and, with that, an important target for the discovery of novel treatments against type 2 diabetes. ß-d-Glucopyranosyl derivatives have provided some of the most potent GP inhibitors discovered to date. In this regard, C-ß-d-glucopyranosyl azole type inhibitors proved to be particularly effective, with 2- and 4-ß-d-glucopyranosyl imidazoles among the most potent designed to date. His377 backbone C=O hydrogen bonding and ion-ion interactions of the protonated imidazole with Asp283 from the 280s loop, stabilizing the inactive state, were proposed as crucial to the observed potencies. Towards further exploring these features, 4-amino-3-(ß-d-glucopyranosyl)-5-phenyl-1H-pyrazole (3) and 3-(ß-d-glucopyranosyl)-4-guanidino-5-phenyl-1H-pyrazole (4) were designed and synthesized with the potential to exploit similar interactions. Binding assay experiments against rabbit muscle GPb revealed 3 as a moderate inhibitor (IC50 = 565 µM), but 4 displayed no inhibition at 625 µM concentration. Towards understanding the observed inhibitions, docking and post-docking molecular mechanics-generalized Born surface area (MM-GBSA) binding free energy calculations were performed, together with Monte Carlo and density functional theory (DFT) calculations on the free unbound ligands. The computations revealed that while 3 was predicted to hydrogen bond with His377 C=O in its favoured tautomeric state, the interactions with Asp283 were not direct and there were no ion-ion interactions; for 4, the most stable tautomer did not have the His377 backbone C=O interaction and while ion-ion interactions and direct hydrogen bonding with Asp283 were predicted, the conformational strain and entropy loss of the ligand in the bound state was significant. The importance of consideration of tautomeric states and ligand strain for glucose analogues in the confined space of the catalytic site with the 280s loop in the closed position was highlighted.

Complementary base pair interactions between different rare tautomers of the second-generation artificial genetic alphabets.

The functionality of a semisynthetic DNA in the biological environment will depend on the base pair nature of its complementary base pairs. To understand this, base pair interactions between complementary bases of recently proposed eight second-generation artificial nucleobases are studied herein by considering their rare tautomeric conformations and a dispersion-corrected density functional theoretic method. It is found that the binding energies of two hydrogen-bonded complementary base pairs are more negative than those of the three hydrogen-bonded base pairs. However, as the former base pairs are endothermic, the semisynthetic duplex DNA would involve the latter base pairs.

Aromatic Character and Relative Stability of Pyrazoloporphyrin Tautomers and Related Protonated Species: Insights into How Pyrazole Changes the Properties of Carbaporphyrinoid Systems.

Pyrazoloporphyrins (PzPs), which are porphyrin analogues incorporating a pyrazole subunit, are examples of carbaporphyrin-type structures with a carbon atom within the macrocyclic cavity. DFT calculations were used to assess a series of 17 PzP tautomers, nine monoprotonated species and four related diprotonated PzP dications. The geometries of the structures were optimized using M06-2X/6-311++G(d,p), and the relative stabilities computed with the cc-PVTZ functional. Nucleus independent chemical shifts, both NICS(0) and NICS(1)zz, were calculated, and the anisotropy of the induced current density (AICD) plots were generated for all of the species under investigation. The results for free base PzPs show that fully aromatic PzP tautomers are not significantly more stable than weakly aromatic cross-conjugated species. In addition, strongly aromatic structures with internal CH2's are much less stable, a feature that is also seen for protonated PzPs. The degree of planarity for the individual macrocycles does not significantly correlate with the stability of these structures. The results allow significant aromatic conjugation pathways to be identified in many cases, and provide insights into the aromatic properties of this poorly studied system. These investigations also complement experimental results for PzPs and emphasize the need for further studies in this area.

A DFT Investigation of the Reactivity of Guanidinium Salts in Tandem aza-Michael Addition/Intramolecular Cyclization.

A proposed mechanism of the reaction of guanidinium chlorides with dimethyl acetylenedicarboxylate in a tandem aza-Michael addition reaction/intramolecular cyclization was investigated by DFT M06-2X and B3LYP computational approaches. The energies of the products were compared against the G3, M08-HX, M11, and wB97xD data or experimentally obtained product ratios. The structural diversity of the products was interpreted by the concurrent formation of different tautomers formed in situ upon deprotonation with a 2-chlorofumarate anion. A comparison of relative energies of the characteristic stationary points along the examined reaction paths indicated that the initial nucleophilic addition is energetically the most demanding process. The overall reaction is strongly exergonic, as predicted by both methods, which is primarily due to methanol elimination during the intramolecular cyclization step producing cyclic amide structures. Formation of a five-membered ring upon intramolecular cyclization is highly favored for the acyclic guanidine, while optimal product structure for the cyclic guanidines is based on a 1,5,7-triaza [4.3.0]-bicyclononane skeleton. Relative stabilities of the possible products calculated by the employed DFT methods were compared against the experimental product ratio. The best agreement was obtained for the M08-HX approach while the B3LYP approach provided somewhat better results than the M06-2X and M11 methods.

Intermolecular Interactions of Edaravone in Aqueous Solutions of Ethaline and Glyceline Inferred from Experiments and Quantum Chemistry Computations.

Edaravone, acting as a cerebral protective agent, is administered to treat acute brain infarction. Its poor solubility is addressed here by means of optimizing the composition of the aqueous choline chloride (ChCl)-based eutectic solvents prepared with ethylene glycol (EG) or glycerol (GL) in the three different designed solvents compositions. The slurry method was used for spectroscopic solubility determination in temperatures between 298.15 K and 313.15 K. Measurements confirmed that ethaline (ETA = ChCl:EG = 1:2) and glyceline (GLE = ChCl:GL = 1:2) are very effective solvents for edaravone. The solubility at 298.15 K in the optimal compositions was found to be equal xE = 0.158 (cE = 302.96 mg/mL) and xE = 0.105 (cE = 191.06 mg/mL) for glyceline and ethaline, respectively. In addition, it was documented that wetting of neat eutectic mixtures increases edaravone solubility which is a fortunate circumstance not only from the perspective of a solubility advantage but also addresses high hygroscopicity of eutectic mixtures. The aqueous mixture with 0.6 mole fraction of the optimal composition yielded solubility values at 298.15 K equal to xE = 0.193 (cE = 459.69 mg/mL) and xE = 0.145 (cE = 344.22 mg/mL) for glyceline and ethaline, respectively. Since GLE is a pharmaceutically acceptable solvent, it is possible to consider this as a potential new liquid form of this drug with a tunable dosage. In fact, the recommended amount of edaravone administered to patients can be easily achieved using the studied systems. The observed high solubility is interpreted in terms of intermolecular interactions computed using the Conductor-like Screening Model for Real Solvents (COSMO-RS) approach and corrected for accounting of electron correlation, zero-point vibrational energy and basis set superposition errors. Extensive conformational search allowed for identifying the most probable contacts, the thermodynamic and geometric features of which were collected and discussed. It was documented that edaravone can form stable dimers stabilized via stacking interactions between five-membered heterocyclic rings. In addition, edaravone can act as a hydrogen bond acceptor with all components of the studied systems with the highest affinities to ion pairs of ETA and GLE. Finally, the linear regression model was formulated, which can accurately estimate edaravone solubility utilizing molecular descriptors obtained from COSMO-RS computations. This enables the screening of new eutectic solvents for finding greener replacers of designed solvents. The theoretical analysis of tautomeric equilibria confirmed that keto-isomer edaravone is predominant in the bulk liquid phase of all considered deep eutectic solvents (DES).

Imidazole-amino acids. Conformational switch under tautomer and pH change.

Replacement of the main chain peptide bond by imidazole ring seems to be a promising tool for the peptide-based drug design, due to the specific prototropic tautomeric as well as amphoteric properties. In this study, we present that both tautomer and pH change can cause a conformational switch of the studied residues of alanine (1-4) and dehydroalanine (5-8) with the C-terminal peptide group replaced by imidazole. The DFT methods are applied and an environment of increasing polarity is simulated. The conformational maps (Ramachandram diagrams) are presented and the stability of possible conformations is discussed. The neutral forms, tautomers τ (1) and π (2), adapt the conformations αRτ (φ, ψ = - 75°, - 114°) and C7eq (φ, ψ = - 75°, 66°), respectively. Their torsion angles ψ differ by about 180°, which results in a considerable impact on the peptide chain conformation. The cation form (3) adapts both these conformations, whereas the anion analogue (4) prefers the conformations C5 (φ, ψ = - 165°, - 178°) and ß2 (φ, ψ ~ - 165°, - 3°). Dehydroamino acid analogues, the tautomers τ (5) and π (6) as well as the anion form (8), have a strong tendency toward the conformations ß2 (φ, ψ = - 179°, 0°) and C5 (φ, ψ = - 180°, 180°). The preferences of the protonated imidazolium form (7) depend on the environment. The imidazole ring, acting as a donor or acceptor of the hydrogen bonds created within the studied residues, has a profound effect on the type of conformation.

Reliable gas-phase tautomer equilibria of drug-like molecule scaffolds and the issue of continuum solvation.

Accurate calculation of relative tautomer energies in different environments is a prerequisite to many parameters of relevance in drug discovery. This work provides a thorough benchmark of the semiempirical methods AM1, PM3 and GFN2-xTB, the force-field OPLS4, Hartree-Fock and HF-3c, the density functionals PBEh-3c, B97-3c, r2SCAN-3c, PBE, PBE0, TPSS, r2SCAN, ω-B97X-V, M06-2X, B3LYP, B2PLYP, and second-order perturbation theory MP2 versus the gold-standard coupled-cluster DLPNO-CCSD(T) using the def2-QZVPP basis set. The outperforming method identified is M06-2X, whereas r2SCAN-3c is the best-perfoming one in the set of cost-optimized methods. Application of the two methods on a challenging subset from the SAMPL2 challenge provides evidence that deviations from experiment are caused by deficiencies of current continuum solvation methods.

Intramolecular Interactions in Derivatives of Uracil Tautomers.

The influence of solvents on intramolecular interactions in 5- or 6-substituted nitro and amino derivatives of six tautomeric forms of uracil was investigated. For this purpose, the density functional theory (B97-D3/aug-cc-pVDZ) calculations were performed in ten environments (1 > ε > 109) using the polarizable continuum model (PCM) of solvation. The substituents were characterized by electronic (charge of the substituent active region, cSAR) and geometric parameters. Intramolecular interactions between non-covalently bonded atoms were investigated using the theory of atoms in molecules (AIM) and the non-covalent interaction index (NCI) method, which allowed discussion of possible interactions between the substituents and N/NH endocyclic as well as =O/−OH exocyclic groups. The nitro group was more electron-withdrawing in the 5 than in the 6 position, while the opposite effect was observed in the case of electron donation of the amino group. These properties of both groups were enhanced in polar solvents; the enhancement depended on the ortho interactions. Substitution or solvation did not change tautomeric preferences of uracil significantly. However, the formation of a strong NO∙∙∙HO intramolecular hydrogen bond in the 5-NO2 derivative stabilized the dienol tautomer from +17.9 (unsubstituted) to +5.4 kcal/mol (substituted, energy relative to the most stable diketo tautomer).

Four pairs of tautomers from rhizomes of Acorus calamus and fruits of Cornus officinalis.

A pair of new guaiane-type sesquiterpene tautomers (1) was isolated from rhizomes of Acorus calamus. Meanwhile, three pairs of known compounds, including a pair of dihydroflavone glycoside tautomers (2), a pair of 4-hydroxybenzoic acid ester glycoside tautomers (3), as well as a pair of secoiridoid glycoside tautomers (4) were isolated from fruits of Cornus officinalis. Their structures were elucidated by extensive spectroscopic and computational methods. Furthermore, the tautomeric mechanisms were discussed.

A review on synthetic strategy, molecular pharmacology of indazole derivatives, and their future perspective.

With different nitrogen-containing heterocyclic moieties, Indazoles earn one of the places among the top investigated molecules in medicinal research. Indazole, an important fused aromatic heterocyclic system containing benzene and pyrazole ring with a chemical formula of C7 H6 N2 , is also called benzopyrazole. Indazoles consist of three tautomeric forms in which 1H-tautomers (indazoles) and 2H-tautomers (isoindazoles) exist in all phases. The tautomerism in indazoles greatly influences synthesis, reactivity, physical and even the biological properties of indazoles. The thermodynamic internal energy calculation of these tautomers points view 1H-indazole as the predominant and stable form over 2H-indazole. The natural source of indazole is limited and exists in alkaloidal nature (i.e., nigellidine, nigeglanine, nigellicine, etc.) found from Nigella plants. Some of the FDA-approved drugs like Axitinib, Entrectinib, Niraparib, Benzydamine, and Granisetron are being used to treat renal cell cancer, non-small cell lung cancer (NSCLC), epithelial ovarian cancer, chronic inflammation, chemotherapy-induced nausea, vomiting, and many more uses. Besides all these advantages regarding its biological activity, the main issue about indazoles is the less abundance in plant sources, and their synthetic derivatives also often face problems with low yield. In this review article, we discuss its chemistry, tautomerism along with their effects, different schematics for the synthesis of indazole derivatives, and their different biological activities.

Solvent Influence on Absorption Spectra and Tautomeric Equilibria of Symmetric Azomethine-Functionalized Derivatives: Structural Elucidation and Computational Studies.

A new series of azomethine-functionalized compounds was synthesized from the condensation of 2-hydroxy-1,3-propanediamine and 2-thienylcarboxaldehydes in the presence of a drying agent. The derivatives were spectroscopically characterized by NMR, LC-MS, UV/Vis, IR and elemental analysis. Variable temperature 1 H-NMR (-60 to +60 °C) was performed to investigate the effect of solvent polarity; the capability of solvent to form H-bond was found to dramatically influencing the tautomerization process of the desired structures. The calculated thermochemical parameters (ΔH298 , ΔG298 and ΔS298 ) at DFT and MP2 levels of theory explained that 3 b exists in equilibrium with two tautomers. The basis of the electronic absorptions was pursued through Time-Dependent Density-Functional Theory (TD-DFT). Analysis of the structural surfaces was inspected and the molecular electrostatic potential (MEP) demonstrated that the three functionalized compounds were relatively analogous in the electronic distributions. Furthermore, the electrophilic and nucleophilic centers lying on the molecular surfaces were probably playing a key-role in stabilizing the compounds through the nonclassical C-H⋅⋅⋅π interactions and hydrogen bonding. The impact of solvent polarity on absorption spectra were investigated via solvatochromic shifts. For instance, compound 3 c displayed a gradual shift of the maximum absorption to the red area when the solvent polarity was increased, recording a 21 nm of bathochromic shift. In contrast, no significant solvent-effect on 3 a and 3 b was observed. The solvation relation was pursued between Gutmann's donicity numbers the experimental λmax ; exhibited almost positive linear performance with a minor oscillation, that ascribe to the possible weak interface between the molecules of solute and designated solvents. The bandgap energy of all products were assessed experimentally using optical absorption spectra following Tauc approach, giving -4.050 (3 a), -3.900 (3 b) and -3.210 (3 c) eV. However, the ΔE were computationally figured out from TD-DFT simulation to be -4.258 (3 a), -4.022 (3 b) and -3.390 (3 c) eV.