Density functional theory study on the ionic liquid pyridinium hydrogen sulfate

Copper-mediated cyclization of thiosemicarbazones leading to 1,3,4-thiadiazoles: Structural elucidation, DFT calculations, in vitro biological evaluation and in silico evaluation studies

Cancer's global impact necessitates innovative and less toxic treatments. Thiosemicarbazones (TSCs), adaptable metal chelators, offer such potential. In this study, we have synthesized N (4)-substituted heterocyclic TSCs from syringaldehyde (TSL1, TSL2), and also report the unexpected copper-mediated cyclization of the TSCs to form thiadiazoles (TSL3, TSL4), expanding research avenues. This work includes extensive characterization and studies such as DNA/protein binding, molecular docking, and theoretical analyses to demonstrate the potential of the as-prepared TSCs and thiadiazoles against different cancer cells. The DFT results depict that the thiadiazoles exhibit greater structural stability and reduced reactivity compared to the corresponding TSCs. The docking results suggest superior EGFR inhibition for TSL3 with a binding constant value of - 6.99 Kcal/mol. According to molecular dynamics studies, the TSL3-EGFR complex exhibits a lower average RMSD (1.39 nm) as compared to the TSL1-EGFR complex (3.29 nm) suggesting that both the thiadiazole and thiosemicarbazone examined here can be good inhibitors of EGFR protein, also that TSL3 can inhibit EGFR better than TSL1. ADME analysis indicates drug-likeness and oral availability of the thiadiazole-based drugs. The DNA binding experiment through absorption and emission spectroscopy discovered that TSL3 is more active towards DNA which is quantitatively calculated with a Kb value of 4.74 × 106 M-1, Kq value of 4.04 × 104 M-1and Kapp value of 5 × 106 M-1. Furthermore, the BSA binding studies carried out with fluorescence spectroscopy showed that TSL3 shows better binding capacity (1.64 × 105 M-1) with BSA protein. All the compounds show significant cytotoxicity against A459-lung, MCF-7-breast, and HepG2-liver cancer cell lines; TSL3 exhibits the best cytotoxicity, albeit less effective than cisplatin. Thiadiazoles demonstrate greater cytotoxicity than the TSCs. Overall, the promise of TSCs and thiadiazoles in cancer research is highlighted by this study. Furthermore, it unveils unexpected copper-mediated cyclization of the TSCs to thiadiazoles.

Meglumine-based Sustainable Three-component Deep Eutectic Solvent Applicable for the Synthesis of Pyrazolo[5,1-b]quinazoline-3-carboxylates as a Sensing Probe for Cu2+ Ions

An unprecedented meglumine-based three-component deep eutectic solvent (3c-DES) (MegPAc) was synthesized using meglumine, p-toluenesulfonic acid (PTSA), and acetic acid as a renewable, and non-toxic solvent. The exploitation of the MegPAc as an eco-friendly reaction media to construct a selective and sensitive small organic molecular sensing probe, namely, pyrazolo[5,1-b]quinazoline-3-carboxylates (PQCs) was executed. Captivatingly, the MegPAc served the dual role of solvent and catalyst, and it delivered the title components with 69-94 % yields within 67-150 minutes. Furthermore, a UV-visible study unfolds the selective detection of Cu2+ ions with our synthetic probe 4 ba and resulted in hypsochromic shift due to electrostatic interactions. Additionally, 1 H NMR titration study and density functional theory (DFT) calculations were performed to attest the binding mechanism of sensing probe 4 ba and Cu2+ ions. Worthy of mention, this protocol unveils the efficacy of meglumine-based 3c-DES for the first time as a bio-renewable system to synthesize the PQCs.

The structure of protic ionic liquids based on sulfuric acid, doped with excess of sulfuric acid or with water

Neutron scattering with isotopic substitution was used to study the structure of concentrated sulfuric acid, and two protic ionic liquids (PILs): a Brønsted-acidic PIL, synthesised using pyridine and excess of sulfuric acid, [Hpy][HSO₄]·H₂SO₄, and a hydrated PIL, in which an equimolar mixture of sulfuric acid and pyridine has been doped with water, [Hpy][HSO₄]·2H₂O. Brønsted acidic PILs are excellent solvents/catalysts for esterifications, driving reaction to completion by phase-separating water and ester products. Water-doped PILs are efficient solvents/antisolvents in biomass fractionation. This study was carried out to provide an insight into the relationship between the performance of PILs in the two respective processes and their liquid structure. It was found that a persistent sulfate/sulfuric acid/water network structure was retained through the transition from sulfuric acid to PILs, even in the presence of 2 moles (∼17 wt%) of water. Hydrogen sulfate PILs have the propensity to incorporate water into hydrogen-bonded anionic chains, with strong and directional hydrogen bonds, which essentially form a new water-in-salt solvent system, with its own distinct structure and physico-chemical properties. It is the properties of this hydrated PIL that can be credited both for the good performance in esterification and beneficial solvent/antisolvent behaviour in biomass fractionation.

