A quantum description of the proton movement in an idealized NHN+ bridge

A series of model calculations was done to analyze the delocalization of the proton in the linking hydrogen bond of the (Dih)(2)H(+) cation (Dih: 4,5-dihydro-1H-imidazole). Standard quantum chemical calculations (B3LYP/D95+(d,p)) predict a low barrier hydrogen bond (LBHB) and thereby a delocalized proton in the NHN(+) hydrogen bridge. Explicit quantum calculations on the proton indicate that the delocalization of the proton does not provide enough energy to stabilize a permanent LBHB. Additional Born-Oppenheimer Molecular Dynamics (BOMD) simulations indicate further that the proton is localized at either side of the NHN(+) bridge and that a central proton position is the result of temporal averaging. The possibility of the proton to tunnel from one side to the other side of the NHN(+) bridge increases with the temperature as the trajectory of the (Dih)(2)H(+) cation runs through regions where the thermal excitation of Dih ring vibrations creates equal bonding opportunities for the proton on both sides of the bridge (vibrationally assisted proton tunneling). The quantum calculations for the proton in (Dih)(2)H(+) suggest further a broad peak for the 1 ← 0 transition with a maximum at 938 cm(-1) similar to that observed for LBHBs. Moreover, the asymmetric NHN(+) bridge in a thermally fluctuating environment is strong enough to create a significant peak at 1828 cm(-1) for the 2 ← 0 transition, while contributions from the 2 ← 1 are expected to be weak for the same reason.     

Survey of the proton affinities of adenine, cytosine, thymine and uracil dideoxyribonucleosides, deoxyribonucleosides and ribonucleosides

The kinetic method was applied to the determination of the proton affinities (PAs) of modified deoxy- and dideoxyribonucleosides. A correlation between the measured PAs and the replacement of one of the three hydroxyl groups of the ribose unit is presented. A PA scale was obtained which shows that the replacement of the primary or of one or both secondary hydroxyl groups of a ribonucleoside with a hydrogen atom induces the lowering or the enhancement of the nucleoside PA, respectively. The scale extends over a very narrow range of approximately 2 kcal mol(-1), thus demonstrating the sensitivity of the kinetic method in the evaluation of small differences in thermodynamic parameters.     

Phonon driven proton transfer in crystals with short strong hydrogen bonds

Recent work on understanding why protons migrate with increasing temperature in short, strong hydrogen bonds is extended here to three more organic, crystalline systems. Inelastic neutron scattering and density functional theory based simulations are used to investigate structure, vibrations, and dynamics of these systems as functions of temperature. The mechanism determined in a previous work on urea phosphoric acid of low frequency vibrations stabilizing average crystal structures, in which the potential energy well of the hydrogen bond has its minimum shifted towards the center of the bond, is found to be valid here. The new feature of the N-H...O hydrogen bonds studied in this work is that the proton is transferred from the donor atom to the acceptor atom. Molecular dynamics simulations show that in an intermediate temperature regime, in which the proton is not completely transferred, the proton is bistable, jumping from one side of the hydrogen bond to the other. In the case of 3,5-pyridine dicarboxylic acid, which has been studied in most detail, specific phonons are identified, which influence the potential energy surface of the proton in the short, strong hydrogen bond.     

The coordination chemistry of the proton

It is well known that an acidic hydrogen atom can form hydrogen bonds to a hydrogen bond acceptor, a Lewis base. It is considerably less known that the proton can coordinate two or more atoms conveniently in bonding modes that cannot be described as hydrogen bonding. Agostic interactions, bridging hydrides, 3-centre-2-electron bonds in boranes, bifurcated hydrogen atoms, they are all elements of the coordination chemistry of the proton and, of course, the hydrogen bond comes in more than one facette as well.     

Hydrogen-bonding partner of the proton-conducting histidine in the influenza m2 proton channel revealed from (1)h chemical shifts

