Electronic properties of hydrogen-bonded complexes of benzene(HCN)(1-4): comparison with benzene(H2O)(1-4)

The electronic properties, specifically, the dipole and quadrupole moments and the ionization energies of benzene (Bz) and hydrogen cyanide (HCN), and the respective binding energies, of complexes of Bz(HCN)(1-4), have been studied through MP2 and OVGF calculations. The results are compared with the properties of benzene-water complexes, Bz(H(2)O)(1-4), with the purpose of analyzing the electronic properties of microsolvated benzene, with respect to the strength of the CH/π and OH/π hydrogen-bond (H-bond) interactions. The linear HCN chains have the singular ability to interact with the aromatic ring, preserving the symmetry of the latter. A blue shift of the first vertical ionization energies (IEs) of benzene is observed for the linear Bz(HCN)(1-4) clusters, which increases with the length of the chain. NBO analysis indicates that the increase of the IE with the number of HCN molecules is related to a strengthening of the CH/π H-bond, driven by cooperative effects, increasing the acidity of the hydrogen cyanide H atom involved in the π H-bond. The longer HCN chains (n ≥ 3), however, can bend to form CH/N H-bonds with the Bz H atoms. These cyclic structures are found to be slightly more stable than their linear counterparts. For the nonlinear Bz(HCN)(3-4) and Bz(H(2)O)(2-4) complexes, an increase of the binding energy with the number of solvent molecules and a decrease of the IE of benzene, relative to the values for the Bz(HCN) and Bz(H(2)O) complexes, respectively, are observed. Although a strengthening of the CH/π and OH/π H-bonds, with increasing n, also takes place for the Bz(H(2)O)(2-4) and Bz(HCN)(3-4) nonlinear complexes, Bz proton donor, CH/O, and CH/N interactions are at the origin of this decrease. Thus CH/π and OH/π H-bonds lead to higher IEs of Bz, whereas the weaker CH/N and CH/O H-bond interactions have the opposite effect. The present results emphasize the importance of both aromatic XH/π (X = C, O) and CH/X (X = N, O) interactions for understanding the structure and electronic properties of Bz(HCN)(n) and Bz(H(2)O)(n) complexes.     

Effect of solvent upon CH...O hydrogen bonds with implications for protein folding

The series of CH...O bonds formed between CF(n)H(4-n) (n = 0-3) and water are studied by quantum calculations under vacuum and in various solvents, including aqueous environment. The results are compared with the OH...O bond of the water dimer in the same solvents. Increasing polarity of the solvent leads in all cases to a lessening of the H-bond interaction energy, in a uniform fashion such that the CH...O bonds all remain weaker than OH...O in any solvent. These H-bond weakenings are coupled to a shortening of the inter-subunit separation. The contraction of the covalent CH bond to the bridging proton is reduced as the solvent becomes more polar, and the blue shift of its stretching vibration is likewise diminished. A process is considered that simulates protein folding by starting from a pair of noninteracting subunits in aqueous solvent and then goes to a H-bonded pair within the confines of a protein environment. This process is found to be energetically more favorable for some of the CH...O H-bonds than for the nominally stronger conventional OH...O H-bond. This finding suggests that CH...O bonds can make important energetic contributions to protein folding, on par with those made by traditional H-bonds.     

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.     

Weak H-bonds. Comparisons of CH···O to NH···O in proteins and PH···N to direct P···N interactions

Whereas CH···O H-bonds are usually weaker than interpeptide NH···O H-bonds, this is not necessarily the case within proteins. The nominally weaker CH···O are surprisingly strong, comparable to, and in some cases stronger than, the NH···O H-bonds in the context of the forces that hold together the adjacent strands in protein β-sheets. The peptide NH is greatly weakened as proton donor in certain conformations of the protein backbone, particularly extended structures, and forms correspondingly weaker H-bonds. The PH group is a weak proton donor, but will form PH···N H-bonds. However, there is a stronger interaction in which P can engage, in which the P atom, not the H, directly approaches the N electron donor to establish a direct P···N interaction. This approach is stabilized by the same sort of electron transfer from the N lone pair to the P-H σ* antibond that characterizes the PH···N H-bond.     

