Directional hydrogen bonding in the MM3 force field: II

Extensive calculations on hydrogen bonded systems were carried out using the improved MM3 directional hydrogen bond potential. The resulting total function was reoptimized. Comparisons of the hydrogen bonding potential function from ab initio calculations (MP2/6-31G**); the original MM3(89); and the reoptimized MM3 force field MM3(96), for a variety of C, N, O, and Cl systems including the formamide dimer and formamide–water complex, are described herein. Hydrogen bonding is shown to be a far more complicated and ubiquitous phenomenon than is generally recognized. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 1001–1016, 1998     

Alcohols,ethers, carbohydrates,and related compounds. III. The 1,2-dimethoxyethane system

Ethylene glycol, its dimethyl ether, and some related compounds have been studied using the MM4 molecular mechanics force field. The MM4 calculated structural and energetic results have been brought into satisfactory agreement with a considerable number of experimental data and MP2/6-311++G(2d,2p) ab initio calculations. The heats of formation of these compounds are also well calculated. The MM4 ethylene glycol conformations in particular are in good agreement, both geometrically and in terms of energy, with those from the ab initio calculations. The corresponding dimethyl ether is of special interest, because it has been suggested that the trans-gauche conformation is unusually stable due to the hydrogen bonding of a hydrogen on a methyl group with the more distant oxygen. It is shown in the present work that while this conformation is more stable than might have been expected, the energy is adequately calculated by MM4 without using any hydrogen bonding between the Cbond;H bond and the oxygen. If such hydrogen bonding occurs, it amounts to no more than about 0.5 kcal/mol in energy, and is too small to detect with certainty. Additionally, energetic relationships in trans-1,2-dimethoxycyclohexane, 1,3,5,7-tetraoxadecalin, and 3-methoxytetrahydropyran have been studied, and the calculated results are compared with experimental information, which is adequately reproduced.     

Amide–Amide Hydrogen Bonding. Semirubin Amides

Summary. Semirubins are analogs for one-half of the bilirubin structure and capable of intramolecular hydrogen bonding. Semirubin amides of ammonia and primary amines are also capable of intramolecular hydrogen bonding. From a combination of spectroscopic methods (¹H NMR, NOE, and VPO), the primary amide is found to engage very effectively in intramolecular hydrogen bonding. The secondary and tertiary amides engage in both intramolecular (i) and intermolecular (ii) hydrogen bonding: N-methyl (i, monomer + ii, dimer), N-tert-butyl (ii, dimer), N,N-diethyl (i, monomer + ii, dimer). With an oxo-group at C(10), all of the amides are monomeric and most engage in intramolecular hydrogen bonding.     

Amide–Amide Hydrogen Bonding. Semirubin Amides

Semirubins are analogs for one-half of the bilirubin structure and capable of intramolecular hydrogen bonding. Semirubin amides of ammonia and primary amines are also capable of intramolecular hydrogen bonding. From a combination of spectroscopic methods (¹H NMR, NOE, and VPO), the primary amide is found to engage very effectively in intramolecular hydrogen bonding. The secondary and tertiary amides engage in both intramolecular (i) and intermolecular (ii) hydrogen bonding: N-methyl (i, monomer + ii, dimer), N-tert-butyl (ii, dimer), N,N-diethyl (i, monomer + ii, dimer). With an oxo-group at C(10), all of the amides are monomeric and most engage in intramolecular hydrogen bonding.     

Young’s modulus calculations for cellulose I<Subscript>β</Subscript> by MM3 and quantum mechanics

