Ab Initio Quantum Mechanics Analysis of Imidazole C-H…O Water Hydrogen Bonding and a Molecular Mechanics Forcefield Correction

Abstract

While it is well established that classical hydrogen bonds play an important role in enzyme structure, function and dynamics, the role of weaker, but ‘activated’ C-H donor hydrogen bonds is poorly understood. The most important such case involves histidine which often plays a direct role in enzyme catalysis and possesses the most acidic C-H donor group of the standard amino acids. In the present study, we obtained optimized geometries and hydrogen bond interaction energies for C-H…O hydrogen bonded complexes between methane, ethylene, benzene, acetylene, and imidazole with water at the MP2-FC/6-31++G(2d,2p) and MP2-FC/aug-cc-pVDZ//MP2-FC/6-31++G(2d,2p) levels of theory. A strong linear relationship is obtained between the stability of the various hydrogen bonded complexes and both separation distances for H…0 and C—O. In general, these calculations indicate that C-H…0 interactions can be classified as hydrogen bonding interactions, albeit significantly weaker than the classical hydrogen bonds, but significantly stronger than just van der Waals interactions. For instance, while the electronic energy of stabilization at the MP2-FC/aug-cc-pVDZ//MP2-FC/6-31++G(2d,2p) level of theory of a water C-H…O water hydrogen bond is 4.36 kcal/mol more stable than the methane C-H…O water interaction, the water-water hydrogen bond is only 2.06 kcal/mol more stable than the imidazole C_e?H…O water hydrogen bond. Neglecting this latter hydrogen bonding interaction is obviously unacceptable. We next compare the potential energy surfaces for the imidazole C_e?H…O water and imidazole N_d?H…O hydrogen bonded complexes computed at the MP2/6-31++G(2d,2p) level of theory with the potential energy surface computed using the AMBER molecular mechanics program and forcefields. While the Weiner et al and Cornell et al AMBER forcefields reasonably account for the imidazole N-H…O water interaction, these forcefields do not adequately account for the imidazole C_e?H…O water hydrogen bond. A forcefield modification is offered that results in excellent agreement between the ab initio and molecular mechanics geometry and energy for this C-H…O hydrogen bonded complex.     

Effects of the position and manner of hydration on the stability of solvated N-methylacetamides and the strength of binding between N-methylacetamide and water clusters: a computational study

This work mainly studies the effects of the position (there are two possible hydrated sites) and the manner (i.e., whether water acts as a proton donor or acceptor) of hydration by various numbers of water molecules on the stability of 14 solvated N-methylacetamide structures, NMA-(H₂O) _n (n = 1–3), as well as the binding strength between the NMA and the water cluster, using molecular dynamics (MD) and B3LYP methods. Natural bond orbital (NBO) analysis is used to explore the origin of these effects. Some novel observations are obtained from the work. Our results show that monohydration at the carbonyl site favors stability and binding strength compared to monohydration at the amino site. Similarly, the preferred hydration at the carbonyl site is observed for dihydrated NMAs when the second water is added as a proton donor to the C=O group or the first water is H-bonded to the C=O group. However, unfavorable hydration at the C=O site occurs if the second water acts as a proton acceptor. Trihydration by a ring cluster of three water molecules at either the carbonyl site or the amino one yields relatively stable complexes, but significantly disfavors binding strength. The other trihydrated NMAs show similar behavior to dihydrated NMAs. In addition, our results show that the C=O and N–H frequencies can still be utilized to examine the H-bond effects of the water cluster.     

Theoretical investigations on the hydrogen bonding of nitrile isomers with H2O,HF, NH3 and H2S

The hydrogen bonds formed by the interaction of nitriles with water, hydrogen fluoride, ammonia and hydrogen sulphide have been studied using B3LYP and second-order Møller–Plesset perturbation (MP2) methods and 6-311+ + G(d,p) basis set. The energies and structures of 80 hydrogen-bonded complexes between nitriles and small molecules were examined systematically using B3LYP and MP2 procedure. Categorisation of the hydrogen bonds involved in the various complexes led to an ordering of hydrogen bond donor and acceptor abilities for some functional groups. The interaction energies have been corrected for the basis set superposition error using Boy's counterpoise correction method. The Morokuma energy decomposition analysis reveals that the strong interactions are due to the attractive contributions from the electrostatic (ES), polarisation (PL) and charge transfer (CT) components. The topological parameters, electron density and Laplacian of electron density show excellent correlation with the hydrogen bond length. Natural bond orbital (NBO) analysis has also been performed to study the CT from proton acceptor to the antibonding orbital of the H–Y bond in the proton donor part of complexes. The frequency analysis of C–H…Y bond in the complexes indicates the blue-shifting nature largely in case of sp² hybridised carbon atom.     

