Estimation of the intramolecular O-H···O═C hydrogen bond energy via the molecular tailoring approach. Part I: aliphatic structures

Molecular Tailoring Approach for the Direct Estimation of Individual Noncovalent Interaction Energies in Molecular Systems

The noncovalent interactions (NCIs) are omnipresent in chemistry, physics, and biology. The study of such interactions offers insights into various physicochemical phenomena. Some indirect approaches proposed in the literature for exploring the NCIs are briefly reviewed in Section 1 of this Perspective. These include: (i) Shift in the stretching frequency of an X-Y bond involved in X-Y···Z interaction. (ii) Topological analysis of molecular electron density. (iii) Empirical equations derived employing experimental and theoretical quantities. However, a direct method for estimating individual intramolecular/intermolecular interaction energies has been conspicuous by its absence from the literature. We have developed a molecular tailoring approach (MTA)-based method enabling a direct and reliable estimation of the energy of intra- as well as intermolecular interactions. This method offers a direct and reliable estimation of these interactions, in particular of the hydrogen bonds (HB) in molecules/weakly bound clusters along with the respective cooperativity contribution. In Section 2, the basis of our method is discussed, along with some illustrative examples. The application of this method to a variety of molecules and clusters, with a special emphasis on estimating the HB energy along with the energy of other NCIs is presented in Section 3. Section 4 discusses some computational strategies for applying our method to large molecular clusters. The last Section provides a summary and a discussion on future developments.

Quantification of hydrogen bond energy based on equations using spectroscopic, structural, QTAIM-based, and NBO-based descriptors which calibrated by the molecular tailoring approach

CONTEXT

Hydrogen bonds critically influence the structure and properties of both organic molecules and biomolecules, as well as supramolecular assemblies. For this reason, the development and elaboration of methods for quantitative assessment of hydrogen bond energy is an urgent challenge. In this study, using a large series of hydroxycarbonyl aliphatic compounds with the O‒H···O = C intramolecular hydrogen bond, a bank of hydrogen bond descriptors was created, including spectroscopic, structural, QTAIM-based, and NBO-based parameters. It was shown that the O‒H vibration frequency, OH chemical shift as the spectroscopic descriptors, the O···H hydrogen bond length, O···O distance, and O‒H covalent bond length as the structural descriptors, the electron density and its Laplacian, electron potential energy density in the hydrogen bond critical point, the electron density at the ring critical point as the QTAIM-based descriptors change in a correlated manner. The same correlation is found in change of the charge transfer energy through a hydrogen bond, the occupancy of the O‒H bond antibonding orbital, the Wiberg indices of the O···H hydrogen bond, and the O‒H covalent bond, as well as the polarization of the O‒H bond, which are the NBO-based descriptors. It was also recognized that the specified descriptors from the spectroscopic, structural, QTAIM-based, and NBO-based categories are functionally related to the values of intramolecular hydrogen bond energy, quantified via the molecular tailoring approach. This allowed one to obtain a system of equations for quantitative estimation of intramolecular hydrogen bond energy based on the spectroscopic, structural, QTAIM, and NBO descriptors, which makes such quantification more dependable and reliable.

METHODS

To obtain the spectroscopic descriptors, the vibrational spectra and shielding constants were calculated using the GIAO method. Structural descriptors were obtained for the equilibrium geometry of molecules, calculated at the MP2(FC)/6-311 +  + (2d,2p) level using the Gaussian 09 program. The QTAIM-based descriptors were calculated using the AIMAll program within the framework of the quantum theory "Atoms in Molecules." The NBO-based descriptors were calculated using the NBO 3.1 program implemented into Gaussian 09. To quantify the energy of intramolecular hydrogen bonds, molecular fragmentation was used within the molecular tailoring approach.

Revealing the Reasons for Degeneration of Resonance-Assisted Hydrogen Bond on the Aromatic Platform: Calculations of Ortho-, Meta-, Para-Disubstituted Benzenes, and (Z)-(E)-Olefins

