Quantitative characterization of the global electrophilicity power of common diene/dienophile pairs in Diels–Alder reactions

DFT Study of Some Copper Complexes and Their Detection Limit

A theoretical investigation was probed to shed light on the correlation between low detection limit (LOD) in AdSV technique and metal trace complexes stability energy. The study was conducted by means of DFT calculations of copper traces complexation by using three different organic molecules as chelating agents, such as: morin, red pyrogallol and thymolphtalexone. The quantum chemistry calculations were carried out at the B3LYP/6-31G(d) level implemented in Gaussian 09 program package. The results of the electrophilicity index ω indicate that all the studied molecules have a tendency to exchange electron with copper. The negative values of free energy G and enthalpy H show that the complexation reactions are spontaneous in nature and exothermic. According to DFT calculations, copper-red pyrogallol complex with better detection limit (0.07 ng•mL-1) has the lowest total energy (-5100.213 a.u.). Thus, there is a very strong relationship between the total energy of the three complexes and their detection limits in AdSV technique. Hence, the more stable complex has the better detection limit value.

Synthesis, virtual screening and computational approach of a quinoxaline derivative as potent anti-HIV agent targeting the reverse transcriptase enzyme

Infection by the human immunodeficiency virus still represents a continuous serious concern and a global threat to human health. Due to the appearance of multi-resistant virus strains and the serious adverse side effects of the antiretroviral therapy administered, there is an urgent need for the development of new treatment agents that are more active, less toxic, and with increased tolerability to mutations. Quinoxaline derivatives are a class of heterocyclic compounds with a wide range of organic and remedial applications. In addition, they are known to significantly inhibit HIV reverse transcriptase (RT) and HIV replication in cell cultures. For these reasons, we are investigating the synthesis and computational studies of quinoxaline derivatives with a focus on their effects on the HIV RT enzyme, and we present here the structure of one such molecule, methyl 2-[(2E)-3-oxo-1,2,3,4-tetrahydroquinoalin-2-ylidene] acetate, which was confirmed by X-ray diffraction studies. In the crystal, N-H···O and C-H···O hydrogen bonds form ribbons whose mean planes are inclined to (111) by 25.69(8)°. The ribbons are formed into stacks by C-H···π(ring) interactions and π-stacking interactions between carbonyl groups. The Hirshfeld surface map allows us to understand the nature of interactions in the contribution to crystal packing. A density functional theory (DFT) calculation was performed to optimize the geometrical parameters and then they were compared with the solid-state phase. The molecular electrostatic potential map displays reactive sites on the surface, which are responsible for intermolecular interaction in the chemical species. Computational molecular docking, in addition to molecular dynamics simulations and MMGB/PBSA binding energy techniques, was used to assess the affinity of the molecule for the HIV reverse transcriptase enzyme. The new quinoxaline derivative is more powerful in terms of binding affinity and binding conformation stability with the HIV reverse transcriptase enzyme, which suggests the molecule is a good candidate for further biological optimization.Communicated by Ramaswamy H. Sarma.

Molecular docking-assisted investigation of Pd(II) complexes carrying "SNS" pincer-type pyridine-thioether ligand as potential drug candidates against COVID-19

