Low-barrier Hydrogen Bond Research Articles

Constrained nuclear-electronic orbital QM/MM approach for simulating complex systems with quantum nuclear delocalization effects incorporated

The hybrid quantum mechanics/molecular mechanics (QM/MM) approach, which combines the accuracy of QM methods with the efficiency of MM methods, is widely used in the study of complex systems. However, past QM/MM implementations often neglect or face challenges in addressing nuclear quantum effects, despite their crucial role in many key chemical and biological processes. Recently, our group developed the constrained nuclear-electronic orbital (CNEO) theory, a cost-efficient approach that accurately addresses nuclear quantum effects, especially quantum nuclear delocalization effects. In this work, we integrate CNEO with the QM/MM approach through the electrostatic embedding scheme and apply the resulting CNEO QM/MM to two hydrogen-bonded complexes. We find that both solvation effects and nuclear quantum effects significantly impact hydrogen bond structures and dynamics. Notably, in the glutamic acid–glutamate complex, which mimics a common low barrier hydrogen bond in biological systems, CNEO QM/MM accurately predicts nearly equal proton sharing between the two residues. With an accurate description of both quantum nuclear delocalization effects and environmental effects, CNEO QM/MM is a promising new approach for simulating complex chemical and biological systems.

A low-barrier proton shared between two aspartates acts as a conformational switch that changes the substrate specificity of the β-lactamase BlaC

Serine β-lactamases inactivate β-lactam antibiotics in a two-step mechanism comprising acylation and deacylation. For the deacylation step, a water molecule is activated by a conserved glutamate residue to release the adduct from the enzyme. The third-generation cephalosporin ceftazidime is a poor substrate for the class A β-lactamase BlaC from Mycobacterium tuberculosis but it can be hydrolyzed faster when the active site pocket is enlarged, as was reported for mutant BlaC P167S. The conformational change in the Ω-loop of the P167S mutant displaces the conserved glutamate (Glu166), suggesting it is not required for deacylation of the ceftazidime adduct. Here, we report the characterization of wild type BlaC and BlaC E166A at various pH values. The presence of Glu166 strongly enhances activity against nitrocefin but not ceftazidime, indicating it is indeed not required for deacylation of the adduct of the latter substrate. At high pH wild type BlaC was found to exist in two states, one of which converts ceftazidime much faster, resembling the open state previously reported for the BlaC mutant P167S. The pH-dependent switch between the closed and open states is caused by the loss at high pH of a low-barrier hydrogen bond, a proton shared between Asp172 and Asp179. These results illustrate how readily shifts in substrate specificity can occur as a consequence of subtle changes in protein structure.

Dynamic interchange between two protonation states is characteristic of active sites of cholinesterases.

Cholinesterases are well-known and widely studied enzymes crucial to human health and involved in neurology, Alzheimer's, and lipid metabolism. The protonation pattern of active sites of cholinesterases influences all the chemical processes within, including reaction, covalent inhibition by nerve agents, and reactivation. Despite its significance, our comprehension of the fine structure of cholinesterases remains limited. In this study, we employed enhanced-sampling quantum-mechanical/molecular-mechanical calculations to show that cholinesterases predominantly operate as dynamic mixtures of two protonation states. The proton transfer between two non-catalytic glutamate residues follows the Grotthuss mechanism facilitated by a mediator water molecule. We show that this uncovered complexity of active sites presents a challenge for classical molecular dynamics simulations and calls for special treatment. The calculated proton transfer barrier of 1.65 kcal/mol initiates a discussion on the potential existence of two coupled low-barrier hydrogen bonds in the inhibited form of butyrylcholinesterase. These findings expand our understanding of structural features expressed by highly evolved enzymes and guide future advances in cholinesterase-related protein and drug design studies.

Solvatochromism mechanism of perylene diimide radicals stabilized by low barrier hydrogen bond and detection of formaldehyde in real food samples

