Intramolecular hydrogen bond energy and cooperative interactions in α-, β-, and γ-cyclodextrin conformers

Accurate estimation of individual intramolecular hydrogen bond (H-bond) energies is an intricate task for multiply H-bonded systems. In such cases, the hydrogen bond strengths could be highly influenced by the cooperative interactions, for example, those between hydroxyl groups in sugars. In this work, we use the recently proposed molecular tailoring approach-based quantification (Deshmukh, Gadre, and Bartolotti, J Phys Chem A 2006, 110, 12519) to the extended systems of cyclodextrins (CDs). Further, the structure and stability of different conformers of α-, β-, and γ-CDs are explained based on the energetics and cooperative contribution to the strength of these H-bonds. The estimated O-H···O H-bond energies in the various CD conformers are found to vary widely from 1.1 to 8.3 kcal mol(-1). The calculated energy contributions to cooperativity toward the H-bond strengths fall in the range of 0.25-2.75 kcal mol(-1).     

Two-way effects between hydrogen bond and intramolecular resonance effect: an ab initio study on complexes of formamide and its derivatives with water

Ab initio calculations up to MP2/aug-cc-pVTZ//MP2/cc-pVTZ level, including natural charge population and natural resonance theory analyses, have been carried out to study the two-way effects between hydrogen bond (H-bond) and the intramolecular resonance effect by using the H-bonded complexes of formamide ( FAO) and its derivatives ( FAXs, X represents the heavy atoms in the substituent groups, CH 2, NH, SiH 2, PH, and S) with water as models. Unlike NH 3 and NH 2CH 3 which prefer being H-bond acceptors ( HA) to form H-bond with water, the amino groups in the six monomers, because of the resonance effect, prefer being H-bond donors ( HD) rather HA. Six monomers can all form HD complexes with water, and only two ( FAC and FASi) with the weakest resonance effect are able to form HA complexes with water. The HD H-bond and resonance effect enhance each other (positive two-way effects) whereas the HA H-bond and resonance effect weaken each other (negative two-way effects). The H-bond energies in the six HD complexes are nearly linearly correlated with the weights of the dipolar resonance in Pauling's model and the N-C bond lengths; the correlation coefficients are 0.91 and 0.93, respectively. The positive two-way effects also happens in FAO-water complex, in which the FAO CO group serves as HA ( HA co ). Interestingly, when the HD and HA co H-bonds are present in FAO H-bond complex simultaneously, the enhancements are much more significant, and the energies of the two types of H-bonds are much larger than those when only one type of H-bond is present, reflecting the cooperative effects. By using the knowledge to the two-way effects, we computationally designed a molecule ( FAO- BH 3 ) to increase H-bond energy. Because of the oxygen lone pair donation to the empty pi orbital of BH 3, FAO- BH 3 has a much stronger resonance effect than FAO. As a result, the H-bond energy (-5.55 kcal/mol) in HD H 2O ... FAO- BH 3 complex is much greater than the -3.30 kcal/mol in the HD H 2O...FAO complex. The two-way effects can be rationalized as follows: the resonance effect leads to intramolecular charge shifts in the monomers which facilitate or prevent the charge donation or acceptation of their H-bond partners. Therefore, the H-bonds are strengthened or weakened. In reverse, the charge donations or acceptations of their H-bond partners facilitate or prevent the intramolecular charge shifts in the monomer moieties, which enhance or weaken the resonance effect. The understanding to the two-way effects may be helpful in drug design and refinement by modulating the H-bond strength and in building empirical H-bond models to study large biological molecules. The study supports Pauling's resonance model.     

Intramolecular hydrogen bond energy in polyhydroxy systems: a critical comparison of molecular tailoring and isodesmic approaches

The intramolecular hydrogen bond (H-bond) energies in several polyhydroxy systems are estimated using an isodesmic/homodesmic reaction approach as well as a molecular tailoring approach (MTA) [Deshmukh, M. M.; Gadre, S. R.; Bartolotti, L. J. J. Phys. Chem. A 2006, 110, 12519]. It is shown that the isodesmic/homodesmic reaction approach as advocated in the literature does not give true H-bond energy but includes the effect of strain energy due to the formation of a ring structure. Such a ring strain is duly accounted for in the MTA method. The isodesmic H-bond energies are found to be smaller than their MTA energy counterparts typically by the strain energy. The MTA is applied to decitol, a system with more than five different H-bonds for which an application of an appropriate isodesmic reaction is extremely difficult. It has been shown that the MTA method is able to predict not only the H-bond energies but also the trends in conformational energies for three different conformers of decitol studied in the present work.     