Refined structure of human carbonic anhydrase II at 2.0 Å resolution

The structure of human erythrocytic carbonic anhydrase II has been refined by constrained and restrained structure–factor least-squares refinement at 2.0 Å resolution. The conventional crystallographic R value is 17.3%. Of 167 solvent molecules associated with the protein, four are buried and stabilize secondary structure elements. The zinc ion is ligated to three histidyl residues and one water molecule in a nearly tetrahedral geometry. In addition to the zinc-bound water, seven more water molecules are identified in the active site. Assuming that Glu-106 is deprotonated at pH 8.5, some of the hydrogen bond donor–acceptor relations in the active site can be assigned and are described here in detail. The Oγ1 atom of Thr-199 donates its proton to the Oε1 atom of Glu-106 and can function as a hydrogen bond acceptor only in additional hydrogen bonds.     

Detection of N-H...N hydrogen bonding in RNA via scalar couplings in the absence of observable imino proton resonances

Hydrogen bond networks stabilize RNA secondary and tertiary structure and are thus essentially important for protein recognition. During structure refinements using either NMR or X-ray techniques, hydrogen bonds were usually inferred indirectly from the proximity of donor and acceptor functional groups. Recently, quantitative heteronuclear J(N,N)-HNN COSY NMR experiments were introduced that allowed the direct identification of donor and acceptor nitrogen atoms involved in hydrogen bonds. However, protons involved in base pairing interactions in nucleic acids are often not observable due to exchange processes. The application of a modified quantitative J(N,N)-HNN COSY pulse scheme permits observation of ^2hJ(N,N) couplings via non-exchangeable protons. This approach allowed the unambiguous identification of the A27·U23 reverse Hoogsteen base pair involved in a U-A·U base triple in the HIV-2 transactivation response element–argininamide complex. Despite a wealth of NOE information, direct evidence for this interaction was lacking due to the rapid exchange of the U23 imino proton. The ability to directly observe hydrogen bonds, even in D₂O and in the presence of rapid exchange, should facilitate structural studies of RNA.     

Study of conformation structure and molecular interactions of rosamycin by infrared spectrum analysis

IR spectra of rosamycin and its solutions in inert (CCl4 and C2Cl4), proton acceptor (tetrahydrofuran, hexametapol and diethylamine) and proton donor (CHCl3 and CH3OD) solvents were studied at various concentrations (0.1 to 0.001 mol/l) and temperatures (20 to 100 degrees C) in the region of the vC = O and vOH absorption bands (1600-1800 and 3200 3650 sm 1). It was found that the absorption bands at 3480 and 3560 sm-1 observed in the spectra of rosamycin diluted solutions in the inert solvents referred to variations of vOH...N of the aminosugar fragment and to vOH...O = C of the ester group of the macrocycle. Bands at 1697 and 1717 sm-1 referred to vC = O of the ketone and aldehyde carbonyl groups and band at 1728 sm-1 referred to vC = O of the ester group whose carbonyl was involved in the C = H...HO intramolecular hydrogen bond. Intensity of vC = O band (1745 sm-1) of the free ester group was nought. However, it increased with using the proton acceptor solvents. OH...N and OH...O = C intramolecular hydrogen bonds stabilized rosamycin molecule conformation. Mechanism of rosamycin interaction with the proton donor and acceptor molecules was elucidated. It was shown that tertiary nitrogen was the center of rosamycin molecule protonation.     

The influences of water solvent on the structures and stabilities of Na+–AD conformers