Formaldehyde (FA) is widely used for food preservation and whitening because of its preservative and bleaching effects. Excessive FA seriously endangers human health, so monitoring FA-treated food is of great significance for food safety. In this paper, perylene diimide radicals with the chirality (BTFPDI) and water-solubility (TFPDI-bPEI) were designed respectively based on our previous studies. Property of solvatochromism induced by these perylene diimide radicals aggregation was demonstrated by the concentration and temperature dependent absorption spectrum and the circular dichroism spectrum. Further the colorimetric monitoring of trace FA in gaseous and liquid was realized by the condensation reaction of FA with amine to induce TFPDI-bPEI aggregation. After the addition of FA, TFPDI-bPEI molecule aggregated in aqueous solution due to the specific reaction of a large number of –NH2 groups with FA, resulting in a blue-shift of the maximum absorption peak from 780 nm to 561 nm, and the detection limit as low as 0.86 μM to FA and the response time of only 20 s. This method could further detect FA in food with good recoveries and detect FA in contaminated food samples.

Electronic Fingerprint of the Protonated Imidazole Dimer Probed by X-ray Absorption Spectroscopy.

Protons in low-barrier superstrong hydrogen bonds are typically delocalized between two electronegative atoms. Conventional methods to characterize such superstrong hydrogen bonds are vibrational spectroscopy and diffraction techniques. We introduce soft X-ray spectroscopy to uncover the electronic fingerprints for proton sharing in the protonated imidazole dimer, a prototypical building block enabling effective proton transport in biology and high-temperature fuel cells. Using nitrogen core excitations as a sensitive probe for the protonation status, we identify the X-ray signature of a shared proton in the solvated imidazole dimer in a combined experimental and theoretical approach. The degree of proton sharing is examined as a function of structural variations that modify the shape of the low-barrier potential in the superstrong hydrogen bond. We conclude by showing how the sensitivity to the quantum distribution of proton motion in the double-well potential is reflected in the spectral signature of the shared proton.

Low-Barrier Hydrogen Bond Determines Target-Binding Affinity and Specificity of the Antitubercular Drug Bedaquiline.

The role of short strong hydrogen bonds (SSHBs) in ligand-target binding remains largely unexplored, thereby hindering a potentially important avenue in rational drug design. Here we investigate the interaction between the antituberculosis drug bedaquiline (Bq) and the mycobacterial ATP synthase to unravel the role of a specific hydrogen bond to a conserved acidic residue in the target affinity and specificity. Our ab initio molecular dynamics simulations reveal that this bond belongs to the SSHB category and accounts for a substantial fraction of the target binding free energy. We also demonstrate that the presence of an extra acidic residue, i.e., aspartic acid at position 32 (D32), found exclusively in mycobacteria, cooperatively enhances the HB strength, ensuring specificity for the mycobacterial target. Consistently, we show that the removal of D32 markedly weakens the affinity, leading to Bq resistance associated with mutations of D32 to nonacidic residues. By designing simple Bq analogs, we then explore the possibility to overcome the resistance and potentially broaden the Bq antimicrobial spectrum by making the SSHB independent of the presence of the extra acidic residue.

Molecular insights into the water dissociation and proton dynamics at the β-TaON (100)/water interface.

Understanding the dynamic nature of the semiconductor-water interface is crucial for developing efficient photoelectrochemical water splitting catalysts, as it governs reactivity through charge and mass transport. In this study, we employ ab initio molecular dynamics simulations to investigate the structural and dynamical properties of water at the β-TaON (100) surface. We observed that a well-defined interface is established through the spontaneous dissociation of water and the reorganization of surface chemical bonds. This leads to the formation of a partially hydroxylated surface, accompanied by a strong network of hydrogen bonds at the TaON-water interface. Consequently, various proton transport routes, including the proton transfer through "low-barrier hydrogen bond" path, become active across the interface, dramatically increasing the overall rate of the proton hopping at the interface. Based on our findings, we propose that the observed high photocatalytic activity of TaON-based semiconductors could be attributed to the spontaneous water dissociation and the resulting high proton transfer rate at the interface.

