σ-Electrons Responsible for Cooperativity and Ring Equalization in Hydrogen-Bonded Supramolecular Polymers

Invited for this month's cover are collaborators from the TheoCheM group of the Vrije Universiteit Amsterdam and the University of Perugia. The cover picture shows a σ-electron traveling through a hydrogen-bonded squaramide linear chain. The charge transfer within the σ-electronic system is the cause for the cooperativity in the investigated urea, deltamide, and squaramide polymers. More information can be found in the Full Paper by Célia Fonseca Guerra, and co-workers.

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.

σ-Electrons Responsible for Cooperativity and Ring Equalization in Hydrogen-Bonded Supramolecular Polymers

We have quantum chemically analyzed the cooperative effects and structural deformations of hydrogen-bonded urea, deltamide, and squaramide linear chains using dispersion-corrected density functional theory at BLYP-D3(BJ)/TZ2P level of theory. Our purpose is twofold: (i) reveal the bonding mechanism of the studied systems that lead to their self-assembly in linear chains; and (ii) rationalize the C-C bond equalization in the ring moieties of deltamide and squaramide upon polymerization. Our energy decomposition and Kohn-Sham molecular orbital analyses reveal cooperativity in all studied systems, stemming from the charge separation within the σ-electronic system by charge transfer from the carbonyl oxygen lone pair donor orbital of one monomer towards the σ* N-H antibonding acceptor orbital of the neighboring monomer. This key orbital interaction causes the C=O bonds to elongate, which, in turn, results in the contraction of the adjacent C-C single bonds that, ultimately, makes the ring moieties of deltamide and squaramide to become more regular. Notably, the π-electron delocalization plays a much smaller role in the total interaction between the monomers in the chain.

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 (Mg²⁺, Ca²⁺, Sr²⁺, and Ba²⁺), and (post‐)transition metal (Zn²⁺, Cd²⁺, Hg²⁺, and Pb²⁺) 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.

A Quantitative Molecular Orbital Perspective of the Chalcogen Bond

Invited for this month's cover are the groups of Prof. Dr. Teodorico C. Ramalho (Federal University of Lavras and University Hradec Kralove) and Prof. Dr. F. Matthias Bickelhaupt (Vrije Universiteit Amsterdam and Radboud University). The cover picture shows the key message of their work, that is, the covalency of the chalcogen bonds, in an elegantly simple and attractive manner. To that end, the chalcogen bonds are represented by schematic 3D structures of the bond donor D₂ Ch and the bond acceptor A^- , and their attractive interaction in green. Then, a colorful molecular orbital (MO) diagram where the HOMO-LUMO mixing is represented by the mixing of red (HOMO) and blue (LUMO) into purple (MO) is presented. Read the full text of their Full Paper at 10.1002/open.202000323.

A Quantitative Molecular Orbital Perspective of the Chalcogen Bond

We have quantum chemically analyzed the structure and stability of archetypal chalcogen-bonded model complexes D2 Ch⋅⋅⋅A- (Ch = O, S, Se, Te; D, A = F, Cl, Br) using relativistic density functional theory at ZORA-M06/QZ4P. Our purpose is twofold: (i) to compute accurate trends in chalcogen-bond strength based on a set of consistent data; and (ii) to rationalize these trends in terms of detailed analyses of the bonding mechanism based on quantitative Kohn-Sham molecular orbital (KS-MO) theory in combination with a canonical energy decomposition analysis (EDA). At odds with the commonly accepted view of chalcogen bonding as a predominantly electrostatic phenomenon, we find that chalcogen bonds, just as hydrogen and halogen bonds, have a significant covalent character stemming from strong HOMO-LUMO interactions. Besides providing significantly to the bond strength, these orbital interactions are also manifested by the structural distortions they induce as well as the associated charge transfer from A- to D2 Ch.

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)₆²⁺, Mn(CO)₆⁺, Cr(CO)₆, V(CO)₆⁻ and Ti(CO)₆²⁻, 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 ΔV_elstat into four fundamental terms, we are able to rationalize why ΔV_elstat gives rise to the nonclassical behavior, leading to new insights into the driving forces behind C−O contraction.

Halogen Bonds in Ligand-Protein Systems: Molecular Orbital Theory for Drug Design

Halogen bonds are highly important in medicinal chemistry as halogenation of drugs, generally, improves both selectivity and efficacy toward protein active sites. However, accurate modeling of halogen bond interactions remains a challenge, since a thorough theoretical investigation of the bonding mechanism, focusing on the realistic complexity of drug–receptor systems, is lacking. Our systematic quantum-chemical study on ligand/peptide-like systems reveals that halogen bonding is driven by the same bonding interactions as hydrogen bonding. Besides the electrostatic and the dispersion interactions, our bonding analyses, based on quantitative Kohn–Sham molecular orbital theory together with energy decomposition analysis, reveal that donor–acceptor interactions and steric repulsion between the occupied orbitals of the halogenated ligand and the protein need to be considered more carefully within the drug design process.

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.

Nature of Intramolecular Resonance Assisted Hydrogen Bonding in Malonaldehyde and Its Saturated Analogue

The nature of resonance-assisted hydrogen bonds (RAHB) is still subject of an ongoing debate. We therefore analyzed the σ and π charge redistributions associated with the formation of intramolecular hydrogen bonds in malonaldehyde (MA) and its saturated analogue 3-hydroxypropanal (3-OH) and addressed the question whether there is a resonance assistance phenomenon in the sense of a synergistic interplay between the σ and π electron systems. Our quantum chemical calculations at the BP86/TZ2P level of theory show that the π charge flow is indeed in line with the Lewis structure as proposed by the RAHB model. This typical rearrangement of charge is only present in the unsaturated system, and not in its saturated analogue. Resonance in the π electron system assists the intramolecular hydrogen bond by reducing the hydrogen bond distance, and by providing an additional stabilizing component to the net bonding energy. The σ orbital interaction plays an important role in the enhanced hydrogen bond strength in MA as well. However, there is no resonance assistance in the sense of an interplay between σ charge transfer and π polarization; σ and π contribute independently from each other.

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.