The influenza M2 protein conducts protons through a critical histidine (His) residue, His37. Whether His37 only interacts with water to relay protons into the virion or whether a low-barrier hydrogen bond (LBHB) also exists between the histidines to stabilize charges before proton conduction is actively debated. To address this question, we have measured the imidazole (1)H(N) chemical shifts of His37 at different temperatures and pH using 2D (15)N-(1)H correlation solid-state NMR. At low temperature, the H(N) chemical shifts are 8-15 ppm at all pH values, indicating that the His37 side chain forms conventional hydrogen bonds (H-bonds) instead of LBHBs. At ambient temperature, the dynamically averaged H(N) chemical shifts are 4.8 ppm, indicating that the H-bonding partner of the imidazole is water instead of another histidine in the tetrameric channel. These data show that His37 forms H-bonds only to water, with regular strength, thus supporting the His-water proton exchange model and ruling out the low-barrier H-bonded dimer model.     

Theoretical study of the symmetry of the (OH...O)- hydrogen bonds in vinyl alcohol-vinyl alcoholate systems

The interactions between substituted vinyl alcohols and vinyl alcoholates (X = NH(2), H, F, Cl, CN) are studied at the B3LYP/6-311++G(d,p) level of theory. In a first step, the conformation of the monomers is investigated and the proton affinities (PA(A(-))) of the enolates are calculated. The enols and enolates are held together by strong (OH...O)(-) hydrogen bonds, the hydrogen bond energies ranging from 19.1 to 34.6 kcal mol(-1). The optimized O...O distances are between 2.414 and 2.549 A and the corresponding OH distances from 1.134 and 1.023 A. The other geometry parameters such as C[double bond]C or CO distances also indicate that, in the minimum energy configuration, the hydrogen bonds are characterized by a double well potential. The Mulliken charges on the different atoms of the proton donors and proton acceptors and the frequencies of the nu(OH) stretching vibrations agree with this statement. All the data indicate that the hydrogen bonds are the strongest in the homomolecular complexes. The transition state for hydrogen transfer is located with the transition barrier estimated to be about zero. Upon addition of the zero-point vibration energies to the total potential energy, the barrier vanishes. This is a characteristic feature of low-barrier hydrogen bonds (LBHBs). The hydrogen bond energies are correlated to the difference 1.5 PA(AH) - PA(A(-)). The correlation predicts different energies for homomolecular hydrogen bonds, in agreement with the theoretical calculations. Our results suggest that a PA (or pK(a)) match is not a necessary condition for forming LBHBs in agreement with recent data on the intramolecular hydrogen bond in the enol form of benzoylacetone (J. Am. Chem. Soc. 1998, 120, 12117).     

Interstitial water and the formation of low barrier hydrogen bonds: A computational model study

The B3LYP/D95+(d,p) analysis of the uncharged low barrier hydrogen bond (LBHB) between 4‐methyl‐1H‐imidazole (Mim) and acetic acid (HAc) shows that uncharged LBHBs can be formed either by adding three water molecules around the cluster or by placing the Mim–HAc pair in a dielectric environment created by a polarizable continuum model with a permittivity larger than 20.7. The permittivity of environment around uncharged LBHB can be lowered significantly by including water molecules into the system. A Mim–HAc LBHB stabilized with one water molecule observed in diethyl ether (ε = 4.34), with two water molecules in toluene (ε = 2.38), and with three water molecules in vacuo (ε = 1). Solvation models with different numbers of water molecules predict average differences in the proton affinities of the hydrogen bonded bases (ΔPA) for stable uncharged LBHB systems in vacuo to be 91.5 kcal/mol being different from the ΔPA values close to zero in charge‐assisted LBHB systems. The results clearly indicate that small amounts of interstitial water molecules at the active site of enzymes do not preclude the existence of LBHBs in biological catalysis. Our results also show that interstitial water molecules provide a useful clue in the search for uncharged LBHBs in an enzymatic environment and the number of water molecules can be used as a relative measure for the polarity around the direct environment of LBHBs. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

Vibrations of a chelated proton in a protonated tertiary diamine

