Chalcogen Noncovalent Interactions between Diazines and Sulfur Oxides in Supramolecular Circular Chains

The noncovalent chalcogen interaction between SO2/SO3 and diazines was studied through a dispersion-corrected DFT Kohn-Sham molecular orbital together with quantitative energy decomposition analyses. For this, supramolecular circular chains of up to 12 molecules were built with the aim of checking the capability of diazine molecules to detect SO2/SO3 compounds within the atmosphere. Trends in the interaction energies with the increasing number of molecules are mainly determined by the Pauli steric repulsion involved in these σ-hole/π-hole interactions. But more importantly, despite the assumed electrostatic nature of the involved interactions, the covalent component also plays a determinant role in its strength in the involved chalcogen bonds. Noticeably, π-hole interactions are supported by the charge transfer from diazines to SO2/SO3 molecules. Interaction energies in these supramolecular complexes are not only determined by the S···N bond lengths but attractive electrostatic and orbital interactions also determine the trends. These results should allow us to establish the fundamental characteristics of chalcogen bonding based on its strength and nature, which is of relevance for the capture of sulfur oxides.

Effects of Chalcogen Atoms on Excited-State Double-Proton Transfer Behavior for 3,6-bis(4,5-Dihydroxyoxazo-2-yl)benzene-1,2-diol Derivatives: A Computational Investigation

The impact of the chalcogen atomic electronegativity (O, S, and Se atoms) of new organic molecules on excited-state dynamical reactions is self-evident. Inspired by this kind of distinguished photochemical characteristic, in this work, we performed a computational investigation of chalcogen-substituted 3,6-bis(4,5-dihydroxyoxazo-2-yl)benzene-1,2-diol (BDYBD) derivatives (i.e., BDYBD-O, BDYBD-S, and BDYBD-Se). In this paper, we pay close attention to characteristic BDYBD derivatives that contain intramolecular double hydrogen bonds (O1-H2···N3 and O4-H5···N6). The main goal of this study was to explore how changes in atomic electronegativity affect the way hydrogen bonds interact and how excited molecules affect transfer protons. We go into further detail in the main text of the paper. By fixing our attention to geometrical variations and infrared (IR) vibrational spectra between the S0 and S1 states, exploring hydrogen bonding behaviors using the core-valence bifurcation (CVB) index, and simulating hydrogen bonding energy (EHB) via the atom in molecule (AIM) method, we clarified the photo-induced strengthened dual hydrogen bonding interactions that facilitate the excited-state dual-proton transfer (ESDPT) behavior of BDYBD derivatives. The reorganization of charge stemming from photoexcitation further verifies the tendencies of ESDPT reactions. We relied on constructing potential energy surfaces (PESs) by adopting a restrictive optimization approach, and herein, we finally clarify the gradual ESDPT mechanism of BDYBD derivatives. Particularly, we confirm that the variation in chalcogen atomic electronegativity has a regulatory effect on the ESDPT behavior of BDYBD derivatives; that is, the lower the atomic electronegativity, the more favorable it is for the gradual ESDPT reaction.

Complementary base pair interactions between different rare tautomers of the second-generation artificial genetic alphabets

The functionality of a semisynthetic DNA in the biological environment will depend on the base pair nature of its complementary base pairs. To understand this, base pair interactions between complementary bases of recently proposed eight second-generation artificial nucleobases are studied herein by considering their rare tautomeric conformations and a dispersion-corrected density functional theoretic method. It is found that the binding energies of two hydrogen-bonded complementary base pairs are more negative than those of the three hydrogen-bonded base pairs. However, as the former base pairs are endothermic, the semisynthetic duplex DNA would involve the latter base pairs.

Solvent effects on the NMR shieldings of stacked DNA base pairs

