A new scheme for determining the intramolecular seven-membered ring N–H⋯OC hydrogen-bonding energies of glycine and alanine peptides

A Critical Check for the Role of Resonance in Intramolecular Hydrogen Bonding

Although resonance-assisted H-bonds (RAHBs) are well recognized, the role of π resonance in RAHBs is controversial, as the seemingly enhanced H-bonds in unsaturated compounds may result from the constraints imposed by the σ skeleton. Herein the block-localized wave function (BLW) method, which can derive optimal yet resonance-quenched structures with related physiochemical properties, was employed to examine the correlation between π resonance and the strength of intramolecular RAHBs. Examination of a series of paradigmatic molecules with RAHBs and their saturated analogues showed that it is inappropriate to compare a conjugated system with its saturated counterpart, as they may have quite different σ frameworks. Nevertheless, comparison between a conjugated system and its resonance-quenched (i.e., electron-localized) state, which have identical σ skeletons, shows that in all studied cases, π resonance unanimously reduces the bonding distance by 0.111-0.477 Å, strengthens the bonding by 40-56 %, and redshifts the D-H vibrational frequency by 104-628 cm^-1 . Furthermore, there is an excellent correlation between hydrogen-bonding strength and the classical Coulomb attraction between the hydrogen-bond donor and the acceptor, which suggests that the dominant role of the electrostatic interaction in H-bonds and RAHBs originates from the charge flow from H-bond donors to acceptors through π conjugation.

The Origin of the Non-Additivity in Resonance-Assisted Hydrogen Bond Systems

The concept of resonance-assisted hydrogen bond (RAHB) has been widely accepted, and its impact on structures and energetics can be best studied computationally using the block-localized wave function (BLW) method, which is a variant of ab initio valence bond (VB) theory and able to derive strictly electron-localized structures self-consistently. In this work, we use the BLW method to examine a few molecules that result from the merging of two malonaldehyde molecules. As each of these molecules contains two hydrogen bonds, these intramolecular hydrogen bonds may be cooperative or anticooperative, depended on their relative orientations, and compared with the hydrogen bond in malonaldehyde. Apart from quantitatively confirming the concept of RAHB, the comparison of the computations with and without π resonance shows that both σ-framework and π-resonance contribute to the nonadditivity in these RAHB systems with multiple hydrogen bonds.

Rapid evaluation of the interaction energies for carbohydrate-containing hydrogen-bonded complexes via the polarizable dipole–dipole interaction model combined with NBO or AM1 charge

The polarizable dipole–dipole interaction model has been developed to rapidly and accurately estimate the hydrogen bond distances and interaction energies for carbohydrate-containing hydrogen-bonded complexes.

ESTIMATION ON THE INTRAMOLECULAR 8- AND 12-MEMBERED RING N–H…O=C HYDROGEN BONDING ENERGIES IN β-PEPTIDES

Computation of accurate hydrogen bonding energies in peptides is of great importance in understanding the conformational stabilities of peptides. In this paper, the intramolecular 8- and 12-membered ring N – H … O = C hydrogen bonding energies in β-peptide structures were evaluated. The optimal structures of the β-peptide conformers were obtained using MP2/6-31G(d) method. The MP2/6-311++G(d,p) calculations were then carried out to evaluate the single-point energies. The results show that the intramolecular 8-membered ring N – H … O = C hydrogen bonding energies in the five β-dipeptide structures β-di, β-di-R1, β-di-R2, β-di-R3, and β-di-R4 are -5.50, -5.40, -7.28, -4.94, and -6.84 kcal/mol with BSSE correction, respectively; the intramolecular 12-membered ring N – H … O = C hydrogen bonding energies in the nine β-tripeptide structures β-tri, β-tri-R1, β-tri-R2, β-tri-R3, β-tri-R4, β-tri-R1', β-tri-R2', β-tri-R3' and β-tri-R4' are -10.23, -10.32, -9.53, -10.30, -10.32, -10.55, -10.09, -10.51, and -9.60 kcal/mol with BSSE correction, respectively. Our calculation results further indicate that for the intramolecular 8-membered ring hydrogen bondings, the structures where the orientation of the side chain methyl group is "a–a" have stronger intramolecular hydrogen bondings than those where the orientation of the side chain methyl group is "e–e", while for the intramolecular 12-membered ring hydrogen bondings, the structures where the orientation of the side chain methyl group is "e–e" have stronger intramolecular hydrogen bondings than those where the orientation of the side chain methyl group is "a–a". The method is also applied to estimate the individual intermolecular hydrogen bonding energies in the dimers of amino-acetaldehyde, 2-amino-acetamide, 2-oxo-acetamide, 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.71, -1.50, -4.67, and -3.22 kcal/mol at the MP2/6-311++G(d,p) level, which are in good agreement with the values of -1.84, -1.72, -4.93, and -3.26 kcal/mol predicted by the supermolecular method.

