Intramolecular Repulsion Visible Through the Electrostatic Potential in Disulfides: Analysis via Varying Iso-density Envelopes

Radial Behavior of Electrostatic Potentials of Atoms and Ions Revisited: Isotropy and Anisotropy

In this paper we revisit earlier work relating to monoatomic atoms and ions published by pioneers in the area of electrostatic potentials. We include plots of the radial distributions of the electrostatic potentials for spherically symmetric atoms and cations, and for singly, doubly and triply negative anions. For atoms with anisotropy in their densities and electrostatic potentials, such as the halonium cations, it is shown how the molecular surface approach for plotting electrostatic potentials complements that achieved by directional radial distributions.

Unveiling the Intermolecular Interactions between Drug 5-Fluorouracil and Watson-Crick/Hoogsteen Base Pairs: A Computational Analysis

The adsorption of 5-fluorouracil (5FU) on Watson-Crick (WC) base pairs and Hoogsteen (HT) base pairs has been studied using the dispersion-corrected density functional theory (DFT). The adsorption, binding energy, and thermochemistry for the drug 5FU on the WC and HT base pairs were determined. The most stable geometries were near planar geometry, and 5FU has a higher preference for WC than HT base pairs. The adsorption energies of 5FU on nucleobase pairs are consistently higher than pristine nucleobase pairs, indicating that nucleobase pair cleavage is less likely during the adsorption of the 5FU drug. The enthalpy change for the formation of 5FU-DNA base pairs is higher than that for the formation of 5FU-nucleobases and is enthalpy-driven. The E gap of AT base pairs is higher, suggesting that their chemical reactivity toward further reaction would be less than that of GC base pairs. The electron density difference (EDD) analysis shows a significant decrease in electron density in aromatic regions on the purine bases (adenine/guanine) compared to the pyrimidine bases. The MESP diagram of the stable 5FU-nucleobase pair complexes shows a directional interaction, with the positive regions in a molecule interacting with the negative region of other molecules. The atoms in molecule analysis show that the ρ(r) values of C=O···H-N are higher than those of N···H/N-H···O. The N···H intermolecular bonds between the base pair/drug and nucleobases are weak, closed shell interactions and are electrostatic in nature. The noncovalent interaction analysis shows that several new spikes are engendered along with an increase in their strength, which indicates that the H-bonding interactions are stronger and play a dominant role in stabilizing the complexes. Energy decomposition analysis shows that the drug-nucleobase pair complex has a marginal increase in the electrostatic contributions compared to nucleobase pair complexes.

Probing intramolecular interactions using molecular electrostatic potentials: changing electron density contours to unveil both attractive and repulsive interactions

We focus on intramolecular interactions, using the electrostatic potential plotted on iso-density surfaces to lead the way. We show that plotting the electrostatic potential on varying iso-density envelopes much closer to the nuclei than the commonly used 0.001 a.u. contour can reveal the driving forces for such interactions, whether they be stabilizing or destabilizing. Our approach involves optimizing the structures of molecules exhibiting intramolecular interactions and then finding the contour of the electronic density which allows the interacting atoms to be separated; we call this the nearly-touching contour. The electrostatic potential allows then to identify the intramolecular interactions as either attractive or repulsive. The discussed 1,5- and 1,6-intramolecular interactions in o-bromophenol and o-nitrophenol are attractive, while the interactions between terminal methyl hydrogens in diethyl disulfides (as shown recently) and those between the closest hydrogens in planar biphenyl and phenanthrene are clearly repulsive in nature. For the attractive 1,4-interactions in trinitromethane and chlorotrinitromethane, and the 1,3-S⋯N and the 1,4-Si⋯N interactions in the ClH2Si(CH2)nNH2 series, the lack of (3,-1) bond critical points has often been cited as reason to not identify such interactions as attractive in nature. Here, by looking at the nearly-touching contours we see that bond critical points are neither necessary nor sufficient for attractive interactions, as others have pointed out, and in some instances also pointing to repulsive interactions, as the examples of planar biphenyl and phenanthrene discussed in this work show.

The Formation of σ-Hole Bonds: A Physical Interpretation

This paper discusses two quite different computational experiments relating to the formation of σ-hole bonds A···B. The first involves looking at the complex at equilibrium and finding the contour X of the electronic density which allows the iso-density envelopes of A and B to be nearly touching. This contour increases, becoming closer to the nuclei, as the strength of the interaction increases. The second experiment involves allowing A and B to approach each other, with the aim of finding the distance at which their 0.001 a.u. iso-density envelopes are nearly merging into one envelope. What is found in the second experiment may be somewhat surprising, in that the ratio of the distance between interacting atoms at this nearly merging point-divided by the sum of the van der Waals radii of these atoms-covers a narrow range, typically between 1.2 and 1.3. It is intriguing to note that for the dataset presented, approaching molecules attracted to each other appear to do so unknowing of the strength of their ultimate interaction. This second experiment also supports the notion that one should expect favorable interactions, in some instances, to have close contacts significantly greater than the sums of the van der Waals radii.