Interplay of two low-barrier hydrogen bonds in long-distance proton-coupled electron transfer for water oxidation.

D1-Tyr161 (TyrZ) forms a low-barrier H-bond with D1-His190 and functions as a redox-active group in photosystem II. When oxidized to the radical form (TyrZ-O•), it accepts an electron from the oxygen-evolving Mn4CaO5 cluster, facilitating an increase in the oxidation state (Sn; n = 0-3). In this study, we investigated the mechanism of how TyrZ-O• drives proton-coupled electron transfer during the S2 to S3 transition using a quantum mechanical/molecular mechanical approach. In response to TyrZ-O• formation and subsequent loss of the low-barrier H-bond, the ligand water molecule at the Ca2+ site (W4) reorients away from TyrZ and donates an H-bond to D1-Glu189 at Mn4 of Mn4CaO5 together with an adjacent water molecule. The H-bond donation to the Mn4CaO5 cluster triggers the release of the proton from the lowest pKa site (W1 at Mn4) along the W1…D1-Asp61 low-barrier H-bond, leading to protonation of D1-Asp61. The interplay of the two low-barrier H-bonds, involving the Ca2+ interface and forming the extended Grotthuss-like network [TyrZ…D1-His190]-[Mn4CaO5]-[W1…D1-Asp61], rather than the direct electrostatic interaction, is likely a basis of the apparent long-distance interaction (11.4 Å) between TyrZ-O• formation and D1-Asp61 protonation.

Selective Electrochemical Oxidation of Benzylic C-H to Benzylic Alcohols with the Aid of Imidazolium Radical Mediators.

Selective oxidation of benzylic C-H to benzylic alcohols is a well-known challenge in the chemical community since benzylic C-H is more prone to be overoxidized to benzylic ketones. In this work, we report the highly selective electro-oxidation of benzylic C-H to benzylic alcohols in an undivided cell in ionic liquid-based solution. As an example, the selectivity toward xanthydrol could be as high as 95.7% at complete conversion of xanthene, a typical benzylic C-H compound, on gram-scale in imidazolium bromide/H2O/DMF. Mechanism investigation reveals that the imidazolium radical generated in situ participants in a proton-coupled electron transfer process and low-barrier hydrogen bonds stabilize the reaction intermediates, together steering the redox equilibrium, favoring benzylic alcohols over benzylic ketones.

Free-Base Corrole Anion.

Free-base corroles have long been known to be acidic, readily undergoing deprotonation by mild bases and in polar solvents. The conjugate base, however, has not been structurally characterized until now. Presented here is a first crystal structure of a free-base corrole anion, derived from tris(p-cyanophenyl)corrole, as the tetrabuylammonium salt. The low-temperature (100 K) structure reveals localized hydrogens on a pair of opposite pyrrole nitrogens. DFT calculations identify such a structure as the global minimum but also point to two cis tautomers only 4-7 kcal/mol above the ground state. In terms of free energy, however, the cis tautomers are above or essentially flush with the trans-to-cis barrier so the cis tautomers are unlikely to exist or be observed as true intermediates. Thus, the hydrogen bond within each dipyrrin unit on either side of the molecular pseudo-C2 axis through C10 (i.e., between pyrrole rings A and B or between C and D) qualifies as or closely approaches a low-barrier hydrogen bond. Proton migration across the pseudo-C2 axis entails much higher activation energies >20 kcal/mol, reflecting the relative rigidity of the molecule along the C1-C19 pyrrole-pyrrole linkage.

Human 2'-Deoxynucleoside 5'-Phosphate N-Hydrolase 1: Mechanism of 2'-Deoxyuridine 5'-Monophosphate Hydrolysis.

