Revealing Intra- and Intermolecular Interactions Determining Physico-Chemical Features of Selected Quinolone Carboxylic Acid Derivatives

The intra- and intermolecular interactions of selected quinolone carboxylic acid derivatives were studied in monomers, dimers and crystals. The investigated compounds are well-recognized as medicines or as bases for further studies in drug design. We employed density functional theory (DFT) in its classical formulation to develop gas-phase and solvent reaction field (PCM) models describing geometric, energetic and electronic structure parameters for monomers and dimers. The electronic structure was investigated based on the atoms in molecules (AIM) and natural bond orbital (NBO) theories. Special attention was devoted to the intramolecular hydrogen bonds (HB) present in the investigated compounds. The characterization of energy components was performed using symmetry-adapted perturbation theory (SAPT). Finally, the time-evolution methods of Car–Parrinello molecular dynamics (CPMD) and path integral molecular dynamics (PIMD) were employed to describe the hydrogen bond dynamics as well as the spectroscopic signatures. The vibrational features of the O-H stretching were studied using Fourier transformation of the autocorrelation function of atomic velocity. The inclusion of quantum nuclear effects provided an accurate depiction of the bridged proton delocalization. The CPMD and PIMD simulations were carried out in the gas and crystalline phases. It was found that the polar environment enhances the strength of the intramolecular hydrogen bonds. The SAPT analysis revealed that the dispersive forces are decisive factors in the intermolecular interactions. In the electronic ground state, the proton-transfer phenomena are not favourable. The CPMD results showed generally that the bridged proton is localized at the donor side, with possible proton-sharing events in the solid-phase simulation of stronger hydrogen bridges. However, the PIMD enabled the quantitative estimation of the quantum effects inclusion—the proton position was moved towards the bridge midpoint, but no qualitative changes were detected. It was found that the interatomic distance between the donor and acceptor atoms was shortened and that the bridged proton was strongly delocalized.     

Intermolecular Interactions and Spectroscopic Signatures of the Hydrogen-Bonded System—n-Octanol in Experimental and Theoretical Studies

n-Octanol is the object of experimental and theoretical study of spectroscopic signatures and intermolecular interactions. The FTIR measurements were carried out at 293 K for n-octanol and its deuterated form. Special attention was paid to the vibrational features associated with the O-H stretching and the isotope effect. Density Functional Theory (DFT) in its classical formulations was applied to develop static models describing intermolecular hydrogen bond (HB) and isotope effect in the gas phase and using solvent reaction field reproduced by Polarizable Continuum Model (PCM). The Atoms in Molecules (AIM) theory enabled electronic structure and molecular topology study. The Symmetry-Adapted Perturbation Theory (SAPT) was used for energy decomposition in the dimers of n-octanol. Finally, time-evolution methods, namely classical molecular dynamics (MD) and Car-Parrinello Molecular Dynamics (CPMD) were employed to shed light onto dynamical nature of liquid n-octanol with emphasis put on metric and vibrational features. As a reference, CPMD gas phase results were applied. Nuclear quantum effects were included using Path Integral Molecular Dynamics (PIMD) and a posteriori method by solving vibrational Schrödinger equation. The latter applied procedure allowed to study the deuterium isotope effect.     

Car-Parrinello simulation of an O-H stretching envelope and potential of mean force of an intramolecular hydrogen bonded system: application to a Mannich base in solid state and in vacuum

Car-Parrinello molecular dynamics (CPMD) study was performed for an anharmonic system-an intramolecularly hydrogen bonded Mannich-base-type compound, 4,5-dimethyl-2(N,N-dimethylaminemethyl)phenol, to investigate the vibrational spectrum associated with the O-H stretching. Calculations were carried out for the solid state and for an isolated molecule. The classical CPMD simulation was performed and then the proton potential snapshots were extracted from the trajectory. The vibrational Schrodinger equation for the snapshots was solved numerically, and the (O-H) envelope was calculated as a superposition of the 0-->1 transitions. The potential of mean force for the proton stretching mode was calculated from the proton vibrational eigenfunctions and eigenvalues incorporating statistical sampling, nuclear quantum effects, and effects of the environment. Perspectives for application of the presented methodology in the computational support of biocatalysis are given in the study.     

