Molecular electrostatic potential analysis of non-covalent complexes

Unraveling pnicogen bonding cooperativity: Insights from molecular electrostatic potential analysis

A theoretical investigation on the cooperativity of a series of binary, ternary, and quaternary complexes interconnected by pnicogen bonds has been conducted using calculations at the M06-2X/aug-cc-pVTZ level of density functional theory. By measuring changes in the molecular electrostatic potential (MESP) at the nucleus of interacting atoms in all of the complexes, it is possible to quantify the substantial reorganization of the electron density triggered by the formation of pnicogen bonds. The positive change in MESP, indicating a loss of electron density from the donor molecule in a dimer, facilitates the acceptance of electron density from a third molecule, resulting in the formation of a ternary complex with a stronger pnicogen bond compared to the one present in the binary complex. Similarly, the acceptor molecule in a dimer with a negative change in MESP showed an enhanced tendency to donate electron density to an electron-deficient third molecule. The MESP analysis provided valuable insights into the donor/acceptor characteristics of pnicogen bonds within the quaternary complexes. The proposed MESP hypotheses are consistent with the positive cooperativity observed in the pnicogen-bonded clusters. To quantify the changes in MESP, both at the donor atom (ΔVdonor ) and the acceptor atom (ΔVacceptor ), for all pnicogen bonds in the cluster, the total change in MESP (ΔΔVn ) was measured as ΔΔVn  = ∑(ΔVdonor )-∑(ΔVacceptor ). Remarkably, ΔΔVn exhibited a strong linear relationship with the sum of the bond energies of the pnicogen bonds in the cluster. This establishes the MESP analysis as a robust approach for understanding the strength and cooperative behavior of pnicogen-bonded clusters. Additionally, the MESP features provided clear evidence of pnicogen bond formation, further supporting the reliability of this approach.

Wolfium bonds in homodimers of MX4Y (M = Mo, W; X = F, Cl, Br; Y = O, S, Se)

The term "wolfium bond" has been recently introduced to describe the noncovalent attraction between an atom of group 6 and a nucleophile via a σ-hole binding site. Crystal structures commonly contain a motif wherein two MX4Y units are arranged in close proximity, where M represents either Mo or W, and X and Y refer to halogen and chalcogen atoms respectively. DFT calculations were thus applied to a wide range of homodimers of these molecules so as to assess their preferred arrangements, and to characterize the types of bonding that are present in each in a systematic manner. The most stable Dual-X configuration is symmetric and contains a pair of equivalent M⋯X bonds. The interaction energies range from -8 to -29 kcal mol-1, and are largest for X = F, Y = O, and M = W. The X electron donor is replaced by Y, and the two wolfium bonds are reduced to one, in the less stable Mono-Y structure, with interaction energies between -2 and -10 kcal mol-1. There is some question as to whether the weaker bonds of this type constitute true wolfium bonds.

DFT Exploration of Metal Ion-Ligand Binding: Toward Rational Design of Chelating Agent in Semiconductor Manufacturing

Chelating agents are commonly employed in microelectronic processes to prevent metal ion contamination. The ligand fragments of a chelating agent largely determine its binding strength to metal ions. Identification of ligands with suitable characteristics will facilitate the design of chelating agents to enhance the capture and removal of metal ions from the substrate in microelectronic processes. This study employed quantum chemical calculations to simulate the binding process between eleven ligands and the hydrated forms of Ni2+, Cu2+, Al3+, and Fe3+ ions. The binding strength between the metal ions and ligands was quantified using binding energy and binding enthalpy. Additionally, we explored the binding interaction mechanisms and explained the differences in binding abilities of the eleven ligands using frontier molecular orbitals, nucleophilic indexes, electrostatic potentials, and energy decomposition calculations based on molecular force fields. Based on our computational results, promising chelating agent structures are proposed, aiming to guide the design of new chelating agents to address metal ion contamination issues in integrated circuit processes.

