Calculated vibrational spectra of weak hydrogen bonds

Potential energy and dipole moment surfaces for the H-bond SH... S in the dimeric methanethiol have been calculated by the SCF-MO-LCGO method, and the vibrational spectrum — transition frequencies and IR absorption intensities at 20 °K — computed. This spectrum is compared with that of the monomeric species and with experimental results. The resulting dimerization energy is 1.4 kcal/mole.     

Comparative study of weak interactions in molecular crystals: H-H bonds vs hydrogen bonds

The crystal structures of tetraphenylphosphonium squarate, bianthrone, and bis(benzophenone)azine are shown to contain a variety of C-H(delta+)...(delta+)H-C interactions, as well as a variety of C-H...O and C-H...C(pi) interactions. Each of these molecules possesses interactions that can possibly be characterized as either H-H bonds or weak hydrogen bonds based on the first four criteria proposed by Koch and Popelier. These interactions have been completely characterized topologically after the multipole refinement of the structures. It appears that weak interactions of the form C-H(delta+)...(delta+)H-C possess certain correlations between the various properties of the electron density at the bond critical points. The coexistence of the three types of interactions makes it possible to establish similarities and differences in the correlations of these weak interactions. This all leads to a better understanding of H-H interactions and how they fit into the hierarchy of weak interactions.     

The role of weak hydrogen bonds in chiral recognition

Chiral recognition has been studied in neutral or ionic weakly bound complexes isolated in the gas phase by combining laser spectroscopy and quantum chemical calculations. Neutral complexes of the two enantiomers of lactic ester derivatives with chiral chromophores have been formed in a supersonic expansion. Their structure has been elucidated by means of IR-UV double resonance spectroscopy in the 3 μm region. In both systems described here, the main interaction ensuring the cohesion of the complex is a strong hydrogen bond between the chromophore and methyl-lactate. However, an additional hydrogen bond of much weaker strength plays a discriminative role between the two enantiomers. For example, the 1:1 heterochiral complex between R-(+)-2-naphthyl-ethanol and S-(+) methyl-lactate is observed, in contrast with the 1:1 homochiral complex which lacks this additional hydrogen bond. On the other hand, the same kind of insertion structures is formed for the complex between S-(±)-cis-1-amino-indan-2-ol and the two enantiomers of methyl-lactate, but an additional addition complex is formed for R-methyl-lactate only. This selectivity rests on the formation of a weak CHπ interaction which is not possible for the other enantiomer. The protonated dimers of Cinchona alkaloids, namely quinine, quinidine, cinchonine and cinchonidine, have been isolated in an ion trap and studied by IRMPD spectroscopy in the region of the ν(OH) and ν(NH) stretch modes. The protonation site is located on the alkaloid nitrogen which acts as a strong hydrogen bond donor in all the dimers studied. While the nature of the intermolecular hydrogen bond is similar in the homochiral and heterochiral complexes, the heterochiral complex displays an additional weak CHO hydrogen bond located on its neutral part, which results in slightly different spectroscopic fingerprints in the ν(OH) stretch region. This first spectroscopic evidence of chiral recognition in protonated dimers opens the way to the study of the complexes of Cinchona alkaloids involved in enantioselective catalysis. These examples show how secondary hydrogen bonds controlled by stereochemical factors govern molecular recognition processes.     

Cation-anion hydrogen bonds: a new class of hydrogen bonds that extends their strength beyond the covalent limit. A theoretical characterization

The existence of O-H···O hydrogen bonds having a strength within the -80 to -210 kcal/mol range, that is, in the range of strength of covalent bonds and well beyond the so-called covalent limit (-50 kcal/mol), is reported on complexes where the O-H proton donor and O acceptor groups are located in ions of opposite sign. A complete analysis of short distance O-H···O hydrogen bonds between charged fragments was performed for cases where the OH and O groups are both located on charged molecules. It shows that these interactions (a) are nonsymmetrical for the O-H and H···O distances, (b) have a noncovalent H···O bond critical point, and (c) have a strong and energetically stable electrostatic component when the OH and O groups are located in oppositely charged molecules. These cation-anion O-H···O interactions are energetically stable, satisfy the usual topology for hydrogen bonds, HBs, and also have the same directionality found in other HBs. Therefore, they should be considered as a new class of HBs, the cation-anion hydrogen bonds.     

Characterizing the strength of individual hydrogen bonds in DNA base pairs

The energies of individual hydrogen bonds (H-bonds) in A-T and G-C Watson-Crick base pairs were calculated according to the natural bond orbital (NBO) analysis of intermolecular interactions. The extent to which individual H-bonds are helpful in holding the two base pairs together was previously investigated quantitatively by a few different approaches, and the results of the present and previous estimations were compared. The method was validated by the determination of the H-bond strength changes in A-T and G-C pairs upon the substitution of the monomer (base) by two cationic substituents; the systems for which the changes were previously anticipated based on the modifications of the H-bonds' distances.     

Investigation of single-molecule kinetics mediated by weak hydrogen bonds within a biological nanopore

