Proton transfer in a short hydrogen bond caused by solvation shell fluctuations: an ab initio MD and NMR/UV study of an (OHO)(-) bonded system

Very Strong Hydrogen Bond in Nitrophthalic Cocrystals

This work presents the studies of a very strong hydrogen bond (VSHB) in biologically active phthalic acids. Research on VSHB comes topical due to its participation in many biological processes. The studies cover the modelling of intermolecular interactions and phthalic acids with 2,4,6-collidine and N,N-dimethyl-4-pyridinamine complexes with aim to obtain a VSHB. The four synthesized complexes were studied by experimental X-ray, IR, and Raman methods, as well as theoretical Car-Parrinello Molecular Dynamics (CP-MD) and Density Functional Theory (DFT) simulations. By variation of the steric repulsion and basicity of the complex' components, a very short intramolecular hydrogen bond was achieved. The potential energy curves calculated by the DFT method were characterized by a low barrier (0.7 and 0.9 kcal/mol) on proton transfer in the OHN intermolecular hydrogen bond for 3-nitrophthalic acid with either 2,4,6-collidine or N,N-dimethyl-4-pyridinamine cocrystals. Moreover, the CP-MD simulations exposed very strong bridging proton dynamics in the intermolecular hydrogen bonds. The accomplished crystallographic and spectroscopic studies indicate that the OHO intramolecular hydrogen bond in 4-nitrophthalic cocrystals is VSHB. The influence of a strong steric effect on the geometry of the studied cocrystals and the stretching vibration bands of the carboxyl and carboxylate groups was elaborated.

Mutual influence of non-covalent interactions formed by imidazole: A systematic quantum-chemical study

Imidazole is a five-membered heterocycle that is part of a number of biologically important molecules such as the amino acid histidine and the hormone histamine. Imidazole has a unique ability to participate in a variety of non-covalent interactions involving the NH group, the pyridine-like nitrogen atom or the π-system. For many biologically active compounds containing the imidazole moiety, its participation in formation of hydrogen bond NH⋯O/N and following proton transfer is the key step of mechanism of their action. In this work a systematic study of the mutual influence of various paired combinations of non-covalent interactions (e.g., hydrogen bonds and π-interactions) involving the imidazole moiety was performed by means of quantum chemistry (PW6B95-GD3/def2-QZVPD) for a series of model systems constructed based on analysis of available x-ray data. It is shown that for considered complexes formation of additional non-covalent interactions can only enhance the proton-donating ability of imidazole. At the same time, its proton-accepting ability can be both enhanced and weakened, depending on what additional interactions are added to a given system. The mutual influence of non-covalent interactions involving imidazole can be classified as weak geometric and strong energetic cooperativity-a small change in the length of non-covalent interaction formed by imidazole can strongly influence its strength. The latter can be used to develop methods for controlling the rate and selectivity of chemical reactions involving the imidazole fragment in larger systems. It is shown that the strong mutual influence of non-covalent interactions involving imidazole is due to the unique ability of the imidazole ring to effectively redistribute electron density in non-covalently bound systems with its participation.

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.

Absorption wavelength along chromophore low-barrier hydrogen bonds

In low-barrier hydrogen bonds (H-bonds), the pK_a values for the H-bond donor and acceptor moieties are nearly equal, whereas the redox potential values depend on the H⁺ position. Spectroscopic details of low-barrier H-bonds remain unclear. Here, we report the absorption wavelength along low-barrier H-bonds in protein environments, using a quantum mechanical/molecular mechanical approach. Low-barrier H-bonds form between Glu46 and p-coumaric acid (pCA) in the intermediate pR_CW state of photoactive yellow protein and between Asp116 and the retinal Schiff base in the intermediate M-state of the sodium-pumping rhodopsin KR2. The H⁺ displacement of only ∼0.4 Å, which does not easily occur without low-barrier H-bonds, is responsible for the ∼50 nm-shift in the absorption wavelength. This may be a basis of how photoreceptor proteins have evolved to proceed photocycles using abundant protons.

•

The low-barrier H-bond formation is a prerequisite for proton transfer

•

How the absorption wavelength changes as H⁺ moves is an open question

•

The H⁺ displacement of ∼0.4 Å leads to the absorption wavelength shift of ∼50 nm

•

The localization of the molecular orbitals plays a key role in the wavelength shift

Alkylammonium Cation Affinities of Nitrogenated Organobases: The Roles of Hydrogen Bonding and Proton Transfer

