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

Multiple hydrogen-bonded dimers: are only the frontier atoms relevant?

Non-frontier atom exchanges in hydrogen-bonded aromatic dimers can induce significant interaction energy changes (up to 6.5 kcal mol-1). Our quantum-chemical analyses reveal that the relative hydrogen-bond strengths of N-edited guanine-cytosine base pair isosteres, which cannot be explained from the frontier atoms, follow from the charge accumulation in the monomers.

Nearest-neighbour parametrization of DNA single, double and triple mismatches at low sodium concentration

DNA mismatches, that is, base pairs different from the canonical AT and CG, are involved in numerous biological processes and can be a problem for technological applications such as PCR amplification. The nearest-neighbour (NN) model is the standard approach for predicting melting temperatures and is used in methods of secondary structure predictions and modelling of hybridization kinetics. However, despite its biological and technological importance, existing NN parameters that include DNA mismatches are incomplete, and those available were obtained from a limited set of melting temperature at high sodium concentration. To our knowledge, there is currently no NN set of parameters for up to three mismatches covering all configurations at low sodium concentrations. Here, we are applying the NN model to a large set of 4096 published melting temperatures, covering all combinations of single, double and triple mismatches. Dealing with such a large set of temperature is challenging in several ways, bringing new methodological problems. Here, optimizing a large number of 252 independent parameters has required the development of a new method where we readjust the seed parameters using the definition of the Gibbs free energy. The new parameters predict the training set within 1.1 °C and the validation set to 2.7 °C.

Strengthened cooperativity of DNA-based cyclic hydrogen-bonded rosettes by subtle functionalization

Cooperative effects cause extra stabilization of hydrogen-bonded supramolecular systems. In this work we have designed hydrogen-bonded rosettes derived from a guanine-cytosine Janus-type motif with the aim of finding a monomer that enhances the synergy of supramolecular systems. For this, relativistic dispersion-corrected density functional theory computations have been performed. Our proposal involves a monomer with three hydrogen-bonds pointing in the same direction, which translates into shorter bonds, stronger donor-acceptor interactions, and more attractive electrostatic interactions, thus giving rise to rosettes with strengthened cooperativity. This newly designed rosette has triple the cooperativity found for the naturally occurring guanine quadruplex.

Structural Insights into the Antiparallel G-Quadruplex in the Presence of K+ and Mg2+ Ions

G-Quadruplex (GQ) is a secondary structural unit of DNA, formed at the telomere region of the chromosome with a high guanine content. It is reported that the GQs can hinder many biological processes. Thus, research thrives to explore the structural stability of GQs. Studies based on circular dichroism (CD) and nuclear magnetic resonance (NMR) experiments established the vital role of cations such as K⁺ and Mg²⁺ in the stability of antiparallel G-quadruplexes (AGQs). However, there is a need to understand how stability in AGQ is attained in the presence of cations. Here, we employed molecular dynamics (MD) simulations, steered MD (SMD) simulations, and QM/MM calculations to understand the biophysical and electronic bases of the stability imparted to AGQs via cation binding. Our results showed that Mg²⁺ prefers to bind in the plane with the guanine tetrad, whereas K⁺ binds in between the AGQ tetrads. Thus, three Mg²⁺ cations or two K⁺ ions are needed to stabilize an AGQ molecule, where each and every tetrad binds to Mg²⁺ more robustly with a higher binding affinity. SMD revealed that the traversal of K⁺ through the AGQ central channel required less force than that of Mg²⁺, illustrating the presence of more strong interactions between Mg²⁺ and AGQ tetrads compared to K⁺. The stabilization in the AGQ tetrads due to cation binding is reassessed by employing ab initio simulations. Mixed QM/MM calculations confirmed that Mg²⁺ binds strongly to AGQ compared to K⁺, and it induces higher interactions between the guanine tetrads. However, K⁺ binding to AGQ induces a higher stabilization energy than Mg²⁺ binding to AGQ tetrads. Despite the higher binding energy, Mg²⁺ binding imparts lower stabilization to AGQ due to its unfavorable fermionic quantum energy.

