σ-Hole and Lone-Pair Hole Interactions in Chalcogen-Containing Complexes: A Comparative Study

Cell-free biodegradable electroactive scaffold for urinary bladder regeneration

Tissue engineering heavily relies on cell-seeded scaffolds to support the complex biological and mechanical requirements of a target organ. However, in addition to safety and efficacy, translation of tissue engineering technology will depend on manufacturability, affordability, and ease of adoption. Therefore, there is a need to develop scalable biomaterial scaffolds with sufficient bioactivity to eliminate the need for exogenous cell seeding. Herein, we describe synthesis, characterization, and implementation of an electroactive biodegradable elastomer for urinary bladder tissue engineering. To create an electrically conductive and mechanically robust scaffold to support bladder tissue regeneration, we developed a phase-compatible functionalization method wherein the hydrophobic conductive polymer poly(3,4-ethylenedioxythiophene) (PEDOT) was polymerized in situ within a similarly hydrophobic citrate-based elastomer poly(octamethylene-citrate-co-octanol) (POCO) film. We demonstrate the efficacy of this film as a scaffold for bladder augmentation in athymic rats, comparing PEDOT-POCO scaffolds to mesenchymal stromal cell-seeded POCO scaffolds. PEDOT-POCO recovered bladder function and anatomical structure comparably to the cell-seeded POCO scaffolds and significantly better than non-cell seeded POCO scaffolds. This manuscript reports: (1) a new phase-compatible functionalization method that confers electroactivity to a biodegradable elastic scaffold, and (2) the successful restoration of the anatomy and function of an organ using a cell-free electroactive scaffold.

Yet another perspective on hole interactions, part II: lp-hole vs. lp-hole interactions

Most of the experimental and theoretical work in hole interactions (HIs) is mainly focused on exploiting the nature and characteristics of σ and π-holes. In this perspective, we focus our attention on understanding the origin and properties of lone-pair holes. These holes are present on an atom opposite to its lone-pair region. Utilizing some new and old examples, such as X₃N/P⋯F^- (X = F/Cl/Br/I), F-Cl/Br/I⋯H₃P⋯NCH and H₃B-NBr₃ along with other molecular systems, we explored to what extent these lp-holes participate in lp-hole interactions, if they participate at all.

Nature and Strength of Weak O⋅⋅⋅O Interactions in Nitryl Halide Dimers

The use of real space functions and molecular graphs has pushed some chemists to wonder: Are interactions between negatively charged oxygen atoms possible? In this contribution we analyze whether there is a real interaction between oxygen atoms in nitryl halide dimers (XNO₂ )₂ (X=F, Cl, Br and I) and in tetranitromethane and derivatives. Based on ab-initio and density functional theories (DFT) methods, we show these complexes are weakly stabilized. Energy decomposition analyses based on local molecular orbitals (LMOEDA) and interacting quantum atoms (IQA) reveal both dispersion and exchange play a crucial role in the stabilization of these complexes. Electron charge density and IQA analyses indicate that the oxygen atoms are connected by privileged exchange channels. In addition, electrostatic interactions between O and N atoms are also vital for the stabilization of the complexes. Finally, a reasonable explanation is given for the dynamic behavior of nitryl groups in tetranitromethane and derivatives.

Synthesis, crystal structure, DFT and molecular docking studies of N-acetyl-2,4-[diaryl-3-azabicyclo[3.3.1]nonan-9-yl]-9-spiro-4'-acetyl-2'-(acetylamino)-4',9-dihydro-[1',3',4']-thiadiazoles: A potential SARS-nCoV-2 Mpro (COVID-19) inhibitor