The influences of water solvent on the structures and stabilities of the complex ion conformers formed by the coordination of alanine dipeptide (AD) and Na⁺ have been investigated using supramolecular and polarizable continuum solvation models at the level of B3LYP/6-311++G**, respectively; 12 monohydrated and 12 dihydrated structures of Na⁺–AD complex ion were obtained after full geometrical optimization. The results showed that H₂O molecules easily bind with Na⁺ of Na⁺–AD complex ion, forming an ion-lone pair interaction with the Na–O bond length of 2.1–2.3 Å. Besides, H₂O molecules also can form hydrogen bonds O_W–H_W···O(1), O_W–H_W···O(2), N(1)–H(1)···O_W or N(2)–H(2)···O_W with O or N groups of the Na⁺–AD backbone. The most stable gaseous bidentate conformer C7AB of Na⁺–AD is still the most stable one in the solvent of water. However, the structure of the most unstable gaseous conformer α′B of Na⁺–AD collapses under the attack of H₂O molecules and changes into C7AB conformation. Computations with IEFPCM solvation model of self-consistent reaction field theory give that aqueous C5A is more stable than C7_eqB and that the stabilization energies of water solvent on monodentate conformers of Na⁺–AD complex ion (about 272–294 kJ/mol) are more than those on bidentate ones (about 243 kJ/mol).     

Theoretical ab Initio Study of the Effects of Methylation on Structure and Stability of G:C Watson-Crick Base Pair

Abstract

Methylation of DNA occurs most readily at N(3), N(7), and O(6) of purine bases and N(3) and O(2) of pyrimidines. Methylated bases are continuously formed through endogenous and exogenous mechanisms. The results of a theoretical ab initio study on the methylation of G:C base pair components are reported. The geometries of the local minima were optimized without symmetry restrictions by the gradient procedure at DFT level of theory and were verified by energy second derivative calculations. The standard 6–31G(d) basis set was used. The single-point calculations have been performed at the MP2/6–31G(d,p), MP2/6–31₊₊G(d,p), and MP2/6–311₊₊G(2d,2p) levels of theory. The geometrical parameters, relative stability and counterpoise corrected interaction energies are reported. Also, using a variation-perturbation energy decomposition scheme we have found the vital contributions to the total interaction energy.     

Electronic structure of alloxan and its dimers: QM/QD simulations and quantum chemical topology analysis

This study aims to identify the origin of the extra stability of alloxan, a biologically active pyrimidine. To achieve this goal, detailed DFT computations and quantum dynamics simulations have been performed to establish the most stable conformation and the global minimum structure on the alloxan potential energy surface. The effects of the solvent, basis set, and DFT method have been examined to validate the theoretical model adopted throughout the work. Two non-covalent intermolecular dimers of alloxan, the H-bonded and dipolar dimers, have been investigated at the ωB97X-D and M06-2X levels of theory using the triple zeta 6-311++G** to establish their relative stability. Quantum chemical topology features and natural bond orbital analysis (NBO) have been performed to identify and characterize the forces that govern the structures and underlie the extra stability of alloxan.     

On the possible influence of different pathways of binding of amides with water on “Pyramidalization” of valence bonds of peptide group nitrogen

With the aim to study solvation effects in peptide structure organization, the behavior of the energy of different types of hydration in simple amines and amides has been analyzed. On the example of quantum-chemical DFT and PM3 calculations of amino derivatives of composition CH₃-(CH₂)₃)-NH₂, (CH3)₂-NH, CH₃-NH₂, NH₃, CH₂=CH-NH₂, H-CC-NH₂, O=C(CH)₃-N(CH₃)₂, O=C(CH₃)-NH(CH₃), O=C(CH₃)-NH₂, O=CH-N(CH₃)H, and O=CH-NH₂ it has been shown that: (1) in the given set of molecules, the proton acceptor N…H-O variant of hydrogen bonding of NH₂ group with a water molecule is dominating only for the simplest amines. Being primordially weaker, the proton donor N-H…OH variant of water H-bonding gradually increases in energy in the given set as the basicity of the compound decreases, and for the case of amides of carboxylic acids it becomes already a significant channel of the hydration; (2) the intermolecular N-H…O=C bonding of trans-N-methylacetamides, which models the peptide hydrogen bonds in proteins, induces “planarization” of its initially nonplanar O=C-NH fragments. However, the addition of water molecules to the complex through the proton acceptor N…H-O variant of binding of N atom not only restores but even strengthens the “pyramidalization” of valence bonds of peptide groups.     