Vibrational spectra of the conjugate acid of Me(2)NCH(2)CH(2)CH(2)CH(2)NMe(2) (N,N,N',N'-tetramethylputrescine) have been examined in the gaseous and crystalline phases using Infrared Multiple Photon Dissociation (IRMPD) spectroscopy, Inelastic Neutron Scattering (INS), and high pressure Raman spectroscopy. A band observed near 530 cm(-1) is assigned to the asymmetric stretch of the bridging proton between the two nitrogens, based on deuterium substitution and pressure dependence. The NN distance measured by X-ray crystallography gives a good match to DFT calculations, and the experimental band position agrees with the value predicted from theory using a 2-dimensional potential energy surface. The reduced dimensionality potential energy surface, which treats the ion as though it possesses a linear NHN geometry, shows low barriers to proton transit from one nitrogen to the other, with zero point levels close to the barrier tops. In contrast, two other related systems containing strong hydrogen bonds do not exhibit the same spectroscopic signature of a low barrier hydrogen bond (LBHB). On the one hand, the IRMPD spectra of the conjugate acid ions of the amino acid N,N,N',N'-tetramethylornithine (in which the two nitrogens have different basicities) show fewer bands and no comparable isotopic shifts in the low frequency domain. On the other hand, the IRMPD spectrum of the shorter homologue Me(2)NCH(2)CH(2)CH(2)NMe(2) (N,N,N',N'-tetramethyl-1,3-propanediamine), for which the NHN bond angle deviates substantially from linearity, displays more than one band in the 1100-1400 cm(-1) domain, which vanish as a consequence of deuteration.     

Electronic basis of the comparable hydrogen bond properties of small H2CO/(H2O)n and H2NO/(H2O)n systems (n = 1, 2)

The electronic and structural properties of dihydronitroxide/water clusters are investigated and compared to the properties of formaldehyde/water clusters. Exploring the stationary points of their potential energy surfaces (structurally, vibrationally, and energetically) and characterizing their hydrogen bonds (by both atoms in molecules and natural bond orbitals methods) clearly reveal the strong similarity between these two kind of molecular systems. The main difference involves the nature of the hydrogen bond taking place between the X-H bond and the oxygen atom of a water molecule. All the properties of the hydrogen bonds occurring in both kind of clusters can be easily interpreted in terms of competition between intermolecular and intramolecular hyperconjugative interactions.     

Hydrogen bonding nature between calix[6]arene and piperidine/triethylamine

We investigated the binding nature of the 1,2,3-alternate calix[6]arene with one piperidine, two piperidines, and two triethyl amines with a special emphasis on the hydrogen bonding networks by density functional theory calculations. The 1,2,3-alternate calix[6]arene strongly binds with piperidines and triethylamines at two different binding sites, exo and endo sites. In the two binding sites, the hydrogen bonding nature shows a characteristic difference. In the exo site, there formed only one hydrogen bond, while in the endo site, two hydrogen bonds except for the triethylamine. The proton transfer within the hydrogen bonding and the hydrogen bonding types, normal hydrogen bonding (NHB), short strong hydrogen bond (SSHB), and low barrier hydrogen bonding (LBHB), will be discussed in detail.     

The Role of Hydrogen Bond in Catalytic Triad of Serine Proteases

In order to investigate the origin of catalytic power for serine proteases, the role of the hydrogen bond in the catalytic triad was studied in the proteolysis process of the peptides chymotrypsin inhibitor 2 (CI2), MCTI-A, and a hexapeptide (SUB), respectively. We first calculated the free energy profile of the proton transfer between His and Asp residues of the catalytic triad in the enzyme-substrate state and transition state by employing QM/MM molecular dynamics simulations. The results show that a low-barrier hydrogen bond (LBHB) only forms in the transition state of the acylation of CI2, while it is a normal hydrogen bond in the acylation of MCTI-A or SUB. In addition, the change of the hydrogen bond strength is much larger in CI2 and SUB systems than in MCTI-A system, which decreases the acylation energy barrier significantly for CI2 and SUB. Clearly, a LBHB formed in the transition state region helps accelerate the acylation reaction. But to our surprise, a normal hydrogen bond can also help to decrease the energy barrier. The key to reducing the reaction barrier is the increment of hydrogen bond strength in the transition state state, whether it is a LBHB or not. Our studies cast new light on the role of the hydrogen bond in the catalytic triad, and help to understand the catalytic triad of serine proteases.     