In this paper, we describe the synthesis and crystal structure analysis of N-acetyl-2,4-[diphenyl-3-azabicyclo[3.3.1]nonan-9-yl]-9-spiro-4′-acetyl-2′-(acetylamino)-4′,9-dihydro-[1′,3′,4′]-thiadiazole (3a) and N-acetyl- 2,4-[bis(p-methoxyphenyl)-3-azabicyclo[3.3.1]nonan-9-yl]-9-spiro-4′-acetyl-2′-(acetylamino)-4′,9-dihydro-[1′,3′,4′]-thiadiazole (3b). The title compounds 3a and 3b are characterized by 1D NMR and single crystal x-ray diffraction analysis. Non-covalent interactions in a molecule were identified by Hirshfeld surface (d_norm contacts and 2D fingerprint plot) analysis. In addition, the existence of chalcogen bond (S•••O bond) in the molecular structures (3a and 3b) are described by NCI-RDG and QTAIM analysis. NBO analysis is employed to describe the orbital interactions and electron transfer between sulfur and oxygen atoms. Molecular docking is carried out for compounds 3a and 3b with COVID-19 viral protein SARS-nCoV-2 M^pro (PDB ID: 6LU7).

Theoretical study of the NO3 radical reaction with CH2ClBr, CH2ICl, CH2BrI, CHCl2Br, and CHClBr2

The potential reaction of the nitrate radical (NO₃), the main nighttime atmospheric oxidant, with five alkyl halides, halons (CH₂ClBr, CH₂ICl, CH₂BrI, CHCl₂Br, and CHClBr₂) has been studied theoretically. The most favorable reaction corresponds to a hydrogen atom transfer. The stationary points on the potential energy surfaces of these reactions have been characterized. The reactions can be classified into two groups based on the number of hydrogen atoms in the halon molecules (1 or 2). The reactions with halons with only one hydrogen atom show more exothermic profiles than those with two hydrogen atoms. In addition, the kinetics of the reaction of NO₃ + CH₂BrI was studied in much higher detail using a multi-well Master Equation solver as a representative example of the nitrate radical reactivity against these halocarbons. These results indicate that the chemical lifetime of the alkyl halides would not be substantially affected by nitrate radical reactions, even in the case of NO₃-polluted atmospheric conditions.

The Pnictogen Bond: The Covalently Bound Arsenic Atom in Molecular Entities in Crystals as a Pnictogen Bond Donor

In chemical systems, the arsenic-centered pnictogen bond, or simply the arsenic bond, occurs when there is evidence of a net attractive interaction between the electrophilic region associated with a covalently or coordinately bound arsenic atom in a molecular entity and a nucleophile in another or the same molecular entity. It is the third member of the family of pnictogen bonds formed by the third atom of the pnictogen family, Group 15 of the periodic table, and is an inter- or intramolecular noncovalent interaction. In this overview, we present several illustrative crystal structures deposited into the Cambridge Structure Database (CSD) and the Inorganic Chemistry Structural Database (ICSD) during the last and current centuries to demonstrate that the arsenic atom in molecular entities has a significant ability to act as an electrophilic agent to make an attractive engagement with nucleophiles when in close vicinity, thereby forming σ-hole or π-hole interactions, and hence driving (in part, at least) the overall stability of the system’s crystalline phase. This overview does not include results from theoretical simulations reported by others as none of them address the signatory details of As-centered pnictogen bonds. Rather, we aimed at highlighting the interaction modes of arsenic-centered σ- and π-holes in the rationale design of crystal lattices to demonstrate that such interactions are abundant in crystalline materials, but care has to be taken to identify them as is usually done with the much more widely known noncovalent interactions in chemical systems, halogen bonding and hydrogen bonding. We also demonstrate that As-centered pnictogen bonds are usually accompanied by other primary and secondary interactions, which reinforce their occurrence and strength in most of the crystal structures illustrated. A statistical analysis of structures deposited into the CSD was performed for each interaction type As···D (D = N, O, S, Se, Te, F, Cl, Br, I, arene’s π system), thus providing insight into the typical nature of As···D interaction distances and ∠R–As···D bond angles of these interactions in crystals, where R is the remainder of the molecular entity.

σ-Hole and LP-Hole Interactions of Pnicogen···Pnicogen Homodimers under the External Electric Field Effect: A Quantum Mechanical Study