Microsolvation effect, hydrogen-bonding pattern, and electron affinity of the uracil-water complexes U-(H2O)n (n = 1, 2, 3)

To achieve a systematic understanding of the influence of microsolvation on the electron accepting behaviors of nucleobases, the reliable theoretical method (B3LYP/DZP++) has been applied to a comprehensive conformational investigation on the uracil-water complexes U-(H(2)O)(n) (n = 1, 2, 3) in both neutral and anionic forms. For the neutral complexes, the conformers of hydration on the O2 of uracil are energetically favored. However, hydration on the O4 atom of uracil is more stable for the radical anions. The electron structure analysis for the H-bonding patterns reveal that the CH...OH(2) type H-bond exists only for di- and trihydrated uracil complexes in which a water dimer or trimer is involved. The electron density structure analysis and the atoms-in-molecules (AIM) analysis for U-(H(2)O)(n) suggest a threshold value of the bond critical point (BCP) density to justify the CH...OH(2) type H-bond; that is, CH...OH(2) could be considered to be a H-bond only when its BCP density value is equal to or larger than 0.010 au. The positive adiabatic electron affinity (AEA) and vertical detachment energy (VDE) values for the uracil-water complexes suggest that these hydrated uracil anions are stable. Moreover, the average AEA and VDE of U-(H(2)O)(n) increase as the number of the hydration waters increases.     

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.     

Molecular basis of the recognition process: hydrogen-bonding patterns in the guanine primary recognition site of ribonuclease T1

Investigation of the intrinsic H-bonding pattern of the guanine complex with a sizable segment (from Asn43 to Glu46) of the primary recognition site (PRS) in RNase T1 at the B3LYP/6-311G(d,p) level of theory enables the electronic density characteristics of the H-bonding patterns of the guanine-PRS complexes to be identified. The perfect H-bonding pattern in the guanine recognition site is achieved through the guanine complex interactions with the large segment of the PRS. Two significant short H-bonds, O epsilon 1...HN1 and O epsilon 2...HN2, have been identified. The similar short H-bond distances found in the anionic GC- base pair and in this study suggest that the short hydrogen-bond distances may be characteristic of the multiple H-bonded anionic nucleobases. The H-bonding energy distribution, the geometric analysis of the H-bonding pattern, and the electron structure characteristics of the H-bonds in the guanine PRS of RNase T1 all suggest that the O epsilon 1...HN1 and O epsilon 2...HN2 side-chain H-bonds dominate the binding at the guanine recognition site of RNase T1. Also, the geometry evidence, the electron structure characteristics, and the properties of the bond critical points of the H-bonds reveal that the side-chain H-bonding and the main-chain H-bonding are mutually intensifying. Thus the positive cooperativity between Asn43 to Tyr45 and Glu46 is proposed.     

Comparison of various types of hydrogen bonds involving aromatic amino acids

Ab initio calculations are used to compare the abilities of the aromatic groups of the Phe, Tyr, Trp, and His amino acids (modeled respectively by benzene, phenol, indole, and imidazole) to form H-bonds of three different types. Strongest of all are the conventional H-bonds (e.g., OH..O and OH..N). His forms the strongest such H-bond, followed by Tyr, and then by Trp. Whereas OH..phi bonds formed by the approach of a proton donor to the pi electron cloud above the aromatic system are somewhat weaker, they nonetheless represent an important class of stabilizing interactions. The strengths of H-bonds in this category follow the trend Trp > His > Tyr approximately Phe. CH.O interactions are weaker still, and only those involving His and Trp are strong enough to make significant contributions to protein structure. A protonated residue such as HisH(+) makes for a very powerful proton donor, such that even its CH..O H-bonds are stronger than the conventional H-bonds formed by neutral groups.     

