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.

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.

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-G_q 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-G_q 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-G_q complex and rather induces an inactive conformation in the absence of the G_q 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.

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.

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 CH_2'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.

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 x_E = 0.158 (c_E = 302.96 mg/mL) and x_E = 0.105 (c_E = 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 x_E = 0.193 (c_E = 459.69 mg/mL) and x_E = 0.145 (c_E = 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).

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).

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 C₇ H₆ N₂ , 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 ¹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 (ΔH₂₉₈, ΔG₂₉₈ and ΔS₂₉₈) 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.

The critical role of QM/MM X-ray refinement and accurate tautomer/protomer determination in structure-based drug design

Conventional protein:ligand crystallographic refinement uses stereochemistry restraints coupled with a rudimentary energy functional to ensure the correct geometry of the model of the macromolecule—along with any bound ligand(s)—within the context of the experimental, X-ray density. These methods generally lack explicit terms for electrostatics, polarization, dispersion, hydrogen bonds, and other key interactions, and instead they use pre-determined parameters (e.g. bond lengths, angles, and torsions) to drive structural refinement. In order to address this deficiency and obtain a more complete and ultimately more accurate structure, we have developed an automated approach for macromolecular refinement based on a two layer, QM/MM (ONIOM) scheme as implemented within our DivCon Discovery Suite and "plugged in" to two mainstream crystallographic packages: PHENIX and BUSTER. This implementation is able to use one or more region layer(s), which is(are) characterized using linear-scaling, semi-empirical quantum mechanics, followed by a system layer which includes the balance of the model and which is described using a molecular mechanics functional. In this work, we applied our Phenix/DivCon refinement method—coupled with our XModeScore method for experimental tautomer/protomer state determination—to the characterization of structure sets relevant to structure-based drug design (SBDD). We then use these newly refined structures to show the impact of QM/MM X-ray refined structure on our understanding of function by exploring the influence of these improved structures on protein:ligand binding affinity prediction (and we likewise show how we use post-refinement scoring outliers to inform subsequent X-ray crystallographic efforts). Through this endeavor, we demonstrate a computational chemistry ↔ structural biology (X-ray crystallography) "feedback loop" which has utility in industrial and academic pharmaceutical research as well as other allied fields.

Factors Governing the Chemical Stability and NMR Parameters of Uracil Tautomers and Its 5-Halogen Derivatives

We report on the density functional theory (DFT) modelling of structural, energetic and NMR parameters of uracil and its derivatives (5-halogenouracil (5XU), X = F, Cl, Br and I) in vacuum and in water using the polarizable continuum model (PCM) and the solvent model density (SMD) approach. On the basis of the obtained results, we conclude that the intramolecular electrostatic interactions are the main factors governing the stability of the six tautomeric forms of uracil and 5XU. Two indices of aromaticity, the harmonic oscillator model of aromaticity (HOMA), satisfying the geometric criterion, and the nuclear independent chemical shift (NICS), were applied to evaluate the aromaticity of uracil and its derivatives in the gas phase and water. The values of these parameters showed that the most stable tautomer is the least aromatic. A good performance of newly designed xOPBE density functional in combination with both large aug-cc-pVQZ and small STO(1M)-3G basis sets for predicting chemical shifts of uracil and 5-fluorouracil in vacuum and water was observed. As a practical alternative for calculating the chemical shifts of challenging heterocyclic compounds, we also propose B3LYP calculations with small STO(1M)-3G basis set. The indirect spin-spin coupling constants predicted by B3LYP/aug-cc-pVQZ(mixed) method reproduce the experimental data for uracil and 5-fluorouracil well.

A New Bias Site for Epigenetic Modifications: How Non-Canonical GC Base Pairs Favor Mechanochemical Cleavage of DNA

Properties of non-canonical GC base pairs and their relations with mechanochemical cleavage of DNA are analyzed. A hypothesis of the involvement of the transient GC wobble base pairs both in the mechanisms of the mechanochemical cleavage of DNA and epigenetic mechanisms involving of 5-methylcytosine, is proposed. The hypothesis explains the increase in the frequency of the breaks of the sugar-phosphate backbone of DNA after cytosines, the asymmetric character of these breaks, and an increase in break frequency in CpG after cytosine methylation. As an alternative hypothesis, probable implication of GC⁺ Hoogsteen base pairs is considered, which now exemplify the best-studied non-canonical GC base pairs in the DNA double helix. Also see the video abstract here https://youtu.be/EUunVWL0ptw.