σ-Hole and lone-pair (lp)-hole interactions within σ-hole···σ-hole, σ-hole···lp-hole, and lp-hole···lp-hole configurations were comparatively investigated on the pnicogen···pnicogen homodimers (PCl3)2, for the first time, under field-free conditions and the influence of the external electric field (EEF). The electrostatic potential calculations emphasized the impressive versatility of the examined PCl3 monomers to participate in σ-hole and lp-hole pnicogen interactions. Crucially, the sizes of σ-hole and lp-hole were enlarged under the influence of the positively directed EEF and decreased in the case of reverse direction. Interestingly, the energetic quantities unveiled more favorability of the σ-hole···lp-hole configuration of the pnicogen···pnicogen homodimers, with significant negative interaction energies, than σ-hole···σ-hole and lp-hole···lp-hole configurations. Quantum theory of atoms in molecules and noncovalent interaction index analyses were adopted to elucidate the nature and origin of the considered interactions, ensuring their closed shell nature and the occurrence of attractive forces within the studied homodimers. Symmetry-adapted perturbation theory-based energy decomposition analysis alluded to the dispersion force as the main physical component beyond the occurrence of the examined interactions. The obtained findings would be considered as a fundamental underpinning for forthcoming studies pertinent to chemistry, materials science, and crystal engineering.

Type I-IV Halogen⋯Halogen Interactions: A Comparative Theoretical Study in Halobenzene⋯Halobenzene Homodimers

In the current study, unexplored type IV halogen⋯halogen interaction was thoroughly elucidated, for the first time, and compared to the well-established types I−III interactions by means of the second-order Møller−Plesset (MP2) method. For this aim, the halobenzene⋯halobenzene homodimers (where halogen = Cl, Br, and I) were designed into four different types, parodying the considered interactions. From the energetic perspective, the preference of scouted homodimers was ascribed to type II interactions (i.e., highest binding energy), whereas the lowest binding energies were discerned in type III interactions. Generally, binding energies of the studied interactions were observed to decline with the decrease in the σ-hole size in the order, C6H5I⋯IC6H5 > C6H5Br⋯BrC6H5 > C6H5Cl⋯ClC6H5 homodimers and the reverse was noticed in the case of type IV interactions. Such peculiar observations were relevant to the ample contributions of negative-belt⋯negative-belt interactions within the C6H5Cl⋯ClC6H5 homodimer. Further, type IV torsional trans → cis interconversion of C6H5X⋯XC6H5 homodimers was investigated to quantify the π⋯π contributions into the total binding energies. Evidently, the energetic features illustrated the amelioration of the considered homodimers (i.e., more negative binding energy) along the prolonged scope of torsional trans → cis interconversion. In turn, these findings outlined the efficiency of the cis configuration over the trans analog. Generally, symmetry-adapted perturbation theory-based energy decomposition analysis (SAPT-EDA) demonstrated the predominance of all the scouted homodimers by the dispersion forces. The obtained results would be beneficial for the omnipresent studies relevant to the applications of halogen bonds in the fields of materials science and crystal engineering.

Stand up for Electrostatics: The Disiloxane Case

The basicity of the simplest silicone, disiloxane (H 3 Si-O-SiH 3 ), is strongly affected by the Si-O-Si angle (α). We use high-level ab initio MP2/aug'-cc-pVTZ calculations and the molecular electrostatic potential (MEP) to analyze the relationship between the increase in basicity and the reduction of α. Our results clearly point out that this increase can be explained through the MEP, as the interactions between oxygen from disiloxane and the acceptors are mostly electrostatic. Furthermore, the effect of α on the tetrel bond between disiloxane and several Lewis bases can again be rationalized using the MEP. Finally, we explore the cooperativity throughout α for ternary complexes where disiloxane simultaneously interacts with a Lewis acid and a Lewis base. Both non-covalent interactions remain cooperative for all α values, although the largest cooperativity effects are not always those maximizing the binding energy in the binary complexes. Overall, the MEP remains a powerful predictor for noncovalent interactions.