The energies of the O-H∙∙∙O=C intramolecular hydrogen bonds were compared quantitatively for the series of ortho-disubstituted benzenes and Z-isomers of olefins via a molecular tailoring approach. It was established that the hydrogen bond energy in the former series is significantly less than that in the latter one. The reason for lowering the hydrogen bond energy in the ortho-disubstituted benzenes compared to the Z-isomers of olefins is the decrease in the π-contribution to the total energy of the complex interaction, in which the hydrogen bond per se is enhanced by the resonance effect. By the example of the para- and meta-disubstituted benzenes, as well as E-isomers of olefins, it was explicitly shown that the aromatic ring is a much poorer conductor of the resonance effect compared to the double bond. The hydrogen bond in the ortho-disubstituted benzenes has a lower energy than a typical resonance-assisted hydrogen bond because the aromatic moiety cannot properly assist the hydrogen bond with a resonance effect. Thus, a hydrogen bond on an aromatic platform should fall into a special category, namely an aromaticity-assisted hydrogen bond, which is closer by nature to a simple hydrogen bond rather than to a resonance-assisted one.

Guide to tuning the chalcone molecular properties based on the push-pull effect energy scale created via the molecular tailoring approach

Using the molecular tailoring approach, a total energy scale for the push-pull effect in the range from -40 to 100 kcal/mol is created for the wide series of neutral, charged and doubly charged compounds on the chalcone platform. Taking into account similar energy scale for hydrogen bonds, the strength of the push-pull effect is ranked in the seven categories, ranging from negative (anti-push-pull) to very weak and very strong push-pull effect. It is demonstrated that the molecular properties of chalcone can be tuned prior synthesis due to the created energy scale for the push-pull effect. The single bonds of the π-spacer in the chalcones are shortened, the double ones are lengthened, and the C=O bond vibrations are red shifted when the push-pull effect is enhanced along the energy scale. The HOMO and LUMO energies change systematically while the HOMO-LUMO energy gap narrows as the strength of the push-pull effect increases.

In vitro cell composition identification of wood decay fungi by Fourier transform infrared spectroscopy

Decay pathogens follow dissimilar metabolic mechanisms to cause irreversible damage to woody tissues. The objective of this study is to perform inter- and intra-species microbial cell structural comparison using attenuated total reflectance Fourier transform infrared spectroscopy. Representative fungi species, causing brown rot and white rot, namely, Postia placenta and Trametes versicolor, respectively, were cultured in laboratory conditions. In vitro spectral measurements were performed at periodic two week intervals of fungal growth. The study shows structural differences for both species of fungi. The prominent presence of protein amide, carbohydrate and carboxyl bands was of interest. Spectral deconvolution of the infrared broadband around approximately 3300 cm-1 produced peaks at four different wavenumbers. The hydrogen bond energy obtained at the four wavenumbers, from deconvolution, varied from approximately 41 kJ mol-1 to approximately 7 kJ mol-1, indicating the presence of strong and weak forces in microbial cell structure. The hydrogen bond distance, obtained at the deconvoluted wavenumbers, varied between 2.7 Å-2.8 Å, indicating the presence of short and long-distance forces within microbial cells. Microscopic observation showed mycelium colonization, hyphal tip and lateral branching.

Interaction-Deletion: A Composite Energy Method for the Optimization of Molecular Systems Selectively Removing Specific Nonbonded Interactions

The complex interactions between different portions of a large molecule can be challenging to analyze through traditional electronic structure calculations. Moreover, standard methods cannot easily quantify the physical consequences of individual pairwise interactions inside a molecule. By creating a set of molecular fragments, we propose a composite energy method to explore changes in a molecule caused by removing selected nonbonded interactions between different molecular portions. Energies and forces are easily obtained with this composite approach, allowing geometry optimizations that lead to chemically meaningful structures that describe how the omitted interactions contribute to changes in the local geometrical minima. We illustrate the application of our new hybrid scheme by computing the influence of intramolecular hydrogen-bonding interactions in two small molecules: 1,6-(tG⁺G⁺TG⁺G⁺g^-)-hexanediol and a cyclic analogue, cis-1,4-cyclohexanediol. The resulting structural and energetic changes are interpreted to yield key physical insights and quantify concepts such as "preparation energy" or "reorganization energy". We demonstrate that the composite method can be extended to larger molecular systems by showing its application on a Si(100) surface model containing interactions between dissociated ammonia molecules on adjacent surface dimers. The scheme's efficacy is also tested by applying it to systems having multiple intramolecular interactions, viz., 3₁₀-polyglycine and H⁺GPGG. Furthermore, the cooperative nature of intramolecular hydrogen bonds is explored by using interaction-deletion in 2-nitrobenzene-1,3-diol.