Alkylammonium cation affinities of 64 nitrogen-containing organobases, as well as the respective proton transfer processes from the alkylammonium cations to the base, have been computed in the gas phase by using DFT methods. The guanidine bases show the highest proton transfer values (191.9-233 kJ mol^-1 ) whereas the cis-2,2'-biimidazole presents the largest affinity towards the alkylammonium cations (>200 kJ mol^-1 ) values. The resulting data have been compared with the experimentally reported proton affinities of the studied nitrogen-containing organobases revealing that the propensity of an organobase for the proton transfer process increases linearly with its proton affinity. This work can provide a tool for designing senors for bioactive compounds containing amino groups that are protonated at physiological pH.

Combined NMR and UV-Vis Spectroscopic Studies of Models for the Hydrogen Bond System in the Active Site of Photoactive Yellow Protein: H-Bond Cooperativity and Medium Effects

Intramolecular hydrogen bonds in aprotic media were studied by combined (simultaneous) NMR and UV-vis spectroscopy. The species under investigation were anionic and featured single or coupled H-bonds between, for example, carboxylic groups and phenolic oxygen atoms (COO···H···OC)^-, among phenolic oxygen atoms (CO···H···OC)^-, and hydrogen bond chains between a carboxylic group and two phenolic oxygen atoms (COO···H···(OC)···H···OC)^-. The last anion may be regarded as a small molecule model for the hydrogen bond system in the active site of wild-type photoactive yellow protein (PYP) while the others mimic the corresponding H-bonds in site-selective mutants. Proton positions in isolated hydrogen bonds and hydrogen bond chains were assessed by calculations for vacuum conditions and spectroscopically for the two media, CD₂Cl₂ and the liquefied gas mixture CDClF₂/CDF₃ at low temperatures. NMR parameters allow for the estimation of time-averaged H-bond geometries, and optical spectra give additional information about geometry distributions. Comparison of the results from the various systems revealed the effects of the formation of hydrogen bond chains and changes of medium conditions on the geometry of individual H-bonds. In particular, the proton in a hydrogen bond to a carboxylic group shifts from the phenolic oxygen atom in the system COO^-···H-OC to the carboxylic group in COO-H···(OC)^-···H-OC as a result of hydrogen bond formation to the additional phenolic donor. Increase in medium polarity may, however, induce the conversion of a structure of a type COO-H···(OC)^-···H-OC to the type COO^-···H-(OC)···H-OC. Application of these results obtained from the model systems to PYP suggests that both cooperative effects within the hydrogen bond chain and a low-polarity protein environment are prerequisites for the stabilization of negative charge on the cofactor and hence for the spectral tuning of the photoreceptor.

Actual Symmetry of Symmetric Molecular Adducts in the Gas Phase, Solution and in the Solid State

This review discusses molecular adducts, whose composition allows a symmetric structure. Such adducts are popular model systems, as they are useful for analyzing the effect of structure on the property selected for study since they allow one to reduce the number of parameters. The main objectives of this discussion are to evaluate the influence of the surroundings on the symmetry of these adducts, steric hindrances within the adducts, competition between different noncovalent interactions responsible for stabilizing the adducts, and experimental methods that can be used to study the symmetry at different time scales. This review considers the following central binding units: hydrogen (proton), halogen (anion), metal (cation), water (hydrogen peroxide).

Sensitivity of 31 P NMR chemical shifts to hydrogen bond geometry and molecular conformation for complexes of phosphinic acids with pyridines

The results of the quantum-chemical investigation of a series of hydrogen-bonded 1:1 acid-base complexes formed by model phosphinic acids, Me₂ POOH, and PhHPOOH, are reported. A series of substituted pyridines (pK_a range from 0.5 to 10) was chosen as proton acceptors. Gradual changes of isotropic ³¹ P nuclear magnetic resonance (NMR) chemical shift, δP, were correlated with the bridging proton position in the intermolecular OHN hydrogen bond, namely, r (OH) distance; the proposed correlation could easily be extended to other phosphinic acids as well. For complexes with pyridine and 2,4,6-trimethylpyridine, we have investigated in more detail several factors influencing the δP values: (1) the proton transfer within the OHN hydrogen bond; (2) the rotation of the pyridine ring around the hydrogen bond axis (associated with the formation/breakage of additional weak PO···H-C hydrogen bond); and (3) the rotation of the phenyl substituent in phenylphosphinic acid around the P-C axis. All these factors appeared to be of similar magnitude, thus masking their individual contributions that have to be independently estimated for a reliable spectral interpretation.

Modeling of Solute-Solvent Interactions Using an External Electric Field-From Tautomeric Equilibrium in Nonpolar Solvents to the Dissociation of Alkali Metal Halides

