Natural bond critical point analysis: Quantitative relationships between natural bond orbital-based and QTAIM-based topological descriptors of chemical bonding

A comprehensive electron wavefunction analysis toolbox for chemists, Multiwfn

Analysis of electron wavefunction is a key component of quantum chemistry investigations and is indispensable for the practical research of many chemical problems. After more than ten years of active development, the wavefunction analysis program Multiwfn has accumulated very rich functions, and its application scope has covered numerous aspects of theoretical chemical research, including charge distribution, chemical bond, electron localization and delocalization, aromaticity, intramolecular and intermolecular interactions, electronic excitation, and response property. This article systematically introduces the features and functions of the latest version of Multiwfn and provides many representative examples. Through this article, readers will be able to fully understand the characteristics and recognize the unique value of Multiwfn. The source code and precompiled executable files of Multiwfn, as well as the manual containing a detailed introduction to theoretical backgrounds and very rich tutorials, can all be downloaded for free from the Multiwfn website (http://sobereva.com/multiwfn).

Natural resonance-theoretic conceptions of extreme electronic delocalization in soft materials

In the broad context of Dalton's atomic hypothesis and subsequent classical vs. quantum understanding of macroscopic materials, we show how Pauling's resonance-type conceptions, as quantified in natural resonance theory (NRT) analysis of modern wavefunctions, can be modified to unify description of interatomic interactions from the Lewis-like limit of localized e-pair covalency in molecules to the extreme delocalized limit of supramolecular "soft matter" aggregation. Such "NRT-centric" integration of NRT bond orders for hard- and soft-matter interactions is illustrated with application to a long-predicted and recently synthesized organometallic sandwich-type complex ("diberyllocene") that exhibits bond orders ranging from the soft limit (bBeC ≈ 0.01) to the typical values (bCC ≈ 1.35) of molecular resonance-covalency in the organic domain, with intermediate value (bBeBe ≈ 0.86) for intermetallic Be⋯Be interaction.

Decoding Chemical Bonds: Assessment of the Basis Set Effect on Overlap Electron Density Descriptors and Topological Properties in Comparison to QTAIM

Quantum chemical bonding descriptors based on the total and overlap density can provide valuable information about chemical interactions in different systems. However, these descriptors can be sensitive to the basis set used. To address this, different numerical treatments of electron density have been proposed to reduce the basis set dependency. In this work, we introduce overlap properties (OPs) obtained through numerical treatment of the electron density and present the topology of overlap density (TOP) for the first time. We compare the basis set dependency of numerical OP and TOP descriptors with their quantum theory of atoms in molecules (QTAIM) counterparts, considering the total electron density. Three single (C-C, C-O, and C-F) bonds in ethane, methanol, and fluoromethane and two double (C═C and C═O) bonds in ethene and formaldehyde were analyzed. Diatomic molecules Li-X with X = F, Cl, and Br were also analyzed. Eight parameters, including QTAIM descriptors and OP/TOP descriptors, are used to assess the basis dependency at the ωB97X-D level of theory using 28 basis sets from three classes: Pople, Ahlrichs, and Dunning. The study revealed that the topological overlap electron density properties exhibit comparatively lesser dependence on the basis set compared to their total electron density counterparts. Remarkably, these properties retain their chemical significance even with reduced basis set dependency. Similarly, numerical OP descriptors show less basis set dependency than their QTAIM counterparts. The excess of polarization functions increases charge concentration in the interatomic region and influences both QTAIM and OP descriptors. The basis sets Def2TZVP, 6-31++G(d,p), 6-311++G(d,p), cc-pVDZ, cc-pVTZ, and cc-pVQZ demonstrate reduced variability for the tested bond classes in this study, with particular emphasis on the triple-ζ quality Ahlrichs' basis set. We recommend against using basis sets with numerous polarization functions, such as augmented Dunning's and Ahlrichs' quadruple-ζ.

Computational Evaluation of N-Based Transannular Interactions in Some Model Fused Medium-Sized Heterocyclic Systems and Implications for Drug Design

As part of a project on fused medium-sized ring systems as potential drugs, we have previously demonstrated the usefulness of Density Functional Theory (DFT) to evaluate amine nitrogen-based transannular interactions across the central 10-membered ring in the bioactive dibenzazecine alkaloid, protopine. A range of related hypothetical systems have been investigated, together with transannular interactions involving ring-embedded imino or azo group nitrogens and atoms or groups (Y) across the ring. Electrostatic potential energies mapped onto electron density surfaces in the different ring conformations were evaluated in order to characterise these conformations. Unexpectedly, the presence of sp² hybridised nitrogen atoms in the medium-sized rings did not influence the conformations appreciably. The strength and type of the N^…Y interactions are determined primarily by the nature of Y. This is also the case when the substituent on the interacting nitrogen is varied from CH₃ (protopine) to H or OH. With Y = BOH, very strong interactions were observed in protopine analogues, as well as in rings incorporating imino or azo groups. Strong to moderate interactions were observed with Y = CS, CO and SO in all ring systems. Weaker interactions were observed with Y = S, O and weaker ones again with an sp³ hybridised carbon (Y = CH₂). The transannular interactions can influence conformational preferencing and shape and change electron distributions at key sites, which theoretically could modify properties of the molecules while providing new or enhanced sites for biological target interactions, such as the H or OH substituent. The prediction of new strong transannular interaction types such as with Y = BOH and CS should be helpful in informing priorities for synthesis and other experimental studies.

A glutamine-based single α-helix scaffold to target globular proteins

The binding of intrinsically disordered proteins to globular ones can require the folding of motifs into α-helices. These interactions offer opportunities for therapeutic intervention but their modulation with small molecules is challenging because they bury large surfaces. Linear peptides that display the residues that are key for binding can be targeted to globular proteins when they form stable helices, which in most cases requires their chemical modification. Here we present rules to design peptides that fold into single α-helices by instead concatenating glutamine side chain to main chain hydrogen bonds recently discovered in polyglutamine helices. The resulting peptides are uncharged, contain only natural amino acids, and their sequences can be optimized to interact with specific targets. Our results provide design rules to obtain single α-helices for a wide range of applications in protein engineering and drug design.

Nitronium salts as mild and inexpensive oxidizing reagents toward designing efficient strategies in organic syntheses; A mechanistic investigation based on the DFT insights

Today, introducing and evaluating the performance of novel reagents are an undeniable part of designing a successful synthetic strategy. Herein, we study the efficiency and mechanism of recently synthesized nitronium salts (e.g., NO₂FSO₃, NO₂CF₃SO₃, NO₂HS₂O₇, NO₂BF₄, NO₂PF₆, and NO₂HSO₄) in the oxidation reaction of ethanol to acetic acid, as a model of the primary alcohol transformations to linear carboxylic acid. An aldehyde molecule is the first produced species in this reaction which is converted to the acetic acid molecule in the presence of in situ-produced nitric acid. Concerning the proposed mechanism, among the studied nitronium salts, two different behaviors can be observed in the transition state of the step in which the aldehyde molecule is formed. The calculated barrier energies of this step have been scrutinized by powerful descriptors such as Quantum Theory of Atoms in Molecules (QTAIM), Natural Bond Orbital (NBO), Electrostatic Potential (ESP) surfaces, and Activation Strain Model (ASM). The outcomes of the studied descriptors illustrate that nitronium salts have different performances in progressing the formation of the aldehyde molecule. Indeed, the likeness of the transition state of this step to the products for NO₂FSO₃, NO₂CF₃SO₃, and NO₂HS₂O₇ species is more significant than the others. Accordingly, these reagents have more potential to apply as oxidizing agents in the primary alcohol transformations to linear carboxylic acid.

Switchable Cycloadditions of Mesoionic Dipoles: Refreshing up a Regioselective Approach to Two Distinctive Heterocycles

Mesoionic rings are among the most versatile 1,3-dipoles, as witnessed recently by their incorporation into bio-orthogonal strategies, and capable of affording unconventional heterocycles beyond the expected scope of Huisgen cycloadditions. Herein, we revisit in detail the reactivity of thiazol-3-ium-4-olates with alkynes, leading to thiophene and/or pyrid-2-one derivatives. A structural variation at the parent mesoionic dipole alters sufficiently the steric outcome, thereby favoring the regioselective formation of a single transient cycloadduct, which undergoes chemoselective fragmentation to either five- or six-membered heterocycles. The synthetic protocol benefits largely from microwave (MW) activation, which enhances reaction rates. The mechanism has been interrogated with the aid of density functional theory (DFT) calculations, which sheds light into the origin of the regioselectivity and points to a predictive formulation of reactivity involving competing pathways of mesoionic cycloadditions.

High-Density “Windowpane” Coordination Patterns of Water Clusters and Their NBO/NRT Characterization

Cluster mixture models for liquid water at higher pressures suggest the need for water clusters of higher coordination and density than those commonly based on tetrahedral H-bonding motifs. We show here how proton-ordered water clusters of increased coordination and density can assemble from a starting cyclic tetramer or twisted bicyclic (Möbius-like) heptamer to form extended Aufbau sequences of stable two-, three-, and four-coordinate “windowpane” motifs. Such windowpane clusters exhibit sharply reduced (~90°) bond angles that differ appreciably from the tetrahedral angles of idealized crystalline ice I_h. Computed free energy and natural resonance theory (NRT) bond orders provide quantitative descriptors for the relative stabilities of clusters and strengths of individual coordinative linkages. The unity and consistency of NRT description is demonstrated to extend from familiar supra-integer bonds of the molecular regime to the near-zero bond orders of the weakest linkages in the present H-bond clusters. Our results serve to confirm that H-bonding exemplifies resonance–covalent (fractional) bonding in the sub-integer range and to further discount the dichotomous conceptions of “electrostatics” for intermolecular bonding vs. “covalency” for intramolecular bonding that still pervade much of freshman-level pedagogy and force-field methodology.

Chain-Length- and Concentration-Dependent Isomerization of Bithiophenyl-Based Diaminotriazine Derivatives in Two-Dimensional Polymorphic Self-Assembly§

Bithiophenyl-based diaminotriazine derivatives (2TDT-n, n = 10, 12, 16, and 18) with different chain lengths display colhex/p6mm mesophases. Their supramolecular self-assembled mechanism is investigated using scanning tunneling microscopy (STM) at the 1-octanoic acid/graphite interface at various concentrations. The chain length effect on the two-dimensional adlayers is observed in this system, and 2TDT-n molecules show a structural phase transition from the four-leaf arrangement to the two-row linear nanostructure accompanied by the emergence of molecular isomerization with the increase of the side-chain length. The self-assembled structure of 2TDT-10 is composed of a four-leaf pattern with uniform s-cis conformers. In 2TDT-12, three kinds of nanostructures (bamboo-like, two-row linear pattern-I, and flower-like) are observed. These nanostructures are randomly constituted by cis and trans conformers, and the ratios of the s-cis conformer in three kinds of patterns are 55.7, 42.3, and 62.5%, respectively. Furthermore, when n = 16 and 18, the ratio of the s-cis conformer further decreases to 19.0 and 4.3%, respectively. Those molecules mainly form linear nanostructures consisting of s-trans conformers. Therefore, it is reasonable to conclude that the side-chain length has a great effect on the self-assembled patterns and the molecular conformation of bithiophenyl-based diaminotriazine derivatives. Density functional theory calculations are applied to optimize molecular conformers and assess their single-point energies, showing that the s-cis conformation has higher energy than the s-trans conformer. We speculate that the ratio of two conformers in nanostructures might be similar to that of the liquid crystalline phase.

Analysis of Conformational Preferences in Caffeine

High level DLPNO−CCSD(T) electronic structure calculations with extended basis sets over B3LYP−D3 optimized geometries indicate that the three methyl groups in caffeine overcome steric hindrance to adopt uncommon conformations, each one placing a C−H bond on the same plane of the aromatic system, leading to the C−H bonds eclipsing one carbonyl group, one heavily delocalized C−N bond constituent of the fused double ring aromatic system, and one C−H bond from the imidazole ring. Deletion of indiscriminate and selective non-Lewis orbitals unequivocally show that hyperconjugation in the form of a bidirectional −CH3 ⇆ aromatic system charge transfer is responsible for these puzzling conformations. The structural preferences in caffeine are exclusively determined by orbital interactions, ruling out electrostatics, induction, bond critical points, and density redistribution because the steric effect, the allylic effect, the Quantum Theory of Atoms in Molecules (QTAIM), and the non-covalent interactions (NCI), all predict wrong energetic orderings. Tiny rotational barriers, not exceeding 1.3 kcal/mol suggest that at room conditions, each methyl group either acts as a free rotor or adopts fluxional behavior, thus preventing accurate determination of their conformations. In this context, our results supersede current experimental ambiguity in the assignation of methyl conformation in caffeine and, more generally, in methylated xanthines and their derivatives.

Ring Vibrations to Sense Anionic Ibuprofen in Aqueous Solution as Revealed by Resonance Raman

We unravel the potentialities of resonance Raman spectroscopy to detect ibuprofen in diluted aqueous solutions. In particular, we exploit a fully polarizable quantum mechanics/molecular mechanics (QM/MM) methodology based on fluctuating charges coupled to molecular dynamics (MD) in order to take into account the dynamical aspects of the solvation phenomenon. Our findings, which are discussed in light of a natural bond orbital (NBO) analysis, reveal that a selective enhancement of the Raman signal due to the normal mode associated with the C-C stretching in the ring, νC=C, can be achieved by properly tuning the incident wavelength, thus facilitating the recognition of ibuprofen in water samples.

A Comprehensive Picture of the Structures, Energies, and Bonding in the Alanine Dimers

High level quantum mechanical computations and extensive stochastic searches of the potential energy surfaces of the Alanine dimers uncover rich and complex structural and interaction landscapes. A total of 416 strongly bound (up 13.4 kcal mol^-1 binding energies at the DLPNO-CCSD(T)/6-311++G(d,p) level corrected by the basis set superposition error and by the zero point vibrational energies over B3LYP-D3 geometries), close energy equilibrium structures were located, bonded via 32 specific types of intermolecular contacts including Y⋅⋅⋅H-X primary and Y⋅⋅⋅H-C secondary hydrogen bonds, H⋅⋅⋅H dihydrogen contacts, and non conventional anti-electrostatic Y δ - ⋯ X δ - interactions. The putative global minimum is triply degenerate, corresponding to the structure of the common dimer of a carboxylic acid. All quantum descriptors of chemical bonding point to a multitude of weak individual interactions within each dimer, whose cumulative effect results in large binding energies and in an attractive fluxional wall of non-covalent interactions in the interstitial region between the monomers.

Thermodynamics and Intermolecular Interactions during the Insertion of Anionic Naproxen into Model Cell Membranes

The insertion process of Naproxen into model dimyristoylphosphatidylcholine (DMPC) membranes is studied by resorting to state-of-the-art classical and quantum mechanical atomistic computational approaches. Molecular dynamics simulations indicate that anionic Naproxen finds an equilibrium position right at the polar/nonpolar interphase when the process takes place in aqueous environments. With respect to the reference aqueous phase, the insertion process faces a small energy barrier of ≈5 kJ mol^-1 and yields a net stabilization of also ≈5 kJ mol^-1. Entropy changes along the insertion path, mainly due to a growing number of realizable microstates because of structural reorganization, are the main factors driving the insertion. An attractive fluxional wall of noncovalent interactions is characterized by all-quantum descriptors of chemical bonding (natural bond orbitals, quantum theory of atoms in molecules, noncovalent interaction, density differences, and natural charges). This attractive wall originates in the accumulation of tiny transfers of electron densities to the interstitial region between the fragments from a multitude of individual intermolecular contacts stabilizing the tertiary drug/water/membrane system.

The CH···HC interaction in biphenyl is a delocalized, molecular-wide and entirely non-classical interaction: Results from FALDI analysis

In this study we aim to determine the origin of the electron density describing a CH···HC interaction in planar and twisted conformers of biphenyl. In order to achieve this, the fragment, atomic, localized, delocalized, intra- and inter-atomic (FALDI) decomposition scheme was utilized to decompose the density in the inter-nuclear region between the ortho-hydrogens in both conformers. Importantly, the structural integrity, hence also topological properties, were fully preserved as no 'artificial' partitioning of molecules was implemented. FALDI-based qualitative and quantitative analysis revealed that the majority of electron density arises from two, non-classical and non-local effects: strong overlap of ortho CH σ-bonds, and long-range electron delocalization between the phenyl rings and ortho carbons and hydrogens. These effects resulted in a delocalized electron channel, that is, a density bridge or a bond path in a QTAIM terminology, linking the H-atoms in the planar conformer. The same effects and phenomena are present in both conformers of biphenyl. We show that the CH···HC interaction is a molecular-wide event due to large and long-range electron delocalization, and caution against approaches that investigate CH···HC interactions without fully taking into account the remainder of the molecule.

Amphipathic 1,3-oxazolidines from N-alkyl glucamines and benzaldehydes: stereochemical and mechanistic studies

Chiral N,O-heterocycles appended to a non-reducing carbohydrate chain, which are valuable synthons, reveal further stereodynamic implications.

Comparison of Chemical and Interpretative Methods: the Carbon-Boron π-Bond as a Test Case*

Quantum chemical calculations and NBO, ETS-NOCV, QTAIM and ELF interpretative approaches have been carried out on C-donor ligand-stabilized dihydrido borenium cations. Numerous descriptors of the C-B π-bond strength obtained from orbital localization, energy partitioning or topological methods as well as from structural and chemical parameters have been calculated for 39 C-donor ligands including N-heterocyclic carbenes and carbones. Comparison of the results allows the identification of relative and absolute descriptors of the π interaction. For both families of descriptors excellent correlations are obtained. This enables the establishment of a π-donation capability scale and shows that the interpretative methods, despite their conceptual differences, describe the same chemical properties. These results also reveal noticeable shortcomings in these popular methods, and some precautions that need to be taken to interpret their results adequately.

Double-Ring Epimerization in the Biosynthesis of Clavulanic Acid

All reaction steps during the biosynthesis of suicidal clavulanic acid (coformulated with β-lactam antibiotics and used to fight bacterial infections) are known, except for the crucial 3S,5S → 3R,5R double epimerization needed to produce a biologically active stereoisomer, for which mechanistic hypothesis is subject to debate. In this work, we provide evidence for a reaction channel for the double inversion of configuration that involves a total of six reaction steps. When mediated by an enzyme with a terminal S-H bond, this highly complex reaction is spontaneous in the absence of solvents. Polarizable continuum models introduce reaction barriers in aqueous environments because of the strong destabilization of the first transition state. Molecular geometries and electronic structures in both cases indicate that solvent-free spontaneity and aqueous medium barriers are both firmly rooted in a substantial reorganization of the electron density right at the onset of the reaction, mostly involving a cyclic evolution/involution of large regions of π delocalization used to stabilize the excess charge left after the initial proton abstraction.

Block deformation analysis: Density matrix blocks as intramolecular deformation density

Block deformation analysis as deformation density of atomic orbitals is introduced to analyze intramolecular interactions. In this respect, density matrix blocks in terms of natural atomic orbitals are employed to find interacting and noninteracting multicenter subsystem and extract the corresponding deformation density. Eigenanalysis of this deformation density is performed to result eigenvalues and eigenorbitals as displaced charge due to the intramolecular interaction and orbital space responsible for charge reorganization, respectively, that possesses advantages of other methods, simultaneously. It is applied to several small molecules, different types of carbon allotropes including zero-, one-, and two-dimensional nanostructures, and challenging systems such as ortho-hydrogen atoms in planar biphenyl. Results highly correlate with delocalization and Wiberg bond indices and show that eigenvalues of block deformation analysis deserved to be considered as bonding index.

Molecular Orbitals Support Energy-Stabilizing "Bonding" Nature of Bader's Bond Paths

Our MO-based findings proved a bonding nature of each density bridge (DB, or a bond path with an associated critical point, CP) on a Bader molecular graph. A DB pinpoints universal physical and net energy-lowering processes that might, but do not have to, lead to a chemical bond formation. Physical processes leading to electron density (ED) concentration in internuclear regions of three distinctively different homopolar H,H atom-pairs as well as classical C-C and C-H covalent bonds were found to be exactly the same. Notably, properties of individual MOs are internuclear-region specific as they (i) concentrate, deplete, or do not contribute to ED at a CP and (ii) delocalize electron-pairs through either in- (positive) or out-of-phase (negative) interference. Importantly, dominance of a net ED concentration and positive e^--pairs delocalization made by a number of σ-bonding MOs is a common feature at a CP. This feature was found for the covalently bonded atoms as well as homopolar H,H atom-pairs investigated. The latter refer to a DB-free H,H atom-pair of the bay in the twisted biphenyl (Bph) and DB-linked H,H atom-pairs (i) in cubic Li₄H₄, where each H atom is involved in three highly repulsive interactions (over +80 kcal/mol), and (ii) in a weak attractive interaction when sterically clashing in the planar Bph.

A detailed look at the bonding interactions in the microsolvation of monoatomic cations

Global and local descriptors of the properties of intermolecular bonding, formally derived from independent methodologies (QTAIM, NCI, NBO, density differences) afford a highly complex picture of the bonding interactions responsible for microsolvation of monoatomic cations. In all cases, the dominant factor dictating geometries and interaction strengths is the electrophilic power of the metal cation. The formal charge disrupts the hydrogen bonding network otherwise present in pristine water clusters, making the hydrogen bonds considerably stronger, even inducing some degree of covalency. All MO interactions are highly ionic, with strengths than in some cases approach that of the reference LiCl bond. Accumulation of electron density in the region connecting MO is observed, thus, ionic bonding in the microsolvation of monoatomic cations is not as simple as an electrostatic interaction between opposing charges.

A Detailed Classification of Three-Centre Two-Electron Bonds

We evaluate the three-centre two-electron (3c-2e) bonds using atoms in molecules (AIM) and natural bond orbital (NBO) theoretical analyses. They have been classified as ‘open (V)’ or ‘closed (Δ)’, depending on how the three centres were bonded. Herein, we show that they could be classified as V, L, Δ, Y, T and I (linear) arrangements depending on the way the three centres are bonded. These different structures are found in B2H6 (V), CH5+ (V), Me-C2H2+ (L), B3+ (Δ), C3H3+ (Δ), H3+ (Y), 2-norbornyl+ (T), SiH5+ (T), and Al2H7− (I). Our results suggest that CH3Li2+ does not contain a 3c-2e bond according to NBO analysis. Therefore, we propose that 3c-2e bonds are classified more accurately as V, L, Δ, Y, T, or I, based on the electron density topology.

Reactions of the Lipid Hydroperoxides With Aminic Antioxidants: The Influence of Stereoelectronic and Resonance Effects on Hydrogen Atom Transfer

Aminic radical-trapping antioxidants (RTAs), as one of the most important antioxidants, have not received sufficient attention yet. But, an increasing number of aminic RTAs have been identified as ferroptosis inhibitors in recent years, which can potentially mediate many pathological states including inflammation, cancer, neurodegenerative disease, as well as ocular and kidney degeneration. This highlights the importance of aminic RTAs in the field of medicine. Herein, we systematically explored the radical scavenging mechanism of aminic RTAs with a quantum chemical method, particularly emphasizing the role of stereoelectronic factors and resonance factors on the transfer of H-atom and the stability to one-electron oxidation. These theoretical results elucidate the diversity of free radical scavenging mechanisms for aminic RTAs, and has significant implications for the rational design of new aminic RTAs.

Aromaticity and Electron Density of Hypericin

The influence of the substituents on the geometry of the central ring system of hypericin has been analyzed. Substitution that causes flattening of the hypericin central rings is connected with introducing the aromatic character of the empty rings. All the hypericin rings have an aromatic character illustrated by the Harmonic Oscillator Measure of Aromaticity (HOMA), Nucleus Independent Chemical Shift (NICS), Fluctuation Index (FLU), and Ellipticity Index (EL) indices. Quantum Theory of Atoms in Molecules (QTAIM) and Natural Bond Orbital (NBO) analyses performed on 7,14-dihydrophenanthro[1,10,9,8-opqra]perylene, its substituted analogues, and hypericin show an influence of this substitution on electron density of the central rings.

Interplay of Ring Puckering and Hydrogen Bonding in Deoxyribonucleosides

The Cremer-Pople ring puckering analysis and the Konkoli-Cremer local mode analysis supported by the topological analysis of the electron density were applied for the first comprehensive analysis of the interplay between deoxyribose ring puckering and intramolecular H-bonding in 2'-deoxycytidine, 2'-deoxyadenosine, 2'-deoxythymidine, and 2'-deoxyguanosine. We mapped for each deoxyribonucleoside the complete conformational energy surface and the corresponding pseudorotation path. We found only incomplete pseudorotation cycles, caused by ring inversion, which we coined as pseudolibration paths. On each pseudolibration path a global and a local minimum separated by a transition state were identified. The investigation of H-bond free deoxyribonucleoside analogs revealed that removal of the H-bond does not restore the full conformational flexibility of the sugar ring. Our work showed that ring puckering predominantly determines the conformational energy; the larger the puckering amplitude, the lower the conformational energy. In contrast no direct correlation between conformational energy and H-bond strength was found. The longest and weakest H-bonds are located in the local minimum region, whereas the shortest and strongest H-bonds are located outside the global and local minimum regions at the turning points of the pseudolibration paths, i.e., H-bonding determines the shape and length of the pseudolibration paths. In addition to the H-bond strength, we evaluated the covalent/electrostatic character of the H-bonds applying the Cremer-Kraka criterion of covalent bonding. H-bonding in the puric bases has a more covalent character whereas in the pyrimidic bases the H-bond character is more electrostatic. We investigated how the mutual orientation of the CH₂OH group and the base influences H-bond formation via two geometrical parameters describing the rotation of the substituents perpendicular to the sugar ring and their tilting relative to the ring center. According to our results, rotation is more important for H-bond formation. In addition we assessed the influence of the H-bond acceptor, the lone pair (N, respectively O), via the delocalization energy. We found larger delocalization energies corresponding to stronger H-bonds for the puric bases. The global minimum conformation of 2'-deoxyguanosine has the strongest H-bond of all conformers investigated in this work with a bond strength of 0.436 which is even stronger than the H-bond in the water dimer (0.360). The application of our new analysis to DNA deoxyribonucleotides and to unnatural base pairs, which have recently drawn a lot of attention, is in progress.

Side chain to main chain hydrogen bonds stabilize a polyglutamine helix in a transcription factor

Polyglutamine (polyQ) tracts are regions of low sequence complexity frequently found in transcription factors. Tract length often correlates with transcriptional activity and expansion beyond specific thresholds in certain human proteins is the cause of polyQ disorders. To study the structural basis of the association between tract length, transcriptional activity and disease, we addressed how the conformation of the polyQ tract of the androgen receptor, associated with spinobulbar muscular atrophy (SBMA), depends on its length. Here we report that this sequence folds into a helical structure stabilized by unconventional hydrogen bonds between glutamine side chains and main chain carbonyl groups, and that its helicity directly correlates with tract length. These unusual hydrogen bonds are bifurcate with the conventional hydrogen bonds stabilizing α-helices. Our findings suggest a plausible rationale for the association between polyQ tract length and androgen receptor transcriptional activity and have implications for establishing the mechanistic basis of SBMA.

Polyglutamine (polyQ) tracts are low-complexity regions and their expansion is linked to certain neurodegenerative diseases. Here the authors combine experimental and computational approaches to find that the length of the androgen receptor polyQ tract correlates with its helicity and show that the polyQ helical structure is stabilized by hydrogen bonds between the Gln side chains and main chain carbonyl groups.

C-H···O Hydrogen Bonding. The Prototypical Methane-Formaldehyde System: A Critical Assessment

Distinguishing the functionality of C-H···O hydrogen bonds (HBs) remains challenging, because their properties are difficult to quantify reliably. Herein, we present a study of the model methane-formaldehyde complex (MFC). Six stationary points on the MFC potential energy surface (PES) were obtained at the CCSD(T)/ANO2 level. The CCSDT(Q)/CBS interaction energies of the conformers range from only -1.12 kcal mol^-1 to -0.33 kcal mol^-1, denoting a very flat PES. Notably, only the lowest energy stationary point (MFC1) corresponds to a genuine minimum, whereas all other stationary points-including the previously studied ideal case of a_e(C-H···O) = 180°-exhibit some degree of freedom that leads to MFC1. Despite the flat PES, we clearly see that the HB properties of MFC1 align with those of the prototypical water dimer O-H···O HB. Each HB property generally becomes less prominent in the higher-energy conformers. Only the MFC1 conformer prominently exhibits (1) elongated C-H donor bonds, (2) attractive C-H···O═C interactions, (3) n(O) → σ*(C-H) hyperconjugation, (4) critical points in the electron density from Bader's method and from the noncovalent interactions method, (5) positively charged donor hydrogen, and (6) downfield NMR chemical shifts and nonzero ²J(C_M-H_M···O_F) coupling constants. Based on this research, some issues merit further study. The flat PES hinders reliable determinations of the HB-induced shifts of the C-H stretches; a similarly difficult challenge is observed for the experiment. The role of charge transfer in HBs remains an intriguing open question, although our BLW and NBO computations suggest that it is relevant to the C-H···O HB geometries. These issues notwithstanding, the prominence of the HB properties in MFC1 serves as clear evidence that the MFC is predominantly bound by a C-H···O HB.

Bond softness sensitive bond-valence parameters for crystal structure plausibility tests

Based on a description of bond valence as a function of valence electron density, a systematic bond softness sensitive approach to determine bond-valence parameters and related quantities such as coordination numbers is elaborated and applied to determine bond-valence parameters for 706 cation-anion pairs. While the approach is closely related to the earlier softBV parameter set, the new softNC1 parameters proposed in this work may be simpler to apply in plausibility checks of crystal structures, as they follow the first coordination shell convention. The performance of this softNC1 bond-valence parameter set is compared with that of the previously derived softBV parameter set that also factors in contributions from higher coordination shells, and with a benchmarking parameter set that has been optimized following the conventional choice of a universal value of the bond-valence parameter b. The results show that a systematic adaptation of the bond-valence parameters to the bond softness leads to a significant improvement in the bond-valence parameters, particularly for bonds involving soft anions, and is safer than individual free refinements of both R₀ and b from a limited number of reference cation environments.

Quantum Chemical Modeling of Hydrogen Bonding in Ionic Liquids

Hydrogen bonding (H-bonding) is an important and very general phenomenon. H-bonding is part of the basis of life in DNA, key in controlling the properties of water and ice, and critical to modern applications such as crystal engineering, catalysis applications, pharmaceutical and agrochemical development. H-bonding also plays a significant role for many ionic liquids (IL), determining the secondary structuring and affecting key physical parameters. ILs exhibit a particularly diverse and wide range of traditional as well as non-standard forms of H-bonding, in particular the doubly ionic H-bond is important. Understanding the fundamental nature of the H-bonds that form within ILs is critical, and one way of accessing this information, that cannot be recovered by any other computational method, is through quantum chemical electronic structure calculations. However, an appropriate method and basis set must be employed, and a robust procedure for determining key structures is essential. Modern generalised solvation models have recently been extended to ILs, bringing both advantages and disadvantages. QC can provide a range of information on geometry, IR and Raman spectra, NMR spectra and at a more fundamental level through analysis of the electronic structure.