The Acetylene−Ammonia Dimer as a Prototypical C−H···N Hydrogen-Bonded System: An Assessment of Theoretical Procedures

Molecular Junctions Inspired by Nature: Electrical Conduction through Noncovalent Nanobelts

Charge transport occurs in a range of biomolecular systems, whose structures have covalent and noncovalent bonds. Understanding from these systems have yet to translate into molecular junction devices. We design junctions which have hydrogen-bonds between the edges of a series of prototype noncovalent nanobelts (NCNs) and vary the number of donor-acceptors to study their electrical properties. From frontier molecular orbitals (FMOs) and projected density of state (DOS) calculations, we found these NCN dimer junctions to have low HOMO-LUMO gaps and states at the Fermi level, suggesting these are metallic-like systems. Their conductance properties were studied with nonequilibrium Green's functions density functional theory (NEGF-DFT) and was found to decrease with cooperative H-bonding, that is, the conductance decreased as the alternating donor-acceptors around the nanobelts attenuates to a uniform distribution in the H-bonding arrays. The latter gave the highest conductance of 51.3 × 10^-6 S and the Seebeck coefficient showed n-type (-36 to -39 μV K^-1) behavior, while the lower conductors with alternating H-bonds are p-type (49.7 to 204 μV K^-1). In addition, the NCNs have appreciable binding energies (19.8 to 46.1 kcal mol^-1), implying they could form self-assembled monolayer (SAM) heterojunctions leading to a polymeric network for long-range charge transport.

Dispersion interactions between neighboring Bi atoms in (BiH3 )2 and Te(BiR2 )2

Triggered by the observation of a short Bi⋯Bi distance and a BiTeBi bond angle of only 86.6° in the crystal structure of bis(diethylbismuthanyl)tellurane quantum chemical computations on interactions between neighboring Bi atoms in Te(BiR₂ )₂ molecules (R = H, Me, Et) and in (BiH₃ )₂ were undertaken. Bi⋯Bi distances atoms were found to significantly shorten upon inclusion of the d shells of the heavy metal atoms into the electron correlation treatment, and it was confirmed that interaction energies from spin component-scaled second-order Møller-Plesset theory (SCS-MP2) agree well with coupled-cluster singles and doubles theory including perturbative triples (CCSD(T)). Density functional theory-based symmetry-adapted perturbation theory (DFT-SAPT) was used to study the anisotropy of the interplay of dispersion attraction and steric repulsion between the Bi atoms. Finally, geometries and relative stabilities of syn-syn and syn-anti conformers of Te(BiR₂ )₂ (R = H, Me, Et) and interconversion barriers between them were computed. © 2018 Wiley Periodicals, Inc.

Photofragment Spectroscopy and Predissociation Dynamics of Weakly Bound Molecules

Photofragment spectroscopy is combined with imaging techniques and time-resolved measurements of photoions and photoelectrons to explore the predissociation dynamics of weakly bound molecules. Recent experimental advances include measurements of pair-correlated distributions, in which energy disposal in one cofragment is correlated with a state-selected level of the other fragment, and femtosecond pump-probe experiments, in some cases with coincidence detection. An application in which coincident measurements are carried out in the molecular frame is also described. To illustrate these state-selective and time-resolved techniques, we review two recent applications: (a) the photoinitiated dissociation of the covalently bound NO dimer on the ground and excited electronic states and the role of state couplings and (b) the state-selected vibrational predissociation of hydrogen-bonded acetylene dimers with HCl (acid) and ammonia (base) and the importance of angular momentum constraints. We highlight the crucial role of theoretical models in interpreting results.

Bonding rearrangements of hydrogen-bonded complexes involving alkynes

Molecules containing a C-C triple bond, such as HC[triple bond]CH, FC[triple bond]CF, and the C[triple bond]CH radical, are allowed to interact with a partner molecule of H2O, NH3, or HF. Quantum chemical calculations show that these C[triple bond]CH...X H-bonded complexes are bound by up to 4 kcal x mol(-1). More importantly, they can rearrange in such a way that the partner molecule adds to the triple bond so as to form a double C=C bond. Whereas this process is strongly exoergic, there is a high-energy barrier to this rearrangement process. On the other hand, when a second water molecule is added to the complex, it can shuttle protons from the donor part of the complex to the acceptor, and thereby greatly reduce the rearrangement energy barrier. In the case of CCH + 2H2O, this barrier is computed to be less than 4 kcal x mol(-1).

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.

Imaging the State-Specific Vibrational Predissociation of the C2H2−NH3Hydrogen-Bonded Dimer