Substrate-Induced Formation of Ribosomal Decoding Center for Accurate and Rapid Genetic Code Translation

Accurate translation of genetic information is crucial for synthesis of functional proteins in all organisms. We use recent experimental data to discuss how induced fit affects accuracy of initial codon selection on the ribosome by aminoacyl transfer RNA in ternary complex ( T₃) with elongation factor Tu (EF-Tu) and guanosine-5'-triphosphate (GTP). We define actual accuracy ([Formula: see text]) of a particular protein synthesis system as its current accuracy and the effective selectivity ([Formula: see text]) as [Formula: see text] in the limit of zero ribosomal binding affinity for T₃. Intrinsic selectivity ([Formula: see text]), defined as the upper thermodynamic limit of [Formula: see text], is determined by the free energy difference between near-cognate and cognate T₃ in the pre-GTP hydrolysis state on the ribosome. [Formula: see text] is much larger than [Formula: see text], suggesting the possibility of a considerable increase in [Formula: see text] and [Formula: see text] at negligible kinetic cost. Induced fit increases [Formula: see text] and [Formula: see text] without affecting [Formula: see text], and aminoglycoside antibiotics reduce [Formula: see text] and [Formula: see text] at unaltered [Formula: see text].

Influence of Hydroxyl Functional Group on the Structure and Stability of Xanthone: A Computational Approach

The present work addresses computational research focused on the energetic and structural properties of four isomers monohydroxyxanthone, using the G3(MP2)//B3LYP method, in order to evaluate the influence of the hydroxyl (-OH moiety) functional group on the xanthone molecule. The combination of these computational results with previous experimental data of these compounds enabled the determination of their enthalpies, entropies and Gibbs energies of formation, in the gaseous phase, and consequently to infer about the relative thermodynamic stability of the four isomers. Other issues were also addressed for the hydroxyxanthone isomers, namely the conformational and the tautomeric equilibrium analysis of the optimized molecular structures, the frontier orbitals, and the electrostatic potential energy maps. Complementarily, an energetic study of the intramolecular O - H ⋯ O hydrogen bond for 1-hydroxanthone was also performed.

A comprehensive analysis of the computed tautomer fractions of the imidazole ring of histidines in Loligo vulgaris

A recently introduced electrostatic-based method to determine the pKa values of ionizable residues and fractions of ionized and tautomeric forms of histidine (His) and acid residues in proteins, at a given fixed pH, is applied here to the analysis of a His-rich protein, namely Loligo vulgaris (pdb id 1E1A), a 314-residue all-β protein. The average tautomeric fractions for the imidazole ring of each of the six histidines in the sequence were computed using an approach that includes, but is not limited to, molecular dynamic simulations coupled with calculations of the ionization states for all 94 ionizable residues of protein 1E1A in water at pH 6.5 and 300 K. The electrostatic-calculated tautomeric fractions of the imidazole ring of His were compared with predictions obtained from an existent NMR-based methodology. Our results indicate that: (i) the averaged electrostatic-based tautomeric predictions for the imidazole ring of all histidines of Loligo vulgaris are dominated by the Nε2-H rather than the Nδ1-H form, although such preferences from the NMR-based methodology are not so well defined; (ii) the computed average absolute difference between the electrostatic- and the NMR-based tautomeric predictions among all six histidines vary among 0% to 17%; (iii) for the His showing the largest fraction of the neutral form (81%), the absolute difference between the NMR- and electrostatic-based computed tautomeric predictions is only 3%; and (iv) the tautomeric predictions for the imidazole ring of His computed with the NMR-based methodology are stable within a certain, well-defined, range of variations of a tautomer-related parameter.

Multinuclear NMR Measurements and DFT Calculations for Capecitabine Tautomeric Form Assignment in a Solution

The molecular structure of capecitabine (a widely applied prodrug of 5-fluorouracil) was studied by multinuclear NMR measurements and DFT quantum mechanical calculations. One or two tautomeric forms in a solution were detected depending on the solvent used. In the organic solvents, a mixture of two forms of capecitabine was observed: carbamate and imine tautomers. In the aqueous solution, only the carbamate form was found. The methylation of capecitabine yields mainly two products in different proportions: N³-methylcapecitabine and N⁷-methylcapecitabine. The protonation of capecitabine in organic solvents with perchloric acid occurs at the N3 nitrogen atom. DFT calculations strongly support the results coming from the analysis of the NMR spectra.