On the Possibility of Using UV Spectroscopy as a Measure of the Low-Barrier Hydrogen Bond

In many enzyme-catalyzed biochemical pathways, a short, strong hydrogen bond between an enzyme and substrate is an important structural feature. These bonds are termed low-barrier hydrogen bonds. In this paper, we show that UV spectra can be used as an experimental technique to determine if a system contains a low-barrier hydrogen bond (LBHB). We simulate, using the time-dependent view of UV spectroscopy, several different UV spectra: absorption, photodissociation, and emission, on systems containing a low-barrier hydrogen bond. We find several distinguishing spectral features in these UV spectra for systems that possess a LBHB.     

Structures,vibrational frequencies,topologies, and energies of hydrogen bonds in cysteine‐formaldehyde complexes

The hydrogen bonding interactions between cysteine (Cys) and formaldehyde (FA) were studied with density functional theory regarding their geometries, energies, vibrational frequencies, and topological features of the electron density. The quantum theory of atoms in molecules and natural bond orbital analyses were employed to elucidate the interaction characteristics in the Cys‐FA complexes. The intramolecular hydrogen bonds (H‐bonds) formed between the hydroxyl and the N atom of cysteine moiety in some Cys‐FA complexes were strengthened because of the cooperativity. Most of intermolecular H‐bonds involve the O atom of cysteine/FA moiety as proton acceptors, while the strongest H‐bond involves the O atom of FA moiety as proton acceptor, which indicates that FA would rather accept proton than providing one. The H‐bonds formed between the CH group of FA and the S atom of cysteine in some complexes are so weak that no hydrogen bonding interactions exist among them. In most of complexes, the orbital interaction of H‐bond is predominant during the formation of complex. The electron density (ρ_b) and its Laplace (?²ρ_b) at the bond critical point significantly correlate with the H‐bond parameter δR, while a linearly relationship between the second‐perturbation energy E(2) and ρ_b has been found as well. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

QTAIM study of strong H-bonds with the O-H...a fragment (A=O, N) in three-dimensional periodical crystals

The relationship between the d(H...A) distance (A=O, N) and the topological properties at the H...A bond critical point of 37 strong (short) hydrogen bonds occurring in 26 molecular crystals are analyzed using the quantum theory of atoms in molecules (QTAIM). Ground-state wave functions of the three-dimensional periodical structures representing the accurate experimental geometries calculated at the B3LYP/6-31G** level of approximation were used to obtain the QTAIM electron density characteristics. The use of an electron-correlated method allowed us to reach the quantitatively correct values of electron density rhob at the H...A bond critical point. However, quite significant differences can appear for small absolute values of the Laplacian (<0.5 au). The difference between the H...O and H...N interactions is described using the rhob versus d(H...A) dependence. It is demonstrated that the values of parameters in this dependence are defined by the nature of the heavy atom forming the H...A bond. An intermediate (or transit) region separating the shared and closed-shell interactions is observed for the H-bonded crystals in which the bridging proton can move from one heavy atom to another. The crystalline environment changes the location of the bridging proton in strong H-bonded systems; however, the d(O-H)/d(H...O) ratio is approximately the same for both the gas-phase complexes and molecular crystals with a linear or near-linear O-H...O bond.     

Substituent effects on two-center three-electron bonds and hydrogen bonds involving unsaturated organic functional groups and an ammonia radical cation--the resonance contribution