Nature of one-dimensional short hydrogen bonding: bond distances, bond energies, and solvent effects

On the basis of recently synthesized calix[4]hydroquinone (CHQ) nanotubes which were self-assembled with infinitely long one-dimensional (1-D) short hydrogen bonds (SHB), we have investigated the nature of 1-D SHB using first-principles calculations for all the systems including the solvent water. The H-bonds relay (i.e., contiguous H-bonds) effect in CHQs shortens the H...O bond distances significantly (by more than 0.2 A) and increases the bond dissociation energy to a large extent (by more than approximately 4 kcal/mol) due to the highly enhanced polarization effect along the H-bond relay chain. The H-bonds relay effect shows a large increase in the chemical shift associated with the SHB. The average binding energies for the infinite 1-D H-bond arrays of dioles and dions increase by approximately 4 and approximately 9 kcal/mol per H-bond, respectively. The solvent effect (due to nonbridging water molecules) has been studied by explicitly adding water molecules in the CHQ tube crystals. This effect is found to be small with slight weakening of the SHB strength; the H...O bond distance increases only by 0.02 A, and the average binding energy decreases by approximately 1 kcal/mol per H-bond. All these results based on the first-principles calculations are the first detailed analysis of energy gain by SHB and energy loss by solvent effect, based on a partitioning scheme of the interaction energy components. These reliable results elucidate not only the self-assembly phenomena based on the H-bond relay but also the solvent effect on the SHB strength.     

Estimation of intramolecular hydrogen bond energy via molecular tailoring approach

A novel method, based on the molecular tailoring approach for estimating intramolecular hydrogen bond energies, is proposed. Here, as a case study, the O-H...O bond energy is directly estimated by addition/subtraction of the single point individual fragment energies. This method is tested on polyhydroxy molecules at MP2 and B3LYP levels of theory. It is seen to be able to distinguish between weak ( approximately 1 kcal mol(-1)) and moderately strong ( approximately 5 kcal mol(-1)) hydrogen bonds in polyhydroxy molecules.     

Charge-assisted intramolecular hydrogen bonds in disubstituted cyclohexane derivatives

In this paper, the N(+)-H···N, N(+)-H···O, and O-H···O(-) charge-assisted intramolecular hydrogen bonds (CAHBs) are investigated using different theoretical approaches. Monocharged cyclohexyldiamines (CHDA), aminocyclohexanols (ACHO), and cyclohexanediols (CHDO) are used as model compounds. Geometry optimizations at the MP2/aug-cc-pVDZ level are used to find the equilibrium structures for all possible H-bonded conformers. CAHBs are characterized geometrically and spectroscopically, and their energy is evaluated by means of homodesmic reactions. By comparison with the neutral forms, the presence of the charge is found to have a deep influence on the geometric and energetic H-bond parameters. In addition, these parameters are strongly dependent on the type of the groups involved as well as on their relative position in the cyclohexyl ring. For the systems under study, the H-bond energies vary from -23 to -113 kJ mol(-1), being classified from moderate to strong H-bonds. These H-bonds are also characterized by the application of the NBO and AIM theories. NBO analysis reveals that the energy corresponding to the charge transfer between the lone-pairs of the electron donor group and the antibonding orbitals of the acceptor group represents an important contribution in the H-bond stabilization. From the application of the AIM theory it is possible to see that these H-bonds possess some covalence which varies according to the type and relative position of the intervenient groups.     

Unique charge-separated intermolecular and eight-membered intramolecular H-bonds in bis-(thio)barbiturates