An implicit account of the solvent effect can be carried out using traditional static quantum chemistry calculations by applying an external electric field to the studied molecular system. This approach allows one to distinguish between the effects of the macroscopic reaction field of the solvent and specific solute-solvent interactions. In this study, we report on the dependence of the simulation results on the use of the polarizable continuum approximation and on the importance of the solvent effect in nonpolar solvents. The latter was demonstrated using experimental data on tautomeric equilibria between the pyridone and hydroxypyridine forms of 2,6-di-tert-butyl-4-hydroxy-pyridine in cyclohexane and chloroform.

Adduct under Field-A Qualitative Approach to Account for Solvent Effect on Hydrogen Bonding

The location of a mobile proton in acid-base complexes in aprotic solvents can be predicted using a simplified Adduct under Field (AuF) approach, where solute–solvent effects on the geometry of hydrogen bond are simulated using a fictitious external electric field. The parameters of the field have been estimated using experimental data on acid-base complexes in CDF₃/CDClF₂. With some limitations, they can be applied to the chemically similar CHCl₃ and CH₂Cl₂. The obtained data indicate that the solute–solvent effects are critically important regardless of the type of complexes. The temperature dependences of the strength and fluctuation rate of the field explain the behavior of experimentally measured parameters.

Elucidating Solvent Effects on Strong Intramolecular Hydrogen Bond: DFT-MD Study of Dibenzoylmethane in Methanol Solution

The dynamic aspect of solvation plays a crucial role in determining properties of strong intramolecular hydrogen bonds since solvent fluctuations modify instantaneous hydrogen-bonded proton transfer barriers. Previous studies pointed out that solvent-solute interactions in the first solvation shell govern the position of the proton but the ability of the electric field due to other solvent molecules to localize the proton remains an important issue. In this work, we examine the structure of the O-H⋅⋅⋅O intramolecular hydrogen bond of dibenzoylmethane in methanol solution by employing density functional theory-based molecular dynamics and quantum chemical calculations. Our computations showed that homogeneous electric fields with intensities corresponding to those found in polar solvents are able to considerably alter the proton transfer barrier height in the gas phase. In methanol solution, the proton position is correlated with the difference in electrostatic potentials on the oxygen atoms of dibenzoylmethane even when dibenzoylmethane-methanol hydrogen bonding is lacking. On a timescale of our simulation, the hydrogen bonding and solvent electrostatics tend to localize the proton on different oxygen atoms. These findings provide an insight into the importance of the solvent electric field on the structure of a strong intramolecular hydrogen bond.

Influence of Hydrogen Bonds in 1:1 Complexes of Phosphinic Acids with Substituted Pyridines on 1H and 31P NMR Chemical Shifts

Two series of 1:1 complexes with strong OHN hydrogen bonds formed by dimethylphosphinic and phenylphosphinic acids with 10 substituted pyridines were studied experimentally by liquid state NMR spectroscopy at 100 K in solution in a low-freezing polar aprotic solvent mixture CDF₃/CDClF₂. The hydrogen bond geometries were estimated using previously established correlations linking ¹H NMR chemical shifts of bridging protons with the O···H and H···N interatomic distances. A new correlation is proposed allowing one to estimate the interatomic distance within the OHN bridge from the displacement of ³¹P NMR signal upon complexation. We show that the values of ³¹P NMR chemical shifts are affected by an additional CH···O hydrogen bond formed between the P═O group of the acid and ortho-CH proton of the substituted pyridines. Breaking of this bond in the case of 2,6-disubstituted bases shifts the ³¹P NMR signal by ca. 1.5 ppm to the high field.

Libration of phenyl groups detected by VT-SSNMR: Comparison with X-ray crystallography

The X-ray crystal structure of 2-benzyl-1H-benzimidazole, 2BnBzIm, was determined at 293 K showing no dynamic phenomena (disorder) of any class. On the other hand, some ¹³ C NMR signals were absent in the CPMAS spectrum (100 MHz, 300 K). We decided to carry out variable-temperature SSNMR and discovered that the missing signals are ortho and meta carbons of the phenyl ring of the benzyl group. Line-shape analysis and the Eyring equation were used to determine the barrier, which was compared with the calculated DFT for the gas phase that it is much lower.

Simplified calculation approaches designed to reproduce the geometry of hydrogen bonds in molecular complexes in aprotic solvents