Evaluation of anticancer potential of tetracene-5,12-dione (A01) and pyrimidine-2,4-dione (A02) via caspase 3 and lactate dehydrogenase cytotoxicity investigations

Cancer stands as a significant global cause of mortality, predominantly arising from the dysregulation of key enzymes and DNA. One strategic avenue in developing new anticancer agents involves targeting specific proteins within the cancer pathway. Amidst ongoing efforts to enhance the efficacy of anticancer drugs, a range of crucial medications currently interact with DNA at the molecular level, exerting profound biological effects. Our study is driven by the objective to comprehensively explore the potential of two compounds: (7S,9S)-7-[(2R,4S,5S,6S)-4-amino-5-hydroxy-6-methyloxan-2-yl]oxy-6,9,11-trihydroxy-9-(2-hydroxyacetyl)-4-methoxy-8,10-dihydro-7H-tetracene-5,12-dione (A01) and 5-fluoro-1H-pyrimidine-2,4-dione (A02). These compounds have demonstrated marked efficacy against breast and cervical cancer cell lines, positioning them as promising anticancer candidates. In our investigation, A01 has emerged as a particularly potent candidate, with its potential bolstered by corroborative evidence from lactate dehydrogenase release and caspase-3 activity assays. On the other hand, A02 has exhibited remarkable anticancer potential. To further elucidate their molecular mechanisms and interactions, we employed computational techniques, including molecular docking and molecular dynamics simulations. Notably, our computational analyses suggest that the A01-DNA complex predominantly interacts via the minor groove, imparting significant insights into its mechanism of action. While earlier studies have also highlighted the anticancer activity of A01, our research contributes by providing a deeper understanding of its binding mechanisms through computational investigations. This knowledge holds potential for designing more effective drugs that target cancer-associated proteins. These findings lay a robust groundwork for future inquiries and propose that derivatives of A01 could be synthesized as potent bioactive agents for cancer treatment. By elucidating the distinctive aspects of our study's outcomes, we address the concern of distinguishing our findings from those of prior research.

A Structural Approach to the Strength Evaluation of Linear Chalcogen Bonds

The experimental structural features of chalcogen bonding (ChB) interactions in over 34,000 linear fragments R–Ch⋯A (Ch = S, Se, Te; R = C, N, O, S, Se, Te; A = N, O, S, Se, Te, F, Cl, Br, I) were analyzed. The bond distances dR–Ch and the interaction distances dCh⋯A were investigated, and the functions δR–Ch and δCh⋯A were introduced to compare the structural data of R–Ch⋯A fragments involving different Ch atoms. The functions δR−Ch and δCh⋯A were calculated by normalizing the differences between the relevant bond dR–Ch and ChB interaction dCh⋯A distances with respect to the sum of the relevant covalent (rcovR + rcovCh) and the van der Waals (vdW) radii (rvdWCh + rvdWA), respectively. A systematic comparison is presented, highlighting the role of the chalcogen involved, the role of the R atoms covalently bonded to the Ch, and the role of the A species playing the role of chalcogen bond acceptor. Based on the results obtained, an innovative approach is proposed for the evaluation and categorization of the ChB strength based on structural data.

An in silico study of the selective adsorption and separation of CO2 from a flue gas mixture (CH4, CO2, N2) by ZnLi5+ clusters

