Low-barrier hydrogen bonds and low fractionation factor bases in enzymatic reactions

Structural Insight on Supramolecular Polyion Salts: Inositol Hexaphosphate Enclosed in Cationic Macrocyclic Clusters

Supramolecular macrocyclic forces have been used to trap phytate, myo-inositol-1,2,3,4,5,6-hexakisphosphate, a key bioanion with multiple roles in metabolic processes. Due to the complex chemistry of six multivalent phosphates surrounding the small, cyclic inositol framework, crystallographic information of simple phytate salts has been elusive. This report represents a combined crystallographic, theoretical, and solution binding investigation of a supramolecular macrocyclic complex of phytate. Together, the results provide significant insight to phytate's intramolecular and intermolecular interactions at the microenvironment level. The macrocycle-phytate aggregates consist of phytate anionic pairs, each partly sandwiched by two 24-membered, amide/amine-based cationic macrocycles. The phytate ion pairs hold the tetrameric macrocyclic array together by six strong intermolecular hydrogen bonds. Both phytates crystallize in 1a5e phosphate conformations (one axial (P2) and five equatorial phosphates). Solution NMR binding studies in 1 : 1 DMSO-d6  : D2 O indicate 2 : 1 macrocycle:phytate associations, suggesting that the sandwich-like nature of the complex holds together in solution. DFT studies indicate the likely occurrence of dynamic intramolecular interchange of phosphate protons, as well as important roles for the axial (P2) phosphate in both intramolecular and intermolecular hydrogen bonding interactions.

Chemical features and machine learning assisted predictions of protein-ligand short hydrogen bonds

There are continuous efforts to elucidate the structure and biological functions of short hydrogen bonds (SHBs), whose donor and acceptor heteroatoms reside more than 0.3 Å closer than the sum of their van der Waals radii. In this work, we evaluate 1070 atomic-resolution protein structures and characterize the common chemical features of SHBs formed between the side chains of amino acids and small molecule ligands. We then develop a machine learning assisted prediction of protein-ligand SHBs (MAPSHB-Ligand) model and reveal that the types of amino acids and ligand functional groups as well as the sequence of neighboring residues are essential factors that determine the class of protein-ligand hydrogen bonds. The MAPSHB-Ligand model and its implementation on our web server enable the effective identification of protein-ligand SHBs in proteins, which will facilitate the design of biomolecules and ligands that exploit these close contacts for enhanced functions.

Brønsted Acids Promote Olefin Oxidations by Bioinspired Nonheme CoIII(PhIO)(OH) Complexes: A Role for Low-Barrier Hydrogen Bonds