Reaction of symmetrical and unsymmetrical (thio)barbituric acids with aldehydes in the presence of triethylamine afforded a new form of bis-(thio)barbiturate containing charge-separated inter- and eight-membered intramolecular H-bonds. The reaction products were obtained as bis-(thio)barbiturates containing eight-membered intramolecular H-bond in the presence of l-(+)-tartaric acid (TA). The intramolecular H-bond strength (kcal/mol) and corresponding pKa value for 4ab′ were estimated to be 37 kcal/mol and ?1.3, respectively.     

Cooperativity between OH...O and CH...O hydrogen bonds involving dimethyl sulfoxide-H2O-H2O complex

The cooperativity between the O-H...O and C-H...O hydrogen bonds has been studied by quantum chemical calculations at the MP2/6-311++G(d,p) level in gaseous phase and at the B3LYP/6-311++G(d,p) level in solution. The interaction energies of the O-H...O and C-H...O H-bonds are increased by 53 and 58%, respectively, demonstrating that there is a large cooperativity. Analysis of hydrogen-bonding lengths, OH bond lengths, and OH stretching frequencies also supports such a conclusion. By NBO analysis, it is found that orbital interaction plays a great role in enhancing their cooperativity. The strength increase of the C-H...O H-bond is larger than that of the O-H...O H-bond due to the cooperativity. The solvent has a weakening effect on the cooperativity.     

Estimation on the intramolecular 10-membered ring N-H...O=C hydrogen-bonding energies in glycine and alanine peptides

Computation of accurate intramolecular hydrogen-bonding energies for peptides is of great importance in understanding the conformational stabilities of peptides and developing a more accurate force field for proteins. We have proposed a method to determine the intramolecular seven-membered ring N-H...O=C hydrogen-bonding energies in glycine and alanine peptides. In this article, the method is further applied to evaluate the intramolecular 10-membered ring N-H...O=C hydrogen-bonding energies in peptides. The optimal structures of the intramolecular 10-membered ring N-H...O=C hydrogen bonds in glycine and alanine tripetide molecules are obtained at the MP2 level with 6-31G(d), 6-31G(d,p), and 6-31+G(d,p) basis sets. The intramolecular 10-membered ring N-H...O=C hydrogen-bonding energies are then evaluated based on our method at the MP2/6-311++G(3df,2p) level with basis set superposition error correction. The intramolecular 10-membered ring N-H...O=C hydrogen-bonding energies are calculated to be in the range of -6.84 to -7.66, -4.44 to -4.98, and -6.95 to -7.88 kcal/mol. The method is also applied to estimate the individual intermolecular hydrogen-bonding energies in the dimers of amino-acetaldehyde, 2-amino-acetamide, formamide, and oxalamide, each dimer having two identical intermolecular hydrogen bonds. According to our method, the individual intermolecular hydrogen-bonding energies in the four dimers are calculated to be -1.77, -1.67, -6.35, and -4.82 kcal/mol at the MP2/6-311++G(d,p) level, which are in good agreement with the values of -1.84, -1.72, -6.23, and -4.93 kcal/mol predicted by the supermolecular method.     

Single-site mutation and secondary structure stability: an isodesmic reaction approach. The case of unnatural amino acid mutagenesis Ala-->Lac