Due to the increasing concentration of CO₂ in the atmosphere and its negative effect on the environment, selective adsorption of CO₂ from flue gas has become significantly important. In this study, we have considered a Zn-doped lithium cluster, ZnLi₅⁺ cluster, featuring a planar pentacoordinate Zn centre, as a potential candidate for selective CO₂ capture and separation from a flue gas mixture (CH₄, CO₂, N₂). The binding energy calculation and non-covalent interaction study showed that CO₂ molecules bind relatively strongly as compared to N₂ and CH₄ molecules. The metal cluster can bind five CO₂, five CH₄, and four N₂ molecules with average binding energies of -9.2, -4.4, and -6.1 kcal mol^-1, respectively. Decomposition of the binding energy through symmetry-adapted perturbation theory analysis reveals that the electrostatic component plays a major role. The cationic cluster may be a promising candidate for selective CO₂ capture and can be used as a pollution-controlling agent. The calculated adsorption energy of H₂S is quite closer to that of CO₂, suggesting competitive adsorption between CO₂ and H₂S. The adsorption energies of H₂O and NH₃ are higher compared to CO₂, indicating that these gases may be a potential threat to CO₂ capture.

Computational insight into a mechanistic overview of water exchange kinetics and thermodynamic stabilities of bis and tris-aquated complexes of lanthanides

A thorough investigation of Ln³⁺ complexes with more than one inner-sphere water molecule is crucial for designing high relaxivity contrast agents (CAs) used in magnetic resonance imaging (MRI). This study accomplished a comparative stability analysis of two hexadentate (H₃cbda and H₃dpaa) and two heptadentate (H₄peada and H₃tpaa) ligands with Ln³⁺ ions. The higher stability of the hexadentate H₃cbda and heptadentate H₄peada ligands has been confirmed by the binding affinity and Gibbs free energy analysis in aqueous solution. In addition, energy decomposition analysis (EDA) reveals the higher binding affinity of the peada^4- ligand than the cbda^3- ligand towards Ln³⁺ ions due to the higher charge density of the peada^4- ligand. Moreover, a mechanistic overview of water exchange kinetics has been carried out based on the strength of the metal-water bond. The strength of the metal-water bond follows the trend Gd-O47 (w) > Gd-O39 (w) > Gd-O36 (w) in the case of the tris-aquated [Gd(cbda)(H₂O)₃] and Gd-O43 (w) > Gd-O40 (w) for the bis-aquated [Gd(peada)(H₂O)₂]^- complex, which was confirmed by bond length, electron density (ρ), and electron localization function (ELF) at the corresponding bond critical points. Our analysis also predicts that the activation energy barrier decreases with the decrease in bond strength; hence k _ex increases. The ¹⁷O and ¹H hyperfine coupling constant values of all the coordinated water molecules were different, calculated by using the second-order Douglas-Kroll-Hess (DKH2) approach. Furthermore, the ionic nature of the bonding in the metal-ligand (M-L) bond was confirmed by the Quantum Theory of Atoms-In-Molecules (QTAIM) and ELF along with energy decomposition analysis (EDA). We hope that the results can be used as a basis for the design of highly efficient Gd(iii)-based high relaxivity MRI contrast agents for medical applications.

The use of electrostatic potential at nuclei in the analysis of halogen bonding

Molecular electrostatic potential data at interacting nuclei provide strong evidence of bond formation in many intermolecular halogen bonded complexes.

A novel sulfamethoxazole derivative as an inhibitory agent against HSP70: A combination of computational with in vitro studies

In the current study, a novel derivative of sulfamethoxazole (a sulfonamide containing anti-biotic) named ZM-093 (IUPAC name: (E)-4-((4-(bis(2-hydroxyethyl)amino)phenyl)diazenyl)-N-(5-methylisoxazole-3-yl)benzenesulfonamide) was synthesized via common diazotization-coupling reactions from sulfamethoxazole and subsequently characterized through NMR/FT-IR spectroscopy. After evaluation, the compound was geometrically optimized at the DFT level of theory with BL3YP method and 6/31++G (d,p) basis set and from the optimized structure, several molecular descriptors important in the biological reactivity of the compound, such as global reactivity parameters, molecular electrostatic potential, average local ionization energy, and drug-likeness features of the compound were computationally analyzed. The experimental in vitro investigations of the interaction between ZM-093 and heat shock protein 70 (HSP70), a protein that is highly expressed in several types of cancers, exhibited a significant inhibitory effect against the chaperone activity of HSP70 for the titled compound (P-value < 0.01) and the comparison between the experimental studies with the mentioned computational analysis, as well as molecular docking, illustrated that ZM-093 may inhibit HSP70 through binding to its substrate-binding domain. Finally, by taking all the previous results into account, a new method for assessing the inhibitory activity of ligand to HSP70 is introduced based on protonography, a recently developed method that is dependent on the catalytic activity of carbonic anhydrase on polyacrylamide gel electrophoresis.