Rapid evaluation of the binding energies in hydrogen-bonded amide-thymine and amide-uracil dimers in gas phase

The binding energies and the equilibrium hydrogen bond distances as well as the potential energy curves of 48 hydrogen-bonded amide-thymine and amide-uracil dimers are evaluated from the analytic potential energy function established in our lab recently. The calculation results show that the potential energy curves obtained from the analytic potential energy function are in good agreement with those obtained from MP2/6-311+G** calculations by including the BSSE correction. For all the 48 dimers, the analytic potential energy function yields the binding energies of the MP2/6-311+G** with BSSE correction within the error limits of 0.50 kcal/mol for 46 dimers, only two differences are larger than 0.50 kcal/mol and the largest one is only 0.60 kcal/mol. The analytic potential energy function produces the equilibrium hydrogen bond distances of the MP2/6-311+G** with BSSE correction within the error limits of 0.050 Å for all the 48 dimers. The analytic potential energy function is further applied to four more complicated hydrogen-bonded amide-base systems involving amino acid side chain and β-sheet. The values of the binding energies and equilibrium hydrogen bond distances obtained from the analytic potential energy function are also in good agreement with those obtained from MP2 calculations with the BSSE correction. These results demonstrate that the analytic potential energy function can be used to evaluate the binding energies in hydrogen-bonded amide-base dimers quickly and accurately.

A scheme for rapid prediction of cooperativity in hydrogen bond chains of formamides, acetamides, and N-methylformamides

A scheme is proposed in this article to predict the cooperativity in hydrogen bond chains of formamides, acetamides, and N-methylformamides. The parameters needed in the scheme are derived from fitting to the hydrogen bonding energies of MP2/6-31+G** with basis set superposition error (BSSE) correction of the hydrogen bond chains of formamides containing from two to eight monomeric units. The scheme is then used to calculate the individual hydrogen bonding energies in the chains of formamides containing 9 and 12 monomeric units, in the chains of acetamides containing from two to seven monomeric units, in the chains of N-methylformamides containing from two to seven monomeric units. The calculation results show that the cooperativity predicted by the scheme proposed in this paper is in good agreement with those obtained from MP2/6-31+G** calculations by including the BSSE correction, demonstrating that the scheme proposed in this article is reasonable. Based on our scheme, a cooperativity effect of almost 240% of the dimer hydrogen bonding energy in long hydrogen bond formamide chains, a cooperativity effect of almost 190% of the dimer hydrogen bonding energy in long hydrogen bond acetamide chains, and a cooperativity effect of almost 210% of the dimer hydrogen bonding energy in long hydrogen bond N-methylformamide chains are predicted. The scheme is further applied to some heterogeneous chains containing formamide, acetamide, and N-methylformamide. The individual hydrogen bonding energies in these heterogeneous chains predicted by our scheme are also in good agreement with those obtained from Møller-Plesset calculations including BSSE correction.

Structural selection of ionic-complementary peptides with electrostatic interactions

