Direct evaluation of individual hydrogen bond energy in situ in intra- and intermolecular multiple hydrogen bonds system

The Role of Hydrogen Bonding in the Raman Spectral Signals of Caffeine in Aqueous Solution

The identification and quantification of caffeine is a common need in the food and pharmaceutical industries and lately also in the field of environmental science. For that purpose, Raman spectroscopy has been used as an analytical technique, but the interpretation of the spectra requires reliable and accurate computational protocols, especially as regards the Resonance Raman (RR) variant. Herein, caffeine solutions are sampled using Molecular Dynamics simulations. Upon quantification of the strength of the non-covalent intermolecular interactions such as hydrogen bonding between caffeine and water, UV-Vis, Raman, and RR spectra are computed. The results provide general insights into the hydrogen bonding role in mediating the Raman spectral signals of caffeine in aqueous solution. Also, by analyzing the dependence of RR enhancement on the absorption spectrum of caffeine, it is proposed that the sensitivity of the RR technique could be exploited at excitation wavelengths moderately far from 266 nm, yet achieving very low detection limits in the quantification caffeine content.

Water-Mediated Oxidation of Guanine by a Repair Enzyme: Simulation Using the ABEEM Polarizable Force Field

The recognition mechanism of oxidative damage in organisms has long been a research hotspot. Water is an important medium in the recognition process, but its specific role remains unknown. There is a need to develop a suitable force field that can adequately describe the electrostatic, hydrogen bond, and other interactions among the molecules in the complex system of the repair enzyme and oxidized base. The developing ABEEM polarizable force field (PFF) has been used to simulate the repaired enzyme hOGG1 and oxidized DNA (PDB ID: 1EBM) in a biological environment, and the corresponding results are better than those of the fixed-charge force fields OPLS/AA and AMBER OL15. 8-Oxo-G is recognized by Gln315 of hOGG1 mainly through hydrogen bonds mediated by continuous exchange of 2 water molecules. Phe319 and Cys253 are stacked on both sides of the π planes of bases to form sandwich structures. The charge polarization effect gives an important signal to drive the exchange of water molecules and maintains the recognition of oxidation bases by enzymes. The mediated main water molecule A and mediated auxiliary water molecule B together pull Gln315 to recognize 8-oxo-G by hydrogen bond interactions. Then, the charge polarization signal of solvent water molecule C with a large absolute charge causes the absolute charge of O atoms in water molecule A or B to increase by approximately 0.2 e, and water molecule A or B leaves Gln315 and 8-oxo-G. The other water molecule and water molecule C synergistically recognize 8-oxo-G with Gln315. Even though the water molecules between Gln315 and 8-oxo-G are removed, the MD simulation results show that water molecules appear between Gln315 and 8-oxo-G in a very short time (<2 ps). The dwell time of each water molecule is approximately 60 ps. The radial distribution function and dwell time support the correctness of the above mechanism. These polarization effects and hydrogen bonding interactions cannot be simulated by a fixed-charge force field.

Experimental and computational diagnosis of the fluxional nature of a benzene-1,3,5-tricarboxamide-based hydrogen-bonded dimer

Precise characterization of the hydrogen bond network present in discrete self-assemblies of benzene-1,3,5-tricarboxamide monomers derived from amino-esters (ester BTAs) is crucial for the construction of elaborated functional co-assemblies. For all ester BTA dimeric structures previously reported, ester carbonyls in the side chain acted as hydrogen bond acceptors, yielding well-defined dimers stabilized by six hydrogen bonds. The ester BTA monomer derived from glycine (BTA Gly) shows a markedly different self-assembly behaviour. We report herein a combined experimental and computational investigation aimed at determining the nature of the dimeric species formed by BTA Gly. Two distinct dimeric structures were characterized by single-crystal X-ray diffraction measurements. Likewise, a range of spectroscopic and scattering techniques as well as molecular modelling were employed to diagnose the nature of dynamic dimeric structures in toluene. Our results unambiguously establish that both ester and amide carbonyls are involved in the hydrogen bond network of the discrete dimeric species formed by BTA Gly. The participation of roughly 4.5 ester carbonyls and 1.5 amide carbonyls per dimer as determined by FT-IR spectroscopy implies that several conformations coexist in solution. Moreover, NMR analysis and modelling data reveal rapid interconversion between these different conformers leading to a symmetric structure on the NMR timescale. Rapid hydrogen bond shuffling between conformers having three (three), two (four), one (five) and zero (six) amide carbonyl groups (ester carbonyl groups, respectively) as hydrogen bond acceptors is proposed to explain the magnetic equivalence of the amide N-H on the NMR timescale. When compared to other ester BTA derivatives in which only ester carbonyls act as hydrogen bond acceptors, the fluxional behaviour of the hydrogen-bonded dimers of BTA Gly likely originates from a larger range of energetically favorable conformations accessible through rotation of the BTA side chains.