Car-Parrinello molecular dynamics study of anharmonic systems: a Mannich base in solution

A Car-Parrinello molecular dynamics study was performed for 4,5-dimethyl-2-(N,N-dimethylaminomethyl)phenol, a Mannich base, to investigate the vibrational properties in solution of its intramolecular hydrogen bond. The dynamic behavior of this hydrogen-bonded system was investigated using an explicit solvent model. Addition of a nonpolar solvent permitted inclusion of delicate environmental effects on the strongly anharmonic system which was studied from first principles. Molecular dynamics and a posteriori quantization of the O-H motion were applied to reproduce the vibrational features of the O-H stretching mode. Consistent application of Car-Parrinello dynamics based on the density functional theory with subsequent solution of the vibrational Schr?dinger equation for the O-H stretching motion offers an effective method for strongly anharmonic systems, and this is supported by the comparison of the results with experimental spectra. As a further element of the intramolecular hydrogen bond study, the effects of deuteration were taken into account and a successful application of the O-H stretching mode quantization technique to the liquid phase is demonstrated. This provides a valuable computational methodology for investigations incorporating nuclear quantum effects in the liquid phase and enzyme active centers and can be used to investigate numerous systems that are not readily susceptible to experimental analysis.     

Investigations of the hydrogen bond in the crystals of tropolone and thiotropolone via car-parrinello and path integral molecular dynamics

Car-Parrinello and path integrals molecular dynamics (CPMD and PIMD) simulations were carried out for the 10π-electron aromatic systems: 2-hydroxy-2,4,6-cycloheptatrien-1-one, commonly known as Tropolone (I) and 2-hydroxy-2,4,6-cycloheptatriene-1-thione, called Thiotropolone (II) in vacuo and in the solid state. The extremely fast proton transfer (FPT) and “prototropy” tautomerism in the keto-enol (thione-enethiol) systems have been analyzed on the basis of CPMD and PIMD methods level. Comparisons of two-dimensional (2D) free-energy landscapes of reaction coordinate δ-parameter and R_O…O or R_O…S distances shows that the OH… tautomer to be more favorable in the Thiotropolone. The hydrogen between the oxygen and the sulfur atoms adopts a starkly asymmetrical position in the double potential well. The values of the energy barriers for the FPT were calculated and suggested a strong hydrogen bond with low barrier for FPT mechanism. These studies and the 2D average index of π-delocalization 〈λ〉 landscape of time evolutions of R_O1…O2 and R_C7O2 or R_C7S1 distances for the both crystals indicate that hydrogen bonds in the crystals of Tropolone (I) and Thiotropolone (II) have characteristic properties for the type of bonding model resonance-assisted hydrogen bonds and also low-barrier hydrogen bonds. In the crystal of the Thiotropolone (II), we found the hydrogen bond O H…S existing without the equilibrium of the two tautomers whereas in the crystal of the Tropolone (I) has been confirmed the hydrogen bond O H…O existing with the equilibrium of the two tautomers. It was also found the significant differences in frequency, speed, and the image of the FPT in the studied crystals. © 2018 Wiley Periodicals, Inc.     

Direct assessment of quantum nuclear effects on hydrogen bond strength by constrained-centroid ab initio path integral molecular dynamics