The enzyme 2'-deoxynucleoside 5'-phosphate N-hydrolase 1 (DNPH1) catalyzes the N-ribosidic bond cleavage of 5-hydroxymethyl-2'-deoxyuridine 5'-monophosphate to generate 2-deoxyribose 5-phosphate and 5-hydroxymethyluracil. DNPH1 accepts other 2'-deoxynucleoside 5'-monophosphates as slow-reacting substrates. DNPH1 inhibition is a promising strategy to overcome resistance to and potentiate anticancer poly(ADP-ribose) polymerase inhibitors. We solved the crystal structure of unliganded human DNPH1 and took advantage of the slow reactivity of 2'-deoxyuridine 5'-monophosphate (dUMP) as a substrate to obtain a crystal structure of the DNPH1:dUMP Michaelis complex. In both structures, the carboxylate group of the catalytic Glu residue, proposed to act as a nucleophile in covalent catalysis, forms an apparent low-barrier hydrogen bond with the hydroxyl group of a conserved Tyr residue. The crystal structures are supported by functional data, with liquid chromatography-mass spectrometry analysis showing that DNPH1 incubation with dUMP leads to slow yet complete hydrolysis of the substrate. A direct UV-vis absorbance-based assay allowed characterization of DNPH1 kinetics at low dUMP concentrations. A bell-shaped pH-rate profile indicated that acid-base catalysis is operational and that for maximum kcat/KM, two groups with an average pKa of 6.4 must be deprotonated, while two groups with an average pKa of 8.2 must be protonated. A modestly inverse solvent viscosity effect rules out diffusional processes involved in dUMP binding to and possibly uracil release from the enzyme as rate limiting to kcat/KM. Solvent deuterium isotope effects on kcat/KM and kcat were inverse and unity, respectively. A reaction mechanism for dUMP hydrolysis is proposed.

Proton-Tunneling Dynamics along Low-Barrier Hydrogen Bonds: A Full-Dimensional Instanton Study of 6-Hydroxy-2-formylfulvene.

Understanding the dynamics of proton transfer along low-barrier hydrogen bonds remains an outstanding challenge of great fundamental and practical interest, reflecting the central role of quantum effects in reactions of chemical and biological importance. Here, we combine ab initio calculations with the semiclassical ring-polymer instanton method to investigate tunneling processes on the ground electronic state of 6-hydroxy-2-formylfulvene (HFF), a prototypical neutral molecule supporting low-barrier hydrogen-bonding. The results emerging from a full-dimensional ab initio instanton analysis reveal that the tunneling path does not pass through the instantaneous transition-state geometry. Instead, the tunneling process involves a multidimensional reaction coordinate with concerted reorganization of the heavy-atom skeletal framework to substantially reduce the donor-acceptor distance and drive the ensuing intramolecular proton-transfer event. The predicted tunneling-induced splittings for HFF isotopologues are in good agreement with experimental findings, leading to percentage deviations of only 20-40%. Our full-dimensional results allow us to characterize vibrational contributions along the tunneling path, highlighting the intrinsically multidimensional nature of the attendant hydron-migration dynamics.

Distinction between low-barrier hydrogen bond and ordinary hydrogen bond: a case study of varying nature of charge assisted hydrogen bonds of diglycine perchlorate crystal

Hydrogen bonding is a complex phenomenon that is a resultant of many energy components like the electrostatic, dispersive, covalent, charge cloud overlap repulsion etc, nature of hydrogen bond (H-bond) depends on which of these components play a dominant role. Low barrier hydrogen bond (LBHB) constitutes a special category of hydrogen bonds characterized by near delocalization of proton between donor and acceptor groups of the H- bond unlike an ordinary hydrogen bond (OHB) having proton clearly localized near the donor group. The significance of LBHBs in macromolecular interactions has been highly controversial, despite may attempts the existence and potential importance of protein LBHBs remains debatable. In order to answer questions like whether or not a distinct class of LBHBs exists and if they do exist under what conditions they are formed and how do they behave differently from OHBs, a detailed study of H-bonding in Diglycine Perchlorate (DGPCl) crystal containing five unique hydrogen bonded glycinium-glycine pairs is undertaken. All O-H–O bonds of DGPCl are between the carboxyl (-COOH) and carboxylate (-COO−) groups with slightly different electron distributions resulting in observable variations in the H-bond geometries, this is an indication of varying strength of these short strong H-bonds. It is found that LBHB nature of the five O-H—O bonds between glycinium-glycine pairs (P1-P5) varies as P1 < P4 < P2 < P3 < P5. This study gives an experimental evidence of the existence of LBHBs and demonstrates that the behaviour of LBHBs is very different from that of strong OHBs.