The structures of the peptides and their assembly are largely modulated by the environment. To discover the physical principles governing the structural modulations of peptides by the environment would be useful for many applications. As the typical examples, the structures of three kinds of ionic-complementary EAK16-family peptides under various environmental conditions are studied with simulations in this work. A model with intermediate resolution is used, in which both the backbone hydrogen bonds and electrostatic interactions are explicitly considered. The thermodynamics of these peptides (including the free energy and heat capacity) are described for various strengths of the electrostatic interactions which reflect the variation of environment. With these results, the phase diagrams of these peptides related to the temperature and the strength of electrostatic interactions are presented and compared. Based on the differences in the phase structures of the peptide, the different aggregation behaviors are explained based on the monomeric structural features of the peptides. Through the analysis on the stability of various secondary structures of these peptides, it is demonstrated that the charge pattern is the basic reason of the different responses of the EAK16-family peptides to the environmental changes. These results provide some examples and insights for the principles of structural selection by environment and may be helpful for further analysis and designs of peptide systems.

Estimation of N-H...O=C intramolecular hydrogen bond energy in polypeptides

The previously proposed molecular tailoring approached (MTA) [Deshmukh, M. M.; Gadre, S. R.; Bartolotti, L. J. J. Phys. Chem. A 2006, 110, 12519] for the estimation of intramolecular O-H...O hydrogen bond energy is extended to that for the N-H...O=C bond within polypeptides. The methodology is initially tested on a tetrapeptide containing two types of N-H...O=C hydrogen bonds and is found to distinguish between them. The estimated values are in good agreement with the trends predicted by the geometrical parameters. Furthermore, this methodology is applied to partially as well as fully substituted, capped polyglycines that contain five glycine residues (acetyl-(gly)(5)-NH(2)) to check the effect of substituents on the energetics of hydrogen bonds. The estimated N-H...O=C bond energy values lie in the range of 4-6 kcal/mol. These estimated values are not only in concurrence with the geometric parameters but also able to reflect the subtle effects of substituents for the substituted polypeptides studied in the present work.

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.

Contributions of NH...O and CH...O hydrogen bonds to the stability of beta-sheets in proteins

Ab initio quantum calculations are applied to both the parallel and the antiparallel arrangements of the beta-sheets of proteins. The energies of the NH...O and CH...O hydrogen bonds present in the beta-sheet are evaluated separately from one another by appropriate modifications of the model systems. The bond energies of these two sorts of hydrogen bonds are found to be very nearly equal in the parallel beta-sheet. The NH...O bonds are stronger than CH...O in the antiparallel geometry but only by a relatively small margin. Moreover, the former NH...O bonds are weakened when placed next to one another, as occurs in the antiparallel beta-sheet. As a result, there is little energetic distinction between the NH...O and CH...O bonds in the full antiparallel beta-sheet, just as in the parallel structure.

Finding the transition state without initial guess: The growing string method for Newton trajectory to isomerization and enantiomerization reaction of alanine dipeptide and poly(15)alanine

We report a new, high-dimensional application of a method for finding a transition state (TS) between a reactant and a product on the potential energy surface: the search of a growing string along a reaction path defined by any Newton trajectory in combination with the Berny method (Quapp, J Chem Phys (2005), 122, 174106; we have provided this algorithm on a web page). Two given minima are connected by a one-dimensional, but usually curvilinear reaction coordinate. It leads to the TS region. The application of the method to alanine dipeptide finds the TS of the isomerisation C(7 ax) --> C(5), some TSs of the enantiomerisation of C(7 ax) from L-form to quasi-D-form, and it finds the TS region of a transition of a partly unfolded, bent structure which turns back into a mainly alpha-helix in the Ac(Ala)(15)NHMe polyalanine (all at the quantum mechanical level B3LYP/6-31G: the growing string calculation is interfaced with the Gaussian03 package). The formation or dissolvation of some alpha- or 3(10)-hydrogen bonds of the helix are discussed along the TS pathway, as well as the case of an enantiomer at the central residue of the helix.