Stacking effects are among the most important effects in DNA. We have recently studied their influence in fragments of DNA through the analysis of NMR magnetic shieldings, firstly in vacuo. As a continuation of this line of research we show here the influence of solvent effects on the shieldings through the application of both explicit and implicit models. We found that the explicit solvent model is more appropriate for consideration due to the results matching better in general with experiments, as well as providing clear knowledge of the electronic origin of the value of the shieldings. Our study is grounded on a recently developed theoretical model of our own, by which we are able to learn about the magnetic effects of given fragments of DNA molecules on selected base pairs. We use the shieldings of the atoms of a central base pair (guanine-cytosine) of a selected fragment of DNA molecules as descriptors of physical effects, like π-stacking and solvent effects. They can be taken separately and altogether. The effect of π-stacking is introduced through the addition of some pairs above and below of the central base pair, and now, the solvent effect is considered including a network of water molecules that consist of two solvation layers, which were fixed in the calculations performed in all fragments. We show that the solvent effects enhance the stacking effects on the magnetic shieldings of atoms that belong to the external N-H bonds. The net effect is of deshielding on both atoms. There is also a deshielding effect on the carbon atoms that belong to CO bonds, for which the oxygen atom has an explicit hydrogen bond (HB) with a solvent water molecule. Solvent effects are found to be no higher than a few percent of the total value of the shieldings (between 1% and 5%) for most atoms, although there are few for which such an effect can be higher. There is one nitrogen atom, the acceptor of the HB between guanine and cytosine, that is more highly shielded (around 15 ppm or 10%) when the explicit solvent is considered. In a similar manner, the most external nitrogen atom of cytosine and the hydrogen atom that is bonded to it are highly deshielded (around 10 ppm for nitrogen and around 3 ppm for hydrogen).

Affinity and Correlation in DNA

A statistical analysis of important DNA sequences and related proteins has been performed to study the relationships between monomers, and some general considerations about these macromolecules can be provided from the results. First, the most important relationship between sites in all the DNA sequences examined is that between two consecutive base pairs. This is an indication of an energetic stabilization due to the stacking interaction of these couples of base pairs. Secondly, the difference between human chromosome sequences and their coding parts is relevant both in the relationships between sites and in some specific compositional rules, such as the second Chargaff rule. Third, the evidence of the relationship in two successive triplets of DNA coding sequences generates a relationship between two successive amino acids in the proteins. This is obviously impossible if all the relationships between the sites are statistical evidence and do not involve causes; therefore, in this article, due to stacking interactions and this relationship in coding sequences, we will divide the concept of the relationship between sites into two concepts: affinity and correlation, the first with physical causes and the second without. Finally, from the statistical analyses carried out, it will emerge that the human genome is uniform, with the only significant exception being the Y chromosome.

B-DNA Structure and Stability: The Role of Nucleotide Composition and Order

We have quantum chemically analyzed the influence of nucleotide composition and sequence (that is, order) on the stability of double-stranded B-DNA triplets in aqueous solution. To this end, we have investigated the structure and bonding of all 32 possible DNA duplexes with Watson-Crick base pairing, using dispersion-corrected DFT at the BLYP-D3(BJ)/TZ2P level and COSMO for simulating aqueous solvation. We find enhanced stabilities for duplexes possessing a higher GC base pair content. Our activation strain analyses unexpectedly identify the loss of stacking interactions within individual strands as a destabilizing factor in the duplex formation, in addition to the better-known effects of partial desolvation. Furthermore, we show that the sequence-dependent differences in the interaction energy for duplexes of the same overall base pair composition result from the so-called "diagonal interactions" or "cross terms". Whether cross terms are stabilizing or destabilizing depends on the nature of the electrostatic interaction between polar functional groups in the pertinent nucleobases.

Matrix effects on hydrogen bonding and proton transfer in fluoropyridine - HCl complexes