The impact of quantum nuclear effects on hydrogen (H-) bond strength has been inferred in earlier work from bond lengths obtained from path integral molecular dynamics (PIMD) simulations. To obtain a direct quantitative assessment of such effects, we use constrained-centroid PIMD simulations to calculate the free energy changes upon breaking the H-bonds in dimers of HF and water. Comparing ab initio simulations performed using PIMD and classical nucleus molecular dynamics (MD), we find smaller dissociation free energies with the PIMD method. Specifically, at 50 K, the H-bond in (HF)(2) is about 30% weaker when quantum nuclear effects are included, while that in (H(2)O)(2) is about 15% weaker. In a complementary set of simulations, we compare unconstrained PIMD and classical nucleus MD simulations to assess the influence of quantum nuclei on the structures of these systems. We find increased heavy atom distances, indicating weakening of the H-bond consistent with that observed by direct calculation of the free energies of dissociation.     

Proton dynamics in the strong chelate hydrogen bond of crystalline picolinic acid N-oxide. A new computational approach and infrared, raman and INS study

Infrared, Raman and INS spectra of picolinic acid N-oxide (PANO) were recorded and examined for the location of the hydronic modes, particularly O-H stretching and COH bending. PANO is representative of strong chelate hydrogen bonds (H-bonds) with its short O...O distance (2.425 A). H-bonding is possibly well-characterized by diffraction, NMR and NQR data and calculated potential energy functions. The analysis of the spectra is assisted by DFT frequency calculations both in the gas phase and in the solid state. The Car-Parrinello quantum mechanical solid-state method is also used for the proton dynamics simulation; it shows the hydron to be located about 99% of time in the energy minimum near the carboxylic oxygen; jumps to the N-O acceptor are rare. The infrared spectrum excels by an extended absorption (Zundel's continuum) interrupted by numerous Evans transmissions. The model proton potential functions on which the theories of continuum formation are based do not correspond to the experimental and computed characteristics of the hydrogen bond in PANO, therefore a novel approach has been developed; it is based on crystal dynamics driven hydronium potential fluctuation. The envelope of one hundred 0 --> 1 OH stretching transitions generated by molecular dynamics simulation exhibits a maximum at 1400 cm-1 and a minor hump at approximately 1600 cm-1. These positions square well with ones predicted for the COH bending and OH stretching frequencies derived from various one- and two-dimensional model potentials. The coincidences with experimental features have to be considered with caution because the CPMD transition envelope is based solely on the OH stretching coordinate while the observed infrared bands correspond to heavily mixed modes as was previously shown by the normal coordinate analysis of the IR spectrum of argon matrix isolated PANO, the present CPMD frequency calculation and the empirical analysis of spectra. The experimental infrared spectra show some unusual characteristics such as large temperature effects on the intensity of some bands, thus presenting a challenge for theoretical band shape treatments. Our calculations clearly show that the present system is characterized by an asymmetric single well potential with no large amplitudes in the hydronium motion, which extends the existence of Zundel-type spectra beyond the established set of hydrogen bonds with large hydronic vibrational amplitudes.     

(1Z,3Z)-1,4-Di(pyridin-2-yl)buta-1,3-diene-2,3-diol: the planar highly conjugated symmetrical enediol with multiple intramolecular hydrogen bonds

1H, (13)C, and (15)N NMR spectral data show that in chloroform solution (1Z,3Z)-1,4-di(pyridin-2-yl)buta-1,3-diene-2,3-diol, OO, is in ca. 9:1 equilibrium with (3Z)-3-hydroxy-1,4-di(pyridin-2-yl)but-3-en-2-one, OK, while no 1,4-di(pyridin-2-yl)-2,3-butanedione, KK, was detected. The species present in the tautomeric mixture were identified by comparing their experimental chemical shifts with those known for similar compounds as well as with the theoretically calculated (GIAO-HF/DFT) values. Ab initio calculations show that OK and especially the highly conjugated OO forms are preferred in the tautomeric mixtures both in vacuo and in chloroform solution. Comparison of experimental (Arrhenius) and calculated (ab initio) energies of OK shows that the MP2/6-31G//RHF/6-31G method gives the most precise results. There are one and two strong O-H.N hydrogen bonds present in the molecules of the former and latter compound, respectively. Other tautomeric forms, e.g., dienaminedione [(1Z,4Z)-1,4-di[pyridin-2(1H)-ylidene]butane-2,3-dione], and their rotamers were found to have higher energies. The single-crystal X-ray diffraction studies show that dienediol OO is the only tautomeric form present in the crystal at 173 K. Its almost perfectly planar molecule is stabilized by two strong intramolecular O-H.N hydrogen bonds.     

Competition of Intra- and Intermolecular Forces in Anthraquinone and Its Selected Derivatives

Intra- and intermolecular forces competition was investigated in the 9,10-anthraquinone (1) and its derivatives both in vacuo and in the crystalline phase. The 1,8-dihydroxy-9,10-anthraquinone (2) and 1,8-dinitro-4,5-dihydroxy-anthraquinone (3) contain Resonance-Assisted Hydrogen Bonds (RAHBs). The intramolecular hydrogen bonds properties were studied in the electronic ground and excited states employing Møller-Plesset second-order perturbation theory (MP2), Density Functional Theory (DFT) method in its classical formulation as well as its time-dependent extension (TD-DFT). The proton potential functions were obtained via scanning the OH distance and the dihedral angle related to the OH group rotation. The topological analysis was carried out on the basis of theories of Atoms in Molecules (AIM—molecular topology, properties of critical points, AIM charges) and Electron Localization Function (ELF—2D maps showing bonding patterns, calculation of electron populations in the hydrogen bonds). The Symmetry-Adapted Perturbation Theory (SAPT) was applied for the energy decomposition in the dimers. Finally, Car–Parrinello molecular dynamics (CPMD) simulations were performed to shed light onto bridge protons dynamics upon environmental influence. The vibrational features of the OH stretching were revealed using Fourier transformation of the autocorrelation function of atomic velocity. It was found that the presence of OH and NO2 substituents influenced the geometric and electronic structure of the anthraquinone moiety. The AIM and ELF analyses showed that the quantitative differences between hydrogen bonds properties could be neglected. The bridged protons are localized on the donor side in the electronic ground state, but the Excited-State Intramolecular Proton Transfer (ESIPT) was noticed as a result of the TD-DFT calculations. The hierarchy of interactions determined by SAPT method indicated that weak hydrogen bonds play modifying role in the organization of these crystal structures, but primary ordering factor is dispersion. The CPMD crystalline phase results indicated bridged proton-sharing in the compound 2.     

Identification of active sites of biomolecules. 1. Methyl-alpha-mannopyranoside and Fe(III)

Car-Parrinello molecular dynamics (CPMD) simulations, DFT chemical reactivity index calculations, and mass spectrometric measurements are combined in an integrated effort to elucidate the details of the coordination of a transition-metal ion to a carbohydrate. The impact of the interaction with the FeIII ion on the glycosidic linkage conformation of methyl-alpha-d-mannopyranoside is studied by classical molecular dynamics (MD) and CPMD simulations. This study shows that FeIII interacts with specific hydroxyl oxygen atoms of the carbohydrate, affecting the ground state carbohydrate conformation. These conformational details are discussed in terms of a set of supporting experiments involving electrospray ionization mass spectrometry, and CPMD simulations clearly indicate that the specific conformational preference is due to intramolecular hydrogen bonding. Classical MD simulations proved insensitive to these important chemical properties. Thus, we demonstrate the importance of chemical reactivity calculations and CPMD simulations in predicting the active sites of biological molecules toward metal cations.     

Convenient synthesis of 5-benzyl-2,3,5-trichloro-4,4-dimethoxycyclopent-2-en-1-one and some its reactions

5-Benzyl-2,3,5-trichloro-4,4-dimethoxycyclopent-2-en-1-one was synthesized by reaction of 1,2,3,4-tetrachloro-5,5-dimethoxycyclopentadiene with benzyl alcohol in methylene chloride in the presence of powdered sodium hydroxide and benzyltrimethylammonium chloride as phase-transfer catalyst. Deprotection of the title compound gave 2-benzyl-2,4,5-trichlorocyclopent-4-ene-1,3-dione which was subjected to unusual intramolecular carbocyclization initiated by molecular iodine.     

Aromaticity-controlled tautomerism and resonance-assisted hydrogen bonding in heterocyclic enaminone-iminoenol systems

The contribution of aromaticity and intramolecular hydrogen bonding to relative stability, for a set of (1H-azahetero-2-ylidene)-acetaldehyde and 2-azahetero-2-yl-ethanol tautomeric pairs, has been investigated by means of quantum chemical DFT and ab initio methods up to the MP4(SDTQ)/AUG-cc-pVDZ and MP2/AUG-cc-pVTZ levels of theory. It is found that the relative energy of the tautomers is governed by the change in the degree of heterocycle aromaticity upon intramolecular hydrogen transfer. An analysis of geometrical parameters of a hydrogen-bonded system reveals a clear relationship between the aromaticity of the heterocycle, the conjugation in a resonant spacer, and the strengths of the intramolecular hydrogen bonds. This allows the conclusion to be drawn that intramolecular N-H...O and O-H...N hydrogen bonds formed are found to be resonance-assisted and their strength is dependent on the pi-donating/accepting properties of the heterocycle. On the basis of the results of the calculations, a simple model describing the mechanism of resonance assistance of hydrogen bonding has been suggested.     

Nuclear quantum and H/D isotope effects on three-centered bonding diborane: Path integral molecular dynamics simulations

Nuclear quantum and H/D isotope effects of bridging and terminal hydrogen atoms of diborane (B₂H₆) molecules were systematically studied by classical ab initio molecular dynamics (CLMD) and ab initio path integral molecular dynamics (PIMD) simulations with BHandHLYP/6-31++G** level of theory at room temperature (298.15 K). Calculated results clearly show that H/D isotope effect appears in the distribution of hydrogen (deuterium) of B₂H₆ (B₂D₆). Geometry of B₂H₆ also plays a significant role in the nuclear quantum effect proved by PIMD simulations, but slightly deviated from its equilibrium structure when simulated via CLMD simulation. The bond lengths between boron atoms R (B1 … B2) and the bridging hydrogen atoms R_HH (H_B1 … H_B2) of the B₂H₆ molecule obtained from PIMD simulations are slightly longer than those of the deuterated form of the diborane (B₂D₆) molecule. The principal component analysis (PCA) was also employed to distinguish the important modes of bridging hydrogen as related to the nuclear quantum and H/D isotope effects. The highest level of contribution obtained from PCA of PIMD simulations is bending, while various mixed vibrations with less contribution were also found. Therefore, the nuclear quantum and H/D isotope effects need to be taken into account for a better understanding of diborane geometry.     

Chiral 5-(diphenylphosphanyl)-1,2,3,4-tetrahydroacridines: new N,P-ligands for asymmetric catalysis

Three chiral 5-(diphenylphosphanyl)-1,2,3,4-tetrahydroacridines, as first representative examples of a new class of chiral N,P-ligands were prepared from (+)-nopinone, (+)-camphor and 5α-androst-2-en-17-one. These ligands have been assessed in the enantioselective palladium-catalysed allylic substitution of 1,3-diphenylprop-2-enyl acetate with dimethyl malonate. Enantioselectivity up to 74% has been obtained.     

Isotope effects in ice Ih: a path-integral simulation

Ice Ih has been studied by path-integral molecular dynamics simulations, using the effective q-TIP4P/F potential model for flexible water. This has allowed us to analyze finite-temperature quantum effects in this solid phase from 25 to 300 K at ambient pressure. Among these effects we find a negative thermal expansion of ice at low temperatures, which does not appear in classical molecular dynamics simulations. The compressibility derived from volume fluctuations gives results in line with experimental data. We have analyzed isotope effects in ice Ih by considering normal, heavy, and tritiated water. In particular, we studied the effect of changing the isotopic mass of hydrogen on the kinetic energy and atomic delocalization in the crystal as well as on structural properties such as interatomic distances and molar volume. For D(2)O ice Ih at 100 K we obtained a decrease in molar volume and intramolecular O-H distance of 0.6% and 0.4%, respectively, as compared to H(2)O ice.     

Monitoring selected hydrogen bonds in crystal hydrates of amino acid salts: combining variable-temperature single-crystal X-ray diffraction and polarized Raman spectroscopy

Predicting behaviour of hydrogen bonds with varying temperature, in particular-correlating donor-acceptor distances in the O-H···O hydrogen bonds with the frequencies of O-H stretching vibrations is important for understanding dynamics of biomolecules and phase transitions in crystals. A commonly used correlation suggested earlier in the literature is based on statistical analysis of different compounds [A. Novak, Structure and Bonding, 1974, 18, 177; K. Nakamoto, M. Margoshes, R. E. Rundle, J. Am. Chem. Soc., 1955, 77, 6480]. The present study is a rare example when correlations between geometry and energy parameters have been found for selected individual hydrogen bonds in the same crystalline compound at multiple temperatures. The properties of several types of O-H···O hydrogen bonds in bis(DL-serinium) oxalate dihydrate and DL-alaninium semi-oxalate monohydrate have been studied by a combination of variable-temperature single-crystal X-ray diffraction and polarized Raman spectroscopy. The changes in the hydrogen bonds geometry could be compared with the changes of the corresponding spectral modes. The correlation suggested by Novak is roughly followed, better for medium and weak, than for short hydrogen bonds. Fine details of spectral changes differ for individual bonds. The way how H-bonds are affected by cooling depends on their environment in the crystal structure. Short O-H···O hydrogen bonds in bis(DL-serinium) oxalate dihydrate expand or remain almost unchanged on cooling, whereas in DL-alaninium semi-oxalate monohydrate all strong H-bonds are compressed under these conditions. The distortion of individual hydrogen bonds on temperature variations is correlated with the anisotropy of lattice strain.     

Comparative estimation of the energies of intramolecular C-H…O,N-H…O,and O-H…O hydrogen bonds according to the QTAIM analysis and NMR spectroscopy data

The energies of intramolecular C-H…O, N-H…O, and O-H…O hydrogen bonds in model compounds are empirically estimated based on the values of the hydrogen bond induced weak-field shift of the bridging hydrogen atom signal in the ¹H NMR spectrum. It is supported by a theoretical estimation of these energies based on the electron density value at the hydrogen bond critical point calculated within the QTAIM method. Good agreement between the empirical and theoretical estimates is found, which gives evidence of their reliability. It is shown that from the standpoint of their strength the intramolecular N-H…O and O-H…O hydrogen bonds can be classified as moderate whereas the intramolecular C-H…O hydrogen bonds must be classified as very weak interactions similar in their energy significance to van der Waals interactions.     

Direct observation of the substitution effects on the hydrogen bridge dynamics in selected Schiff bases--a comparative molecular dynamics study

We have studied substituent effects on the properties of the intramolecular hydrogen bond of some ortho-hydroxy Schiff bases using density functional theory (DFT) based first-principle molecular dynamics (FPMD) and path integral molecular dynamics. The studied compounds possess a strong intramolecular hydrogen bond (r((O???N)) ≤ 2.6 A?), which can be tuned by substitution to either (i) enhance the basicity of the acceptor moiety by induction effects or (ii) decrease the hydrogen bond length through steric repulsion. DFT calculations and FPMD were employed to investigate structural and dynamical properties of the selected molecules, while quantum effects on the structural properties were assessed using path integral FPMD. The simulations were performed in vacuo and in the solid state to study the influence of the environment on the hydrogen bond and spectroscopic properties. We give computational support to the suggestion that induction effects are less effective to tune the intramolecular hydrogen bond properties of the discussed ortho-hydroxy Schiff bases than the steric or the environmental effects.