Energetic vs Synergetic Stability: A Theoretical Study

Can we quantitatively evaluate the mutual impacts of intramolecular metal-ligand bonds the same as intermolecular noncovalent bonds?

In this paper, we have reviewed several equations for calculating the cooperative energy of two chemical bonds between three fragments/species, regardless of whether they are atoms, ions or molecules, and whether the bonds between them are intra- or intermolecular. It is emphasized that two chemical bonds upon cooperation in a new compound change the bond dissociation energy of each other exactly by the same quantitative value, their cooperative energy, regardless of the nature of the bonds or whether one bond is very weak and another one is very strong. However, the final benefit/drawback of weak bonds from this cooperation can be considerably larger than that of strong bonds. The above statements are supported by a computational study on the various types of inter- and intramolecular chemical bonds.

Cooperation/Competition between Halogen Bonds and Hydrogen Bonds in Complexes of 2,6-Diaminopyridines and X-CY3 (X = Cl, Br; Y = H, F)

The DFT calculations have been performed on a series of two-element complexes formed by substituted 2,6-diaminopyridine (R−PDA) and pyridine (R−Pyr) with X−CY3 molecules (where X = Cl, Br and Y = H, F). The primary aim of this study was to examine the intermolecular hydrogen and halogen bonds in the condition of their mutual coexistence. Symmetry/antisymmetry of the interrelation between three individual interactions is addressed. It appears that halogen bonds play the main role in the stabilization of the structures of the selected systems. However, the occurrence of one or two hydrogen bonds was associated with the favourable geometry of the complexes. Moreover, the impact of different substituent groups attached in the para position to the aromatic ring of the 2,6-diaminopyridine and pyridine on the character of the intermolecular hydrogen and halogen bonds was examined. The results indicate that the presence of electron-donating substituents strengthens the bonds. In turn, the presence of electron-withdrawing substituents reduces the strength of halogen bonds. Additionally, when hydrogen and halogen bonds lose their leading role in the complex formation, the nonspecific electrostatic interactions between dipole moments take their place. Analysis was based on geometric, energetic, and topological parameters of the studied systems.

A detailed hydrogen bonding analysis on the compositions of H2SO4/HNO3/H2O ternary systems: A computational study

Hydrogen bonding properties of H₂SO₄/HNO₃/H₂O ternary molecular clusters have been studied by means of structural, energetic, topological, and spectroscopic perspectives. The roles of the hydrogen bonds in the formation of these clusters are considered according to the molecule positions (proton donor or proton acceptor) in the clusters. 33 stable conformers were identified on the potential energy surface. The global minimum one was obtained when HNO₃ donates a proton to H₂SO₄, however, the cooperativity effect contribution was found to be significant for the cluster where HNO₃ acts as a proton donor to an H₂O molecule.

New equation for calculating total interaction energy in one noncyclic ABC triad and new insights into cooperativity of noncovalent bonds

In this work, a new equation consist of A⋅⋅⋅B, B⋅⋅⋅C, A⋅⋅⋅BC, and AB⋅⋅⋅C interactions is proposed for calculating the total interaction energy of noncyclic ABC triads. New equations are also proposed for calculating the changes in values of A⋅⋅⋅B and B⋅⋅⋅C interactions on the formation of triad from the corresponding dyads. The advantages of equations proposed here in comparison with many-body interaction energy approach are discussed. All proposed equations were tested in F₃ MLi⋅⋅⋅NCH⋅⋅⋅HLH and F₃ MLi⋅⋅⋅HLH⋅⋅⋅HCN (M = C, Si; L = Be, Mg) as well as H₃ N⋅⋅⋅XY⋅⋅⋅HF (X, Y = F, Cl, Br) noncyclic A⋅⋅⋅B⋅⋅⋅C triads. The data show that the total cooperativity of triad correlates well with the sum of the changes in values of A⋅⋅⋅B and B⋅⋅⋅C interactions calculated through new equations proposed here. © 2016 Wiley Periodicals, Inc.

Hydrogen and Dihydrogen Bonds in the Reactions of Metal Hydrides

The dihydrogen bond-an interaction between a transition-metal or main-group hydride (M-H) and a protic hydrogen moiety (H-X)-is arguably the most intriguing type of hydrogen bond. It was discovered in the mid-1990s and has been intensively explored since then. Herein, we collate up-to-date experimental and computational studies of the structural, energetic, and spectroscopic parameters and natures of dihydrogen-bonded complexes of the form M-H···H-X, as such species are now known for a wide variety of hydrido compounds. Being a weak interaction, dihydrogen bonding entails the lengthening of the participating bonds as well as their polarization (repolarization) as a result of electron density redistribution. Thus, the formation of a dihydrogen bond allows for the activation of both the MH and XH bonds in one step, facilitating proton transfer and preparing these bonds for further transformations. The implications of dihydrogen bonding in different stoichiometric and catalytic reactions, such as hydrogen exchange, alcoholysis and aminolysis, hydrogen evolution, hydrogenation, and dehydrogenation, are discussed.