A recent intermezzo at the Ribosome Club

Two sets of ribosome structures have recently led to two different interpretations of what limits the accuracy of codon translation by transfer RNAs. In this review, inspired by this intermezzo at the Ribosome Club, we briefly discuss accuracy amplification by energy driven proofreading and its implementation in genetic code translation. We further discuss general ways by which the monitoring bases of 16S rRNA may enhance the ultimate accuracy (d-values) and how the codon translation accuracy is reduced by the actions of Mg²⁺ ions and the presence of error inducing aminoglycoside antibiotics. We demonstrate that complete freezing-in of cognate-like tautomeric states of ribosome-bound nucleotide bases in transfer RNA or messenger RNA is not compatible with recent experiments on initial codon selection by transfer RNA in ternary complex with elongation factor Tu and GTP. From these considerations, we suggest that the sets of 30S subunit structures from the Ramakrishnan group and 70S structures from the Yusupov/Yusupova group may, after all, reflect two sides of the same coin and how the structurally based intermezzo at the Ribosome Club may be resolved simply by taking the dynamic aspects of ribosome function into account.

This article is part of the themed issue ‘Perspectives on the ribosome’.

Stabilization and Structure of the cis Tautomer of a Free-Base Porphyrin

Single-crystal X-ray analysis of the β-heptakis(trifluoromethyl)-meso-tetrakis(p-fluorophenyl)porphyrin, H₂ [(CF₃ )₇ TpFPP], has revealed the first example of a stable cis tautomer of a free-base porphyrin, the long-postulated intermediate of porphyrin tautomerism. The stability of the unique molecule appears to reflect a dual origin: a strongly saddled porphyrin skeleton, which alleviates electrostatic repulsion between the two NH protons, and two polarization-enhanced, transannular N-H⋅⋅⋅O-H⋅⋅⋅N hydrogen bond chains, each involving a molecule of water. DFT calculations suggest that the observed tautomer has a lower energy than the alternative, doubly hydrated trans tautomer by some 8.3 kcal mol^-1 . A fascinating prospect thus exists that H₂ [(CF₃ )₇ TpFPP]⋅2 H₂ O and cognate structures may act as supramolecular synthons, which, given their chirality, may even be amenable to resolution into optically pure enantiomers.

Getting the chemistry right: protonation, tautomers and the importance of H atoms in biological chemistry

There are more H atoms than any other type of atom in an X-ray crystal structure of a protein-ligand complex, but as H atoms only have one electron they diffract X-rays weakly and are `hard to see'. The positions of many H atoms can be inferred by our chemical knowledge, and such H atoms can be added with confidence in `riding positions'. For some chemical groups, however, there is more ambiguity over the possible hydrogen placements, for example hydroxyls and groups that can exist in multiple protonation states or tautomeric forms. This ambiguity is far from rare, since about 25% of drugs have more than one tautomeric form. This paper focuses on the most common, `prototropic', tautomers, which are isomers that readily interconvert by the exchange of an H atom accompanied by the switch of a single and an adjacent double bond. Hydrogen-exchange rates and different protonation states of compounds (e.g. buffers) are also briefly discussed. The difference in heavy (non-H) atom positions between two tautomers can be small, and careful refinement of all possible tautomers may single out the likely bound ligand tautomer. Experimental methods to determine H-atom positions, such as neutron crystallography, are often technically challenging. Therefore, chemical knowledge and computational approaches are frequently used in conjugation with experimental data to deduce the bound tautomer state. Proton movement is a key feature of many enzymatic reactions, so understanding the orchestration of hydrogen/proton motion is of critical importance to biological chemistry. For example, structural studies have suggested that, just as a chemist may use heat, some enzymes use directional movement to protonate specific O atoms on phosphates to catalyse phosphotransferase reactions. To inhibit `wriggly' enzymes that use movement to effect catalysis, it may be advantageous to have inhibitors that can maintain favourable contacts by adopting different tautomers as the enzyme `wriggles'.

XModeScore: a novel method for accurate protonation/tautomer-state determination using quantum-mechanically driven macromolecular X-ray crystallographic refinement