We report an extensive computational and spectroscopic study of several fluoropyridine-HCl complexes, and the parent, pyridine-HCl system. Matrix-IR spectra for pentafluoropyridine-HCl, 2,6-difluororpyridine-HCl, and 3,5-difluororpyridine-HCl in solid neon exhibit shifts for the H-Cl stretching band that parallel the effects of fluorination on hydrogen-bond strength. Analogous spectral shifts observed across various host environments (solid neon, argon, and nitrogen) for pentafluoropyridine-HCl and 2,6-difluororpyridine-HCl convey a systematically varying degree of matrix stabilization on the hydrogen bonds in these complexes. An extended quantum-chemical study of pyridine-HCl and eight fluorinated analogs, including 2-, 3-, and 4-fluoropyridine-HCl, 2,6- and 3,5-difluororpyridine-HCl, 2,4,6- and 3,4,5-trifluropyridine-HCl, as well as pentafluoropyridine-HCl, was also performed. Equilibrium structures and binding energies for the gas-phase complexes illustrate two clear trends in how fluorine substitution affects hydrogen bond strength; increasing fluorination weakens these interactions, yet substitution at the 2- and 6-positions has the most pronounced effect. Bonding analyses for a select subset of these systems reveal shifts in electron density that accompany hydrogen bonding, and most notably, the values of the electron density at the N-H bond critical points among the stronger systems in this subset significantly exceed those typical for moderately strong hydrogen-bonds. We also explored the effects of dielectric media on the structural and bonding properties of these systems. For pyridine-HCl, 3-fluoropyridine-HCl, and 3,5-difluororpyridine-HCl, a transition to proton transfer-type structures is observed at ε-values of 1.2, 1.5, and 2.0, respectively. This is signaled by key structural changes, as well as an increase in the negative charge on the chorine, and dramatic shifts in topological properties of the H-Cl and N-H bonds. In the case of pentafluoropyridine-HCl, and 2,6-difluororpyridine-HCl, we do not predict proton transfer in dielectric media up to ε = 20.0. However, there are clear indications that the media enhance hydrogen-bond strength, and moreover, these observations are completely consistent with the experimental IR spectra.

Unveiling the complex pattern of intermolecular interactions responsible for the stability of the DNA duplex

Herein, we provide new insights into the intermolecular interactions responsible for the intrinsic stability of the duplex structure of a large portion of human B-DNA by using advanced quantum mechanical methods. Our results indicate that (i) the effect of non-neighboring bases on the inter-strand interaction is negligibly small, (ii) London dispersion effects are essential for the stability of the duplex structure, (iii) the largest contribution to the stability of the duplex structure is the Watson-Crick base pairing - consistent with previous computational investigations, (iv) the effect of stacking between adjacent bases is relatively small but still essential for the duplex structure stability and (v) there are no cooperativity effects between intra-strand stacking and inter-strand base pairing interactions. These results are consistent with atomic force microscope measurements and provide the first theoretical validation of nearest neighbor approaches for predicting thermodynamic data of arbitrary DNA sequences.

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.

Can Polarity-Inverted Surfactants Self-Assemble in Nonpolar Solvents?

We investigate the self-assembly process of a surfactant with inverted polarity in water and cyclohexane using both all-atom and coarse-grained hybrid particle-field molecular dynamics simulations. Unlike conventional surfactants, the molecule under study, proposed in a recent experiment, is formed by a rigid and compact hydrophobic adamantane moiety, and a long and floppy triethylene glycol tail. In water, we report the formation of stable inverted micelles with the adamantane heads grouping together into a hydrophobic core and the tails forming hydrogen bonds with water. By contrast, microsecond simulations do not provide evidence of stable micelle formation in cyclohexane. Validating the computational results by comparison with experimental diffusion constant and small-angle X-ray scattering intensity, we show that at laboratory thermodynamic conditions the mixture resides in the supercritical region of the phase diagram, where aggregated and free surfactant states coexist in solution. Our simulations also provide indications as to how to escape this region to produce thermodynamically stable micellar aggregates.

Revisiting the Potential Energy Surface of the Stacked Cytosine Dimer: FNO-CCSD(T) Interaction Energies, SAPT Decompositions, and Benchmarking