A method is described to evaluate backbone interactions in proteins via computational unnatural amino acid mutagenesis. Several N-acetyl polyalanyl amides (AcA(n)NH(2)) were optimized in the representative helical (3(10)-, 4(13)-, and a "hybrid" kappa-helix, n = 7, 9, 10, 14) and hairpin (two- and three-stranded antiparallel beta-sheets with type I turns betaalphaalphaepsilon, n = 6, 9, 10) conformations, and extended conformers of N-acetyl polyalanyl methylamides (n = 2, 3) were used to derive multistranded beta-sheet fragments. Subsequently, each residue of every model structure was substituted, one at a time, with l-lactic acid. The resulting mutant structures were again optimized, and group-transfer energies DeltaE(GT) were obtained as heats of the isodesmic reactions: AcA(n)NHR + AcOMe --> AcA(x)LacA(y)NHR + AcNHMe (R = H, CH(3)). These group-transfer energies correlate with the degree of charge polarization of the substituted peptide linkages as measured by the difference Deltae in H and O Mulliken populations in HN-C=O and with the H-bond distances in the "wild-type" structures. A good correlation obtains for the HF/3-21G and B3LYP/6-31G* group-transfer energies. The destabilization effects are interpreted in terms of loss of interstrand and intrastrand H-bonds, decrease in Lewis basicity of the C=O group, and O...O repulsion. On the basis of several comparisons of Ala --> Lac DeltaE(GT)'s with heats of the NH --> CH(2) substitutions, the latter contribution is estimated (B3LYP/6-31G*) to range between 1.5 and 2.4 kcal mol(-1), a figure close to the recent experimental DeltaDeltaG(o) value of 2.6 kcal mol(-1) (McComas, C. C.; Crowley, B. M.; Boger, D. L. J. Am.Chem. Soc. 2003, 125, 9314). The partitioning yields the following maximum values of the electronic association energy of H-bonds in the examined sample of model structures (B3LYP/6-31G* estimates): 3(10)-helix D(e) = -1.7 kcal mol(-1), alpha-helix D(e) = -3.8 kcal mol(-1), beta-sheet D(e) = -6.1 kcal mol(-1). The premise of experimental evaluations of the backbone-backbone H-bonding that Ala --> Lac substitution in proteins is isosteric (e.g., Koh, J. T.; Cornish, V. W.; Schultz, P. G. Biochemistry 1997, 36, 11314) is often but not always corroborated. Examination of the integrity of H-bonding pattern and phi(i), psi(i) distribution identified several mutants with significant distortions of the "wild-type" structure resulting inter alia from the transitions between i, i + 3 and i, i + 4 H-bonding in helices, observed previously in the crystallographic studies of depsipeptides (Ohyama, T.; Oku, H.; Hiroki, A.; Maekawa, Y.; Yoshida, M.; Katakai, R. Biopolymers 2000, 54, 375; Karle, I. L.; Das, C.; Balaram, P. Biopolymers 2001, 59, 276). Thus, the isodesmic reaction approach provides a simple way to gauge how conformation of the polypeptide chain and dimensions of the H-bonding network affect the strength of backbone-backbone C=O...HN bonds. The results indicate that the stabilization provided by such interactions increases on going from 3(10)-helix to alpha-helix to beta-sheet.     

Vibrational spectra, co-operative intramolecular hydrogen bonding and conformations of calix[4]arene and thiacalix[4]arene molecules and their para-tert-butyl derivatives

The IR and Raman spectra and conformations of calix[4]arene, thiacalix[4]arene and their p-tert-butyl derivatives have been analysed within the framework of scaled quantum mechanics (SQM). It is shown that the introduction of four p-tert-Bu groups into the calixarene molecules influences the relative energies of their conformers and the enthalpy of the cooperative intramolecular H-bonding (DeltaH(intra)) almost negligibly. DeltaH(intra), evaluated from Iogansen's rule, amounts to approximately 26-28 kcal mol(-1) for the calixarenes and approximately 20-21 kcal mol(-1) for the thiacalixarenes, which essentially exceeds the enthalpies of non-cooperative H-bonds formed by related phenols. As a result of this strong bonding, bands of stretching, bending and torsion vibrations of an eight-membered cyclic system (OH...)4 arise, e.g., two delta(OH)4 bands are observed in the IR spectra of the most highly symmetric C4 cone conformations of calix[4]arene and thiacalix[4]arene. The "duplication" of the number of OH infrared bands is a good new indicator of cooperativity of intramolecular H-bonding of the calixarenes.     

How strong can the bend be on a DNA helix from cisplatin? DFT and MP2 quantum chemical calculations of cisplatin-bridged DNA purine bases