Unusual chalcogen⋯chalcogen interactions in like⋯like and unlike YCY⋯YCY complexes (Y = O, S, and Se)

Chalcogen⋯chalcogen interactions were investigated within four types of like⋯like and unlike YCY⋯YCY complexes (where Y = O, S, or Se). A plethora of quantum mechanical calculations, including molecular electrostatic potential (MEP), surface electrostatic potential extrema, point-of-charge (PoC), quantum theory of atoms in molecules (QTAIM), noncovalent interaction (NCI), and symmetry-adapted perturbation theory-based energy decomposition analysis (SAPT-EDA) calculations, were executed. The energetic findings revealed a preferential tendency of the studied chalcogen-bearing molecules to engage in type I, II, III, or IV chalcogen⋯chalcogen interactions. Notably, the selenium-bearing molecules exhibited the most potent ability to favorably participate in all the explored chalcogen⋯chalcogen interactions. Among like⋯like complexes, type IV interactions showed the most favorable negative binding energies, whereas type III interactions exhibited the weakest binding energies. Unexpectedly, oxygen-containing complexes within type IV interactions showed an alien pattern of binding energies that decreased along with an increase in the chalcogen atomic size level. QTAIM analysis provided a solo BCP, via chalcogen⋯chalcogen interactions, with no clues as to any secondary ones. SAPT-EDA outlined the domination of the explored interactions by the dispersion forces and indicated the pivotal shares of the electrostatic forces, except type III σ-hole⋯σ-hole and di-σ-hole interactions. These observations demonstrate in better detail all the types of chalcogen⋯chalcogen interactions, providing persuasive reasons for their more intensive use in versatile fields related to materials science and drug design.

R•-hole interactions of group IV-VII radical-containing molecules: A comparative study

For the first time, the potentiality of the sp²-hybridized group IV-VII radical (R^•)-containing molecules to participate in R^•-hole interactions was comparatively assessed using ^•SiF_3,^•POF₂, ^•SO₂F, and ^•ClO₃ models in the trigonal pyramidal geometry. In that spirit, a plethora of quantum mechanical calculations was performed at the MP2/aug-cc-pVTZ level of theory. According to the results, all the investigated R^•-containing molecules exhibited potent versatility to engage in R^•-hole … Lewis base interactions with significant negative binding energies for the NCH-based complexes. The strength of R^•-hole interactions was perceived to obey the ^•ClO₃ … > ^•SO₂F … > ^•POF₂ … > ^•SiF₃ … Lewis base order, outlining an inverse correlation between the binding energy and the atomic size of the R^•-hole donor. Benchmarking of the binding energy at the CCSD/CBS(T) computational level was executed for all the explored interactions and addressed an obvious similarity between the MP2 and CCSD energetic findings. QTAIM analysis critically unveiled the closed-shell nature of the explored R^•-hole interactions. SAPT-EDA proclaimed the reciprocal contributions of electrostatic and dispersion forces to the total binding energy. These observations demonstrate in better detail the nature of R^•-hole interactions, leading to a convincing amelioration for versatile fields relevant to materials science and drug design.

Effect of External Electric Field on Tetrel Bonding Interactions in (FTF3···FH) Complexes (T = C, Si, Ge, and Sn)