Cocrystals; basic concepts, properties and formation strategies

Cocrystallization is an old technique and remains the focus of several research groups working in the field of Chemistry and Pharmacy. This technique is basically in field for improving physicochemical properties of material which can be active pharmaceutical ingredients (APIs) or other chemicals with poor profile. So this review article has been presented in order to combine various concepts for scientists working in the field of chemistry, pharmacy or crystal engineering, also it was attempt to elaborate concepts belonging to crystal designing, their structures and applications. A handsome efforts have been made to bring scientists together working in different fields and to make chemistry easier for a pharmacist and pharmacy for chemists pertaining to cocrystals. Various aspects of chemicals being used as co-formers have been explored which predict the formation of co-crystals or molecular salts and even inorganic cocrystals.

Quantification of secondary electrostatic interactions in H-bonded complexes

The H-bonding properties of compounds that contain multiple functional groups are difficult to predict, because there are through-bond polarisation effects and long-range secondary electrostatic interactions that have significant effects on the interactions with solvents and other molecules. Here we use experimental measurements of association constants for formation of 1 : 1 H-bonded complexes that contain a single well-defined H-bond and a single well-defined secondary electrostatic interaction to quantify the magnitude of this effect. The results were used to develop a computational method for calculating functional group H-bond parameters that accurately reproduce the magnitudes of both primary H-bonding interaction and secondary electrostatic interactions. The effects of secondary electrostatic interactions are observed in calculations of ab initio Molecular Electrostatic Potential (MEP) values, but at the van der Waals surface, the magnitude of the effect is highly overestimated. MEP values calculated on electron density isosurfaces that lie closer to the nuclei provide a more accurate description of the experimental observations. H-bond parameters calculated using this approach successfully account for the properties of arrays of multiple H-bond donor and acceptor groups in different configurations. The results provide insight into the factors that govern the interaction properties of molecules that contain multiple functional groups and provide an accurate method for prediction of solution phase complexation free energies based on gas phase calculations of individual molecules.

How the Chalcogen Atom Size Dictates the Hydrogen‐Bond Donor Capability of Carboxamides, Thioamides, and Selenoamides

The amino groups of thio‐ and selenoamides can act as stronger hydrogen‐bond donors than of carboxamides, despite the lower electronegativity of S and Se. This phenomenon has been experimentally explored, particularly in organocatalysis, but a sound electronic explanation is lacking. Our quantum chemical investigations show that the NH₂ groups in thio‐ and selenoamides are more positively charged than in carboxamides. This originates from the larger electronic density flow from the nitrogen lone pair of the NH₂ group towards the lower‐lying π*_C=S and π*_C=Se orbitals than to the high‐lying π*_C=O orbital. The relative energies of the π* orbitals result from the overlap between the chalcogen np and carbon 2p atomic orbitals, which is set by the carbon‐chalcogen equilibrium distance, a consequence of the Pauli repulsion between the two bonded atoms. Thus, neither the electronegativity nor the often‐suggested polarizability but the steric size of the chalcogen atom determines the amide's hydrogen‐bond donor capability.

Highly-efficient PVDF adsorptive membrane filtration based on chitosan@CNTs-COOH simultaneous removal of anionic and cationic dyes

An adsorptive membrane filtration based on polyvinylidene fluoride (PVDF) with chitosan (CS) and carboxylated carbon nanotubes (CNTs-COOH) is prepared by method of phase conversion, and the PVDF-CS@CNTs-COOH membranes can effectively separate anionic and cationic dye wastewater. Compared to pure PVDF membranes, PVDF-CS@CNTs-COOH increases pure water flux from 36.39 (L·m^-2·h^-1) to 85.25 (L·m^-2·h^-1), an increase of nearly 230%. The membrane exhibits excellent rejection performance in the filtration of six types of dye wastewater. The modified membranes also performed well in terms of rejection of mixed anionic and cationic dyes and also have a high performance in recycling, with a flux of over 94% for both anionic and cationic dyes. In addition, the adsorption curve fitting results showed that the adsorption process was more consistent with the pseudo-second-order adsorption kinetic model and Langmuir mode.