The impact of the environment onto the geometry of hydrogen bonds can be critically important for the properties of the questioned molecular system. The paper reports on the design of calculation approaches capable to simulate the effect of aprotic polar solvents on the geometric and NMR parameters of intermolecular hydrogen bonds. A hydrogen fluoride and pyridine complex has been used as the main model system because the experimental estimates of these parameters are available for it. Specifically, F-H, F⋯N, and H-N distances, the values of ¹⁵N NMR shift, and spin-spin coupling constants ¹J(¹⁹F¹H), ^1hJ(¹H¹⁵N), and ^2hJ(¹⁹F¹⁵N) have been analyzed. Calculation approaches based on the gas-phase and the Polarizable Continuum Model (PCM) approximations and their combinations with geometric constraints and additional noncovalent interactions have been probed. The main result of this work is that the effect of an aprotic polar solvent on the geometry of a proton-donor⋯H⋯proton-acceptor complex cannot be reproduced under the PCM approximation if no correction for solvent-solute interactions is made. These interactions can be implicitly accounted for using a simple computational protocol.

Refinement of labile hydrogen positions based on DFT calculations of 1H NMR chemical shifts: comparison with X-ray and neutron diffraction methods

Numerous gas phase electron diffraction, ultra-fast electron diffraction, X-ray and neutron diffraction experiments on β-dicarbonyl compounds exhibiting enol-enol tautomeric equilibrium, with emphasis on acetylacetone and dibenzoylmethane, have so far been reported with conflicting results on the structural details of the O-HO intramolecular hydrogen bond and resulted in alternative hypotheses on the intramolecular hydrogen bond potential function either a double minimum potential corresponding to two tautomeric forms in equilibrium or a single symmetrical one. We demonstrate herein, firstly, that the DFT calculated OH ¹H NMR chemical shifts of acetylacetone and dibenzoylmethane exhibit a strong linear dependence on the computed OO hydrogen bond length of ∼-50 ppm Å^-1 and as a function of the O-HO bond angle of ∼1 ppm per degree, upon the transfer of the hydrogen atom from the ground state toward the transition state. Secondly, the refinement of labile hydrogen atomic positions in intramolecular hydrogen bonds based on the root-mean-square deviation between experimentally determined and DFT calculated ¹H NMR chemical shifts in solution can provide high resolution structures of O-H and O(H)O bond lengths and O-HO bond angles with an accuracy of ∼10^-2 Å and ∼0.5°, respectively. Thirdly, the calculated ¹H NMR chemical shifts in solution of the two ground state tautomers in equilibrium of acetylacetone and dibenzoylmethane are in excellent agreement with the experimental value, even for moderate basis sets for energy minimization. In contrast, the single symmetrical structure in a strongly delocalized system is a transition state with calculated ¹H NMR chemical shifts which strongly deviate from the experimental value. Fourth, the DFT calculated ground state O-H bond lengths of acetylacetone and dibenzoylmethane are in quantitative agreement with the literature data which take into account the effect of quantum nuclear motion. The DFT structural results are critically discussed with respect to the state-of-the-art variable temperature X-ray and neutron diffraction methods.

Polar solvent fluctuations drive proton transfer in hydrogen bonded complexes of carboxylic acid with pyridines: NMR, IR and ab initio MD study

We study a series of intermolecular hydrogen-bonded 1 : 1 complexes formed by chloroacetic acid with 19 substituted pyridines and one aliphatic amine dissolved in CD₂Cl₂ at low temperature by ¹H and ¹³C NMR and FTIR spectroscopy. The hydrogen bond geometries in these complexes vary from molecular (O-HN) to zwitterionic (O^-H-N⁺) ones, while NMR spectra show the formation of short strong hydrogen bonds in intermediate cases. Analysis of C[double bond, length as m-dash]O stretching and asymmetric CO₂^- stretching bands in FTIR spectra reveal the presence of proton tautomerism. On the basis of these data, we construct the overall proton transfer pathway. In addition to that, we also study by use of ab initio molecular dynamics the complex formed by chloroacetic acid with 2-methylpyridine, surrounded by 71 CD₂Cl₂ molecules, revealing a dual-maximum distribution of hydrogen bond geometries in solution. The analysis of the calculated trajectory shows that the proton jumps between molecular and zwitterionic forms are indeed driven by dipole-dipole solvent-solute interactions, but the primary cause of the jumps is the formation/breaking of weak CHO bonds from solvent molecules to oxygen atoms of the carboxylate group.

Hydration of the Carboxylate Group in Anti-Inflammatory Drugs: ATR-IR and Computational Studies of Aqueous Solution of Sodium Diclofenac