Nucleobase stacking interactions are crucial for the stability of nucleic acids. This study investigates base stacking energies of the cytosine homodimer in different configurations, including intermolecular separation plots, detailed twist dependence, and displaced structures. Highly accurate ab initio quantum chemical single point energies using an energy function based on MP2 complete basis set extrapolation ([6 → 7]ZaPa-NR) and a CCSD(T)/cc-pVTZ-F12 high-level correction are presented as new reference data, providing the most accurate stacking energies of nucleobase dimers currently available. Accurate SAPT2+(3)δMP2 energy decomposition is used to obtain detailed insights into the nature of base stacking interactions at varying vertical distances and twist values. The ab initio symmetry adapted perturbation theory (SAPT) energy decomposition suggests that the base stacking originates from an intricate interplay between dispersion attraction, short-range exchange-repulsion, and Coulomb interaction. The interpretation of the SAPT data is a complex issue as key energy terms vary substantially in the region of optimal (low energy) base stacking geometries. Thus, attempts to highlight one leading stabilizing SAPT base stacking term may be misleading and the outcome strongly depends on the used geometries within the range of geometries sampled in nucleic acids upon thermal fluctuations. Modern dispersion-corrected density functional theory (among them DSD-BLYP-D3, ωB97M-V, and ωB97M-D3BJ) is benchmarked and often reaches up to spectroscopic accuracy (below 1 kJ/mol). The classical AMBER force field is benchmarked with multiple different sets of point-charges (e.g. HF, DFT, and MP2-based) and is found to produce reasonable agreement with the benchmark data.

Structure and bonding in triorganotin chlorides: a perspective from energy decomposition analysis

The Sn-Cl chemical bond of four organotin halides (Me₃SnCl, Et₃SnCl, Bu₃SnCl, and Ph₃SnCl) was studied by using relativistic density functional theory in combination with a quantitative energy decomposition analysis to explain the formation of charged species. The σ orbital is the dominant contributor to the stabilization of the Sn-Cl bond, and the π-orbital interactions also have a significant contribution to the stabilization of Ph₃Sn⁺ cation when the aromatic groups are bonded to the tin atom. The aromaticity of the phenyl groups delocalizes the positive charge, donating electrons to tin atom by conjugation. Although Me₃SnCl and Ph₃SnCl are constituted by groups which the size of the substituents is different, the interaction energies obtained with the energy decomposition analysis present similar values, which also occur with the thermodynamic parameters. Graphical abstract Organotin compounds have widely studied as a potential antitumoral agent. The mechanism in triorganotin compounds includes the formation of cation species, R₃Sn⁺. This article studies the influence of the R groups on the rupture of Sn-Cl bond using the fragment analysis and quantitative energy decomposition analysis.

Molecular Mechanisms of DNA Replication and Repair Machinery: Insights from Microscopic Simulations

Reproduction, the hallmark of biological activity, requires making an accurate copy of the genetic material to allow the progeny to inherit parental traits. In all living cells, the process of DNA replication is carried out by a concerted action of multiple protein species forming a loose protein-nucleic acid complex, the replisome. Proofreading and error correction generally accompany replication but also occur independently, safeguarding genetic information through all phases of the cell cycle. Advances in biochemical characterization of intracellular processes, proteomics and the advent of single-molecule biophysics have brought about a treasure trove of information awaiting to be assembled into an accurate mechanistic model of the DNA replication process. In this review, we describe recent efforts to model elements of DNA replication and repair processes using computer simulations, an approach that has gained immense popularity in many areas of molecular biophysics but has yet to become mainstream in the DNA metabolism community. We highlight the use of diverse computational methods to address specific problems of the fields and discuss unexplored possibilities that lie ahead for the computational approaches in these areas.