A quantum chemical study was accomplished on the σ-hole interactions of the barely explored group IV elements, for the first time, in the absence and presence of the positively and negatively directed external electric field (EEF). The analyses of molecular electrostatic potential addressed the occurrence of the σ-hole on all the inspected tetrel atoms, confirming their salient versatility to engage in σ-hole interactions. MP2 energetic findings disclosed the occurrence of favorable σ-hole interactions within the tetrel bonding complexes. The tetrel bonding interactions became stronger in the order of C < Si < Ge < Sn for F-T-F3···FH complexes with the largest interaction energy amounting to -19.43 kcal/mol for the optimized F-Sn-F3···FH complex under the influence of +0.020 au EEF. The interaction energy conspicuously evolved by boosting the magnitude of the positively directed EEF value and declining the negatively directed EEF one. The decomposition analysis for the interaction energies was also executed in terms of symmetry-adapted perturbation theory, illuminating the dominant electrostatic contribution to all the studied complexes' interactions except carbon-based interactions controlled by dispersion forces. The outcomes that emerged from the current work reported significantly how the direction and strength of the EEF affect the tetrel-bonding interactions, leading to further improvements in the forthcoming studies of supramolecular chemistry and materials science.

On the Potentiality of X-T-X3 Compounds (T = C, Si, and Ge, and X = F, Cl, and Br) as Tetrel- and Halogen-Bond Donors

The versatility of the X-T-X3 compounds (where T = C, Si, and Ge, and X = F, Cl, and Br) to participate in tetrel- and halogen-bonding interactions was settled out, at the MP2/aug-cc-pVTZ level of theory, within a series of configurations for (X-T-X3)2 homodimers. The electrostatic potential computations ensured the remarkable ability of the investigated X-T-X3 monomers to participate in σ-hole halogen and tetrel interactions. The energetic findings significantly unveil the favorability of the tetrel···tetrel directional configuration with considerable negative binding energies over tetrel···halogen, type III halogen···halogen, and type II halogen···halogen analogs. Quantum theory of atoms in molecules and noncovalent interaction analyses were accomplished to disclose the nature of the tetrel- and halogen-bonding interactions within designed configurations, giving good correlations between the total electron densities and binding energies. Further insight into the binding energy physical meanings was invoked through using symmetry-adapted perturbation theory-based energy decomposition analysis, featuring the dispersion term as the most prominent force beyond the examined interactions. The theoretical results were supported by versatile crystal structures which were characterized by the same type of interactions. Presumably, the obtained findings would be considered as a solid underpinning for future supramolecular chemistry, materials science, and crystal engineering studies, as well as a fundamental linchpin for a better understanding of the biological activities of chemicals.

Study of Beryllium, Magnesium, and Spodium Bonds to Carbenes and Carbodiphosphoranes

The aim of this article is to present results of theoretical study on the properties of C⋯M bonds, where C is either a carbene or carbodiphosphorane carbon atom and M is an acidic center of MX2 (M = Be, Mg, Zn). Due to the rarity of theoretical data regarding the C⋯Zn bond (i.e., the zinc bond), the main focus is placed on comparing the characteristics of this interaction with C⋯Be (beryllium bond) and C⋯Mg (magnesium bond). For this purpose, theoretical studies (ωB97X-D/6-311++G(2df,2p)) have been performed for a large group of dimers formed by MX2 (X = H, F, Cl, Br, Me) and either a carbene ((NH2)2C, imidazol-2-ylidene, imidazolidin-2-ylidene, tetrahydropyrymid-2-ylidene, cyclopropenylidene) or carbodiphosphorane ((PH3)2C, (NH3)2C) molecule. The investigated dimers are characterized by a very strong charge transfer effect from either the carbene or carbodiphosphorane molecule to the MX2 one. This may even be over six times as strong as in the water dimer. According to the QTAIM and NCI method, the zinc bond is not very different than the beryllium bond, with both featuring a significant covalent contribution. However, the zinc bond should be definitely stronger if delocalization index is considered.

Cospatial σ-Hole and Lone Pair Interactions of Square-Pyramidal Pentavalent Halogen Compounds with π-Systems: A Quantum Mechanical Study