Catalyst-Free, Room-Temperature Accessible Regioselective Synthesis of Spiroquinolines and Their Antioxidant Study

An efficient, regioselective, and environmentally benign approach was established using the multicomponent reaction-based synthesis of novel antioxidant spiroquinoline derivatives such as spiro[dioxolo[4,5-g]quinoline], spiro[dioxino[2,3-g]quinoline], and spiro[pyrazolo[4,3-f]quinoline] by reaction of aryl aldehyde, Meldrum's acid, and amine derivatives under an additive-free reaction in aqueous ethanol. Here, two asymmetric carbon centers, three new C-C bonds, and one C-N bond are developed in the final motif. This synthetic methodology offers excellent yields with an easy workup procedure, high diastereoselectivity [d.r. >50:1 (cis/trans)], admirable atom economy, and low E-factor values. Synthesized spiro compounds were investigated for their in vitro antioxidant activity by 2,2-diphenyl-1-picryl-hydrazyl-hydrate (DPPH) and 2,2'-azino-bis(3-ethylbenzothiazoline-6-sulfonic acid) (ABTS) radical scavenging assays. In the ABTS radical scavenging assay, compounds 4d, 4f, and 4l exhibit excellent potency, and in the DPPH radical scavenging assay, compounds 4a, 4d, 4f, and 4g, exhibit excellent potency.

Synthesis, crystal structure and in silico studies of novel 2,4-dimethoxy-tetrahydropyrimido[4,5-b]quinolin-6(7H)-ones

Herein, acetic acid mediated multicomponent synthesis of novel 2,4-dimethoxy-tetrahydropyrimido[4,5-b]quinolin-6(7H)-one (2,4-dimethoxy-THPQs) was reported. Single-crystal XRD analysis of four newly developed crystals of 2,4-dimethoxy-THPQs and their DFT study were also reported. The structure of all molecules was optimized using DFT B3LYP/6-31G(d) level and compared with the corresponding single-crystal XRD data. As a result, the theoretical and experimental geometrical parameters (bond lengths and bond angles) were found to be in good agreement. Frontier molecular orbital (FMO) and molecule electrostatic potential (MEP) analyses were used to investigate the physicochemical properties and relative reactivity of 2,4-dimethoxy-THPQs. The formation of strong C–H⋯O and N–H⋯O interaction was investigated by Hirshfeld analysis. Furthermore, electronic charge density concentration in 2,4-dimethoxy-THPQs was analysed by the Mulliken atomic charges which helps to predict the ability of 2,4-dimethoxy-THPQs to bind in the receptor. The molecular docking of the crystal structure of 2,4-dimethoxy-THPQs in the main protease (M^pro) of SARS-CoV-2 suggested that all four 2,4-dimethoxy-THPQs efficiently docked in M^pro. Furthermore, 2,4-dimethoxy-THPQs with a 3-chloro substitution in the phenyl ring have the highest binding affinity because of the additional formation of halogen bonds and highest dipole moment.

Single-crystal XRD analysis of 2,4-dimethoxy THPQs and their relative reactivity with properties were investigated using DFT calculation. Molecular docking studies show they effectively docked with main protease of SARS-CoV-2.

Cation-anion interaction effect on the nonlinear optical behavior of pyridinium-based ionic liquids

Influence of the cation-anion interaction (expressed as internal hydrogen bonds) on the nonlinear optical properties (dipole moment, polarizability, anisotropy of polarizability and first-order hyperpolarizability) for the ionic liquids pyridinium hydrogen sulfate (PHS) and pyridinium bis(dihydrogen phosphate) (P2HP) was investigated in this paper. It was achieved analyzing vibrational modes (FT-IR) and electronic transitions (UV-vis) in PHS and P2HP. Molecular geometry of the title compounds was optimized by means of density functional theory calculations at B3LYP/6-311++G(2d,2p) level. Structural analysis exposed stronger internal hydrogen bonds (H⁷···O³ and H¹⁸···O⁴) between cationic and anionic fragments in the pyridinium hydrogen sulfate than these (H²⁶···O⁴ and H²⁵···O¹²) in pyridinium bis(dihydrogen phosphate). These observations are confirmed by FT-IR and UV-vis methods, where: (i) the infrared bands due to H···O and N-C interactions in P2HP appear at lower frequencies than the same ones in the spectrum of PHS and (ii) the n → π* electron transitions in PHS are more pronounced in comparison with these in P2HP. Based on the weaker cation-anion interaction in P2HP with respect to PHS, the former one is characterized with a remarkably higher first-order hyperpolarizability value. Hence, the ionic liquid P2HP was referred as a better nonlinear optical material than PHS.