σ-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.

Hydrogen and Halogen Bonding in Homogeneous External Electric Fields: Modulating the Bond Strengths

Recent years have witnessed various fascinating phenomena arising from the interactions of noncovalent bonds with homogeneous external electric fields (EEFs). Here we performed a computational study to interpret the sensitivity of intrinsic bond strengths to EEFs in terms of steric effect and orbital interactions. The block-localized wavefunction (BLW) method, which combines the advantages of both ab initio valence bond (VB) theory and molecular orbital (MO) theory, and the subsequent energy decomposition (BLW-ED) approach were adopted. The sensitivity was monitored and analyzed using the induced energy term, which is the variation in each energy component along the EEF strength. Systems with single or multiple hydrogen (H) or halogen (X) bond(s) were also examined. It was found that the X-bond strength change to EEFs mainly stems from the covalency change, while generally the steric effect rules the response of H-bonds to EEFs. Furthermore, X-bonds are more sensitive to EEFs, with the key difference between H- and X-bonds lying in the charge transfer interaction. Since phenylboronic acid has been experimentally used as a smart linker in EEFs, switchable sensitivity was scrutinized with the example of the phenylboronic acid dimer, which exhibits two conformations with either antiparallel or parallel H-bonds, thereby, opposite or consistent responses to EEFs. Among the studied systems, the quadruple X-bonds in molecular capsules exhibit remarkable sensitivity, with its interaction energy increased by -95.2 kJ mol^-1 at the EEF strength 0.005 a.u.

Cooperative Self-Assembly in Linear Chains Based on Halogen Bonds

Cooperative properties of halogen bonds were investigated with computational experiments based on dispersion-corrected relativistic density functional theory. The bonding mechanism in linear chains of cyanogen halide (X-CN), halocyanoacetylene (X-CC-CN), and 4-halobenzonitrile (X-C6 H4 -CN) were examined for X = H, Cl, Br, and I. Our energy decomposition and Kohn-Sham molecular-orbital analyses revealed the bonding mechanism of the studied systems. Cyanogen halide and halocyanoacetylene chains possess an extra stabilizing effect with increasing chain size, whereas the 4-halobenzonitrile chains do not. This cooperativity can be traced back to charge separation within the σ-electronic system by charge-transfer between the lone-pair orbital of the nitrogen (σHOMO ) on one unit and the acceptor orbital of the C-X (σ*LUMO ) on the adjacent unit. As such, the HOMO-LUMO gap in the σ-system decreases, and the cooperativity increases with chain length revealing the similarity in the bonding mechanisms of hydrogen and halogen bonds.

Nature of Alkali- and Coinage-Metal Bonds versus Hydrogen Bonds

We have quantum chemically studied the structure and nature of alkali- and coinage-metal bonds (M-bonds) versus that of hydrogen bonds between A-M and B- in archetypal [A-M⋅⋅⋅B]- model systems (A, B=F, Cl and M=H, Li, Na, Cu, Ag, Au), using relativistic density functional theory at ZORA-BP86-D3/TZ2P. We find that coinage-metal bonds are stronger than alkali-metal bonds which are stronger than the corresponding hydrogen bonds. Our main purpose is to understand how and why the structure, stability and nature of such bonds are affected if the monovalent central atom H of hydrogen bonds is replaced by an isoelectronic alkali- or coinage-metal atom. To this end, we have analyzed the bonds between A-M and B- using the activation strain model, quantitative Kohn-Sham molecular orbital (MO) theory, energy decomposition analysis (EDA), and Voronoi deformation density (VDD) analysis of the charge distribution.