Many-body energies during proton transfer in an aqueous system

The energetics of the mechanism of proton transfer from a hydronium ion to one of the water molecules in its first solvation shell are studied using density functional theory and the Møller–Plesset perturbation (MP2) method. The potential energy surface of the proton transfer mechanism is obtained at the B3LYP and MP2 levels with the 6-311++G** basis set. Many-body analysis is applied to the proton transfer mechanism to obtain the change in relaxation energy, two-body, three-body and four-body energies when proton transfer occurs from the hydronium ion to one of the water molecules in its first solvation shell. It is observed that the binding energy (BE) of the complex decreases during the proton transfer process at both levels of theory. During the proton transfer process, the % contribution of the total two-body energy to the binding energy of the complex increases from 62.9 to 68.09% (39.9 to 45.95%), and that of the total three-body increases from 25.9 to 27.09% (24.16 to 26.17%) at the B3LYP/6-311++G** (MP2/ 6-311++G**) level. There is almost no change in the water–water–water three-body interaction energy during the proton transfer process at both levels of theory. The contribution of the relaxation energy and the total four-body energy to the binding energy of the complex is greater at the MP2 level than at the B3LYP level. Significant differences are found between the relaxation energies, the hydronium–water interaction energies and the four-body interaction energies at the B3LYP and MP2 levels.     

Computational study of the intramolecular proton transfer reactions of 3-hydroxytropolone (2,7-dihydroxycyclohepta-2,4,6-trien-1-one) and its dimers

The proton transfer reaction and dimerization processes of 3-hydroxytropolone (3-OHTRN) have been investigated using density functional theory (DFT) at the B3LYP/6–31+G** level. The influence of the solvent on the proton transfer reaction of 3-OHTRN was examined using the self-consistent isodensity polarized continuum model (SCI-PCM) with different dielectric constants (ε?=?4.9, CHCI₃; ε?=?32.63, CH₃OH; ε?=?78.39, H₂O). The intramolecular proton transfer reaction occurs more readily in the gas phase than in solution. Results also show that the stability of 3-OHTRN dimers in the gas phase is directly affected by the hydrogen bond length in the dimer structure.     

Dissecting non-canonical interactions in frameshift-stimulating mRNA pseudoknots

A variety of powerful NMR experiments have been introduced over the last few years that allow for the direct identification of different combinations of donor and acceptor atoms involved in hydrogen bonds in biomolecules. This ability to directly observe tertiary structural hydrogen bonds in solution tremendously facilitates structural studies of nucleic acids. We show here that an adiabatic HNN-COSY pulse scheme permits observation and measurement of J(N,N) couplings for nitrogen sites that are separated by up to 140 ppm in a single experiment at a proton resonance frequency of 500 MHz. Crucial hydrogen bond acceptor sites in nucleic acids, such as cytidine N3 nitrogens, can be unambiguously identified even in the absence of detectable H41 and H42 amino protons using a novel triple-resonance two-dimensional experiment, denoted H5(C5C4)N3. The unambiguous identification of amino nitrogen donor and aromatic nitrogen acceptor sites associated with both major groove as well as minor groove triple base pairs reveal the details of hydrogen bonding networks that stabilize the complex architecture of frameshift-stimulating mRNA pseudoknots. Another key tertiary interaction involving a 2′-OH hydroxyl proton that donates a hydrogen bond to an aromatic nitrogen acceptor in a cis Watson–Crick/sugar edge interaction can also be directly detected using a quantitative J(H,N) ¹H,¹⁵N-HSQC experiment.     

Influence of Cooperativity on Hydrogen Bond Networks

Abstract

Cooperative effects are known to strongly affect the geometrical, energetic and vibrational properties of hydrogen bonded systems. In particular, such effects strongly favor molecular arrangements where each molecule is simultaneously a donor and an acceptor of hydrogen bonds (HBs), regardless of the chemical nature of the monomer subunits. In the particular case of water systems, it has been shown that the more a molecule is a proton donor in HBs, the more the HBs where it is a proton acceptor are reinforced. Such a property could be at the origin of the equilibrium between the two species of hydrogen bonded water molecules in liquid water (one with a strong hydrogen bonding character, and one with a weaker one), as experimentally evidenced and as a molecular dynamic study of the small (H₂O)₂₄ cluster clearly suggests.     