Competition between Inter and Intramolecular Tetrel Bonds: Theoretical Studies Complemented by CSD Survey

Crystal structures document the ability of a TF₃ group (T=Si, Ge, Sn, Pb) situated on a naphthalene system to engage in an intramolecular tetrel bond (TB) with an amino group on the adjoining ring. Ab initio calculations evaluate the strength of this bond and evaluate whether it can influence the ability of the T atom to engage in a second, intermolecular TB with another nucleophile. A very strong CN^- anionic base can approach the T either along the extension of a T-C or T-F bond and form a strong TB with an interaction energy approaching 100 kcal/mol, although this bond is weakened a bit by the presence of the internal T⋅⋅⋅N bond. The much less potent NCH base engages in a correspondingly longer and weaker TB, less than 10 kcal/mol. Such an intermolecular TB is weakened by the presence of the internal TB, to the point that it only occurs for the two heavier tetrel atoms Sn and Pb.

Crystallographic and Theoretical Evidences of Anion⋅⋅⋅Anion Interaction

Planar (HgCl₃ )^- anions are stacked fairly closely together in a slipped parallel arrangement within several crystal structures. Quantum chemical analysis shows evidence of strong noncovalent spodium bonds between the Hg π-hole of one unit and the Cl atom of an adjacent unit. Anion⋅⋅⋅anion spodium bonds work in tandem with crystal packing forces.

Design and DFT study of nitrogen-rich donor systems for improved photovoltaic performance in dye-sensitized solar cells

Donor modifications, especially through N-annulation, for enhancing the structure–performance relationship of D–π–A systems for DSSC applications.

Computational Study on Metal-Ion-Decorated Prismane Molecules for Selective Adsorption of CO2 from Flue Gas Mixtures

Selective adsorption of CO₂ from flue gas is extremely significant because of its increasing concentration in air and its deleterious effect on the environment. In this work, we have explored metal-ion-bound prismane molecules for selective CO₂ adsorption from the flue gas mixture. The Ca²⁺-bound prismane complex exhibits superior CO₂ selectivity and adsorption capacity. The calculated binding energy and molecular electrostatic potential (MESP) analysis showed that the rectangular face of prismane binds strongly with metal ions as compared to its triangular face. The CBS-QB3 and density functional theory-based functional M06-2X/6-311+G(d) calculations show that the prismane molecule can bind to one Li⁺, K⁺, Mg²⁺, and Ca²⁺ ion with favorable binding energy. The metal-ion-bound prismane complexes have been examined for their CO₂, N₂, and CH₄ adsorption capacity. Prismane-Ca²⁺ can bind with six CO₂ molecules strongly with an average binding energy of -18.1 kcal/mole as compared to six N₂ (-12.6) and five CH₄ (-13.4) gas molecules. The gravimetric density calculated for the CO₂-adsorbed prismane-Ca²⁺ complex has been found to be 69.1 wt %. The discrete hydrocarbon structure for selective separation of CO₂ is rare in the literature and can have potential applications for cost-effective CO₂ capture from the flue gas mixture.

Dominance of unique Pπ phosphorus bonding with π donors: evidence using matrix isolation infrared spectroscopy and computational methodology

