Peptide bond formation via glycine condensation in the gas phase.

Four unique gas phase mechanisms for peptide bond formation between two glycine molecules have been mapped out with quantum mechanical electronic structure methods. Both concerted and stepwise mechanisms, each leading to a cis and trans glycylglycine product (four mechanisms total), were examined with the B3LYP and MP2 methods and Gaussian atomic orbital basis sets as large as aug-cc-pVTZ. Electronic energies of the stationary points along the reaction pathways were also computed with explicitly correlated MP2-F12 and CCSD(T)-F12 methods. The CCSD(T)-F12 computations indicate that the electronic barriers to peptide bond formation are similar for all four mechanisms (ca. 32-39 kcal mol(-1) relative to two isolated glycine fragments). The smallest barrier (32 kcal mol(-1)) is associated with the lone transition state for the concerted mechanism leading to the formation of a trans peptide bond, whereas the largest barrier (39 kcal mol(-1)) was encountered along the concerted pathway leading to the cis configuration of the glycylglycine dipeptide. Two significant barriers are encountered for the stepwise mechanisms. For both the cis and trans pathways, the early electronic barrier is 36 kcal mol(-1) and the subsequent barrier is approximately 1 kcal mol(-1) lower. A host of intermediates and transition states lie between these two barriers, but they all have very small relative electronic energies (ca. ± 4 kcal mol(-1)). The isolated cis products (glycylglycine + H2O) are virtually isoenergetic with the isolated reactants (within -1 kcal mol(-1)), whereas the trans products are about 5 kcal mol(-1) lower in energy. In both products, however, the water can hydrogen bond to the dipeptide and lower the energy by roughly 5-9 kcal mol(-1). This study indicates that the concerted process leading to a trans configuration about the peptide bond is marginally favored both thermodynamically (exothermic by ca. 5 kcal mol(-1)) and kinetically (barrier height ≈ 32 kcal mol(-1)) according to the CCSD(T)-F12/haTZ electronic energies. The other pathways have slightly larger barrier heights (by 4-8 kcal mol(-1)).

Streptococcal hyaluronate lyase reveals the presence of a structurally significant C-H⋠⋠⋠O hydrogen bond.

Carbon-donated hydrogen bonds (CDHBs) are weak forms of hydrogen bonding (0.5-1.0 kcal mol(-1) ) that are difficult to detect, and thus their roles in the structure and functionality of chemical systems often go unrecognized. Utilizing a computational approach, the existence of a structurally significant CDHB in the medically relevant protein Streptococcus pneumoniae hyaluronate lyase (SpnHL) is affirmed. The structure of a tetrapeptide fragment model containing the CDHB was optimized with second-order perturbation theory. From this, a CDHB with bond distance and angle consistent with previously discovered CDHBs and comparable to neighboring traditional HBs in the fragment model was found. The CDHB competes with another donor T253 OH, whereby the two alternate in strength between protein conformations, imbuing αHelix 3 appreciable flexibility. The CDHB seems to exist in spite of torsional and steric strain on the donor methyl group. It is postulated that the CDHB could aid in either counteracting the macrodipole of αHelix 3 or protecting the A249 CO from destabilizing interactions with the adjacent solvent. Employing the energy gradients from the optimization, the torque generated by the fragment model was computed, which accurately predicts the direction of rotation of αHelix 3 observed from experiment. A strongly correlated motion between αHelix 3 and αHelices 2, 4, and 5 was noted, which the interactions of the fragment model help drive by generating a torque much larger than necessary to rotate just αHelix 3. Considering these results, we conclude that CDHBs should be considered as possible beneficial components of chemical and biological phenomena.

Computational evidence for methyl-donated hydrogen bonds and hydrogen-bond networking in 1,2-ethanediol-dimethyl sulfoxide.

The 1:1 complex of 1,2-ethanediol with dimethyl sulfoxide was studied using density functional theory. A network of three hydrogen bonds holds the complex together, including two in which each methyl group donates to the same hydroxyl oxygen. Four lines of evidence support the existence of methyl-donated hydrogen bonds. The interaction energy is 36 +/- 5 kJ/mol using Becke's three parameter hybrid theory with the 1991 nonlocal correlation functional of Perdew and Wang, and a moderately large basis set (B3PW91/6-311++G**//B3PW91/6-31+G**). To determine the energy of each hydrogen bond, a relaxed potential energy scan was performed in a smaller basis set to break the weaker hydrogen bonds by forced systematic rotation of the methyl groups. Two cross-checking analyses show cooperative effects that cause individual hydrogen bond energies in the network to be nonadditive. When one methyl hydrogen bond is broken, the remaining interactions stabilize the complex by storing an additional 2-3 kJ/mol. With all hydrogen bonds intact, the O[bond]H...O[bond]S hydrogen bond contributes 26 +/- 2 kJ/mol stability, and each weak methyl bond stores 5 +/- 2 kJ/mol.