The B3LYP/6-31G(d) level of theory was used for the optimization of [Pt(NH(3))(4)](2+), [Pt(NH(3))(3)(H(2)O)](2+), cis-[Pt(NH(3))(2)(H(2)O)(2)](2+), and related platinum complexes. In addition, water or ammonium ligands were replaced by DNA purine bases so that finally cis-diammineplatinum with two bases (Pt-bridged complexes) is obtained. Single point calculations using the MP2/6-31+G(d) method were performed on the obtained reference geometries and were utilized for estimating bond dissociation energies (BDEs) and stabilization energies, and for electron density analyses. After reoptimization, IR spectra were determined from HF second derivatives. It was found that replacement of both water and ammonium by the DNA base is an exothermic process (20-50 kcal/mol depending on the ligands present in the complex). Asymmetric structures with one interbase H-bond were obtained for cis-diammine[bond](N(7),N(7)'-diadenine)[bond]platinum and mixed cis-diammine[bond](N(7)-adenine)[bond](N(7)-guanine)[bond]platinum complexes. In the case of the diguanine Pt-bridge, a symmetrical complex with two ammonium...O(6) H-bonds was found. The higher stabilization energy of the di-guanine complex is linked to a larger component of the Coulombic interaction. However, the BDE of Pt[bond]N(7)(G) is smaller in this complex than the BDE of Pt[bond]N(7)(G) from the mixed Pt[bond]AG complex. Also, steric repulsion of the ligands is about 10 kcal/mol smaller for the asymmetrical Pt[bond]AA and Pt[bond]AG bridges. The influence of the trans effect on DBE can be clearly seen. Adenine exhibits the largest trans effect, followed by guanine, ammonium, and water. The strength of the H-bond can be determined from the IR spectra. The strongest H-bond is the interbase H-bridge between adenine and guanine in the mixed Pt[bond]AG complex; otherwise, the H-bonds of adenine complexes are weaker than in guanine complexes. BDE can be traced in the guanine-containing complexes. The nature of the covalent bonding is analyzed in terms of partial charges and MO. A general explanation of the lower affinity of transition metals to oxygen than nitrogen can be partially seen in the less favorable geometrical orientation of lone electron pairs of oxygen.     

Equilibrium thermodynamic studies in water: reactions of dihydrogen with rhodium(III) porphyrins relevant to Rh-Rh, Rh-H, and Rh-OH bond energetics

Aqueous solutions of rhodium(III) tetra p-sulfonatophenyl porphyrin ((TSPP)Rh(III)) complexes react with dihydrogen to produce equilibrium distributions between six rhodium species including rhodium hydride, rhodium(I), and rhodium(II) dimer complexes. Equilibrium thermodynamic studies (298 K) for this system establish the quantitative relationships that define the distribution of species in aqueous solution as a function of the dihydrogen and hydrogen ion concentrations through direct measurement of five equilibrium constants along with dissociation energies of D(2)O and dihydrogen in water. The hydride complex ([(TSPP)Rh-D(D(2)O)](-4)) is a weak acid (K(a)(298 K) = (8.0 +/- 0.5) x 10(-8)). Equilibrium constants and free energy changes for a series of reactions that could not be directly determined including homolysis reactions of the Rh(II)-Rh(II) dimer with water (D(2)O) and dihydrogen (D(2)) are derived from the directly measured equilibria. The rhodium hydride (Rh-D)(aq) and rhodium hydroxide (Rh-OD)(aq) bond dissociation free energies for [(TSPP)Rh-D(D(2)O)](-4) and [(TSPP)Rh-OD(D(2)O)](-4) in water are nearly equal (Rh-D = 60 +/- 3 kcal mol(-1), Rh-OD = 62 +/- 3 kcal mol(-1)). Free energy changes in aqueous media are reported for reactions that substitute hydroxide (OD(-)) (-11.9 +/- 0.1 kcal mol(-1)), hydride (D(-)) (-54.9 kcal mol(-1)), and (TSPP)Rh(I): (-7.3 +/- 0.1 kcal mol(-1)) for a water in [(TSPP)Rh(III)(D(2)O)(2)](-3) and for the rhodium hydride [(TSPP)Rh-D(D(2)O)](-4) to dissociate to produce a proton (9.7 +/- 0.1 kcal mol(-1)), a hydrogen atom (approximately 60 +/- 3 kcal mol(-1)), and a hydride (D(-)) (54.9 kcal mol(-1)) in water.     

Mechanism of OH formation from ozonolysis of isoprene: a quantum-chemical study

The formation and unimolecular reactions of primary ozonides and carbonyl oxides arising from the O(3)-initiated reactions of isoprene have been investigated using density functional theory and ab initio molecular orbital calculations. The activation energies of O(3) cycloaddition to the two double bonds of isoprene are found to be comparable (3.3-3.4 kcal mol(-1)), implying that the initial two O(3) addition pathways are nearly equally accessible. The reaction energies of O(3) addition to isoprene are between -47 and -48 kcal mol(-1). Cleavage of primary ozonides to form carbonyl oxides occurs with a barrier of 11-16 kcal mol(-1) above the ground state of the primary ozonide, and the decomposition energies range from -5 to -13 kcal mol(-1). OH formation is shown to occur primarily via decomposition of the carbonyl oxides with the syn-positioned methyl (alkyl) group, which is more favorable than isomerization to form dioxirane (by 1.1-3.3 kcal mol(-1)). Using the transition-state theory and master equation formalism, we determine an OH yield of 0.25 from prompt and thermal decomposition of the carbonyl oxides.     

Estimates of the energy of intramolecular hydrogen bonds

A method for the estimation of the energy of intramolecular hydrogen bonds in conjugated systems existing in a variety of conformations is presented. The method is applied to determine the intramolecular hydrogen bond energy in 3-aminopropenal and 3-aminopropenthial. According to the proposed estimation scheme, the intramolecular H-bond energies are found to be of the order of 5-7 kcal/mol. These results are compared with those obtained by using other estimation schemes as well as with the recent results by other authors. Also, the H-bond energies in dimers and trimers of the two molecules are calculated and compared with the corresponding data for internally hydrogen-bonded monomers. This comparison shows that the bond equalization effect is primarily due to proton donor-proton acceptor proximity. In comparison with intermolecular hydrogen bonds, the rigidity of the chelate skeleton enhances this proximity effect. The same effect can be seen in systems with intermolecular hydrogen bonds, although its magnitude is diminished because of the absence of additional forces which pull the proton donor and proton acceptor groups toward each other. No specific resonance-assisted origin of the intramolecular hydrogen bond energy seems to be needed to elucidate the energetics of these bonds.     

Hydrogen bonding of the nucleobase mimic 2-pyridone to fluorobenzenes: an ab initio investigation

The hydrogen-bonded complexes of the nucleobase mimic 2-pyridone (2PY) with seven different fluorinated benzenes (1-, 1,2-, 1,4-, 1,2,3-, 1,3,5-, 1,2,3,4-, and 1,2,4,5-fluorobenzene) are important model systems for investigating the relative importance of hydrogen bonding versus pi-stacking interactions in DNA. We have shown by supersonic-jet spectroscopy that these dimers are hydrogen bonded and not pi-stacked at low temperature (Leist, R.; Frey, J. A.; Leutwyler, S. J. Phys. Chem. A 2006, 110, 4180). Their geometries and binding energies D(e) were calculated using the resolution of identity (RI) M?ller-Plesset second-order perturbation theory method (RIMP2). The most stable dimers are bound by antiparallel N-H...F-C and C-H...O=C hydrogen bonds. The binding energies are extrapolated to the complete basis set (CBS) limit, , using the aug-cc-pVXZ basis set series. The CBS binding energies range from -D(e,CBS) = 6.4-6.9 kcal/mol and the respective dissociation energies from -D(0,CBS) = 5.9-6.3 kcal/mol. In combination with experiment, the latter represent upper limits to the dissociation energies of the pi-stacked isomers (which are not observed experimentally). The individual C-H...O=C and N-H...F-C contributions to D(e) can be approximately separated. They are nearly equal for 2PY.fluorobenzene; each additional F atom strengthens the C-H...O=C hydrogen bond by approximately 0.5 kcal/mol and weakens the C-F...H-N hydrogen bond by approximately 0.3 kcal/mol. The single H-bond strengths and lengths correlate with the gas-phase acid-base properties of the C-H and C-F groups of the fluorobenzenes.     

Quantum chemical evaluation of protein control over heme ligation: CO/O2 discrimination in myoglobin

Control of O2 versus CO binding in myoglobin (Mb) is tuned by a distal histidine residue through steric and H-bonding interactions. These interactions have been evaluated via Car-Parrinello DFT calculations, whose efficiency allows full quantum mechanical treatment of the 13 closest residues surrounding the heme. The small (8 degrees ) deviation of the Fe-C-O bond angle from linearity results from the steric influence of a distal valine residue and not the distal histidine. H-bond energies were evaluated by replacing the distal histidine with the non-H-bonding residue isoleucine. Binding energies for CO and O2 decreased by 0.8 and 4.1 kcal/mol for MbCO and MbO2, in good agreement with experimental H-bond estimates. Ligand discrimination is dominated by distal histidine H-bonding, which is also found to stabilize a metastable side-on isomer of MbO2 that may play a key role in MbO2 photodynamics.     

Cooperative 4-pyridone H-bonds with extraordinary stability. A DFT molecular orbital study

The N-H...O H-bonding enthalpy between 4-pyridones connected in a chain of H-bonds can achieve 23 kcal/mol for the most central H-bonds, while that between two 4-pyridones is 9.90 kcal/mol based upon DFT calculations on the counterpoise-corrected potential energy surfaces. That the range of enthalpies for N-H...O H-bonds can vary from as little as 2 to as much 23 kcal/mol depends primarily upon the polarizability of whatever internally connects the N-H and C=O within the H-bonding molecule, which are two parallel -C=C- entities in 4-pyridone. The contribution of covalent or charge-transfer interactions between the pi-systems of adjacent 4-pyridones is small.     

Dissection of complex molecular recognition interfaces

The synthesis of a family of zinc porphyrins and pyridine ligands equipped with peripheral H-bonding functionality has provided access to a wide range of closely related supramolecular complexes featuring between zero and four intramolecular H-bonds. An automated UV/vis titration system was used to characterize 120 different complexes, and these data were used to construct a large of number of different chemical double mutant cycles to quantify the intramolecular H-bonding interactions. The results probe the quantitative structure-activity relationship that governs cooperativity in the assembly of complex molecular recognition interfaces. Specifically, variations in the chemical structures of the complexes have allowed us to change the supramolecular architecture, conformational flexibility, geometric complementarity, the number and nature of the H-bond interactions, and the overall stability of the complex. The free energy contributions from individual H-bonds are additive, and there is remarkably little variation with architecture in the effective molarity for the formation of intramolecular interactions. Intramolecular H-bonds are not observed in complexes where they are geometrically impossible, but there are no cases where excellent geometric complementarity leads to very high affinities. Similarly, changes in conformational flexibility seem to have limited impact on the values of effective molarity (EM). The major variation that was found for all of the 48 intramolecular interactions that were examined using double mutant cycles is that the values of EM for intramolecular carboxylate ester-phenol H-bonds (200 mM) are an order of magnitude larger than those found for phosphonate diester-phenol H-bonds (30 mM). The corresponding intermolecular phosphonate diester-phenol H-bonds are 2 orders of magnitude more stable than carboxylate ester-phenol H-bonds, and the large differences in EM may be due to some kind of compensation effect, where the stronger H-bond is harder to make, because it imposes tighter constraints on the geometry of the complex.     

Analysis of multiple H-bond interactions in self-assembly between polyurethane with pendent carboxyl and poly(4-vinylpyridine)

As an extension study, FTIR and molecular simulation methods were combined in the present paper to analyze the H-bond interactions resulting from multiple donors and acceptors that have led to self-assembly based on segmented polyurethane with carboxyl (PUc) and poly(4-vinylpyridine) (P4VP) in our previous work. Of them, FTIR was used to analyze the H-bonding types and interactions as well as their changes before and after self-assembly; molecular mechanics (MM/COMPASS) was used to study the effect of possible conformations on the H-bonds involved and analyze the most probable H-bond patterns; quantum mechanics (QM/B3LYP) was used to help confirm the experimental FTIR band assignments and calculate the H-bond energy. It was found that two types of H-bonds exist, namely, COOH...P4VP (type I) and (OCO)NH...P4VP (type II), based on OH and NH as the strong donors in the interaction between PUc and P4VP. Strong evidence has been obtained for a type II H-bond, which is the specialty in PUc/P4VP assembly. The type I and type II H-bonding energies are -11.293 and -7.150 kcal/mol, respectively. The forming probability of the type I H-bond accounts for 95.87%, while that of the type II H-bond is 4.13%, showing the primary driving force for the assembly based on PUc and P4VP is still the H-bond between COOH and P4VP, yet the H-bonds based on NH and pyridyl in P4VP cannot be ignored.