Ab initio quantum mechanics analysis of imidazole C-H...O water hydrogen bonding and a molecular mechanics forcefield correction

Antagonistic Interplay Between an Intermolecular CH···O and an Intramolecular OH···O Hydrogen Bond in a 1:1 Complex Between 1,2-Cyclohexanedione and Chloroform: A Combined Matrix Isolation Infrared and Quantum Chemistry Study

Matrix isolation infrared spectra of a weak C-H···O hydrogen-bonded complex between the keto-enol form of 1,2-cyclohexanedione (HCHD) and chloroform have been measured. The spectra reveal that the intramolecular O-H···O H-bond of HCHD is weakened as a result of complex formation, manifesting in prominent blue shift (∼23 cm^-1) of the ν_O-H band and red shifts (∼7 cm^-1) of ν_C═O bands of the acceptor (HCHD). The ν_C-H band of donor CHCl₃ undergoes a large red shift of ∼33 cm^-1. Very similar spectral effects are also observed for formation of the complex in CCl₄ solution at room temperature. Our analysis reveals that out of several possible iso-energetic conformational forms of the complex, the one involving antagonistic interplay between the two hydrogen bonds (intermolecular C-H···O and intramolecular O-H···O) is preferred. The combined experimental and calculated data presented here suggest that in condensed media, conformational preferences are guided by directional hyperconjugative charge transfer interactions at the C-H···O hydrogen bonding site of the complex.

Water bridges anchored by a C–H⋯O hydrogen bond: the role of weak interactions in molecular solvation

Hydrogen-bonded water bridges are re-directed from a polar NH bond to a weakly activated C(2)–H bond upon N-methylation. Infrared spectra, supported by ab initio calculations, provide direct evidence of the role of the C(2)–H donor in the solvation of the imidazole ring.

Conformational preferences of monohydrated clusters of imidazole derivatives revisited

IR-UV double resonance spectroscopy was used to identify the conformers of monohydrated benzimidazole andN-methylbenzimidazole in a supersonic jet. A new OH–N bound conformer relevant to histidine containing proteins was discovered. The long standing differences in the literature about the relative energies and abundance of the monohydrated imidazole derivatives have also been resolved.

Communication: Where does the first water molecule go in imidazole?

Supersonic jet FTIR spectroscopy supplemented by (18)O substitution shows unambiguously that water prefers to act as an O-H···N hydrogen bond donor towards imidazole, instead of acting as a N-H···O acceptor. Previous matrix isolation, helium droplet, and aromatic substitution experiments had remained ambiguous, as are standard quantum chemical calculations. The finding is supported by a study of the analogous methanol complexes and by higher level quantum chemical calculations.

Evaluating how discrete water molecules affect protein-DNA π-π and π(+)-π stacking and T-shaped interactions: the case of histidine-adenine dimers

Changes in the magnitude of (M06-2X/6-31+G(d,p)) π-π stacking and T-shaped (nucleobase-edge and amino acid-edge) interactions between (neutral or protonated) histidine (His) and adenine (A) dimers upon microsolvation with up to four discrete water molecules were determined. A variety of histidine-water interactions were considered including conventional (N-H···O, N···H-O, C-H···O) hydrogen bonding and nonconventional (X-H···π (neutral His) or lone-pair···π (protonated His)) contacts. Overall, the effects of discrete His-H(2)O interactions on the neutral histidine-adenine π-π contacts are negligible (<3 kJ mol(-1) or 15%) regardless of the type of water binding, the number of water molecules bound, or the His-A dimer (stacked or (amino acid- or nucleobase-edge) T-shaped) configuration. This suggests that previously reported gas-phase binding strengths for a variety of neutral amino acid-nucleobase dimers are likely relevant for a wide variety of (microsolvated) environments. In contrast, the presence of water decreases the histidine-adenine π(+)-π interaction by up to 15 kJ mol(-1) (or 30%) for all water binding modes and orientations, as well as different stacked and T-shaped His(+)-A dimers. Regardless of the larger effect of discrete histidine-water interactions on the magnitude of the π(+)-π compared with π-π interactions, the π(+)-π binding strengths remain substantially larger than the corresponding π-π contacts. These findings emphasize the distinct nature of π(+)-π and π-π interactions and suggest that π(+)-π contacts can provide significant stabilization in biological systems relative to π-π contacts under many different environmental conditions.

Synthesis and Spectroscopic Investigation of Aggregation through Cooperative π−π and C−H···O Interactions in a Novel Pyrene Octaaldehyde Derivative

[reaction: see text] Synthesis of a pyrene octaaldehyde derivative and its aggregation through pi-pi and C-H...O interactions in solution and in the solid state probed by fluorescence emission and other spectroscopic methods are reported. The effect of solvent, concentration, and temperature on the aggregation is investigated.

One and two dimensional 1H and 13C high resolution NMR investigation of lariat ethers and their alkali metal ionic complexes: a more tangible evidence for the presence of less common C-H***O hydrogen bonds