Molecular Tailoring Approach for the Estimation of Intramolecular Hydrogen Bond Energy

Hydrogen bonds (HBs) play a crucial role in many physicochemical and biological processes. Theoretical methods can reliably estimate the intermolecular HB energies. However, the methods for the quantification of intramolecular HB (IHB) energy available in the literature are mostly empirical or indirect and limited only to evaluating the energy of a single HB. During the past decade, the authors have developed a direct procedure for the IHB energy estimation based on the molecular tailoring approach (MTA), a fragmentation method. This MTA-based method can yield a reliable estimate of individual IHB energy in a system containing multiple H-bonds. After explaining and illustrating the methodology of MTA, we present its use for the IHB energy estimation in molecules and clusters. We also discuss the use of this method by other researchers as a standard, state-of-the-art method for estimating IHB energy as well as those of other noncovalent interactions.

O-H stretching frequency red shifts do not correlate with the dissociation energies in the dimethylether and dimethylsulfide complexes of phenol derivatives

In this perspective, we present a comprehensive report on the spectroscopic and computational investigations of the hydrogen bonded (H-bonded) complexes of Me2O and Me2S with seven para-substituted H-bond donor phenols. The salient finding was that although the dissociation energies, D0, of the Me2O complexes were consistently higher than those of the analogous Me2S complexes, the red-shifts in phenolic O-H frequencies, Δν(O-H), showed the exactly opposite trend. This is in contravention of the general perception that the red shift in the X-H stretching frequency in the X-HY hydrogen bonded complexes is a reliable indicator of H-bond strength (D0), a concept popularly known as the Badger-Bauer rule. This is also in contrast to the trend reported for the H-bonded complexes of H2S/H2O with several para substituted phenols of different pKa values wherein the oxygen centered hydrogen bonded (OCHB) complexes consistently showed higher Δν(O-H) and D0 compared to those of the analogous sulfur centered hydrogen bonded (SCHB) complexes. Our effort was to understand these intriguing observations based on the spectroscopic investigations of 1 : 1 complexes in combination with a variety of high level quantum chemical calculations. Ab initio calculations at the MP2 level and the DFT calculations using various dispersion corrected density functionals (including DFT-D3) were performed on counterpoise corrected surfaces to compute the dissociation energy, D0, of the H-bonded complexes. The importance of anharmonic frequency computations is underscored as they were able to correctly reproduce the observed trend in the relative OH frequency shifts unlike the harmonic frequency computations. We have attempted to find a unified correlation that would globally fit the observed red shifts in the O-H frequency with the H-bonding strength for the four bases, namely, H2S, H2O, Me2O, and Me2S, in this set of H-bond donors. It was found that the proton affinity normalized Δν(O-H) values scale very well with the H-bond strength.

Finding the direct energy-structure correlations in intramolecular aromaticity assisted hydrogen bonding (AAHB)

A predictive model for intramolecular hydrogen bond energy (E_HB) calculation of polyaromatic ortho-hydroxyaldehydes based on a set of small, functionalized hydrocarbons is developed. The complete data set of 18 compounds was used for this study. The model is based on one of four optional categories of molecular descriptors: geometric, spectroscopic, bond order and topological indices. The model of Wiberg bond indices (WBIs) as descriptors of the CC involved bond based on stepwise regression has acceptable prediction abilities for 14 structures of ortho-hydroxyformylobenzo[a]pyrene derivatives already at the semiempirical level. The presented correlation enables a significantly more rapid and quantitative description of the hydrogen bonding strength than the much more time-consuming MTA method. Thus, WBIs are shown to provide a reliable means for fast prescreening of the energy of chelate hydrogen bonds potentially for any polyaromatic derivatives.