Introduction of Brønsted acids into biomimetic nonheme reactions promotes the oxidative ability of metal-oxygen complexes significantly. However, the molecular machinery of the promoted effects is missing. Herein, a comprehensive investigation of styrene oxidation by a cobalt(III)-iodosylbenzene complex, [(TQA)Co^III(OIPh)(OH)]²⁺ (1, TQA = tris(2-quinolylmethyl)amine), in the presence and absence of triflic acid (HOTf) was performed using density functional theory calculations. Results revealed for the first time that there is a low-barrier hydrogen bond (LBHB) between HOTf and the hydroxyl ligand of 1, which forms two valence-resonance structures [(TQA)Co^III(OIPh)(HO^---HOTf)]²⁺ (1_LBHB) and [(TQA)Co^III(OIPh)(H₂O--OTf^-)]²⁺ (1'_LBHB). Due to the oxo-wall, these complexes (1_LBHB and 1'_LBHB) cannot convert to high-valent cobalt-oxyl species. Instead, styrene oxidation by these oxidants (1_LBHB and 1'_LBHB) shows novel spin-state selectivity, i.e., on the ground closed-shell singlet state, styrene is oxidized to an epoxide, whereas on the excited triplet and quintet states, an aldehyde product, phenylacetaldehyde, is formed. The preferred pathway is styrene oxidation by 1'_LBHB, which is initiated by a rate-limiting bond-formation-coupled electron transfer process with an energy barrier of 12.2 kcal mol^-1. The nascent PhIO-styrene-radical-cation intermediate undergoes an intramolecular rearrangement to produce an aldehyde. The halogen bond between the OH^-/H₂O ligand and the iodine of PhIO modulates the activity of the cobalt-iodosylarene complexes 1_LBHB and 1'_LBHB. These new mechanistic findings enrich our knowledge of nonheme chemistry and hypervalent iodine chemistry and will play a positive role in the rational design of new catalysts.

Common Mechanism of Activated Catalysis in P-loop Fold Nucleoside Triphosphatases-United in Diversity

To clarify the obscure hydrolysis mechanism of ubiquitous P-loop-fold nucleoside triphosphatases (Walker NTPases), we analysed the structures of 3136 catalytic sites with bound Mg-NTP complexes or their analogues. Our results are presented in two articles; here, in the second of them, we elucidated whether the Walker A and Walker B sequence motifs-common to all P-loop NTPases-could be directly involved in catalysis. We found that the hydrogen bonds (H-bonds) between the strictly conserved, Mg-coordinating Ser/Thr of the Walker A motif ([Ser/Thr]WA) and aspartate of the Walker B motif (AspWB) are particularly short (even as short as 2.4 ångströms) in the structures with bound transition state (TS) analogues. Given that a short H-bond implies parity in the pKa values of the H-bond partners, we suggest that, in response to the interactions of a P-loop NTPase with its cognate activating partner, a proton relocates from [Ser/Thr]WA to AspWB. The resulting anionic [Ser/Thr]WA alkoxide withdraws a proton from the catalytic water molecule, and the nascent hydroxyl attacks the gamma phosphate of NTP. When the gamma-phosphate breaks away, the trapped proton at AspWB passes by the Grotthuss relay via [Ser/Thr]WA to beta-phosphate and compensates for its developing negative charge that is thought to be responsible for the activation barrier of hydrolysis.

Small Molecules Promote Selective Denaturation and Degradation of Tubulin Heterodimers through a Low-Barrier Hydrogen Bond

Here, we report a novel mechanism to selectively degrade target proteins. 3-(3-Phenoxybenzyl)amino-β-carboline (PAC), a tubulin inhibitor, promotes selective degradation of αβ-tubulin heterodimers. Biochemical studies have revealed that PAC specifically denatures tubulin, making it prone to aggregation that predisposes it to ubiquitinylation and then degradation. The degradation is mediated by a single hydrogen bond formed between the pyridine nitrogen of PAC and βGlu198, which is identified as a low-barrier hydrogen bond (LBHB). In contrast, another two tubulin inhibitors that only form normal hydrogen bonds with βGlu198 exhibit no degradation effect. Thus, the LBHB accounts for the degradation. We then screened for compounds capable of forming an LBHB with βGlu198 and demonstrated that BML284, a Wnt signaling activator, also promotes tubulin heterodimer degradation through the LBHB. Our study provided a unique example of LBHB function and identified a novel approach to obtain tubulin degraders.

Neutron crystallography reveals mechanisms used by Pseudomonas aeruginosa for host-cell binding

The opportunistic pathogen Pseudomonas aeruginosa, a major cause of nosocomial infections, uses carbohydrate-binding proteins (lectins) as part of its binding to host cells. The fucose-binding lectin, LecB, displays a unique carbohydrate-binding site that incorporates two closely located calcium ions bridging between the ligand and protein, providing specificity and unusually high affinity. Here, we investigate the mechanisms involved in binding based on neutron crystallography studies of a fully deuterated LecB/fucose/calcium complex. The neutron structure, which includes the positions of all the hydrogen atoms, reveals that the high affinity of binding may be related to the occurrence of a low-barrier hydrogen bond induced by the proximity of the two calcium ions, the presence of coordination rings between the sugar, calcium and LecB, and the dynamic behaviour of bridging water molecules at room temperature. These key structural details may assist in the design of anti-adhesive compounds to combat multi-resistance bacterial infections.

Pseudomonas aeruginosa employs lectins to bind to its host cells, and is known to be the major cause of lung infections. Lectin B (LecB) from Pseudomonas aeruginosa binds specifically to galactose and fucose and is important for pathogenicity, adhesion and biofilm formation. In this work, the neutron crystal structure (1.9 Å) of the deuterated LecB/Ca/fucose complex is reported. The structure, in combination with perdeuteration of the ligand and the receptor allowed the observation of hydrogen atoms, protonation states and hydrogen bonds involved in the interaction between pathogenic bacteria and host cells. Thus the study provides structural insights into the mechanism of high affinity binding of LecB to its targets.

Beyond the Plateau: pL Dependence of Proton Inventories as a Tool for Studying Ribozyme and Ribonuclease Catalysis

Acid/base catalysis is an important catalytic strategy used by ribonucleases and ribozymes; however, understanding the number and identity of functional groups involved in proton transfer remains challenging. The proton inventory (PI) technique analyzes the dependence of the enzyme reaction rate on the ratio of D₂O to H₂O and can provide information about the number of exchangeable sites that produce isotope effects and their magnitude. The Gross-Butler (GB) equation is used to evaluate H/D fractionation factors from PI data typically collected under conditions (i.e., a "plateau" in the pH-rate profile) assuming minimal change in active site residue ionization. However, restricting PI analysis to these conditions is problematic for many ribonucleases, ribozymes, and their variants due to ambiguity in the roles of active site residues, the lack of a plateau within the accessible pL range, or cooperative interactions between active site functional groups undergoing ionization. Here, we extend the integration of species distributions for alternative enzyme states in noncooperative models of acid/base catalysis into the GB equation, first used by Bevilacqua and colleagues for the HDV ribozyme, to develop a general population-weighted GB equation that allows simulation and global fitting of the three-dimensional relationship of the D₂O ratio (n) versus pL versus k_n/k₀. Simulations using the GPW-GB equation of PI results for RNase A, HDVrz, and VSrz illustrate that data obtained at multiple selected pL values across the pL-rate profile can assist in the planning and interpreting of solvent isotope effect experiments to distinguish alternative mechanistic models.

Mechanism of ArcLight derived GEVIs involves electrostatic interactions that can affect proton wires

The genetically encoded voltage indicators ArcLight and its derivatives mediate voltage-dependent optical signals by intermolecular, electrostatic interactions between neighboring fluorescent proteins (FPs). A random mutagenesis event placed a negative charge on the exterior of the FP, resulting in a greater than 10-fold improvement of the voltage-dependent optical signal. Repositioning this negative charge on the exterior of the FP reversed the polarity of voltage-dependent optical signals, suggesting the presence of “hot spots” capable of interacting with the negative charge on a neighboring FP, thereby changing the fluorescent output. To explore the potential effect on the chromophore state, voltage-clamp fluorometry was performed with alternating excitation at 390 nm followed by excitation at 470 nm, resulting in several mutants exhibiting voltage-dependent, ratiometric optical signals of opposing polarities. However, the kinetics, voltage ranges, and optimal FP fusion sites were different depending on the wavelength of excitation. These results suggest that the FP has external, electrostatic pathways capable of quenching fluorescence that are wavelength specific. One mutation to the FP (E222H) showed a voltage-dependent increase in fluorescence when excited at 390 nm, indicating the ability to affect the proton wire from the protonated chromophore to the H222 position. ArcLight-derived sensors may therefore offer a novel way to map how conditions external to the β-can structure can affect the fluorescence of the chromophore and transiently affect those pathways via conformational changes mediated by manipulating membrane potential.

On the Case of the Misplaced Hydrogens

Few other elements play a more central role in biology than hydrogen. The interactions, bonding and movement of hydrogen atoms are central to biological catalysis, structure and function. Yet owing to the elusive nature of a single hydrogen atom few experimental and computational techniques can precisely determine its location. This is exemplified in short hydrogen bonds (SHBs) where the location of the hydrogen atom is indicative of the underlying strength of the bonds, which can vary from 1-5 kcal/mol in canonical hydrogen bonds, to an almost covalent nature in single-well hydrogen bonds. Owing to the often-times inferred position of hydrogen, the role of SHBs in biology has remained highly contested and debated. This has also led to discrepancies in computational, biochemical and structural studies of proteins thought to use SHBs in performing chemistry and stabilizing interactions. Herein, we discuss in detail two distinct examples, namely the conserved catalytic triad and the photoreceptor, photoactive yellow protein, where studies of these SHB-containing systems have permitted contextualization of the role these unique hydrogen bonds play in biology.

Low barrier hydrogen bonds in protein structure and function

Low-barrier hydrogen bonds (LBHBs) are a special type of short hydrogen bond (HB) that is characterized by the equal sharing of a hydrogen atom. The existence and catalytic role of LBHBs in proteins has been intensely contested. Advancements in X-ray and neutron diffraction methods has revealed delocalized hydrogen atoms involved in potential LBHBs in a number of proteins, while also demonstrating that short HBs are not necessarily LBHBs. More importantly, a series of experiments on ketosteroid isomerase (KSI) have suggested that LBHBs are significantly stronger than standard HBs in the protein microenvironment in terms of enthalpy, but not free energy. The discrepancy between the enthalpy and free energy of LBHBs offers clues to the challenges, and potential solutions, of the LBHB debate, where the unique strength of LBHBs plays a special role in the kinetic processes of enzyme function and structure, together with other molecular forces in a pre-organized environment.

Intramolecular Hydrogen Bond Energy and Its Decomposition—O–H∙∙∙O Interactions

The method to calculate the energy of intramolecular hydrogen bond is proposed and tested for a sample of malonaldehyde and its fluorine derivatives; the corresponding calculations were performed at the ωB97XD/aug-cc-pVTZ level. This method based on relationships found for related intermolecular hydrogen bonds is compared with other approaches which may be applied to estimate the intramolecular hydrogen bond energy. Particularly, methods based on the comparison of the system that contains the intramolecular hydrogen bond compared with corresponding conformations where such interaction does not occur are discussed. The function-based energy decomposition analysis, FB-EDA, of the intramolecular hydrogen bonds is also proposed here.

Inverse Solvent Isotope Effects in Enzyme-Catalyzed Reactions

Solvent isotope effects have long been used as a mechanistic tool for determining enzyme mechanisms. Most commonly, macroscopic rate constants such as k_cat and k_cat/K_m are found to decrease when the reaction is performed in D₂O for a variety of reasons including the transfer of protons. Under certain circumstances, these constants are found to increase, in what is termed an inverse solvent kinetic isotope effect (SKIE), which can be a diagnostic mechanistic feature. Generally, these phenomena can be attributed to an inverse solvent equilibrium isotope effect on a rapid equilibrium preceding the rate-limiting step(s). This review surveys inverse SKIEs in enzyme-catalyzed reactions by assessing their underlying origins in common mechanistic themes. Case studies for each category are presented, and the mechanistic implications are put into context. It is hoped that readers may find the illustrative examples valuable in planning and interpreting solvent isotope effect experiments.

17O NMR studies of organic and biological molecules in aqueous solution and in the solid state

This review describes the latest developments in the field of ¹⁷O NMR spectroscopy of organic and biological molecules both in aqueous solution and in the solid state. In the first part of the review, a general theoretical description of the nuclear quadrupole relaxation process in isotropic liquids is presented at a mathematical level suitable for non-specialists. In addition to the first-order quadrupole interaction, the theory also includes additional relaxation mechanisms such as the second-order quadrupole interaction and its cross correlation with shielding anisotropy. This complete theoretical treatment allows one to assess the transverse relaxation rate (thus the line width) of NMR signals from half-integer quadrupolar nuclei in solution over the entire range of motion. On the basis of this theoretical framework, we discuss general features of quadrupole-central-transition (QCT) NMR, which is a particularly powerful method of studying biomolecules in the slow motion regime. Then we review recent advances in ¹⁷O QCT NMR studies of biological macromolecules in aqueous solution. The second part of the review is concerned with solid-state ¹⁷O NMR studies of organic and biological molecules. As a sequel to the previous review on the same subject [G. Wu, Prog. Nucl. Magn. Reson. Spectrosc. 52 (2008) 118-169], the current review provides a complete coverage of the literature published since 2008 in this area.

Substrate binding mode and catalytic mechanism of human heparan sulfate d-glucuronyl C5 epimerase

Heparan sulfate (HS) is a linear, complex polysaccharide that modulates the biological activities of proteins through binding sites made by a series of Golgi-localized enzymes. Of these, glucuronyl C5-epimerase (Glce) catalyzes C5-epimerization of the HS component, d-glucuronic acid (GlcA), into l-iduronic acid (IdoA), which provides internal flexibility to the polymer and forges protein-binding sites to ensure polymer function. Here we report crystal structures of human Glce in the unbound state and of an inactive mutant, as assessed by real-time NMR spectroscopy, bound with a (GlcA-GlcNS)_n substrate or a (IdoA-GlcNS)_n product. Deep infiltration of the oligosaccharides into the active site cleft imposes a sharp kink within the central GlcNS-GlcA/IdoA-GlcNS trisaccharide motif. An extensive network of specific interactions illustrates the absolute requirement of N-sulfate groups vicinal to the epimerization site for substrate binding. At the epimerization site, the GlcA/IdoA rings are highly constrained in two closely related boat conformations, highlighting ring-puckering signatures during catalysis. The structure-based mechanism involves the two invariant acid/base residues, Glu499 and Tyr578, poised on each side of the target uronic acid residue, thus allowing reversible abstraction and readdition of a proton at the C5 position through a neutral enol intermediate, reminiscent of mandelate racemase. These structures also shed light on a convergent mechanism of action between HS epimerases and lyases and provide molecular frameworks for the chemoenzymatic synthesis of heparin or HS analogs.

Charge-assisted hydrogen bonding as a cohesive force in soil organic matter: water solubility enhancement by addition of simple carboxylic acids

Weak bonds between molecular segments and between separate molecules of natural organic matter (OM) govern OM solubility, adsorption, supramolecular association in solution, and complexation with metal ions and oxides. We tested the hypothesis that especially strong hydrogen bonds, known as (negative) charge-assisted hydrogen bonds, (-)CAHB, contribute significantly to OM cohesion and play a role in the water solubility of solid-phase OM. The (-)CAHB, exemplified by structures such as (-CO2HO2C-)- and (-CO2HO-)-, may form between weak acids with similar proton affinity, and is shorter, more covalent, and much stronger than ordinary hydrogen bonds. Using a high-organic reference soil, we show that (-)CAHBs within the solid OM phase (intra-OM) are disrupted by solutions of aliphatic and aromatic acids, resulting in enhanced solubility of OM. The aromatic acids included naturally occurring plant exudate compounds. At constant pH and ionic strength, OM solubility increased with added organic acid concentration and molecular weight. Polar compounds incapable of forming (-)CAHBs, such as alkanols, acetonitrile, and dimethyl sulfoxide, were ineffective. Solubility enhancement showed behavior consistent with (-)CAHB theory and published observations-namely, (i) that formate is more effective than acetate due to its tendency to form stronger (-)CAHBs; (ii) that solubility enhancement peaks at pH ∼5-6, where the product of the concentration of the interacting carboxylate ions reaches a maximum; and (iii) that elution of acetate or formate through soil columns releases hydroxide ion, consistent with formation of (-)CAHBs between added acid and free weak acid groups on the solid OM. The results support the hypothesis that the (-)CAHB contributes to the cohesion of OM in the solid state.

Spectral Signatures of Proton-Transfer Dynamics at the Cusp of Low-Barrier Hydrogen Bonding

Despite their importance in diverse chemical and biochemical processes, low-barrier hydrogen bonds remain elusive targets to classify and interpret spectroscopically. Here the correlated nature of hydrogen bonding and proton transfer in the low-barrier regime has been probed for the ground and excited electronic states of 6-hydroxy-2-formylfulvene by acquiring jet-cooled fluorescence spectra of the parent and monodeuterated isotopologs. While excited-state profiles reveal regular vibronic patterns devoid of obvious dynamical signatures, their ground-state counterparts display a radically altered energy landscape characterized by spectral bifurcations comparable in magnitude to typical vibrational spacings (>100 cm^-1). Quantitative analyses yield unusual deuterium kinetic isotope effects that straddle limiting values attributed to above-barrier vibration and below-barrier tunneling of the proton adjoining donor/acceptor sites. Our findings provide compelling experimental evidence for ultrafast hydron-migration events commensurate with the onset of low-barrier hydrogen bonding and afford a trenchant glimpse of molecular phenomena taking place at the "tipping point" between disparate dynamical regimes.

N-oxide Derivatives: Car-Parrinello Molecular Dynamics and Electron Localization Function Study on Intramolecular Hydrogen Bonds

First-principle molecular dynamics simulations (CPMD) and electron localization function (ELF) topological analysis were performed for two N-oxides, 2-[( N, N-dimethylamino- N-oxy)methyl]-4,6-dimethylphenol (I) and 5,5'-dibromo-3-[(diethylamino)methyl]-2,2'-biphenol N-oxide (II). Special attention was paid to hydrogen bond dynamics and spectroscopic features. Simulations were carried out in vacuo and in the crystalline phase. It was found that for compound I proton transfer phenomena do not occur spontaneously. An opposite situation was noticed for one of the two intramolecular hydrogen bonds of compound II. The influence of the neighboring molecules (including microsolvation) and the crystal field on the hydrogen bond properties cannot be neglected. Hydrogen bond spectroscopic signatures were quite well reproduced by power spectra of atomic velocity. Metric and spectroscopic theoretical findings were compared with experimental data available, especially in a view of the fact that experimental X-ray diffraction and FT-IR data of compound II indicate strong delocalization of the N-oxide bridge proton. ELF analysis showed that for compound II the N-oxide bridge belongs to the topological category of strong hydrogen bonds with the separate ELF basin for the bridged hydrogen atom.

Solid-State 17 O NMR Reveals Hydrogen-Bonding Energetics: Not All Low-Barrier Hydrogen Bonds Are Strong

While NMR and IR spectroscopic signatures and structural characteristics of low-barrier hydrogen bond (LBHB) formation are well documented in the literature, direct measurement of the LBHB energy is difficult. Here, we show that solid-state ¹⁷ O NMR spectroscopy can provide unique information about the energy required to break a LBHB. Our solid-state ¹⁷ O NMR data show that the HB enthalpy of the O⋅⋅⋅H⋅⋅⋅N LBHB formed in crystalline nicotinic acid is only 7.7±0.5 kcal mol^-1 , suggesting that not all LBHBs are particularly strong.

Short Carboxylic Acid-Carboxylate Hydrogen Bonds Can Have Fully Localized Protons

Short hydrogen bonds (H-bonds) have been proposed to play key functional roles in several proteins. The location of the proton in short H-bonds is of central importance, as proton delocalization is a defining feature of low-barrier hydrogen bonds (LBHBs). Experimentally determining proton location in H-bonds is challenging. Here, bond length analysis of atomic (1.15-0.98 Å) resolution X-ray crystal structures of the human protein DJ-1 and its bacterial homologue, YajL, was used to determine the protonation states of H-bonded carboxylic acids. DJ-1 contains a buried, dimer-spanning 2.49 Å H-bond between Glu15 and Asp24 that satisfies standard donor-acceptor distance criteria for a LBHB. Bond length analysis indicates that the proton is localized on Asp24, excluding a LBHB at this location. However, similar analysis of the Escherichia coli homologue YajL shows both residues may be protonated at the H-bonded oxygen atoms, potentially consistent with a LBHB. A Protein Data Bank-wide screen identifies candidate carboxylic acid H-bonds in approximately 14% of proteins, which are typically short [⟨d_O-O⟩ = 2.542(2) Å]. Chemically similar H-bonds between hydroxylated residues (Ser/Thr/Tyr) and carboxylates show a trend of lengthening O-O distance with increasing H-bond donor pK_a. This trend suggests that conventional electronic effects provide an adequate explanation for short, charge-assisted carboxylic acid-carboxylate H-bonds in proteins, without the need to invoke LBHBs in general. This study demonstrates that bond length analysis of atomic resolution X-ray crystal structures provides a useful experimental test of certain candidate LBHBs.

Less stable tautomers form stronger hydrogen bonds: the case of water complexes

Hydrogen bonding in cyclic complexes of water with tautomeric pairs of molecules M⁰ and M¹ is calculated to be stronger by more than 25% for the less stable tautomer M¹ in all cases where the energy gap between the two tautomers is large (ΔE(M⁰ − M¹) > 10 kcal mol⁻¹).

Proton-Coupled Electron Transfer and a Tyrosine-Histidine Pair in a Photosystem II-Inspired β-Hairpin Maquette: Kinetics on the Picosecond Time Scale

Photosystem II (PSII) and ribonucleotide reductase employ oxidation and reduction of the tyrosine aromatic ring in radical transport pathways. Tyrosine-based reactions involve either proton-coupled electron transfer (PCET) or electron transfer (ET) alone, depending on the pH and the pKa of tyrosine's phenolic oxygen. In PSII, a subset of the PCET reactions are mediated by a tyrosine-histidine redox-driven proton relay, YD-His189. Peptide A is a PSII-inspired β-hairpin, which contains a single tyrosine (Y5) and histidine (H14). Previous electrochemical characterization indicated that Peptide A conducts a net PCET reaction between Y5 and H14, which have a cross-strand π-π interaction. The kinetic impact of H14 has not yet been explored. Here, we address this question through time-resolved absorption spectroscopy and 280-nm photolysis, which generates a neutral tyrosyl radical. The formation and decay of the neutral tyrosyl radical at 410 nm were monitored in Peptide A and its variant, Peptide C, in which H14 is replaced by cyclohexylalanine (Cha14). Significantly, both electron transfer (ET, pL 11, L = lyonium) and PCET (pL 9) were accelerated in Peptide A and C, compared to model tyrosinate or tyrosine at the same pL. Increased electronic coupling, mediated by the peptide backbone, can account for this rate acceleration. Deuterium exchange gave no significant solvent isotope effect in the peptides. At pL 9, but not at pL 11, the reaction rate decreased when H14 was mutated to Cha14. This decrease in rate is attributed to an increase in reorganization energy in the Cha14 mutant. The Y5-H14 mechanism in Peptide A is reminiscent of proton- and electron-transfer events involving YD-H189 in PSII. These results document a mechanism by which proton donors and acceptors can regulate the rate of PCET reactions.

Sequence-Dependent Conformational Heterogeneity and Proton-Transfer Reactivity of the Fluorescent Guanine Analogue 6-Methyl Isoxanthopterin (6-MI) in DNA

The local conformations of individual nucleic acid bases in DNA are important components in processes fundamental to gene regulation. Fluorescent nucleic acid base analogues, which can be substituted for natural bases in DNA, can serve as useful spectroscopic probes of average local base conformation and conformational heterogeneity. Here we report excitation-emission peak shift (EES) measurements of the fluorescent guanine (G) analogue 6-methyl isoxanthoptherin (6-MI), both as a ribonucleotide monophosphate (NMP) in solution and as a site-specific substituent for G in various DNA constructs. Changes in the peak positions of the fluorescence spectra as a function of excitation energy indicate that distinct subpopulations of conformational states exist in these samples on time scales longer than the fluorescence lifetime. Our pH-dependent measurements of the 6-MI NMP in solution show that these states can be identified as protonated and deprotonated forms of the 6-MI fluorescent probe. We implement a simple two-state model, which includes four vibrationally coupled electronic levels to estimate the free energy change, the free energy of activation, and the equilibrium constant for the proton transfer reaction. These parameters vary in single-stranded and duplex DNA constructs, and also depend on the sequence context of flanking bases. Our results suggest that proton transfer in 6-MI-substituted DNA constructs is coupled to conformational heterogeneity of the probe base, and can be interpreted to suggest that Watson-Crick base pairing between 6-MI and its complementary cytosine in duplex DNA involves a "low-barrier-hydrogen-bond". These findings may be important in using the 6-MI probe to understand local base conformational fluctuations, which likely play a central role in protein-DNA and ligand-DNA interactions.

Investigation of the resonance-assisted hydrogen bond in model β-diketones through localized molecular orbital analysis of the spin-spin coupling constants related to the O-H···O hydrogen bond

The resonance-assisted hydrogen bond (HB) phenomenon has been studied theoretically by a localized molecular orbital (LMO) decomposition of the spin-spin coupling constants between atoms either involved or close to the O-H···O system of some β-diketones and their saturated counterparts. The analysis, carried out at the level of the second-order polarization propagator approximation, shows that the contributions in terms of LMO to the paramagnetic spin orbital and the spin dipolar Ramsey terms proof the importance of the delocalized π-electron structure supporting the idea of the existence of the resonance-assisted HB phenomenon phenomenon. The LMO contributions to the Fermi contact term indicate mainly the presence of the HB that may or not be linked to the π-electrons.

Isotopomeric polymorphism in a “doubly-polymorphic” multi-component molecular crystal

Isotopomeric polymorphism is observed in complexes of isonicotinamide with oxalic acid, highly unusual here in that each isotopic complex is itself polymorphic, a situation of “double polymorphism”. The four polymorphic forms exhibit different degrees of hydron transfer.

Intramolecular charge-assisted hydrogen bond strength in pseudochair carboxyphosphate

Carboxyphosphate, a suspected intermediate in ATP-dependent carboxylases, has not been isolated nor observed directly by experiment. Consequently, little is known concerning its structure, stability, and ionization state. Recently, carboxyphosphate as either a monoanion or dianion has been shown computationally to adopt a novel pseudochair conformation featuring an intramolecular charge-assisted hydrogen bond (CAHB). In this work, additive and subtractive correction schemes to the commonly employed open–closed method are used to estimate the strength of the CAHB. Truhlar’s Minnesota M06-2X functional with Dunning’s aug-cc-pVTZ basis set has been used for geometry optimization, energy evaluation, and frequency analysis. The CHARMM force field has been used to approximate the Pauli repulsive terms in the closed and open forms of carboxyphosphate. From our additive correction scheme, differential Pauli repulsion contributions between the pseudochair (closed) and open conformations of carboxyphosphate are found to be significant in determining the CAHB strength. The additive correction modifies the CAHB prediction (ΔE_{closed–open}) of −14 kcal/mol for the monoanion and −12 kcal/mol for the dianion to −22.9 and −18.4 kcal/mol, respectively. Results from the subtractive technique reinforce those from our additive procedure, where the predicted CAHB strength ranges from −17.8 to −25.4 kcal/mol for the monoanion and from −15.7 to −20.9 kcal/mol for the dianion. Ultimately, we find that the CAHB in carboxyphosphate meets the criteria for short-strong hydrogen bonds. However, carboxyphosphate has a unique energy profile that does not result in the symmetric double-well behavior of low-barrier hydrogen bonds. These findings provide deeper insight into the pseudochair conformation of carboxyphosphate, and lead to an improved mechanistic understanding of this intermediate in ATP-dependent carboxylases.

Using unnatural amino acids to probe the energetics of oxyanion hole hydrogen bonds in the ketosteroid isomerase active site

Hydrogen bonds are ubiquitous in enzyme active sites, providing binding interactions and stabilizing charge rearrangements on substrate groups over the course of a reaction. But understanding the origin and magnitude of their catalytic contributions relative to hydrogen bonds made in aqueous solution remains difficult, in part because of complexities encountered in energetic interpretation of traditional site-directed mutagenesis experiments. It has been proposed for ketosteroid isomerase and other enzymes that active site hydrogen bonding groups provide energetic stabilization via "short, strong" or "low-barrier" hydrogen bonds that are formed due to matching of their pKa or proton affinity to that of the transition state. It has also been proposed that the ketosteroid isomerase and other enzyme active sites provide electrostatic environments that result in larger energetic responses (i.e., greater "sensitivity") to ground-state to transition-state charge rearrangement, relative to aqueous solution, thereby providing catalysis relative to the corresponding reaction in water. To test these models, we substituted tyrosine with fluorotyrosines (F-Tyr's) in the ketosteroid isomerase (KSI) oxyanion hole to systematically vary the proton affinity of an active site hydrogen bond donor while minimizing steric or structural effects. We found that a 40-fold increase in intrinsic F-Tyr acidity caused no significant change in activity for reactions with three different substrates. F-Tyr substitution did not change the solvent or primary kinetic isotope effect for proton abstraction, consistent with no change in mechanism arising from these substitutions. The observed shallow dependence of activity on the pKa of the substituted Tyr residues suggests that the KSI oxyanion hole does not provide catalysis by forming an energetically exceptional pKa-matched hydrogen bond. In addition, the shallow dependence provides no indication of an active site electrostatic environment that greatly enhances the energetic response to charge accumulation, consistent with prior experimental results.

Are there really low-barrier hydrogen bonds in proteins? The case of photoactive yellow protein

For a long time, low-barrier hydrogen bonds (LBHBs) have been proposed to exist in many enzymes and to play an important role in their catalytic function, but the proof of their existence has been elusive. The transient formation of an LBHB in a protein system has been detected for the first time using neutron diffraction techniques on a photoactive yellow protein (PYP) crystal in a study published in 2009 (Yamaguchi, S.; et al. Proc. Natl. Acad. Sci. U.S.A. 2009, 106, 440-444). However, very recent theoretical studies based on electronic structure calculations and NMR resonance experiments on PYP in solution (Saito, K.; et al. Proc. Natl. Acad. Sci. U.S.A. 2012, 109, 167-172) strongly indicate that there is not such an LBHB. By means of electronic structure calculations combined with the solution of the nuclear Schrödinger equation, we analyze here under which conditions an LBHB can exist in PYP, thus leading to a more reasonable and conciliating understanding of the above-mentioned studies.

Car–Parrinello and path integral molecular dynamics study of the intramolecular hydrogen bonds in the crystals of benzoylacetone and dideuterobenzoylacetone

The dynamics of the intramolecular short hydrogen bond in the molecular crystal of benzoylacetone and its deuterated analogue are investigated using ab initio molecular dynamics simulations.

Exploring Histidine Conformations in the M2 Channel Lumen of the Influenza A Virus at Neutral pH via Molecular Simulations

The pH-regulated M2 proton channel from the influenza A virus has a His-tetrad in its transmembrane (TM) domain that is essential for proton conduction. At neutral pH, the tetrad has been observed in two distinct configurations, the "His-box" and "dimer-of-dimers", with similar backbone structures suggesting competing models for proton conduction. Here, we propose that both conformations can play a role. In support of this hypothesis, we used molecular dynamics simulations based on density functional theory to simulate the M2-TM domain and force-field-based simulations to estimate the relevant free-energy barriers. Both configurations are stable on accessible simulation time scales, and transitions between them occur faster than the millisecond time scale of proton conduction. Moreover, the deprotonation energy is too high for spontaneous conduction, consistent with their occurrence in the low-current regime. Our computations support a multiconfiguration model with different population levels, thereby connecting experimental data obtained under different conditions.

STUDY OF INTRAMOLECULAR HYDROGEN BONDING IN ORTHO-SUBSTITUTED ACETANILIDE COMPOUNDS USING COMPUTATIONAL METHODS

The intramolecular hydrogen bonding formation in ortho-substituted compounds of Acetanilide, ortho-hydroxy Acetanilide and ortho-nitro Acetanilide, was investigated using Density Functional Theory (DFT), Møller-Plesset second-order (MP2) method and "Atoms in Molecules (AIM)" theory. It was found that in each case, the cis isomer is more stable than the trans isomer and ortho-nitro Acetanilide forms a stronger hydrogen bond than ortho-hydroxy Acetanilide. The effects of hydrogen bonding on structural parameters of the considered systems were studied using Becke's functional (B3LYP) and at the ab initio MP2 level in conjunction with different basis sets and suitable structural factors. The results are in agreement with the results of AIM theory.

Kinetic, mechanistic, and structural modeling studies of truncated wild-type leucine-rich repeat kinase 2 and the G2019S mutant

Leucine-rich repeat kinase 2 (LRRK2), a large and complex protein that possesses two enzymatic properties, kinase and GTPase, is one of the major genetic factors in Parkinson's disease (PD). Here, we characterize the kinetic and catalytic mechanisms of truncated wild-type (t-wt) LRRK2 and its most common mutant, G2019S (t-G2019S), with a structural interpretation of the kinase domain. First, the substitution of threonine with serine in the LRRKtide peptide results in a much less efficient substrate as demonstrated by a 26-fold decrease in k(cat) and a 6-fold decrease in binding affinity. The significant decrease in k(cat) is attributed to a slow chemical transfer step as evidenced by the inverse solvent kinetic isotope effect in the proton inventory and pL (pH or pD)-dependent studies. The shape of the proton inventory and pL profile clearly signals the involvement of a general base (pK(a) = 7.5) in the catalysis with a low fractionation factor in the ground state. We report for the first time that the increased kinase activity of the G2019S mutant is substrate-dependent. Homology modeling of the kinase domain (open and closed forms) and structural analysis of the docked peptide substrates suggest that electrostatic interactions play an important role in substrate recognition, which is affected by G2019S and may directly influence the kinetic properties of the enzyme. Finally, the GTPase activity of the t-G2019S mutant was characterized, and the mutation modestly decreases GTPase activity without significantly affecting GTP binding affinity.

Proton coupled electron transfer and redox-active tyrosine Z in the photosynthetic oxygen-evolving complex

Proton coupled electron transfer (PCET) reactions play an essential role in many enzymatic processes. In PCET, redox-active tyrosines may be involved as intermediates when the oxidized phenolic side chain deprotonates. Photosystem II (PSII) is an excellent framework for studying PCET reactions, because it contains two redox-active tyrosines, YD and YZ, with different roles in catalysis. One of the redox-active tyrosines, YZ, is essential for oxygen evolution and is rapidly reduced by the manganese-catalytic site. In this report, we investigate the mechanism of YZ PCET in oxygen-evolving PSII. To isolate YZ(•) reactions, but retain the manganese-calcium cluster, low temperatures were used to block the oxidation of the metal cluster, high microwave powers were used to saturate the YD(•) EPR signal, and YZ(•) decay kinetics were measured with EPR spectroscopy. Analysis of the pH and solvent isotope dependence was performed. The rate of YZ(•) decay exhibited a significant solvent isotope effect, and the rate of recombination and the solvent isotope effect were pH independent from pH 5.0 to 7.5. These results are consistent with a rate-limiting, coupled proton electron transfer (CPET) reaction and are contrasted to results obtained for YD(•) decay kinetics at low pH. This effect may be mediated by an extensive hydrogen-bond network around YZ. These experiments imply that PCET reactions distinguish the two PSII redox-active tyrosines.

Three-dimensional structure and catalytic mechanism of cytosine deaminase

Cytosine deaminase (CDA) from E. coli is a member of the amidohydrolase superfamily. The structure of the zinc-activated enzyme was determined in the presence of phosphonocytosine, a mimic of the tetrahedral reaction intermediate. This compound inhibits the deamination of cytosine with a K(i) of 52 nM. The zinc- and iron-containing enzymes were characterized to determine the effect of the divalent cations on activation of the hydrolytic water. Fe-CDA loses activity at low pH with a kinetic pK(a) of 6.0, and Zn-CDA has a kinetic pK(a) of 7.3. Mutation of Gln-156 decreased the catalytic activity by more than 5 orders of magnitude, supporting its role in substrate binding. Mutation of Glu-217, Asp-313, and His-246 significantly decreased catalytic activity supporting the role of these three residues in activation of the hydrolytic water molecule and facilitation of proton transfer reactions. A library of potential substrates was used to probe the structural determinants responsible for catalytic activity. CDA was able to catalyze the deamination of isocytosine and the hydrolysis of 3-oxauracil. Large inverse solvent isotope effects were obtained on k(cat) and k(cat)/K(m), consistent with the formation of a low-barrier hydrogen bond during the conversion of cytosine to uracil. A chemical mechanism for substrate deamination by CDA was proposed.