Insights into N-Heterocyclic Carbene-Catalyzed Oxidative α-C(sp3)-H Activation of Aliphatic Aldehydes and Cascade [2 + 2] Cycloaddition with Ketimines

Predicting the chemoselectivity of [2 + 2] cyclizations is an important challenge in organic chemistry. Herein, we provided a valuable case for this issue. Density functional theory calculations were performed to systematically study the possible mechanisms and origin of selectivities for the N-heterocyclic carbene (NHC)-catalyzed oxidative α-C(sp³)-H activation of aliphatic aldehydes and the cascade [2 + 2] cycloaddition with ketimines. The [2 + 2] cycloaddition of azolium enolate intermediates to the C═N bond, rather than the C═O bond of ketimine, is revealed to be determined by chemo- and stereoselectivity. By comparing the energy gap between the frontier molecular orbitals (FMOs) of the two reacting parts involved in the [2 + 2] cycloaddition transition states, we propose a new strategy to determine the origin of the reaction chemoselectivity. Moreover, the local nucleophilic index can efficiently predict the active site of ketimines. Further analyses illustrate that NHC can increase the nucleophilicity of aldehydes and the acidity of the α-C(sp³)-H bond, and 3,3',5,5'-tetra- tert-butyl diphenoquinone (DQ) acts as an oxidant and promotes α-C(sp³)-H bond deprotonation. This work is useful not only for understanding the NHC-catalyzed oxidative [2 + 2] annulation but also for developing new applications of the FMO theory in organocatalytic cyclizations.

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.

Prediction on the origin of chemoselectivity in Lewis base-mediated competition cyclizations between allenoates and chalcones: a computational study

The HOMO_TSs and p-orbital contributions of the center atoms were used to predict the origin of chemoselectivity in this work.

RNA versus DNA G-Quadruplex: The Origin of Increased Stability

DNA quadruplexes have been the subject of investigation because of their biological relevance and because of their potential application in supramolecular chemistry. Similarly, RNA quadruplexes are now gaining increasing attention. Although DNA and RNA quadruplexes are structurally very similar, the latter show higher stability. In this study we report dispersion‐corrected density functional theory (DFT‐D) quantum chemical calculations that were undertaken to understand the difference in stabilities of RNA and DNA quadruplexes. The smallest meaningful model of a stack of quartets, interacting with alkali metal cations, was simulated in an aqueous environment. The energy decomposition analysis allows for in‐depth examination of the interaction energies, emphasising the role of noncovalent interactions and better electrostatics in determining RNA‐GQs higher stabilities, particularly pinpointing the role of the extra 2′‐OH groups. Furthermore, our computations present new insights on why the cation is required for self‐assembly: unexpectedly the cation is not necessary to relieve the repulsion between the oxygen atoms in the central cavity, but it is needed to overcome the entropic penalty.

Nucleophilic Substitution in Solution: Activation Strain Analysis of Weak and Strong Solvent Effects

We have quantum chemically studied the effect of various polar and apolar solvents on the shape of the potential energy surface (PES) of a diverse collection of archetypal nucleophilic substitution reactions at carbon, silicon, phosphorus, and arsenic by using density functional theory at the OLYP/TZ2P level. In the gas phase, all our model SN 2 reactions have single-well PESs, except for the nucleophilic substitution reaction at carbon (SN 2@C), which has a double-well energy profile. The presence of the solvent can have a significant effect on the shape of the PES and, thus, on the nature of the SN 2 process. Solvation energies, charges on the nucleophile or leaving group, and structural features are compared for the various SN 2 reactions in a spectrum of solvents. We demonstrate how solvation can change the shape of the PES, depending not only on the polarity of the solvent, but also on how the charge is distributed over the interacting molecular moieties during different stages of the reaction. In the case of a nucleophilic substitution at three-coordinate phosphorus, the reaction can be made to proceed through a single-well [no transition state (TS)], bimodal barrier (two TSs), and then through a unimodal transition state (one TS) simply by increasing the polarity of the solvent.