Density functional study of the proton transfer effect on vibrations of strong (short) intermolecular O-H...N/O-...H-N+ hydrogen bonds in aprotic solvents

Combined X-ray Crystallographic, IR/Raman Spectroscopic, and Periodic DFT Investigations of New Multicomponent Crystalline Forms of Anthelmintic Drugs: A Case Study of Carbendazim Maleate

Synthesis of multicomponent solid forms is an important method of modifying and fine-tuning the most critical physicochemical properties of drug compounds. The design of new multicomponent pharmaceutical materials requires reliable information about the supramolecular arrangement of molecules and detailed description of the intermolecular interactions in the crystal structure. It implies the use of a combination of different experimental and theoretical investigation methods. Organic salts present new challenges for those who develop theoretical approaches describing the structure, spectral properties, and lattice energy Elatt. These crystals consist of closed-shell organic ions interacting through relatively strong hydrogen bonds, which leads to Elatt > 200 kJ/mol. Some technical problems that a user of periodic (solid-state) density functional theory (DFT) programs encounters when calculating the properties of these crystals still remain unsolved, for example, the influence of cell parameter optimization on the Elatt value, wave numbers, relative intensity of Raman-active vibrations in the low-frequency region, etc. In this work, various properties of a new two-component carbendazim maleate crystal were experimentally investigated, and the applicability of different DFT functionals and empirical Grimme corrections to the description of the obtained structural and spectroscopic properties was tested. Based on this, practical recommendations were developed for further theoretical studies of multicomponent organic pharmaceutical crystals.

Simplified calculation approaches designed to reproduce the geometry of hydrogen bonds in molecular complexes in aprotic solvents

The impact of the environment onto the geometry of hydrogen bonds can be critically important for the properties of the questioned molecular system. The paper reports on the design of calculation approaches capable to simulate the effect of aprotic polar solvents on the geometric and NMR parameters of intermolecular hydrogen bonds. A hydrogen fluoride and pyridine complex has been used as the main model system because the experimental estimates of these parameters are available for it. Specifically, F-H, F⋯N, and H-N distances, the values of ¹⁵N NMR shift, and spin-spin coupling constants ¹J(¹⁹F¹H), ^1hJ(¹H¹⁵N), and ^2hJ(¹⁹F¹⁵N) have been analyzed. Calculation approaches based on the gas-phase and the Polarizable Continuum Model (PCM) approximations and their combinations with geometric constraints and additional noncovalent interactions have been probed. The main result of this work is that the effect of an aprotic polar solvent on the geometry of a proton-donor⋯H⋯proton-acceptor complex cannot be reproduced under the PCM approximation if no correction for solvent-solute interactions is made. These interactions can be implicitly accounted for using a simple computational protocol.

Organic co-crystals of 1,3-bis(4-pyridyl)azulene with a series of hydrogen-bond donors

1,3-Bis(4-pyridyl)azulene has been employed as a hydrogen bond acceptor to construct two-component organic cocrystals.

NMR Study of Solvation Effect on the Geometry of Proton-Bound Homodimers of Increasing Size

Hydrogen bond geometries in the proton-bound homodimers of quinoline and acridine derivatives in an aprotic polar solution have been experimentally studied using ¹H NMR at 120 K. The reported results show that an increase of the dielectric permittivity of the medium results in contraction of the N···N distance. The degree of contraction depends on the homodimer's size and its substituent-specific solvation features. Neither of these effects can be reproduced using conventional implicit solvent models employed in computational studies. In general, the N···N distance in the homodimers of pyridine, quinoline, and acridine derivatives decreases in the sequence gas phase > solid state > polar solvent.

New multi-component solid forms of anti-cancer drug Erlotinib: role of auxiliary interactions in determining a preferred conformation