Intramolecular hydrogen bonding and cooperative interactions in carbohydrates via the molecular tailoring approach

In spite of many theoretical and experimental attempts for understanding intramolecular hydrogen bonding (H-bonding) in carbohydrates, a direct quantification of individual intramolecular H-bond energies and the cooperativity among the H-bonded networks has not been reported in the literature. The present work attempts, for the first time, a direct estimation of individual intramolecular O-H...O interaction energies in sugar molecules using the recently developed molecular tailoring approach (MTA). The estimated H-bond energies are in the range of 1.2-4.1 kcal mol(-1). It is seen that the OH...O equatorial-equatorial interaction energies lie between 1.8 and 2.5 kcal mol(-1), with axial-equatorial ones being stronger (2.0-3.5 kcal mol(-1)). The strongest bonds are nonvicinal axial-axial H-bonds (3.0-4.1 kcal mol(-1)). This trend in H-bond energies is in agreement with the earlier reports based on the water-water H-bond angle, solvent-accessible surface area (SASA), and (1)H NMR analysis. The contribution to the H-bond energy from the cooperativity is also estimated using MTA. This contribution is seen to be typically between 0.1 and 0.6 kcal mol(-1) when H-bonds are a part of a relatively weak equatorial-equatorial H-bond network and is much higher (0.5-1.1 kcal mol(-1)) when H-bonds participate in an axial-axial H-bond network.     

Relative strengths of NH..O and CH..O hydrogen bonds between polypeptide chain segments

Correlated ab initio calculations are used to compare the energetics when the CH and NH groups of the model dipeptide CHONHCH2CONH2 are each allowed to form a H-bond with the proton acceptor O of a peptide group. When the dipeptide is in its C7 conformation, the NH..O H-bond energy is found to be 7.4 kcal/mol, as compared to only 2.8 kcal/mol for the CH..O interaction. On the other hand, the situation reverses, and the CH..O H-bond becomes stronger than NH..O, when the dipeptide adopts a C5 structure. This reversal is important as C5 is nearly equal in stability to C7 for the dipeptide, and is representative of the commonly observed beta-sheet structure in a protein. Immersing the dipeptide-peptide pair in a model solvent weakens both sorts of H-bonds, and in a fairly uniform manner. Consequently, the trends observed in the in vacuo situation retain their validity in either aqueous solution or the protein interior. Likewise, the desolvation penalty, suffered by removing a H-bonded complex from water and placing it in the less polar interior of a protein, is quite similar for the NH..O and CH..O bonds.     

The strength with which a peptide group can form a hydrogen bond varies with the internal conformation of the polypeptide chain

The strength of the H-bond formed between a dipeptide and a proton acceptor molecule is assessed by correlated ab initio quantum calculations for a broad range of different conformations of the dipeptide. The H-bond energy is very sensitive to the internal (phi,psi) angles, even when the geometry of the H-bond does not vary significantly from one conformation to another. This result indicates that the peptide NH is a much less potent proton donor in certain conformations than in others. In particular, extended conformations of a polypeptide are capable of only weak H-bonds. Thus, the interstrand NH...O H-bonds in parallel and antiparallel beta-sheets are expected to be significantly weaker than those found in other conformations, such as helices, ribbons, and beta-bends, even if the H-bond geometries are similar.     

Electronic properties of multifurcated bent hydrogen bonds CH3...Y and CH2...Y