The possible existence of less common hydrogen bonds in three lariat ethers and their alkali-metal ionic complexes have been investigated with one- and two-dimensional (1D and 2D) proton and carbon-13 high resolution liquid state NMR spectroscopy. The occurrence of hydrogen-bonding induced by the addition of metal ions has been identified with the observation of indirect dipolar coupling between the coupling partners involved in the hydrogen-bonding. The addition of metal ions, moreover, causes appreciable change of chemical shift of several protons and carbons. The chemical shift change depends on the ion radius, larger ions causing smaller change. Moreover, the change of chemical shift is in coincidence with the occurrence of hydrogen-bonding. The values of the coupling constants have been obtained for each of these hydrogen bonds and were used for evaluating the hydrogen-bond strength. An intriguing and surprising observation is that a C-H***O hydrogen bond identified in solution by this work was not found in the previous study with X-ray diffraction or other methods.

Kinase inhibitors and the case for CH...O hydrogen bonds in protein-ligand binding

Although the hydrogen bond is known to be an important mediator of intermolecular interactions, there has yet to be an analysis of the role of CH...O hydrogen bonds in protein-ligand complexes. In this work, we present evidence for such nonstandard hydrogen bonds from a survey of aromatic ligands in 184 kinase crystal structures and 358 high-resolution structures from the Protein Data Bank. CH groups adjacent to the positively charged nitrogen of nicotinamide exhibit geometric preferences strongly suggestive of hydrogen bonding interactions, as do heterocyclic CH groups in kinase ligands, while other aromatic CH groups do not exhibit these characteristics. Ab initio calculations reveal a considerable range of CH...O hydrogen bonding potentials among different aromatic ring systems, with nicotinamide and heterocycles preferred in kinase inhibitors showing particularly favorable interactions. These results provide compelling evidence for the existence of CH...O hydrogen bonds in protein-ligand interactions, as well as information on the relative strength of various aromatic CH donors. Such knowledge will be of considerable value in protein modeling, ligand design, and structure-activity analysis.

Comparison of various types of hydrogen bonds involving aromatic amino acids

Ab initio calculations are used to compare the abilities of the aromatic groups of the Phe, Tyr, Trp, and His amino acids (modeled respectively by benzene, phenol, indole, and imidazole) to form H-bonds of three different types. Strongest of all are the conventional H-bonds (e.g., OH..O and OH..N). His forms the strongest such H-bond, followed by Tyr, and then by Trp. Whereas OH..phi bonds formed by the approach of a proton donor to the pi electron cloud above the aromatic system are somewhat weaker, they nonetheless represent an important class of stabilizing interactions. The strengths of H-bonds in this category follow the trend Trp > His > Tyr approximately Phe. CH.O interactions are weaker still, and only those involving His and Trp are strong enough to make significant contributions to protein structure. A protonated residue such as HisH(+) makes for a very powerful proton donor, such that even its CH..O H-bonds are stronger than the conventional H-bonds formed by neutral groups.

CH...O hydrogen bonds at protein-protein interfaces

For the first time, a statistical potential has been developed to quantitatively describe the CH.O hydrogen bonding interaction at the protein-protein interface. The calculated energies of the CH.O pair interaction show a favorable valley at approximately 3.3 A, exhibiting a feature typical of an H-bond and similar to the ab initio quantum calculation result (Scheiner, S., Kar, T., and Gu, Y. (2001) J. Biol. Chem. 276, 9832-9837). The potentials have been applied to a set of 469 protein-protein complexes to calculate the contribution of different types of interactions to each protein complex: the average energy contribution of a conventional H-bond is approximately 30%; that of a CH.O H-bond is 17%; and that of a hydrophobic interaction is 50%. In some protein-protein complexes, the contribution of the CH.O H-bond can reach as high as approximately 40-50%, indicating the importance of the CH.O H-bond at the protein interface. At the interfaces of these complexes, C(alpha)H.O H-bonds frequently occur between adjacent strands in both parallel and antiparallel orientations, having the obvious structural motif of bifurcated H-bonds. Our study suggests that the weak CH.O H-bond makes an important contribution to the association and stability of protein complexes and needs more attention in protein-protein interaction studies.

Evidence for the preservation of specific intermolecular interactions in gaseous protein-oligosaccharide complexes

Arrhenius parameters, obtained with the blackbody infrared radiative dissociation technique, are reported for the dissociation of a series of gaseous protonated complexes composed of mutants of a single chain variable fragment (scFv) of the monoclonal antibody Se155-4 and three trisaccharide ligands. Hydrogen bonding between a ligand and a particular amino acid residue is identified from a comparison of activation energies measured for the complex of the unmodified scFv and the corresponding mutant. It is shown that the specific hydrogen bond between His(101H) and Man C-4 OH is preserved in the gaseous complex.