Albeit the first account of hypervalentπ interactions has been reported with halogenπ interactions, the feasibility of their extension to other hypervalent atoms as possible Lewis acids is still open. In this work, the role of phosphorus as an acceptor from the π electron cloud (Pπ pnicogen or phosphorus bonding) in PCl3-C2H2 and PCl3-C2H4 heterodimers is explored, by combining matrix isolation infrared spectroscopy with ab initio and DFT computational methodologies. The respective potential energy surfaces of the PCl3-C2H2 and PCl3-C2H4 heterodimers reveal unique minima stabilized by a concert of reasonably strong to weak interactions, of which Pπ phosphorus bonding was energetically dominant. Heterodimers, trimers and tetramers bound primarily by this unique phosphorus bond were generated at low temperatures. The dominance of phosphorus bonding in the PCl3-C2H2 and PCl3-C2H4 heterodimers over other interactions (such as Hπ, HCl, HP, Clπ and lone pair-π interactions) was confirmed and substantiated using extended quantum theory of atoms in molecules, natural bond orbital, electrostatic potential mapping and energy decomposition analyses. The following inferences in correlation with results from non-covalent-interaction analysis offer a complete understanding of the nature of the Pπ phosphorus bonding interactions. The significance of electrostatic forces kinetically favoring the formation of phosphorus bonded heterodimers, in addition to thermodynamic stabilization, is demonstrated experimentally.

Microhydration of a Carbonyl Group: How does the Molecular Electrostatic Potential (MESP) Impact the Formation of (H2O)n:(R2C═O)Complexes?

The presence of a carbonyl group in a molecule usually leads to the identification of a π-hole on the molecular electrostatic potential (MESP) of the species. How does this electrophilic site influence the formation of microhydrated complexes? To address this point, a panel of R₂CO solutes with various MESPs was selected, and we identified the structures and properties of several complexes containing one, two, three and six water molecules. The following solutes were considered in the present study: H₂CO, F₂CO, Cl₂CO,(NC)₂CO and H₂C═CO. Geometry optimizations and frequency calculations were carried out at the LC-ωPBE/6-311++G(d,p) level, with the GD3BJ empirical correction for dispersion. For a number of n water molecules around the R₂CO solute, the structure and the features of the most stable (H₂O)_n:(R₂CO) complexes are highly dependent on the MESP of the isolated R₂CO solute. The formation of pi-hole bondings appears to play a decisive role in the initiation of a three-dimensional organization of water molecules around the solute.

Probing the bent bonds in cyclopropane systems for gas storage and separation process: A computational study

The hydrogen, carbon dioxide, and carbon monoxide gas adsorption and storage capacity of lithium-decorated cyclopropane ring systems were examined with quantum chemical calculations at density functional theory, DFT M06-2X functional using 6-31G(d) and cc-pVDZ basis sets. To examine the reliability of M06-2X DFT functional, a few representative systems are also examined with complete basis set CBS-QB3 method and CCSD-aug-cc-pVTZ level of theory. The cyclopropane systems can bind to one Li⁺ ion; however, the corresponding the methylated systems can bind with two Li⁺ ions. The cyclopropane systems can adsorb six hydrogen molecules with an average binding energy of 3.8 kcal/mol. The binding free energy (ΔG) values suggest that the hydrogen adsorption process is feasible at 273.15 K. The calculation of desorption energies indicates the recyclable property of gas adsorbed complexes. The same number of CO₂ and CO gas molecules can also be adsorbed with an average binding energy of -14.4 kcal/mol and -10.7 kcal/mol, respectively. The carbon dioxide showed ~3-4 kcal/mol better binding energy as compared to carbon monoxide and hence such designed systems can function as a potential candidate for the separation of these flue gas molecules. The nature of interactions in complexes was examined with atoms in molecules analysis revealed the electrostatic nature for the interaction of Li⁺ ion with cyclopropane rings. The chemical hardness and electrophilicity calculations showed that the gas adsorbed complexes are rigid and therefore robust as gas storage materials.