Quantum mechanics (QM) and molecular mechanics (MM) calculations were performed to elucidate Young’s moduli for a series of cellulose I_β models. Computations using the second generation empirical force field MM3 with a disaccharide cellulose model, 1,4′-O-dimethyl-β-cellobioside (DMCB), and an analogue, 2,3,6,2′,3′,6′-hexadeoxy-1,4′-O-dimethyl-β-cellobioside (DODMCB), that cannot make hydrogen bonds reveal a considerable contribution of intramolecular hydrogen bonding to the molecular stiffness of cellulose I_β; the moduli for DMCB and DODMCB being 85.2 and 37.6 GPa, respectively. QM calculations confirm this contribution with modulus values of 99.7 GPa for DMCB and 33.0 GPa for DODMCB. However, modulus values for DMCB were considerably lower than values previously reported for cellulose I_β. MM calculations with extended cellulose chains (10–40 glucose units) resulted in modulus values, 126.0–147.5 GPa, more akin to the values reported for cellulose I_β. Comparison of the cellodecaose model, 1,4′-O-dimethyl-β-cellodecaoside (DMCD), modulus with that of its hydrogen bonding-deficient analogue, 2,3,6,2′,3′,6′-hexadeoxy-1,4′-O-dimethyl-β-cellodecaoside (DODMCD), corroborates the observed stiffness conferred by intramolecular hydrogen bonds; the moduli for DMCD and DODMCD being 126.0 and 63.3 GPa, respectively. Additional MM3 determinations revealed that modulus values were not strongly affected by intermolecular hydrogen bonding, with multiple strand models providing values similar to the single strand models; 87.5 GPa for a 7-strand DMCB model and 129.5 GPa for a 7 strand DMCD model.     

Correlation effects on the structure and dynamics of the H3O+ hydrate: B3LYP/MM and MP2/MM MD simulations

Two combined quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) simulations, namely B3LYP/MM and MP2/MM, have been performed to investigate the possible influence of electron correlation on the structure and dynamics of the H(3)O(+) hydrate. In comparison to the previously published HF/MM results, both B3LYP/MM and MP2/MM simulations clearly reveal stronger H(3)O(+)-water hydrogen bond interactions, which are reflected in a slightly greater compactness of the H(3)O(+) hydrate. However, the B3LYP/MM simulation, although providing structural details very close to the MP2/MM data, shows an artificially slow dynamic nature of some first shell water molecules as a consequence of the formation of a long-lived H(3)O(+)···H(2)O hydrogen bonding structure.     

Nature of the Fe-O2 bonding in oxy-myoglobin: effect of the protein

The nature of the Fe-O2 bonding in oxy-myoglobin was probed by theoretical calculations: (a) QM/MM (hybrid quantum mechanical/molecular mechanical) calculations using DFT/MM and CASSCF/MM methods and (b) gas-phase calculations using DFT (density functional theory) and CASSCF (complete active space self-consistent field) methods. Within the protein, the O2 is hydrogen bonded by His64 and the complex feels the bulk polarity of the protein. Removal of the protein causes major changes in the complex. Thus, while CASSCF/MM and DFT/MM are similar in terms of state constitution, degree of O2 charge, and nature of the lowest triplet state, the gas-phase CASSCF(g) species is very different. Valence bond (VB) analysis of the CASSCF/MM wave function unequivocally supports the Weiss bonding mechanism. This bonding arises by electron transfer from heme-Fe(II) to O2 and the so formed species coupled then to a singlet state Fe(III)-O2(-) that possesses a dative sigma(Fe-O) bond and a weakly coupled pi(Fe-O2) bond pair. The bonding mechanism in the gas phase is similar, but now the sigma(Fe-O) bond involves higher back-donation from O2(-) to Fe(III), while the constituents of pi(Fe-O2) bond pair have greater delocalization tails. The protein thus strengthens the Fe(III)-O2(-) character of the complex and thereby affects its bonding features and the oxygen binding affinity of Mb. The VB model is generalized, showing how the protein or the axial ligand of the oxyheme complex can determine the nature of its bonding in terms of the blend of the three bonding models: Weiss, Pauling, and McClure-Goddard.     

Magnetic nonequivalence in NMR spectroscopy—II. The role of hydrogen bonding in 2,2,4-trimethylpentane-1,3-diol and its monoisobutyrates

The methylene protons in 2,2,4-trimethylpentane-1,3-diol ( 1 ), its 1-isobutyrate ( 2 ) and 3-isobutyrate ( 3 ) monoesters, are magnetically nonequivalent with respect to chemical shift. Under the same experimental conditions, the chemical shift differences are 14·4 Hz, 46·2 Hz and 18·2 Hz for 1 , 2 and 3 , respectively. These nonequivalences were examined as a function of temperature and concentration. The large nonequivalence in 2 is attributed to a highly restricted conformation resulting from intramolecular hydrogen bonding. Less favorable steric interactions result in a substantially smaller contribution of intramolecular hydrogen bonding in 3 , while in 1 intermolecular hydrogen bonding occurs exclusively. Infrared data were obtained in support of these conclusions.     

