Investigations of an O-H...S hydrogen bond via Car-Parrinello and path integral molecular dynamics

Unraveling the Nature of Hydrogen Bonds of "Proton Sponges" Based on Car-Parrinello and Metadynamics Approaches

The nature of intra- and intermolecular non-covalent interactions was studied in four naphthalene derivatives commonly referred to as "proton sponges". Special attention was paid to an intramolecular hydrogen bond present in the protonated form of the compounds. The unsubstituted "proton sponge" served as a reference structure to study the substituent influence on the hydrogen bond (HB) properties. We selected three compounds substituted by methoxy, amino, and nitro groups. The presence of the substituents either retained the parent symmetry or rendered the compounds asymmetric. In order to reveal the non-covalent interaction properties, the Hirshfeld surface (HS) was computed for the crystal structures of the studied compounds. Next, quantum-chemical simulations were performed in vacuo and in the crystalline phase. Car-Parrinello molecular dynamics (CPMD), Path Integral Molecular Dynamics (PIMD), and metadynamics were employed to investigate the time-evolution changes of metric parameters and free energy profile in both phases. Additionally, for selected snapshots obtained from the CPMD trajectories, non-covalent interactions and electronic structure were studied. Quantum theory of atoms in molecules (QTAIM) and the Density Overlap Regions Indicator (DORI) were applied for this purpose. It was found based on Hirshfeld surfaces that, besides intramolecular hydrogen bonds, other non-covalent interactions are present and have a strong impact on the crystal structure organization. The CPMD results obtained in both phases showed frequent proton transfer phenomena. The proton was strongly delocalized in the applied time-scale and temperature, especially in the PIMD framework. The use of metadynamics allowed for tracing the free energy profiles and confirming that the hydrogen bonds present in "proton sponges" are Low-Barrier Hydrogen Bonds (LBHBs). The electronic and topological analysis quantitatively described the temperature dependence and time-evolution changes of the electronic structure. The covalency of the hydrogen bonds was estimated based on QTAIM analysis. It was found that strong hydrogen bonds show greater covalency, which is additionally determined by the proton position in the hydrogen bridge.

Sensitivity of Intra- and Intermolecular Interactions of Benzo[h]quinoline from Car-Parrinello Molecular Dynamics and Electronic Structure Inspection

The O-H...N and O-H...O hydrogen bonds were investigated in 10-hydroxybenzo[h]quinoline (HBQ) and benzo[h]quinoline-2-methylresorcinol complex in vacuo, solvent and crystalline phases. The chosen systems contain analogous donor and acceptor moieties but differently coupled (intra- versus intermolecularly). Car–Parrinello molecular dynamics (CPMD) was employed to shed light onto principle components of interactions responsible for the self-assembly. It was applied to study the dynamics of the hydrogen bonds and vibrational features as well as to provide initial geometries for incorporation of quantum effects and electronic structure studies. The vibrational features were revealed using Fourier transformation of the autocorrelation function of atomic velocity and by inclusion of nuclear quantum effects on the O-H stretching solving vibrational Schrödinger equation a posteriori. The potential of mean force (Pmf) was computed for the whole trajectory to derive the probability density distribution and for the O-H stretching mode from the proton vibrational eigenfunctions and eigenvalues incorporating statistical sampling and nuclear quantum effects. The electronic structure changes of the benzo[h]quinoline-2-methylresorcinol dimer and trimers were studied based on Constrained Density Functional Theory (CDFT) whereas the Electron Localization Function (ELF) method was applied for all systems. It was found that the bridged proton is localized on the donor side in both investigated systems in vacuo. The crystalline phase simulations indicated bridged proton-sharing and transfer events in HBQ. These effects are even more pronounced when nuclear quantization is taken into account, and the quantized Pmf allows the proton to sample the acceptor area more efficiently. The CDFT indicated the charge depletion at the bridged proton for the analyzed dimer and trimers in solvent. The ELF analysis showed the presence of the isolated proton (a signature of the strongest hydrogen bonds) only in some parts of the HBQ crystal simulation. The collected data underline the importance of the intramolecular coupling between the donor and acceptor moieties.

The use of variable temperature 13 C solid-state MAS NMR and GIPAW DFT calculations to explore the dynamics of diethylcarbamazine citrate

Experimental ¹³ C solid-state magic-angle spinning (MAS) Nuclear Magnetic Resonance (NMR) as well as Density-Functional Theory (DFT) gauge-including projector augmented wave (GIPAW) calculations were used to probe disorder and local mobility in diethylcarbamazine citrate, (DEC)⁺ (citrate)^- . This compound has been used as the first option drug for the treatment of filariasis, a disease endemic in tropical countries and caused by adult worms of Wuchereria bancrofti, which is transmitted by mosquitoes. We firstly present 2D ¹³ C─¹ H dipolar-coupling-mediated heteronuclear correlation spectra recorded at moderate spinning frequency, to explore the intermolecular interaction between DEC and citrate molecules. Secondly, we investigate the dynamic behavior of (DEC)⁺ (citrate)^- by varying the temperature and correlating the experimental MAS NMR results with DFT GIPAW calculations that consider two (DEC)⁺ conformers (in a 70:30 ratio) for crystal structures determined at 293 and 235 K. Solid-state NMR provides insights on slow exchange dynamics revealing conformational changes involving particularly the DEC ethyl groups.

N-oxide Derivatives: Car-Parrinello Molecular Dynamics and Electron Localization Function Study on Intramolecular Hydrogen Bonds