Interplay between Beryllium Bonds and Anion-π Interactions in BeR2:C6X6:Y- Complexes (R = H, F and Cl, X = H and F, and Y = Cl and Br)

A theoretical study of the beryllium bonds in BeR₂:C₆X₆ (R = H, F, Cl and X = H and F) has been carried out by means of MP2/aug′-cc-pVDZ computational methods. In addition, the ternary complexes BeR₂:C₆X₆:Y⁻ (Y = Cl and Br) have been analyzed. Geometric, energetic and electronic aspects of the complexes have been taken into account. All the parameters analyzed provide a clear indication of favorable cooperativity in both interactions observed, beryllium bond and aromatic ring:anion interaction.

Reconciling experiment and theory in the use of aryl-extended calix[4]pyrrole receptors for the experimental quantification of chloride-π interactions in solution

In this manuscript we consider from a theoretical point of view the recently reported experimental quantification of anion–π interactions (the attractive force between electron deficient aromatic rings and anions) in solution using aryl extended calix[4]pyrrole receptors as model systems. Experimentally, two series of calix[4]pyrrole receptors functionalized, respectively, with two and four aryl rings at the meso positions, were used to assess the strength of chloride–π interactions in acetonitrile solution. As a result of these studies the contribution of each individual chloride–π interaction was quantified to be very small (<1 kcal/mol). This result is in contrast with the values derived from most theoretical calculations. Herein we report a theoretical study using high-level density functional theory (DFT) calculations that provides a plausible explanation for the observed disagreement between theory and experiment. The study reveals the existence of molecular interactions between solvent molecules and the aromatic walls of the receptors that strongly modulate the chloride–π interaction. In addition, the obtained theoretical results also suggest that the chloride-calix[4]pyrrole complex used as reference to dissect experimentally the contribution of the chloride–π interactions to the total binding energy for both the two and four-wall aryl-extended calix[4]pyrrole model systems is probably not ideal.

Novel CX⋯π halogen bonds in complexes of acetylene and its derivatives of Na and MPH3 (M=Cu, Ag, Au) with XCCF (X=Cl, Br, I)

Ab initio calculations have been carried out for a variety of model systems with a T-shaped CX⋯π motif. The CX⋯π interaction of acetylene with the halogen donor molecule XCCF (X=Cl, Br, I) is invariably found to be weak with the interaction energy less than 11kJ/mol in magnitude. Substitution of the two protons in acetylene with more electron-donating sodium atoms increases the π electron density in the CC bond and leads to a substantial enhancement in its interaction with the halogen donor. The calculated interaction energies increase to as much as 73kJ/mol in the case of C2Na2-ICCF. The interaction of XCCF with a model coinage metal ethynide, H3PMCCMPH3 (M=Cu, Ag, Au), is intermediate between these two extremes, and the interaction energy is related to the nature of coinage metals. The CX⋯π halogen bonds have been analyzed with natural bond orbital, atoms in molecules, and energy decomposition.

The mutual interactions based on amphipathic tetraoxacalix[2]arene[2]triazine: recognition cases of anion and cation investigated by a computational study

This work provides insights into the interaction nature of amphiphilic tetraoxacalix[2]arene[2]triazine with the system (J. Am. Chem. Soc., 2013, 135, 892) as well as the recognition cases of anion and cation.

Aromatic rings in chemical and biological recognition: energetics and structures

This review describes a multidimensional treatment of molecular recognition phenomena involving aromatic rings in chemical and biological systems. It summarizes new results reported since the appearance of an earlier review in 2003 in host-guest chemistry, biological affinity assays and biostructural analysis, data base mining in the Cambridge Structural Database (CSD) and the Protein Data Bank (PDB), and advanced computational studies. Topics addressed are arene-arene, perfluoroarene-arene, S⋅⋅⋅aromatic, cation-π, and anion-π interactions, as well as hydrogen bonding to π systems. The generated knowledge benefits, in particular, structure-based hit-to-lead development and lead optimization both in the pharmaceutical and in the crop protection industry. It equally facilitates the development of new advanced materials and supramolecular systems, and should inspire further utilization of interactions with aromatic rings to control the stereochemical outcome of synthetic transformations.

Can lone pair-π and cation-π interactions coexist? A theoretical study