In the spirit of the mounting interest in noncovalent interactions, the present study was conducted to scrutinize a special type that simultaneously involved both σ-hole and lone pair (lp) interactions with aromatic π-systems. Square-pyramidal pentavalent halogen-containing molecules, including X-Cl-F4, F-Y-F4, and F-I-X4 compounds (where X = F, Cl, Br, and I and Y = Cl, Br, and I) were employed as σ-hole/lp donors. On the other hand, benzene (BZN) and hexafluorobenzene (HFB) were chosen as electron-rich and electron-deficient aromatic π-systems, respectively. The investigation relied upon a variety of quantum chemical calculations that complement each other. The results showed that (i) the binding energy of the X-Y-F4···BZN complexes increased (i.e., more negative) as the Y atom had a larger magnitude of σ-hole, contrary to the pattern of X-Y-F4···HFB complexes; (ii) the interaction energies of X-Y-F4···BZN complexes were dominated by both dispersion and electrostatic contributions, while dispersive interactions dominated X-Y-F4···HFB complexes; and (iii) the X4 atoms in F-I-X4···π-system complexes governed the interaction energy pattern: the larger the X4 atoms were, the greater the interaction energies were, for the same π-system. The results had illuminating facets in regard to the rarely addressed cases of the σ-hole/lp contradictory scene.

Comparison of ±σ-hole and ±R˙-hole interactions formed by tetrel-containing complexes: a computational study

For the first time, unconventional ±R˙-hole interactions were unveiled in tetrel-containing complexes. The nature and characteristics of ±R˙-hole interactions were explored relative to their ±σ-hole counterparts for ˙TF3⋯ and W-T-F3⋯B/R˙/A complexes (where T = C, Si, and Ge, W = H and F, B = Lewis bases, R˙ = free radicals, and A = Lewis acids). In an effort to thoroughly investigate such interactions, a plethora of quantum mechanical calculations, including molecular electrostatic potential (MEP), maximum positive electrostatic potential (V s,max), point-of-charge (PoC), interaction energy, symmetry adapted perturbation theory (SAPT), and reduced density gradient-noncovalent interaction (RDG-NCI) calculations, were applied. The most notable findings to emerge from this study are that (i) from the electrostatic perspective, the molecular stabilization energies of ˙TF3 and W-T-F3 monomers became more negative as the Lewis basicity increased, (ii) the most stable complexes were observed for the ones containing Lewis bases, forming -σ-hole and -R˙-hole interactions, and the interaction energies systematically increased in the order H-T-F3⋯B < ˙TF3⋯B < F-T-F3⋯B, (iii) contrariwise, the +σ-hole and +R˙-hole interactions with Lewis acids are more energetically favorable in the order F-T-F3⋯A < ˙TF3⋯A < H-T-F3⋯A, and (iv) generally, the dispersion force plays a key role in stabilizing the tetrel-containing complexes, jointly with the electrostatic and induction forces for the interactions with Lewis bases and acids, respectively. Concretely, the findings presented in this paper add to our understanding of the characteristics and nature of such intriguing interactions.

External electric field effects on the σ-hole and lone-pair hole interactions of group V elements: a comparative investigation

σ-hole and lone-pair (lp) hole interactions of trivalent pnicogen-bearing (ZF3) compounds were comparatively scrutinized, for the first time, under field-free and external electric field (EEF) conditions. Conspicuously, the sizes of the σ-hole and lp-hole were increased by applying an EEF along the positive direction, while the sizes of both holes decreased through the reverse EEF direction. The MP2 energetic calculations of ZF3⋯FH/NCH complexes revealed that σ-holes exhibited more impressive interaction energies compared to the lp-holes. Remarkably, the strengths of σ-hole and lp-hole interactions evolved with the increment of the positive value of the considered EEF; i.e., the interaction energy increased as the utilized EEF value increased. Unexpectedly, under field-free conditions, nitrogen-bearing complexes showed superior strength for their lp-hole interactions than phosphorus-bearing complexes. However, the reverse picture was exhibited for the interaction energies of nitrogen- and phosphorus-bearing complexes interacting within lp-holes by applying the high values of a positively directed EEF. These results significantly demonstrate the crucial influence of EEF on the strength of σ-hole and lp-hole interactions, which in turn leads to an omnipresent enhancement for variable fields, including biological simulations and material science.