The state-to-state vibrational predissociation (VP) dynamics of the hydrogen-bonded ammonia-acetylene dimer were studied following excitation in the asymmetric CH stretch. Velocity map imaging (VMI) and resonance-enhanced multiphoton ionization (REMPI) were used to determine pair-correlated product energy distributions. Following vibrational excitation of the asymmetric CH stretch fundamental, ammonia fragments were detected by 2 + 1 REMPI via the B1E'' <-- X1A1' and C'1A1' <-- X1A1' transitions. The fragments' center-of-mass (c.m.) translational energy distributions were determined from images of selected rotational levels of ammonia with one or two quanta in the symmetric bend (nu2 umbrella mode) and were converted to rotational-state distributions of the acetylene co-fragment. The latter is always generated with one or two quanta of bending excitation. All the distributions could be fit well when using a dimer dissociation energy of D0 = 900 +/- 10 cm(-1). Only channels with maximum translational energy <150 cm(-1) are observed. The rotational excitation in the ammonia fragments is modest and can be fit by temperatures of 150 +/- 50 and 50 +/- 20 K for 1nu2 and 2nu2, respectively. The rotational distributions in the acetylene co-fragment pair-correlated with specific rovibrational states of ammonia appear statistical as well. The vibrational-state distributions, however, show distinct state specificity among channels with low translational energy release. The predominant channel is NH3(1nu2) + C2H2(2nu4 or 1nu4 + 1nu5), where nu4 and nu5 are the trans- and cis-bend vibrations of acetylene, respectively. A second observed channel, with much lower population, is NH3(2nu2) + C2H2(1nu4). No products are generated in which the ammonia is in the vibrational ground state or the asymmetric bend (1nu4) state, nor is acetylene ever generated in the ground vibrational state or with CC stretch excitation. The angular momentum (AM) model of McCaffery and Marsh is used to estimate impact parameters in the internal collisions that give rise to the observed rotational distributions. These calculations show that dissociation takes place from bent geometries, which can also explain the propensity to excite fragment bending levels. The low recoil velocities associated with the observed channels facilitate energy exchange in the exit channel, which results in statistical-like fragment rotational distributions.

Resonance-assisted hydrogen bonds: a critical examination. Structure and stability of the enols of beta-diketones and beta-enaminones

The characteristics of the intramolecular hydrogen bond (IMHB) for a series of 40 different enols of beta-diketones and their nitrogen counterparts have been systematically analyzed at the B3LYP/6-311+G(3df,2p)//B3LYP/6-311+G(d,p) level of theory. In some cases, two tautomers may exist which are interconnected by a hydrogen shift through the IMHB. In tautomer a the HB donor group (YH) is attached to the six-membered ring, while in tautomer b the HB acceptor (X) is the one that is attached to the six-membered ring. We found that changing an O to a N favors the a tautomer when the atom is endo and the contrary when it is exo, while the presence of a double bond favors the a tautomers. As expected, the OH group behaves as a better HB donor than the NH2 group and the C=NH group as a better HB acceptor than the C=O group, although the first effect clearly dominates. Accordingly, the expected IMHB strength follows the [donor, acceptor] trend: [OH, C=NH] > [OH, C=O] > [NH2, C=NH] > [NH2, C=O]. For all those compounds in which the functionality exhibiting the IMHB is unsaturated (I-type), the IMHB is much stronger than in their saturated counterparts (II-type). However, when the systems of the II-type subset, which are saturated, are constrained to have the HB donor and the HB acceptor lying in the same plane and at the same distance as in the corresponding unsaturated analogue, the IMHB is of similar or even larger strength. Hence, we conclude that, at least for this series of unsaturated compounds, the resonance-assisted hydrogen bond effect is not the primary reason behind the strength of their IMHBs, which is simply a consequence of the structure of the sigma-skeleton of the system that keeps the HB donor and the HB acceptor coplanar and closer to each other.

Effects of Hydrogen-Bonding and Stacking Interactions with Amino Acids on the Acidity of Uracil

The effects of hydrogen-bonding interactions with amino acids on the (N1) acidity of uracil are evaluated using (B3LYP) density functional theory. Many different binding arrangements of each amino acid to three uracil binding sites are considered. The effects on the uracil acidity are found to significantly depend upon the nature of the amino acid and the binding orientation, but weakly depend on the binding site. Our results reveal that in some instances small models for the amino acids can be used, while for other amino acids larger models are required to properly describe the binding to uracil. The gas-phase acidity of uracil is found to increase by up to approximately 60 kJ mol(-1) due to discrete hydrogen-bonding interactions. Although (MP2) stacking interactions with aromatic amino acids decrease the acidity of uracil, unexpected increases in the acidity are found when any of the aromatic amino acids, or the backbone, hydrogen bond to uracil. Consideration of enzymatic and aqueous environments leads to decreases in the effects of the amino acids on the acidity of uracil. However, we find that the magnitude of the decrease varies with the nature of the molecule bound, as well as the (gas-phase) binding orientations and strengths, and therefore solvation effects should be considered on a case-by-case basis in future work. Nevertheless, the effects of amino acid interactions within enzymatic environments are as much as approximately 35 kJ mol(-1). The present study has general implications for understanding the nature of active site amino acids in enzymes, such as DNA repair enzymes, that catalyze reactions involving anionic nucleobase intermediates.

Pentachlorocyclopropane/Base Complexes: Matrix Isolation Infrared Spectroscopic and Density Functional Study of C−H- - -N Hydrogen Bonds