The study of factors essential for protein-peptide interactions and protein pore-mediated peptide transport are of particular relevance in biology. Wild-type α-hemolysin was adopted as a "nanoreactor" in which perturbations of the current through a protein containing a lumen-residing, aryl-capped antimicrobial peptide were seen for the first time and studied at the single-molecule level. Energy and steric considerations hint that Met-aryl interactions between aromatic residues placed at a peptide's extremities and any of the methionines lining the α-hemolysin constriction region may be the primary cause of peptide stabilization within the lumen and may be particularly important to the peptide-α-hemolysin interaction.     

Comparison of the performance of dispersion-corrected density functional theory for weak hydrogen bonds

Potential energy curves for five complexes with weak to medium strong hydrogen bonds have been computed with dispersion corrected DFT methods. The electronic density based vdW-DF2 and VV10 van der Waals density functionals have been tested, as well as an atom pair-wise correction method (DFT-D3). The short-range exchange-correlation components BLYP and rPW86-PBE together with the extended aug-cc-pVQZ basis sets have been employed. Reference data have been computed at the estimated CCSD(T)/CBS(aQ-a5) level of theory. The investigated systems are CH(4)·NH(3), Cl(3)CH·NH(3), NH(3)·NH(3), CH(3)F·C(2)H(2) and CH(3)F·H(2)O with binding energies ranging from -0.7 kcal mol(-1) to -5.5 kcal mol(-1). We find that all dispersion corrected methods perform reasonably well for these hydrogen bonds, but also observe distinct differences. The BLYP-D3 method provides the best results for three out of five systems. For the fluorinated complexes, the VV10 method gives remarkably good results. The vdW-DF2 method yields good interaction energies similar to the other methods (mean average deviation of 0.2-0.3 kcal mol(-1)), but fails to provide accurate equilibrium separations. Based on these results and previous experience with the computation of non-covalent interactions, for large-scale applications we can recommend DFT-D3 based structure optimizations with subsequent checking of interaction energies by single-point VV10 computations. Comparison of the DFT-D3 and VV10 results leads to the conclusion that the short-range exchange-correlation functional and not the dispersion correction mainly determines the achievable accuracy.     

A new method for quick predicting the strength of intermolecular hydrogen bonds

A new method is proposed to quick predict the strength of intermolecular hydrogen bonds. The method is employed to produce the hydrogen-bonding potential energy curves of twenty-nine hydrogen-bonded dimers. The calculation results show that the hydrogen-bonding potential energy curves obtained from this method are in good agreement with those obtained from MP2/6-31+G** calculations by including the BSSE correction, which demonstrate that the method proposed in this work can be used to calculate the hydrogen-bonding interactions in peptides. Supported by the National Natural Science Foundation of China (Grants Nos. 20573049 and 20633050) and the research fund of the Educational Department of Liaoning Province (2004C019, 20060469)     

Estimation on the individual hydrogen-bond strength in molecules with multiple hydrogen bonds

A simple atom-replacement approach is proposed for estimating the individual contributions of each intermolecular hydrogen bond (HB) in multiple hydrogen-bonded systems. The approach is validated by calculations on the homodimer of formylformamide and then applied to nucleic acid base pairs (adenine-thymine and guanine-cytosine) and some quadruply hydrogen-bonded dimers. With the help of this method, it is easy to distinguish the relative strength of each HB, and identify the main factors contributing to the total binding energies of multiple HBs.     

Azabicyclo[3.3.1]nonanone: a case when weak interactions are preferred over strong hydrogen bonds

Derivatives of azabicyclo[3.3.1]nonanone tend to prefer for weak interactions in the crystal over strong N–H···O hydrogen bonds. The main stabilizing forces in the investigated azatricyclo[7.3.1.0^2,7]trideca-trienone derivatives are C–H···O, N–H···π and C–H···π interactions, leading to interesting structural patterns. The azabicyclo[3.3.1]nonanone ring adopts chair-envelope conformation having exo-C2,C4-aromatic substituents. Amino NH is in trigonal pyramidal configuration. The interesting stereochemistry of azabicyclo[3.3.1]nonanone, driving exceptional preference for weaker interactions over strong hydrogen bonds serves a useful example toward engineering and design strategy, and structure prediction methodologies.     

Enhancement of the solubilities of polycyclic aromatic hydrocarbons by weak hydrogen bonds with water

Summary The thermodynamics of mobile order is applied to predict the aqueous solubility of liquid and solid aliphatic and polycyclic aromatic hydrocarbons. The solubility values are mainly determined by the magnitude of the hydrophobic effect. However, contrary to the solubilities of the alkanes, the solubilities of polycyclic aromatic hydrocarbons in water predicted in absence of solute-solvent hydrogen (H) bonds are systematically too low. This shows the contribution of weak specific interactions between the OH groups and the electrons of the aromatic substances. According to the theory, these interactions are characterized by a stability bility constant K_o which can be derived from solubility data. At 25°C, this constant amounts to 80 cm³/mol, the order of magnitude of which can be explained by the competition of these intermolecular bonds with the rather weak self-association bonds in the secondary chains of water.     

Non-conventional hydrogen bonds: pterins-metal anions