ABEEM/MM OH- Models for OH-(H2O)n Clusters and Aqueous OH-: Structures, Charge Distributions, and Binding Energies

Based on the atom-bond electronegativity equalization method fused into molecular mechanics (ABEEM/MM), two fluctuating charge models of OH^--water system were proposed. The difference between these two models is whether there is charge transfer between OH^- and its first-shell water molecules. The structures, charge distributions, charge transfer, and binding energies of the OH^-(H₂O)_n (n = 1-8, 10, 15, 23) clusters were studied by these two ABEEM/MM models, the OPLS/AA force field, the OPLS-SMOOTH/AA force field, and the QM methods. The results demonstrate that two ABEEM/MM models can search out all stable structures just as the QM methods, and the structures and charge distributions agree well with those from the QM calculations. The structures, the charge transfer, and the strength of hydrogen bonds in the first hydration shell are closely related to the coordination number of OH^-. Molecular dynamics simulations on the aqueous OH^- solution are performed at 298 and 278 K using ABEEM/MM-I model. The MD results show that the populations of three-, four-, and five-coordinated OH^- are 29.6%, 67.1%, and 3.4% at 298 K, respectively, and those of two-, three-, four-, and five-coordinated OH^- are 10.8%, 44.9%, 39.2%, and 4.9% at 278 K, respectively; the average hydrogen bond lengths and the hydrogen bond angle in the first shell increase with the temperature decreasing.

Reaction mechanism of NO with hydrolysates of NAMI-A: an MD simulation by combining the QM/MM(ABEEM) with the MD-FEP method

Nitrosylation reaction mechanisms of the hydrolysates of NAMI-A and hydrolysis reactions of ruthenium nitrosyl complexes were investigated in the triplet state and the singlet state. Activation free energies were calculated by combining the QM/MM(ABEEM) method with free energy perturbation theory, and the explicit solvent environment was simulated by an ABEEMσπ polarizable force field. Our results demonstrate that nitrosylation reactions of the hydrolysates of NAMI-A occur in both the triplet and the singlet states. The Ru-N-O angle of the triplet ruthenium nitrosyl complexes is in the range of 132.0°-138.2°. However, all the ruthenium nitrosyl complexes at the singlet state show an almost linear Ru-N-O angle. The nitrosylation reaction happens prior to the hydrolysis reaction for the first-step hydrolysates. The activation free energies of the nitrosylation reactions show that the H₂ O-NO exchange reaction of [RuCl₄ (Im)(H₂ O)] in the singlet spin sate is the most likely one. Comparing with the activation free energies of the hydrolysis reactions of the ruthenium nitrosyl complexes, the results indicate that the rate of the DMSO-H₂ O exchange reaction of [RuCl₃ (NO)(Im)(DMSO)] is faster than that of [RuCl₃ (H₂ O)(Im)(DMSO)] in both the triplet spin state and the singlet spin state. © 2018 Wiley Periodicals, Inc.

Combined covalent-electrostatic model of hydrogen bonding improves structure prediction with Rosetta

Interactions between polar atoms are challenging to model because at very short ranges they form hydrogen bonds (H-bonds) that are partially covalent in character and exhibit strong orientation preferences; at longer ranges the orientation preferences are lost, but significant electrostatic interactions between charged and partially charged atoms remain. To simultaneously model these two types of behavior, we refined an orientation dependent model of hydrogen bonds [Kortemme et al. J. Mol. Biol. 2003, 326, 1239] used by the molecular modeling program Rosetta and then combined it with a distance-dependent Coulomb model of electrostatics. The functional form of the H-bond potential is physically motivated and parameters are fit so that H-bond geometries that Rosetta generates closely resemble H-bond geometries in high-resolution crystal structures. The combined potentials improve performance in a variety of scientific benchmarks including decoy discrimination, side chain prediction, and native sequence recovery in protein design simulations and establishes a new standard energy function for Rosetta.