Understanding chemical reactivity using the activation strain model

Contains fulltext : 216198.pdf (publisher's version ) (Closed access

The Nature of Hydrogen Bonds: A Delineation of the Role of Different Energy Components on Hydrogen Bond Strengths and Lengths

Hydrogen bonds are a complex interplay between different energy components, and their nature is still subject of an ongoing debate. In this minireview, we therefore provide an overview of the different perspectives on hydrogen bonding. This will be done by discussing the following individual energy components: 1) electrostatic interactions, 2) charge-transfer interactions, 3) π-resonance assistance, 4) steric repulsion, 5) cooperative effects, 6) dispersion interactions and 7) secondary electrostatic interactions. We demonstrate how these energetic factors are essential in a correct description of the hydrogen bond, and discuss several examples of systems whose energetic and geometrical features are not captured by easy-to-use predictive models

Hydrogen-Bond Strength of CC and GG Pairs Determined by Steric Repulsion:Electrostatics and Charge Transfer Overruled

Theoretical and experimental studies have elucidated the bonding mechanism in hydrogen bonds as an electrostatic interaction, which also exhibits considerable stabilization by charge transfer, polarization, and dispersion interactions. Therefore, these components have been used to rationalize the differences in strength of hydrogen-bonded systems. A completely new viewpoint is presented, in which the Pauli (steric) repulsion controls the mechanism of hydrogen bonding. Quantum chemical computations on the mismatched DNA base pairs CC and GG (C=cytosine, G=guanine) show that the enhanced stabilization and shorter distance of GG is determined entirely by the difference in the Pauli repulsion, which is significantly less repulsive for GG than for CC. This is the first time that evidence is presented for the Pauli repulsion as decisive factor in relative hydrogen-bond strengths and lengths

A Quantitative Molecular Orbital Perspective of the Chalcogen Bond

Contains fulltext : 232227.pdf (Publisher’s version ) (Open Access

Secondary Electrostatic Interaction Model Revised: Prediction Comes Mainly from Measuring Charge Accumulation in Hydrogen-Bonded Monomers

The secondary electrostatic interaction (SEI) model is often used to predict and explain relative hydrogen bond strengths of self-assembled systems. The SEI model oversimplifies the hydrogen-bonding mechanisms by viewing them as interacting point charges, but nevertheless experimental binding strengths are often in line with the model's predictions. To understand how this rudimentary model can be predictive, we computationally studied two tautomeric quadruple hydrogen-bonded systems, DDAA-AADD and DADA-ADAD. Our results reveal that when the proton donors D (which are electron-donating) and the proton acceptors A (which are electron-withdrawing) are grouped together as in DDAA, there is a larger accumulation of charge around the frontier atoms than when the proton donor and acceptor groups are alternating as in DADA. This accumulation of charge makes the proton donors more positive and the proton acceptors more negative, which enhances both the electrostatic and covalent interactions in the DDAA dimer. The SEI model is thus predictive because it provides a measure for the charge accumulation in hydrogen-bonded monomers. Our findings can be understood from simple physical organic chemistry principles and provide supramolecular chemists with meaningful understanding for tuning hydrogen bond strengths and thus for controlling the properties of self-assembled systems

Tuning the Binding Strength of even and Uneven Hydrogen-Bonded Arrays with Remote Substituents

We investigated the tunability of hydrogen bond strength by altering the charge accumulation around the frontier atoms with remote substituents. For pyridine···H2O with NH2 and CN substituted at different positions on pyridine, we find that the electron-withdrawing CN group decreases the negative charge accumulation around the frontier atom N, resulting in weakening of the hydrogen bond, whereas the electron-donating NH2 group increases the charge accumulation around N, resulting in strengthening of the hydrogen bond. By applying these design principles on DDAA-AADD, DADA-ADAD, DAA-ADD, and ADA-DAD hydrogen-bonded dimers, we find that the effect of the substituent is delocalized over the whole molecular system. As a consequence, systems with an equal number of hydrogen bond donor (D) and acceptor (A) atoms are not tunable in a predictable way because of cancellation of counteracting strengthening and weakening effects. Furthermore, we show that the position of the substituent and long-range electrostatics can play an important role as well. Overall, the design principles presented in this work are suitable for monomers with an unequal number of donor and acceptor atoms and can be exploited to tune the binding strength of supramolecular building blocks

The Nature of Nonclassical Carbonyl Ligands Explained by Kohn–Sham Molecular Orbital Theory

When carbonyl ligands coordinate to transition metals, their bond distance either increases (classical) or decreases (nonclassical) with respect to the bond length in the isolated CO molecule. C−O expansion can easily be understood by π-back-donation, which results in a population of the CO's π*-antibonding orbital and hence a weakening of its bond. Nonclassical carbonyl ligands are less straightforward to explain, and their nature is still subject of an ongoing debate. In this work, we studied five isoelectronic octahedral complexes, namely Fe(CO)62+, Mn(CO)6+, Cr(CO)6, V(CO)6− and Ti(CO)62−, at the ZORA-BLYP/TZ2P level of theory to explain this nonclassical behavior in the framework of Kohn–Sham molecular orbital theory. We show that there are two competing forces that affect the C−O bond length, namely electrostatic interactions (favoring C−O contraction) and π-back-donation (favoring C−O expansion). It is a balance between those two terms that determines whether the carbonyl is classical or nonclassical. By further decomposing the electrostatic interaction ΔVelstat into four fundamental terms, we are able to rationalize why ΔVelstat gives rise to the nonclassical behavior, leading to new insights into the driving forces behind C−O contraction

Front cover artwork - How Divalent Cations Interact with the Internal Channel Site of Guanine Quadruplexes

The front cover artwork is provided by the TheoCheM group of the Vrije Universiteit Amsterdam. The image shows that guanine quadruplexes have a stronger binding affinity for divalent cations than monovalent cations. Read the full text of the Article at 10.1002/cphc.202100529

How Divalent Cations Interact with the Internal Channel Site of Guanine Quadruplexes

The formation of guanine quadruplexes (GQ) in DNA is crucial in telomere homeostasis and regulation of gene expression. Pollution metals can interfere with these DNA superstructures upon coordination. In this work, we study the affinity of the internal GQ channel site towards alkaline earth metal (Mg2+, Ca2+, Sr2+, and Ba2+), and (post-)transition metal (Zn2+, Cd2+, Hg2+, and Pb2+) cations using density functional theory computations. We find that divalent cations generally bind to the GQ cavity with a higher affinity than conventional monovalent cations (e. g. K+). Importantly, we establish the nature of the cation-GQ interaction and highlight the relationship between ionic and nuclear charge, and the electrostatic and covalent interactions. The covalent interaction strength plays an important role in the cation affinity and can be traced back to the relative stabilization of cations’ unoccupied atomic orbitals. Overall, our findings contribute to a deeper understanding of how pollution metals could induce genomic instability