XModeScore determines the correct protomeric/tautomeric state or mode of active-site residues along with any bound ligand(s) using quantum-mechanics-based X-ray refinement followed by post-refinement scoring based on a combination of energetic strain (or ligand strain) and rigorous difference electron-density analysis.

Gaining an understanding of the protein–ligand complex structure along with the proper protonation and explicit solvent effects can be important in obtaining meaningful results in structure-guided drug discovery and structure-based drug discovery. Unfortunately, protonation and tautomerism are difficult to establish with conventional methods because of difficulties in the experimental detection of H atoms owing to the well known limitations of X-ray crystallography. In the present work, it is demonstrated that semiempirical, quantum-mechanics-based macromolecular crystallographic refinement is sensitive to the choice of a protonation-state/tautomer form of ligands and residues, and can therefore be used to explore potential states. A novel scoring method, called XModeScore, is described which enumerates the possible protomeric/tautomeric modes, refines each mode against X-ray diffraction data with the semiempirical quantum-mechanics (PM6) Hamiltonian and scores each mode using a combination of energetic strain (or ligand strain) and rigorous statistical analysis of the difference electron-density distribution. It is shown that using XModeScore it is possible to consistently distinguish the correct bound protomeric/tautomeric modes based on routine X-ray data, even at lower resolutions of around 3 Å. These X-ray results are compared with the results obtained from much more expensive and laborious neutron diffraction studies for three different examples: tautomerism in the acetazolamide ligand of human carbonic anhydrase II (PDB entries 3hs4 and 4k0s), tautomerism in the 8HX ligand of urate oxidase (PDB entries 4n9s and 4n9m) and the protonation states of the catalytic aspartic acid found within the active site of an aspartic protease (PDB entry 2jjj). In each case, XModeScore applied to the X-ray diffraction data is able to determine the correct protonation state as defined by the neutron diffraction data. The impact of QM-based refinement versus conventional refinement on XModeScore is also discussed.

Geometric consequences of electron delocalization for adenine tautomers in aqueous solution