Thermal dissociation of protein-oligosaccharide complexes in the gas phase: mapping the intrinsic intermolecular interactions

Blackbody infrared radiative dissociation (BIRD) and functional group replacement are used to map the location and strength of hydrogen bonds between an antibody single chain fragment (scFv) and its natural trisaccharide receptor, alpha-D-Galp (1-->2)[alpha-D-Abep (1-->3)]alpha-D-Manp1-->OMe (1), in the gaseous, multiply protonated complex. Arrhenius activation parameters (E(a) and A) are reported for the loss of 1 and a series of monodeoxy trisaccharide congeners (5-8 identical with tri) from the (scFv + tri + 10H)(+10) complex. The energetic contribution of the specific oligosaccharide OH groups to the stability of the (scFv + 1 + 10H)(+10) complex is determined from the differences in E(a) measured for the trisaccharide analogues and 1 (55.2 kcal/mol). A decrease of 6 to 11 kcal/mol in E(a), measured for the monodeoxy trisaccharides, indicates that the deleted OH groups interact strongly with the scFv and that they account for a majority of the stabilizing intermolecular interactions. A partial map of the hydrogen bond donor/acceptor groups of 1 and the strength of the interactions is presented for the protonated +10 complex. A comparison of the gas-phase map with the crystal structure indicates that significant structural differences exist. The hydroxyl groups located outside of the binding pocket, and exposed to solvent in solution, participate in new protein-oligosaccharide hydrogen bonds in the gas phase. The decrease in kinetic and energetic stability of the (scFv + 2 + nH)(n)()(+) complex with increasing charge-state is attributed to conformational differences in the binding region induced by electrostatic repulsion. The similarity in the Arrhenius parameters for the +9 and +10 charge states suggests that repulsion effects on the structure of the binding region are negligible below +11.

Effects of high hydrostatic pressures on living cells: a consequence of the properties of macromolecules and macromolecule-associated water

Sixty percent of the Earth's biomass is found in the sea, at depths greater than 1000 m, i.e., at hydrostatic pressures higher than 100 atm. Still more surprising is the fact that living cells can reversibly withstand pressure shifts of 1000 atm. One explanation lies in the properties of cellular water. Water forms a very thin film around macromolecules, with a heterogeneous structure that is an image of the heterogeneity of the macromolecular surface. The density of water in contact with macromolecules reflects the physical properties of their different domains. Therefore, any macromolecular shape variations involving the reorganization of water and concomitant density changes are sensitive to pressure (Le Chatelier's principle). Most of the pressure-induced changes to macromolecules are reversible up to 2000 atm. Both the effects of pressure shifts on living cells and the characteristics of pressure-adapted species are opening new perspectives on fundamental problems such as regulation and adaptation.

Aspects of the mechanism of catalysis of glucose oxidase: a docking, molecular mechanics and quantum chemical study

The complex structure of glucose oxidase (GOX) with the substrate glucose was determined using a docking algorithm and subsequent molecular dynamics simulations. Semiempirical quantum chemical calculations were used to investigate the role of the enzyme and FAD co-enzyme in the catalytic oxidation of glucose. On the basis of a small active site model, substrate binding residues were determined and heats of formation were computed for the enzyme substrate complex and different potential products of the reductive half reaction. The influence of the protein environment on the active site model was estimated with a point charge model using a mixed QM/MM method. Solvent effects were estimated with a continuum model. Possible modes of action are presented in relation to experimental data and discussed with respect to related enzymes. The calculations indicate that the redox reaction of GOX differs from the corresponding reaction of free flavins as a consequence of the protein environment. One of the active site histidines is involved in substrate binding and stabilization of potential intermediates, whereas the second histidine is a proton acceptor. The former one, being conserved in a series of oxidoreductases, is also involved in the stabilization of a C4a-hydroperoxy dihydroflavin in the course of the oxidative half reaction.

A role for CH...O interactions in protein-DNA recognition

The concept of CH...O hydrogen bonds has recently gained much interest, with a number of reports indicating the significance of these non-classical hydrogen bonds in stabilizing nucleic acid and protein structures. Here, we analyze the CH...O interactions in the protein-DNA interface, based on 43 crystal structures of protein-DNA complexes. Surprisingly, we find that the number of close intermolecular CH...O contacts involving the thymine methyl group and position C5 of cytosine is comparable to the number of protein-DNA hydrogen bonds involving nitrogen and oxygen atoms as donors and acceptors. A comprehensive analysis of the geometries of these close contacts shows that they are similar to other CH...O interactions found in proteins and small molecules, as well as to classical NH...O hydrogen bonds. Thus, we suggest that C5 of cytosine and C5-Met of thymine form relatively weak CH...O hydrogen bonds with Asp, Asn, Glu, Gln, Ser, and Thr, contributing to the specificity of recognition. Including these interactions, in addition to the classical protein-DNA hydrogen bonds, enables the extraction of simple structural principles for amino acid-base recognition consistent with electrostatic considerations.