Origins of Lewis acid acceleration in nickel-catalysed C–H, C–C and C–O bond cleavage

The effects of charge transfer, Pauli repulsion and electrostatics/polarization are identified as dominant factors for Lewis acid accelerations in Ni-catalyzed C–X (X = H, C and O) bond cleavages.

Large Non-planar Conjugated Molecule with Strong Intermolecular Interactions Achieved with Homoleptic Zn(II) Complex of Di(5-quinolylethynyl)-tetraphenylazadipyrromethene

Zinc(II) complexes of tetraphenylazadipyrromethenes are potential non-planar n-type conjugated materials. To tune the properties, we installed 5-quinolylethynyl groups at the pyrrolic positions. Compared to the complex with 1-napthylethynyl, we found evidence for stronger intermolecular interactions in the new complex, including much higher overlap integrals in crystals. X-ray analysis revealed unconventional C-H···N hydrogen bonding between two quinolyls of neighboring molecules, pointing to a new strategy for the development of non-planar molecular semiconductors with stronger intermolecular interactions.

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.

Anatomy of Base Pairing in DNA by Interacting Quantum Atoms

The formation of purine and pyrimidine base pairs (BPs), which contributes to shaping of the canonical and noncanonical 3D structures of nucleic acids, is one the most investigated phenomena in chemistry and life sciences. In this contribution, the anatomy of the bond energy (BDE) of the base-pairing interaction in 39 different arrangements found experimentally or predicted for DNA structures containing the four common nucleobases (A, C, G, T) in their neutral or protonated forms is described in light of the theory of interacting quantum atoms within the context of the quantum theory of atoms in molecules. The interplay of individual energy components involved in the three stages of the bond formation process (structural deformation, electron-density promotion, and intermolecular interaction) is studied. We recognized that for the neutral BPs, variations in the kinetic and electrostatic contributions to the BDE are rather negligible, leaving the exchange-correlation energy as the main stabilizing component. It is shown that the contribution of the exchange-correlation term can be recovered by including atoms that are formally assumed to be hydrogen bonded (primary interaction). In contrast, to recover the electrostatic component of interaction, one must consider both the primary and secondary (formally nonbonded atoms) interatomic interactions. The results of our study were employed to design new types of BPs with altered bonding anatomy. We demonstrate that improving the electrostatic characteristics of the BPs does not necessarily result in greater interaction energies if weak secondary hydrogen bonding is destroyed. However, the main tuning factor for systems with conserved interacting faces (primary interactions) is the electrostatic component of the interaction energy resulting from the secondary atom-atom electrostatics.

Melting temperature measurement and mesoscopic evaluation of single, double and triple DNA mismatches

Unlike the canonical base pairs AT and GC, the molecular properties of mismatches such as hydrogen bonding and stacking interactions are strongly dependent on the identity of the neighbouring base pairs. As a result, due to the sheer number of possible combinations of mismatches and flanking base pairs, only a fraction of these have been studied in varying experiments or theoretical models. Here, we report on the melting temperature measurement and mesoscopic analysis of contiguous DNA mismatches in nearest-neighbours and next-nearest neighbour contexts. A total of 4032 different mismatch combinations, including single, double and triple mismatches were covered. These were compared with 64 sequences containing all combinations of canonical base pairs in the same location under the same conditions. For a substantial number of single mismatch configurations, 15%, the measured melting temperatures were higher than the least stable AT base pair. The mesoscopic calculation, using the Peyrard-Bishop model, was performed on the set of 4096 sequences, and resulted in estimates of on-site and nearest-neighbour interactions that can be correlated to hydrogen bonding and base stacking. Our results confirm many of the known properties of mismatches, including the peculiar sheared stacking of tandem GA mismatches. More intriguingly, it also reveals that a number of mismatches present strong hydrogen bonding when flanked on both sites by other mismatches. To highlight the applicability of our results, we discuss a number of practical situations such as enzyme binding affinities, thymine DNA glycosylase repair activity, and trinucleotide repeat expansions.