A theoretical investigation of the substituent effects on the two-center, three-electron (2c-3e) bond involved between unsaturated functional groups and an amine nitrogen is presented. The competitive hydrogen-bonded complexes are also studied. In both cases, the bond energies are found to be in the range of 20-30 kcal mol(-1). The variation of these energies is discussed with respect to the electron-donating effect of the substitutents, as well as with respect to the alpha-bonded atom of the organic functional group (O, S, NH). For the 2c-3e bonds, the results are discussed on the basis of the differences of the ionization potential (delta IP) of the separated fragments and can be rationalized through the valence bond theory. For the hydrogen bonds, the substituent influence is discussed by using the differences of the proton affinities (delta PA) of the substrates. The resonating nature of the hydrogen bond in these cationic systems is investigated and is found to account for most of the binding energy. Marcus theory is compared with the proposed resonating model.     

Theoretical studies of the proton transfer behaviors in molecular complexes analogous to catalytic triad of serine protease: Toward understanding the existence and significance of the low-barrier hydrogen-bond in enzymatic catalysis

A representative acetate-(5-methylimidazole)-methanol system has been employed as a model of catalytic triad in serine protease to validate the formation processes of low-barrier H-bonds (LBHB) at the B3LYP/6-311++G** level of theory, and variable H-bonding characters from conventional ones to LBHBs have been represented along with the proceedings of proton transfer. Solvent effect is an important factor in modulation of the existence of an LBHB, where an LBHB (or a conventional H-bond) in the gas phase can be changed into a non-LBHB (an LBHB) upon solvation. The origin of the additional stabilization energy arising from the LBHB may be attributed to the H-bonding energy difference before and after proton transfer because the shared proton can freely move between the proton donor and proton acceptor. Most importantly, the order of magnitude of the stabilization energy depends on the studied systems. Furthermore, the nonexistence of LBHBs in the catalytic triad of serine proteases has been verified in a more sophisticated model treated using the ONIOM method. As a result, only the single proton transfer mechanism in the catalytic triad has been confirmed and the origin of the powerful catalytic efficiency of serine proteases should be attributed to other factors rather than the LBHB. Supported by the National Natural Science Foundation of China (Grant Nos. 20633060 & 20573063), the Natural Science Foundation of Shandong Province (Grant No. Y2007B23), the Scientific Research Foundation of Qufu Normal University (Grant Nos. Bsqd2007003 and Bsqd2007008), and the State Key Laboratory of Physical Chemistry of Solid Surfaces (Xiamen University)     

Theoretical study on the bromomethane-water 1:2 complexes

Bromomethane-water 1:2 complexes have been theoretically studied to reveal the role of hydrogen bond and halogen bond in the formation of different aggregations. Four stable structures exist on the potential energy surface of the CH3Br(H2O)2 complex. The bromine atom acts mainly as proton acceptor in the four studied structures. It is also capable of participating in the formation of the halogen bond. The properties and characteristics of the hydrogen bond and the halogen bond are investigated employing several different quantum chemical analysis methods. Cooperative effects for the pure hydrogen bonds or the mixed hydrogen bonds with halogen bonds and the possibility of describing cooperative effects in terms of the topological analysis of the electronic density or the charge-transfer stabilization energy are discussed in detail. An atoms-in-molecules study of the hydrogen bond or the halogen bond in the bromomethane-water 1:2 complexes suggests that the electronic density topology of the hydrogen bond or the halogen bond is insensitive to the cooperative effect. The charge-transfer stabilization energy is proportional to the cooperative effect, which indicates the donor-acceptor electron density transfer to be mainly responsible for the trimer nonadditive effect.     

Theoretical study of (CH...C)- hydrogen bonds in CH(4-n)X(n) (X = F, Cl; n = 0, 1, 2) systems complexed with their homoconjugate and heteroconjugate carbanions