The Realm of Unconventional Noncovalent Interactions in Proteins: Their Significance in Structure and Function.

Proteins and their assemblies are fundamental for living cells to function. Their complex three-dimensional architecture and its stability are attributed to the combined effect of various noncovalent interactions. It is critical to scrutinize these noncovalent interactions to understand their role in the energy landscape in folding, catalysis, and molecular recognition. This Review presents a comprehensive summary of unconventional noncovalent interactions, beyond conventional hydrogen bonds and hydrophobic interactions, which have gained prominence over the past decade. The noncovalent interactions discussed include low-barrier hydrogen bonds, C5 hydrogen bonds, C-H···π interactions, sulfur-mediated hydrogen bonds, n → π* interactions, London dispersion interactions, halogen bonds, chalcogen bonds, and tetrel bonds. This Review focuses on their chemical nature, interaction strength, and geometrical parameters obtained from X-ray crystallography, spectroscopy, bioinformatics, and computational chemistry. Also highlighted are their occurrence in proteins or their complexes and recent advances made toward understanding their role in biomolecular structure and function. Probing the chemical diversity of these interactions, we determined that the variable frequency of occurrence in proteins and the ability to synergize with one another are important not only for ab initio structure prediction but also to design proteins with new functionalities. A better understanding of these interactions will promote their utilization in designing and engineering ligands with potential therapeutic value.

Symmetry of Hydrogen Bonds: Application of NMR Method of Isotopic Perturbation and Relevance of Solvatomers.

Short, strong, symmetric, low-barrier hydrogen bonds (H-bonds) are thought to be of special significance. We have been searching for symmetric H-bonds by using the NMR technique of isotopic perturbation. Various dicarboxylate monoanions, aldehyde enols, diamines, enamines, acid-base complexes, and two sterically encumbered enols have been investigated. Among all of these, we have found only one example of a symmetric H-bond, in nitromalonamide enol, and all of the others are equilibrating mixtures of tautomers. The nearly universal lack of symmetry is attributed to the presence of these H-bonded species as a mixture of solvatomers, meaning isomers (or stereoisomers or tautomers) that differ in their solvation environment. The disorder of solvation renders the two donor atoms instantaneously inequivalent, whereupon the hydrogen attaches to the less well solvated donor. We therefore conclude that there is no special significance to short, strong, symmetric, low-barrier H-bonds. Moreover, they have no heightened stability or else they would have been more prevalent.

Identification of the protonation and oxidation states of the oxygen-evolving complex in the low-dose X-ray crystal structure of photosystem II.

In photosystem II (PSII), the O3 and O4 sites of the Mn4CaO5 cluster form hydrogen bonds with D1-His337 and a water molecule (W539), respectively. The low-dose X-ray structure shows that these hydrogen bond distances differ between the two homogeneous monomer units (A and B) [Tanaka etal., J. Am Chem. Soc. 2017, 139, 1718]. We investigated the origin of the differences using a quantum mechanical/molecular mechanical (QM/MM) approach. QM/MM calculations show that the short O4-OW539 hydrogen bond (~2.5 Å) of the B monomer is reproduced when O4 is protonated in the S1 state. The short O3-NεHis337 hydrogen bond of the A monomer is due to the formation of a low-barrier hydrogen bond between O3 and doubly-protonated D1-His337 in the overreduced states (S-1 or S-2). It seems plausible that the oxidation state differs between the two monomer units in the crystal.