Hydrogen-bonded complexes of pentachlorocyclopropane with the bases acetonitrile, ammonia, monomethylamine, and dimethylamine have been isolated and characterized for the first time in argon matrices at 16 K. Coordination of the proton of pentachlorocyclopropane (Pccp) to the electron donor (N) of the base was evidenced by red shifts of the CH stretching mode. These shifts, which range from 22 to 170 cm(-1), increase in the order CH3CN, NH3, (CH3)NH2, and (CH3)2NH. Density functional theory (DFT) calculations at the B3LYP level agree well with experiment and support the formation of 1:1 complexes of Pccp/base. Distinct changes were observed in ring modes as well as CCl and CCl2 modes. The hydrogen bond energy of the complexes varies from 2.95 to 4.22 kcal/mol and is stronger than our previously studied bromocyclopropane-ammonia complex (2.35 kcal/mol, MP2).

The Interaction of Monosubstituted Benzenes with the Stationary Liquid in Gas Liquid Chromatography

In gas liquid chromatography (GLC), the relative retention values log gamma was mainly expressed by van der Waals energy (the sum of the dispersion E(dis) and repulsive E(rep) energies) to the interactions between monosubstituted benzene derivatives and the nonpolar stationary liquid as squalane. The single exception was that of anilines, and it was corrected by the electrostatic energy (E(ES)) due to C-H/pi hydrogen bond. When the stationary liquid changed from the nonpolar to polar, log gamma was estimated by the inductive interaction energy (included in E(ES)) in addition to the sum of E(dis) and E(rep). In the benzene solution, the relative equilibrium values log K/K(o) introduced from the interactions between phenol and substituted benzene derivatives were estimated by E(ES). The E(ES) of COCH(3), CO(2)C(2)H(5) groups is especially originated in the excited dipole moments micro(e). The relative frequency values log nu/nu(o) derived from O-H or O-D stretching vibration of phenol or methanol-D gave the correlation to E(ES) as well as log K/K(o). That of anilines-methanol-D however had been out of a linear relation to E(ES). The cause is concluded that the aniline-methanol-D is making the proton transfer structure from the discussion about the proton affinity (PA) of the base.

N-methylformamide-benzene complex as a prototypical peptide N-H...pi hydrogen-bonded system: density functional theory and MP2 studies

Although the existence of peptide N-H...pi hydrogen bonds recently has been reported in protein structures, little is known about their strength and binding nature and, therefore, the relative importance of the interaction. To shed light on this binding, the N-methylformamide-benzene complex, as a model of the peptide N-H...pi hydrogen bonding, was studied by using density functional theory and Møller-Plesset second-order perturbation (MP2) methods. The geometry of the complex was fully optimized at the B3LYP/6-31G(d,p) and MP2/6-31G(d,p) levels. The optimized interaction distances are about 3.6 and 3.2 A, respectively, at the two levels. The binding energy corrected by basis set superposition error with the MP2/cc-pVTZ method based on the MP2/6-31G geometry is -4.37 kcal/mol, which is as strong as the conventional hydrogen bonding. The calculated results suggest that the peptide N-H...pi hydrogen bonding is of sufficient strength to play an important role in the stabilization of protein structures. These results are helpful to better understand the characteristics and nature of the peptide N-H...pi interaction as well as to modify current force fields to better represent this special interaction.

Strength of the Calpha H..O hydrogen bond of amino acid residues

Although the peptide C(alpha)H group has historically not been thought to form hydrogen bonds within proteins, ab initio quantum calculations show it to be a potent proton donor. Its binding energy to a water molecule lies in the range between 1.9 and 2.5 kcal/mol for nonpolar and polar amino acids; the hydrogen bond (H-bond) involving the charged lysine residue is even stronger than a conventional OH..O interaction. The preferred H-bond lengths are quite uniform, about 3.32 A. Formation of each interaction results in a downfield shift of the bridging hydrogen's chemical shift and a blue shift in the C(alpha)H stretching frequency, potential diagnostics of the presence of such an H-bond within a protein.

Short-strong hydrogen bonds and a low barrier transition state for the proton transfer reaction in RNase A catalysis: a quantum chemical study

There is growing evidence that some enzymes catalyze reactions through the formation of short-strong hydrogen bonds as first suggested by Gerlt and Gassman. Support comes from several experimental and quantum chemical studies that include correlation energies on model systems. In the present study, the process of proton transfer between hydroxyl and imidazole groups, a model of the crucial step in the hydrolysis of RNA by the enzymes of the RNase A family, is investigated at the quantum mechanical level of density functional theory and perturbation theory at the MP2 level. The model focuses on the nature of the formation of a complex between the important residues of the protein and the hydroxyl group of the substrate. We have also investigated different configurations of the ground state that are important in the proton transfer reaction. The nature of bonding between the catalytic unit of the enzyme and the substrate in the model is investigated by Bader's atoms in molecule theory. The contributions of solvation and vibrational energies corresponding to the reactant, the transition state and the product configurations are also evaluated. Furthermore, the effect of protein environment is investigated by considering the catalytic unit surrounded by complete proteins--RNase A and Angiogenin. The results, in general, indicate the formation of a short-strong hydrogen bond and the formation of a low barrier transition state for the proton transfer model of the enzyme.