Geometric consequences of electron delocalization were studied for all possible adenine tautomers in aqueous solution by means of ab initio methods {PCM(water)//DFT(B3LYP)/6-311+G(d,p)} and compared to those in the gas phase {DFT(B3LYP)/6-311+G(d,p)}. To measure the consequences of any type of resonance conjugation (π-π, n-π, and σ-π), the geometry-based harmonic oscillator model of electron delocalization (HOMED) index, recently extended to the isolated (DFT) and hydrated (PCM//DFT) molecules, was applied to the molecular fragments (imidazole, pyrimidine, 4-aminopyrimidine, and purine) and also to the whole tautomeric system. For individual tautomers, the resonance conjugations and consequently the bond lengths strongly depend on the position of the labile protons. The HOMED indices are larger for tautomers (or their fragments) possessing the labile proton(s) at the N rather than C atom. Solvent interactions with adenine tautomers slightly increase the resonance conjugations. Consequently, they slightly shorten the single bonds and lengthen the double bonds. When going from the gas phase to water solution, the HOMED indices increase (by less than 0.15 units). There is a good relation between the HOMED indices estimated in water solution and those in the gas phase for the neutral and ionized forms of adenine. Subtle effects, being a consequence of intramolecular interactions between the neighboring groups, are so strongly reduced by solvent that the relation between the HOMED indices and the relative energies for the neutral adenine tautomers seems to be better in water solution than in the gas phase.

GIAO/DFT studies on 1,2,4-triazole-5-thiones and their propargyl derivatives

Density functional theory (DFT)/Becke-Lee-Yang-Parr (B3LYP) and gauge-including atomic orbital (GIAO) calculations were performed on a number of 1,2,4-triazole derivatives, and the optimized structural parameters were employed to ascertain the nature of their predominant tautomers. (13)C and (15)N NMR chemical shifts of 3-substituted 1,2,4-triazole-5-thiones and their propargylated derivatives were calculated via GIAO/DFT approach at the B3LYP level of theory with geometry optimization using a 6-311++G** basis set. A good agreement between theoretical and experimental (13)C and (15)N NMR chemical shifts could be found for the systems investigated. The data generated were useful in predicting (15)N chemical shifts of all the nitrogen atoms of the triazole ring, some of which could not be obtained in solution state (15)N HMBC/HSQC NMR measurements. The energy profile computed for the dipropargylated derivatives was found to follow the product distribution profile of regioisomers formed during propargylation of 1,2,4-triazole thiones.

SERS, XPS, and DFT Study of Adenine Adsorption on Silver and Gold Surfaces

The adsorption of adenine on silver and gold surfaces has been investigated combining density functional theory calculations with surface-enhanced Raman scattering and angle-resolved X-ray photoelectron spectroscopy measurements, obtaining useful insight into the orientation and interaction of the nucleobase with the metal surfaces.

Bound anionic states of adenine

Anionic states of nucleic acid bases are involved in DNA damage by low-energy electrons and in charge transfer through DNA. Previous gas phase studies of free, unsolvated nucleic acid base parent anions probed only dipole-bound states, which are not present in condensed phase environments, but did not observe valence anionic states, which for purine bases are thought to be adiabatically unbound. Contrary to this expectation, we have demonstrated that some thus far ignored tautomers of adenine, which result from enamine-imine transformations, support valence anionic states with electron vertical detachment energies as large as 2.2 eV, and at least one of these anionic tautomers is adiabatically bound. Moreover, we predict that the new anionic tautomers should also dominate in solutions and should be characterized by larger values of electron vertical detachment energy than the canonical valence anion. All of the newfound anionic tautomers might be formed in the course of dissociative electron attachment followed by a hydrogen atom attachment to a carbon atom, and they might affect the structure and properties of DNA and RNA exposed to low-energy electrons. The new valence states observed here, unlike the dipole-bound state, could exist in condensed phases and might be relevant to radiobiological damage. The discovery of these valence anionic states of adenine was facilitated by the development of (i) an experimental method for preparing parent anions of nucleic acid bases for photoelectron experiments, and (ii) a combinatorial/quantum chemical approach for identification of the most stable tautomers of organic molecules.

The guanine cation radical: investigation of deprotonation states by ESR and DFT

This work reports ESR studies that identify the favored site of deprotonation of the guanine cation radical (G*+) in an aqueous medium at 77 K. Using ESR and UV-visible spectroscopy, one-electron oxidized guanine is investigated in frozen aqueous D2O solutions of 2'-deoxyguanosine (dGuo) at low temperatures at various pHs at which the guanine cation radical, G*+ (pH 3-5), singly deprotonated species, G(-H)* (pH 7-9), and doubly deprotonated species, G(-2H)*- (pH > 11), are found. C-8-deuteration of dGuo to give 8-D-dGuo removes the major proton hyperfine coupling at C-8. This isolates the anisotropic nitrogen couplings for each of the three species and aids our analyses. These anisotropic nitrogen couplings were assigned to specific nitrogen sites by use of 15N-substituted derivatives at N1, N2, and N3 atoms in dGuo. Both ESR and UV-visible spectra are reported for each of the species: G*+, G(-H)*, and G(-2H)*-. The experimental anisotropic ESR hyperfine couplings are compared to those obtained from DFT calculations for the various tautomers of G(-H)*. Using the B3LYP/6-31G(d) method, the geometries and energies of G*+ and its singly deprotonated state in its two tautomeric forms, G(N1-H)* and G(N2-H)*, were investigated. In a nonhydrated state, G(N2-H)* is found to be more stable than G(N1-H)*, but on hydration with seven water molecules G(N1-H)* is found to be more stable than G(N2-H)*. The theoretically calculated hyperfine coupling constants (HFCCs) of G*+, G(N1-H)*, and G(-2H)*- match the experimentally observed HFCCs best on hydration with seven or more waters. For G(-2H)*-, the hyperfine coupling constant (HFCC) at the exocyclic nitrogen atom (N2) is especially sensitive to the number of hydrating water molecules; good agreement with experiment is not obtained until nine or 10 waters of hydration are included.

Diastereoselective Synthesis of a Simplified Core of Rishirilide B

A route enabling the synthesis of the stereo-triad of rishirilide B (1) from 2-hydroxy-3-methylnaphthalene-1,4-dione, is reported. Key transformations include the regioselective 1,2-Grignard addition to a tautomeric mixture of o- and p-quinones, regioselective carbamoylation of a tautomeric mixture, and a synopsis of the methods explored to convert various terminal vinyl ethers into the corresponding carboxylic acid by cleavage.