Investigation on electronic structure, vibrational spectra, NBO analysis, and molecular docking studies of aflatoxins and selected emerging mycotoxins against wild-type androgen receptor

The geometry, frontier molecular orbitals (FMOs), vibrational, NBO analysis, and molecular docking simulations of aflatoxins (B1, B2, M1, M2, G1, G2), zearalenone (ZEA) emodin (EMO), alternariol (AOH), alternariol monoethyl ether (AMME), and tenuazonic acid (TeA) mycotoxins have been extensively theoretically studied and discussed based on quantum density functional theory calculations using Gaussian 16 software package. The theoretical computation for the geometry optimization, NBOs, and the molecular docking interaction was conducted using Density Functional Theory with B3LYP/6-31+G(d,p), NBO program, and AutoDock Vina tools respectively. Charge delocalization patterns and second-order perturbation energies of the most interacting natural bond orbitals (NBOs) of these mycotoxins have also been computed and predicted. Interestingly, among the mycotoxins investigated, aflatoxin G1 is seen to give the strongest stabilization energy while Zearalenone shows the highest tendency to accept electron(s) and emodin, an emerging mycotoxin gave the best binding pose within the androgen receptor pocket with a mean binding affinity of -7.40 kcal/mol.

Theoretical investigation of the stability, reactivity, and the interaction of methyl-substituted peridinium-based ionic liquids

This research work focuses on the reactivity, stability, and electronic interaction of pyridinium hydrogen nitrate (PHN)-based ionic liquids and the influence of methyl substituent on this class of ionic liquids: Ortho- (O-MPHN), meta- (M-MPHN), and para- (P-MPHN) substitution. Natural bond orbital (NBO) calculations were performed at the density functional theory (DFT) with Becke’s Lee Yang and Parr functional (B3LYP) methods and DFT/B3LYP/6-311++G(d,p) as basis set using GAUSSIAN 09W and GAUSSVIEW 6.0 software and the most important interaction between donor (Filled Lewis-type NBO’s) and the acceptor (vacant non-Lewis NBOs) were observed. From our natural bond orbital (NBO) result, it could be deduced that the higher the stabilization energy value, the greater the interaction between the donor and acceptor NBOs. The stability of the studied compounds is said to follow the order from O-MPHN > PHN > P-MPHN > M-MPHN based on the hyperconjugative interaction (stabilization energy) of the most significant interaction. The result of the highest occupied molecular orbital (HOMO), shows that PHN has the highest HOMO while the substituted derivatives have similar HOMO values between −7.70 and −7.98 eV thus PHN complex is the best electron donor while the substituted derivatives act as electron acceptors due to the presence of methyl group substituent which is observed to be electron deficient as a result of its withdrawal effect from the aromatic ring. Furthermore, the electron density, real space functions such as energy density and Laplacian of electron density at bond critical point (BCP) of the hydrogen bond interaction of the studied compounds were analyzed using Multifunctional Wavefunction analyzer software version 3.7 and it was observed that the hydrogen at position 6 and oxygen at position 11 (H6–O11) of M-methyl pyridinium nitrate with bond distance of 4.59 (Å) gave binding energy with the strongest electrostatic interaction between the cation and anion of the compounds under investigation. We also observed from our results that, substitution at the ortho position enhances the stability and strengthen the extent of charge transfer. This therefore implies that substitution at ortho position is more favorable for inter- and intramolecular interactions resulting to stabilization of the studied molecules.

A Systematic Study of DFT Performance for Geometry Optimizations of Ionic Liquid Clusters