The coronavirus disease 2019 (COVID-19) pandemic has posed a threat to public health throughout the world since the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) was discovered in late 2019. Since the beginning of the pandemic, scientists have done a tremendous amount of work in this area. However, among these studies, the investigation of the effect of newly synthesized compounds against coronavirus is rather weak. Examining the newly synthesized compounds with a computer-assisted molecular docking study provides quite an advantage in terms of the estimation and analysis of the biochemical activity and binding affinity of existing synthesized compounds against a biological target in a labor, time, and cost-saving way. In this study, the SNS pincer type 2,6-bis[[(4-methylphenyl)thio]methyl]pyridine ligand(L) (1) and its novel Pd(II) complexes ([Pd(κ2-L)(OAc)2]·3H2O (2) and [Pd(κ2-L)Cl2]·3H2O (3)) were synthesized and characterized by using FT-IR, UV-Vis, NMR, mass and elemental analysis techniques. The synthesized Pd complexes exhibited a square planar structure. The compounds were found to have non-electrolytic behavior. In the meantime, in silico investigations have defined and justified interaction processes between these molecules and Pd(II) at the atomic level. Furthermore, using molecular docking against target proteins of SARS-CoV-2, the efficiency of the SNS pincer type ligand and its Pd (II) complexes produced was studied and discussed for the first time. The experimental data has been supported and illuminated using computational visual methods and molecular docking, and the findings produced indicate compatibility. The binding energy values of the relevant compounds on the four protein model structures of SARS-CoV-2 (Main Protease, Papain-like protease, RdRp without RNA, and RdRp with RNA) are represented. Compound 2 ([Pd(κ2-L)(OAc)2]·3H2O) is the structure that exhibits the highest biochemical activity. According to all of the docking studies, Papain-like protease is the SARS-CoV-2 protein with which the three compounds exhibit mutual interaction. The compound 2 structure, in particular, is the most effective in terms of structural and interaction with the targets, as well as binding orientations.

Mechanistic aspects of the synthesis of seven-membered internal nitronates via stepwise [4 + 3] cycloaddition involving conjugated nitroalkenes: Molecular Electron Density Theory computational study

The density functional theory computational study indicates the possibility of the synthesis of seven-membered internal nitronates via cycloaddition reactions involving Z-C-aryl-N-methylnitrones and E-2-aryl-1-cyano-1-nitroethenes. The detailed exploration of the reaction paths indicates a polar, stepwise reaction mechanism through the zwitterionic intermediate. Using bonding evolution theory (BET), we have deciphered the molecular mechanism of the [4 + 3] cycloaddition reaction between E-2-phenyl-1-cyano-1-nitroethene and Z-C-phenyl-N-methylnitrone. The BET study has revealed that the formation of two CO single bonds takes place in the same way, through the depopulation of NC and CC bonding regions and monosynaptic basins, respectively. The first O1C7 single bond was formed in the sixth phase, while the second C3O4 bond was formed in the last ninth phase.

Virtual Prospection of Marine Cyclopeptides as Therapeutics by Means of Conceptual DFT and Computational ADMET

Bioactive peptides are chemical compounds created through the covalent bonding of amino acids, known as amide or peptide bonds. Due to their unusual chemistry and various biological effects, marine bioactive peptides have garnered considerable research. The effectiveness of a bioactive marine peptide is attributed to its structural features, such as amino acid content and sequence, which vary depending on the degree of action. Cyclic peptides combine several favorable properties such as good binding affinity, target selectivity and low toxicity that render them an attractive modality for the development of therapeutics. The apratoxins are a class of molecules formed by a series of cyclic depsipeptides with potent cytotoxic activities. The objective of this research is to pursue a computational prospection of the molecular structures and properties of several cylopeptides of marine origin with potential therapeutic applications. The methodology will be based on the determination of the chemical reactivity descriptors of the studied molecules through the consideration of the Conceptual DFT model and validation of a particular model chemistry, MN12SX/Def2TZVP/H2O. These studies will be complemented by a determination of the pharmacokinetics and ADMET parameters by resorting to certain cheminformatics tools.

Identification of a 3-(5-methyl-2-thiazolylamino)phthalide as a new minor groove agent