H-bonding angle angleYHX has an important effect on the electronic properties of the H-bond Y...HX, such as intra- and intermolecular hyperconjugations and rehybridization, and topological properties of electron density. We studied the multifurcated bent H-bonds of the proton donors H3CZ (Z = F, Cl, Br), H2CO and H2CF2 with the proton acceptors Cl(-) and Br(-) at the four high levels of theory: MP2/6-311++G(d,p), MP2/6-311++G(2df,2p), MP2/6-311++G(3df,3pd) and QCISD/6-311++G(d,p), and found that they are all blue-shifted. These complexes have large interaction energies, 7-12 kcal mol(-1), and large blue shifts, delta r(HC) = -0.0025 --0.006 A and delta v(HC) = 30-90 cm(-1). The natural bond orbital analysis shows that the blue shifts of these H-bonds Y...HnCZ are mainly caused by three factors: rehybridization; indirect intermolecular hyperconjugation n(Y) -->sigma*(CZ), in that the electron density from n(Y) of the proton acceptor is transferred not to sigma*(CH), but to sigma*(CZ) of the donor; intramolecular hyperconjugation n(Z) -->sigma*(CH), in that the electron density in sigma*(CH) comes back to n(Z) of the donor such that the occupancy in sigma*(CH) decreases. The topological properties of the electron density of the bifurcated H-bonds Y...H2CZ are similar to those of the usual linear H-bonds, there is a bond critical point between Y and each hydrogen, and a ring critical point inside the tetragon YHCH. However, the topological properties of electron density of the trifurcated H-bonds Y...H3CZ are essentially different from those of linear H-bonds, in that the intermolecular bond critical point, which represents a closed-shell interaction, is not between Y and hydrogen, but between Y and carbon.     

Preferential formation of the different hydrogen bonds and their effects in tetrahydrofuran and tetrahydropyran microhydrated complexes

The role of cycloether-water (c-w) and water-water (w-w) hydrogen bonds (H-bonds) on the stability of the tetrahydrofuran THF/(H(2)O)(n) and the tetrahydropyran THP/(H(2)O)(n) complexes with n = 1-4 was investigated herein using the density functional and ab initio methods and the atoms in molecules theory. Geometry optimizations for these complexes were carried out with various possible initial guess structures. It was revealed that the major contributions of the mono and dihydrated complexes came from c-w H-bonds. A competition between c-w and w-w H-bonds contribution was observed for trihydrated complexes. For most of tetrahydrated complexes, the inter-water H-bonds provided the greatest contribution, whereas the c-w contributions were small but not negligible. It was confirmed that to produce a hydrophobic hydration of cycloethers, the C-H···O(w) H-bond should be associated with a network of H-bonds that connects both portions of the solute, through the formation of a bifunctional H-bond. A linear correlation is obtained for the sum of electron density at the bond critical points (ρ(b)) with the interaction energy (ΔE) and with the solute-solvent interaction energy (ΔE(s-w)) of the microhydrated complexes. In addition, a new way to estimate the energetic contribution as well as the preferential formation of the different H-bonds based completely on ρ(b) was found. Even more, it allows to differentiate the contribution from c-w interactions in both hydrophilic and hydrophobic contributions, it is therefore a useful tool for studying the hydration of large biomolecules. The analysis of the modifications in the atomic and group properties brought about by successive addition of H(2)O molecules allowed to pinpoint the atoms or molecular groups that undergo the greatest changes in electron population and energetic stabilization. It was identified that the remarkable stabilization of the water oxygen atoms is crucial for the stabilization of the complexes.     

Hydrogen bonding interactions between <Emphasis Type="Italic">N</Emphasis>,<Emphasis Type="Italic">N</Emphasis>-dimethylformamide and cysteine: DFT studies of structures,properties, and topologies

The hydrogen bonding interactions between cysteine and N,N-dimethylformamide (DMF) were studied at the extended hybrid functional DFT-X3LYP/6-311++G(d,p) level regarding their geometries, energies, vibrational frequencies, and topological features of the electron density. The quantum theory of atoms in molecules (QTAIM) and natural bond orbital (NBO) analyses were employed to elucidate the interaction characteristics in the complexes. The results show that two intermolecular hydrogen bonds (H-bonds) are formed in one complex except few complexes with one intermolecular H-bond. The H-bonds involving O atom of DMF as H-bond acceptor usually are red-shifting H-bonds, while the blue-shifting H-bond usually involve methyl of DMF or methenyl of cysteine moiety as H-bond donors. Both hydrogen bonding interaction and structural deformation play important roles in the relative stabilities of the complexes. Due to the π-bond cooperativity, the strongest H-bond is formed between hydroxyl of cysteine moiety and O atom of DMF, however, the serious deformation counteract the hydrogen bonding interaction to a great extent. The complex involves a stronger hydrogen bonding interaction as well as the smaller deformation is the most stable one. The electron density (ρ_b) as well as its Laplacian (∇²ρ_b) at the H-bond critical point predicted by QTAIM is strongly correlated with the H-bond structural parameter (δR _H···Y) and the second-perturbation energies E(2) in the NBO scheme.     

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.     