This work deals with a theoretical study of the (CH...C)- hydrogen bonds in CH4, CH3X, and CH2X2 (X = F, Cl) complexed with their homoconjugate and heteroconjugate carbanions. The properties of the complexes are calculated with the B3LYP method using the 6-311++G(d,p) or 6-311++G(2df,2p) basis sets. The deprotonation enthalpies (DPE) of the CH bond or the proton affinities of the carbanions (PA(C-) are calculated as well. All the systems with the exception of the CH4...CHCl2(-) one are characterized by a double minimum potential. In some of the complexes, the (CH(b)...C)- hydrogen bond is linear. In other systems, such as CH3F...CH2F- and CH3F...CHF2(-), there is a large departure from linearity, the systems being stabilized by electrostatic interactions between the nonbonded H of the neutral molecule and the F atom of the carbanion. In the transition state, the (CH(b)...C)- bond is linear, and there is a large contraction of the intermolecular C...C distance. The binding energies vary within a large range, from -1.4 to -11.1 kcal mol(-1) for the stable complexes and -8.6 to -44.1 kcal mol(-1) for the metastable complexes. The energy barriers to proton transfer are between 5 and 20 kcal mol(-1) for the heteroconjugate systems and between 3.8 and 8.3 kcal mol(-1) for the homoconjugate systems. The binding energies of the linear complexes depend exponentially on 1.5DPE - PA(C-), showing that the proton donor is more important than the proton acceptor in determining hydrogen bond strength. The NBO analysis indicates an important electronic reorganization in the two partners. The elongations of the CH bond resulting from the interaction with the carbanion depend on the occupation of the sigma*(CH(b)) antibonding orbitals and on the hybridization of the C bonded to H(b). The frequency shifts of the nu(CH)(A1) stretching vibration range between 15 and 1150 cm(-1). They are linearly correlated to the elongation of the CH(b) bond.     

Hydrogen-bond mediated catalysis: the aminolysis of 6-chloropyrimidine as catalyzed by derivatives of uracil.

The aminolysis of 6-chloropyrimidine and 2-amino-6-chloropyrimidine has been examined by using density functional theory. Relative to the aminolysis of 6-chloropyrimidine, the addition of an electron-donating NH(2) group to C(2) increases the barrier to aminolysis, indicating that the third hydrogen bond does not play a catalytic role but introduces additional rigidity into the system. However, the computations suggest that there is an interesting correlation between the barrier to aminolysis and the proton affinity of the species that interacts with the incoming NH(3). To extend the range of proton affinities, the aminolysis of 6-chloropyrimidine was examined by using fluoro, imine, and thioketo derivatives of the uracil-derived bases. The proton affinity of the moiety that hydrogen bonds with NH(3) is decreased by fluoro substitution, and thus the aminolysis barriers are increased. Similarly, imine substitution enhances the PA of the moiety, which is reflected in a decrease in the aminolysis barriers. The same correlation exists for the thioketo-derived bases, whose PAs are intermediate between the fluoro and imine derivatives. Thus, the aminolysis of 6-chloropryimidine and 2-amino-6-chloropyrimidine demonstrates the importance of a well-chosen proton acceptor and the catalytic possibilities associated with the formation of multiple hydrogen bonds.     

Special pair dance and partner selection: elementary steps in proton transport in liquid water

Conditional and time-dependent radial distribution functions reveal the details of the water structure surrounding the hydronium during the proton mobility process. Using this methodology for classical multistate empirical valence bond (MS-EVB) and ab initio molecular dynamics trajectories, as well as quantal MS-EVB trajectories, we supply statistical proof that proton hops in liquid water occur by a transition from the H3O+[3H2O] Eigen-complex, via the H5O2+ Zundel-complex, to a H3O+[3H2O] centered on a neighboring water molecule. In the "resting period" before a transition, there is a distorted hydronium with one of its water ligands at a shorter distance and another at a longer distance than average. The identity of this "special partner" interchanges rapidly within the three first-shell water ligands. This is coupled to cleavage of an acceptor-type hydrogen bond. Just before the transition, a partner is selected by an additional translation of the H3O+ moiety in its direction, possibly enabled by loosening of donor-type hydrogen bonds on the opposite side. We monitor the transition in real time, showing how the average structure is converted to a distorted H5O2+ cation constituting the transitional complex for proton hopping between water molecules.