Inter- vs. Intra-Molecular Hydrogen Bond in Complexes of Nitrophthalic Acids with Pyridine.

This study covers the analysis of isomeric forms of nitrophthalic acids with pyridine. This work dwells on the complementary experimental (X-ray, IR and Raman) and theoretical (Car-Parrinello Molecular Dynamics (CPMD) and Density Functional Theory (DFT)) studies of the obtained complexes. The conducted studies showed that steric repulsion between the nitro group in ortho-position and the carboxyl group causes significant isomeric changes. Modeling of the nitrophthalic acid-pyridine complex yielded a short strong intramolecular hydrogen bond (SSHB). The transition energy from the isomeric form with an intermolecular hydrogen bond to the isomeric form with an intramolecular hydrogen bond was estimated.

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)CoIII(OIPh)(OH)]2+ (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)CoIII(OIPh)(HO---HOTf)]2+ (1LBHB) and [(TQA)CoIII(OIPh)(H2O--OTf-)]2+ (1'LBHB). Due to the oxo-wall, these complexes (1LBHB and 1'LBHB) cannot convert to high-valent cobalt-oxyl species. Instead, styrene oxidation by these oxidants (1LBHB 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-/H2O ligand and the iodine of PhIO modulates the activity of the cobalt-iodosylarene complexes 1LBHB 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.

Unraveling the Nature of Hydrogen Bonds of "Proton Sponges" Based on Car-Parrinello and Metadynamics Approaches.

The nature of intra- and intermolecular non-covalent interactions was studied in four naphthalene derivatives commonly referred to as "proton sponges". Special attention was paid to an intramolecular hydrogen bond present in the protonated form of the compounds. The unsubstituted "proton sponge" served as a reference structure to study the substituent influence on the hydrogen bond (HB) properties. We selected three compounds substituted by methoxy, amino, and nitro groups. The presence of the substituents either retained the parent symmetry or rendered the compounds asymmetric. In order to reveal the non-covalent interaction properties, the Hirshfeld surface (HS) was computed for the crystal structures of the studied compounds. Next, quantum-chemical simulations were performed in vacuo and in the crystalline phase. Car-Parrinello molecular dynamics (CPMD), Path Integral Molecular Dynamics (PIMD), and metadynamics were employed to investigate the time-evolution changes of metric parameters and free energy profile in both phases. Additionally, for selected snapshots obtained from the CPMD trajectories, non-covalent interactions and electronic structure were studied. Quantum theory of atoms in molecules (QTAIM) and the Density Overlap Regions Indicator (DORI) were applied for this purpose. It was found based on Hirshfeld surfaces that, besides intramolecular hydrogen bonds, other non-covalent interactions are present and have a strong impact on the crystal structure organization. The CPMD results obtained in both phases showed frequent proton transfer phenomena. The proton was strongly delocalized in the applied time-scale and temperature, especially in the PIMD framework. The use of metadynamics allowed for tracing the free energy profiles and confirming that the hydrogen bonds present in "proton sponges" are Low-Barrier Hydrogen Bonds (LBHBs). The electronic and topological analysis quantitatively described the temperature dependence and time-evolution changes of the electronic structure. The covalency of the hydrogen bonds was estimated based on QTAIM analysis. It was found that strong hydrogen bonds show greater covalency, which is additionally determined by the proton position in the hydrogen bridge.

Nuclear quantum and H/D isotope effects on intramolecular hydrogen bond in curcumin.

Curcumin and its derivatives possess intramolecular low-barrier hydrogen bonds for intramolecular proton transfer. The π-delocalization in the OCCCO framework of the hydrogen bond in these compounds is reorganized concomitantly with the proton transfer. To characterize the hydrogen bond and π-delocalization, we performed path integral molecular dynamics simulations, revealing that although the proton migration and reorganization of the π-delocalized structure showed a positive correlation, the correlation was weak.