Erlotinib is a BCS (biopharmaceutical classification system) class II drug used for the treatment of non-small cell lung cancer. There is an urgent need to obtain new solid forms of higher solubility to improve the bioavailability of the API (active pharmaceutical ingredient). In this context, cocrystals with urea, succinic acid, and glutaric acid and salts with maleic acid, adipic acid, and saccharin were prepared via wet granulation and solution crystallizations. Crystal structures of the free base (Z' = 2), cocrystals of erlotinib-urea (1:1), erlotinib-succinic acid monohydrate (1:1:1), erlotinib-glutaric acid monohydrate (1:1:1) and salts of erlotinib-adipic acid adipate (1:0.5:0.5) are determined and their hydrogen-bonding patterns are analyzed. Self recognition via the (amine) N-H...N (pyridine) hydrogen bond between the API molecules is replaced by several heterosynthons such as acid-pyridine, amide-pyridine and carboxylate-pyridinium in the new binary systems. Auxiliary interactions play an important role in determining the conformation of the API in the crystal. FT-IR spectroscopy is used to distinguish between the salts and cocrystals in the new multi-component systems. The new solid forms are characterized by powder X-ray diffraction (PXRD) and differential scanning calorimetry (DSC) to confirm their unique phase identity.

Quantification of Conventional and Nonconventional Charge-Assisted Hydrogen Bonds in the Condensed and Gas Phases

Charge-assisted hydrogen bonds (CAHBs) play critical roles in many systems from biology through to materials. In none of these areas has the role and function of CAHBs been explored satisfactorily because of the lack of data on the energy of CAHBs in the condensed phases. We have, for the first time, quantified three types of CAHBs in both the condensed and gas phases for 1-(2'-hydroxylethyl)-3-methylimidazolium acetate ([C2OHmim][OAc]). The energy of conventional OH···[OAc](-) CAHBs is ∼10 kcal·mol(-1), whereas nonconventional C(sp2)H···[OAc](-) and C(sp3)H···[OAc](-) CAHBs are weaker by ∼5-7 kcal·mol(-1). In the gas phase, the strength of the nonconventional CAHBs is doubled, whereas the conventional CAHBs are strengthened by <20%. The influence of cooperativity effects on the ability of the [OAc](-) anion to deprotonate the imidazolium cation is evaluated. The ability to quantify CAHBs in the condensed phase on the basis of easier accessible gas-phase estimates is highlighted.

The applicability of the dimeric heterosynthon concept to molecules with equivalent binding sites. A DFT study of crystalline urea–H2O2

The limited applicability of the dimeric heterosynthon concept to a two-component urea–H₂O₂crystal is reported.

Measurement of proton tunneling in short hydrogen bonds in single crystals of 3,5 pyridinedicarboxylic acid using nuclear magnetic resonance spectroscopy

In this Letter, we present NMR spin-lattice and relaxometry data for proton transfer in one of the shortest known N-H⋯O hydrogen bonds in a single crystal of 3,5 pyridinedicarboxylic acid (35PDCA). It is widely believed that proton transfer by quantum tunneling does not occur in short hydrogen bonds since the ground state energy level lies above the potential barrier, yet these data show a temperature independent, proton tunneling rate below 77 K and a clear deviation from classical dynamics below 91 K. This study therefore suggests that proton tunneling occurs in all hydrogen bonds at low temperature and the crossover temperature to classical hopping must be determined when evaluating whether proton tunneling persists at higher temperature, for example in enzyme catalysis under physiological conditions.

Hydrogen-bonding patterns in a series of multi-component molecular solids formed by 2,3,5,6-tetramethylpyrazine with selected carboxylic acids

Hydrogen-bonding patterns in a series of cocrystals of 2,3,5,6-tetramethylpyrazine with selected carboxylic acids are discussed.

How strong is hydrogen bonding in ionic liquids? Combined X-ray crystallographic, infrared/Raman spectroscopic, and density functional theory study

Hydrogen bonding in ionic liquids based on the 1-(2'-hydroxylethyl)-3-methylimidazolium cation ([C₂OHmim](+)) and various anions ([A](-)) of differing H-bond acceptor strength, viz. hexafluorophosphate [PF6](-), tetrafluoroborate [BF₄](-), bis(trifluoromethanesulfonimide) [Tf₂N](-), trifluoromethylsulfonate [OTf](-), and trifluoroacetate [TFA](-), was studied by a range of spectroscopic and computational techniques and, in the case of [C₂OHmim][PF6], by single crystal X-ray diffraction. The first quantitative estimates of the energy (E(HB)) and the enthalpy (-ΔH(HB)) of H-bonds in bulk ILs were obtained from a theoretical analysis of the solid-state electron-density map of crystalline [C₂OHmim][PF6] and an analysis of the IR spectra in crystal and liquid samples. E(HB) for OH···[PF6](-) H-bonds amounts to ~3.4-3.8 kcal·mol(-1), whereas weaker H-bonds (2.8-3.1 kcal·mol(-1)) are formed between aromatic C2H group of imidazolium ring and the [PF6](-) anion. The enthalpy of the OH···[A](-) H-bonds follows the order: [PF6] (2.4 kcal·mol(-1)) < [BF₄] (3.3 kcal·mol(-1)) < [Tf₂N] (3.4 kcal·mol(-1)) < [OTf] (4.7 kcal·mol(-1)l) < [TFA] (6.2 kcal·mol(-1)). The formation of aggregates of self-associated [C₂OHmim](+) cations is present in liquid [C₂OHmim][PF6], [C₂OHmim][BF₄], and [C₂OHmim][Tf₂N], with the energy of the OH···OH H-bonds amounting to ~6 kcal·mol(-1). Multiple secondary interactions in the bulk ILs influence their structure, vibrational spectra, and H-bond strength. In particular, these interactions can blue-shift the stretching frequencies of the CH groups of the imidazolium ring in spite of red-shifting CH···[A](-) H-bonds. They also weaken the H-bonding in the IL relative to the isolated ion pairs, with these anticooperative effects amounting to ca. 50% of the E(HB) value.

Is there a simple way to reliable simulations of infrared spectra of organic compounds?

To assess the ability of the quantum-chemical computations to reproduce the experimental relative intensities in the infrared (IR) spectra of both the gas- and condensed-phase systems, the hybrid DFT functional B3LYP has been applied to simulation of IR spectra for species containing from three to twelve first- or second-row atoms, both in the gas phase and in CCl4 solutions. The results demonstrate that B3LYP, combined with the highly compact double-ζ basis set 6-31+G* and "scaled quantum mechanics" techniques, offers excellent quantitative performance in the calculations of relative IR intensities and frequencies (ν ≤ 2200 cm(-1)) for the bands of vibrations of medium-size isolated molecules, whereas it produces unsatisfactory results for the solutions of the same species. Neither larger basis sets nor implicit treatment of the media effects improve the agreement of the simulated spectra with the condensed-phase experiment. At the same time, some preliminary results suggest that explicit modeling of media effects could offer better quality of the IR spectral simulations for the condensed-phase systems.

Energy Analysis of Competing Non-Covalent Interaction in 1:1 and 1:2 Adducts of Collidine with Benzoic Acids by Means of X-Ray Diffraction

The hydrogen bond pattern and the types of non-covalent interactions in the crystals of the 1:1 and 1:2 adducts of 2,4,6-trimethylpyridine and benzoic acids are studied using high-resolution X-ray diffraction. The geometries of the hydrogen bonds are estimated using a combined XRD/DFT approach that provides the geometrical parameters within the margin of error of neutron diffraction studies. The energies of the non-covalent interactions are estimated on the base of the experimental electron density distribution function. It is shown that the structures of the adducts are governed by the NOH and OHO hydrogen bonds. In turn, C-H...O contacts and stacking interactions define the packing of the adducts in the crystal. On the other hand, it is important to note that the latter interactions affect the competition of the former hydrogen bonds in some 1:2 adducts.

Acridine – a Promising Fluorescence Probe of Non-Covalent Molecular Interactions

Fluorescence and absorption spectral parameters of acridine have been studied in solution at room temperature in the presence of different proton donors aiming to inspect whether or not acridine can be used as a spectroscopic probe suitable to measure the geometry of hydrogen bonds under different conditions. It has been shown that the most appropriate spectral parameter is the position of fluorescence maximum that changes heavily upon a contraction of the N···H distance. Presumably, also the intensity of the maximum strongly depends on the hydrogen bond geometry. These two parameters can be used to establish two independent, mutually complementary correlations connecting the spectral manifestations and the geometry of hydrogen bond.

Intermolecular hydrogen bond energies in crystals evaluated using electron density properties: DFT computations with periodic boundary conditions

The hydrogen bond (H-bond) energies are evaluated for 18 molecular crystals with 28 moderate and strong O-H···O bonds using the approaches based on the electron density properties, which are derived from the B3LYP/6-311G** calculations with periodic boundary conditions. The approaches considered explore linear relationships between the local electronic kinetic G(b) and potential V(b) densities at the H···O bond critical point and the H-bond energy E(HB). Comparison of the computed E(HB) values with the experimental data and enthalpies evaluated using the empirical correlation of spectral and thermodynamic parameters (Iogansen, Spectrochim. Acta Part A 1999, 55, 1585) enables to estimate the accuracy and applicability limits of the approaches used. The V(b)-E(HB) approach overestimates the energy of moderate H-bonds (E(HB) < 60 kJ/mol) by ~20% and gives unreliably high energies for crystals with strong H-bonds. On the other hand, the G(b)-E(HB) approach affords reliable results for the crystals under consideration. The linear relationship between G(b) and E(HB) is basis set superposition error (BSSE) free and allows to estimate the H-bond energy without computing it by means of the supramolecular approach. Therefore, for the evaluation of H-bond energies in molecular crystals, the G(b) value can be recommended to be obtained from both density functional theory (DFT) computations with periodic boundary conditions and precise X-ray diffraction experiments.

Evaluation of coupling terms between intra- and intermolecular vibrations in coarse-grained normal-mode analysis: does a stronger acid make a stiffer hydrogen bond?

Using theory of harmonic normal-mode vibration analysis, we developed a procedure for evaluating the anisotropic stiffness of intermolecular forces. Our scheme for coarse-graining of molecular motions is modified so as to account for intramolecular vibrations in addition to relative translational/rotational displacement. We applied this new analytical scheme to four carboxylic acid dimers, for which coupling between intra- and intermolecular vibrations is crucial for determining the apparent stiffness of the intermolecular double hydrogen bond. The apparent stiffness constant was analyzed on the basis of a conjunct spring model, which defines contributions from true intermolecular stiffness and molecular internal stiffness. Consequently, the true intermolecular stiffness was in the range of 43-48 N m(-1) for all carboxylic acids studied, regardless of the molecules' acidity. We concluded that the difference in the apparent stiffness can be attributed to differences in the internal stiffness of the respective molecules.

H-bond network in amino acid cocrystals with H2O or H2O2. The DFT study of serine-H2O and serine-H2O2

The structure, IR spectrum, and H-bond network in the serine-H(2)O and serine-H(2)O(2) crystals were studied using DFT computations with periodic boundary conditions. Two different basis sets were used: the all-electron Gaussian-type orbital basis set and the plane wave basis set. Computed frequencies of the IR-active vibrations of the titled crystals are quite different in the range of 10-100 cm(-1). Harmonic approximation fails to reproduce IR active bands in the 2500-2800 frequency region of serine-H(2)O and serine-H(2)O(2). The bands around 2500 and 2700 cm(-1) do exist in the anharmonic IR spectra and are caused by the first overtone of the OH bending vibrations of H(2)O and a combination vibration of the symmetric and asymmetric bendings of H(2)O(2). The quantum-topological analysis of the crystalline electron density enables us to describe quantitatively the H-bond network. It is much more complex in the title crystals than in a serine crystal. Appearance of water leads to an increase of the energy of the amino acid-amino acid interactions, up to ~50 kJ/mol. The energy of the amino acid-water H-bonds is ~30 kJ/mol. The H(2)O/H(2)O(2) substitution does not change the H-bond network; however, the energy of the amino acid-H(2)O(2) contacts increases up to 60 kJ/mol. This is caused by the fact that H(2)O(2) is a much better proton donor than H(2)O in the title crystals.

Geometry and spectral properties of the protonated homodimer of pyridine in the liquid and solid states. A combined NMR, X-ray diffraction and inelastic neutron scattering study

The structure and spectral signatures of the protonated homodimer of pyridine in its complex with a poorly coordinating anion have been studied in solution in CDF(3)/CDClF(2) down to 120 K and in a single crystal. In both phases, the hydrogen bond is asymmetric. In the solution, the proton is involved in a fast reversible transfer that determines the multiplicity of NMR signals and the sign of the primary H/D isotope effect of --0.95 ppm. The proton resonates at 21.73 ppm that is above any value reported in the past and is indicative of a very short hydrogen bond. By combining X-ray diffraction analysis with model computations, the position of the proton in the crystal has been defined as d(N-H) = 1.123 Å and d(H···N) = 1.532 Å. The same distances have been estimated using a (15)N NMR correlation. The frequency of the protonic out-of-plane bending mode is 822 cm(-1) in agreement with Novak's correlation.

A model proton-transfer system in the condensed phase: NH4(+)OOH(-), a crystal with short intermolecular H-bonds

The crystal structure of NH(4)(+)OOH(-) is determined from single-crystal x-ray data obtained at 150 K. The crystal belongs to the space group P2(1)/c and has four molecules in a unit cell. The structure consists of discrete NH(4)(+) and OOH(-) ions. The OOH(-) ions are linked by short hydrogen bonds (2.533 Å) to form parallel infinite chains. The ammonium ions form links between these chains (the N⋯O distances vary from 2.714 to 2.855 Å) giving a three-dimensional network. The harmonic IR spectrum and H-bond energies are computed at the Perdew-Burke-Ernzerhof (PBE)/6-31G(∗∗) level with periodic boundary conditions. A detailed analysis of the shared (bridging) protons' dynamics is obtained from the CPMD simulations at different temperatures. PBE functional with plane-wave basis set (110 Ry) is used. At 10 K the shared proton sits near the oxygen atom, only a few proton jumps along the chain are detected at 70 K while at 270 K numerous proton jumps exist in the trajectory. The local-minimum structure of the space group Cc is localized. It appears as a result of proton transfer along a chain. This process is endothermic (∼2 kJ/mol) and is described as P2(1)/c↔2Cc. The computed IR spectrum at 10 K is close to the harmonic one, the numerous bands appear at 70 K while at 270 K it shows a very broad absorption band that covers frequencies from about 1000 to 3000 cm(-1). The advantages of the NH(4)(+)OOH(-) crystal as a promising model for the experimental and DFT based molecular dynamics simulation studies of proton transfer along the chain are discussed.