Diclofenac (active ingredient of Voltaren) has a significant, multifaceted role in medicine, pharmacy, and biochemistry. Its physical properties and impact on biomolecular structures still attract essential scientific interest. However, its interaction with water has not been described yet at the molecular level. In the present study, we shed light on the interaction between the steric hindrance (the intramolecular N-H···O bond, etc.) carboxylate group (-CO₂^-) with water. Aqueous solution of sodium declofenac is investigated using attenuated total reflection-infrared (ATR-IR) and computational approaches, i.e., classical molecular dynamics (MD) simulations and density functional theory (DFT). Our coupled classical MD simulations, DFT calculations, and ATR-IR spectroscopy results indicated that the -CO₂^- group of the diclofenac anion undergoes strong specific interactions with the water molecules. The combined experimental and theoretical techniques provide significant insights into the spectroscopic manifestation of these interactions and the structure of the hydration shell of the -CO₂^- group. Moreover, the developed methodology for the theoretical analysis of the ATR-IR spectrum could serve as a template for the future IR/Raman studies of the strong interaction between the steric hindrance -CO₂^- group of bioactive molecules with the water molecules in dilute aqueous solutions.

Ab-Initio Molecular Dynamics Simulations and Calculations of Spectroscopic Parameters in Hydrogen-Bonding Liquids in Confinement (Project 8)

We investigate the effect of several nanoscale confinements on structural and dynamical properties of liquid water and binary aqueous mixtures. By means of molecular dynamics simulations based on density functional theory and atomistic force fields. Our main focus is on the dependence on the structure and the hydrogen-bonding-network of the liquids near the confinement interface at atomistic resolution. As a complementary aspect, spatially resolved profiles of the proton NMR chemical shift values are used to quantify the local strength of the hydrogen-bond-network.

Hydrogen Atomic Positions of O-H···O Hydrogen Bonds in Solution and in the Solid State: The Synergy of Quantum Chemical Calculations with ¹H-NMR Chemical Shifts and X-ray Diffraction Methods

The exact knowledge of hydrogen atomic positions of O-H···O hydrogen bonds in solution and in the solid state has been a major challenge in structural and physical organic chemistry. The objective of this review article is to summarize recent developments in the refinement of labile hydrogen positions with the use of: (i) density functional theory (DFT) calculations after a structure has been determined by X-ray from single crystals or from powders; (ii) ¹H-NMR chemical shifts as constraints in DFT calculations, and (iii) use of root-mean-square deviation between experimentally determined and DFT calculated ¹H-NMR chemical shifts considering the great sensitivity of ¹H-NMR shielding to hydrogen bonding properties.

Symmetry and dynamics of FHF− anion in vacuum, in CD2Cl2 and in CCl4. Ab initio MD study of fluctuating solvent–solute hydrogen and halogen bonds

Asymmetric solvation of FHF⁻ by halogen- and hydrogen-bonding solvents breaks the symmetry of the anion.

Some Brief Notes on Theoretical and Experimental Investigations of Intramolecular Hydrogen Bonding

A review of selected literature data related to intramolecular hydrogen bonding in ortho-hydroxyaryl Schiff bases, ortho-hydroxyaryl ketones, ortho-hydroxyaryl amides, proton sponges and ortho-hydroxyaryl Mannich bases is presented. The paper reports on the application of experimental spectroscopic measurements (IR and NMR) and quantum-mechanical calculations for investigations of the proton transfer processes, the potential energy curves, tautomeric equilibrium, aromaticity etc. Finally, the equilibrium between the intra- and inter-molecular hydrogen bonds in amides is discussed.

Selective recognition of fluoride salts by vasarenes: a key role of a self-assembled in situ dimeric entity via an exceptionally short [O-H-O](-) H-bond

A self-assembled supramolecular dimeric entity via an exceptionally short (2.404 Å) and strong (22.9 kcal mol(-1)) [O-H-O](-) hydrogen bond is the key to the special reactivity of vasarenes with fluoride salts. Vasarene is a self-assembled, vase-shaped compound, obtained by the reaction between ninhydrin and phloroglucinol. Analogous compounds are prepared by replacing the phloroglucinol with other polyhydroxy aromatics. Vasarenes show special affinity towards compounds of the type M(+)F(-), where M being a large monovalent cation, producing ion-pair-vasarene adducts. The first step in the proposed mechanism is the dissociation of the M(+)F(-) salt releasing F(-) to the solution, which may provide an explanation as to why only MF salts, which include large monovalent cations, undergo this reaction. From a practical point of view, the ease of their preparation and their special affinity towards fluoride salts make vasarenes potential means for salt separation. The readily formed dimeric structure with the very short [O-H-O](-) negative charge-assisted H-bond (-CAHB) can also be further used as a model in theoretical studies of such systems and understanding their role in biological processes.