Clusters of two ion pairs of imidazolium-based ionic liquids were optimized with 43 different levels of theory, including DFT functionals and MP2-based methods combined with varying Dunning's basis sets, and added dispersion corrections. Better preforming DFT functionals were then applied to clusters consisting of four ion pairs. Excellent performance of some DFT functionals for the two ion pair clusters did not always match that of the four ion-paired clusters despite interionic distances remaining constant between the optimized two and four ion-paired clusters of the same ionic liquid. Combinations of DFT functional and basis set such as ωB97X-D/cc-pVDZ, M06-2X/aug-cc-pVDZ, B3LYP-D3/cc-pVTZ, and TPSS-D3/cc-pVTZ gave excellent results for geometry optimization of two ion-paired clusters of imidazolium ionic liquids but gave larger deviations when applied to the four ion-paired clusters of varying ionic liquids. Empirical dispersion corrections were seen to be crucial in correctly capturing correlation effects in the studied ionic liquid clusters, becoming more important in larger clusters. Dunning's double-ζ basis set, cc-pVDZ, is associated with the smallest root mean squared deviations for geometries; however, it also produces the largest deviations in total electronic energies. ωB97X-D and M06-2X produced the best performance with the augmented version of this basis set. The triple-ζ basis set, cc-pVTZ, leads to the best performance of most of the DFT functionals (especially the dispersion-corrected ones) used, whereas its augmented version, aug-cc-pVTZ, was not seen to improve results. The combinations of functional and basis set that gave the best geometry and energetics in both two and four ion-paired clusters were PBE-D3/cc-pVTZ, ωB97X-D/aug-cc-pVDZ, and BLYP-D3/cc-pVTZ. All three combinations are recommended for geometry optimizations of larger clusters of ionic liquids. PBE-D3/cc-pVTZ performed the best with an average deviation of 2.3 kJ mol^-1 and a standard deviation of 3.4 kJ mol^-1 for total electronic energy when applied to four ion-paired clusters. Geometries optimized with FMO2-SRS-MP2/cc-pVTZ produced total energy within 2.0 kJ mol^-1 off the benchmark in two ion-paired clusters, with the cc-pVDZ basis set performing unsurprisingly poorly with the same method. The error increased to 4.8 kJ mol^-1 on average in four ion-paired clusters, with the smallest RMSD deviations in geometries when compared to the benchmark ones. This study is the first report that investigated the performance of DFT functionals for two and four ion-paired clusters of a wide range of ionic liquids consisting of commonly used cations such as pyrrolidinium, imidazolium, pyridinium, and ammonium. It also identified the importance of assessing the performance of quantum chemical methods for ionic liquids on a variety of cation-anion combinations.

A combined DFT and FT-IR study on the surface interactions in alumina supported ionic liquid [H-Pyr]+[HSO4]. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2020;226:117545. [PMID: 31710889 DOI: 10.1016/j.saa.2019.117545]

The surface interactions in [H-Pyr]⁺[HSO₄]^-/γ-Al₂O₃ system (prepared by wetness impregnation method) are investigated theoretically and experimentally. For that purpose, atoms in molecules theory, natural bond orbital and natural charge population analyses in a combination with a vibrational spectroscopy (FT-IR) are used. It is established that a bond formation between the hydrogen sulfate anion and the alumina influences the IL-support interaction in a great extent. However, a support-cation and a cation-anion interaction in the immobilized [H-Pyr]⁺[HSO₄]^- present as well. A comparative analysis between the experimental and the calculated vibrational modes is carried out and a significant number of infrared bands are assigned. The results indicate a good correlation between the experimental and the theoretical IR frequencies. It is found that the aforementioned interactions affected the vibrational frequencies in the supported IL.

Mechanistic investigation of molecular geometry, intermolecular interactions and spectroscopic properties of pyridinium nitrate

The molecular structure, vibrational frequencies of the fundamental modes and electronic transitions of the ionic liquid pyridinium nitrate ([H-Pyr]⁺[NO₃]^-) have been determined by means of density functional theory (DFT) at B3LYP/6-311++G (2d,2p) level. The chemical bonds nature and the intermolecular interactions in the title compound were investigated using the Quantum theory of atoms in molecules and Hirshfeld surface analysis. Natural bond orbital analysis has been performed in order to elucidate the hybridization and delocalization of electron density within the ion pair [H-Pyr]⁺[NO₃]^-. A detailed vibrational spectral analysis was carried out and the assignments of the observed bands have been proposed on the basis of potential energy distribution. A good correlation between experimental and theoretical IR frequencies was observed. To study the charge transfers occurred in the title molecule, UV-vis analysis was conducted.

Electronic structure, reactivity, and Hirshfeld surface analysis of carvone

The density functional theory (at the B3LYP level using 6-311++G(2d,2p) basis set) was used for the investigation of the geometry and electronic properties of the carvone. The electronic properties and chemical activity of the titled compound were investigated by means of several theoretical approaches, molecular electrostatic potential surface, natural bond orbital, and frontier molecular orbital analyses. It was established that the oxygen atom in the structure characterized the electrophilic reactivity; the positive regions are localized on the hydrogen atoms, which can be considered as possible sites for nucleophilic attack. A detailed analysis of the intermolecular interactions via Hirshfeld surface analysis and fingerprint plots revealed that the carvone structure is stabilized mainly by the formation of O. . .H/H. . .O hydrogen bonds. However, close contacts were established between C. . .H/H. . .C and H. . .H contacts.