A new 3-(5-methyl-2-thiazolylamino)phthalide molecule, 3-((5-methylthiazol-2-yl)amino)isobenzofuran-1(3H)-one, was synthesized and characterized experimentally by FT-IR, NMR, UV-Vis, and single-crystal X-ray analysis and theoretically by quantum chemical calculations. The single-crystal X-ray studies revealed that the compound crystallizes in the monoclinic space group P-2_1/c with unit-cell parameters a = 8.0550(6) Å, b = 6.1386(3) Å, c = 23.3228(18) Å, β = 97.724(6)° and Z = 4. Optimized geometries and the vibrational frequencies were studied at the density functional theory (DFT) level by using the hybrid functional B3LYP with a 6-311 G (d,p) basis set. The title compound was evaluated for its anti-quorum sensing (anti-QS) activity on Chromobacterium violaceum 12472 and additionally for its antibacterial activity against Staphylococcus aureus 29213, Staphylococcus epidermidis 12228, Pseudomonas aeruginosa 27853, Escherichia coli 25922, and Proteus mirabilis 14153. The lowest MIC value was 0.24 μg/mL for S. aureus 29213 and the highest MIC value was 30.75 μg/mL for E. coli 25922. While anti-bacterial activity was observed in those other than the S. epidermidis and P. Mirabilis, anti-QS activity wasn't detected. Investigations on dsDNA binding affinity indicate that the title compound binds to dsDNA via the groove binding mode. Molecular docking calculations and molecular dynamics simulations results showed also that the title compound prefers binding to the minor groove of dsDNA and remains stable in the minor groove throughout the molecular dynamics simulation.Communicated by Ramaswamy H. Sarma.

Computational Chemistry Study of Natural Apocarotenoids and Their Synthetic Glycopeptide Conjugates as Therapeutic Drugs

The objective of the research to be presented in the chapter is the determination of the chemical reactivity properties of some natural apocarotenoids and their synthetic glycopeptide conjugates that could have the ability to inhibit SARS-CoV-2 replication. The study will be based on the consideration of the Conceptual DFT branch of Density Functional Theory (DFT) through the consideration of particular successful model chemistry which has been demonstrated as satisfying the Janak and Ionization Energy theorems within Generalized Gradient Approximation (GGA) theory. The research will be complemented by a report of the ADMET and pharmacokinetic properties hoping that this information could be of help in the development of new pharmaceutical drugs for fighting COVID-19.

A DFT analysis of electronic, reactivity, and NLO responses of a reactive orange dye: the role of Hartree-Fock exchange corrections

An experimental and theoretical study based on DFT/TD-DFT approximations is presented to understand the nature of electronic excitations, reactivity, and nonlinear optical (NLO) properties of reactive orange 16 dye (RO16), an azo chromophore widely used in textile and pharmacological industries. The results show that the solvent has a considerable influence on the electronic properties of the material. According to experimental results, the absorption spectrum is formed by four intense transitions, which have been identified as [Formula: see text] states using TD-DFT calculations. However, the TD-DFT results reveal a weak [Formula: see text] in the low-lying spectral region. Continuum models of solvation indicate that these states suffer from bathochromic (ca. 15 nm) and hypsochromic shifts (ca. 4 nm), respectively. However, the expected blue shift for the absorption [Formula: see text] is only described using long-range or dispersion-corrected DFT methods. RO16 is classified as a strong electrophilic system, with electrophilicity ω > 1.5 eV. Concerning the nucleophilicity parameter (N), from vacuum to solvent, the environment is active and changes the nucleophilic status from strong to moderate nucleophile (2.0 ≤ N ≤ 3.0 eV). The results also suggest that all electrical constants are strongly dependent on long-range and Hartree-Fock exchange contributions, and the absence of these interactions gives results far from reality. In particular, the results for the NLO response show that the chromophore presents a potential application in this field with a low refractive index and first hyperpolarizability ca. 214 times bigger than the value usually reported for urea (β = 0.34 × 10^- 30 esu), which is a standard NLO material. Concerning the solvent effects, the results indicate that the polarizability increases [Formula: see text] esu from gas to solvent while the first hyperpolarizability is calculated as [Formula: see text] esu, ca. 180%, regarding the vacuum. The results suggest RO16 is a potential compound in NLO applications. Graphical Abstract The frontier molecular orbitals, and the inverse relation between the energy-gap (E_gap) and the first hyperpolarizability (β).

New 1,2,3-Triazoles from (R)-Carvone: Synthesis, DFT Mechanistic Study and In Vitro Cytotoxic Evaluation