A theoretical study on the hydrogen-bonding interactions between flavonoids and ethanol/water

Ethanol and water are the solvents most commonly used to extract flavonoids from propolis. Do hydrogen-bonding interactions exist between flavonoids and ethanol/water? In this work, this question was addressed by using density functional theory (DFT) to provide information on the hydrogen-bonding interactions between flavonoids and ethanol/water. Chrysin and Galangin were chosen as the representative flavonoids. The investigated complexes included chrysin–H₂O, chrysin–CH₃CH₂OH, galangin–H₂O and galangin–CH₃CH₂OH dyads. Molecular geometries, hydrogen-bond binding energies, charges of monomers and dyads, and topological analysis were studied at the B3LYP/M062X level of theory with the 6?31++G(d,p) basis set. The main conclusions were: (1) nine and ten optimized hydrogen-bond geometries were obtained for chrysin–H₂O/CH₃CH₂OH and galangin–H₂O/CH₃CH₂OH complexes, respectively. (2) The hydrogen atoms except aromatic H1 and H5 and all of the oxygen atoms can form hydrogen-bonds with H₂O and CH₃CH₂OH. Ethanol and water form strong hydrogen-bonds with the hydroxyl, carbonyl and ether groups in chrysin/galangin and form weak hydrogen-bonds with aromatic hydrogen atoms. Except in structures labeled A and B, chrysin and galangin interact more strongly with H₂O than CH₃CH₂OH. (3) When chrysin and galangin form hydrogen-bonds with H₂O and CH₃CH₂OH, charge transfers from the hydrogen-bond acceptor (H₂O and CH₃CH₂OH in structures A, B, G, H, I, J) to the hydrogen-bond donor (chrysin and galangin in structure A, B, G, H, I, J). The stronger hydrogen-bond makes the hydrogen-bond donor lose more charge (A> B> G> H> I> J). (4) Most of the hydrogen-bonds in chrysin/galangin?H₂O/CH₃CH₂OH complexes may be considered as electrostatic dominant, while C?O2···H in structures labeled E and C?O5···H in structures labeled J are hydrogen-bonds combined of electrostatic and covalent characters. H9, H7, and O4 are the preferred hydrogen-bonding sites.     

猴头菌丝多糖降血糖作用研究

目的:研究猴头菌丝多糖的降血糖作用。方法:以液体发酵生产的猴头菌丝体为原料,经热水浸提、浓缩、酒精沉淀获得菌丝粗多糖;以常规降糖药物格列本脲为阳性治疗对照,通过四氧嘧啶诱发小鼠糖尿病的预防试验,比较猴头菌丝多糖各剂量与格列本脲的降血糖效果。结果:猴头菌丝多糖得率为7.14%,粗多糖再经Sevage法去除蛋白质,获得猴头菌丝精多糖(HMP),得率为10.92%;猴头多糖高、中、低三个剂量均能有效的对抗四氧嘧啶诱发的高血糖;其中,高剂量的降血糖作用与格列本脲相比,差异极显著。结论:猴头菌丝多糖对四氧嘧啶型高血糖模型小鼠有降血糖作用,作用效果优于格列本脲,对糖尿病小鼠的胰腺具有一定的保护作用。     

Theoretical investigation on the structure and thermodynamic properties of the 2,4-dinitroimidazole complex with methanol

The structure and thermodynamic properties of the 2, 4-dinitroimidazole complex with methanol were investigated using the B3LYP and MP2(full) methods with the 6-31++G(2d,p) and 6-311++G(3df,2p) basis sets. Four types of hydrogen bonds [N–H?O, C–H?O, O–H?O (nitro oxygen) and O–H?π] were found. The hydrogen-bonded complex having the highest binding energy had a N–H?O hydrogen bond. Analyses of natural bond orbital (NBO) and atoms-in-molecules (AIM) revealed the nature of the intermolecular hydrogen-binding interaction. The changes in thermodynamic properties from monomers to complexes with temperatures ranging from 200.0 to 800.0 K were investigated using the statistical thermodynamic method. Hydrogen-bonded complexes of 2,4-dinitroimidazole with methanol are fostered by low temperatures.