Intramolecular Hydrogen Bond Energy and Its Decomposition—O–H∙∙∙O Interactions

The method to calculate the energy of intramolecular hydrogen bond is proposed and tested for a sample of malonaldehyde and its fluorine derivatives; the corresponding calculations were performed at the ωB97XD/aug-cc-pVTZ level. This method based on relationships found for related intermolecular hydrogen bonds is compared with other approaches which may be applied to estimate the intramolecular hydrogen bond energy. Particularly, methods based on the comparison of the system that contains the intramolecular hydrogen bond compared with corresponding conformations where such interaction does not occur are discussed. The function-based energy decomposition analysis, FB-EDA, of the intramolecular hydrogen bonds is also proposed here.

Approximating the Strength of the Intramolecular Hydrogen Bond in 2-Fluorophenol and Related Compounds: A New Application of a Classic Technique

Fluorinated organic compounds are ubiquitous in the pharmaceutical and agricultural industries. To better discern the mode of action of these compounds, it is critical to understand the strengths of hydrogen bonds involving fluorine. While established techniques can determine these strengths for intermolecular complexes, there is no analogous scheme for intramolecular hydrogen bonds. This work uses ¹H nuclear magnetic resonance spectroscopy to measure the strength of intramolecular hydrogen bonds in ortho-substituted phenols. Titration of each phenol with DMSO in CCl₄ yields a free energy of binding (ΔG). Subtraction of this value from the ΔG of binding of the standard, 4-fluorophenol, is shown to give the difference in ΔG for the cis and trans isomers of the ortho-substituted phenols. This difference is conventionally taken to be approximately equal to the ΔG of the intramolecular hydrogen bond. These data complement theoretical methods, which yield slightly larger ΔGs. Both theory and experiment point to a weak intramolecular hydrogen bond in 2-fluorophenol. The other 2-X-phenols have stronger hydrogen bonds, following the order F < Cl ≈ Br < OCH₃. The methodology developed here can be readily applied to other systems with intramolecular hydrogen bonds.

Quantitative decomposition of resonance-assisted hydrogen bond energy in β-diketones into resonance and hydrogen bonding (π- and σ-) components using molecular tailoring and function-based approaches

Using the molecular tailoring and function-based approaches allows one to divide the energy of the O─H⋯O═C resonance-assisted hydrogen bond in a series of the β-diketones into resonance and hydrogen bonding components. The magnitude of the resonance component is assessed as about 6 kcal mol^-1 . This value increases by ca. 1 kcal mol^-1 on going from the weak to strong resonance-assisted hydrogen bonding. The magnitude of the hydrogen bonding component varies in the wide range from 2 to 20 kcal mol^-1 depending on the structure of the β-diketone in question.

A molecular tailoring approach – a new guide to quantify the energy of push–pull effects: a case study on (E)-3-(1H-pyrrol-2-yl)prop-2-enones

We report a new approach to quantify the push–pull effect in molecules with internal charge transfer.

Fragmentation method reveals a wide spectrum of intramolecular hydrogen bond energies in antioxidant natural products

Very strong and weak IHBs in curcumin.

Cis/Trans Energetics in Epoxide, Thiirane, Aziridine and Phosphirane Containing Cyclopentanols: Effects of Intramolecular OH⋯O, S, N and P Contacts