Aseries of novel 1,4-disubstituted 1,2,3-triazoles were synthesized from an (R)-carvone terminal alkyne derivative via a Cu (I)-catalyzed azide–alkyne cycloaddition reaction using CuSO₄,5H₂O as the copper (II) source and sodium ascorbate as a reducing agent which reduces Cu (II) into Cu (I). All the newly synthesized 1,2,3-triazoles 9a–h were fully identified on the basis of their HRMS and NMR spectral data and then evaluated for their cell growth inhibition potential by MTS assay against HT-1080 fibrosarcoma, A-549 lung carcinoma, and two breast adenocarcinoma (MCF-7 and MDA-MB-231) cell lines. Compound 9d showed notable cytotoxic effects against the HT-1080 and MCF-7 cells with IC₅₀ values of 25.77 and 27.89 µM, respectively, while compound 9c displayed significant activity against MCF-7 cells with an IC₅₀ value of 25.03 µM. Density functional calculations at the B3LYP/6-31G* level of theory were used to confirm the high reactivity of the terminal alkyne as a dipolarophile. Quantum calculations were also used to investigate the mechanism of both the uncatalyzed and copper (I)-catalyzed azide–alkyne cycloaddition reaction (CuAAC). The catalyzed reaction gives complete regioselectivity via a stepwise mechanism streamlining experimental observations. The calculated free-energy barriers 4.33 kcal/mol and 29.35 kcal/mol for the 1,4- and 1,5-regioisomers, respectively, explain the marked regioselectivity of the CuAAC reaction.

Conceptual DFT-Based Computational Peptidology, Pharmacokinetics Study and ADMET Report of the Veraguamides A–G Family of Marine Natural Drugs

As a continuation of our research on the chemical reactivity, pharmacokinetics and ADMET properties of cyclopeptides of marine origin with potential therapeutic abilities, in this work our already presented integrated molecular modeling protocol has been used for the study of the chemical reactivity and bioactivity properties of the Veraguamides A–G family of marine natural drugs. This protocol results from the estimation of the conceptual density functional theory (CDFT) chemical reactivity descriptors together with several chemoinformatics tools commonly considered within the process of development of new therapeutic drugs. CP-CDFT is a branch of computational chemistry and molecular modeling dedicated to the study of peptides, and it is a protocol that allows the estimation with great accuracy of the CDFT-based reactivity descriptors and the associated physical and chemical properties, which can aid in determining the ability of the studied peptides to behave as potential useful drugs. Moreover, the superiority of the MN12SX density functional over other long-range corrected density functionals for the prediction of chemical and physical properties in the presence of water as the solvent is clearly demonstrated. The research was supplemented with an investigation of the bioactivity of the molecular systems and their ADMET (absorption, distribution, metabolism, excretion, and toxicity) parameters, as is customary in medicinal chemistry. Some instances of the CDFT-based chemical reactivity descriptors’ capacity to predict the pKas of peptides as well as their potential as AGE inhibitors are also shown.

A machine learning approach for predicting the nucleophilicity of organic molecules

Nucleophilicity provides important information about the chemical reactivity of organic molecules. Experimental determination of the nucleophilicity parameter is a tedious and resource-intensive approach. Herein, we present a novel machine learning protocol that uses key structural descriptors to predict the nucleophilicities of organic molecules, which agree well with the experimental values. A data driven approach was used where quantum mechanical molecular and thermodynamic descriptors from a wide range of structurally diverse nucleophiles and relevant solvents were extracted and modelled using advanced algorithms against the experimentally available nucleophilicity values. Despite the structural diversity of nucleophiles, we are able to achieve statistically robust models with a high predictive power using tree-based and neural network algorithms trained on an in-house developed unique dataset consisting of 752 nucleophilicity values and 27 molecular descriptors.

Modelling the DFT structural and reactivity study of feverfew and evaluation of its potential antiviral activity against COVID-19 using molecular docking and MD simulations

Abstract