A theoretical study on unusual intermolecular T-shaped X–H...π interactions between the singlet state HB=BH and HF, HCl, HCN or H2C2

The unusual T-shaped X–H...π hydrogen bonds are found between the B=B double bond of the singlet state HB=BH and the acid hydrogen of HF, HCl, HCN and H₂C₂ using MP2 and B3LYP methods at 6-311++G(2df,2p) and aug-cc-pVTZ levels. The binding energies follow the order of HB=BH...HF>HB=BH...HCl>HB=BH...HCN>HB=BH...H₂C₂. The hydrogen-bonded interactions in HB=BH...HX are found to be stronger than those in H₂C=CH₂...HX and OCB≡BCO...HX. The analyses of natural bond orbital (NBO) and the electron density shifts reveal that the nature of the T-shaped X–H...π hydrogen-bonded interaction is that much of the lost density from the π-orbital of B=B bond is shifted toward the hydrogen atom of the proton donor, leading to the electron density accumulation and the formation of the hydrogen bond. The atoms in molecules (AIM) theory have also been applied to characterize bond critical points and confirm that the B=B double bond can be a potential proton acceptor. The unusual T-shaped X–H...π hydrogen bonds are found between the B=B double bond of the singlet state HB=BH and the acid hydrogen of HF, HCl, HCN and H₂C₂     

Quantum-chemical analysis of C-H...O and C-H...N interactions in RNA base pairs--H-bond versus anti-H-bond pattern

Geometries and interaction energies of unusual UU and AA base pairs with one standard hydrogen bond (H-bond) and additional C-H...O or C-H...N contacts have been determined by quantum-chemical methods taking into account electron correlation. Whereas the C-H bond length in the UU C-H...O contact increases upon complex formation (H-bond pattern), the C-H bond of the AA C-H....N interaction is shortened (anti-H-bond pattern). The same properties are found for model complexes between U or A and formaldehyde that have intermolecular C-H...acceptor contacts but no standard H-bonds. Both the C-H...acceptor H-bond and anti-H-bond interactions are attractive. A possible influence of the donor CH group charge distribution on the interaction pattern is discussed.     

Successful prediction of the hydrogen bond network of the 3-oxo-12alpha-hydroxy-5beta-cholan-24-oic acid crystal from resolution of the crystal structure of deoxycholic acid and its three 3,12-dihydroxy epimers

Crystal structures of p-xylene-crystallized deoxycholic acid (3alpha,12alpha-dihydroxy-5beta-cholan-24-oic acid) and its three epimers (3beta,12alpha-; 3alpha,12beta-; and 3beta,12beta-) have been solved. Deoxycholic acid forms a crystalline (P21) complex with the solvent with a 2:1 stoichiometry whereas crystals of the three epimers do not form inclusion compounds. Crystals of the 3beta,12beta-epimer are hexagonal, whereas the 3alpha,12beta-and 3beta,12alpha-epimers crystallize in the P2(1)2(1)2(1) orthorhombic space group. The three hydrogen bond sites (two hydroxy groups, i. e. O3-H, and O12-H, and the carboxylic acid group of the side chain, O24bO24a-H) simultaneously act as hydrogen bond donors and acceptors. The hydrogen bond network in the crystals was analyzed and the following sequences have been observed: two chains (abcabc... or acbacb... ) and two rings (abc or acb), which constitute a complete set of all the possible sequences which can be drawn for an intermolecular hydrogen bond network formed by three hydrogen bond donor/acceptor sites forming crossing hydrogen bonds. The orientation of O3-H (alpha or beta) determines the sequence of the acceptor and the donor groups involved in the pattern: O24a --> O12 --> O3 --> O24b when it is alpha and O24a --> O3 --> O12--> O24B when it is beta. These observations were used to predict the hydrogen bond network of p-xylene-crystallized 3-oxo,12alpha-hydroxy-5beta-cholan-24-oic acid. This compound has two hydrogen bond donor and three potential hydrogen bond acceptor sites. According to the previous sequence set, this compound should crystallize in the monoclinic P21 system, should form a complex with the solvent, O24b should not participate in the hydrogen bond network, and the chain sequence O24a --> O12 --> O3 would be followed. All predictions were confirmed experimentally.