A recent computational analysis of the stabilizing intramolecular OH⋯O contact in 1,2-dialkyl-2,3-epoxycyclopentanol diastereomers has been extended to thiiriane, aziridine and phosphirane analogues. Density functional theory (DFT), second-order Møller-Plesset perturbation theory (MP2) and CCSD(T) coupled-cluster computations with simple methyl and ethyl substituents indicate that electronic energies of the cis isomers are lowered by roughly 3 to 4 kcal mol⁻¹ when the OH group of these cyclopentanol systems forms an intramolecular contact with the O, S, N or P atom on the adjacent carbon. The results also suggest that S and P can participate in these stabilizing intramolecular interactions as effectively as O and N in constrained molecular environments. The stabilizing intramolecular OH⋯O, OH⋯S, OH⋯N and OH⋯P contacts also increase the covalent OH bond length and significantly decrease the OH stretching vibrational frequency in every system with shifts typically on the order of −41 cm⁻¹.

Estimating the energy of intramolecular hydrogen bonds from 1H NMR and QTAIM calculations

The values of the downfield chemical shift of the bridge hydrogen atom were estimated for a series of compounds containing an intramolecular hydrogen bond O-HO, O-HN, O-HHal, N-HO, N-HN, C-HO, C-HN and C-HHal. Based on these values, the empirical estimation of the hydrogen bond energy was obtained by using known relationships. For the compounds containing an intramolecular hydrogen bond, the DFT B3LYP/6-311++G(d,p) method was used both for geometry optimization and for QTAIM calculations of the topological parameters (electron density ρ^BCP and the density of potential energy V in the critical point of the hydrogen bond). The calculated geometric and topological parameters of hydrogen bonds were also used to evaluate the energy of the hydrogen bond based on the equations from the literature. Comparison of calibrating energies from the ¹H NMR data with the energies predicted by calculations showed that the most reliable are the linear dependence on the topological ρ^BCP and V parameters. However, the correct prediction of the hydrogen bond energy is determined by proper fitting of the linear regression coefficients. To obtain them, new linear relationships were found between the calculated ρ^BCP and V parameters and the hydrogen bond energies obtained from empirical ¹H NMR data. These relationships allow the comparison of the energies of different types of hydrogen bonds for various molecules and biological ensembles.

Energy of Intramolecular Hydrogen Bonding in ortho-Hydroxybenzaldehydes, Phenones and Quinones. Transfer of Aromaticity from ipso-Benzene Ring to the Enol System(s)

Intramolecular hydrogen bonding (HB) is one of the most studied noncovalent interactions of molecules. Many physical, spectral, and topological properties of compounds are under the influence of HB, and there are many parameters used to notice and to describe these changes. Hitherto, no general method of measurement of the energy of intramolecular hydrogen bond (EHB) has been put into effect. We propose the molecular tailoring approach (MTA) for EHB calculation, modified to apply it to Ar-O-H∙∙∙O=C systems. The method, based on quantum calculations, was checked earlier for hydroxycarbonyl-saturated compounds, and for structures with resonance-assisted hydrogen bonding (RAHB). For phenolic compounds, the accuracy, repeatability, and applicability of the method is now confirmed for nearly 140 structures. For each structure its aromaticity HOMA indices were calculated for the central (ipso) ring and for the quasiaromatic rings given by intramolecular HB. The comparison of calculated HB energies and values of estimated aromaticity indices allowed us to observe, in some substituted phenols and quinones, the phenomenon of transfer of aromaticity from the ipso-ring to the H-bonded ring via the effect of electron delocalization.

Assessment of the Presence and Strength of H-Bonds by Means of Corrected NMR

The downfield shift of the NMR signal of the bridging proton in a H-bond (HB) is composed of two elements. The formation of the HB causes charge transfer and polarization that lead to a deshielding. A second factor is the mere presence of the proton-accepting group, whose electron density and response to an external magnetic field induce effects at the position of the bridging proton, exclusive of any H-bonding phenomenon. This second positional shielding must be subtracted from the full observed shift in order to assess the deshielding of the proton caused purely by HB formation. This concept is applied to a number of H-bonded systems, both intramolecular and intermolecular. When the positional shielding is removed, the remaining chemical shift is in much better coincidence with other measures of HB strength.

Analysis of Different Fragmentation Strategies on a Variety of Large Peptides: Implementation of a Low Level of Theory in Fragment-Based Methods Can Be a Crucial Factor