The unavailability of a proper drug against SARS-CoV-2 infections and the emergence of various variants created a global crisis. In the present work, we have studied the antiviral behavior of feverfew plant in treating COVID-19. We have reported a systematic in silico study with the antiviral effects of various phytoconstituents Borneol (C₁₀H₁₈O), Camphene (C₁₀H₁₆), Camphor (C₁₀H₁₆O), Alpha-thujene (C₁₀H₁₆), Eugenol (C₁₀H₁₄O), Carvacrol (C₁₀H₁₄O) and Parthenolide (C₁₅H₂₀O₃) of feverfew on the viral protein of SARS-CoV-2. Parthenolide shows the best binding affinity with both main protease (M^pro) and papain-like protease (PL^pro). The molecular electrostatic potential and Mulliken atomic charges of the Parthenolide molecule shows the high chemical reactivity of the molecule. The docking of Parthenolide with PL^pro give score of −8.0 kcal/mol that validates the good binding of Parthenolide molecule with PL^pro. This complex was further considered for molecular dynamics simulations. The binding energy of the complex seems to range in between −3.85 to −11.07 kcal/mol that is high enough to validate the stability of the complex. Free energy decomposition analysis have been also performed to understand the contribution of residues that reside into the binding site. Good binding affinity and reactivity response suggested that Parthenolide can be used as a promising drug against the COVID-19.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1007/s11696-022-02067-6.

A Theoretical Evaluation of the Behavior of Nitrosoamidine upon Reacting with Methoxy Butadiene, as Potential Heterodiene or Heterodienophile

Abstract:

The molecular mechanism of experimentally observed regio- and chemoselectivity of the cycloaddition reaction of nitrosoamidine 1 and 1-methoxy butadiene 2 has been investigated using DFT calculations at M06-2X/cc-pVDZ level. Accordingly, the possible reaction pathways and factors, which govern selectivity are investigated systematically. Analysis of the calculated results show that the most favorable cyclization reaction takes place through the [2+4] endo-proximal pathway, which is under kinetically and thermodynamically controls. Moreover, analysis of the global and local reactivity indices correctly explains the source of the experimentally observed regio- and chemoselectivity. The electron localization function (ELF) analysis of some selected points along the IRC profile of the most preferred pathway suggests that the reaction takes place via a two-stage one-step mechanism. NCI topological analysis of the possible pathways of [2+4] cycloaddition reaction of 1-E and 2-Z reveals the roles of the attractive interactions between reaction sites, the weak noncovalent interactions in the endo approaches, and the repulsive interactions in the regio- and stereoselectivity of the reaction.

Synergistic effect of GO/SrFe12O19 as magnetic hybrid nanocatalyst for regioselective ring-opening of epoxides with amines under eco-friendly conditions

Highly efficient magnetically separable hybrid GO/SrFe12O19 nanocomposite was synthesized, as catalyst for epoxide ring-opening, via dispersing M-type strontium hexaferrite (SrFe12O19) on graphene oxide (GO) sheets.

Unveiling the regioselective synthesis of antiviral 5-isoxazol-5-yl-2´-deoxyuridines from the perspective of a molecular electron density theory

The regioselective synthesis of a potent antiviral sugar nucleoside isoxazole analogue in the [3+2] cycloaddition (32CA) reaction of acetonitrile- -N-oxide (ANO) and acetyl-protected 5-ethynyl-2?-deoxyuridine (EDU) has been studied at the MPWB1K/6-311G(d,p) level within perspective of the molecular electron density theory (MEDT). From an electron localization function (ELF) analysis, ANO is classified as a zwitterionic species devoid of any pseudoradical or carbenoid centre. The ortho regioisomer is energetically preferred over the meta one by the activation enthalpy of 21.7?24.3 kJ mol-1, suggesting complete regioselectivity in agreement with the experiment. The activation enthalpy increases from 53.9 kJ mol-1 in the gas phase to 71.5 kJ mol-1 in water, suggesting more facile reaction in low polar solvents. The minimal global electron density transfer (GEDT) at the TSs suggests non-polar character and the formation of new covalent bonds has not been started at the located TSs, showing non-covalent intermolecular interactions from an atoms-in- -molecules (AIM) study and in the independent gradient model (IGM) isosurfaces. The AIM analysis shows more accumulation of electron density at the C?C interacting region relative to the C?O one, and earlier C?C bond formation is predicted from a bonding evolution theory (BET) study.