Hydrogen bond design principles

Hydrogen bonding principles are at the core of supramolecular design. This overview features a discussion relating molecular structure to hydrogen bond strengths, highlighting the following electronic effects on hydrogen bonding: electronegativity, steric effects, electrostatic effects, π-conjugation, and network cooperativity. Historical developments, along with experimental and computational efforts, leading up to the birth of the hydrogen bond concept, the discovery of nonclassical hydrogen bonds (C-H…O, O-H…π, dihydrogen bonding), and the proposal of hydrogen bond design principles (e.g., secondary electrostatic interactions, resonance-assisted hydrogen bonding, and aromaticity effects) are outlined. Applications of hydrogen bond design principles are presented.

C5-Substituted 2-Selenouridines Ensure Efficient Base Pairing with Guanosine; Consequences for Reading the NNG-3' Synonymous mRNA Codons

5-Substituted 2-selenouridines (R5Se2U) are post-transcriptional modifications present in the first anticodon position of transfer RNA. Their functional role in the regulation of gene expression is elusive. Here, we present efficient syntheses of 5-methylaminomethyl-2-selenouridine (1, mnm5Se2U), 5-carboxymethylaminomethyl-2-selenouridine (2, cmnm5Se2U), and Se2U (3) alongside the crystal structure of the latter nucleoside. By using pH-dependent potentiometric titration, pKa values for the N3H groups of 1–3 were assessed to be significantly lower compared to their 2-thio- and 2-oxo-congeners. At physiological conditions (pH 7.4), Se2-uridines 1 and 2 preferentially adopted the zwitterionic form (ZI, ca. 90%), with the positive charge located at the amino alkyl side chain and the negative charge at the Se2-N3-O4 edge. As shown by density functional theory (DFT) calculations, this ZI form efficiently bound to guanine, forming the so-called “new wobble base pair”, which was accepted by the ribosome architecture. These data suggest that the tRNA anticodons with wobble R5Se2Us may preferentially read the 5′-NNG-3′ synonymous codons, unlike their 2-thio- and 2-oxo-precursors, which preferentially read the 5′-NNA-3′ codons. Thus, the interplay between the levels of U-, S2U- and Se2U-tRNA may have a dominant role in the epitranscriptomic regulation of gene expression via reading of the synonymous 3′-A- and 3′-G-ending codons.

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.

Understanding chemical reactivity using the activation strain model

Understanding chemical reactivity through the use of state-of-the-art computational techniques enables chemists to both predict reactivity and rationally design novel reactions. This protocol aims to provide chemists with the tools to implement a powerful and robust method for analyzing and understanding any chemical reaction using PyFrag 2019. The approach is based on the so-called activation strain model (ASM) of reactivity, which relates the relative energy of a molecular system to the sum of the energies required to distort the reactants into the geometries required to react plus the strength of their mutual interactions. Other available methods analyze only a stationary point on the potential energy surface, but our methodology analyzes the change in energy along a reaction coordinate. The use of this methodology has been proven to be critical to the understanding of reactions, spanning the realms of the inorganic and organic, as well as the supramolecular and biochemical, fields. This protocol provides step-by-step instructions-starting from the optimization of the stationary points and extending through calculation of the potential energy surface and analysis of the trend-decisive energy terms-that can serve as a guide for carrying out the analysis of any given reaction of interest within hours to days, depending on the size of the molecular system.

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.

The staining efficiency of cyanine dyes for single-stranded DNA is enormously dependent on nucleotide composition