The interplay between two important noncovalent interactions involving different aromatic rings is studied by means of ab initio calculations (MP2/6-31++G**) computing the non-additivity energies. In this study we demonstrate the existence of cooperativity effects when cation-π and lone pair-π interactions coexist in the same system. These effects are studied theoretically using energetic and geometric features of the complexes. In addition we use Bader’s theory of atoms-in-molecules and Molecular Interaction Potential with polarization (MIPp) partition scheme to characterize the interactions. Experimental evidence for this combination of interactions has been obtained from the Cambridge Structural Database.

Effects of the biological backbone on stacking interactions at DNA-protein interfaces: the interplay between the backbone···π and π···π components

The (gas-phase) MP2/6-31G*(0.25) π···π stacking interactions between the five natural bases and the aromatic amino acids calculated using (truncated) monomers composed of conjugated rings and/or (extended) monomers containing the biological backbone (either the protein backbone or deoxyribose sugar) were previously compared. Although preliminary energetic results indicated that the protein backbone strengthens, while the deoxyribose sugar either strengthens or weakens, the interaction calculated using truncated models, the reasons for these effects were unknown. The present work explains these observations by dissecting the interaction energy of the extended complexes into individual backbone···π and π···π components. Our calculations reveal that the total interaction energy of the extended complex can be predicted as a sum of the backbone···π and π···π components, which indicates that the biological backbone does not significantly affect the ring system through π-polarization. Instead, we find that the backbone can indirectly affect the magnitude of the π···π contribution by changing the relative ring orientations in extended dimers compared with truncated dimers. Furthermore, the strengths of the individual backbone···π contributions are determined to be significant (up to 18 kJ mol(-1)). Therefore, the origin of the energetic change upon model extension is found to result from a balance between an additional (attractive) backbone···π component and differences in the strength of the π···π interaction. In addition, to understand the effects of the biological backbone on the stacking interactions at DNA-protein interfaces in nature, we analyzed the stacking interactions found in select DNA-protein crystal structures, and verified that an additive approach can be used to examine the strength of these interactions in biological complexes. Interestingly, although the presence of attractive backbone···π contacts is qualitatively confirmed using the quantum theory of atoms in molecules (QTAIM), QTAIM electron density analysis is unable to quantitatively predict the additive relationship of these interactions. Most importantly, this work reveals that both the backbone···π and π···π components must be carefully considered to accurately determine the overall stability of DNA-protein assemblies.

A preliminary investigation of the additivity of pi-pi or pi+-pi stacking and T-shaped interactions between natural or damaged DNA nucleobases and histidine

Previous computational studies have examined pi-pi and pi(+)-pi stacking and T-shaped interactions in nucleobase-amino acid dimers, yet it is important to investigate how additional amino acids affect these interactions since simultaneous contacts often appear in nature. Therefore, this paper investigates the geometries and binding strengths of amino acid-nucleobase-amino acid trimers, which are compared to the corresponding nucleobase-amino acid dimer interactions. We concentrate on systems containing the natural nucleobase adenine or its (cationic) damaged counterpart, 3-methyladenine, and the aromatic amino acid histidine, in both the neutral and protonated forms. This choice of molecules provides information about pi-pi and pi(+)-pi stacking and T-shaped interactions in asymmetric, biologically relevant systems. We determined that both stacked and T-shaped interactions, as well as both pi-pi and pi(+)-pi interactions, exhibit geometric additivity. To investigate the energetic additivity in our trimers, the synergy (E(syn)) and the additivity (E(add)) energy were examined. E(add) reveals that it is important to consider the interaction between the two amino acids when examining the additivity of nucleobase-amino acid interactions. Additionally, E(syn) and E(add) indicate that pi(+)-pi interactions are quite different from pi-pi interactions. The magnitude of E(add) is generally less than 2 kJ mol(-1), which suggests that these interactions are additive. However, the interaction energy analysis does not provide information about the individual interactions in the trimers. Therefore, the quantum theory of atoms in molecules (QTAIM) was implemented. We find inconsistent conclusions from our QTAIM analysis and interaction energy evaluation. However, the magnitudes of the differences between the dimer and trimer critical point properties are extremely small and therefore may not be able to yield conclusive descriptions of differences (if any) between the dimer and trimer interactions. We hypothesize that, due to the limited number of investigations of this type, it is currently unclear how QTAIM can improve our understanding of pi-pi and pi(+)-pi dimers and trimers. Therefore, future work must systematically alter the pi-pi or pi(+)-pi system to definitively determine how the geometry, symmetry, and system size alter the QTAIM analysis, which can then be used to understand biologically relevant complexes.