Computational study of quinoline-based thiadiazole compounds as potential antileishmanial inhibitors

Leishmaniasis is a severe disease caused by protozoan parasites of the genus Leishmania and it is accountable for sizable morbidity and mortality worldwide.

Unveiling the intramolecular [3 + 2] cycloaddition reactions of C,N-disubstituted nitrones from the molecular electron density theory perspective

The zwitterionic type intramolecular [3 + 2] cycloaddition reactions of nitrones show remarkable substituent effects and varying reactivities in the cases of cyclic and acyclic alkenes.

Three decades of unveiling the complex chemistry of C-nitroso species with computational chemistry

This review revisits the complex reactivity of C-nitroso derivatives through the synergistic combination of computational and synthetic organic chemistry, with an emphasis on the rationalization of mechanisms and selectivities. 

DFT and MD simulation investigation of favipiravir as an emerging antiviral option against viral protease (3CLpro) of SARS-CoV-2

As per date, around 20 million COVID-19 cases reported from across the globe due to a tiny 125 nm sized virus: SARS-CoV-2 which has created a pandemic and left an unforgettable impact on our world. Besides vaccine, medical community is in a race to identify an effective drug, which can fight against this disease effectively. Favipiravir (F) has recently attracted too much attention as an effective repurposed drug against COVID-19. In the present study, the pertinency of F has been tested as an antiviral option against viral protease (3CLpro) of SARS-CoV-2 with the help of density functional theory (DFT) and MD Simulation. Different electronic properties of F such as atomic charges, molecular electrostatic properties (MEP), chemical reactivity and absorption analysis have been studied by DFT. In order to understand the interaction and stability of inhibitor F against viral protease, molecular docking and MD simulation have been performed. Various output like interaction energies, number of intermolecular hydrogen bonding, binding energy etc. have established the elucidate role of F for the management of CoV-2 virus for which there is no approved therapies till now. Our findings highlighted the need to further evaluate F as a potential antiviral against SARS-CoV-2.

Integrated experimental and theoretical insights for Malachite Green Dye adsorption from wastewater using low cost adsorbent

Water pollution is one of the problems that threaten humanity, and to confront it with only experimental procedures is not enough. It is necessary to integrate both practical methods and theoretical calculations to achieve decontamination with the most accurate interpretation. Hence, discussing the experimental mechanism study of Malachite Green (MG) dye adsorption with the help of the application of density functional theory (DFT) calculations is the main goal of this article. The experimental results affirmed that the preparation of γ-Al₂O₃ by precipitation method using (NH₄)₂CO₃ improved the porosity, the surface capability, and the adsorbent capacities (q_max = 210 mg/g) at optimum condition compared with the previous studies. Kinetic and equilibrium studies showed that the adsorption follows the pseudo-second-order model and Freundlich isotherm model, respectively. Also, the calculated and observed thermodynamic parameters exerted positive values of ΔH° and ΔS°, which translates into an endothermic process with increasing disorder of the system. Theoretical calculations at DFT- B3LYP/6-31G (d,P) level of theory were calculated to show the selectivity of using the cationic form of MG in the experimental measurements to find the interaction mechanism. The electronic structure and intramolecular charge transfer of MG, its cationic form and the complex of MG-Al were investigated theoretically at the B3LYP/6-31 G (d,p) level of theory. The equilibrium geometries of MG, its cationic form and the complex of MG-Al were determined, and it was found that these geometries are non-planar. The E_HOMO and E_LUMO energies can be used to calculate the global properties; chemical hardness (η), softness (S) and electronegativity (χ). The calculated non-linear optical parameters (NLO) of the studied compounds, the electronic dipole moment (μ), first-order hyperpolarizability (β), the hyper-Rayleigh scattering (β_HRS) and the depolarization ratio (DR), showed promising optical properties. Finally, the computational and the experimental results indicated that the adsorption efficiency of MG from wastewater was directly associated with the dye electrophilicity power.