In this paper, we present an analysis of the interaction of metal ions (Cu, Ag and Au) with three different pterins (pterin, isoxanthopterin and sepiapterin) to provide insights concerning the formation of conventional and non-conventional H bonds. Density functional theory calculations were performed in order to reveal the optimized structures of pterin molecules, dimers and tetramers compounds, both with and without metal anions (M). The interaction with small metal clusters (M(3)) is also considered. The formation of different systems is characterized in terms of the structural parameters and hydrogen binding energies (HBE). The HBE values for pterin-M systems presented in this study lie between 22 and 60 kcal mol(-1) and can therefore be classified as strong conventional and strong non-conventional hydrogen bonds. The HBE with small metal clusters (pterin-M(3)) are smaller than the HBE with metal atoms. Vertical electron detachment energies (VEDEs) are also reported in order to analyze the influence of the hydrogen bond on electronic properties. A direct correlation between VEDEs and HBE was found for pterin-M and pterin-M(3) complexes; i.e. as the VEDEs increase, the HBE also augment. The only exception is with Ag(3). The main conclusion derived from this study is that the strong non-conventional hydrogen bonds formed between pterins, dimers and tetramers do not affect the formation of conventional hydrogen bonds between pterins but they do influence the VEDEs.     

Topological properties of the electrostatic potential in weak and moderate N...H hydrogen bonds

The topological analyses of the electrostatic potential phi(r) and the electron density distribution rho(r) have been performed for a set of 20 neutral complexes with weak and moderate N...H bonds. In all cases, a zero flux surface of the electrostatic potential containing a saddle point analogous to the bond critical point of the electron density distribution is observed. These surfaces define an equivalent of the atomic basin of rho(r) for the electrostatic potential, which exhibits zero net charge and can be regarded as an electrostatically isolated region if its volume is finite. The phi(r) and rho(r) zero flux surfaces divide the hydrogen-bonding region in three parts, being the central one related to the electrostatic interaction between donor and acceptor. This central region exhibits a relative size of approximately 13-14% of the N...H distance dNH, it belongs to the outermost shell of the nitrogen and is mainly associated with its lone pair. Topological properties of both rho(r) and phi(r), as well as the electron kinetic (G) and potential (V) energy densities, show similar dependences with dNH at both bond critical points (phi-BCP and rho-BCP). Phenomenological proportionalities between the rho(r) curvatures and G and V are also found at the electrostatic potential critical point. The curvatures of the electrostatic potential, which are interpreted in terms of the electrostatic forces in the bonding region, present the same exponential dependency as the electron density distribution, to which they are related by Poisson's equation.     

C-H...O and other weak hydrogen bonds. From crystal engineering to virtual screening

Hydrogen bonds, X-H...A, formed by weak donors (X = C) and acceptors (A = pi system) were generally dismissed as being of little consequence before and even during the 1970s. This situation changed in the early 1980s, and during the two following decades they were implicated as being significant in many small molecule crystal structures, and also in solution. Today, knowledge gained about these interactions is being used to understand the structure of biomolecules with implications for structure based drug design.     

Temperature-dependent intensity anomalies in amino acid esters: weak hydrogen bonds in protected glycine, alanine and valine

Esters of glycine, alanine and valine are investigated by FTIR and Raman spectroscopy in supersonic jets as gas phase model systems for the neutral peptide N-terminus. The NH-stretching vibrations exhibit very large temperature- and substitution-dependent intensity anomalies which are related to weak, bifurcated intramolecular hydrogen bonds to the carbonyl group. Comparison to theory is only satisfactory at low temperature. Spectral NH aggregation shifts are small or even negligible and the associated IR intensity is remarkably low. In the case of valine, chirality recognition effects are nevertheless detected and rationalized. Comparison to quantum-chemical calculations for dimers shows that dispersion interactions are essential. It also rules out cooperative hydrogen bond topologies and points at deficiencies in standard harmonic treatments with the linear dipole approximation.     

Influence of substituents on the strength of aryl C-H...anion hydrogen bonds

[structure: see text]When electron-withdrawing substituents are present, aryl C-H groups become powerful hydrogen bond donors, forming stronger complexes than obtained with conventional O-H and N-H groups.     

Some measures for making halogen bonds stronger than hydrogen bonds in H2CS-HOX (X = F, Cl, and Br) complexes

The properties and applications of halogen bonds are dependent greatly on their strength. In this paper, we suggested some measures for enhancing the strength of the halogen bond relative to the hydrogen bond in the H(2)CS-HOX (X = F, Cl, and Br) system by means of quantum chemical calculations. It has been shown that with comparison to H(2)CO, the S electron donor in H(2)CS results in a smaller difference in strength for the Cl halogen bond and the corresponding hydrogen bond, and the Br halogen bond is even stronger than the hydrogen bond. The Li atom in LiHCS and methyl group in MeHCS cause an increase in the strength of halogen bonding and hydrogen bonding, but the former makes the halogen bond stronger and the latter makes the hydrogen bond stronger. In solvents, the halogen bond in the Br system is strong enough to compete with the hydrogen bond. The interaction nature and properties in these complexes have been analyzed with the natural bond orbital theory.