Density functional theory and topological analysis on the hydrogen bonding interactions in N-protonated adrenaline–DMSO complexes

In this study, complexes formed via hydrogen bond interactions between N-protonated adrenaline (AdH⁺) and DMSO have been studied by density functional theory (DFT). The relevant geometries, energies, and IR characteristics of the hydrogen bonds (H-bonds) have been systematically investigated. The natural bond orbital (NBO) and the quantum theory of atoms in molecule (QTAIM) analysis have also been applied to understand the nature of the hydrogen bonding interactions in complexes. The H-bonds involving amino or hydroxyls as H-donor are dominant H-bonds in complexes and are attributed to strong H-bonds. The weak H-bonds, such as π H-bonds and H-bonds involving methyl (DMSO) or methenyls (C2H6 and C5H7 of AdH⁺) as H-acceptors, were found in complexes as well. The complexes in which the dominant H-bond involves amino of AdH⁺ as H-donor are more stable than those with the dominant H-bond involving hydroxyls as H-donor. Some relationships between various properties of QTAIM, NBO, geometry as well as frequency were also investigated.     

Hydrogen-bonding interaction in a complex of amino acid with urea studied by DFT calculations

Theoretical studies on hydrogen-bonded complexes between amino acids (glycine, alanine, and leucine) and urea in gas phase have been carried out using density functional theory (DFT) and ab initio methods at the B3LYP/6-311++g** and MP2/6-311++g** theory levels. The structures, binding energy, Chelpg (charges from electrostatic potentials using a grid-based method) charge distribution, and bond characteristics of the mentioned complexes were calculated. Urea is a good H-bond donor and an excellent receptor for highly electronegative atoms like O and N, through the formation of two or more hydrogen bonds. The NH₂ and COOH groups of amino acids can form several different types of H-bonds with urea molecular, as well as C^αH and alkyl side chains. The calculated high binding energy also suggests multiple H-bonds formed in one complex. The OH···O contact is the strongest hydrogen bond interaction with H···O separation around 1.65 Å and its relevant angle close to 176°. The closely linear amide H-bonds NH···O and OH···N strongly stabilize the amino acid–urea complex with H···O separation between 1.89 and 2.38 Å. The weaker CH···O/N H-bonds are also discussed as significant interaction in biological systems involving amino acids.     

Influence of H-bond strength on chelate cooperativity

