Dihydrogen Bonding-Seen through the Eyes of Vibrational Spectroscopy.

In this work, we analyzed five groups of different dihydrogen bonding interactions and hydrogen clusters with an H3+ kernel utilizing the local vibrational mode theory, developed by our group, complemented with the Quantum Theory of Atoms-in-Molecules analysis to assess the strength and nature of the dihydrogen bonds in these systems. We could show that the intrinsic strength of the dihydrogen bonds investigated is primarily related to the protonic bond as opposed to the hydridic bond; thus, this should be the region of focus when designing dihydrogen bonded complexes with a particular strength. We could also show that the popular discussion of the blue/red shifts of dihydrogen bonding based on the normal mode frequencies is hampered from mode-mode coupling and that a blue/red shift discussion based on local mode frequencies is more meaningful. Based on the bond analysis of the H3+(H2)n systems, we conclude that the bond strength in these crystal-like structures makes them interesting for potential hydrogen storage applications.

On the formation of CN bonds in Titan's atmosphere-a unified reaction valley approach study.

In this work, we investigated the formation of protonated hydrogen cyanide HCNH+ and methylene amine cation CH[Formula: see text] (both identified in Titan's upper atmosphere) from three different pathways which stem from the interaction between CH4 and N+(3P). As a mechanistic tool, we used the Unified Reaction Valley Approach (URVA) complemented with the Local Mode Analysis (LMA) assessing the strength of the CN bonds formed in these reactions. Our URVA studies could provide a comprehensive overview on bond formation/cleavage processes relevant to the specific mechanism of eight reactions R1- R8 that occur across the three pathways. In addition, we could explain the formation of CH[Formula: see text] and the appearance of HCNH+ and CHNH[Formula: see text] along these paths. Although only smaller molecules are involved in these reactions including isomerization, hydrogen atom abstraction, and hydrogen molecule capture, we found a number of interesting features, such as roaming in reaction R3 or the primary interaction of H2 with the carbon atom in HCNH+ in reaction R8 followed by migration of one of the H2 hydrogen atoms to the nitrogen which is more cost effective than breaking the HH bond first; a feature often found in catalysis. In all cases, charge transfer between carbon and nitrogen could be identified as a driving force for the CN bond formation. As revealed by LMA, the CN bonds formed in reactions R1-R8 cover a broad bond strength range from very weak to very strong, with the CN bond in protonated hydrogen cyanide HCNH+ identified as the strongest of all molecules investigated in this work. Our study demonstrates the large potential of both URVA and LMA to shed new light into these extraterrestrial reactions to help better understand prebiotic processes as well as develop guidelines for future investigations involving areas of complex interstellar chemistry. In particular, the formation of CN bonds as a precursor to the extraterrestrial formation of amino acids will be the focus of future investigations. Formation of CN bonds in Titan's atmosphere visualized via the reaction path curvature.

Hydrogen Bonding in Natural and Unnatural Base Pairs-A Local Vibrational Mode Study.

In this work hydrogen bonding in a diverse set of 36 unnatural and the three natural Watson Crick base pairs adenine (A)-thymine (T), adenine (A)-uracil (U) and guanine (G)-cytosine (C) was assessed utilizing local vibrational force constants derived from the local mode analysis, originally introduced by Konkoli and Cremer as a unique bond strength measure based on vibrational spectroscopy. The local mode analysis was complemented by the topological analysis of the electronic density and the natural bond orbital analysis. The most interesting findings of our study are that (i) hydrogen bonding in Watson Crick base pairs is not exceptionally strong and (ii) the N-H⋯N is the most favorable hydrogen bond in both unnatural and natural base pairs while O-H⋯N/O bonds are the less favorable in unnatural base pairs and not found at all in natural base pairs. In addition, the important role of non-classical C-H⋯N/O bonds for the stabilization of base pairs was revealed, especially the role of C-H⋯O bonds in Watson Crick base pairs. Hydrogen bonding in Watson Crick base pairs modeled in the DNA via a QM/MM approach showed that the DNA environment increases the strength of the central N-H⋯N bond and the C-H⋯O bonds, and at the same time decreases the strength of the N-H⋯O bond. However, the general trends observed in the gas phase calculations remain unchanged. The new methodology presented and tested in this work provides the bioengineering community with an efficient design tool to assess and predict the type and strength of hydrogen bonding in artificial base pairs.

Interplay of Ring Puckering and Hydrogen Bonding in Deoxyribonucleosides.

The Cremer-Pople ring puckering analysis and the Konkoli-Cremer local mode analysis supported by the topological analysis of the electron density were applied for the first comprehensive analysis of the interplay between deoxyribose ring puckering and intramolecular H-bonding in 2'-deoxycytidine, 2'-deoxyadenosine, 2'-deoxythymidine, and 2'-deoxyguanosine. We mapped for each deoxyribonucleoside the complete conformational energy surface and the corresponding pseudorotation path. We found only incomplete pseudorotation cycles, caused by ring inversion, which we coined as pseudolibration paths. On each pseudolibration path a global and a local minimum separated by a transition state were identified. The investigation of H-bond free deoxyribonucleoside analogs revealed that removal of the H-bond does not restore the full conformational flexibility of the sugar ring. Our work showed that ring puckering predominantly determines the conformational energy; the larger the puckering amplitude, the lower the conformational energy. In contrast no direct correlation between conformational energy and H-bond strength was found. The longest and weakest H-bonds are located in the local minimum region, whereas the shortest and strongest H-bonds are located outside the global and local minimum regions at the turning points of the pseudolibration paths, i.e., H-bonding determines the shape and length of the pseudolibration paths. In addition to the H-bond strength, we evaluated the covalent/electrostatic character of the H-bonds applying the Cremer-Kraka criterion of covalent bonding. H-bonding in the puric bases has a more covalent character whereas in the pyrimidic bases the H-bond character is more electrostatic. We investigated how the mutual orientation of the CH2OH group and the base influences H-bond formation via two geometrical parameters describing the rotation of the substituents perpendicular to the sugar ring and their tilting relative to the ring center. According to our results, rotation is more important for H-bond formation. In addition we assessed the influence of the H-bond acceptor, the lone pair (N, respectively O), via the delocalization energy. We found larger delocalization energies corresponding to stronger H-bonds for the puric bases. The global minimum conformation of 2'-deoxyguanosine has the strongest H-bond of all conformers investigated in this work with a bond strength of 0.436 which is even stronger than the H-bond in the water dimer (0.360). The application of our new analysis to DNA deoxyribonucleotides and to unnatural base pairs, which have recently drawn a lot of attention, is in progress.