The staining of nucleic acids with fluorescent dyes is one of the most fundamental technologies in relevant areas of science. For reliable and quantitative analysis, the staining efficiency of the dyes should not be very dependent on the sequences of the specimens. However, this assumption has not necessarily been confirmed by experimental results, especially in the staining of ssDNA (and RNA). In this study, we found that both SYBR Green II and SYBR Gold did not stain either homopyrimidines or ssDNA composed of only adenine (A) and cytosine (C). However, these two dyes emit strong fluorescence when the ssDNA contains both guanine (G) and C (and/or both A and thymine (T)) and form potential Watson-Crick base pairs. Interestingly, SYBR Gold, but not SYBR Green II, strongly stains ssDNA consisting of G and A (or G and T). Additionally, we found that the secondary structure of ssDNA may play an important role in DNA staining. To obtain reliable results for practical applications, sufficient care must be paid to the composition and sequence of ssDNA.

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.

Designing Self-Assembled Rosettes: Why Ammeline is a Superior Building Block to Melamine

In supramolecular chemistry, the rational design of self-assembled systems remains a challenge. Herein, hydrogen-bonded rosettes of melamine and ammeline have been theoretically examined by using dispersion-corrected density functional theory (DFT-D). Our bonding analyses, based on quantitative Kohn-Sham molecular orbital theory and corresponding energy decomposition analyses (EDA), show that ammeline is a much better building block than melamine for the fabrication of cyclic complexes based on hydrogen bonds. This superior capacity is explained by both stronger hydrogen bonding and the occurrence of a strong synergy.

Why do A·T and G·C self-sort? Hückel aromaticity as a driving force for electronic complementarity in base pairing

Computations reveal that the potential for aromaticity gain and loss in nucleobases play key roles in modulating base pairing strengths.

Factors Controlling the Diels-Alder Reactivity of Hetero-1,3-Butadienes

We have quantum chemically explored the Diels-Alder reactivities of a systematic series of hetero-1,3-butadienes with ethylene by using density functional theory at the BP86/TZ2P level. Activation strain analyses provided physical insight into the factors controlling the relative cycloaddition reactivity of aza- and oxa-1,3-butadienes. We find that dienes with a terminal heteroatom, such as 2-propen-1-imine (NCCC) or acrolein (OCCC), are less reactive than the archetypal 1,3-butadiene (CCCC), primarily owing to weaker orbital interactions between the more electronegative heteroatoms with ethylene. Thus, the addition of a second heteroatom at the other terminal position (NCCN and OCCO) further reduces the reactivity. However, the introduction of a nitrogen atom in the backbone (CNCC) leads to enhanced reactivity, owing to less Pauli repulsion resulting from polarization of the diene HOMO in CNCC towards the nitrogen atom and away from the terminal carbon atom. The Diels-Alder reactions of ethenyl-diazene (NNCC) and 1,3-diaza-butadiene (NCNC), which contain heteroatoms at both the terminal and backbone positions, are much more reactive due to less activation strain compared to CCCC.

NMR 1H-Shielding Constants of Hydrogen-Bond Donor Reflect Manifestation of the Pauli Principle

NMR spectroscopy is one of the most useful methods for detection and characterization of hydrogen bond (H-bond) interactions in biological systems. For H bonds X-H···Y, where X and Y are O or N, it is generally believed that a decrease in ¹H-shielding constants relates to a shortening of H-bond donor-acceptor distance. Here we investigated computationally the trend of ¹H-shielding constants for hydrogen-bonded protons in a series of guanine C8-substituted GC pair model compounds as a function of the molecular structure. Furthermore, the electron density distribution around the hydrogen atom was analyzed with the Voronoi deformation density (VDD) method. Our findings demonstrate that ¹H-shielding values of the hydrogen bond are determined by the depletion of charge around the hydrogen atom, which stems from the fact that electrons obey the Pauli exclusion principle.

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.

The magnetic fingerprint of dithiazolyl-based molecule magnets

The ferromagnetic fingerprint of dithiazolyl-based molecule materials is uncovered. Interestingly geometrical rather than electronic structure factors play the leading role.

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.