Competition between Intra and Intermolecular Triel Bonds. Complexes between Naphthalene Derivatives and Neutral or Anionic Lewis Bases

: A TrF2 group (Tr = B, Al, Ga, In, Tl) is placed on one of the α positions of naphthalene, and its ability to engage in a triel bond (TrB) with a weak (NCH) and strong (NC-) nucleophile is assessed by ab initio calculations. As a competitor, an NH2 group is placed on the neighboring Cα, from which point it forms an intramolecular TrB with the TrF2 group. The latter internal TrB reduces the intensity of the π-hole on the Tr atom, decreasing its ability to engage in a second external TrB. The intermolecular TrB is weakened by a factor of about two for the smaller Tr atoms but is less severe for the larger Tl. The external TrB can be quite strong nonetheless; it varies from a minimum of 8 kcal/mol for the weak NCH base, up to as much as 70 kcal/mol for CN-. Likewise, the appearance of an external TrB to a strong base like CN- lessens the ability of the Tr to engage in an internal TrB, to the point where such an intramolecular TrB becomes questionable.

On the ability of pnicogen atoms to engage in both σ and π-hole complexes. Heterodimers of ZF2C6H5 (Z = P, As, Sb, Bi) and NH3

When bound to a pair of F atoms and a phenyl ring, a pyramidal pnicogen (Z) atom can form a pnicogen bond wherein an NH₃ base lies opposite one F atom. In addition to this σ-hole complex, the ZF₂C₆H₅ molecule can distort in such a way that the NH₃ approaches on the opposite side to the lone pair on Z, where there is a so-called π-hole. The interaction energies of these π-hole dimers are roughly 30 kcal/mol, much larger than the equivalent quantities for the σ-hole complexes, which are only 4-13 kcal/mol. On the other hand, this large interaction energy is countered by the considerable deformation energy required for the Lewis acid to adopt the geometry necessary to form the π-hole complex. The overall energetics of the complexation reaction are thus more exothermic for the σ-hole dimers than for the π-hole dimers.

Formation of large clusters of CO2 around anions: DFT study reveals cooperative CO2 adsorption

The cooperative O⋯C secondary interactions compensate for the diminishing effect of primary anion⋯C interactions in anionic clusters of CO₂ molecules.

Electrostatics for probing lone pairs and their interactions

The value of the molecular electrostatic potential minimum (V_min ) and its topographical features (position, as well as the eigenvalues and eigenvectors of the corresponding Hessian matrix) are recently proposed as the criteria for characterizing a lone pair (Kumar A. et al., J. Phys. Chem. 2014, A118, 526). This electrostatic characterization of lone pairs is examined for a large number of small molecules employing MP4/6-311++G(d,p)//MP2/6-311++G(d,p) theory. The eigenvector of the Hessian matrix corresponding to its largest eigenvalue (λ_max ), is found to be directed toward the lone pair-bearing-atom, with λ_max showing a strong linear correlation with V_min . Large magnitudes of V_min and λ_max indicate a charge-dense lone pair. The topographical features of V_min are seen to provide insights into the interactive behavior of the molecules with model electrophiles, viz. HF, CO₂ , and Li⁺ . In all the complexes of HF and majority of the other complexes, the interaction energy (E_int ) correlates well with the respective V_min, value, but for some deviations occurring due to other competing secondary interactions. The electrostatic interactions are found to be highly directional in nature as the orientation of interacting atom correlates strongly to the position of lone pair. In summary, the present study on a large number of test molecules shows that electrostatics is able to probe lone pairs in molecules and offers a simple interpretation of chemical reactivity. © 2017 Wiley Periodicals, Inc.

Toward a uniform description of hydrogen bonds and halogen bonds: correlations of interaction energies with various geometric, electronic and topological parameters

Correlations between interaction energies and various structural parameters were established to reveal the differences between hydrogen bonds and halogen bonds.