Ab initio studies of molecules with N---C---C=O units. Part 2. 1-Amino-2-propanone, 2-methylaminoethanal and 2-aminopropanal

A complete conformational analysis of all the monomethylated derivatives of 2-aminoethanal (2AE) was carried out using MO ab initio with the 6-31G^** basis set, and it was inferred that methylation produces an increase in stability with respect to the initial compound 2AE. Geometric tendencies related to the existence of intramolecular hydrogen bonding and anomeric effects are discussed. Finally, the differences between the ab initio results and those produced by the current consistent MM3(92) force field are analysed.     

Theoretical study on the hydration of hydrogen peroxide in terms of ab initio method and atom-bond electronegativity equalization method fused into molecular mechanics

In this paper, the interaction between hydrogen peroxide (HP) and water were systemically studied by atom-bond electronegativity equalization method fused into molecular mechanics (ABEEM/MM) and ab initio method. The results show that the optimized geometries, interaction energies and dipole moments of hydrated HP clusters HP(H₂O) _n (n = 1–6) calculated by ABEEM/MM model are fairly consistent with the MP2/aug-cc-pVTZ//MP2/aug-cc-pVDZ results. The ABEEM/MM results indicate that n = 4 is the transition state structure from 2D planar structure to 3D network structure. The variations of the average hydrogen bond length with the increasing number of water molecules given by ABEEM/MM model agree well with those of ab initio studies. Moreover, the radial distribution functions (RDFs) of water molecule around HP in HP aqueous solution have been analyzed in detail. It can be confirmed that HP is a good proton donor and poor proton acceptor in aqueous solution by analysis of the RDFs.     

Simulations of liquid ammonia based on the combined quantum mechanical/molecular mechanical (QM/MM) approach

Two combined quantum mechanics/molecular mechanics (QM/MM) molecular dynamics simulations, namely, HF/MM and B3LYP/MM, have been performed to investigate the local structure and dynamics of liquid ammonia. The most interesting region, a sphere containing a central reference molecule and all its nearest surrounding molecules (first coordination shell), was treated by the Hartree-Fock (HF) and hybrid density functional B3LYP methods, whereas the rest of the system was described by the classical pair potentials. On the basis of both HF and B3LYP methods, it is observed that the hydrogen bonding in this peculiar liquid is weak. The structure and dynamics of this liquid are suggested to be determined by the steric packing effects, rather than by the directional hydrogen bonding interactions. Compared to previous empirical as well as Car-Parrinello (CP) molecular dynamics studies, our QM/MM simulations provide detailed information that is in better agreement with experimental data.     

Hydrogen‐Bonding Toughened Hydrogels and Emerging CO2‐Responsive Shape Memory Effect

A double hydrogen bonding (DHB) hydrogel is constructed by copolymerization of 2‐vinyl‐4,6‐diamino‐1,3,5‐triazine (hydrophobic hydrogen bonding monomer) and N,N‐dimethylacrylamide (hydrophilic hydrogen bonding monomer) with polyethylene glycol diacrylates. The DHB hydrogels demonstrate tunable robust mechanical properties by varying the ratio of hydrogen bonding monomer or crosslinker. Importantly, because of synergistic energy dissipating mechanism of strong diaminotriazine (DAT) hydrogen bonding and weak amide hydrogen bonding, the DHB hydrogels exhibit high toughness (up to 2.32 kJ m⁻²), meanwhile maintaining 0.7 MPa tensile strength, 130% elongation at break, and 8.3 MPa compressive strength. Moreover, rehydration can help to recover the mechanical properties of the cyclic loaded–unloaded gels. Attractively, the DHB hydrogels are responsive to CO₂ in water, and demonstrate unprecedented CO₂‐triggered shape memory behavior owing to the reversible destruction and reconstruction of DAT hydrogen bonding upon passing and degassing CO₂ without introducing external acid. The CO₂ triggering mechanism may point out a new approach to fabricate shape memory hydrogels.     