We have investigated the performance of two classes of fragmentation methods developed in our group (Molecules-in-Molecules (MIM) and Many-Overlapping-Body (MOB) expansion), to reproduce the unfragmented MP2 energies on a test set composed of 10 small to large biomolecules. They have also been assessed to recover the relative energies of different motifs of the acetyl(ala)18NH2 system. Performance of different bond-cutting environments and the use of Hartree-Fock and different density functionals (as a low level of theory) in conjunction with the fragmentation strategies have been analyzed. Our investigation shows that while a low level of theory (for recovering long-range interactions) may not be necessary for small peptides, it provides a very effective strategy to accurately reproduce the total and relative energies of larger peptides such as the different motifs of the acetyl(ala)18NH2 system. Employing M06-2X as the low level of theory, the calculated mean total energy deviation (maximum deviation) in the total MP2 energies for the 10 molecules in the test set at MIM(d=3.5Å), MIM(η=9), and MOB(d=5Å) are 1.16 (2.31), 0.72 (1.87), and 0.43 (2.02) kcal/mol, respectively. The excellent performance suggests that such fragment-based methods should be of general use for the computation of accurate energies of large biomolecular systems.

Resonance-Assisted Hydrogen Bonding as a Driving Force in Synthesis and a Synthon in the Design of Materials

Resonance-assisted hydrogen bonding (RAHB), a concept introduced by Gilli and co-workers in 1989, concerns a kind of intramolecular H-bonding strengthened by a conjugated π-system, usually in 6-, 8-, or 10-membered rings. This Review highlights the involvement of RAHB as a driving force in the synthesis of organic, coordination, and organometallic compounds, as a handy tool in the activation of covalent bonds, and in starting moieties for synthetic transformations. The unique roles of RAHB in molecular recognition and switches, E/Z isomeric resolution, racemization and epimerization of amino acids and chiral amino alcohols, solvatochromism, liquid-crystalline compounds, and in synthons for crystal engineering and polymer materials are also discussed. The Review can provide practical guidance for synthetic chemists that are interested in exploring and further developing RAHB-assisted synthesis and design of materials.

Topological insights into the 1/1 diacetyl/water complex gained using a new methodological approach

The 1/1 diacetyl/water complex is of atmospheric relevance. Previous experimental and theoretical studies have focused on two isomeric forms, and geometry optimizations were carried out on them. Herein, we propose a six-step methodological approach based on topological properties to search for and characterize all of the isomeric forms of the 1/1 noncovalent diacetyl/water complex: (1) a molecular electrostatic potential (MESP) study to get an overview of the V min and V max regions on the molecular surfaces of the separate molecules (diacetyl and water); (2) a topological (QTAIM and ELF) study allowing thorough characterization of the electron densities (QTAIM) and irreducible ELF basins of the separate molecules; (3) full optimization of the predicted structures based on the interaction between complementary reaction sites; (4) energetic characterization based on a symmetry-adapted perturbation theory (SAPT) analysis; (5) topological characterization of the optimized complexes; (6) analysis of the complexes in terms of orbital overlaps (natural bond orbitals, NBO analysis). Using this approach, in addition to achieving the topological characterization of the two isomeric forms already reported, a third possible isomer was identified and characterized. Graphical Abstract Topological tools to study monohydrated complexes.

Intramolecular O-H···O═C hydrogen bond energy via the molecular tailoring approach to RAHB structures

A method for the calculation of the intramolecular hydrogen bond (HB) energy (EHB) by molecular tailoring approach for hydroxycarbonyl aliphatic compounds has been used for compounds with resonance-assisted hydrogen bonding (RAHB). The intramolecular hydrogen bond energies estimated for 229 structures (of 186 compounds) range from 8.2 to 26.3 kcal/mol and show correlation with the geometry descriptors of hydrogen bonds, with the calculated frequencies as well as with topological parameters obtained from the atoms in molecules (AIM) theory. These correlations differ significantly from obtained formerly for saturated nonenolizable structures and prove the special character of the resonance-assisted hydrogen-bonded systems.