Intermolecular complexes formed between metalloporphyrins and pyridine ligands equipped with multiple H-bond donors and acceptors have been used to measure the free energy contributions due to intramolecular ether-phenol H-bonding in the 24 different supramolecular architectures using chemical double mutant cycles in toluene. The ether-phenol interactions are relatively weak, and there are significant populations of partially bound states where between zero and four intramolecular H-bonds are made in addition to the porphyrin-ligand coordination interaction. The complexes were analyzed as ensembles of partially bound states to determine the effective molarities for the intramolecular interactions by comparison with the corresponding intermolecular ether-phenol H-bonds. The properties of the ether-phenol interactions were compared with phosphonate diester-phenol interactions in a closely related ligand system, which has more powerful H-bond acceptor oxygens positioned at the same location on the ligand framework. This provides a comparison of the properties of weak and strong H-bonds embedded in the same 24 supramolecular architectures. When the product of the intermolecular association constant and the effective molarity KEM > 1, there is a linear increase in the free energy contribution due to H-bonding with log EM, because the intramolecular interactions contribute fully to the stability of the complex. When KEM < 1, the H-bonded state is not significantly populated, and there is no impact on the overall stability of the complex. Intermolecular phosphonate diester-phenol H-bonds are 2 orders of magnitude stronger than ether-phenol H-bonds in toluene, so for the phosphonate diester ligand system, 23 of the 24 supramolecular architectures make intramolecular H-bonds. However, only 8 of these architectures lead to detectable H-bonding in the ether ligand system. The other 15 complexes have a suitable geometry for formation of H-bonds, but the ether-phenol interaction is not strong enough to overcome the reorganization costs associated with making intramolecular contacts, i.e., KEM < 1 for the ether ligands, and KEM > 1 for the phosphonate diester ligands. The values of EM measured for two different types of H-bond acceptor are linearly correlated, which suggests that EM is a property of the supramolecular acrchitecture. However, the absolute value of EM for an intramolecular phosphonate diester H-bond is about 4 times lower than the corresponding value for an intramolecular ether-phenol interaction embedded in the same supramolecular framework, which suggests that there may be some interplay of K and EM.     

Protein–ligand interfaces are polarized: discovery of a strong trend for intermolecular hydrogen bonds to favor donors on the protein side with implications for predicting and designing ligand complexes

Understanding how proteins encode ligand specificity is fascinating and similar in importance to deciphering the genetic code. For protein–ligand recognition, the combination of an almost infinite variety of interfacial shapes and patterns of chemical groups makes the problem especially challenging. Here we analyze data across non-homologous proteins in complex with small biological ligands to address observations made in our inhibitor discovery projects: that proteins favor donating H-bonds to ligands and avoid using groups with both H-bond donor and acceptor capacity. The resulting clear and significant chemical group matching preferences elucidate the code for protein-native ligand binding, similar to the dominant patterns found in nucleic acid base-pairing. On average, 90% of the keto and carboxylate oxygens occurring in the biological ligands formed direct H-bonds to the protein. A two-fold preference was found for protein atoms to act as H-bond donors and ligand atoms to act as acceptors, and 76% of all intermolecular H-bonds involved an amine donor. Together, the tight chemical and geometric constraints associated with satisfying donor groups generate a hydrogen-bonding lock that can be matched only by ligands bearing the right acceptor-rich key. Measuring an index of H-bond preference based on the observed chemical trends proved sufficient to predict other protein–ligand complexes and can be used to guide molecular design. The resulting Hbind and Protein Recognition Index software packages are being made available for rigorously defining intermolecular H-bonds and measuring the extent to which H-bonding patterns in a given complex match the preference key.     

Attractive strain: the disadvantages of rigid multiple H-bond donors and acceptors. A theoretical analysis of the hydrogen-bonding interactions in complexes of tetraazaanthracenedione with pyridylureas

We report density functional theory calculations at the B3LYP/D95(d,p) level on the hydrogen-bonding complexes of tetraazaanthracenedione, 1, with N-(pyridin-2-yl)urea, 2H, or N-(6-aminopyridin-2-yl)urea, 2N. The interaction energy of the 1-2H complex exceeds that of 1-2N, despite the fact that 1-2N contains a strong N-H...O interaction in place of a weak C-H...O interaction in 1-2H. We show that the 1-2N interaction is weaker than the sum of the four normal individual H-bonding interactions because the steric constraints of the complex prevent the H-bonding donors and acceptors from optimally approaching each other to form the two central H-bonds. This steric phenomenon, which we call attractive strain, is likely present to some extent in most H-bonding systems that contain more than two H-bonds between rigid monomers. Attractive strain is unusually important in 1-2N. Attractive strain can be conceived of as an enthalpic cost for the entropic benefits of freezing the dihedral angles of the multiple H-bond donors and acceptors by designing rigid systems.