First-principle molecular dynamics simulations (CPMD) and electron localization function (ELF) topological analysis were performed for two N-oxides, 2-[( N, N-dimethylamino- N-oxy)methyl]-4,6-dimethylphenol (I) and 5,5'-dibromo-3-[(diethylamino)methyl]-2,2'-biphenol N-oxide (II). Special attention was paid to hydrogen bond dynamics and spectroscopic features. Simulations were carried out in vacuo and in the crystalline phase. It was found that for compound I proton transfer phenomena do not occur spontaneously. An opposite situation was noticed for one of the two intramolecular hydrogen bonds of compound II. The influence of the neighboring molecules (including microsolvation) and the crystal field on the hydrogen bond properties cannot be neglected. Hydrogen bond spectroscopic signatures were quite well reproduced by power spectra of atomic velocity. Metric and spectroscopic theoretical findings were compared with experimental data available, especially in a view of the fact that experimental X-ray diffraction and FT-IR data of compound II indicate strong delocalization of the N-oxide bridge proton. ELF analysis showed that for compound II the N-oxide bridge belongs to the topological category of strong hydrogen bonds with the separate ELF basin for the bridged hydrogen atom.

Resonance-assisted stabilisation of hydrogen bonds probed by NMR spectroscopy and path integral molecular dynamics

Path integral molecular dynamics and experimental NMR data are used to investigate resonance-assisted hydrogen bonds (RAHBs). When nuclear delocalisation is included in chemical shift calculations, the agreement with experiment is excellent, while static calculations show very poor performance. The results support the concept of RAHB, which has recently been questioned.

N-H⋯O versus O-H⋯O: density functional calculation and first principle molecular dynamics study on a quinoline-2-carboxamide N-oxide

N-oxide-type compounds are the object of current research interest due to the presence of resonance-assisted N–H⋯O hydrogen bonds. Here, the metric and spectroscopic parameters of N-methyl-quinoline-2-carboxamide 1-oxide were computed on the basis of density functional theory and Car-Parrinello molecular dynamics. Computations were performed in vacuo and in solid state; for both phases additional simulations with Grimme’s dispersion correction were carried out. The approaches used were able to reproduce correctly the structural aspects of the studied compound and shed more light on the hydrogen bonding with special focus on bridge proton mobility. Proton transfer phenomena were found not to occur in the investigated compound, and the bridge proton was localized to the donor site. This observation is in agreement with the classical theory of the acidity of donor–acceptor sites. The presence of hydrogen bonding was confirmed using atoms-in-molecules theory. The computational results were compared with available experimental data.

Flexibility of a biotinylated ligand in artificial metalloenzymes based on streptavidin--an insight from molecular dynamics simulations with classical and ab initio force fields

In the field of enzymatic catalysis, creating activity from a non catalytic scaffold is a daunting task. Introduction of a catalytically active moiety within a protein scaffold offers an attractive means for the creation of artificial metalloenzymes. With this goal in mind, introduction of a biotinylated d(6)-piano-stool complex within streptavidin (SAV) affords enantioselective artificial transfer-hydrogenases for the reduction of prochiral ketones. Based on an X-ray crystal structure of a highly selective hybrid catalyst, displaying significant disorder around the biotinylated catalyst [eta(6)-(p-cymene)Ru(Biot-p-L)Cl], we report on molecular dynamics simulations to shed light on the protein-cofactor interactions and contacts. The results of these simulations with classical force field indicate that the SAV-biotin and SAV-catalyst complexes are more stable than ligand-free SAV. The point mutations introduced did not affect significantly the overall behavior of SAV and, unexpectedly, the P64G substitution did not provide additional flexibility to the protein scaffold. The metal-cofactor proved to be conformationally flexible, and the S112K or P64G mutants proved to enhance this effect in the most pronounced way. The network of intermolecular hydrogen bonds is efficient at stabilizing the position of biotin, but much less at fixing the conformation of an extended biotinylated ligand. This leads to a relative conformational freedom of the metal-cofactor, and a poorly localized catalytic metal moiety. MD calculations with ab initio potential function suggest that the hydrogen bonds alone are not sufficient factors for full stabilization of the biotin. The hydrophobic biotin-binding pocket (and generally protein scaffold) maintains the hydrogen bonds between biotin and protein.

Proton-transfer dynamics in the (HCO3 −)2 dimer of KHCO3 from Car–Parrinello and path-integrals molecular dynamics calculations

The proton motion in the (HCO_3^-)2 dimer of KHCO3 at 298 K has been studied with Car–Parrinello molecular dynamics (CPMD) and path-integrals molecular dynamics (PIMD) simulations. According to earlier neutron diffraction studies at 298 K hydrogen is disordered and occupies two positions with an occupancy ratio of 0.804/0.196. A simulation with only one unit cell is not sufficient to reproduce the disorder of the protons found in the experiments. The CPMD results with four cells, 0.783/0.217, are in close agreement with experiment. The motion of the two protons along the O...O bridge is highly correlated inside one dimer, but strongly uncoupled between different dimers. The present results support a mechanism for the disorder which involves proton transfer from donor to acceptor and not orientational disordering of the entire dimer. The question of simultaneous or successive proton transfer in the two hydrogen bonds in the dimer remains unanswered. During the simulation situations with almost simultaneous proton transfer with a time gap of around 1 fs were observed, as well as successive processes where first one proton is transferred and then the second one with a time gap of around 20 fs. The calculated vibrational spectrum is in good agreement with the experimental IR spectrum, but a slightly different assignment of the bands is indicated by the present simulations.