An ONIOM study of the QA site semiquinone in the Rhodobacter sphaeroides photosynthetic reaction centre

ONIOM (QM/MM) calculations are performed to investigate the spin density distribution for the ubisemiquinone anion radical in the Q_A binding site of the photosynthetic bacterium Rhodobacter sphaeroides. The calculated spin density in the Q_A site model suggests that differential hydrogen bonding strength to the O1 and O4 oxygen atoms of the radical results in an asymmetric spin density distribution in the semiquinone anion free radical form. The origin of the spin density asymmetry is attributed to the presence of the divalent iron or zinc ion situated between the Q_A and Q_B sites.     

Importance of selecting proper basis set in quantum mechanical studies of potential energy surfaces of carbohydrates

An extensive quantum mechanical study of a water dimer suggests that the introduction of a diffuse function into the basis set, which significantly reduces the basis set superposition error (BSSE) in the hydrogen bonding energy calculation, is the key to better calculations of the potential energy surfaces of carbohydrates. This article examines the potential energy surfaces of selected d -aldo- and d -ketohexoses (a total of 82 conformers) by quantum mechanics (QM) and molecular mechanics (MM) methods. In contrast to the results with a smaller basis set (B3LYP/6-31G** 5d), we found at the higher level calculation (B3LYP/6-311++G(2d,2p)//B3LYP/6-31G** 5d) that, in most cases, the furanose forms are less stable than the pyranose forms. These discrepancies are mainly due to the fact that intramolecular hydrogen bonding energies are overestimated in the lower level calculations. The higher level QM calculations of the potential energy surfaces of d -aldo- and d -ketohexoses now are more comparable to the MM3 results. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 1593–1603, 1999     

Investigation of Steric Influences on Hydrogen‐Bonding Motifs in Cyclic Ureas by Using X‐Ray,Neutron, and Computational Methods

A series of urea‐derived heterocycles, 5N‐substituted hexahydro‐1,3,5‐triazin‐2‐ones, has been prepared and their structures have been determined for the first time. This family of compounds only differ in their substituent at the 5‐position (which is derived from the corresponding primary amine), that is, methyl ( 1 ), ethyl ( 2 ), isopropyl ( 3 ), tert‐butyl ( 4 ), benzyl ( 5 ), N,N‐(diethyl)ethylamine ( 6 ), and 2‐hydroxyethyl ( 7 ). The common heterocyclic core of these molecules is a cyclic urea, which has the potential to form a hydrogen‐bonding tape motif that consists of self‐associative (8) dimers. The results from X‐ray crystallography and, where possible, Laue neutron crystallography show that the hydrogen‐bonding motifs that are observed and the planarity of the hydrogen bonds appear to depend on the steric hindrance at the α‐carbon atom of the N substituent. With the less‐hindered substituents, methyl and ethyl, the anticipated tape motif is observed. When additional methyl groups are added onto the α‐carbon atom, as in the isopropyl and tert‐butyl derivatives, a different 2D hydrogen‐bonding motif is observed. Despite the bulkiness of the substituents, the benzyl and N,N‐(diethyl)ethylamine derivatives have methylene units at the α‐carbon atom and, therefore, display the tape motif. The introduction of a competing hydrogen‐bond donor/acceptor in the 2‐hydroxyethyl derivative disrupts the tape motif, with a hydroxy group interrupting the N? H???O?C interactions. The geometry around the hydrogen‐bearing nitrogen atoms, whether planar or non‐planar, has been confirmed for compounds 2 and 5 by using Laue neutron diffraction and rationalized by using computational methods, thus demonstrating that distortion of O‐C‐N‐H torsion angles occurs to maintain almost‐linear hydrogen‐bonding interactions.     

Detection of halogen bond formation by correlation of proton solvent shifts. 1. Haloforms in n-electron donor solvents

The solvent shifts of haloformic protons, (Cl₃CH, Br₃CH, I₃CH), have been measured in 24 n-electron donor solvents consisting of halogenated hydrocarbons, esters, ketones, ethers and amines. Deviations of Δ_Br and Δ₁ from linear dependence with Δ_Cl are indicative of the presence of halogen bond formation competitive with hydrogen bonding interactions. Bromoform interacts predominantly by hydrogen bonding, halogen bonding being detected to a small extent in chlorinated hydrocarbons and amines. Iodoform shows halogen bonding interactions which increase in relative importance to hydrogen bonding with solvent basicity. Halogen bonding is predominant for solutions of iodoform in amines.     

Alcohols,ethers, carbohydrates,and related compounds. IV. Carbohydrates

Ab initio calculations [B3LYP/6-311++G(2d,2p)] have been carried out on 84 conformations of 12 different sugars (hexoses), in both pyranose and furanose forms, with the idea of generating a data base for carbohydrate structural energies that may be used for developing the predictive value of molecular mechanics calculations for carbohydrates. The average value for the apparent gas phase anomeric effect for a series of 31 pairs of pyranose conformations was found to be 1.83 kcal/mol (vs. 2.67 kcal/mol with a smaller basis set used in earlier calculations). In developing MM4 to reproduce these data, it was necessary first to have good energies for simple alcohols and ethers, together with an adequate treatment of hydrogen bonding, and then to include the anomeric effect, and the ethylene glycol type system, as was previously recognized. It was also found that the so-called delta-2 effect, long recognized in carbohydrates, must be explicitly included, in order to obtain acceptable results. When a force field that included all of these items as developed from the small molecules based on the MM4 hydrocarbon force field was applied without any parameter adjustment to the set of hexopyranose and furanose conformations mentioned earlier, the E(beta) - E(alpha) was found to have an average value of 1.88 kcal/mol, versus 1.74 for the quantum calculations. The signed average and RMS deviations of the MM4 from the QM results were +0.15 and 0.87 kcal/mol.     

Conformational analysis of 2-substituted alkylphosphoryl compounds. Part 3: (2-hydroxypentyl)diphenylphosphine oxide and its acetate

Conformational analysis of (2-hydroxypentyl)diphenylphosphine oxide 1 and its acetate 2 is described. The NMR, X-ray, IR, and molecular mechanics (MM) modeling studies indicated that phosphine oxide 1 favors different conformers in the solid state and in solution and that conformational preferences are strongly influenced by the nature of the hydrogen bonding. Vicinal proton-proton coupling constants and MM modeling solvation studies indicated that there are deviations from perfectly staggered conformers. Conformational analysis based on the twisted staggered conformers for phosphine oxide 1 made marked changes to the estimated conformational populations. For the acetate 2 , NMR spectroscopy established that the position of the conformational equilibrium is strongly dependent on the polarity of the medium. © 1996 John Wiley & Sons, Inc.     

Whether proton transition to the triphosphate tail of ATP occurs at protein kinase environment: A Car‐Parrinello ab initio molecular dynamics study

Protonation state of the triphosphate tail of ATP (adenosine triphosphate) in protein environment is a fundamental issue, which has significant impact on the mechanism investigation of biochemical processes with ATP involved. Proton transition from surroundings (water molecule coordinating to magnesium, H_W; amino group of Lys, H_L) to the ATP tail in the catalytic core of protein kinase found recently disproved the commonly accepted deprotonation state of ATP tail. In this account, Car‐Parrinello ab initio molecular dynamics (CP‐AIMD) method has been employed to examine whether the proton transition occurs. To provide a comparison basis for the dynamics simulations, static quantum mechanics (QM), and combined quantum mechanics and molecular mechanics (QM/MM) calculations have also been carried out. Consistent results have been obtained that complete transition of hydrogen from the surroundings to the triphosphate tail of ATP is not allowed. The most dominant conformations correspond to the ones with H_W bonding to O(W) and H‐bonding to O(ATP), [O(W)‐H_W···O(ATP)], H_L bonding to N(Lys) and H‐bonding to O(ATP), [N(Lys)‐H_L···O(ATP)]. Metastable structures with one hydrogen atom bonding with two heavy atoms (hydrogen acceptors) were also located by our dynamic simulations. This bonding mode can satisfy the hungering for hydrogen of the two heavy atoms simultaneously. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2008